[go: up one dir, main page]
More Web Proxy on the site http://driver.im/
You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
Sign in to use this feature.

Years

Between: -

Subjects

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Journals

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Article Types

Countries / Regions

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Search Results (936)

Search Parameters:
Keywords = groundwater variations

Order results
Result details
Results per page
Select all
Export citation of selected articles as:
25 pages, 12496 KiB  
Article
Impact of Future Climate Change on Groundwater Storage in China’s Large Granary: A Study Based on LSTM and CMIP6 Models
by Haiqing Wang, Peng Qi, Moran Xu, Yao Wu and Guangxin Zhang
Water 2025, 17(3), 315; https://doi.org/10.3390/w17030315 - 23 Jan 2025
Viewed by 301
Abstract
Northeast China, as a primary grain-producing region, has long drawn attention for its intensive groundwater extraction for irrigation. However, previous studies on the future spatiotemporal changes of groundwater storage (GWS) are lacking. Utilizing the Global Land Data Assimilation System Version 2.2 (GLDAS-2.2), which [...] Read more.
Northeast China, as a primary grain-producing region, has long drawn attention for its intensive groundwater extraction for irrigation. However, previous studies on the future spatiotemporal changes of groundwater storage (GWS) are lacking. Utilizing the Global Land Data Assimilation System Version 2.2 (GLDAS-2.2), which simulates groundwater storage (as Equivalent Water Height) using the Catchment Land Surface Model (CLSM-F2.5) and calibrates it with terrestrial water storage data from the GRACE satellite, we analyzed the spatiotemporal variations of GWS in northeast China and employed a Long Short-Term Memory (LSTM) neural network model to quantify the responses of GWS to future climate change. Maintaining current socio–economic factors and combining climate factors from four scenarios (SSP126, SSP245, SSP370, and SSP585) under the CMIP6 model, we predicted GWS from 2022 to 2100. The results indicate that historically, groundwater storage exhibits a decreasing trend in the south and an increasing trend in the north, with a 44° N latitude boundary. Under the four scenarios, the predicted GWS increments in northeast China are 0.08 ± 0.09 mm/yr in SSP126, 0.11 ± 0.08 mm/yr in SSP245, 0.12 ± 0.09 mm/yr in SSP370, and 0.20 ± 0.07 mm/yr in SSP585. Although overall groundwater storage has slightly increased and the model projections indicate a continued increase, the southern part of the region may not return to past levels and faces water stress risks. This study provides an important reference for the development of sustainable groundwater management strategies. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Spatial distribution of mean GWS (as Equivalent Water Height, representing the spatial pattern of groundwater storage) from February 2003 to December 2022 based on GLDAS-2.2 [<a href="#B29-water-17-00315" class="html-bibr">29</a>].</p>
Full article ">Figure 2
<p>Comparison of GLDAS and CSR GRACE GWSA Data for northeast China, with a NSE of 0.758.</p>
Full article ">Figure 3
<p>Plot of correlation analysis between current and corrected data for climate influences from 1950 to 2014, from left to right, for precipitation, potential evapotranspiration, and temperature.</p>
Full article ">Figure 4
<p>LSTM network structure.</p>
Full article ">Figure 5
<p>Spatial distributions of GWS interannual spacing values in the northeast China, 2004–2022.</p>
Full article ">Figure 6
<p>Spatial trend of the GWS in northeast China, 2003.02–2022.12. (<b>a</b>) The spatial trend of the GWS in northeast China, (<b>b</b>) spatial distribution of <span class="html-italic">p</span>-value corresponding to the trend.</p>
Full article ">Figure 7
<p>Spatial distribution of monthly-scale groundwater storage anomaly (GWSA) in northeast China.</p>
Full article ">Figure 8
<p>Spatial distributions of monthly trends in GWS in northeast China.</p>
Full article ">Figure 9
<p>Spatial distribution of correlation coefficients between GWS and corresponding influencing factors in northeast China (inter-annual data), with no arable land in the Daxinganling region.</p>
Full article ">Figure 10
<p>Spatial distribution of LSTM model fit scores. The data from 2003.02 to 2016.09 were used as the training set, and the data from 2016.10 to 2019.12 were used as the test set. (<b>a</b>,<b>b</b>) were the NSE and RMSE spatial distribution plots of the model test set and the status quo data, respectively.</p>
Full article ">Figure 11
<p>Spatial distribution map of the future GWS. (<b>a</b>) Spatial trend of GWS in different contexts over time and its plotting. (<b>b</b>) Spatial variation of the mean GWS in different contexts over time and its plotting against the mean GWS in the historical period.</p>
Full article ">Figure 12
<p>Temporal changes in GWS in the northeast, with shading indicating the range of fluctuations in GWS under different scenarios. 2004–2019 is a plot of interannual changes in GWS for the historical period, and 2002–2100 is a plot of simulated interannual changes in GWS for the future.</p>
Full article ">Figure 13
<p>Box plots of monthly trends in mean values of future scenarios of groundwater storage in the northeast and changes in GWS. (<b>a</b>,<b>b</b>) Changes in inter-monthly trends for different periods. (<b>c</b>,<b>d</b>) Comparisons of inter-monthly mean values of GWS for different periods with historical GWS. (<b>a</b>,<b>c</b>) The period in 2022–2060; (<b>b</b>,<b>d</b>) the period in 2061–2100.</p>
Full article ">Figure 14
<p>Fitted plots for the four sampled data points. In the figure, the yellow line represents the historical data, the blue line represents the prediction result of the training set, and the green line represents the prediction result of the testing set. (<b>a</b>,<b>b</b>) LSTM model; (<b>c</b>,<b>d</b>) Conv-LSTM model.</p>
Full article ">Figure 15
<p>Schematic diagram of LSTM model.</p>
Full article ">Figure 16
<p>Graph of temporal trends in future climate factors.</p>
Full article ">
19 pages, 2189 KiB  
Article
Study on Alfalfa Water Use Efficiency and Optimal Irrigation Strategy in Agro-Pastoral Ecotone, Northwestern China
by Xiangyang Miao, Guoshuai Wang, Bing Xu, Ruiping Li, Delong Tian, Jie Ren, Zekun Li, Ting Fan, Zisen Zhang and Qiyu Xu
Agronomy 2025, 15(2), 258; https://doi.org/10.3390/agronomy15020258 - 21 Jan 2025
Viewed by 261
Abstract
Agro-pastoral ecotone is an important livestock production area in the north of China, and alfalfa is the main pasture crop in this area. Aiming to address the issues of groundwater overexploitation in the area with water demand, we assessed the consumption pattern, irrigation [...] Read more.
Agro-pastoral ecotone is an important livestock production area in the north of China, and alfalfa is the main pasture crop in this area. Aiming to address the issues of groundwater overexploitation in the area with water demand, we assessed the consumption pattern, irrigation scheduling, and water usage efficiency of alfalfa under subsurface drip irrigation. Alfalfa was used as the research object in this study. A DSSAT model was used to simulate the soil moisture, yield, and other alfalfa grow characteristics during a two-year in situ observation study and provide information on the best irrigation techniques and the water-use efficiency of alfalfa in the agro-pastoral ecotone of Northwestern China. The results showed that the ARE, nRMSE, and R2 values of the alfalfa soil water content, leaf area index, and yield varied between 3.82% and 5.57%, 4.81% and 8.06%, and 0.86 and 0.93, respectively, the accuracy of the calibrated and validated parameters were acceptable, and the model could be applied to this study. The water consumption of alfalfa ranged from 395.6 mm to 421.8 mm during the whole year, and the critical water consumption period was the branching stage and the bud stage. During the branching stage and the bud stage, water consumption was 30–31% and 31–33% of the total water consumption, and the water consumption intensity averaged 2.97–3.04 mm/d and 4.23–4.97 mm/d. The variations of WUE and IWUE were 11.74–14.39 kg·m−3 and 7.12–9.31 kg·m−3. Irrigation increased the water productivity of rain-fed alfalfa by 49.48–64.70% and increased the yield of alfalfa by 17.87–34.72%. With the highest yield as the goal, the recommended irrigation volumes for normal and dry flow years were 200 mm and 240 mm; with the goal of the highest utilization of groundwater resources, the recommended irrigation volumes for normal and dry flow years were 160 mm and 192 mm. The results of this study are expected to provide scientific and technological support for the rational utilization of groundwater and the scientific improvement of alfalfa yields in the agro-pastoral ecotone of Northwestern China. Full article
Show Figures

Figure 1

Figure 1
<p>Location of the study area.</p>
Full article ">Figure 2
<p>Experimental design drawing.</p>
Full article ">Figure 3
<p>Growth stages of different crops.</p>
Full article ">Figure 4
<p>Validation of the soil water content, leaf area, yield of the alfalfa.</p>
Full article ">Figure 5
<p>Soil water dynamics during the depth of 0–60 cm soil profile. Note: SMC stands for Soil Moisture Content. The dot is the measured value, and the curve is the simulation value.</p>
Full article ">Figure 6
<p>Leaf area Change trends of alfalfa. Note: Lal—Leaf area index.</p>
Full article ">Figure 7
<p>Yield change trend of alfalfa.</p>
Full article ">Figure 8
<p>Rainfall frequency curve.</p>
Full article ">Figure 9
<p>Rainfall during alfalfa growth stage from 1991 to 2023. Note: The red line is the trend line.</p>
Full article ">Figure 10
<p>Effects of irrigation schemes on alfalfa yield in different hydrological years.</p>
Full article ">
15 pages, 1727 KiB  
Article
Oxygen Consumption and Carbon Budget in Groundwater-Obligate and Surface-Dwelling Diacyclops Species (Crustacea Copepoda Cyclopoida) Under Temperature Variability
by Tiziana Di Lorenzo, Agostina Tabilio Di Camillo, Sanda Iepure, Diana M. P. Galassi, Nataša Mori and Tatjana Simčič
Environments 2025, 12(1), 32; https://doi.org/10.3390/environments12010032 - 20 Jan 2025
Viewed by 566
Abstract
This study explores the metabolic response and carbon budget of two cyclopoid copepod species, Diacyclops belgicus Kiefer, 1936 (a stygobitic, groundwater-adapted species) and Diacyclops crassicaudis crassicaudis (Sars G.O., 1863) (a stygophilic, predominantly surface-associated species). We measured oxygen consumption rates (OCRs), carbon requirements (CRs), [...] Read more.
This study explores the metabolic response and carbon budget of two cyclopoid copepod species, Diacyclops belgicus Kiefer, 1936 (a stygobitic, groundwater-adapted species) and Diacyclops crassicaudis crassicaudis (Sars G.O., 1863) (a stygophilic, predominantly surface-associated species). We measured oxygen consumption rates (OCRs), carbon requirements (CRs), ingestion (I) rates, and egestion (E) rates at 14 °C and 17 °C, representing current and predicted future conditions in the collection habitats of the two species. Diacyclops belgicus displayed OCRs (28.15 and 18.32 µL O2/mg DW × h at 14 and 17 °C, respectively) and carbon budget (CR: 0.14 and 0.10 µg C/mg × d at 14 and 17 °C) lower than those of D. crassicaudis crassicaudis (OCR: 55.67 and 47.93 µL O2/mg DW × h at 14 and 17 °C; CR: 0.3 and 0.27 µg C/mg × d at 14 and 17 °C). However, D. belgicus exhibited metabolic rates and carbon requirements comparable to those of other epigean species, challenging the assumption that low metabolic rates are universal among stygobitic species. Temperature variations did not significantly affect the metabolic responses and carbon requirements of the two species, suggesting that they may cope with moderate temperature increases. Full article
Show Figures

