The Study of Drought in Future Climate Scenarios in the Huang-Huai-Hai Region

: In the context of global warming, agricultural production and social and economic development are signiﬁcantly affected by drought. The future change of climate conditions is uncertain; thus, it is of great importance to clarify the aspects of drought in order to deﬁne local and regional drought adaptation strategies. In this study, the meteorological data from 1976 to 2005 was used as a historical reference, and nine Global Climate Models (GCMs), downscaling to meteorological stations from 2039 to 2089, were used as future climate data. Based on Penman–Monteith, the reference crop Evapotranspiration (ET 0 ) and Standardized Precipitation Evapotranspiration Index (SPEI) of the reference crop in three emission scenarios of RCP2.6, RCP4.5, and RCP8.5, under future climate conditions, were calculated. A non-parameter Mann–Kendall trend test was performed on temperature, precipitation, ET 0 , and SPEI to analyze the drought spatiotemporal distribution traits under upcoming climate scenarios. The results showed that, under future climate conditions, SPEI values in most areas of the Huang-Huai-Hai region would continuously increase year by year, and drought would be alleviated to some extent at the same pace. However, with the increase of greenhouse gas concentration in the emission scenarios, SPEI values continued to decline. In the RCP8.5 scenario, the area of severe drought was large. To sum up, in the future climate scenario, the degree of drought in the Huang-Huai-Hai region will be alleviated to some extent with the increase of rainfall, but with the increase of greenhouse gas concentration, the degree of drought will be further intensi-ﬁed, posing a huge challenge to agricultural water use in the region. This study provides a theoretical foundation for alleviating drought in the Huang-Huai-Hai region in future climate scenarios.


Introduction
Under the general trend of global warming, extreme climate events occur frequently [1].The increase of temperature will accelerate water evaporation, causing drought, especially extreme drought, to become one of the most serious environmental issues in the world [2].According to statistics, the global economic losses caused by drought are as high as USD 6-8 billion every year, severely affecting the sustainable development of global agriculture [3,4].At the same time, drought will also lead to a shortage of water, which will have a consequential impact on water supply and agriculture [5].Meteorological hazards are the most serious natural disasters affecting agriculture in China, and drought accounts for Water 2021, 13, 3474 2 of 17 approximately 50% of these [6,7].The severe drought in the southwest China from 2009 to 2010 caused huge damage to local agriculture, society, the economy, and the ecological environment [8].The Huang-Huai-Hai region is one of the largest main grain producing areas in China.There have been several major droughts in this region, including the once-in-a-century drought in 1997, in which 137,000 hm 2 of summer grain crops dried up and died, and 6700 hm 2 of farmland failed to yield, severely affecting food safety and ecological balance in the region [9].Studies suggest that the range of extreme drought may increase from 1% of the current global area to 30% by the end of the 21st century [10].As global warming continues to intensify in the future, the negative impact of drought on food safety and agricultural water management will become more serious.Therefore, it is of a great significance to define the drought characteristics under future climate conditions and suggest reasonable countermeasures [11][12][13][14].
Most studies classify drought into five types: meteorological drought, agricultural drought, hydrological drought, socioeconomic drought, and groundwater drought.Meteorological drought is defined as the water deficit caused by the imbalance between total precipitation and surface evaporation in the study area over a certain period.Meteorological drought can lead to other types of droughts.Thus, accurate detection of meteorological drought can effectively prevent and alleviate drought [15,16].Drought index can measure, compare, and analyze the degree of drought, which plays an important role in drought research.Qu et al. [17] examined the spatiotemporal distribution aspects of drought intensity in different regions of Inner Mongolia based on the MCI index.Wei et al. [18] found that the accuracy of Pa, in describing drought intensity, was lower than that of PDSI in most areas of China.Zuo et al. [19] discussed the feasibility of using short-term precipitation series in China to calculate SPI, based on data from 2416 precipitation observation stations during the period 1961−2019.Tang et al. [20] used SPEI and Standardized Runoff Index (SRI), which was developed on the basis of SPI, to analyze the evolution of hydrological and meteorological droughts in southwest China during the period 1968−2018.However, SPI calculation only considers the precipitation impact on drought, disregarding changes caused by the temperature rise.Due to different calculation focuses of various drought indexes, there will be deviations in evaluation results [21].Yao et al. [22] found that drought severity was mainly affected by evapotranspiration in regions with substantial evapotranspiration changes in Xinjiang, China.SPEI considers, not only precipitation and temperature, but also evapotranspiration changes.Meanwhile, the calculation of SPEI is relatively simple and has more suitable application prospects [23].Tan et al. [24] used SPEI and SPI to study drought in Ningxia, China, and found that SPEI could describe the drought more accurately.Tirivarombo et al. [25] established that SPEI reflected a more comprehensive drought situation in the Kafue Basin.
The size of potential reference crop evapotranspiration (ET 0 ) affects the drought index SPEI.By studying the impacts of precipitation and evapotranspiration on global drought in the 21st century, Cook et al. [26] discovered that the increase of global drought area in the 2080s was mainly due to the increase of ET 0 .Therefore, it is very important to determine ET 0 , using appropriate methods to accurately evaluate the drought characteristics.At first, many scholars used the Thornthwaite equation to estimate evapotranspiration when calculating SPEI [27].The Thornthwaite equation is a method expressing the potential evapotranspiration by temperature, but evapotranspiration is also affected by precipitation, wind speed, humidity, and other factors besides temperature; thus, there will be some deviation when using this equation to assess evapotranspiration [28][29][30].In 1990, The Food and Agriculture Organization (FAO) included the Penman-Monteith (PM) equation as the standard calculation for evapotranspiration.This equation is more accurately based on physics, but it considers larger parameters and is more complicated to use [29].Bae et al. [31] calculated SPEI based on the Thornthwaite and Penman-Monteith equations and conducted comparative analysis.The results indicated that, when the meteorological data were the same, the drought degree calculated using the PM equation was more intense.
First, this study examined the precipitation and temperature changes in the Huang-Huai-Hai region, based on meteorological data in future scenarios, and then calculated the ET 0 and SPEI indexes under RCP2.6,RCP4.5, and RCP8.5 scenarios, based on the PM equation, and conducted MK trend tests.Finally, the spatiotemporal distribution and evolution of drought in the Huang-Huai-Hai region were quantified for future scenarios.This study can provide the basis for decision-making on reducing drought risk in the Huang-Huai-Hai region under future meteorological conditions.

