the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Detecting the resilience of soil moisture dynamics to drought periods as function of soil type and climatic region
Abstract. Abrupt changes in climatic conditions and land management can cause permanent shifts in soil hydraulic response to climatic inputs, impacting soil functions and established soil–climate interactions. To quantify the resilience of soil water content dynamics after abrupt changes in environmental conditions, we present a model framework combining a neural network with seasonal trend analysis (STL). Using data from a series of lysimeters from the TERrestrial ENvironmental Observatories (TERENO) – SOILCan lysimeter network, we identified changes in soil water content responses after an extremely hot and dry summer in Germany in 2018. The model incorporates meteorological variables decomposed into seasonal and long-term components along with a categorical indicator of current moisture conditions. It is trained on data from a reference site with stable soil water content response and applied to lysimeters from multiple origins exposed to contrasting climates. By analysing annual residual patterns—particularly mean bias over time—soil water content state dynamics is classified as ‘stable’, ‘resilient’, or ‘changed’, reflecting whether the system maintains, recovers, or diverges from its original state. We found that soils preserve the response function to environmental forcing under typical conditions but exhibit structural change when relocated to new environments, even when soil texture remains constant. The proposed method offers a scalable and non-invasive tool for tracking changes in the response of soil water content to climatic change and provides early indicators of changes in essential soil functions and soil health status.
- Preprint
(1822 KB) - Metadata XML
-
Supplement
(817 KB) - BibTeX
- EndNote
Status: final response (author comments only)
-
RC1: 'Comment on egusphere-2025-5141', Anonymous Referee #1, 19 Nov 2025
-
AC1: 'Reply on RC1', Nedal Aqel, 06 Feb 2026
The comment was uploaded in the form of a supplement: https://egusphere.copernicus.org/preprints/2025/egusphere-2025-5141/egusphere-2025-5141-AC1-supplement.pdf
-
AC1: 'Reply on RC1', Nedal Aqel, 06 Feb 2026
-
EC1: 'Handling editor comment on egusphere-2025-5141', Nunzio Romano, 12 Jan 2026
Dear Authors,
As the handling editor, I would like to offer a few comments that should be discussed, especially for the benefit of readers, potential discussants, and assigned reviewers.
First, let's consider the term "resilience" employed in the current submission. This concept is used in a variety of ways in different contexts. Here, the authors classify soil-water content dynamics as “stable”, “resilient”, and “changed’, depending on changes to the soil-water content profile from its original state (Lines 21-24). I suggest framing this classification more consistently with existing literature, because using the term “resilient” alongside “stable” and “changed” could be misleading. In the hydrological and environmental sciences, resilience is generally accepted as the ability of a system (a lysimeter in this paper) to withstand disturbance, to maintain its functions, whether potential or actual, (which a bit reflects the concepts of resistance or vulnerability), and the system’s recovery rate following disturbance (e.g., Lal, 1997; Pimm, 1984; Gunderson and Holling, 2002). Using the available data in the present submission, it would be interesting to quantify the rate of recovery. Overall, the concepts of resistance, recovery time, and recovery rate should be addressed, I guess.
Anyhow, at line 342, the authors correctly argue that "... it is still possible that the response function may recover in the future ...". However, this important statement is then lost as the text progresses, whereas it should be reiterated at least in the conclusions.
The other comment refers to the authors’ statement on Line 653, I quote: “… This response depends on soil hydraulic properties …”. Of course, the impact that different local climatic conditions, soil being equal, exerts on the dynamics of soil-water content in the lysimeter also depends on changes in soil organic carbon and microbial activity. These latter two aspects seem to be overlooked in this submission. Since the authors claimed to have measured both soil-water content and matric potential, it would be interesting to examine how soil-water retention and hydraulic conductivity functions changed over the investigated period. It is plausible that any changes are mainly attributable to changes in soil organic matter. Another related issue is demonstrating that changes in soil particle arrangement did not occur due to transportation from one site to another.
References cited
Gunderson, L.H., C.S. Holling (eds.), 2002. Panarchy. Washington (DC), 507 pp., Island Press.
Lal, R. 1997. Phil. Trans. R. Soc. Lond. B 352, 997-1010.
Pimm, S.L. 1984. Nature 307, 321-326.
Citation: https://doi.org/10.5194/egusphere-2025-5141-EC1 -
AC3: 'Reply on EC1', Nedal Aqel, 06 Feb 2026
The comment was uploaded in the form of a supplement: https://egusphere.copernicus.org/preprints/2025/egusphere-2025-5141/egusphere-2025-5141-AC3-supplement.pdf
-
AC3: 'Reply on EC1', Nedal Aqel, 06 Feb 2026
-
RC2: 'Comment on egusphere-2025-5141', Anonymous Referee #2, 17 Jan 2026
In this paper the authors use a neural network trained on meteorological variables, a seasonal trend analysis, and categorical definitions of soil conditions to quantify and evaluate soil water response across soil types and climate regions after an extremely hot summer in Germany.
Overall, I found the study and content of this paper compelling as well as the author’s motivations for incorporating increasingly available remotely sensed products into soil stability analyses. However, I would like to see more discussion of the evaluation of the developed model and reasoning behind decisions, such as using Nash–Sutcliffe efficiency.
Further, as the authors identify, performance struggles at the site on which the model was not trained. This paper could be made far more robust through the reverse analysis—training the model in wetter conditions—to see how comparable the results are. Or training the model on both a wet- and dry-regime lysometer simultaneously.
Specific comments
- In Line 75, I’d like to see additional citations to back up the claim of ‘experimental studies’, implying multiple, showing these shifts in soil water content dynamics following extreme events. Or rephrase.
- Clarifying question, beginning at the end of Line 79, are you stating that soil texture remains constant through land-use changes? Or in your modeling approach? Or something else?
- Line 118: While the authors state that predictive modeling providing accurate predictions of soil water dynamics is "not the focus in this study,” I do think that claim requires substantial citations. Alternatively, they could discuss their choice to use the NSE classification here as well as add more information about the pros and cons of the NSE method as this is missing and needed, when the strength of the paper relies on NSE performing as described. Additionally, this sentence (Line 115-119) is quite unwieldy, and it isn’t clear what point the authors are trying to make.
- Figure 4: Considering the study focuses on the summer of 2018, would highlighting that data be evocative rather than the other groupings? It is not clear why the different time periods have been selected as they are not all the same length. Is this arbitrary or is it related to the mean bias results? Regardless, I think some minor explanation is needed. If the purpose of Figure 4 is only to motivate the selection of a particular lysemeter, the spread of values in the contrasting lysimeter is more than enough, without the additional colors).
- Line 252: This may be a naive comment, but is there a reasoning behind the numerical values chosen to represent each categorical condition? For example, why is ‘wet’ 1 rather than 10?
- Relatedly, on what basis are you defining wet, moderate, and dry using the 30th and 70th percentiles? Why not 25th and 75th? Is this arbitrary? Further, have you considered looking at ‘extremely’ wet or dry categories? Considering you call 2018 an ‘extreme’ year throughout the manuscript, I’d be curious to know how incorporating the 95th or 99th percentiles for drying and wetting would be reflected.
- In the mean bias aggregation, did you also look at a seasonal (wet v dry season) mean bias, or only annual? I would be curious if the response function is stable across seasons as well (to ensure there isn’t variability between wet and dry that are cancelling each other out, for instance). This seems to be most evident in Figure 7, where the predicted water content does not capture both the driest and wettest water content periods. In later years the mean bias improves, but the prediction only appears to have a slightly better fit.
- Line 423: How do you conclude that this is a result for coarse soil rather than all/many/some compositions of topsoil? When you say ‘this coarse soil’ are you simply referring to that translocated lysimeter? Or making a statement about coarse soils in general? If the latter, I don’t think this is adequately supported.
- I’d recommend reviewing the language in section 3.2 for clarity, when discussing transposed soil and ‘changed’ climate states. For example in line 447 the use of ‘unchanged climate’ is confusing. Perhaps consider consistent language such as ‘original’ climate, rather than (un)changed to avoid confusion.
- I would caution against describing the differences in climate between the sites as ‘climatic shifts’ (e.g., Line 518). To me this implies a step change in the climate at a given location, not transposing soil between two sites. The timescale of translocating soil is more abrupt than most/all climatic natural and/or anthropogenic climate change-forced shifts.
- Were additional training methods explored for this analysis? (e.g., Is one lysometer enough to train this model? Does fit improve or plateau when trained on multiple lysometers? Would training the model on lysometers from each climate perform better than one lysometer from a single climate?)
- Line 600: what evidence for the physical soil response do you have? I don’t see how this could have been determined from the model or information presented prior. Was the autumn and winter of 2019 comparable to conditions in prior years? Would we necessarily expect a rapid rebound to previous year’s conditions?
- Line 650: I would not classify the top 10cm of soil moisture evaluated in this paper as defining water availability to plants. Most plants have root systems far deeper than 10cm. It seems a stretch to connect these results to agricultural impacts as surface vegetation was not varied or evaluated.
- 673: ‘extreme climatic events’ seems overly broad, as you’re only focusing on a single intense drought, you cannot draw conclusions against other extreme events (e.g. extreme precipitation).
- 675: how are you defining ‘past’ here? Did you evaluate how unique the 2018 drought was in these locations/for these soils?
Technical corrections
- Inconsistent spacing in em-dash usage. There should not be any spaces between words separated by em-dashes (e.g. use is correct on line 22, but not correct on lines 54, 55, 77, 78.
- Line 75: possibly missing an ‘a’ or time scale should be plural.
- Line 101-102: ‘for a set of lysimeters’ I’m not sure ‘for’ is the correct preposition, perhaps ‘with a set of lysimeters’? Unless I’m misunderstanding the function.
- Line 114: for clarity, I’d recommend adding commas around ‘for example’ to set off the clause.
- Line 116: appears to be an extra space between ‘which’ and ‘the’
- Line 160: I would recommend offsetting the place names with commas (e.g., “including both, a wetter site, Selhausen, and a drier site, Bad Lauchstädt, …”). Alternatively, you could use parenthesis.
- Please align tense, particularly in the methods section, where both present and past tense alternate quite frequently.
- Line 215 is stating the same information as Line 130. I recognize the authors are re-introducing the reader to the input variables, but I would recommend restructuring as to not imply that water balance has not been discussed previously. My suggestion would be to either re-introduce all inputs in this section or, better yet, move this discussion up near Line 130. Moving the input variable discussion higher up would also reduce a bit of confusion with the title of subsection 2.2.1 as you aren’t discussing the STD until the second paragraph. In that same vein, the beginning of the 2.2.1 subsection implies that the input variables were selected after data cleaning, which does not seem to be the case (as the input variables were selected for cleaning/alignment).
- Line 298: Should the Adam optimizer be cited?
- Line 322, 670: I’m not sure the 2 hyphen is needed in Nash–Sutcliffe efficiency, nor should the E be capitalized. It may also be more clear to state: ‘we use the Nash-Sutcliffe efficiency coefficient (NSE, see eq. 5) as a general descriptor…”
- Line 444: use ‘which’ or ‘that’ rather than ‘who’
- Lines 519, 520: hyphens used instead of em-dashes.
- Line 615: missing a period.
- Lines 639, 641: hyphens used instead of em-dashes.
- Line 669, 672: the preceding bullets end with a period, whereas these two do not.
- Line 676, 676: include only one does before the (i) and make require plural.
Citation: https://doi.org/10.5194/egusphere-2025-5141-RC2 -
AC2: 'Reply on RC2', Nedal Aqel, 06 Feb 2026
The comment was uploaded in the form of a supplement: https://egusphere.copernicus.org/preprints/2025/egusphere-2025-5141/egusphere-2025-5141-AC2-supplement.pdf
Viewed
| HTML | XML | Total | Supplement | BibTeX | EndNote | |
|---|---|---|---|---|---|---|
| 337 | 79 | 39 | 455 | 33 | 16 | 20 |
- HTML: 337
- PDF: 79
- XML: 39
- Total: 455
- Supplement: 33
- BibTeX: 16
- EndNote: 20
Viewed (geographical distribution)
| Country | # | Views | % |
|---|
| Total: | 0 |
| HTML: | 0 |
| PDF: | 0 |
| XML: | 0 |
- 1
This paper applies statistical modelling techniques (combining a neural network model with seasonal trend analysis) to a comprehensive dataset of soil water contents and pressure potentials measured at 10 cm depth in lysimeters moved to two different locations, in order to identify shifts in hydrological responses to climate forcing.
The shifts in these “in situ” water retention curves (WRC) are intriguing and really quite dramatic (e.g. figures 4, 7 and 8). But I do wonder about the mechanisms and underlying processes. The authors are rather vague about the causes, suggesting that they are due to changes in soil structure tiggered by climate (lines 674-676). I’m not fully convinced about this interpretation, not least because the largest changes seem to occur in the very dry range of the WRC where structure should not play such a large role.
(in this respect, I think the WRC curves should be plotted with matric potential on a log-axis for improved readability. On a linear scale, we can’t really see what is happening close to saturation, which is where most of the structural changes would be expected).
There may be alternative explanations for the observations, including (slowly reversible) swell-shrink behaviour and preferential (non-equilibrium) flow. I would encourage the authors to try to strengthen the discussion and interpretation of the data with respect to the underlying mechanisms, including the above-mentioned processes. Nevertheless, although the responses to climate of apparent WRC observed by the authors “in situ” seem stronger than I would expect (especially in the dry range), I am aware of two previous large-scale (regional-continental) statistical analyses of water retention curves measured in the laboratory that have shown significant impacts of climatic factors on the structural pore space (Hirmas, D. et al. 2018, Nature 561, 100-103; Klöffel, T., et al., 2024. Geoderma, 442, 116772). These studies could be mentioned as they would give support to the authors’ inferences and interpretations.
Specific comments