Figure 1

Figure 1
<p>Phraetobiological net used to collect <span class="html-italic">Diacyclops belgicus</span> (<b>a</b>) and <span class="html-italic">Diacyclops crassicaudis crassicaudis</span> (<b>b</b>).</p>
Full article ">Figure 2
<p>The respirometric setup for oxygen consumption measurements of CV copepodids of <span class="html-italic">Diacyclops belgicus</span> and <span class="html-italic">Diacyclops crassicaudis crassicaudis</span> at 14 °C and 17 °C. One hour after collection, the specimens were kept in darkness for 21 days, transitioning through three media: 100% bore water (black beaker), a 50% bore and standard water mix (green beaker), and 100% standard water (blue beaker). Individual CV copepodids were placed in 80 μL glass wells with oxygen sensor spots, housed in a microplate, and monitored for oxygen levels over 18 h. At the end of the measurements, the specimens were measured, and their body volume was computed based on body dimensions.</p>
Full article ">Figure 3
<p>A heatmap of oxygen consumption rates (OCRs; μL O<sub>2</sub>/mg DW × h) of freshwater cyclopoid species across a temperature gradient: (<b>a</b>) all stages; (<b>b</b>) nauplii; (<b>c</b>) copepodids; (<b>d</b>) adults. Colour intensity indicates OCR magnitude. Mvi: <span class="html-italic">Megacyclops viridis</span> (Jurine, 1820); Mbr: <span class="html-italic">Mesocyclops brasilianus</span> Kiefer, 1933; Ese: <span class="html-italic">Eucyclops serrulatus serrulatus</span> (Fischer, 1851); Eag: <span class="html-italic">Eucyclops agilis agilis</span> (Koch, 1838); Dcr: <span class="html-italic">Diacyclops crassicaudis crassicaudis</span> (Sars G.O., 1863); Dbi: <span class="html-italic">Diacyclops bicuspidatus bicuspidatus</span> (Claus, 1857); Mva: <span class="html-italic">Microcyclops varicans varicans</span> (Sars G.O., 1863); Cvi: <span class="html-italic">Cyclops vicinus vicinus</span> Uljanin, 1875; Dbe: <span class="html-italic">Diacyclops belgicus</span> Kiefer, 1936.</p>
Full article ">
16 pages, 4408 KiB  
Article
Dynamic Heat Transfer Modeling and Validation of Super-Long Flexible Thermosyphons for Shallow Geothermal Applications
by Jianhua Liu, Yanghuiqin Ding, Hao Liu, Liying Zheng, Xiaoyuan Wang and Yuezhao Zhu
Energies 2025, 18(2), 433; https://doi.org/10.3390/en18020433 - 20 Jan 2025
Viewed by 343
Abstract
In comparison to borehole heat exchangers that rely on forced convection, super-long thermosyphons offer a more efficient approach to extracting shallow geothermal energy. This work conducted field tests on a super-long flexible thermosyphon (SFTS) to evaluate its heat transfer characteristics. The tests investigated [...] Read more.
In comparison to borehole heat exchangers that rely on forced convection, super-long thermosyphons offer a more efficient approach to extracting shallow geothermal energy. This work conducted field tests on a super-long flexible thermosyphon (SFTS) to evaluate its heat transfer characteristics. The tests investigated the effects of cooling water temperature and the inclination angle of the condenser on the start-up characteristics and steady-state heat transfer performance. Based on the field test results, the study proposed a dynamic heat transfer modeling method for SFTSs using the equivalent thermal conductivity (ETC) model. Furthermore, a full-scale 3D CFD model for geothermal extraction via SFTS was developed, taking into account weather conditions and groundwater advection. The modeling validation showed that the simulation results aligned well with the temperature and heat transfer power variations observed in the field tests when the empirical coefficient in the ETC model was specified as 2. This work offers a semi-empirical dynamic heat transfer modeling method for geothermal thermosyphons, which can be readily incorporated into the overall simulation of a geothermal system that integrates thermosyphons. Full article
(This article belongs to the Section J1: Heat and Mass Transfer)
Show Figures

Figure 1

Figure 1
<p>(<b>a</b>) Structure of the SFTS, (<b>b</b>) Schematic view of the experimental setup for the super-long thermosyphon, (<b>c</b>) Arrangement of temperature monitoring points.</p>
Full article ">Figure 2
<p>Picture of the experimental setup.</p>
Full article ">Figure 3
<p>Variations of <span class="html-italic">T</span><sub>e</sub>, <span class="html-italic">T</span><sub>c</sub> and <span class="html-italic">Q</span><sub>p</sub> in the test with inclination angle at 0°.</p>
Full article ">Figure 4
<p>Experimental total thermal resistance under varying <span class="html-italic">T</span><sub>cooling</sub> and inclination angles.</p>
Full article ">Figure 5
<p>Variations in thermal parameters from beginning to the steady state with an inclination angle at 0°: (<b>a</b>) outer wall temperature; (<b>b</b>) temperature difference between <span class="html-italic">T</span><sub>e</sub> and <span class="html-italic">T</span><sub>c</sub>; (<b>c</b>) <span class="html-italic">Q</span><sub>p</sub> and (<b>d</b>) <span class="html-italic">Kt</span>.</p>
Full article ">Figure 6
<p>Data fitting for the equivalent thermal conductivity of the SFTS at a steady state.</p>
Full article ">Figure 7
<p>Geometry and the meshing structure of the computational domain.</p>
Full article ">Figure 8
<p>Boundary conditions of the computational domain.</p>
Full article ">Figure 9
<p>Comparison between experimental results and CFD predictions at a steady state: (<b>a</b>) <span class="html-italic">Q</span><sub>p</sub>, and (<b>b</b>) <span class="html-italic">T</span><sub>c</sub>.</p>
Full article ">Figure 10
<p>Comparison between experimental results and CFD predictions in start-up: (<b>a</b>) <span class="html-italic">Q</span><sub>p</sub>, and (<b>b</b>) <span class="html-italic">T</span><sub>c</sub>.</p>
Full article ">
20 pages, 10254 KiB  
Article
Discernible Orientation for Tortuosity During Oxidative Precipitation of Fe(II) in Porous Media: Laboratory Experiment and Micro-CT Imaging
by Wenran Cao, Ekaterina Strounina, Harald Hofmann and Alexander Scheuermann
Minerals 2025, 15(1), 91; https://doi.org/10.3390/min15010091 - 19 Jan 2025
Viewed by 520
Abstract
In the mixing zone, where submarine groundwater carrying ferrous iron [Fe(II)] meets seawater with dissolved oxygen (DO), the oxidative precipitation of Fe(II) occurs at the pore scale (nm~μm), and the resulting Fe precipitation significantly influences the seepage properties at the Darcy scale (cm~m). [...] Read more.
In the mixing zone, where submarine groundwater carrying ferrous iron [Fe(II)] meets seawater with dissolved oxygen (DO), the oxidative precipitation of Fe(II) occurs at the pore scale (nm~μm), and the resulting Fe precipitation significantly influences the seepage properties at the Darcy scale (cm~m). Previous studies have presented a challenge in upscaling fluid dynamics from a small scale to a large scale, thereby constraining our understanding of the spatiotemporal variations in flow paths as porous media evolve. To address this limitation, this study simulated subsurface mixing by injecting Fe(II)-rich freshwater into a DO-rich saltwater flow within a custom-designed syringe packed with glass beads. Micro-computed tomography imaging at the representative elementary volume scale was utilized to track the development of Fe precipitates over time and space. Experimental observations revealed three distinct stages of Fe hydroxides and their effects on the flow dynamics. Initially, hydrous Fe precipitates were characterized by a low density and exhibited mobility, allowing temporarily clogged pathways to intermittently reopen. As precipitation progressed, the Fe precipitates accumulated, forming interparticle bonding structures that redirected the flow to bypass clogged pores and facilitated precipitate flushing near the syringe wall. In the final stage, a notable reduction in the macroscopic capillary number from 3.0 to 0.05 indicated a transition from a viscous- to capillary-dominated flow, which led to the construction of ramified, tortuous flow channels. This study highlights the critical role of high-resolution imaging techniques in bridging the gap between pore-scale and continuum-scale analyses of multiphase flows in hydrogeochemical processes, offering valuable insights into the complex groundwater–seawater mixing. Full article
(This article belongs to the Special Issue Mineral Dissolution and Precipitation in Geologic Porous Media)
Show Figures

Figure 1