Study Area
The Huang-Huai-Hai region, located in the east of China, is one of its main grain producing areas.The study area of this article mainly includes Henan, Hebei, Shandong, and Anhui, covering most provinces in the Huang-Huai-Hai region (Figure 1).This location is in the middle latitude zone of the monsoon climate zone in eastern China.In summer (June to August), the long-term mean temperature is 24.6 • C and the precipitation is 122.9 mm.In winter (December to February), the mean temperature is 2.3 • C and the precipitation is 9.8 mm.The annual mean temperature and precipitation decreased with the northern and southern latitudes increases [32].Considering global warming, the climate in the Huang-Huai-Hai region shows a warm and dry trend.The environment of insufficient water resources and increased accumulated temperature has brought many negative effects on local agricultural production.
Water 2021, 13, x FOR PEER REVIEW 3 of 18 equations and conducted comparative analysis.The results indicated that, when the meteorological data were the same, the drought degree calculated using the PM equation was more intense.
First, this study examined the precipitation and temperature changes in the Huang-Huai-Hai region, based on meteorological data in future scenarios, and then calculated the ET0 and SPEI indexes under RCP2.6,RCP4.5, and RCP8.5 scenarios, based on the PM equation, and conducted MK trend tests.Finally, the spatiotemporal distribution and evolution of drought in the Huang-Huai-Hai region were quantified for future scenarios.This study can provide the basis for decision-making on reducing drought risk in the Huang-Huai-Hai region under future meteorological conditions.

Study Area
The Huang-Huai-Hai region, located in the east of China, is one of its main grain producing areas.The study area of this article mainly includes Henan, Hebei, Shandong, and Anhui, covering most provinces in the Huang-Huai-Hai region (Figure 1).This location is in the middle latitude zone of the monsoon climate zone in eastern China.In summer (June to August), the long-term mean temperature is 24.6 °C and the precipitation is 122.9 mm.In winter (December to February), the mean temperature is 2.3 °C and the precipitation is 9.8 mm.The annual mean temperature and precipitation decreased with the northern and southern latitudes increases [32].Considering global warming, the climate in the Huang-Huai-Hai region shows a warm and dry trend.The environment of insufficient water resources and increased accumulated temperature has brought many negative effects on local agricultural production.

Data Source
In this study, the historical reference period is 1976-2005.The daily meteorological data of the Huang-Huai-Hai region were obtained from the China Meteorological Data Website (http://data.cma.cn/en,accessed on 6 October 2021).It mainly includes daily maximum and minimum temperature, near-surface relative humidity, near-surface wind speed, sunlight, and other factors of selected meteorological stations in the study area.Taking 2039-2089 as the study period for future scenarios, the method of downscaling was used [33] to obtain the meteorological data of nine climate models (Table 1) in the future, downscaling to station in three emission scenarios (RCP2.6,RCP4.5, and RCP8.5).The RCP is a series of comprehensive concentration and emission scenarios that are used as input parameters of the climate change predictive model, under the influence of 21st century human activities, to describe the emissions of greenhouse gases, reactive gases, and aerosols as well as concentration of atmospheric components under upcoming changes in population, socioeconomics, science and technology, energy consumption, and land use [34][35][36].The RCP2.6 represents the total radiative forcing to stabilize at 2.6 Wm −2 after 2100, which is a very low greenhouse gas concentration model.The RCP4.5 represents the total radiative forcing to stabilize at 4.5 Wm −2 after 2100, which is the climate scenario under government intervention.The RCP8.5 represents the total radiative forcing to stabilize at 8.5 Wm −2 after 2100, which is a baseline scenario in the absence of climate change policy intervention.In this study, the PM equation recommended by FAO was used to calculate ET 0 , and the calculation process is as follows [37,38]: where ET 0 denotes the reference crop evapotranspiration (mm/day), R n denotes the net radiation of crop surface (MJ/(m 2 .d)),G denotes soil heat flux (MJ/(m 2 .d)),T denotes the daily mean temperature at a height of 2 m ( • C), U 2 denotes the wind speed at a height of 2 m (m/s), e s denotes the saturated vapor pressure (kPa/ • C), and γ denotes the hygrometer constant (kPa/ • C).

Calculation of Standardized Precipitation Evapotranspiration Index (SPEI)
SPEI is a drought index that reflects water deficit at different time scales, based on ET 0 and precipitation.Therefore, SPEI can describe the aspects of drought caused by changes in meteorological factors such as precipitation and temperature.Table 2 shows the drought degree corresponding to the SPEI index in the different ranges in this article.For a specific month, water deficit, D, and rainfall, Pe, minus reference crop evapotranspiration, ET 0 .The equation is as follows: Monthly cumulative water deficit for a given year, i, with a time scale of 12 months is calculated as follows [23]: where X k i,j denotes the cumulative water deficit in month j in year i, k denotes the selected time scale, and k is 12 months in this study.
The probability distribution of the three-parameter log-logistic distribution was used to calculate SPEI, and the distribution function is [23]: where x is the value in the acceptable range (γ < x < ∞); the equations of α, β, and γ is as follows [23]: ) where i denotes the ascending sequence number of cumulative water deficit (X 1 ≤ X 2 . . .≤ X n ), N denotes the number of samples, and Γ(β) denotes the gamma function of β.Pe was defined as the part that exceeds the determined D, and SPEI was calculated [23]:

