pith. sign in

arxiv: 2511.04074 · v2 · pith:2B5SNLGUnew · submitted 2025-11-06 · ⚛️ physics.geo-ph

Insights on Numerical Damping Formulations Gained from Calibrating Two-Dimensional Ground Response Analyses at Downhole Array Sites

Nishkarsha Dawadi , Kami Mohammadi , Mohamad M. Hallal , Brady R. Cox This is my paper

Pith reviewed 2026-05-21 20:27 UTC · model grok-4.3

classification ⚛️ physics.geo-ph
keywords ground response analysisnumerical dampingsmall-strain dampingRayleigh dampingempirical transfer functionsdownhole arraysseismic site responsevelocity contrast
0
0 comments X

The pith

Two-dimensional ground response analyses match empirical transfer functions best when using Rayleigh Mass damping with an inflated small-strain damping multiplier.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper examines numerical damping in two-dimensional ground response analyses to better capture seismic wave attenuation at four downhole array sites. It tests different damping formulations and shows that inflating the small-strain damping ratio and applying Rayleigh Mass damping reproduces the observed transfer functions more accurately than conventional methods. A sympathetic reader would care because accurate ground response modeling is essential for predicting earthquake shaking and designing safer structures. The results also indicate that the needed inflation depends on the velocity contrast at the site.

Core claim

Using an inflated D_min multiplier m in 2D GRAs with Rayleigh Mass damping consistently achieved the closest match to empirical transfer functions at three of the four sites while offering faster computational performance. The appropriate m correlates with the site's velocity contrast. Full Rayleigh and Maxwell damping with inflation overdamped higher modes.

What carries the argument

Rayleigh Mass damping formulation applied with a calibrated multiplier m to the small-strain damping ratio D_min in two-dimensional ground response analyses

If this is right

  • Inflated D_min accounts for unmodeled attenuation such as diffraction in 2D models.
  • Rayleigh Mass damping provides better predictions than frequency-independent formulations.
  • The m value correlates with velocity contrast, allowing site-specific adjustments.
  • This combination is computationally faster than alternative damping methods.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • This method may allow simpler 2D models to substitute for more complex 3D simulations in many practical cases.
  • Validating the velocity contrast correlation at new sites would test broader applicability.
  • Adjustments to damping could influence estimates of nonlinear soil response during large earthquakes.

Load-bearing premise

The m multipliers calibrated to match transfer functions at these four specific sites can be generalized to other sites using the observed correlation with velocity contrast.

What would settle it

Conducting the calibration process at an additional downhole array site and verifying if the Rayleigh Mass damping with the corresponding m still provides the closest match to the empirical transfer functions.

Figures

Figures reproduced from arXiv: 2511.04074 by Brady R. Cox, Kami Mohammadi, Mohamad M. Hallal, Nishkarsha Dawadi.

Figure 3
Figure 3. Figure 3: The maximum resolvable frequencies (fs) from the 3D model developed for each site are also shown in [PITH_FULL_IMAGE:figures/full_fig_p011_3.png] view at source ↗
Figure 6
Figure 6. Figure 6 [PITH_FULL_IMAGE:figures/full_fig_p017_6.png] view at source ↗
read the original abstract

Accurately modeling seismic wave attenuation is critical for ground response analyses (GRAs), which aim to replicate local site effects in ground motions. However, theoretical transfer functions (TTFs) from GRAs often overestimate empirical transfer functions (ETFs) when the small-strain damping ratio ($D_{\text{min}}$) is set equal to laboratory measurements. Prior studies addressed this by inflating $D_{\text{min}}$ in one-dimensional (1D) GRAs to account for apparent damping mechanisms such as diffraction and mode conversions that cannot be captured in 1D. Although this approach improved fundamental-mode predictions, it often overdamped higher modes. This study explores more direct modeling of apparent damping using two-dimensional (2D) GRAs at four downhole array sites: Delaney Park (DPDA), I-15 (I15DA), Treasure Island (TIDA), and Garner Valley (GVDA). At each site, three numerical damping formulations, Full Rayleigh, Maxwell, and Rayleigh Mass, were implemented using both conventional $D_{\text{min}}$ and an inflated $D_{\text{min}}$ ($m \times D_{\text{min}}$) obtained from site-specific calibration. Results show that the appropriate $D_{\text{min}}$ multiplier ($m$) correlates with the site's velocity contrast. Using inflated $D_{\text{min}}$, Full Rayleigh and Maxwell damping systematically overdamped higher modes, with Maxwell damping also shifting modal peaks. In contrast, Rayleigh Mass damping consistently achieved the closest match to ETFs at three of the four sites while offering faster computational performance. These findings demonstrate that inflated $D_{\text{min}}$ can represent unmodeled attenuation in 2D GRAs, particularly at sites with low velocity contrast, and that frequency-dependent formulations such as Rayleigh Mass damping can more accurately predict site response than traditional frequency-independent approaches.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript examines three numerical damping formulations (Full Rayleigh, Maxwell, and Rayleigh Mass) in 2D ground response analyses at four downhole array sites (DPDA, I15DA, TIDA, GVDA). It calibrates a site-specific multiplier m on the small-strain damping ratio D_min to improve agreement with empirical transfer functions (ETFs), reports that m correlates with velocity contrast, and concludes that Rayleigh Mass damping with inflated D_min (m × D_min) yields the closest ETF match at three sites while being computationally faster; this is interpreted as evidence that inflated D_min captures unmodeled attenuation, especially at low-contrast sites.