Figure 1
<p>(<b>a</b>) Schematic of the experimental setup. Fe(II)-rich freshwater and DO-rich saltwater are co-injected using peristaltic pumps, resulting in a Fe precipitation zone near the inlet. Four cylinders are connected to the syringe via flexible tubes, and their weight readings are recorded using a Raspberry Pi. Outflow was collected in a volumetric container placed on a digital scale. (<b>b</b>) Diagram of the seepage test under a constant head (<span class="html-italic">h</span>) using a Marriot bottle. It is used to determine the sample permeability after dismantling the experiment. Purple arrows represent the flow direction.</p>
Full article ">Figure 2
<p>Segmentation of micro-CT images: (<b>a</b>) 2D raw image of the sample in grayscale, (<b>b</b>) histogram of grayscale values in micro-CT images, and (<b>c</b>) segmented image of the sample (pore space in a specified color).</p>
Full article ">Figure 3
<p>(<b>a</b>) The 3D microstructure showing the pore space (blue) and solid particles (gray). (<b>b</b>) The 3D pore network extraction of the sample (pore space in blue).</p>
Full article ">Figure 4
<p>(<b>a</b>) Spatial distribution of Fe precipitation near the injection point on days 0, 5, 10, 15, 20, and 25; (<b>b</b>) total amount of Fe precipitates in the syringe. (<b>c</b>) Experimental outflow rate over time. Purple triangles represent experimental data; the red dashed line shows the trend. (<b>d</b>) Observed reduction in total permeability. The vertical axis indicates the ratio of permeability <span class="html-italic">k</span> to <span class="html-italic">k</span><sub>0</sub>.</p>
Full article ">Figure 5
<p>(<b>a</b>) Deterministic REV (dREV) analysis. Cubic domains in yellow, purple, and red illustrate increasing volumes with side length increments. (<b>b</b>) Porosity variation across the cubic domain. The dashed line at the 5 mm side length represents the sample sizes used in previous studies on multiphase flow imaging.</p>
Full article ">Figure 6
<p>Schematic diagram of (<b>a</b>) experimental sample setup; (<b>b</b>) early stage—no or temporary pore clogging, with reopened channels in the central sample; (<b>c</b>) middle stage—partial clogging in the central sample redirects flow towards the walls; and (<b>d</b>) final stage—complete clogging in the center, with ramified flow paths along the walls. The purple dashed rectangle indicates the REV used for analysis, the red line represents the tortuous flow path, the blue arrows denote streamlines, and the dashed blue line shows the general flow pattern.</p>
Full article ">Figure 7
<p>Two-dimensional (2D) slices of three-dimensional (3D) images for the transversal cross-section of the sample over time. (<b>a</b>) Initial sample before the experiment. (<b>b</b>) Sample on day 5. (<b>c</b>) Sample on day 15. (<b>d</b>) Sample on day 25. Note that solid particles are depicted in bright gray; saltwater-filled pores are shown in dark blue. The syringe wall is also shown in blue since high-density polythene (HDPE) has a similar density to seawater. Hydrous Fe precipitates are illustrated in cyan, while solid Fe precipitates are shown in orange.</p>
Full article ">Figure 8
<p>Temporal evolution of (<b>a</b>) hydraulic gradient (<span class="html-italic">i</span>) and (<b>b</b>) capillary number (Ca) during Fe precipitation. Purple circles or rectangles represent experimental data; the blue or red dashed line indicates the predicted trend.</p>
Full article ">
27 pages, 4621 KiB  
Article
Analyzing the Impact of Geoenvironmental Factors on the Spatiotemporal Dynamics of Forest Cover via Random Forest
by Hendaf N. Habeeb and Yaseen T. Mustafa
Earth 2025, 6(1), 3; https://doi.org/10.3390/earth6010003 - 14 Jan 2025
Viewed by 573
Abstract
Understanding the dynamic relationships between geoenvironmental factors and forest vegetation cover is crucial for effective forest management and planning. This study investigates the spatiotemporal dynamics of forest cover in the Duhok District in the Kurdistan Region of Iraq over a decade (2013–2023), emphasizing [...] Read more.
Understanding the dynamic relationships between geoenvironmental factors and forest vegetation cover is crucial for effective forest management and planning. This study investigates the spatiotemporal dynamics of forest cover in the Duhok District in the Kurdistan Region of Iraq over a decade (2013–2023), emphasizing the impact of geoenvironmental factors via Random Forest algorithms and Landsat data. This research integrates datasets including fractional vegetation cover (FVC), groundwater levels, climate data, topography, and soil moisture data, offering a comprehensive analysis of the factors influencing forest cover. The results show that in 2013, altitude and rainfall were the primary factors influencing FVC, with areas of higher altitudes and adequate rainfall exhibiting up to 30% denser forest cover. By 2023, soil moisture and groundwater levels had emerged as the dominant factors, with soil moisture levels accounting for 25% of the variation in FVC. This shift underscores the increasing importance of water management strategies to maintain forest health. The Random Forest model demonstrated high predictive accuracy, achieving an R2 value of 0.918 (RMSE of 0.016 and MAE of 0.013) for 2013 and 0.916 (RMSE of 0.018 and MAE of 0.014) for 2023, underscoring the model’s robustness in handling nonlinear ecological processes. This study’s insights are crucial for guiding sustainable forest management practices and assisting decision-makers in formulating strategies for resource management, environmental preservation, and future planning. This study underscores the necessity of adaptive management strategies that consider evolving climatic and hydrological conditions, emphasizing continuous monitoring and advanced technologies to ensure the resilience of forest ecosystems. Full article
Show Figures

Figure 1

Figure 1
<p>Flowchart of the adopted methodology.</p>
Full article ">Figure 2
<p>Study area map: (<b>a</b>) Iraq; (<b>b</b>) Kurdistan Region of Iraq; (<b>c</b>) Duhok District with overlaid Landsat OLI imagery; (<b>d</b>) example plot with sample points; (<b>e</b>) example plot with sample points.</p>
Full article ">Figure 3
<p>Geoenvironmental factors in forest cover analysis for Duhok: FVC ((<b>a</b>): 2013; (<b>b</b>): 2023), SWL ((<b>c</b>): 2013; (<b>d</b>): 2023), T ((<b>e</b>): 2013; (<b>f</b>): 2023), R ((<b>g</b>): 2013; (<b>h</b>): 2023), SMI ((<b>i</b>): 2013; (<b>j</b>): 2023), AL (<b>k</b>), AS (<b>l</b>), and S (<b>m</b>), all of which are used as inputs in the RF model.</p>
Full article ">Figure 4
<p>Heatmaps and Variance Inflation Factor values showing the correlation and multicollinearity of geoenvironmental factors (T: temperature; R: rainfall; AL: altitude; AS: aspect; S: slope; SMI: soil moisture index) for the years (<b>a</b>,<b>c</b>) 2013 and (<b>b</b>,<b>d</b>) 2023. The heatmaps (<b>a</b>,<b>b</b>) display the Pearson correlation coefficients between the factors, whereas the VIF values (<b>c</b>,<b>d</b>) help identify multicollinearity.</p>
Full article ">Figure 5
<p>Correlation and multicollinearity analyses of geoenvironmental factors after excluding temperature (T) for (<b>a</b>,<b>c</b>) 2013 and (<b>b</b>,<b>d</b>) 2023. The heatmaps (<b>a</b>,<b>b</b>) illustrate the correlation structure of the remaining factors, whereas the VIF values (<b>c</b>,<b>d</b>) show reduced multicollinearity after temperature exclusion.</p>
Full article ">Figure 6
<p>Scatter plots comparing the observed and predicted fractional vegetation cover via the Random Forest model for (<b>a</b>) 2013 and (<b>b</b>) 2023. Each dot represents a predicted FVC value compared with its corresponding observed value, with the dashed diagonal line indicating perfect agreement (1:1 correlation). The R<sup>2</sup>, RMSE, and MAE values quantify the model’s predictive accuracy for each year.</p>
Full article ">Figure 7
<p>Bar plots showing the relative importance of geoenvironmental factors influencing fractional vegetation cover as determined by the Random Forest model for (<b>a</b>) 2013 and (<b>b</b>) 2023. The height of each bar represents the contribution of each factor to the model’s predictive accuracy, indicating the dominant factors affecting FVC in the Duhok region for both years.</p>
Full article ">Figure 8
<p>Comparative analysis of geoenvironmental factors and their influence on fractional vegetation cover for 2013 and 2023. This figure highlights how the relationship between environmental factors and forest cover has changed over the decade, with a particular focus on shifts in altitude, precipitation, and soil moisture distributions.</p>
Full article ">
28 pages, 15707 KiB  
Article
Characterizing Groundwater Level Response to Precipitation at Multiple Timescales in the Lubei Plain Region Using Transfer Function Analysis
by Lewei Xu, Huili Gong, Beibei Chen, Chaofan Zhou, Xueting Zhong, Ziyao Ma and Dexin Meng
Remote Sens. 2025, 17(2), 208; https://doi.org/10.3390/rs17020208 - 8 Jan 2025
Viewed by 443
Abstract
Groundwater is essential for ecosystem stability and climate adaptation, with precipitation variations directly affecting groundwater levels (GWLs). Human activities, particularly groundwater exploitation, disrupt the recharge mechanism and the regional water cycle. In this study, we propose a new research framework: On the basis [...] Read more.
Groundwater is essential for ecosystem stability and climate adaptation, with precipitation variations directly affecting groundwater levels (GWLs). Human activities, particularly groundwater exploitation, disrupt the recharge mechanism and the regional water cycle. In this study, we propose a new research framework: On the basis of analyzing the spatiotemporal variability characteristics of precipitation and shallow GWL, we used transfer function analysis (TFA) to quantify the multi-timescale characteristics of precipitation–GWL response under the effects of climate change and human activities. In addition, we evaluated the GWL seasonality and seasonal response while also considering apportionment entropy. We applied this framework to the Lubei Plain (LBP), and the findings indicated the following: (1) Annual precipitation in the LBP decreased from southeast to northwest, with July and August contributing 51.5% of total rainfall; spatial autocorrelation of GWL was high and was influenced by geological conditions and cropland irrigation. (2) The coherence between GWL and precipitation was 0.96 in the high-precipitation areas but was only 0.6 in overexploited areas, and sandy soils enhanced the effective groundwater recharge, with a gain of 1.65 and a lag time of 2.1 months. (3) Over interannual scales, GWL response was driven by precipitation distribution and aquifer characteristics, while shorter timescales (4 months) were significantly affected by human activities, with a longer lag time in overexploited areas, which was nearly 60% longer than areas that were not overexploited. (4) Groundwater exploitation reduced the seasonality of GWL, and irrigation reduced the coherence between GWL and precipitation (0.5), with a gain of approximately 0.5, while a coherence of 0.8 and a gain of 3.5 were observed in the non-irrigation period. This study clarified the multi-timescale characteristics of the precipitation–GWL response, provided a new perspective for regional research on groundwater response issues, and proposed an important basis for the short-term regulation and sustainable development of water resources. Full article
Show Figures