MK Test
The Mann-Kendall (MK) trend test is a non-parametric method [21], which has a wider test range, and the required data does not need to follow a specific distribution.It is commonly used to examine the variation trend of a data series with time [39].This study used this method to perform a test for annual mean temperature, annual mean Water 2021, 13, 3474 6 of 17 precipitation, annual mean ET 0 , and SPEI in the Huang-Huai-Hai region.The specific steps are as follows [40]: where X k and X l denote the values in the kth and lst year in the time series, respectively, n denotes the length of time series data, t p denotes the binding value corresponding to the pth digit, and Z c denotes the variation trend of time series data.
For a given confidence level, , this indicates rejection of the null hypothesis, suggesting that data series have a significant rising or descending trend; otherwise, the trend is not significant.When the standard statistical variable is Z c > 0, this indicates that the time series has a rising trend.When it is Z c < 0, this indicates that the time series has a decreasing trend [21].

Temperature Change
The spatiotemporal distribution of annual mean temperature in the Huang-Huai-Hai region is shown in Figure 2. The range of annual mean temperature ranges from 6.5 • C to 20.8 • C. On the whole, the annual mean temperature in this area gradually decreases from the south to the north, and the temperature in the south and east is higher, while that in the northwest is lower.In general, the temperature under the RCP8.5 scenario is higher than that under RCP2.6 and RCP4.5.Under the RCP2.6 scenario, the annual mean temperature is 6.5-17.2• C, 6.8-17.4• C, and 6.5-17.3• C in the 2030s, 2050s, and 2080s, respectively.Under RCP4.5, the mean temperature is 6.5−17.1 • C, 7.3-17.9• C, and 8.1-18.5 • C in the 2030s, 2050s, and 2080s, respectively.Under the RCP8.5 scenario, the mean temperature is 6.8−17.4• C, 8.0-18.5 • C, and 10.5-20.8• C in the 2030s, 2050s, and 2080s, respectively.The temperature difference between the three scenarios in the 2030s is smaller, and the spatial distribution pattern of temperature is similar.In the 2050s, the temperature of the three scenarios will increase by approximately 0.5 • C. In the 2080s, the temperature shows an obvious increase in the three scenarios.The annual mean temperature of RCP4.5 is 1.2-1.6 • C higher than that of RCP2.6, and that of RCP8.5 is 2.3 • C higher than that of RCP4.5.With the increase of time, the area with the higher temperature in the south of the region sees a gradual rise in temperature, while the area with lower temperature in the northwest of the Huang-Huai-Hai region sees a gradual decline in temperature.This trend is more evident under the RCP8.5 scenario.
the three scenarios will increase by approximately 0.5 °C.In the 2080s, the temperature shows an obvious increase in the three scenarios.The annual mean temperature of RCP4.5 is 1.2-1.6 °C higher than that of RCP2.6, and that of RCP8.5 is 2.3 °C higher than that of RCP4.5.With the increase of time, the area with the higher temperature in the south of the region sees a gradual rise in temperature, while the area with lower temperature in the northwest of the Huang-Huai-Hai region sees a gradual decline in temperature.This trend is more evident under the RCP8.5 scenario.Based on the MK trend test, the variation of mean temperature in future climate scenarios can be obtained, as shown in Figure 3.Among them, the variation trend of mean temperature in emission scenarios with different concentrations greatly differ.In the RCP2.6 scenario, the Zc value range of annual mean temperature is between −1.5 and 2.5, and the annual mean temperature in the western Huang-Huai-Hai region shows a downward trend.In the RCP4.5 scenario, the Zc value range of annual mean temperature is 7.8-8.8,and the annual mean temperature in the whole region shows a substantial increasing trend.In the RCP8.5 scenario, the Zc value range of annual mean temperature is 9.9-10.4.At this time, the annual mean temperature in Huang-Huai-Hai region shows an extremely significant increasing trend, especially in the northeast region, where the Zc value range of annual mean temperature is between 10.3 and 10.4.Based on the MK trend test, the variation of mean temperature in future climate scenarios can be obtained, as shown in Figure 3.Among them, the variation trend of mean temperature in emission scenarios with different concentrations greatly differ.In the RCP2.6 scenario, the Z c value range of annual mean temperature is between −1.5 and 2.5, and the annual mean temperature in the western Huang-Huai-Hai region shows a downward trend.In the RCP4.5 scenario, the Z c value range of annual mean temperature is 7.8-8.8,and the annual mean temperature in the whole region shows a substantial increasing trend.In the RCP8.5 scenario, the Z c value range of annual mean temperature is 9.9-10.4.At this time, the annual mean temperature in Huang-Huai-Hai region shows an extremely significant increasing trend, especially in the northeast region, where the Z c value range of annual mean temperature is between 10.