Significance. If the performance ranking and correlation prove robust beyond the calibration set, the work would supply actionable guidance for choosing damping models in 2D GRAs, favoring frequency-dependent options such as Rayleigh Mass for better higher-mode fidelity and efficiency. It would also strengthen the case for treating apparent damping via D_min inflation in low velocity-contrast settings.

major comments (2)
  1. Calibration and evaluation procedure (Methods/Results sections): The multiplier m is obtained by fitting each formulation to the same empirical transfer functions later used to rank performance and declare 'closest match.' This renders the reported superiority of Rayleigh Mass damping with inflated D_min partly tautological rather than an independent test, directly affecting the central claim that the results demonstrate inflated D_min represents unmodeled attenuation.
  2. Correlation and generalization statements (Results/Discussion): The claimed correlation between calibrated m and velocity contrast, together with the performance ranking, rests on observations from only four sites with one exception. No out-of-sample validation, statistical quantification of the trend, or sensitivity checks to layering/ETF quality are reported, so the interpretation that inflated D_min is particularly useful at low-contrast sites cannot yet be generalized.
minor comments (2)
  1. Abstract: Specify which of the four sites was the exception to the 'three of four' claim and briefly note the reason, to avoid ambiguity for readers.
  2. Notation and tables: Add a summary table listing the calibrated m values, velocity contrasts, and key modal-match metrics for each site and formulation to improve traceability.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the thoughtful and constructive comments on our manuscript. We address each major comment below, providing clarifications on our methodology and acknowledging limitations where appropriate. We believe these responses strengthen the presentation of our findings from the available downhole array data.

read point-by-point responses

  1. Referee: Calibration and evaluation procedure (Methods/Results sections): The multiplier m is obtained by fitting each formulation to the same empirical transfer functions later used to rank performance and declare 'closest match.' This renders the reported superiority of Rayleigh Mass damping with inflated D_min partly tautological rather than an independent test, directly affecting the central claim that the results demonstrate inflated D_min represents unmodeled attenuation.

    Authors: We appreciate the referee drawing attention to the calibration procedure. For each damping formulation, m is calibrated independently by minimizing the misfit between the theoretical transfer function and the empirical transfer function at that site. The subsequent comparison then evaluates which formulation, when allowed this same site-specific adjustment to D_min, achieves the smallest residual misfit to the observations. This is a standard approach for comparing model classes under equivalent parameterization flexibility rather than a circular test. The closer match obtained with Rayleigh Mass damping indicates that its frequency-dependent character better accommodates the apparent damping effects present in the data. The interpretation that inflated D_min captures unmodeled attenuation (e.g., diffraction and scattering not explicit in the 2D mesh) is grounded in the systematic pattern that larger m values are required at low velocity-contrast sites, where such effects are expected to be more prominent relative to the modeled 1D-like response. We will revise the Methods and Discussion sections to state this rationale more explicitly and to clarify that the ranking reflects relative model performance under consistent calibration. revision: partial

  2. Referee: Correlation and generalization statements (Results/Discussion): The claimed correlation between calibrated m and velocity contrast, together with the performance ranking, rests on observations from only four sites with one exception. No out-of-sample validation, statistical quantification of the trend, or sensitivity checks to layering/ETF quality are reported, so the interpretation that inflated D_min is particularly useful at low-contrast sites cannot yet be generalized.

    Authors: We agree that the dataset is limited to the four downhole array sites possessing both high-quality empirical transfer functions and detailed velocity profiles. The reported correlation between m and velocity contrast is an empirical observation across these sites (with the noted exception at one site), and we did not conduct formal statistical tests or out-of-sample validation because additional independent sites with comparable data are not currently available. In the revised manuscript we will (i) expand the Discussion to explicitly qualify the correlation as preliminary and hypothesis-generating, (ii) report any sensitivity checks we performed regarding ETF frequency range and layering assumptions, and (iii) temper generalization statements to indicate that the utility of D_min inflation at low-contrast sites is suggested by the present results and warrants testing with future datasets. revision: yes