Figure 1

Figure 1
<p>Location maps of the study area. (<b>a</b>) The location of Shandong Province in China. (<b>b</b>) Geographic location of Lubei Plain. (<b>c</b>) Distribution of Lubei Plain. The figure includes the distribution of groundwater monitoring wells and the overexploitation areas in the study area.</p>
Full article ">Figure 2
<p>Overview map of the study area. (<b>a</b>) The distribution map of groundwater monitoring wells; (<b>b</b>) distribution of hydrogeological zoning and faults (inferred) in the LBP, in which “Brackish” represents the brackish water, while “Fresh” represents freshwater in the aquifer; (<b>c</b>) distribution of clay content within the aquifer; (<b>d</b>) distribution of sandy soil content within the aquifer. The clay and sandy soil content data come from National Earth System Science Data Center (<a href="https://www.geodata.cn/" target="_blank">https://www.geodata.cn/</a> (accessed on 11 April 2024)).</p>
Full article ">Figure 3
<p>Precipitation distribution in LBP. (<b>a</b>) Spatial distribution of multi-year average annual precipitation; (<b>b</b>) spatiotemporal distribution of monthly precipitation; (<b>c</b>) verification of multi-source remote sensing precipitation dataset and annual precipitation data of meteorological stations, where the gray column represents the annual precipitation recorded by the meteorological stations, whose positions are shown by the black triangle in (<b>a</b>).</p>
Full article ">Figure 4
<p>Deep and shallow groundwater extraction and groundwater level fluctuations in the regions of NJ (<b>a</b>), WD–YX (<b>b</b>), and BX–GR (<b>c</b>). (SDE is short for shallow groundwater exploitation; DGE is short for deep groundwater exploitation; NJ is short for Ningjin County; WD is short for Wudi County; YX is short for Yangxin County; BX is short for Boxing County; GR is short for Guangrao County).</p>
Full article ">Figure 5
<p>The overall idea and workflow of this paper.</p>
Full article ">Figure 6
<p>Conceptual illustration of transfer function analysis method. (<b>a</b>) Based on the forcing signal <math display="inline"><semantics> <mrow> <mi>X</mi> </mrow> </semantics></math> (i.e., the precipitation series) and the response signal <math display="inline"><semantics> <mrow> <mi>Y</mi> </mrow> </semantics></math> (i.e., the groundwater level series), spectral estimates were obtained using Fast Fourier Transform to extract the major periods and signal components of the two signal resonances. (<b>b</b>) In the time domain, the response relationship between the three main signal components of the forcing signal <math display="inline"><semantics> <mrow> <mi>X</mi> </mrow> </semantics></math> and the response signal <math display="inline"><semantics> <mrow> <mi>Y</mi> </mrow> </semantics></math> is illustrated. In the frequency domain, in order to reflect the frequency-dependent TFA parameters (i.e., coherence, gain, and lag time), four time scales are divided based on the frequency range, which are annual, 6-month, 4-month, and less than 4-month scales, respectively. In particular, the coherence spectrum shows the correlation between the resonances of the two signals; the gain is the ratio of the amplitude of the forcing signal to the amplitude of the response signal at the component with a significant relationship (i.e., <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>A</mi> <mi>m</mi> <mi>p</mi> </mrow> <mrow> <mi>x</mi> </mrow> </msub> <mo>/</mo> <msub> <mrow> <mi>A</mi> <mi>m</mi> <mi>p</mi> </mrow> <mrow> <mi>y</mi> </mrow> </msub> </mrow> </semantics></math>), which represents the response strength of the response signal to a unit forcing signal; and the lag time represents the response time of the response signal to the forcing signal at the corresponding frequency.</p>
Full article ">Figure 7
<p>The semivariogram function of the GWL in the LBP (Gaussian model).</p>
Full article ">Figure 8
<p>Semivariogram function theoretical model and related parameters of groundwater level in the LBP.</p>
Full article ">Figure 9
<p>Spatial distribution of the groundwater level in the LBP (colder colors represent shallower groundwater burial, and warmer colors represent deeper groundwater burial).</p>
Full article ">Figure 10
<p>Spatial distribution of the characteristics of shallow groundwater level response to precipitation in the LBP. The coherence (<b>a</b>), gain (<b>b</b>), and lag time (<b>c</b>) between precipitation and groundwater in the LBP. The blue boxes in (<b>a</b>) represent the two typical zones selected in this study; while the four red circles represent the area a–d.</p>
Full article ">Figure 11
<p>Characteristics of precipitation–groundwater response at monitoring wells within the typical zone. Figures (<b>a</b>–<b>d</b>) correspond to typical areas a–d shown in <a href="#remotesensing-17-00208-f010" class="html-fig">Figure 10</a>b, respectively. The location of the wells is shown in <a href="#remotesensing-17-00208-f002" class="html-fig">Figure 2</a>a.</p>
Full article ">Figure 12
<p>Coherence, gain, and lag time across multiple timescales (annual (<b>a</b>–<b>c</b>), 6-month (<b>d</b>–<b>f</b>), 4-month (<b>g</b>–<b>i</b>), and less than 4-month (<b>j</b>–<b>l</b>) scales) in the LBP.</p>
Full article ">Figure 13
<p>Distribution of precipitation and fluctuation in groundwater level within a year.</p>
Full article ">Figure 14
<p>The distribution of groundwater level seasonality represented by apportionment entropy and land use classification in the LBP. (<b>a</b>) Spatial autocorrelation distribution of groundwater level seasonality. High AE value means low seasonality of GWL; low AE value means significant seasonality. (<b>b</b>) Distribution of land use types and distribution of irrigation points in typical areas in LBP.</p>
Full article ">Figure 15
<p>Comparison of response characteristics between irrigated and non-irrigated periods by transfer function analysis. The location of those wells is shown in <a href="#remotesensing-17-00208-f014" class="html-fig">Figure 14</a>b.</p>
Full article ">
29 pages, 7741 KiB  
Article
Groundwater Storage Estimation in the Saskatchewan River Basin Using GRACE/GRACE-FO Gravimetric Data and Machine Learning
by Mohamed Hamdi, Anas El Alem and Kalifa Goita
Atmosphere 2025, 16(1), 50; https://doi.org/10.3390/atmos16010050 - 6 Jan 2025
Viewed by 726
Abstract
Climate change is having a significant impact on groundwater storage, affecting water resources in many parts of the world. To characterize this impact, remote sensing and machine learning are essential tools to analyze the data accurately and efficiently. This study aims to predicting [...] Read more.
Climate change is having a significant impact on groundwater storage, affecting water resources in many parts of the world. To characterize this impact, remote sensing and machine learning are essential tools to analyze the data accurately and efficiently. This study aims to predicting the variations of groundwater storage (GWS) using GRACE/GRACE-FO and multi-source remote sensing data, combined with machine learning techniques. The approach was applied over the Canadian Prairies region. The study area was classified into three zones of different aquifer potentials (low, medium, and high) using a combination of remote sensing data and the Classification and Regression Trees (CART) approach. The prediction model was developed using a machine-learning approach based on multiple linear regression to estimate GWS variations as a function of various environmental parameters. The results showed that the developed model was able to predict GWS variations with satisfactory accuracy (up to 95% of the explained variance) and good robustness (96% success rate). They also provided a better understanding of the variations in groundwater storage in the Canadian Prairies. Therefore, this work provides a promising method for predicting GWS, which could eventually be applied to other similar environmental conditions. Full article
(This article belongs to the Special Issue The Impact of Climate Change on Water Resources (2nd Edition))
Show Figures

Figure 1

Figure 1
<p>Geographic localization of the Saskatchewan River Basin (SRB).</p>
Full article ">Figure 2
<p>Geological cross section of the SRB [<a href="#B21-atmosphere-16-00050" class="html-bibr">21</a>].</p>
Full article ">Figure 3
<p>General adopted methodology. AHP is for analytic hierarchy process, TWS is for total water storage, GWS is for groundwater storage, GRD is for grid, GPLa is for groundwater piezometric level anomalies. SWAT–MODFLOW is a coupled modeling framework that integrates the Soil and Water Assessment Tool (SWAT) with the MODFLOW groundwater flow model. NARR stands for the North American Regional Reanalysis, a high-resolution dataset produced by the National Centers for Environmental Prediction (NCEP). CART refers to Classification and Regression Trees, a machine-learning algorithm used for both classification and regression tasks.</p>
Full article ">Figure 4
<p>Factors conditioning the groundwater potential used for this study: (<b>a</b>) drainage density, (<b>b</b>) LULC, (<b>c</b>) mean rainfall (2019), (<b>d</b>) slope, and (<b>e</b>) soil map.</p>
Full article ">Figure 5
<p>Groundwater potential map for the SRB: (<b>a</b>) zoom on very low zones and (<b>b</b>) zoom on very high zones.</p>
Full article ">Figure 6
<p>Correlation between the TWS from GRACE/GRACE-FO data and TRMM precipitation data: for zones with (<b>a</b>) low groundwater potential, (<b>b</b>) medium groundwater potential, (<b>c</b>) high groundwater potential, (<b>d</b>) CART classification between TWS and Log (TRMM precipitation) for low groundwater potential zones, (<b>e</b>) Kmeans classification between TWS and Log (TRMM precipitation) for medium groundwater potential zones, and (<b>f</b>) Kmeans classification between TWS and Log (TRMM precipitation) for high groundwater potential zones.</p>
Full article ">Figure 7
<p>(<b>a</b>–<b>c</b>) Histograms showing the frequency distributions of log-transformed precipitation values for Groups 1, 2, and 3. (Columns 2–4) Scatter plots comparing TWS observed versus TWS estimated for Groups 1, 2, and 3, along with performance metrics (R<sup>2</sup>, Nash, Bias, and RMSE). These scatter plots evaluate the model’s performance in estimating TWS across the three hydrological groups.</p>
Full article ">Figure 8
<p>Scatter plots comparing observed total water storage (TWS<sub>Observed</sub>) from GRACE/GRACE-FO with machine learning-based estimated TWS (TWS<sub>Estimated</sub>) for the study area. Colors represent different groups: magenta for <b>Group 1</b> (low TWS and precipitation), blue for <b>Group 2</b> (moderate TWS and precipitation), and yellow for <b>Group 3</b> (high TWS and precipitation). Performance metrics (R<sup>2</sup>, Nash, Bias, and RMSE) are provided for each group ((<b>a</b>) groupe 1, (<b>b</b>) groupe 2 and (<b>c</b>) groupe 3).</p>
Full article ">Figure 9
<p>Variation in the contribution coefficient of ETP and ST NARR (snow depth). REC stands for Relative Explained Contribution. It represents the percentage contribution of each variable to the variability or explanation of the target variable, in this case, groundwater storage (GWS), as determined by the machine learning model. (<b>a</b>,<b>d</b>): Group 1, (<b>b</b>,<b>e</b>): Group 2 and (<b>c</b>,<b>f</b>): Group 3.</p>
Full article ">Figure 10
<p>K-fold cross validation results. The in situ GWS was calculated using groundwater level data from observation wells and specific yield values for the aquifer system. The change in GWS (in cm) (ΔGWS) was derived by multiplying changes in water levels (Δh) with the specific yield (Sy). These values were aggregated spatially to provide a regional estimate of in situ GWS, which was compared with GWS estimates from GRACE/GRACE-FO and machine-learning predictions to assess model performance ((<b>a</b>) groupe 1, (<b>b</b>) groupe 2 and (<b>c</b>) groupe 3).</p>
Full article ">
22 pages, 12094 KiB  
Article
Identification and Analysis on Surface Deformation in the Urban Area of Nanchang Based on PS-InSAR Method
by Mengping Zhang, Jiayi Pan, Peifeng Ma and Hui Lin
Remote Sens. 2025, 17(1), 157; https://doi.org/10.3390/rs17010157 - 5 Jan 2025
Viewed by 541
Abstract
Interferometric Synthetic Aperture Radar (InSAR) technology has emerged as a vital tool for monitoring surface deformation due to its high accuracy and spatial resolution. With the rapid economic development of Nanchang, extensive infrastructure development and construction activities have significantly altered the urban landscape. [...] Read more.
Interferometric Synthetic Aperture Radar (InSAR) technology has emerged as a vital tool for monitoring surface deformation due to its high accuracy and spatial resolution. With the rapid economic development of Nanchang, extensive infrastructure development and construction activities have significantly altered the urban landscape. Underground excavation and groundwater extraction in the region are potential contributors to surface deformation. This study utilized Sentinel-1 satellite data, acquired between September 2018 and May 2023, and applied the Permanent Scatterer Interferometric Synthetic Aperture Radar (PS-InSAR) technique to monitor surface deformation in Nanchang’s urban area. The findings revealed that surface deformation rates in the study area range from −10 mm/a to 6 mm/a, with the majority of regions remaining relatively stable. Approximately 99.9% of the monitored points exhibited deformation rates within −5 mm/a to 5 mm/a. However, four significant subsidence zones were identified along the Gan River and its downstream regions, with a maximum subsidence rate reaching 9.7 mm/a. Historical satellite imagery comparisons indicated that certain subsidence areas are potentially associated with construction activities. Further analysis integrating subsidence data, monthly precipitation, and groundwater depth revealed a negative correlation between surface deformation in Region A and rainfall, with subsidence trends aligning with groundwater level fluctuations. However, such a correlation was not evident in the other three regions. Additionally, water level data from the Xingzi Station of Poyang Lake showed that only Region A’s subsidence trend closely corresponds with water level variations. We conducted a detailed analysis of the spatial distribution of soil types in Nanchang and found that the soil types in areas of surface deformation are primarily Semi-hydromorphic Soils and Anthropogenic Soils. These soils exhibit high compressibility, making them prone to compaction and significantly influencing surface deformation. This study concludes that localized surface deformation in Nanchang is primarily driven by urban construction activities and the compaction of artificial fill soils, while precipitation also has an impact in certain areas. Full article
(This article belongs to the Section Environmental Remote Sensing)
Show Figures

Figure 1

Figure 1
<p>(<b>a</b>) Map of China, highlighting Nanchang’s location. (<b>b</b>) Map of Jiangxi Province, indicating where Nanchang is situated. (<b>c</b>) Map of Nanchang, showing its geographic features.</p>
Full article ">Figure 2
<p>PS-InSAR technical workflow.</p>
Full article ">Figure 3
<p>Temporal and spatial baseline diagram, with the central image being the master image and the others being the slave images.</p>
Full article ">Figure 4
<p>Surface deformation rate map of Nanchang City along the satellite line of sight from 2018 to 2023. Area A is located in Zhongxu Village, Nanchang County; Area B is situated along the shoreline of Xiazhuang Lake in Xinjian District; Area C is located Along Jiangzhong Avenue in Xihu District and at Sunshine Lighting Plaza; Area D is near the Shiqi Resettlement Housing in Nanchang County.</p>
Full article ">Figure 5
<p>Surface deformation rate distribution.</p>
Full article ">Figure 6
<p>Left (<b>a</b>) is the distribution map of Area A, and right (<b>b</b>) is the time-series subsidence line chart of three points in Area A.</p>
Full article ">Figure 7
<p>Left (<b>a</b>) is the distribution map of Area B, and right (<b>b</b>) is the time-series subsidence line chart of three points in Area B.</p>
Full article ">Figure 8
<p>(<b>a</b>–<b>c</b>) show the Google Earth images of Region B. (<b>a</b>) represents 2 May 2014, (<b>b</b>) represents 14 February 2017, and (<b>c</b>) represents 16 November 2019. Label 1 indicates the location of Baojie Machinery Company and Aonong Central China Science and Technology Park, Label 2 represents the location of Huihua Industrial Company.</p>
Full article ">Figure 9
<p>Left (<b>a</b>) is the distribution map of Area C, and right (<b>b</b>) is the time-series subsidence line chart of three points in Area C.</p>
Full article ">Figure 10
<p>(<b>a</b>–<b>c</b>) show the Google Earth images of Region C. (<b>a</b>) represents 16 November 2019, (<b>b</b>) represents 15 November 2020, and (<b>c</b>) represents 3 March 2022. The red polygon indicates the location of the Oupengwan project.</p>
Full article ">Figure 11
<p>Left (<b>a</b>) is the distribution map of Area D, and right (<b>b</b>) is the time-series subsidence line chart of three points in Area D.</p>
Full article ">Figure 12
<p>(<b>a</b>–<b>c</b>) show the Google Earth images of Region D. (<b>a</b>) represents 2 May 2014, (<b>b</b>) represents 27 March 2017, and (<b>c</b>) represents 15 November 2020. Labels 1, 2, and 3 mark the areas of subsidence corresponding to the three points in Region D.</p>
Full article ">Figure 13
<p>Spatial distribution map of soil types in the Nanchang area.</p>
Full article ">Figure 14
<p>(<b>a</b>–<b>d</b>) Relationship between subsidence at various points and monthly cumulative precipitation. (<b>a</b>) subsidence-precipitation relationship in Region A, (<b>b</b>) subsidence-precipitation relationship in Region B, (<b>c</b>) subsidence-precipitation relationship in Region C, (<b>d</b>) subsidence-precipitation relationship in Region D.</p>
Full article ">Figure 15
<p>(<b>a</b>–<b>d</b>) Relationship between precipitation in four regions and the average depth to groundwater in the Poyang Lake Plain. (<b>a</b>) precipitation-groundwater depth relationship in Region A, (<b>b</b>) precipitation-groundwater depth relationship in Region B, (<b>c</b>) precipitation-groundwater depth relationship in Region C, (<b>d</b>) precipitation-groundwater depth relationship in Region D.</p>
Full article ">Figure 16
<p>(<b>a</b>–<b>d</b>) Relationship between water level of Xingzi Station and the subsidence at various points. (<b>a</b>) subsidence-water level relationship in Region A, (<b>b</b>) subsidence- water level relationship in Region B, (<b>c</b>) subsidence- water level relationship in Region C, (<b>d</b>) subsidence- water level relationship in Region D.</p>
Full article ">
21 pages, 10913 KiB  
Article
Impact Assessment of Beach Nourishment on Hot Spring Groundwater on Ibusuki Port Coast
by Nobuyuki Ono, Takatomo Miyake, Kenki Kasamo, Kenji Ishimoto and Toshiyuki Asano
Coasts 2025, 5(1), 1; https://doi.org/10.3390/coasts5010001 - 31 Dec 2024
Viewed by 362
Abstract
This study investigated the thermo-hydrodynamic groundwater environment of a sandy beach where a unique sand bathing method attracts many visitors. The discussed temperatures covered a wide range, from the normal to the boiling temperature of water. We, at first, examined the feasible conditions [...] Read more.
This study investigated the thermo-hydrodynamic groundwater environment of a sandy beach where a unique sand bathing method attracts many visitors. The discussed temperatures covered a wide range, from the normal to the boiling temperature of water. We, at first, examined the feasible conditions for sand bathing and found that the volumetric water content was the crucial factor. Comprehensive field observations were implemented to elucidate two physical quantities: the groundwater flow and the temperature in the sand layer. The latter one was found to be governed by the groundwater level and tidal fluctuations. The characteristics obtained were found to be consistent with the feasible conditions in the landward area. While in the offshore area, the temperature was proved to have suddenly dropped. These results strongly suggest that the underground heat source is distributed in specific spots. A numerical model to describe the groundwater flows and the heat transfer mechanism was developed based on a saturated/unsaturated seepage flow model. The computational results were found to adequately reproduce the observed spatial temperature distribution. The reproduction ability of the model was found to be limited in terms of temporal variations; it was good for the groundwater level, but not for the temperature in the sand. Full article
Show Figures

Figure 1

Figure 1
<p>Location of Ibusuki Port Coast in Japan.</p>
Full article ">Figure 2
<p>Natural spring sand bath (NSSB).</p>
Full article ">Figure 3
<p>An aerial photo of the Ibusuki Port Coast before the implementation of the coastal improvement project (2013/09/17, Geospatial Information Authority of Japan).</p>
Full article ">Figure 4
<p>(<b>a</b>) The topography of the present study site and the regional setting for the numerical analysis (the red, black, and blue lines in (<b>b</b>) delineate the broad domain, intermediate domain, and narrow domain, respectively; the black and red circles in (<b>c</b>) indicate the boring points and boring points including water level measurements, respectively).</p>
Full article ">Figure 5
<p>An image of the completed integrated shore protection project on the Ibusuki Port Coast.</p>
Full article ">Figure 6
<p>Conceptual illustration of hot spring–-groundwater flow in unconfined coastal aquifer.</p>
Full article ">Figure 7
<p>Thermal images of discharging hot spring water at seepage face during ebb tide period.</p>
Full article ">Figure 8
<p>Observation points traversing NSSB area: (<b>a</b>) plane view and (<b>b</b>) side view (number indicates setting points of thermal and water pressure sensors).</p>
Full article ">Figure 9
<p>(<b>a</b>) Installation of measuring devices in observation well. (<b>b</b>,<b>c</b>) Observation poles and thermometer set-up.</p>
Full article ">Figure 10
<p>Geological modeling (<span class="html-italic">n<sub>z</sub></span>: division number of vertical grid; <span class="html-italic">k</span>*: hydraulic conductivity [m/s]).</p>
Full article ">Figure 11
<p>Characteristic curves of volumetric water contents <span class="html-italic">θ</span>; (<b>a</b>) Measured relationship between <span class="html-italic">θ</span> and suction head <span class="html-italic">h</span> (black solid line indicate a fitted curve of VG model); (<b>b</b>) VG model relationship between <span class="html-italic">θ</span> and <span class="html-italic">h</span> (red line shows the relationship <span class="html-italic">θ</span> and relative hydraulic conductivity <span class="html-italic">K<sub>r</sub></span>); (<b>c</b>) Measured relationship between <span class="html-italic">θ</span> and vertical distance from groundwater level <span class="html-italic">z<sub>GWL</sub></span> to sampling point level <span class="html-italic">z<sub>sample</sub></span>.</p>
Full article ">Figure 12
<p>Alongshore distribution of water temperature observations at shoreline: (<b>a</b>) entire measurements (15–17 November 2016) and (<b>b</b>) detailed measurements of area surrounding NSSB site (10–11 May 2017).</p>
Full article ">Figure 13
<p>Cross-sectional distribution of observed TISL: (<b>a</b>) at low tide (26 September 2022 14:00); and (<b>b</b>) at high tide (26 September 2022 18:00) (number indicates setting points of thermal and water pressure sensors, encircled number indicates the points where the sand bathing is performed).</p>
Full article ">Figure 13 Cont.
<p>Cross-sectional distribution of observed TISL: (<b>a</b>) at low tide (26 September 2022 14:00); and (<b>b</b>) at high tide (26 September 2022 18:00) (number indicates setting points of thermal and water pressure sensors, encircled number indicates the points where the sand bathing is performed).</p>
Full article ">Figure 14
<p>Cross-sectional grid system and model setting for heat source depth (white line).</p>
Full article ">Figure 15
<p>Cross-sectional distribution of computed TISL: (<b>a</b>) at low tide (<b>b</b>) at high tide.</p>
Full article ">Figure 16
<p>Comparisons between observed and computed GWL fluctuations: (<b>a</b>) St.-1; (<b>b</b>) St.-3; and (<b>c</b>) St.-5.</p>
Full article ">Figure 17
<p>Comparisons between observed and computed TISL fluctuations: (<b>a</b>) St.-1; (<b>b</b>) St.-2; and (<b>c</b>) St.-3.</p>
Full article ">Figure 18
<p>The profiles of the GSL on the status quo beach and the planned nourishment beaches, and the profiles of the computed GWL for each beach during the low tide condition (the yellow band indicates the feasible region for an NSSB on the status quo beach).</p>
Full article ">Figure 19
<p>The predicted profiles of the TISL after the completion of the beach nourishment (the region enclosed by the red dotted line indicates the feasible region for an NSSB on the status quo beach).</p>
Full article ">
18 pages, 5264 KiB  
Article
Evaluation of Water Replenishment in the Northern Segment of the Yellow River Within the Beijing–Hangzhou Grand Canal, China
by Jingwen Du, Yourang Zan, Qingshun Meng, Hongwei Zhang, Feifei Sun, Chunyan Zhang and Chengpeng Lu
Water 2025, 17(1), 48; https://doi.org/10.3390/w17010048 - 28 Dec 2024
Viewed by 532
Abstract
Recently, frequent water shortages and reductions in water flow have been observed in the northern segment of the Yellow River within the Beijing–Hangzhou Grand Canal. In response, a water replenishment program has been initiated. This study is focused on the section of the [...] Read more.
Recently, frequent water shortages and reductions in water flow have been observed in the northern segment of the Yellow River within the Beijing–Hangzhou Grand Canal. In response, a water replenishment program has been initiated. This study is focused on the section of the Grand Canal north of the Yellow River and uses the GSFLOW model to examine interaction between surface water and groundwater, as well as the effect of water replenishment. The results indicate that, after the water replenishment, the efficiency of water replenishment was highest in the Xiao Canal (64.30%), followed by the Wei Canal (39.09%), the South Canal (12.11%), and the North Canal, which exhibited the lowest efficiency (5.75%). This variation can be attributed to greater water loss with increasing distance from the replenishment source, leading to lower replenishment efficiency. Surface water recharge to groundwater was extended by 32 days, with replenishment effects persisting even after the water supply ceased. The maximum influence distance on either side of the canal reached 5.73 km, with an average impact distance of 1.48 km, resulting in a total affected area of 974.7 km2, accounting for 2.2% of the study area. Water replenishment positively influenced the recovery of groundwater levels along the Grand Canal. Full article
(This article belongs to the Special Issue Advances in Surface Water and Groundwater Simulation in River Basin)
Show Figures

Figure 1

Figure 1
<p>Overview of the northern segment of the Yellow River in the Beijing–Hangzhou Grand Canal.</p>
Full article ">Figure 2
<p>Map of soil types and distribution.</p>
Full article ">Figure 3
<p>(<b>a</b>) Distribution of shallow groundwater types in the study area. (<b>b</b>) Geological cross-section of the North China Plain (east–west direction). Modified from Chen et al. [<a href="#B39-water-17-00048" class="html-bibr">39</a>] and Cao et al. [<a href="#B40-water-17-00048" class="html-bibr">40</a>].</p>
Full article ">Figure 4
<p>Generalization of the study area: (<b>a</b>) Generalization of boundary conditions. (<b>b</b>) Generalized diagram of the river network.</p>
Full article ">Figure 5
<p>Fitting of observed and simulated runoff in the Beijing–Hangzhou Grand Canal using the PRMS model: (<b>a</b>) The calibration period (1 January 2013~31 December 2015) and validation period (1 January 2016~31 December 2017) daily simulation of the model. (<b>b</b>) The calibration period (January 2013~December 2015) and validation period (January 2016~December2017) monthly simulation of the model.</p>
Full article ">Figure 6
<p>Contour fit of groundwater levels in 2017.</p>
Full article ">Figure 7
<p>Process of groundwater storage variation from 2013 to 2017.</p>
Full article ">Figure 8
<p>Hydrograph of surface water–groundwater exchange processes from 2013 to 2017.</p>
Full article ">Figure 9
<p>(<b>a</b>) Process lines of flow rates for each river section during the water replenishment period. (<b>b</b>) Process lines of water replenishment efficiency for each river section.</p>
Full article ">Figure 10
<p>Water exchange relationship between surface water and groundwater in 2022.</p>
Full article ">Figure 11
<p>Scope of influence of water replenishment on the Beijing–Hangzhou Grand Canal in 2022.</p>
Full article ">
14 pages, 2338 KiB  
Article
Effects of Redox Condition on Bacteria-Mediated Hydrochemical Processes and Bacterial Community During Managed Aquifer Recharge
by Mengjie Qin, Haichi You, Weijie Zhang, Longyun Liu, Jinhui Liu and Lu Xia
Sustainability 2025, 17(1), 64; https://doi.org/10.3390/su17010064 - 26 Dec 2024
Viewed by 653
Abstract
During the process of managed aquifer recharge (MAR), when the aerobic surface water is recharged into the reductive aquifer, the redox environment changes along the water pathway. MAR practice can reshape the initial groundwater bacterial community, and further induce variations in the bacteria-mediated [...] Read more.
During the process of managed aquifer recharge (MAR), when the aerobic surface water is recharged into the reductive aquifer, the redox environment changes along the water pathway. MAR practice can reshape the initial groundwater bacterial community, and further induce variations in the bacteria-mediated hydrochemical reactions. In this study, laboratory-scale column experiments were conducted to simulate the processes of aerobic/anaerobic recharge to aquifer. The results showed that the concentration of DO during the aerobic recharge was higher than that of the anaerobic recharge, and ORP showed a similar trend. Active nitrogen transformation was observed during the simulated MAR processes. In the early stages of both the aerobic and anaerobic recharges, nitrate reduction occurred due to denitrification and DNRA. However, in the late stages, nitrification might happen in the aerobic column, and nitrate reduction remained the major process in the anaerobic column. For the bacterial community, Massilia, Ralstonia, Legionella, and Curvibacter predominated under the aerobic recharge. Comparatively, Cedecea, Cupriavidus, and Ralstonia maintained high relative abundances under the anaerobic recharge. Our study provides essential information about the characteristics of bacterial-mediated hydrochemical reactions during the MAR process. The result would enhance understanding of MAR activities and provide valuable insights into the groundwater resources’ sustainable development and management. Full article
Show Figures

Figure 1

Figure 1
<p>Schematic diagram of (<b>a</b>) experimental device and (<b>b</b>) experimental design of percolation experiments. A-MAR and An-MAR represent the simulated aerobic MAR and anaerobic MAR, respectively.</p>
Full article ">Figure 2
<p>Hydrochemical parameters of groundwater during the simulated aerobic MAR (A-MAR) and anaerobic MAR (An-MAR) processes. (<b>a</b>) DO, (<b>b</b>) ORP, (<b>c</b>) DOC, (<b>d</b>) NO<sub>3</sub>-N, (<b>e</b>) NO<sub>2</sub>-N, (<b>f</b>) NH<sub>4</sub>-N.</p>
Full article ">Figure 3
<p>Bacterial community composition on (<b>a</b>) phylum, (<b>b</b>) class, and (<b>c</b>) genus levels, respectively. A1–A3 and An1–An3 represented samples collected from the first layers of the aerobic and anaerobic columns, receptively.</p>
Full article ">Figure 4
<p>Redundancy analysis between groundwater hydrochemical parameters and bacterial communities on (<b>a</b>) phylum and (<b>b</b>) genus levels under different recharge conditions. A1–A3 and An1–An3 represented samples collected from the first layers of the aerobic and anaerobic columns, receptively.</p>
Full article ">Figure 5
<p>Pearson correlation heatmap Circos plot between environmental factors and top 20 dominant genera. Significant differences were indicated as follows: * <span class="html-italic">p</span> &lt; 0.05, ** <span class="html-italic">p</span> &lt; 0.01 and *** <span class="html-italic">p</span> &lt; 0.001.</p>
Full article ">
25 pages, 6044 KiB  
Article
Application of Pumping Tests to Estimate Hydraulic Parameters of Volcanic Aquifers in Lake Tana Basin, Ethiopia
by Fenta Nigate, Alemu Yenehun, Ashebir Sewale Belay, Desale Kidane Asmamaw and Kristine Walraevens
Water 2025, 17(1), 9; https://doi.org/10.3390/w17010009 - 24 Dec 2024
Viewed by 613
Abstract
The purpose of this study was to enhance the understanding and sustainable groundwater management of volcanic aquifer systems by estimating key hydrogeological parameters. The transmissivity of a volcanic aquifer system was estimated using analytical solutions based on 68 constant rate and recovery data [...] Read more.
The purpose of this study was to enhance the understanding and sustainable groundwater management of volcanic aquifer systems by estimating key hydrogeological parameters. The transmissivity of a volcanic aquifer system was estimated using analytical solutions based on 68 constant rate and recovery data sets collected from various sources. A combination of hydro-lithostratigraphy and diagnostic plots was employed to identify the aquifer types and flow conditions, which facilitated model selection. Transmissivity of the confined aquifer was modeled using both Theis and Cooper–Jacob methods, with the Theis residual drawdown solution utilized for estimation. For the unconfined aquifer, the Neuman method was used, and the Hantush/Jacob method was employed for leaky aquifers. The results showed that the transmissivity of the Tertiary basalt varied from 0.38 m2/d to 860 m2/d, while the Quaternary aquifer system ranged from 2.33 m2/d to 1.8 × 104 m2/d, indicating an increase in transmissivity with younger volcanic flows. Specific capacity (SC) was estimated for 74 wells and the values ranged from 0.62 to 5860 m2/d. This wide variation of specific capacity and transmissivity showed significant heterogeneity within the volcanic aquifers. This study introduces the innovative application of derivative diagnostic plots in groundwater research, offering an efficient approach for analyzing and interpreting pumping test data to characterize aquifer systems in various hydrogeologic units. This study focuses on aquifer characterization in hard rock formation, demonstrating methods that can be applied to similar geological environments globally. For the Blue Nile basin in general and for the Lake Tana basin in particular, the study result of aquifer characterization will contribute to exploration, development, and improved groundwater management in the region. Full article
Show Figures

Figure 1

Figure 1
<p>Location map of the study area showing its geographical extent and major boundaries (elevation: above mean sea level).</p>
Full article ">Figure 2
<p>(<b>a</b>) Geological map showcasing the prominent towns within the study area (modified from Nigate et al. [<a href="#B28-water-17-00009" class="html-bibr">28</a>]). (<b>b</b>) Hydrogeological map showing borehole locations in the Tana basin (the locations are represented by number only; the prefix TBH- is omitted here to avoid too much text in the map).</p>
Full article ">Figure 3
<p>Lithologic descriptions and well design for (<b>a</b>,<b>b</b>) Quaternary basaltic aquifer systems and (<b>c</b>) the Tertiary basaltic aquifer.</p>
Full article ">Figure 4
<p>Diagnostic plots and interpretation of constant pumping rate regarding Tertiary volcanic aquifer for (<b>a</b>) TBH-23, (<b>b</b>) TBH-29, (<b>c</b>) TBH-39), and (<b>d</b>) TBH-60 (log-log); Quaternary volcanic aquifer for (<b>e</b>) TBH-67 (log-log) and (<b>f</b>) TBH-84 (semi-log); and Tertiary aquifer with different flow types for (<b>g</b>) TBH-52 (log-log) and (<b>h</b>) TBH-56 (log-log). In the time–drawdown graphs, drawdown is represented by (▯) and log-derivative is represented by (+).</p>
Full article ">Figure 5
<p>Time–drawdown graphs of single well-pumping test data obtained from (<b>a</b>) TBH-29, (<b>b</b>) TBH-39, (<b>c</b>) TBH-23, and (<b>d</b>) TBH-60 for constant rate test in Tertiary hard rock. Theis [<a href="#B30-water-17-00009" class="html-bibr">30</a>] model-type curve is represented by the solid line and drawdown data are represented by (▯).</p>
Full article ">Figure 6
<p>Transmissivity estimation from recovery data for boreholes: (<b>a</b>) TBH-23, (<b>b</b>) TBH-29, (<b>c</b>) TBH-39, and (<b>d</b>) TBH-60 using the Theis residual drawdown method.</p>
Full article ">Figure 7
<p>Time–drawdown graphs of single well-pumping test data obtained from (<b>a</b>) TBH-67 and (<b>b</b>) TBH-84 for constant rate test in the Quaternary hard rock formation. Theis [<a href="#B30-water-17-00009" class="html-bibr">30</a>] model-type curve is represented by a solid line and drawdown data are represented by (▯).</p>
Full article ">Figure 8
<p>Transmissivity estimation from recovery data for boreholes: (<b>a</b>) TBH-67 and (<b>b</b>) TBH-84. (<b>b</b>) using the Theis residual drawdown method [<a href="#B42-water-17-00009" class="html-bibr">42</a>].</p>
Full article ">Figure 9
<p>Spatial distribution of (<b>a</b>) transmissivity and (<b>b</b>) specific capacity in the basin (see <a href="#water-17-00009-f002" class="html-fig">Figure 2</a>b for the legend explaining the hydrogeologic unit and geological structures).</p>
Full article ">Figure 10
<p>(<b>a</b>) frequency distribution of log transmissivity and (<b>b</b>) frequency distribution of log-specific capacity in a log-normal diagram; (<b>c</b>) arithmetic plot of transmissivity versus specific capacity; (<b>d</b>) log-log plot of transmissivity versus specific capacity of the observed data.</p>
Full article ">Figure 11
<p>Scale factor vs. power coefficient relationship for different aquifer units analyzed by different researchers around the globe and for the Lake Tana basin volcanic aquifer.</p>
Full article ">
19 pages, 3765 KiB  
Article
Integrating Satellite Observations and Hydrological Models to Unravel Large TROPOMI Methane Emissions in South Sudan Wetlands
by Yousef A. Y. Albuhaisi, Ype van der Velde, Sudhanshu Pandey and Sander Houweling
Remote Sens. 2024, 16(24), 4744; https://doi.org/10.3390/rs16244744 - 19 Dec 2024
Viewed by 529
Abstract
This study presents a comprehensive investigation of Methane (CH4) emissions in the wetlands of South Sudan, employing an integrated approach that combines TROPOMI satellite data, river altimetry, and hydrological model outputs. TROPOMI data show a strong increase in CH4 concentrations [...] Read more.
This study presents a comprehensive investigation of Methane (CH4) emissions in the wetlands of South Sudan, employing an integrated approach that combines TROPOMI satellite data, river altimetry, and hydrological model outputs. TROPOMI data show a strong increase in CH4 concentrations over the Sudd wetlands from 2018 to 2022. We quantify CH4 emissions using these data. We find a twofold emission increase from 2018 to 2019 (9.2 ± 2.4 Tg yr−1) to 2020 to 2022 (16.3 ± 3.3 Tg yr−1). River altimetry data analysis elucidates the interconnected dynamics of river systems and CH4 emissions. We identify correlations and temporal alignments across South Sudan wetlands catchments. Our findings indicate a clear signature of ENSO driving the wetland dynamics and CH4 emissions in the Sudd by altering precipitation patterns, hydrology, and temperature, leading to variations in anaerobic conditions conducive to CH4 production. Significant correlations are found between CH4 emissions and PCR-GLOBWB-simulated soil moisture dynamics, groundwater recharge, and surface water parameters within specific catchments, underscoring the importance of these parameters on the catchment scale. Lagged correlations were found between hydrological parameters and CH4 emissions, particularly with PCR-GLOBWB-simulated capillary rise. These correlations shed light on the temporal dynamics of this poorly studied and quantified source of CH4. Our findings contribute to the current knowledge of wetland CH4 emissions and highlight the urgency of addressing the complex interplay between hydrology and carbon dynamics in these ecosystems that play a critical role in the global CH4 budget. Full article
(This article belongs to the Section Environmental Remote Sensing)
Show Figures