Precipitation Change
The spatiotemporal distribution of precipitation in the future climate scenarios shown in Figure 4.In terms of spatial distribution, precipitation successively increa from the north to the south, with less precipitation in the northwest of the Huang-Hu  In the future climate scenarios, the MK trend change of annual mean precipitation in the Huang-Huai-Hai region is shown in Figure 5, and the Zc value range of annual mean precipitation is between −0.73 and 5.3.Except for the decreasing trend of precipitation in some regions under the RCP2.6 scenario, almost all ranges in the Huang-Huai-Hai region show an increasing trend under the other two scenarios.In the RCP2.6 scenario, the Zc value range of annual mean precipitation is −0.73-1.6, and the southeast of the region shows a downward trend, especially at the junction of Henan and Anhui.Under the RCP4.5 scenario, the Zc value range of annual mean precipitation is −0.3-3.9.Except for a small part of Hebei, the precipitation in the Huang-Huai-Hai region presents a rising trend, and the increase amplitude in the western territory is greater than that in the eastern.In the RCP8.5 scenario, the Zc value range of annual mean precipitation is 1.6-5.3,and the precipitation in the whole Huang-Huai-Hai region presents an increasing trend.The distribution trend of precipitation is similar to that in RCP4.5, and the increase amplitude in the western region is greater than that in the eastern region.From the perspective of spatial alternation, the variation trend of precipitation in the region shows the obvious east-west difference.In possible events of all three emissions, the precipitation in the west In the future climate scenarios, the MK trend change of annual mean precipitation in the Huang-Huai-Hai region is shown in Figure 5, and the Z c value range of annual mean precipitation is between −0.73 and 5.3.Except for the decreasing trend of precipitation in some regions under the RCP2.6 scenario, almost all ranges in the Huang-Huai-Hai region show an increasing trend under the other two scenarios.In the RCP2.6 scenario, the Z c value range of annual mean precipitation is −0.73-1.6, and the southeast of the region shows a downward trend, especially at the junction of Henan and Anhui.Under the RCP4.5 scenario, the Z c value range of annual mean precipitation is −0.3-3.9.Except for a small part of Hebei, the precipitation in the Huang-Huai-Hai region presents a rising trend, and the increase amplitude in the western territory is greater than that in the eastern.In the RCP8.5 scenario, the Z c value range of annual mean precipitation is 1.6-5.3,and the precipitation in the whole Huang-Huai-Hai region presents an increasing trend.
The distribution trend of precipitation is similar to that in RCP4.5, and the increase amplitude in the western region is greater than that in the eastern region.From the perspective of spatial alternation, the variation trend of precipitation in the region shows the obvious east-west difference.In possible events of all three emissions, the precipitation in the west has a larger increasing trend.

Evolution of ET0 in the Huang-Huai-Hai Region under Future Climate Change Scenarios
In this study, the PM equation was used to calculate the reference crop evapotransp ration, ET0, in the Huang-Huai-Hai region under future climate scenarios.The spatiotem poral distribution of the annual mean ET0 in the region is shown in Figure 6.The spati distribution of annual mean ET0 in the Huang-Huai-Hai region is similar under the thr scenarios.The spatial distribution is higher in the east, especially in Shandong.It is rel tively lower in the northwest and the southwest, especially in northern Hebei and weste Henan.In the RCP2.6 scenario, the annual mean ET0 in the Huang-Huai-Hai region vari from 1611 mm to 2017 mm in the 2030s, from 1670 mm to 2058 mm in the 2050s, and fro 1618 mm to 2079 mm in the 2080s.ET0 continues to increase over time, and the increa from the 2030s to the 2050s is larger than that from the 2050s to the 2080s.However, th spatial distribution of ET0 in this scenario changes insignificantly between years.In th RCP4.5 scenario, the annual mean ET0 in the Huang-Huai-Hai region ranges from 16 mm to 1989 mm in the 2030s, from 1674 mm to 2094 mm in the 2050s, and from 1733 m to 2125 mm in the 2080s.From the 2050s to the 2080s, the increase amplitude of ET0 slower than that of the previous decade, but is larger than that in the RCP2.6 scenario.the RCP8.5 scenario, the annual mean ET0 in the Huang-Huai-Hai region varies from 16 mm to 2022 mm in the 2030s, from 1678 mm to 2104 mm in the 2050s, and from 1820 m to 2217 mm in the 2080s.The increase amplitude of ET0 during the period 2050s−2080s larger than that in the previous decade.In this context, the spatial distribution of ET0 va ies significantly between different decades.With the increase of decade, the area with th larger ET0 is constantly expanding; especially in the 2080s, ET0 in almost the entire Huan Huai-Hai region is above 2000 mm.By comparing the three scenarios, it can be seen th the area with the smaller ET0 in the three scenarios is continuously shrinking, while th