Circularity Check

1 steps flagged

Site-specific calibration of m to ETFs renders reported 'closest match' partly tautological

specific steps
  1. fitted input called prediction [Abstract]
    "an inflated D_min (m × D_min) obtained from site-specific calibration. Results show that the appropriate D_min multiplier (m) correlates with the site's velocity contrast. Using inflated D_min, Full Rayleigh and Maxwell damping systematically overdamped higher modes... In contrast, Rayleigh Mass damping consistently achieved the closest match to ETFs at three of the four sites"

    m is calibrated directly to the ETFs at each site to improve the match; the subsequent claim that Rayleigh Mass damping with this m 'achieved the closest match' is therefore achieved by construction of the calibration rather than through an independent test or prediction.

full rationale

The paper obtains the D_min multiplier m via site-specific calibration to the empirical transfer functions (ETFs) at each of the four arrays, then reports that Rayleigh Mass damping with this inflated D_min produces the closest match to those same ETFs at three sites. This evaluation step reduces to the fitting process by construction. The observed correlation between m and velocity contrast is an empirical observation from the calibration set but does not constitute an independent prediction. With N=4 and no out-of-sample validation, the ranking of damping formulations rests on in-sample fits rather than falsifiable forecasts. This qualifies as partial circularity under the fitted-input-called-prediction pattern, though the cross-formulation comparison and velocity-contrast trend retain some independent descriptive content.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The central claim depends on the validity of the numerical damping implementations and on the assumption that the four chosen sites are representative. The multiplier m is a fitted parameter with no independent physical derivation.

free parameters (1)
  • D_min multiplier m
    Site-specific scalar obtained by calibration to match empirical transfer functions; value varies with velocity contrast.
axioms (1)
  • domain assumption Numerical damping formulations accurately represent small-strain material damping when properly parameterized.
    Invoked when comparing TTFs to ETFs after applying the formulations.

pith-pipeline@v0.9.0 · 5885 in / 1389 out tokens · 37612 ms · 2026-05-21T20:27:34.440754+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Effects of Varying Incident Wave Inclination and Azimuthal Angles on Multi-Dimensional Ground Response Analyses at the Delaney Park Downhole Array Site

    physics.geo-ph 2026-04 unverdicted novelty 4.0

    Inclination angles up to 15 degrees in 2D and 3D ground response models at Delaney Park produce only minor reductions in transfer function amplitudes near the fundamental frequency with limited improvement over vertic...

Reference graph

Works this paper leans on

2 extracted references · 2 canonical work pages · cited by 1 Pith paper

  1. [1]

    https://doi.org/10.1785/0120200214 Itasca Consulting Group, Inc. (2023). FLAC3D — Fast Lagrangian Analysis of Continua in Three- Dimensions, Ver. 9.0. Itasca. Jackson, T. S. (2024). Evaluating 1D and 2D Small-Strain Ground Response Analyses at the I-15 Downhole Array Using Recorded Aftershocks from the M5.7 2020 Magna, Utah Earthquake. [M.S. Thesis, Utah ...

    work page doi:10.1785/0120200214 2023

  2. [2]

    P., Meijers, P., Bommer, J

    https://doi.org/10.1061/(ASCE)1090-0241(2000)126:5(472) Rodriguez‐Marek, A., Kruiver, P. P., Meijers, P., Bommer, J. J., Dost, B., van Elk, J., & Doornhof, D. (2017). A Regional Site‐Response Model for the Groningen Gas Field. Bulletin of the Seismological Society of America, 107(5), 2067–2077. https://doi.org/10.1785/0120160123 Shibuya, S., Mitachi, T., ...

    work page doi:10.1061/(asce)1090-0241(2000)126:5(472 2000