Figure 1

Figure 1
<p>(<b>a</b>) South Sudan study region. Data layers from Stamen Terrain-USA/OSM and OpenStreetMap Humanitarian Data Model. (<b>b</b>) South Sudan river streams within South Sudan borders. The red dots are Hydroweb river altimetry measurement points used in the study.</p>
Full article ">Figure 2
<p>The TROPOMI XCH<sub>4</sub> annual average enhancement over the South Sudan Wetlands Region (SSWR) from 2018 to 2022 at 0.1° × 0.1° resolution. The SSWR, indicated by the black rectangle (latitude 5–10° N and longitude 28–34.5° E), encompasses the area of interest.</p>
Full article ">Figure 3
<p>Monthly CH<sub>4</sub> emissions for the South Sudan Wetlands Region (SSWR) derived from TROPOMI data from 2018 to 2022. The Y-axis represents CH<sub>4</sub> emissions anomalies in Tg CH<sub>4</sub>. The anomalies in CH<sub>4</sub> emissions are calculated as the deviation from the multi-year monthly mean, standardized by dividing the difference by the standard deviation (SD) of the observed values for that particular month across all years. Positive values indicate higher-than-average emissions, while negative values represent lower-than-average emissions.</p>
Full article ">Figure 4
<p>The TROPOMI-inferred CH<sub>4</sub> enhancements (black), river altimetry (pink), GRACE equivalent water thickness (green), and ENSO index (blue and red) for the South Sudan Wetland Region (SSWR) highlighting the interconnected dynamics. The anomalies in CH<sub>4</sub> emissions are calculated as the deviation from the multi-year monthly mean, standardized by dividing the difference by the standard deviation (SD) of the observed values for that particular month across all years. Positive values indicate higher-than-average emissions, while negative values represent lower-than-average emissions.</p>
Full article ">Figure 5
<p>(<b>a</b>) South Sudan Wetlands Region (SSWR) catchments, and (<b>b</b>) TROPOMI CH<sub>4</sub> emission in comparison to river altimetry measurements per catchment.</p>
Full article ">Figure 6
<p>Correlation between PCR-GLOBWB sub-surface parameters anomalies and TROPOMI CH₄ emission anomalies.</p>
Full article ">Figure 7
<p>Correlation between PCR-GLOBWB surface parameters anomalies and TROPOMI CH₄ emission anomalies.</p>
Full article ">Figure 8
<p>Correlation between PCR-GLOBWB river parameters anomalies and TROPOMI CH₄ emission anomalies.</p>
Full article ">Figure 9
<p>Comparison of normalized TROPOMI CH₄ emissions from the SSWR (black) alongside normalized water levels of Lake Victoria (green). The bars indicate the phases of the El Niño Southern Oscillation (ENSO), with El Niño events represented in blue and La Niña events in red, highlighting potential correlations between climatic phenomena and CH₄ emissions.</p>
Full article ">
18 pages, 3901 KiB  
Article
Assessing the Effects of Wheat Planting on Groundwater Under Climate Change: A Quantitative Adaptive Sliding Window Detection Strategy
by Lingling Fan, Shi Chen, Lang Xia, Yan Zha and Peng Yang
Atmosphere 2024, 15(12), 1501; https://doi.org/10.3390/atmos15121501 - 16 Dec 2024
Viewed by 521
Abstract
Climate change has led to changes in precipitation patterns, exacerbating the overextraction of groundwater for wheat irrigation. Although many studies have examined the effects of wheat cultivation on groundwater storage (GWS), few studies have directly assessed the effects of wheat planting on GWS. [...] Read more.
Climate change has led to changes in precipitation patterns, exacerbating the overextraction of groundwater for wheat irrigation. Although many studies have examined the effects of wheat cultivation on groundwater storage (GWS), few studies have directly assessed the effects of wheat planting on GWS. We proposed a wheat subsiding effect detection (WSED) strategy using time-series remote sensing image to assess the effect of wheat area on GWS across China. The subsiding magnitude of the WSED is calculated as the GWS difference between the wheat area and adjacent nonwheat area in the self-adaptive moving window (the size and position of the sliding window can be automatically adjusted based on the characteristics of the data at the central pixel location). The effects of the wheat area on groundwater storage differ greatly among the change types of wheat area and planting regionalization, characterized by the strong subsiding effect in the wheat stable area, gain area, and Huanghuaihai zone (HWW, the most important wheat-producing region in China mainly includes the provinces and municipalities of Beijing, Tianjin, Henan, Hebei, Shandong, Anhui, and Jiangsu). Nearly 80% of the wheat area in the stable and gain regions had lower groundwater depth than nonwheat areas with significant differences (p < 0.05), resulting in a clear declining groundwater trend of approximately −1 cm/year. This study provides quantitative evidence for the effects of wheat planting on GWS regarding agricultural production and climate change adaptations. Full article
(This article belongs to the Special Issue Observation of Climate Change and Cropland with Satellite Data)
Show Figures