Evolution of ET 0 in the Huang-Huai-Hai Region under Future Climate Change Scenarios
In this study, the PM equation was used to calculate the reference crop evapotranspiration, ET 0 , in the Huang-Huai-Hai region under future climate scenarios.The spatiotemporal distribution of the annual mean ET 0 in the region is shown in Figure 6.The spatial distribution of annual mean ET 0 in the Huang-Huai-Hai region is similar under the three scenarios.The spatial distribution is higher in the east, especially in Shandong.It is relatively lower in the northwest and the southwest, especially in northern Hebei and western Henan.In the RCP2.6 scenario, the annual mean ET 0 in the Huang-Huai-Hai region varies from 1611 mm to 2017 mm in the 2030s, from 1670 mm to 2058 mm in the 2050s, and from 1618 mm to 2079 mm in the 2080s.ET 0 continues to increase over time, and the increase from the 2030s to the 2050s is larger than that from the 2050s to the 2080s.However, the spatial distribution of ET 0 in this scenario changes insignificantly between years.In the RCP4.5 scenario, the annual mean ET 0 in the Huang-Huai-Hai region ranges from 1622 mm to 1989 mm in the 2030s, from 1674 mm to 2094 mm in the 2050s, and from 1733 mm to 2125 mm in the 2080s.From the 2050s to the 2080s, the increase amplitude of ET 0 is slower than that of the previous decade, but is larger than that in the RCP2.6 scenario.In the RCP8.5 scenario, the annual mean ET 0 in the Huang-Huai-Hai region varies from 1627 mm to 2022 mm in the 2030s, from 1678 mm to 2104 mm in the 2050s, and from 1820 mm to 2217 mm in the 2080s.The increase amplitude of ET 0 during the period 2050s−2080s is larger than that in the previous decade.In this context, the spatial distribution of ET 0 varies significantly between different decades.With the increase of decade, the area with the larger ET 0 is constantly expanding; especially in the 2080s, ET 0 in almost the entire Huang-Huai-Hai region is above 2000 mm.By comparing the three scenarios, it can be seen that the area with the smaller ET 0 in the three scenarios is continuously shrinking, while the area with larger ET 0 is continuously expanding.This change trend is more evident in the RCP8.5 scenario.

Evolution of SPEI in the Huang-Huai-Hai Region under Future Climate Change
Figure 7 shows the spatial variation of SPEI under different scenarios in the decenniums (2030s, 2050s, and 2080s) in the Huang-Huai-Hai region.As can be the illustration of the RCP2.6 scenario, the annual mean SPEI in the Huang-Hu gion varies from −2.3 to 1.8 in the 2030s, from −2.5 to 2.4 in the 2050s, and from in the 2080s.In the RCP4.5 scenario, the annual mean SPEI in the Huang-Huaivaries from −1.8 to 2.1 in the 2030s, from −2.8 to 2.5 in the 2050s, and from −2.7 t 2080s.In the RCP8.5 scenario, the annual mean SPEI in the Huang-Huai-Hai reg from -3.0−2.7 in the 2030s, from −2.7 to 2.7 in the 2050s, and from −2.0 to 2.1 in In the RCP8.5 scenario, the drought degree reduces with the increase of de drought degree in the other two scenarios expands over time but is alleviate extent in the 2080s.From the perspective of spatial distribution, the area of seve in the north is larger, and the area of severe drought in the same period increase expansion of the concentration in the emission scenario, and this variation tren

Evolution of SPEI in the Huang-Huai-Hai Region under Future Climate Change Scenarios
Figure 7 shows the spatial variation of SPEI under different scenarios in the next three decenniums (2030s, 2050s, and 2080s) in the Huang-Huai-Hai region.As can be seen from the illustration of the RCP2.6 scenario, the annual mean SPEI in the Huang-Huai-Hai region varies from −2.3 to 1.8 in the 2030s, from −2.5 to 2.4 in the 2050s, and from −1.8 to 2.3 in the 2080s.In the RCP4.5 scenario, the annual mean SPEI in the Huang-Huai-Hai region varies from −1.8 to 2.1 in the 2030s, from −2.8 to 2.5 in the 2050s, and from −2.7 to 2.2 in the 2080s.In the RCP8.5 scenario, the annual mean SPEI in the Huang-Huai-Hai region varies from −3.0-2.7 in the 2030s, from −2.7 to 2.7 in the 2050s, and from −2.0 to 2.1 in the 2080s.In the RCP8.5 scenario, the drought degree reduces with the increase of decade.The drought degree in the other two scenarios expands over time but is alleviated to some extent in the 2080s.From the perspective of spatial distribution, the area of severe drought in the north is larger, and the area of severe drought in the same period increases with the expansion of the concentration in the emission scenario, and this variation trend is more apparent in the RCP2.6 scenario.As can be seen from Figure 8, the variation range of SPEI is smaller under RC the value of Zc ranging from −1.4 to 1.4.The variation range of SPEI is large under RC and RCP8.5, the value of Zc ranging from −2 to 2.1 and −2.8 to 2, respectively.Wit growth of emission scenario concentration, the area of increasing SPEI gradually expa Under the RCP2.6 scenario, most regions of Shandong and Anhui show a down trend, while Henan and Hebei show a slight upward trend.In the RCP4.5 scenario, in southern Anhui and the central Hebei decreases slightly, while that in the ce Huang-Huai-Hai region increases.Under the RCP8.5 scenario, the northwest China southern Anhui indicate an obvious trend of increasement, while Shandong, southern nan, and northern Anhui show an obvious trend of decline.As can be seen from Figure 8, the variation range of SPEI is smaller under RCP2.6, the value of Z c ranging from −1.4 to 1.4.The variation range of SPEI is large under RCP4.5 and RCP8.5, the value of Z c ranging from −2 to 2.1 and −2.8 to 2, respectively.With the growth of emission scenario concentration, the area of increasing SPEI gradually expands.Under the RCP2.6 scenario, most regions of Shandong and Anhui show a downward trend, while Henan and Hebei show a slight upward trend.In the RCP4.5 scenario, SPEI in southern Anhui and the central Hebei decreases slightly, while that in the central Huang-Huai-Hai region increases.Under the RCP8.5 scenario, the northwest China and southern Anhui indicate an obvious trend of increasement, while Shandong, southern Henan, and northern Anhui show an obvious trend of decline.As can be seen from Figure 8, the variation range of SPEI is smaller under RCP2.6, the value of Zc ranging from −1.4 to 1.4.The variation range of SPEI is large under RCP4.5 and RCP8.5, the value of Zc ranging from −2 to 2.1 and −2.8 to 2, respectively.With the growth of emission scenario concentration, the area of increasing SPEI gradually expands.Under the RCP2.6 scenario, most regions of Shandong and Anhui show a downward trend, while Henan and Hebei show a slight upward trend.In the RCP4.5 scenario, SPEI in southern Anhui and the central Hebei decreases slightly, while that in the central Huang-Huai-Hai region increases.Under the RCP8.5 scenario, the northwest China and southern Anhui indicate an obvious trend of increasement, while Shandong, southern Henan, and northern Anhui show an obvious trend of decline.Compared with the baseline level (1976−2005), the change of ET 0 level of representative stations in each province under future climate scenarios is shown in Figure 9.The results indicate that, compared with the baseline scenario (1976−2005), the ET 0 level in the 2030s, 2050s, and 2080s shows an upward trend.In the 2080s, the station with the largest growth of ET 0 is in Henan in the RCP2.6 scenario, with a growth rate of 10.7%.In the RCP4.5 scenario, Hebei has the largest growth rate of 13.5%, while in the RCP8.5 scenario, Anhui has the largest growth rate of 21%.The results of the MK trend test for the annual mean ET 0 under future climate scenarios are presented in Figure 9.It can be seen that the Z c values of ET 0 are all greater than zero, between 0.6 and 8.2, and the increase amplitude rises with the increase of concentration in emission scenarios.The Z c value of the ET 0 range is 0.6-4 in the RCP2.6 scenario.ET 0 in the eastern Huang-Huai-Hai Region increases importantly, especially in Anhui, Henan, and eastern Shandong.In the RCP4.5 scenario, the range of Z c value of ET 0 is 3-6, and ET 0 increases greatly in the northwest of the Huang-Huai-Hai region, especially in Hebei.In the RCP8.5 scenario, the range of Z c value of ET 0 is 5-8.2, and that in the southeast of Huang-Huai-Hai region notably increases, especially in Anhui.

