- College of Tobacco Science, Henan Agricultural University, Zhengzhou, China
Introduction: Henan Province plays a pivotal role in ensuring China’s food security through its extensive winter wheat–summer maize rotation system. However, excessive groundwater extraction has become a critical constraint on sustainable agricultural development.
Methods: This study quantitatively analyzed the spatiotemporal dynamics of crop water use and groundwater depletion across Henan from 1961 to 2020 using a water balance framework that integrates crop evapotranspiration (ETc), effective precipitation (Re), irrigation water requirement (Iwr), and net groundwater consumption (NGWC).
Results and Discussion: Results showed distinct and opposite trends in ETc between the two crops: winter wheat exhibited increasing ETc under climate warming, while summer maize experienced a decline driven by improved water-use efficiency. Effective precipitation decreased during key growth stages, leading to increased irrigation demand for winter wheat but a reduction for summer maize. Winter wheat remained the dominant driver of groundwater depletion, consuming three to four times more groundwater annually than maize. However, since 2010, maize groundwater use has risen sharply, intensifying water stress in the central–northern plains. To alleviate overexploitation, a dual strategy combining winter wheat fallowing and summer maize water conservation is proposed. These findings provide a scientific foundation for optimizing cropping patterns and improving water resource sustainability in major grain-producing regions of North China.
1 Introduction
Henan Province, located in the core of the North China Plain, is one of China’s most important grain production regions, supporting national food security through its extensive winter wheat–summer maize rotation system (Liu et al., 2013). The high multiple-cropping index and intensive use of fertilizer and irrigation have substantially increased land productivity, but they also rely heavily on groundwater abstraction to buffer the pronounced mismatch between crop water demand and uneven precipitation in this semi-humid to semi-arid environment (Wang and Shi, 2024b). Over the past several decades, continuous overdraft has led to widespread groundwater level decline, land subsidence, river flow reduction and degradation of associated ecosystems, making groundwater depletion a key constraint on the long-term sustainability of agricultural production in Henan and across the broader North China Plain (Du et al., 2024; Famiglietti, 2014).
In response to these challenges, considerable research in the field of agricultural water management has focused on quantifying crop water use, irrigation demand and groundwater stress in intensively cultivated regions. At the field and regional scales, water balance models, remote-sensing–driven approaches and water footprint assessments have been used to estimate reference evapotranspiration, crop evapotranspiration (ETc), effective precipitation (Re), and irrigation water requirement (Iwr) for major cereal crops, and to evaluate their implications for blue-water consumption and water scarcity (D’odorico et al., 2020; Li et al., 2024b; Wang et al., 2023; Li et al., 2023). In the North China Plain, numerous studies have reported increasing evaporative demand for winter wheat under climate warming, contrasted with declining or stabilized ETc for summer maize due to the combined effects of the evaporation paradox and improved varieties and management (Ti et al., 2018). Parallel work using hydrological and groundwater-flow models has revealed large-scale depletion of aquifers driven by irrigation, and has highlighted the need for more efficient cropping patterns and irrigation strategies in major grain-producing regions (Cao et al., 2013; Gleeson et al., 2020; Tan et al., 2023).
Despite this progress, several important knowledge gaps remain for Henan’s winter wheat–summer maize rotation system. First, most existing studies have treated winter wheat and summer maize separately, or have only analyzed a single crop or limited time periods, making it difficult to capture the coupled and potentially compensating effects of the two crops on groundwater over multi-decadal scales (Mo et al., 2017; Liu et al., 2017). Second, while ETc, Re or Iwr have been widely examined, relatively few studies have integrated these components into a unified water balance framework that explicitly links crop water consumption to net groundwater consumption (NGWC) and its spatial redistribution within double-cropping systems (Zhao et al., 2021). Third, the seasonal shift of groundwater stress within the rotation—particularly the evolving contribution of summer maize in the context of policy-driven yield improvement and climate variability—has not been systematically quantified using long time series and stage-resolved indicators. Finally, although winter wheat fallowing has been proposed as a promising measure to alleviate groundwater depletion in high-risk zones, there is still a lack of empirical evidence that simultaneously evaluates the hydrological consequences of winter wheat fallowing and summer maize water conservation within a consistent analytical framework (Cao et al., 2020; Kang et al., 2017; Pereira et al., 2002).
Addressing these gaps is crucial for designing rotation-specific water-saving strategies that reconcile groundwater protection with food production targets. By jointly considering both crops over several decades, and by explicitly quantifying irrigation water demand and NGWC at different growth stages and spatial scales, it is possible to reveal whether groundwater stress is merely concentrated in winter wheat, or is progressively transferring and amplifying within the summer maize season under changing climatic and management conditions. Such insights are essential for refining fallow policies, optimizing crop structure and improving irrigation scheduling in Henan Province and similar double-cropping regions worldwide.
Therefore, this study focuses on the winter wheat–summer maize rotation system in Henan Province from 1961 to 2020 with three specific objectives: (1) to quantify the long-term spatiotemporal patterns of ETc, Re or Iwr for winter wheat and summer maize under the rotation system; (2) to estimate stage-based and decadal NGWC and annual net groundwater consumption (ANGWC) for both crops and to reveal how groundwater depletion has evolved across space and time; and (3) to assess the implications of alternative cropping and irrigation strategies, including winter wheat fallowing and summer maize water conservation, for relieving groundwater stress and supporting sustainable agricultural water management in this major grain-producing region.
2 Materials and methods
2.1 Study area and regional characteristics
Henan Province, located in central-eastern China on the North China Plain, is one of the nation’s most important agricultural regions, with extensive wheat–maize double-cropping systems forming the backbone of local food production. The province spans 18 prefecture-level cities and covers diverse climatic and soil conditions, ranging from warm temperate semi-humid zones in the north to subtropical monsoon regions in the south (Figure 1).
Figure 1. Geographical distribution of meteorological stations and altitude in Henan Province, China.
Annual precipitation generally increases from northwest to southeast, while temperature and evaporation display inverse gradients. Soils are dominated by alluvial, brown, and yellow-brown types, providing a favorable environment for both summer maize and winter wheat cultivation. This climatic and edaphic heterogeneity underpins Henan’s role as a major grain production base, while also shaping regional differences in irrigation demand and groundwater dependence.
2.2 Data sources
This study utilized a 60-year dataset (1961–2020) integrating meteorological, agricultural, and crop-specific data. Meteorological data were obtained from 17 stations of the China Meteorological Network (http://data.cma.cn), encompassing variables such as daily temperature, wind speed, sunshine duration, precipitation, and humidity. Crop data, including growth-period information and crop coefficients, were derived from the China Meteorological Yearbook and the Food and Agriculture Organization of the United Nations (FAO) (China Meteorological Yearbook Editorial Department, 1961-2020; Allen et al., 1998). County-level sown area data were extracted from municipal statistical yearbooks to estimate spatial distributions of crop cultivation (China Agricultural Yearbook Editorial Committee, 1961-2020).
2.3 Calculation of crop water requirements
Reference crop evapotranspiration (ET0) was calculated using the FAO-recommended Penman–Monteith model (Allen et al., 1998) (Equation 1):
where ETc is reference evapotranspiration (mm·day-1); Re is net radiation (MJ·m-2·day-1), calculated based on solar radiation derived from sunshine duration data; G is soil heat flux, considered negligible (G ≈ 0) for daily calculation steps; T is air temperature (°C); u2 is wind speed at 2 m (m·s-1); es is the saturation vapor pressure (kPa), calculated from air temperature; ea is the actual vapor pressure (kPa), derived from relative humidity data; Δ is the slope of the vapor-pressure curve (kPa·°C-1); and γ is the psychrometric constant (kPa·°C-1).
Crop evapotranspiration (ETc) was derived as Equation 2:
where Kc is the crop coefficient.
Effective precipitation (Re) represents the portion of rainfall retained within the root zone sufficient for crop use (Kang and Cai, 1996) (Equation 3):
where Pi is total precipitation (mm) and αi is the effective utilization coefficient:
The fixed coefficient method employed for estimating effective precipitation (Re), while operationally straightforward and widely applied in large-scale, long-term hydrological assessments, possesses inherent limitations, as shown in Equation 4. This simplified approach does not explicitly account for the dynamic influences of soil hydraulic properties (e.g., texture, structure, water-holding capacity), topography, or rainfall intensity and duration on the actual partitioning of rainfall into infiltration, runoff, and deep percolation. Consequently, this method may introduce systematic biases in Re estimation across different pedo-climatic zones. For instance, in regions dominated by sandy soils with high infiltration rates and low water-holding capacity, the fixed coefficients might overestimate the fraction of precipitation retained within the root zone available for crop use. Conversely, in areas with clayey soils characterized by higher water retention but potentially slower initial infiltration, the method could underestimate Re, particularly for low-intensity rainfall events. These potential inaccuracies in Re would propagate directly into the subsequent calculations of irrigation water requirement (Iwr) and NGWC, likely leading to an overestimation of Iwr and NGWC in sandy soil areas and an underestimation in clayey soil regions. Despite these limitations, the method provides a consistent and computationally efficient framework for analyzing long-term spatiotemporal trends across the extensive and diverse study area. The findings and conclusions of this study should therefore be interpreted with an understanding of this methodological simplification, which highlights a valuable direction for future research to incorporate more process-based soil water balance models.
Irrigation water requirement (Iwr) was then estimated as Equation 5:
Positive Iwr indicates irrigation demand exceeding rainfall; negative values imply rainfall sufficiency.
2.4 Inverse distance weighted interpolation
Spatial interpolation of meteorological-station data to surface scale was achieved using the inverse distance weighting method (Rolim et al., 2011):
where Z is the interpolated value, Zi is the observed value at station i, di is the distance to the interpolation point, n is the number of stations (typically 3), and k = 2 is the distance-decay parameter (Equation 6).
2.5 Estimation of net groundwater consumption
Net groundwater consumption quantifies the average groundwater volume consumed per unit area by crops over their growth cycle, as shown in Equation 7:
where Re,i,j is effective precipitation, ETc,i,j is crop water requirement, m is the number of years, and n is the number of growth stages. Positive values denote groundwater recharge; negative values indicate net depletion.
Annual net groundwater consumption (ANGWC) represents the total volume of groundwater consumed by a specific crop within a given region over a defined time period. To obtain stable and representative results at the county scale, and to minimize the potential influence of extreme climatic events (such as severe droughts or floods) on crop water use in a single year, a decadal mean approach was adopted for the calculation of ANGWC.
In this method, the decadal average NGWC of summer maize and winter wheat was first calculated for each ten-year period (e.g., 1981–1990, 1991–2000, 2001–2010, and 2011–2020). The decadal mean NGWC was then multiplied by the sown area of the corresponding mid-point year of each decade (e.g., 1985, 1995, 2005, and 2015) to estimate the representative ANGWC at the county level. This approach can be expressed as Equation 8:
where denotes the annual net groundwater consumption of crop c (winter wheat or summer maize) in county k; represents the ten-year average net groundwater consumption (mm) for crop c; and is the cultivated area (ha) of the mid-point year of the corresponding decade.
By linking decadal mean groundwater consumption intensity with the representative cropping structure of each period, this approach effectively reduces the interference of short-term climate variability and interannual fluctuations. It provides a more robust and realistic representation of the long-term dynamics of groundwater consumption in Henan’s wheat–maize rotation system, ensuring that the derived county-level ANGWC values reflect both hydrological and agronomic conditions in a balanced and scientifically consistent manner.
2.6 Mann–Kendall Trend test and Sen’s slope estimation
2.6.1 Mann–Kendall trend test
The Mann–Kendall (M-K) test is a non-parametric method for detecting monotonic trends in time-series data (Kendall, 1957) (Equation 9):
With Equation 10
The variance is given by Equation 11:
and the standardized statistic Zs is computed as Equation 12:
At a 5% significance level (α = 0.05), trends are significant when |Zs| > 1.96.
2.6.2 Sen’s slope estimation
Sen’s method estimates the magnitude of a trend by computing the median of all pairwise slopes (Sen, 1968) (Equation 13):
where xj and xz are data at time points j and z (j > z), and b represents the annual change rate. The approach is robust to outliers and widely used to quantify changes in climatic and hydrological variables.
2.7 Statistical analysis
All meteorological and agricultural data were processed using RStudio and Microsoft Excel 2021. Spatial analyses and visualizations were performed in ArcGIS, Origin 2021, and Excel 2019, enabling clear representation of the spatiotemporal evolution of evapotranspiration, irrigation demand, and groundwater consumption for both crops.
3 Results
3.1 Spatiotemporal distribution of crop evapotranspiration
During 1961–2020, crop evapotranspiration (ETc) for summer maize and winter wheat in Henan Province exhibited pronounced spatiotemporal differentiation and strikingly opposite temporal trends (Figures 2, 3). Summer maize showed a consistent downward trend in ETc across all growth stages, with the most substantial decline occurring during the mid-growth stage (M3ETc), where Sen’s slope reached −1.236 mm·year-1. Significant decreases were also observed during the development (M2ETc) and early growth (M1ETc) stages, suggesting a gradual reduction in atmospheric evaporative demand or improved water-use efficiency.
Figure 2. Spatiotemporal variation of ETc (mm) for summer maize across different growth stages in Henan Province, 1961–2020.
Figure 3. Spatiotemporal variation of ETc (mm) for winter wheat across different growth stages in Henan Province, 1961–2020.
In contrast, winter wheat exhibited a distinct increase in ETc, particularly during the development (W2ETc) and mid-growth (W3ETc) stages. In some central and northern regions, this increase was statistically significant (|Z| > 1.96), indicating regionally robust upward trends. This opposite pattern of ETc evolution between the two crops highlights their differing responses to climate change and irrigation practices.
Spatially, both crops displayed peak water requirements during their mid-growth stages, but the spatial patterns and magnitudes differed markedly. Summer maize ETc reached a maximum of 280 mm, with high-value zones concentrated in the central and northern plains, including Zhengzhou, Xuchang, and Luohe, areas characterized by relatively high solar radiation and strong evapotranspiration potential. Winter wheat reached its peak ETc during the development stage (max: 301 mm), with high values distributed across the central, western and northern plains, where temperature and radiation levels rise sharply during spring. The broader range of ETc variation in winter wheat indicates its greater sensitivity to local meteorological fluctuations.
Notably, the increasing ETc trend observed for winter wheat coincided with a declining trend in effective precipitation, indicating a growing mismatch between water demand and natural supply. This imbalance has intensified irrigation dependence during the critical stages of tillering and heading. Conversely, the declining ETc in summer maize, combined with a moderate rise in early-season precipitation, reflects a relatively improved water supply–demand relationship. These opposite patterns underscore a seasonal shift in water stress within the rotation system.
3.2 Spatiotemporal distribution of effective precipitation
Effective precipitation (Re) during the growth stages of summer maize and winter wheat exhibited complex spatiotemporal dynamics across Henan Province (Figures 4, 5). For summer maize, Re displayed distinct phased characteristics. During the early growth stage (M1Re), a significant upward trend was observed across most of the province, suggesting an improvement in moisture conditions during the seedling establishment phase. However, as crop development progressed, spatial divergence in Re became increasingly evident. By the mid-growth (M3Re) and maturity (M4Re) stages, a clear and consistent downward trend dominated the spatial pattern, with the most significant reduction during the mid-growth stage (Sen’s slope = −1.012 mm·year-1). This indicates that rainfall during the main water consumption period of maize has gradually become less effective.
Figure 4. Spatiotemporal variation of effective precipitation (Re, mm) for summer maize across different growth stages in Henan Province, 1961–2020.
Figure 5. Spatiotemporal variation of effective precipitation (Re, mm) for winter wheat across different growth stages in Henan Province, 1961–2020.
For winter wheat, Re variations showed a generally declining trend across all growth stages. The maturity stage (W4Re) exhibited the most pronounced and consistent decrease (Sen’s slope = −0.965 mm·year-1), followed by the mid-growth stage (W3Re). During the early (W1Re) and development (W2Re) stages, spatial differentiation was evident, with localized increases in southern areas but declines dominating the overall pattern.
Spatially, summer maize Re peaked during the development stage at 169 mm, followed by the mid-growth stage at 149 mm, exhibiting strong spatial heterogeneity throughout. High-value zones were mainly distributed in southern Henan. In contrast, winter wheat Re values were generally lower, with a maximum of 105 mm during the development stage and a minimum of 50 mm during the early stage. The relatively low and uneven precipitation during winter wheat’s growing period underscores its vulnerability to moisture deficits and explains its heavy reliance on groundwater irrigation.
Overall, these spatiotemporal variations in effective precipitation constitute a primary driver of the contrasting irrigation requirements and groundwater consumption patterns between summer maize and winter wheat.
3.3 Spatiotemporal distribution of irrigation water requirement
The irrigation water requirements (Iwr) for summer maize and winter wheat exhibited markedly divergent trends between 1961 and 2020 (Figures 6, 7). For summer maize, Iwr displayed a persistent downward trend across all growth stages, consistent with the overall reduction in crop evapotranspiration and the moderate rise in effective precipitation during early growth stages. The most substantial decline occurred during the mid-growth stage (M3Iwr), with a Sen’s slope of −1.273 mm·year-1. Similar downward tendencies were also evident during the early (M1Iwr), development (M2Iwr), and maturity (M4Iwr) stages. Spatially, summer maize Iwr was highly concentrated during the mid-growth stage, reaching a maximum of 252 mm, while remaining relatively low during the other stages. High-value zones were clustered mainly in the central and northern plains.
Figure 6. Spatiotemporal variation of irrigation water requirement (Iwr, mm) for summer maize across different growth stages in Henan Province, 1961–2020.
Figure 7. Spatiotemporal variation of irrigation water requirement (Iwr, mm) for winter wheat across different growth stages in Henan Province, 1961–2020.
In contrast, winter wheat exhibited more complex and spatially heterogeneous patterns. During the early growth stage (W1Iwr), most regions showed a declining trend, reflecting generally favorable moisture conditions during early winter. However, at the critical development stage (W2Iwr), the pattern diverged sharply: while several regions, especially in central Henan, experienced a significant increase (|Z| > 1.96), others exhibited notable decreases. This spatial polarization reflects growing disparities in water supply–demand relationships across different subregions. During the mid-growth stage (W3Iwr), the overall trend shifted towards increasing irrigation demand, implying a sustained intensification of water stress during this critical phase.
Spatially, winter wheat Iwr peaked during the development stage, reaching 298 mm, and exhibited pronounced spatial heterogeneity throughout all stages. High-value zones were concentrated in the central and northern plains, which correspond closely to areas of severe groundwater overexploitation.
Overall, the contrasting trajectories of Iwr between summer maize and winter wheat reveal a seasonal redistribution of irrigation pressure: maize irrigation demand is gradually easing, while wheat irrigation requirements continue to intensify.
3.4 Spatiotemporal distribution of net groundwater consumption
3.4.1 Stage-based net groundwater consumption
Winter wheat and summer maize exhibited fundamentally different effects on groundwater dynamics and depletion trajectories (Figures 8, 9). For winter wheat, NGWC values remained negative throughout all growth stages, indicating a persistent net consumption of groundwater across its entire life cycle. Interannual analyses revealed clear downward trends in NGWC during the mid-growth stage (W3NGWC) and maturity stage (W4NGWC), with Sen’s slopes reaching −1.144 mm·year-1 and −1.167 mm·year-1, respectively. These negative slopes reflect a progressive intensification of groundwater withdrawal during these critical phenological phases. The developmental stage (W2NGWC) displayed spatially mixed trends, with both positive and negative shifts, suggesting heterogeneous management practices and localized hydrological differences across regions.
Figure 8. Spatiotemporal variation of net groundwater consumption (NGWC, mm) for summer maize across different growth stages in Henan Province, 1961–2020.
Figure 9. Spatiotemporal variation of net groundwater consumption (NGWC, mm) for winter wheat across different growth stages in Henan Province, 1961–2020.
In contrast, summer maize exhibited distinct phased behavior in NGWC evolution. During the initial stage (M1NGWC) and development stages (M2NGWC), most NGWC values were positive, indicating that effective precipitation not only satisfied crop water demands but also contributed net recharge to groundwater. However, by the mid-growth stage (M3NGWC), all NGWC values turned negative, identifying this phase as the critical period of net groundwater depletion for summer maize. Particularly noteworthy is the statistically significant upward trend in NGWC during the early growth stage (|Z| > 1.96), signifying a long-term transformation from “net recharge” to “net depletion”.
Spatially, winter wheat exhibited its greatest groundwater depletion during the development stage, reaching a minimum NGWC of −279.98 mm, with high-consumption zones concentrated in the central–northern plains. Similarly, summer maize showed high depletion levels in these same areas but displayed positive (recharge) values in the south, where precipitation and soil water retention are higher.
These results indicate that winter wheat is a net groundwater consumer throughout its entire growth period, whereas summer maize consumption is concentrated during the mid-growth stage. This difference in phased consumption patterns suggests that implementing precision irrigation management tailored to different crops and growth stages holds significant potential for optimizing water resource utilization.
3.4.2 Decadal dynamics of annual net groundwater consumption
To minimize the influence of extreme climatic years on groundwater assessment, decadal mean NGWC values were used to estimate county-level ANGWC (Figures 10, 11). The results reveal contrasting decadal trajectories between the two crops. The ANGWC of summer maize remained relatively stable before 2010, fluctuating around –19.49 T (1 T = 106 m3), but surged sharply to –45.36 T during 2011–2020—an increase exceeding 130%. This acceleration corresponds to both expanded cultivation and rising evapotranspiration deficits, indicating growing groundwater dependence in the post-2010 period. In contrast, winter wheat consistently dominated total groundwater depletion. Its ANGWC increased steadily from –131.25 T during 1991–2000 to –162.02 T in 2011–2020, remaining three to four times higher than that of summer maize.
Figure 10. Spatiotemporal variation in annual net groundwater consumption (ANGWC, 106 m3) of summer maize at the county scale in Henan Province, 1981–2020.
Figure 11. Spatiotemporal variation in annual net groundwater consumption (ANGWC, 106 m3) of winter wheat at the county scale in Henan Province, 1981–2020.
Spatially, the central-northern plains emerged as persistent depletion hotspots throughout all decades, with overlapping high-consumption zones for both crops. These areas correspond to the region’s densest double-cropping belts and most intensive irrigation use, leading to long-term overextraction. The southern counties, characterized by higher precipitation and smaller irrigation reliance, displayed relatively balanced or even slightly positive ANGWC, indicating episodic groundwater recharge.
Overall, the decadal-scale evolution of ANGWC demonstrates that while winter wheat remains the primary driver of chronic groundwater depletion, the rapid post-2010 intensification in summer maize consumption has aggravated the total extraction pressure on aquifers. The superimposed effect of both crops in the central-northern region has transformed this zone into a critical area of groundwater stress.
4 Discussion
This study employs a water balance–based analytical framework to systematically elucidate the mechanisms underlying groundwater depletion within Henan Province’s winter wheat–summer maize rotation system over the past six decades. The results not only reaffirm the widespread agricultural water scarcity across the North China Plain but also reveal—through the opposite temporal trajectories of the two crops—a dynamic process of groundwater stress transfer and amplification within the rotation system itself (Piao et al., 2010). Furthermore, the post-2010 surge in groundwater depletion, particularly for summer maize, underscores the need for an integrated analysis that combines climatic, agronomic, and increasingly important policy and management drivers. These findings deepen the understanding of how climate variability and agronomic adaptation jointly shape groundwater consumption patterns in intensively cultivated double-cropping systems.
4.1 Drivers of changes in crop water requirements
This study demonstrates that the evapotranspiration (ETc) of winter wheat increased significantly during critical growth stages, whereas ETc of summer maize generally decreased. These divergent trajectories reflect the combined influences of climate change and agronomic adaptation (Zhang et al., 2013).
For winter wheat, the rising ETc is consistent with findings across the North China Plain, where increased air temperature and solar radiation during spring and early summer have elevated reference evapotranspiration (ET0) and thereby enhanced ETc (Mo et al., 2005). The spatial significance of this trend confirms that climate warming has systematically raised the evaporative potential of winter wheat fields (Tan et al., 2022).
Conversely, the decline in ETc for summer maize reflects both climatic moderation and technological improvement. The “evaporation paradox”, characterized by declining wind speed and reduced sunshine duration, has moderated ET0 (Wang et al., 2007). Concurrently, advances in maize breeding have introduced drought-tolerant, short-duration, and high-efficiency cultivars with optimized canopy structures and improved transpiration efficiency (Chaves et al., 2003; Cao et al., 2015). These improvements, together with adaptive irrigation and soil management, have enhanced field-level water-use efficiency, offsetting the potential ETc increase from climatic warming (Zhang et al., 2017, 2024).
4.2 Drivers of effective precipitation variation
The evolution of effective precipitation (Re) in Henan Province reflects changing precipitation regimes and surface hydrological conditions. For summer maize, the early-growth-stage increase in Re is linked to interannual fluctuations in the East Asian summer monsoon (Ding et al., 2008). However, the subsequent decline during mid- and late growth stages indicates a structural change, where increasing rainfall concentration and extreme events promote runoff and percolation losses, reducing the fraction of rainfall effectively available to crops (Cao et al., 2013; Sun et al., 2022).
For winter wheat, the widespread decline in Re across all phenological stages corresponds to the long-term decrease in winter–spring precipitation and snowfall observed throughout North China. Spatial analysis further indicates that soil texture and structure exert a strong control over the spatial heterogeneity of Re. The brown and cinnamon soils dominating the central–northern plains of Henan, although possessing favorable tillage properties, have relatively low water-holding capacity compared with the yellow-brown and paddy soils prevalent in the southern region. Under persistently limited rainfall conditions, these coarse-textured soils exhibit a higher risk of desiccation and reduced infiltration efficiency, thereby amplifying the effects of meteorological drought and diminishing the effective utilization of precipitation (Kang et al., 2017).
This spatial coupling between soil hydro-physical properties and precipitation effectiveness explains why regions dominated by brown and cinnamon soils coincide with zones of severe water deficit in winter wheat production (Zhao et al., 2024). The results indicate that precipitation pattern shifts and soil-mediated infiltration limitations together govern the spatial distribution of water supply–demand imbalances across Henan’s wheat–maize rotation system (Wang et al., 2024a).
4.3 Drivers of irrigation water requirement variations
Variations in irrigation water requirement (Iwr), the crucial link between surface processes and groundwater extraction, reflect the combined influences of climate forcing, agronomic adaptation, and water management capacity.
For summer maize, the consistent decline in Iwr across all growth stages results from the joint effects of reduced ETc and increased early-season Re. This finding aligns with reports that water-saving technologies substantially mitigate maize irrigation requirements in the North China Plain (Kang et al., 2017).
Conversely, winter wheat exhibits increasing Iwr during its critical growth phases, accompanied by pronounced spatial heterogeneity. This upward trend arises from the widening gap between water demand and supply—driven by higher ETc and lower Re—and is further complicated by regional disparities in irrigation infrastructure and groundwater accessibility. In areas with historically abundant groundwater and well-developed irrigation systems, farmers tend to increase irrigation frequency to stabilize yields under rising evaporative demand, leading to sharp increases in Iwr. In contrast, in water-scarce regions where groundwater levels have already declined or strict pumping restrictions are enforced, farmers often reduce irrigation intensity, resulting in a decline in Iwr but also potential yield losses (Zhai et al., 2005). This spatial differentiation in irrigation behavior reveals that groundwater overexploitation is not solely a hydrological problem but is deeply rooted in socio-economic disparities and uneven adaptation capacity. Consequently, one-size-fits-all macro-level policies may prove insufficient to address region-specific irrigation challenges (Liu et al., 2018a, 2018).
4.4 Impact of different cropping patterns on groundwater consumption
This study quantitatively confirms winter wheat as the dominant contributor to groundwater depletion in Henan’s double-cropping system, with an ANGWC three to four times greater than that of summer maize (Liu et al., 2017). Seasonal fallowing of winter wheat in high-risk zones thus represents a direct and effective measure to reduce abstraction and facilitate aquifer recovery (Bian et al., 2024).
However, the sharp post-2010 increase in summer maize NGWC signals a critical shift. This transformation is driven by: (1) Policy-driven expansion and adoption of high-yielding, water-intensive cultivars under grain security initiatives (Yang et al., 2013); (2) Irrigation management shifts, including increased frequency and practices like no-till direct seeding that can disrupt water demand-rainfall synchrony (Scanlon et al., 2012); and (3) Climatic factors like more frequent late-summer droughts. Consequently, winter wheat fallowing policies implemented without parallel regulation of summer maize irrigation may see groundwater savings offset by increased maize water use (Li et al., 2024a).
To ensure sustainable groundwater management, a regionally differentiated dual strategy is proposed. This entails targeted winter wheat fallowing coupled with strict irrigation quotas and water-saving technologies for summer maize in the critically overexploited central-northern plains, while optimizing irrigation scheduling in the wetter south. Effective policy frameworks must integrate fallow subsidies with quota-based water allocation (Pereira et al., 2002; Cao et al., 2020) and promote efficient irrigation technologies and cultivars (Duan et al., 2018). Addressing socio-economic challenges, such as farm income losses and adoption costs, through tailored incentives is paramount for balancing groundwater sustainability with food security and livelihoods (Gleick, 2003).
4.5 Research limitations
This study retains certain limitations that point to future research directions. The spatial interpolation using the Inverse Distance Weighting (IDW) method, while effective for capturing spatial patterns, may introduce uncertainty in sparsely sampled western mountainous and peripheral areas, as it cannot quantify prediction uncertainty like geostatistical methods (e.g., Kriging) (Vicente-Serrano et al., 2003). Furthermore, the fixed coefficient method for effective precipitation estimation, though operationally simple and suitable for long-term analysis, may introduce systematic biases by not fully accounting for the dynamic impacts of soil hydraulic properties and rainfall intensity on infiltration and runoff (Kim and Stricker, 1996).
5 Conclusions
This study investigated the spatiotemporal dynamics of crop water use and groundwater depletion in Henan Province’s winter wheat–summer maize rotation system from 1961 to 2020 using a water balance framework integrating ETc, Re, Iwr, and NGWC. Results reveal strong temporal and spatial contrasts between the two crops. Winter wheat ETc increased due to climate warming during spring and early summer, while summer maize ETc decreased, reflecting improved water-use efficiency and adaptive management. Effective precipitation declined in key growth stages, with rainfall becoming more concentrated and infiltration efficiency reduced, particularly in coarse-textured soils of northern Henan. Consequently, Irrigation water demand decreased for maize but increased for winter wheat, indicating widening water supply–demand disparities. Groundwater depletion is dominated by winter wheat, whose ANGWC remains three to four times that of maize. However, the post-2010 increase in summer maize groundwater use signals a shift toward dual-season water stress, especially in the central–northern plains. To ensure sustainable groundwater management, this study recommends a dual strategy integrating winter wheat fallowing and summer maize water conservation. Coordinated policies—combining targeted fallow subsidies, irrigation quotas, and high-efficiency irrigation technologies—are essential for balancing agricultural production with groundwater sustainability. These findings provide a scientific basis for designing adaptive water management strategies in China’s major grain-producing regions.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Author contributions
ZZ: Data curation, Writing – original draft. JW: Data curation, Writing – original draft. JN: Conceptualization, Writing – original draft. JT: Funding acquisition, Project administration, Writing – review & editing.
Funding
The author(s) declared that financial support was received for this work and/or its publication. This research was funded by the Key Scientific Research Project of Henan Higher Education Institutions, grant number 23A210021.
Conflict of interest
The authors declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declared that generative AI was used in the creation of this manuscript. During the revision process, the authors used DeepSeek-V3.2 for polishing the English grammar. The authors have reviewed and edited the output and take full responsibility for the content of this publication.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Glossary
ETc: Crop Evapotranspiration
Re: Effective Precipitation
Iwr: Irrigation Water Requirement
NGWC: Net Groundwater Consumption
ANGWC: Annual Net Groundwater Consumption
ET0: Reference Crop Evapotranspiration
FAO: Food and Agriculture Organization of the United Nations
IDW: Inverse Distance Weighting
M-K Test: Mann-Kendall Test
M1: Initial Stage (Summer Maize)
M2: Development Stage (Summer Maize)
M3: Mid-Season Stage (Summer Maize)
M4: Late Season Stage (Summer Maize)
W1: Initial Stage (Winter Wheat)
W2: Development Stage (Winter Wheat)
W3: Mid-Season Stage (Winter Wheat)
W4: Late Season Stage (Winter Wheat)
M1ETc: ETc at the Initial Stage of Summer Maize
W2Re: Re at the Development Stage of Winter Wheat
M3Iwr: Iwr at the Mid-Season Stage of Summer Maize
W4NGWC: NGWC at the Late Season Stage of Winter Wheat
M1981-1990: ANGWC for the period 1981-1990 (Summer Maize)
M1991-2000: ANGWC for the period 1991-2000 (Summer Maize)
M2001-2010: ANGWC for the period 2001-2010 (Summer Maize)
M2011-2020: ANGWC for the period 2011-2020 (Summer Maize)
W1981-1990: ANGWC for the period 1981-1990 (Winter Wheat)
W1991-2000: ANGWC for the period 1991-2000 (Winter Wheat)
W2001-2010: ANGWC for the period 2001-2010 (Winter Wheat)
W2011-2020: ANGWC for the period 2011-2020 (Winter Wheat)
References
Allen, R. G., Pereira, L. S., Raes, D., and Smith, M. (1998). Crop evapotranspiration-Guidelines for computing crop water requirements-FAO Irrigation and drainage paper 56 Vol. 300 (Rome: Fao), D05109.
Bian, J., Hu, X., Shi, L., Min, L., Zhang, Y., Shen, Y., et al. (2024). Evapotranspiration partitioning by integrating eddy covariance, micro-lysimeter and unmanned aerial vehicle observations: A case study in the North China Plain. Agr. Water Manage. 295, 108735. doi: 10.1016/j.agwat.2024.108735
Cao, X. C., Wang, Y. B., Wu, P. T., and Zhao, X. N. (2015). Water productivity evaluation for grain crops in irrigated regions of China. Ecol. Indic. 55, 107–117. doi: 10.1016/j.ecolind.2015.03.003
Cao, X. C., Zeng, W., Wu, M. Y., Guo, X. P., and Wang, W. G. (2020). Hybrid analytical framework for regional agricultural water resource utilization and efficiency evaluation. Agr. Water Manage. 231, 106027. doi: 10.1016/j.agwat.2020.106027
Cao, G. L., Zheng, C. M., Scanlon, B. R., Liu, J., and Li, W. P. (2013). Use of flow modeling to assess sustainability of groundwater resources in the North China Plain. Water Resour. Res. 49, 159–175. doi: 10.1029/2012wr011899
Chaves, M. M., Maroco, J. P., and Pereira, J. S. (2003). Understanding plant responses to drought - from genes to the whole plant. Funct. Plant Biol. 30, 239–264. doi: 10.1071/fp02076
China Agricultural Yearbook Editorial Committee (1961-2020). China agricultural yearbook. (Beijing: China Agriculture Press).
China Meteorological Yearbook Editorial Department (1961-2020). Meteorological yearbook of China. (Beijing: China Meteorological Press).
D’odorico, P., Chiarelli, D. D., Rosa, L., Bini, A., Zilberman, D., and Rulli, M. C. (2020). The global value of water in agriculture. Proc. Natl. Acad. Sci. U.S.A. 117, 21985–21993. doi: 10.1073/pnas.2005835117
Ding, Y. H., Wang, Z. Y., and Sun, Y. (2008). Inter-decadal variation of the summer precipitation in East China and its association with decreasing Asian summer monsoon. Part I: Observed evidences. Int. J. Climatol. 28, 1139–1161. doi: 10.1002/joc.1615
Du, J., Laghari, Y., Wei, Y.-C., Wu, L., He, A.-L., Liu, G.-Y., et al. (2024). Groundwater depletion and degradation in the North China Plain: Challenges and mitigation options. Water 16, 354. doi: 10.3390/w16020354
Duan, J. Z., Wu, Y. P., Zhou, Y., Ren, X. X., Shao, Y. H., Feng, W., et al. (2018). Approach to higher wheat yield in the huang-huai plain: improving post-anthesis productivity to increase harvest index. Front. Plant Sci. 9. doi: 10.3389/fpls.2018.01457
Famiglietti, J. S. (2014). The global groundwater crisis. Nat. Clim. Change 4, 945–948. doi: 10.1038/nclimate2425
Gleeson, T., Cuthbert, M., Ferguson, G., and Perrone, D. (2020). Global groundwater sustainability, resources, and systems in the anthropocene. Annu. Rev. Earth Planet. Sci. 48, 431–463. doi: 10.1146/annurev-earth-071719-055251
Gleick, P. H. (2003). Global freshwater resources: Soft-path solutions for the 21st century. Science 302, 1524–1528. doi: 10.1126/science.1089967
Kang, S. and Cai, H. (1996). Agricultural water management science. (Beijing: China Agriculture Press).
Kang, S. Z., Hao, X. M., Du, T. S., Tong, L., Su, X. L., Lu, H. N., et al. (2017). Improving agricultural water productivity to ensure food security in China under changing environment: From research to practice. Agr. Water Manage. 179, 5–17. doi: 10.1016/j.agwat.2016.05.007
Kim, C. P. and Stricker, J. N. M. (1996). Influence of spatially variable soil hydraulic properties and rainfall intensity on the water budget. Water Resour. Res. 32, 1699–1712. doi: 10.1029/96WR00603
Li, H., Li, C., Xing, K., Lei, Y., and Shen, Y. (2024a). Surface temperature adjustment in METRIC model for monitoring crop water consumption in North China Plain. Agr. Water Manage. 291, 108654. doi: 10.1016/j.agwat.2023.108654
Li, Y., Yin, Y., and Zhang, W. (2023). Water footprint assessment of major crops in Henan Province and reduction suggestions. Water 15, 1135. doi: 10.3390/w15061135
Li, M., Zhou, S., Shen, S., Wang, J., Yang, Y., Wu, Y., et al. (2024b). Climate-smart irrigation strategy can mitigate agricultural water consumption while ensuring food security under a changing climate. Agr. Water Manage. 292, 108663. doi: 10.1016/j.agwat.2023.108663
Liu, W. X., Ma, G., Wang, C. Y., Wang, J. R., Lu, H. F., Li, S. S., et al. (2018a). Irrigation and nitrogen regimes promote the use of soil water and nitrate nitrogen from deep soil layers by regulating root growth in wheat. Front. Plant Sci. 9. doi: 10.3389/fpls.2018.00032
Liu, W. X., Wang, J. R., Wang, C. Y., Ma, G., Wei, Q. R., Lu, H. F., et al. (2018b). Root growth, water and nitrogen use efficiencies in winter wheat under different irrigation and nitrogen regimes in north China plain. Front. Plant Sci. 9. doi: 10.3389/fpls.2018.01798
Liu, J. G., Yang, H., Gosling, S. N., Kummu, M., Flörke, M., Pfister, S., et al. (2017). Water scarcity assessments in the past, present, and future. Earth’s Future 5, 545–559. doi: 10.1002/2016ef000518
Liu, J. G., Zang, C. F., Tian, S. Y., Liu, J. G., Yang, H., Jia, S. F., et al. (2013). Water conservancy projects in China: Achievements, challenges and way forward. Glob. Environ. Change 23, 633–643. doi: 10.1016/j.gloenvcha.2013.02.002
Mo, X. G., Hu, S., Lin, Z. H., Liu, S. X., and Xia, J. (2017). Impacts of climate change on agricultural water resources and adaptation on the North China Plain. Adv. Clim. Change Res. 8, 93–98. doi: 10.1016/j.accre.2017.05.007
Mo, X. G., Liu, S. X., Lin, Z. H., Xu, Y. Q., Xiang, Y., and Mcvicar, T. R. (2005). Prediction of crop yield, water consumption and water use efficiency with a SVAT-crop growth model using remotely sensed data on the North China Plain. Ecol. Modell. 183, 301–322. doi: 10.1016/j.ecolmodel.2004.07.032
Pereira, L. S., Oweis, T., and Zairi, A. (2002). Irrigation management under water scarcity. Agr. Water Manage. 57, 175–206. doi: 10.1016/s0378-3774(02)00075-6
Piao, S. L., Ciais, P., Huang, Y., Shen, Z. H., Peng, S. S., Li, J. S., et al. (2010). The impacts of climate change on water resources and agriculture in China. Nature 467, 43–51. doi: 10.1038/nature09364
Rolim, J., Catalao, J., and Teixeira, J. (2011). The influence of different methods of interpolating spatial meteorological data on calculated irrigation requirements. Appl. Eng. Agric. 27, 979–989. doi: 10.13031/2013.40625
Scanlon, B. R., Faunt, C. C., Longuevergne, L., Reedy, R. C., Alley, W. M., Mcguire, V. L., et al. (2012). Groundwater depletion and sustainability of irrigation in the US High Plains and Central Valley. Proc. Natl. Acad. Sci. U.S.A. 109, 9320–9325. doi: 10.1073/pnas.1200311109
Sen, P. K. (1968). Estimates of the regression coefficient based on Kendall’s tau. J. Am. Stat. Assoc. 63, 1379–1389. doi: 10.2307/2285891
Sun, S., Shi, P. J., Zhang, Q., Wang, J. A., Wu, J. G., and Chen, D. L. (2022). Evolution of future precipitation extremes: Viewpoint of climate change classification. Int. J. Climatol. 42, 1220–1230. doi: 10.1002/joc.7298
Tan, Q., Liu, Y., Pan, T., Song, X., and Li, X. (2022). Changes and determining factors of crop evapotranspiration derived from satellite-based dual crop coefficients in North China Plain. Crop J. 10, 1496–1506. doi: 10.1016/j.cj.2022.07.013
Tan, L., Zhang, X., Qi, J., Sun, D., Marek, G. W., Feng, P., et al. (2023). Assessment of the sustainability of groundwater utilization and crop production under optimized irrigation strategies in the North China Plain under future climate change. Sci. Total Environ. 899, 165619. doi: 10.1016/j.scitotenv.2023.165619
Ti, J., Yang, Y., Yin, X., Liang, J., Pu, L., Jiang, Y., et al. (2018). Spatio-temporal analysis of meteorological elements in the North China District of China during 1960–2015. Water 10, 789. doi: 10.3390/w10060789
Vicente-Serrano, S. M., Saz-Sánchez, M. A., and Cuadrat, J. M. (2003). Comparative analysis of interpolation methods in the middle Ebro Valley (Spain): application to annual precipitation and temperature. Clim. Res. 24, 161–180. doi: 10.3354/cr024161
Wang, Y., Jiang, T., Bothe, O., and Fraedrich, K. (2007). Changes of pan evaporation and reference evapotranspiration in the Yangtze River basin. Theor. Appl. Climatol. 90, 13–23. doi: 10.1007/s00704-006-0276-y
Wang, M. and Shi, W. (2024b). The annual dynamic dataset of high-resolution crop water use in China from 1991 to 2019. Sci. Data 11, 1373. doi: 10.1038/s41597-024-04185-0
Wang, B., van Dam, J., Yang, X., Ritsema, C., Du, T., and Kang, S. (2023). Reducing water productivity gap by optimizing irrigation regime for winter wheat–summer maize system in the North China Plain. Agr. Water Manage. 280, 108229. doi: 10.1016/j.agwat.2023.108229
Wang, B., Wang, G., van Dam, J., Yang, X., Ritsema, C., Siddique, K. H. M., et al. (2024a). Diversified crop rotations improve crop water use and subsequent cereal crop yield through soil moisture compensation. Agr. Water Manage. 294, 108721. doi: 10.1016/j.agwat.2024.108721
Yang, X., Gao, W., Shi, Q., Chen, F., and Chu, Q. (2013). Impact of climate change on the water requirement of summer maize in the Huang-Huai-Hai farming region. Agr. Water Manage. 124, 20–27. doi: 10.1016/j.agwat.2013.03.017
Zhai, P. M., Zhang, X. B., Wan, H., and Pan, X. H. (2005). Trends in total precipitation and frequency of daily precipitation extremes over China. J. Clim. 18, 1096–1108. doi: 10.1175/jcli-3318.1
Zhang, X., Ding, B., Hou, Y., Feng, P., Li Liu, D., Srinivasan, R., et al. (2024). Assessing the feasibility of sprinkler irrigation schemes and their adaptation to future climate change in groundwater over-exploitation regions. Agr. Water Manage. 292, 108674. doi: 10.1016/j.agwat.2024.108674
Zhang, B., Liu, Y., Xu, D., Zhao, N., Lei, B., Rosa, R. D., et al. (2013). The dual crop coefficient approach to estimate and partitioning evapotranspiration of the winter wheat–summer maize crop sequence in North China Plain. Irrigation Sci. 31, 1303–1316. doi: 10.1007/s00271-013-0405-1
Zhang, Z. Y., Xiong, S. P., Wei, Y. H., Meng, X. D., Wang, X. C., and Ma, X. M. (2017). The role of glutamine synthetase isozymes in enhancing nitrogen use efficiency of N-efficient winter wheat. Sci. Rep. 7, 1000. doi: 10.1038/s41598-017-01071-1
Zhao, J., Yang, J., Huang, R., Xie, H., Qin, X., and Hu, Y. (2024). Estimating evapotranspiration and drought dynamics of winter wheat under climate change: A case study in Huang-Huai-Hai region, China. Sci. Total Environ. 949, 175114. doi: 10.1016/j.scitotenv.2024.175114
Keywords: double-cropping system, groundwater depletion, irrigation water requirement, sustainable agriculture, water balance framework
Citation: Zhang Z, Wang J, Ning J and Ti J (2025) Spatiotemporal dynamics of crop water use and groundwater depletion under the winter wheat–summer maize rotation system in Henan Province, China. Front. Plant Sci. 16:1728535. doi: 10.3389/fpls.2025.1728535
Received: 20 October 2025; Accepted: 28 November 2025; Revised: 28 November 2025;
Published: 11 December 2025.
Edited by:
Zhuanyun Si, Chinese Academy of Agricultural Sciences, ChinaReviewed by:
Xiaojun Shen, Tianjin Agricultural University, ChinaXiquan Wang, Inner Mongolia Agricultural University, China
Copyright © 2025 Zhang, Wang, Ning and Ti. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jinsong Ti, dGlqaW5zb25nQGhlbmF1LmVkdS5jbg==
Zhao Zhang