the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
| 04 Dec 2025
Soil thermal memory regulates event-scale precipitation recycling lag in a dryland environment
Abstract. Precipitation events in dryland environments generate sharp but uneven adjustments in surface and atmospheric conditions. While the atmosphere recovers rapidly from rainfall-induced cooling, the soil retains a substantial portion of the cooling anomaly, creating a land-surface memory. Using multi-year, multi-layer observations from five stations in Ningxia, China, and ERA5 reanalysis, this study investigates how this soil thermal memory timescale (τrec) modulates the timing of recycled moisture return.
Analysis of 112 events reveals a consistent "cold-humid pulse" with rapid atmospheric recovery but slow soil recovery, whose persistence we quantify as τrec (20–80 hours) using an exponential-decay framework. ERA5 diagnostics show the recycled moisture signal peaks 20–40 hours after rainfall, defining a recycling lag (τRR). Event-wise analysis of ten long-duration events reveals a systematic positive correlation (R ≈ 0.57) between τrec and τRR.
Longer soil memory consistently predicts a more delayed recycling peak. We show this relationship is mediated by enhanced moisture-heat feedback (H), where persistent cold soils slow boundary-layer recovery and postpone the reactivation of evaporation. These results identify soil thermal memory as an active regulator; the timing of recycled moisture is not solely an atmospheric process but is partially land-controlled. This work establishes a novel "coupling-memory-recycling" pathway, providing a new mechanism for understanding and modeling dryland precipitation dynamics.
- Preprint
(8506 KB) - Metadata XML
- BibTeX
- EndNote
Status: open (until 03 May 2026)
-
RC1: 'Comment on egusphere-2025-5752', Anonymous Referee #1, 20 Mar 2026
reply
-
AC1: 'Reply on RC1', Ruolin Li, 25 Mar 2026
reply
We are grateful to Referee #1 for a detailed and critical reading. The review led us to uncover an implementation error in the event-scale τ_rec fitting — one that materially affects the central correlation — and prompted us to reconsider how the manuscript frames its contribution. We address each comment below and describe the corresponding revisions.
RC1.1: The methods section is largely based on bullet points… It’s not clear how W_LOC and W_ADV connect to P_ADV/P_LOC. Does P_ADV + P_LOC = P measured at the stations?
AC: The bullet-point structure was a presentational mistake; it fragmented what is fundamentally a continuous derivation. The revised manuscript rewrites this section as prose, with the governing equations laid out explicitly.
The two reservoirs W_LOC and W_ADV track column-integrated water vapour of local and advected origin, respectively, and evolve as:
dW_LOC/dt = E − ρ·P − ∇·(u W_LOC)
dW_ADV/dt = 0 − (1−ρ)·P − ∇·(u W_ADV)
where ρ = W_LOC / (W_LOC + W_ADV) is the instantaneous recycling ratio. Precipitation partitions as P_LOC = ρ·P and P_ADV = (1−ρ)·P, so P_LOC + P_ADV = P at every timestep by construction. Horizontal advection uses an upwind finite-difference scheme applied separately to the u and v components; a CFL-compliant timestep check enforces numerical stability, and non-negativity of both reservoirs is maintained throughout.
On the ERA5 vs. station comparison: we will add a supplementary validation table giving ERA5 total precipitation against station-observed accumulations for each of the ten events, including bias and peak-timing offset. For the domain-averaged values used in the recycling calculation, ERA5 and station observations agree to within 15% in total accumulation for 8 of 10 events. The two exceptions (case2, case7) are the advection-dominated events discussed in RC1.3.
RC1.2: Fig. 4 suggests that the seasonal cycle of temperature is mainly responsible for τ_REC, which isn’t a major surprise.
AC: This is a fair point. The seasonal pattern of τ_rec predominantly reflects temperature-controlled recovery rates — warmer soils rewarm faster after rain; cooler soils sustain the anomaly longer. We do not think this warrants a novelty claim, and Section 4.2 will be reframed accordingly: the correspondence between τ_rec and the expected temperature-controlled seasonal envelope serves as physical validation that the metric is capturing a genuine land-surface signal rather than a fitting artefact.
To separate temperature-independent variability, we will include a partial correlation analysis reporting the residual variance in τ_rec after removing the seasonal T2m signal. This clarifies whether event-scale τ_rec carries information beyond background seasonality.
RC1.3: Claims of a consistent temporal structure of ρ are not the case: 2/10 are negative and 3/10 are greater than 50 hours.
AC: The reviewer is right, and we withdraw the characterisation of τ_RR as showing “consistent temporal structure.” The ten events span −72 to +110 h; two yield negative τ_RR values, meaning the recycling ratio peaks before the precipitation maximum.
We interpret these negative-τ_RR cases as advection-dominated: pre-existing moisture convergence elevates the recycling ratio ahead of the local rainfall peak rather than in response to it. The revised Section 4.3.1 will introduce this classification explicitly — recycling-dominated (τ_RR ≥ 0, n = 8) versus advection-dominated (τ_RR < 0, n = 2) — and use it to stratify all subsequent event-scale analysis. The classification itself is a physically meaningful result: even within a single region and season, individual events can be controlled by fundamentally different moisture pathways.
RC1.4: 6/7 of the τ_REC timescale are equal to 5 — the correlation is therefore totally spurious.
AC: The reviewer has identified a genuine implementation error, and we thank them for catching it. We traced the problem to the upper bound constraint in the optimisation.
The original code fitted the exponential decay to the autocorrelation function (ACF) of T0cm anomalies, with an upper bound of 10 × t_max, where t_max was the maximum ACF lag. Because the ACF was evaluated over a sub-hourly window, t_max ≈ 0.5 h — collapsing the effective upper bound to ~5 h. The optimiser converged to this ceiling for nearly all events. This is a fitting-bound artefact, not a physical result.
We have corrected this by fitting the exponential decay model directly to the post-event T0cm anomaly time series (0–72 h window) with a physically motivated fixed upper bound of 240 h — an approach more directly consistent with the definition in Eq. 8 of the manuscript. The corrected event-scale τ_rec values range from 2 to 202 h, which is consistent with the 20–80 h range obtained independently from the composite-based analysis in Section 4.2, and no longer shows pathological clustering.
With corrected τ_rec, the event-scale Pearson correlation with τ_RR is r = 0.09 (p = 0.80, n = 10; r = 0.14, p = 0.73 for τ_RR ≥ 0 events). This is considerably weaker than originally reported. We note that with n = 10, statistical significance at p < 0.05 requires |r| > 0.63 — a threshold that constrains any event-scale correlation in this dataset regardless of the underlying physical mechanism. In the revised manuscript, the event-scale analysis will be repositioned as exploratory evidence, with corrected τ_rec and τ_RR values presented in a revised Table 2 alongside explicit uncertainty estimates. Among the tested predictors, the moisture–heat feedback index H retains the strongest association with τ_RR (r = 0.49), consistent with the proposed physical pathway; it will be retained as a mechanistic hypothesis to be evaluated with larger event samples in future work.
Given the corrected analysis and the reviewer’s broader comments, we propose to reorganise the central contribution around three findings that remain robust:
- Composite analysis across 112 events establishes the characteristic thermodynamic structure of the post-rainfall surface response in this dryland setting in Section 4.1. This finding is independent of event-scale τ_rec estimation.
- The τ_rec can be robustly estimated from composite T0cm recovery, exhibits a physically interpretable seasonal structure, and provides a meaningful characterisation of land-surface memory following rainfall.
- The τ_rec → H → τ_RR pathway is presented as a physically motivated hypothesis supported by the composite structure and event-scale case analysis — not a statistically confirmed relationship. This positioning more accurately reflects what ten events can and cannot demonstrate, and defines a testable hypothesis for future observational and modelling work.
We believe this framing is more honest about the evidence and more clearly situates the manuscript within the land–atmosphere coupling literature.
Citation: https://doi.org/10.5194/egusphere-2025-5752-AC1
-
AC1: 'Reply on RC1', Ruolin Li, 25 Mar 2026
reply
Viewed
| HTML | XML | Total | BibTeX | EndNote | |
|---|---|---|---|---|---|
| 217 | 179 | 22 | 418 | 25 | 40 |
- HTML: 217
- PDF: 179
- XML: 22
- Total: 418
- BibTeX: 25
- EndNote: 40
Viewed (geographical distribution)
| Country | # | Views | % |
|---|
| Total: | 0 |
| HTML: | 0 |
| PDF: | 0 |
| XML: | 0 |
- 1
This paper uses a novel observational dataset to probe the relationship between land surface conditions and subsequent precipitation recycling. However, I find that there are some major flaws in the methodology and I am concerned about the authors' interpretation of their results.
First, the methods section is largely based on bullet points and itemized lists, which makes it difficult to follow. This is particularly true for the TRR parameter, which is calculated with ERA5 in a way that is very unclear. The reader is never presented with information on how ERA compares to the station observations. It's not clear how the WLOC and WADV parameters connect to the PADV/PLOC time series. Is it a simple scaling or is there something more complex going on? Does the PADV + PLOC = P measured at the stations?
The results presented in Figs 2 and 3 are interesting, but the plot of Fig. 4 suggests that the seasonal cycle of temperature is mainly responsible for the TREC parameter, which isn't a major surprise because temperature is a strong control on surface evaporation.
Section 4.3.1: claims that there is a consistent temporal structure of the ρ parameter, but this is patently not the case, as (2/10) are negative and (3/10) are greater than 50 hours. So I'm not sure where this claim is coming from.
Section 4.3.2: This claims to be the central finding of this study, but 6/7 of the TREC timescale are equal to 5 (how this is possible is never explained) and the correlation is therefore totally spurious.
Given these errors, the overall results of the study are not convincing to me, making the conclusions drawn from the subsequent analysis not supported by the evidence.