Impact of the Future Climate Change on ET0
Compared with the baseline level (1976−2005), the change of ET0 level of representative stations in each province under future climate scenarios is shown in Figure 9.The results indicate that, compared with the baseline scenario (1976−2005), the ET0 level in the 2030s, 2050s, and 2080s shows an upward trend.In the 2080s, the station with the largest growth of ET0 is in Henan in the RCP2.6 scenario, with a growth rate of 10.7%.In the RCP4.5 scenario, Hebei has the largest growth rate of 13.5%, while in the RCP8.5 scenario, Anhui has the largest growth rate of 21%.The results of the MK trend test for the annual mean ET0 under future climate scenarios are presented in Figure 9.It can be seen that the Zc values of ET0 are all greater than zero, between 0.6 and 8.2, and the increase amplitude rises with the increase of concentration in emission scenarios.The Zc value of the ET0 range is 0.6-4 in the RCP2.6 scenario.ET0 in the eastern Huang-Huai-Hai Region increases importantly, especially in Anhui, Henan, and eastern Shandong.In the RCP4.5 scenario, the range of Zc value of ET0 is 3-6, and ET0 increases greatly in the northwest of the Huang-Huai-Hai region, especially in Hebei.In the RCP8.5 scenario, the range of Zc value of ET0 is 5-8.2, and that in the southeast of Huang-Huai-Hai region notably increases, especially in Anhui.Apparently, in different emission scenarios, the annual mean ET0 in most areas of the Huang-Huai-Hai region increases with the increase of decade.At the same time, the increased amplitude of annual mean ET0 in this area becomes higher with the increase of emission concentration in the scenarios.Li et al. [41] also predicted that ET0 in China would manifest an obvious growth gradient from the northeast to the southwest from 2011 to 2099.Liu et al. [42] studied the changes of ET0 in areas with different elevations over the Tibetan Plateau during the period 1961−2017 and found that the annual mean ET0 showed an upward trend within the study period, and mean temperature and wind speed played an important role in ET0 change.Wang et al. [43] found that ET0 in northwest China increased significantly with the evident increase of temperature.In this study, the variation trend of annual mean temperature in the Huang-Huai-Hai region under future meteorological conditions was compared.It was found that, with the increase of greenhouse gas concentration in the emission scenarios, both the mean temperature and the annual mean ET0 showed an increasing trend, but the Zc value spatial variation trend of temperature was not completely consistent with that of ET0.For example, under the RCP8.5 sce- Apparently, in different emission scenarios, the annual mean ET 0 in most areas of the Huang-Huai-Hai region increases with the increase of decade.At the same time, the increased amplitude of annual mean ET 0 in this area becomes higher with the increase of emission concentration in the scenarios.Li et al. [41] also predicted that ET 0 in China would manifest an obvious growth gradient from the northeast to the southwest from 2011 to 2099.Liu et al. [42] studied the changes of ET 0 in areas with different elevations over the Tibetan Plateau during the period 1961−2017 and found that the annual mean ET 0 showed an upward trend within the study period, and mean temperature and wind speed played an important role in ET 0 change.Wang et al. [43] found that ET 0 in northwest China increased significantly with the evident increase of temperature.In this study, the variation trend of annual mean temperature in the Huang-Huai-Hai region under future meteorological conditions was compared.It was found that, with the increase of greenhouse gas concentration in the emission scenarios, both the mean temperature and the annual mean ET 0 showed an increasing trend, but the Z c value spatial variation trend of temperature was not completely consistent with that of ET 0 .For example, under the RCP8.5 scenario, the mean temperature of Hebei in the north of the Huang-Huai-Hai region increases the most, while the annual mean ET 0 of Anhui in the south of the region essentially increases.This also indicates that temperature knowledge is one of the main factors affecting the change of ET 0 .According to the PM equation, we can also see that the main meteorological factors affecting ET 0 include air humidity, temperature, wind speed, and solar radiation.Zhang et al. [44] found that, in arid and semi-arid regions of China, ET 0 was more sensitive to air humidity, followed by solar radiation, temperature, and wind speed.Gao et al. [45] found that the change rates of air humidity, solar radiation, temperature, and wind speed were very small during the study period based on meteorological data from 2001 to 2050, thus there was no significant change of ET 0 in the RCP2.6 and R8.5 scenarios.Roderick and Farquhar [46] believed that the evaporation change rate of the northern hemisphere decreased steadily at 2-4 mm/a 2 in the past 50a, resulting in the phenomenon of "evaporation paradox".