Figure 1

Figure 1
<p>Spatial distribution of Chinese wheat in 2019, with the background map indicating the topography. We divided China into four wheat planting zones (Northern Spring Wheat zone, NSW; Northwest Winter–Spring Wheat zone, NWSW; Huanghuaihai Winter Wheat zone, HWW; Southern Winter Wheat zone, SWW) based on the Chinese wheat planting regionalization obtained from a previous study [<a href="#B23-atmosphere-15-01501" class="html-bibr">23</a>].</p>
Full article ">Figure 2
<p>The schematic representation of the detection method for the effects of wheat planting area on GWS. (<b>a</b>) The mechanism framework of the effects of spatial nearby wheat planting area on GWS within a searching window. (<b>b</b>) The schematic representation of the relationship between the wheat planting area and GWS. The “W”, “N”, and “G” represent wheat pixel, nonwheat pixel, and groundwater storage.</p>
Full article ">Figure 3
<p>Flowchart of the wheat subsiding effect detection (WSED) method. The default parameter values were set as L<sub>0</sub> = 150 km, E = 50 m, and L<sub>max</sub> = 300 km.</p>
Full article ">Figure 4
<p>Temporal (<b>a</b>–<b>c</b>) and spatial (<b>d</b>) changes in wheat area in China for 2000–2019. The line charts (<b>a</b>,<b>b</b>) represent the annual wheat area in China, HWW, NSW, NWSW, and SWW, and the dotted lines represent their fitted trends for the study period. The bar charts represent the annual wheat area in terms of different wheat production regions, and the lines present their annual fractions of the total wheat area (<b>c</b>). The right map shows the changing trend of the wheat area from 2000 to 2019, indicated by the slope of the linear fitted lines of the wheat area in each grid (5 × 5 km) (<b>d</b>).</p>
Full article ">Figure 5
<p>Stacked frequency distributions (<b>a</b>,<b>d</b>,<b>g</b>) and difference (mean ± SD) between wheat and nonwheat areas (<b>b</b>,<b>e</b>,<b>h</b>) and the trends of GWS from 2004 to 2016 (<b>c</b>,<b>f</b>,<b>i</b>) in different change regions of wheat area. A negative (positive) ΔGWS indicates a subsiding (rising) effect. Significance levels of the difference (mean ± SD) in the slope of ΔGWS between wheat and nonwheat areas have been marked (* <span class="html-italic">p</span> &lt; 0.05, ** <span class="html-italic">p</span> &lt; 0.01, *** <span class="html-italic">p</span> &lt; 0.001). NS refers to not significant (<b>b</b>,<b>e</b>,<b>h</b>). The dark red dots and line indicate the annual value and fitted line of the annual GWS, and the dark and light red shaded areas indicate the confidence interval and prediction interval of the GWS from 2004 to 2016 (<b>c</b>,<b>f</b>,<b>i</b>).</p>
Full article ">Figure 6
<p>The averages of annual ΔGWS in the comparison experiments with different values of L<sub>0</sub> (100 km and 200 km, (<b>a</b>,<b>b</b>)), N (3000 and 4000, (<b>c</b>,<b>d</b>)), and E (40 m and 60 m, (<b>e</b>,<b>f</b>)) for a given year (taking the wheat map in 2019 for experiments). The insets show the correlation between the ΔGWS of the comparison experiments and that of the initial experiment (L<sub>0</sub> = 150 km, N = 3500, E = 50 m). The red dotted line is the 1:1 line, and the red asterisk indicates the significance level (* <span class="html-italic">p</span> &lt; 0.05, ** <span class="html-italic">p</span> &lt; 0.01).</p>
Full article ">Figure 7
<p>The averages of annual ΔGWS in the comparison experiments with different values of L<sub>0</sub> (100 km and 200 km, (<b>a</b>,<b>b</b>)), N (3000 and 4000, (<b>c</b>,<b>d</b>)), and E (40 m and 60 m, (<b>e</b>,<b>f</b>)) for the wheat stable area. The insets show the correlation between the ΔGWS of the comparison experiments and that of the initial experiment (L<sub>0</sub> = 150 km, N = 3500, E = 50 m). The red dotted line is the 1:1 line, and the red asterisk indicates the significance level (** <span class="html-italic">p</span> &lt; 0.01).</p>
Full article ">Figure 8
<p>The monthly average groundwater storage (GWS) from April 2002 to March 2021 in China. Blank pixels represent outliers that are replaced by null values.</p>
Full article ">Figure 9
<p>Comparison of the GWS in the wheat and nonwheat areas for the wheat stable area in China. (<b>a</b>) The monthly stacked frequency distributions of the difference (ΔGWS) between wheat and nonwheat areas. (<b>b</b>) Comparison of the GWS values between the wheat and nonwheat areas from January to April. (<b>c</b>) Comparison of the GWS values between the wheat and nonwheat areas from May to August. (<b>d</b>) Comparison of the GWS values between the wheat and nonwheat areas from September to December. Asterisks indicate a significant level (*** <span class="html-italic">p</span> &lt; 0.001, ** <span class="html-italic">p</span> &lt; 0.01, * <span class="html-italic">p</span> &lt; 0.05).</p>
Full article ">
Back to TopTop