The Drought Occurrence in Future Climate Scenarios
As can be seen from Figure 8, there are important differences in the changes of SPEI in different emission scenarios, suggesting that the aspects of drought in future climate scenarios are related to the concentration in the emission path.In the case of the high emission scenario (RCP8.5),SPEI in the Huang-Huai-Hai region shows a decreasing trend.SPEI in the same emission scenario also has significant spatial differences with time.Most droughts occur in the northern Huang-Huai-Hai region.The area of drought that occurred in THE RCP2.6 and RCP4.5 scenarios increases slightly with time, and that in the 2080s decreases significantly.In the RCP8.5 scenario, the area of drought has been decreasing over time.This might be due to the gradual increase of precipitation over time, which alleviates drought conditions in some areas.As the emission concentration increases, the temperature keeps rising, leading to the increase of ET 0 .As a result, the precipitation and ET 0 difference lessens more and more.Therefore, in the RCP8.5 scenario, the drought state in the Huang-Huai-Hai region is the most severe in the 2030s, while there is more of the moist situation in the region in the 2080s in the RCP2.6 scenario.Zhao et al. [47] used the SPEI index to analyze the drought in China from 2020 to 2099, and the conclusion was that there were large differences in the frequency of drought and moist conditions at various time scales.In this study, the percentages of stations with severe drought (Z c < −1.5) and severe moist (Z c > 1.5) in the lower Huang-Huai-Hai region in the 2030s, 2050s, and 2080s, under different emission scenarios, were estimated.It can be seen that, under the same emission scenario, the frequency of drought and moist conditions in this region also varies greatly at different time scales (Figure 10).In the RCP2.6 and RCP8.5 scenarios, there are more stations with severe drought in the 2030s.In the RCP4.5 scenario, there are more stations with severe drought in the 2050s.In the RCP8.5 scenario, there are more stations experiencing severe drought in the 2030s (Figure 10).According to the MK trend test results (Figure 8), SPEI decreases in a larger area in the RCP8.5 scenario, indicating that there are more drought areas, which is similar to the findings of Wang et al. [48].Li et al. [49] also predicted that the frequency, duration, severity and intensity of drought in the Huang-Huai-Hai plain would increase under the RCP8.5 scenario.Gao et al. [45] found that 12-month SPEI in Loess Plateau manifested an increasing trend under the RCP8.5 scenario, and the frequency of drought duration decreased significantly in the first half of the 21st century, which was generally consistent with the research results.The area of drought is relatively large in the 2030s, and decreases after the 2050s (Figures 7 and 10).

Conclusions
In future climate scenarios, the temperature and precipitation in the Huang-Huai-Hai region will apparently change.In most of this region, temperature will increase over time, and this trend is more significant in the RCP8.5 scenario.With the increase of emission concentration, areas with higher temperature gradually expand from the south to the north.In the RCP2.6 scenario, the area with the higher temperature is in Anhui.In the RCP4.5 scenario, areas with higher temperature are distributed in Anhui, Shandong, and eastern Henan.In the RCP8.5 scenario, areas with higher temperature include Anhui, Shandong, Henan, and southern Hebei.Most of the Huang-Huai-Hai region has higher temperature, especially in the 2080s.Precipitation also increases with the expansion of greenhouse gas concentration in the emission scenarios.Generally speaking, the growth trend in the west is greater than that in the east.Areas with higher precipitation are mostly distributed in the southern Huang-Huai-Hai region.From the perspective of spatial variation, the areas with higher ET0 in the three emission scenarios gradually increase from the east to the west with time.MK trend test results indicated that, with the increase of greenhouse gas concentration in the emission scenarios, the growth trend of ET0 would be more evident, and the growth trend in RCP8.5 would be much larger than the other two scenarios.
Under future climate conditions, due to the influence of temperature, precipitation, and other meteorological factors, drought in the Huang-Huai-Hai region presents various spatiotemporal distribution characteristics in different scenarios.The severity and area of drought in the RCP8.5 are the largest in the 2030s.This might be due to the continuous increase of precipitation and SPEI with time, and the drought is alleviated to some extent with the increase of decades.According to the MK trend test, SPEI has a great fluctuation with the increase of greenhouse gas concentration in emission scenarios.In the RCP8.5 scenario, the decrease amplitude and extent of SPEI are greater than those in the RCP2.6 and RCP4.5, indicating that drought is more severe in the high emission scenario.This might be caused by the decrease of water resources due to the increase of greenhouse gas concentration, temperature, and ET0 in different emission scenarios.
According to the results of this study, in future climate scenarios, the degree of drought in the Huang-Huai-Hai region will be alleviated to some extent due to the increase of rainfall.However, if greenhouse gas emissionsare not controlled, with the increase of greenhouse gas, the continuously increasing temperature will lead to water shortage and aggravate drought.The government should determine measures to reduce

Conclusions
In future climate scenarios, the temperature and precipitation in the Huang-Huai-Hai region will apparently change.In most of this region, temperature will increase over time, and this trend is more significant in the RCP8.5 scenario.With the increase of emission concentration, areas with higher temperature gradually expand from the south to the north.In the RCP2.6 scenario, the area with the higher temperature is in Anhui.In the RCP4.5 scenario, areas with higher temperature are distributed in Anhui, Shandong, and eastern Henan.In the RCP8.5 scenario, areas with higher temperature include Anhui, Shandong, Henan, and southern Hebei.Most of the Huang-Huai-Hai region has higher temperature, especially in the 2080s.Precipitation also increases with the expansion of greenhouse gas concentration in the emission scenarios.Generally speaking, the growth trend in the west is greater than that in the east.Areas with higher precipitation are mostly distributed in the southern Huang-Huai-Hai region.From the perspective of spatial variation, the areas with higher ET 0 in the three emission scenarios gradually increase from the east to the west with time.MK trend test results indicated that, with the increase of greenhouse gas concentration in the emission scenarios, the growth trend of ET 0 would be more evident, and the growth trend in RCP8.5 would be much larger than the other two scenarios.
Under future climate conditions, due to the influence of temperature, precipitation, and other meteorological factors, drought in the Huang-Huai-Hai region presents various spatiotemporal distribution characteristics in different scenarios.The severity and area of drought in the RCP8.5 are the largest in the 2030s.This might be due to the continuous increase of precipitation and SPEI with time, and the drought is alleviated to some extent with the increase of decades.According to the MK trend test, SPEI has a great fluctuation with the increase of greenhouse gas concentration in emission scenarios.In the RCP8.5 scenario, the decrease amplitude and extent of SPEI are greater than those in the RCP2.6 and RCP4.5, indicating that drought is more severe in the high emission scenario.This might be caused by the decrease of water resources due to the increase of greenhouse gas concentration, temperature, and ET 0 in different emission scenarios.
According to the results of this study, in future climate scenarios, the degree of drought in the Huang-Huai-Hai region will be alleviated to some extent due to the increase of rainfall.However, if greenhouse gas emissionsare not controlled, with the increase of greenhouse gas, the continuously increasing temperature will lead to water shortage and aggravate drought.The government should determine measures to reduce greenhouse gas emissions, develop water-saving agriculture, carefully allocate water resources to alleviate

Figure 2 .
Figure 2. Spatiotemporal distribution of temperature under future climate scenarios.

Figure 2 .
Figure 2. Spatiotemporal distribution of temperature under future climate scenarios.

Figure 3 .
Figure 3.The MK trend test of temperature under future climate scenarios.

Figure 3 .
Figure 3.The MK trend test of temperature under future climate scenarios.

Figure 4 .
Figure 4. Spatiotemporal distribution of precipitation under future climate scenarios.

Figure 4 .
Figure 4. Spatiotemporal distribution of precipitation under future climate scenarios.

Figure 5 .
Figure 5. MK trend test of precipitation under future climate scenarios.

Figure 5 .
Figure 5. MK trend test of precipitation under future climate scenarios.

Figure 6 .
Figure 6.Spatiotemporal distribution of ET0 under future climate scenarios.

Figure 6 .
Figure 6.Spatiotemporal distribution of ET 0 under future climate scenarios.

Figure 8 .
Figure 8. MK trend test of SPEI under future climate scenarios.

Figure 8 .
Figure 8. MK trend test of SPEI under future climate scenarios.Figure 8. MK trend test of SPEI under future climate scenarios.

Figure 8 .
Figure 8. MK trend test of SPEI under future climate scenarios.Figure 8. MK trend test of SPEI under future climate scenarios.

4 . Discussion 4 . 1 .
Impact of the Future Climate Change on ET 0

Figure 9 .
Figure 9. Zc values of ET0 in the MK trend test and the rate of decadal ET0 variation in major stations in different RCP scenarios.

Figure 9 .
Figure 9. Z c values of ET 0 in the MK trend test and the rate of decadal ET 0 variation in major stations in different RCP scenarios.

Figure 10 .
Figure 10.The variation trend of drought and moist under RCP scenarios in the Huang-Huai-Hai region (red represents severe drought and blue stands for severe moist; colors from light to dark represent 2030s, 2050s, and 2080s).

Figure 10 .
Figure 10.The variation trend of drought and moist under RCP scenarios in the Huang-Huai-Hai region (red represents severe drought and blue stands for severe moist; colors from light to dark represent 2030s, 2050s, and 2080s).

Table 1 .
Climate models used in this study.