Effects on Net Irrigation Water Requirement of Joint Distribution of Precipitation and Reference Evapotranspiration

: To establish the uncertain inﬂuence that the joint distribution of precipitation and reference evapotranspiration has on net irrigation water requirement, a Copula function–Monte Carlo method (CFMC) was proposed to calculate the probability of irrigation water requirement. Taking the Jingdian Irrigation District in Northwest China as an example, the distribution laws of precipitation and reference evapotranspiration were studied. Furthermore, ﬁve typical years under different crop planting structure conditions were selected, and the variation characteristics of net irrigation water requirement in each typical year under the conditions of climate uncertainty were analyzed. The results revealed the optimal distribution functions of precipitation and reference evapotranspiration to be gamma distribution and lognormal distribution. The probability density map of the joint distribution of precipitation and reference evapotranspiration has a “saddle” shape; that is, irrigation water requirement and reference evapotranspiration are usually inversely related. As the probability of the irrigation water requirement increases, the net irrigation water requirement in the irrigation area also increases. The CFMC method can determine the design value of the net irrigation water requirement under a speciﬁc probability for typical years under different crop planting structure conditions, which can provide a reference for agricultural water resource allocation in irrigation areas.


Introduction
The irrigated area in Asia has increased steadily during the last 50 years [1,2]. The use of agricultural water has increased significantly as the irrigated area has increased, leading to increasing water scarcity in a number of regions [3][4][5]. Some developing countries use over 90 percent of their water for irrigation [6]. As a large agricultural country, the average water consumption per unit area in China is only 50% of the global average level [7], and the issue of irrigation water shortage has become increasingly obvious [8,9]. Due to the greater demand for irrigation, the contradiction between the supply and demand of agricultural water resources has intensified [10][11][12][13]. Crop water demand and net irrigation water requirement are important bases for the management planning and design of irrigation water. Therefore, the rational design of net irrigation water requirement can serve to improve irrigation water use efficiency [14][15][16].
Many studies have been conducted on calculating the net irrigation water requirement. Liu et al. used the crop coefficient method and empirical formula for the calculation of crop water demand, effective precipitation, and net irrigation water requirement of major crops in China before analyzing their spatial distribution characteristics [17]. Nkomozepi et al. used a sensitivity analysis method for studying the impact that climate change has on the net irrigation water requirement in the natural agro-ecological areas of Zimbabwe [18]. Wang et al. used meteorological data and crop data for calculating the water requirement of maize in the Chengde area, analyzing the influence meteorological factors have on the water demand of spring maize [19]. Liu et al. calculated the net irrigation water requirement for cotton in Hebei Province, China to be an average of 190.6 mm by using the water balance method, and the irrigation demand index was found to be 0.29 [20]. Wang et al. conducted research on the net irrigation water requirement of rice, finding precipitation to be the primary factor that affects the net irrigation water requirement, followed by crop water demand [21]. Liu et al. used the path analysis method to determine the causes of water requirements of staple crops in Henan Province, China [22]. Xie et al. decomposed the driving factors of the net irrigation water requirement in Wuwei City using the LMDI method, also calculating the effect values of the influencing factors of net irrigation water requirement [23]. Feng et al. used the water balance method for the calculation of the net irrigation water requirement in the Sanyizhai Yellow River Diversion Irrigation District, analyzing its temporal and spatial variation [24]. In addition, several studies have attempted the use of crop models including CROPWAT for the estimation of the irrigation needs of a variety of crops in different agro-ecological units [25,26]. The CROPWAT model estimates crop water demands during different stages and uses input climate, soil, and crop data for the generation of irrigation schedules [27].
The majority of the aforementioned studies use the crop coefficient method and Penman-Monteith formula for the calculation of crop water demand and the empirical formula for the estimation of effective precipitation. The net irrigation water requirement is generally calculated by a water balance equation [28,29]. Precipitation and reference evapotranspiration are two of the key factors in the design and implementation of irrigation schedules, and the relationship between them has important implications for the study of the uncertainty of net irrigation water requirements [30]. The Penman-Monteith formula [31] only requires simple meteorological data for the calculation of reference evapotranspiration. This formula has the advantage of high calculation accuracy and is therefore widely used [32][33][34][35]. As they are both affected by meteorological factors, crop water demand and precipitation are uncertain. The frequency distribution curve is a simple one-way frequency analysis method [36] that is widely used in hydrological research [37][38][39][40]. Due to the correlation between meteorological factors, the design of irrigation schedules must be based on the joint probability distribution of precipitation and crop water demand [41]. Zhang et al. demonstrated that a copula function can be used for the accurate description of the joint probability distribution characteristics of precipitation and reference evapotranspiration and for determining the joint return period of both [30]. However, when the joint return period is determined by a copula function, it is difficult to apply it directly to the design of the net irrigation water requirement due to the infinite possible combinations of precipitation and reference evapotranspiration during a given return period. Therefore, for the utilization efficiency of agricultural water resources to be improved, a calculation method for net irrigation water requirement where uncertainty is considered is urgently required.
In this study, a copula function and Monte Carlo simulation were combined to construct a net irrigation water requirement calculation method considering the joint probability distribution of precipitation and reference evapotranspiration, which can be used to determine the net irrigation water requirement with different probabilities of irrigation. Taking the Jingdian Irrigation District in Gansu Province, China as an example, five typical years under different crop planting structure conditions were selected, and the variation characteristics of net irrigation water requirement in each typical year from 2000 to 2020 were analyzed. This study will provide a reference for the calculation of water demand in irrigation areas under the condition of climate uncertainty, which is of great significance to the efficient use of water resources in irrigation areas.

Calculation of Net Irrigation Water Requirement
The water balance method is widely used for the calculation of the net irrigation water requirement. When groundwater depth is more than 5 m, the impact groundwater recharge has can generally be ignored [34], and the calculation formula that is used for the net irrigation water requirement during the growing period is expressed as [17,42]: where IR w is the net irrigation water requirement of the whole growing period, mm; EP w is the effective precipitation of the whole growing period, mm; and ETc w is the crop water demand of the whole growing period, mm. In Equation (1), effective precipitation refers to the part of precipitation that actually enters the root layer and can be used by crops, excluding surface runoff and deep percolation. When Equation (1) is used for the calculation of the net irrigation water requirement, only effective precipitation during the growth period and crop water demand need to be calculated.
In arid regions with less snowfall, when Equation (1) is used for the calculation of net irrigation water requirement during the entire growth period, the influence of soil moisture changes can generally be ignored. However, when the net irrigation water requirement is calculated for short periods (including daily or monthly scales), changes in soil moisture must be considered, and monthly net irrigation water requirement is then expressed as [34,43]: where IR i is the net irrigation water requirement in month i, mm; ETc i is the crop water demand in month i, mm; EP i is the effective precipitation in month i, mm; and S i is the soil water storage at the end of month i, mm. Monthly effective precipitation is calculated using the empirical formula [44,45]: where AP i represents cumulated precipitation in month i, mm, which is converted from the daily precipitation data of the weather station. Monthly crop water demand is calculated based on the daily crop water demand, and the daily crop water demand is calculated using the crop coefficient method: where n i is the number of days in the month i; Kc j is the crop coefficient on day j; and ETo j is the reference evapotranspiration on day j, mm, which can be calculated using the Penman-Monteith formula that is recommended by FAO [31]: where R n is the net radiation on the crop surface, MJ/(m 2 ·d); G is the soil heat flux, MJ/(m 2 ·d); T is the average temperature at a height of 2 m, • C; u 2 is the wind speed at a height of 2 m, m/d; e s is the saturated water pressure, kPa; e a is the actual vapor pressure, kPa; ∆ is the slope of the saturated water pressure and temperature curve, kPa/°C; and λ is the psychrometric constant, kPa/ • C. The crop coefficient refers to the ratio of the crop water demand to the reference evapotranspiration in a certain period of time. Different crops have different crop coefficients in where n k is the number of crop types in the study area; Kc k is the crop coefficient of the k-th crop; A k is the planting area of the k-th crop, km 2 ; and A is the total planting area, km 2 . The monthly effective precipitation and monthly reference evapotranspiration can be calculated using Equations (3) and (4); however, according to Equation (2), it is necessary to calculate the changes of soil water storage in each month to obtain the monthly net irrigation water demand. According to the principle of water balance, the monthly soil water storage can be calculated by the following equation: According to the above equations, as long as the initial soil water storage is determined, the soil water storage in any month can be calculated. In this study, the monthly effective precipitation one month before the crop growth period is used as the initial soil water storage.

Joint Probability Distribution Based on a Copula Function
Effective precipitation and reference evapotranspiration are both key variables that affect net irrigation water requirement. The two variables are affected by meteorological conditions and have a certain correlation. A copula function is an effective tool for the construction of multivariate joint distribution and can completely retain the correlation information between variables. It also possesses strong flexibility and excellent applicability [46]. Therefore, a copula function is used in this study for establishing the joint probability distribution model of the monthly effective precipitation and reference evapotranspiration. The use of a copula function for establishing the joint probability distribution is divided into two stages. The first is the determination of the optimal distribution of monthly effective precipitation and reference evapotranspiration, and the second is the determination of copula function parameters.
Gamma distribution, lognormal distribution, and generalized extreme value (GEV) distribution are the marginal distribution types that are most commonly used for hydrological analysis [47]. The probability density functions of the three distributions are expressed as: (1) Gamma distribution: where α and β are shape and scale parameters.
(2) Lognormal distribution: where µ and σ are the mean and standard deviation of the logarithm of the variable.
(3) GEV distribution: where k, σ, and µ are shape, scale, and position parameters. The K-S test method (Kolmogorov-Smirnov) is used for testing the hypothesis of marginal distribution, and the Akaike information criterion (AIC) is used for selecting the optimal marginal distribution. The model with the smallest AIC value should be selected from the alternative functions as the optimal model. The expression of the AIC value is [47]: where N is the total number of independent parameters, and L is the maximum likelihood value.
Assuming that X and Y represent effective precipitation and reference evapotranspiration, the corresponding design values are x and y, and their marginal distributions are F X (x) and F Y (y). According to Sklar's theorem, the joint distribution function of X and Y is expressed using a two-dimensional copula function [48]: Copula functions are divided into elliptic, quadratic, and Archimedean types, and Archimedean copula functions are those that are most widely used in the field of hydrology. These include the Gumbel copula, Clayton copula, and Frank copula functions. The first two are generally used for positive correlation variables, whereas the Frank copula is suitable for negative correlation variables. A negative correlation generally exists between precipitation and reference evapotranspiration. Therefore, the Frank copula function was chosen for the construction of the joint distribution function, which is expressed in the following way [49]: where u and v represent the cumulative probability of marginal distribution for two variables; and θ is a shape parameter.

Calculation of Probability of Irrigation Water Requirement
The probability of irrigation water requirement refers to the probability that the water consumption of the irrigation area can be fully satisfied. Usually, it is expressed as the percentage of the number of years in which the design irrigation water consumption is fully satisfied to the total number of years [43]. The probability of irrigation water requirement is usually determined by the typical year method. The selection of the typical year usually only considers the amount of precipitation and does not consider the uneven distribution of precipitation in time. The problem of the uneven distribution of precipitation can only be solved when there are enough precipitation samples. The Monte Carlo method is based on the theory of probability and statistics and is used for the analysis of a variety of uncertain problems. Monte Carlo methods are often used to solve uncertain problems. The method is to generate possible situations through a large number of simulations so as to obtain different results and probability distributions under a certain probability. Then, by performing statistical analysis on a large number of sample values, a result that meets a certain accuracy can be obtained. Therefore, Monte Carlo simulation is a powerful tool for uncertain and risky problems. When using the Monte Carlo method, it is necessary to input a given distribution function and parameters and then use a computer to generate a series of random numbers that satisfy the distribution function. The more random samples there are, the more likely it is to cover all situations. The given distribution function can be a univariate distribution function or a multivariate distribution function. The Monte Carlo method is used in this paper for evaluating the impact that sample uncertainty has on the net irrigation water requirement. Combining the Monte Carlo method and a copula function, the Copula function-Monte Carlo method (CFMC) was established as shown in Figure 1. From Figure 1, the following steps are taken for determining the probability of the irrigation water requirement: Step 1: Calculate the reference evapotranspiration according to the meteorological data and the Penman-Monteith formula.
Step 2: Determine the optimal marginal distribution model and parameters of monthly precipitation and reference evapotranspiration based on the Akaike information criterion.
Step 3: Determine the joint distribution parameters of precipitation and reference evapotranspiration based on the Frank copula function, use the joint distribution model to represent the population, randomly generate 1000 sets of cumulative probability samples of the joint distribution function, and then use the marginal distribution model of precipitation and reference evapotranspiration. Calculate the corresponding precipitation and reference evapotranspiration.
Step 4: Select a typical year under a specific crop planting structure and use Equation (6) to calculate the crop coefficient. Take 1000 groups of randomly simulated precipitation and reference evapotranspiration data as input items, and calculate the net irrigation water demand corresponding to 1000 groups of samples.
Step 5: Calculate the probability of irrigation water requirement according to the 1000 groups of net irrigation water requirements so that the net irrigation water requirements corresponding to different probabilities of irrigation water requirements can be determined. The formula for calculating the probability of irrigation water requirement is: From Figure 1, the following steps are taken for determining the probability of the irrigation water requirement: Step 1: Calculate the reference evapotranspiration according to the meteorological data and the Penman-Monteith formula.
Step 2: Determine the optimal marginal distribution model and parameters of monthly precipitation and reference evapotranspiration based on the Akaike information criterion.
Step 3: Determine the joint distribution parameters of precipitation and reference evapotranspiration based on the Frank copula function, use the joint distribution model to represent the population, randomly generate 1000 sets of cumulative probability samples of the joint distribution function, and then use the marginal distribution model of precipitation and reference evapotranspiration. Calculate the corresponding precipitation and reference evapotranspiration.
Step 4: Select a typical year under a specific crop planting structure and use Equation (6) to calculate the crop coefficient. Take 1000 groups of randomly simulated precipitation and reference evapotranspiration data as input items, and calculate the net irrigation water demand corresponding to 1000 groups of samples.
Step 5: Calculate the probability of irrigation water requirement according to the 1000 groups of net irrigation water requirements so that the net irrigation water require-Agriculture 2022, 12, 801 7 of 16 ments corresponding to different probabilities of irrigation water requirements can be determined. The formula for calculating the probability of irrigation water requirement is: where P(Q) is the probability of the irrigation water requirement corresponding to the net irrigation water requirement of Q, %; n is the total number of years, where n = 1000; and m is the number of years when the net irrigation water requirement is greater than Q.
Step 6: Select different typical years and repeat steps 4 to 5 to obtain the relationship between the net irrigation water requirement and the probability of the irrigation water requirement under different crop planting structures.

Study Area
The Jingdian Irrigation District is situated in the middle of Gansu Province, China ( Figure 2). Its north side is bordered by the Tengger Desert, its east side is in close proximity to the Yellow River, and at its south side are the Changling mountains. According to raw National Meteorological data for Jingtai County, the annual average precipitation in the irrigation area is 186 mm, the annual variation range is 103.5-298.4 mm, and the annual water surface evaporation intensity is 2123 mm (E601 evaporating dish). The annual average temperature is 8.8 • C. Figure 3 shows the monthly average precipitation and temperature. The maximum and minimum temperatures in the Jingdian Irrigation District occur in July and January, respectively, and the maximum and minimum precipitation occur in August and December. The total area of the Jingdian Irrigation District is 1496 km 2 , and the irrigation area is 607.7 km 2 , which includes 201.1 km 2 in phase I and 406.6 km 2 in phase II. The main crop types that are found in the irrigation area include spring wheat, barley, corn, soybean, flax, potato, melon, medlar, and forest fruit. From the test results of soil samples, it was found that the soil texture of the root layer in the irrigation area is mainly sandy loam. where P(Q) is the probability of the irrigation water requirement corresponding to the net irrigation water requirement of Q, %; n is the total number of years, where n = 1000; and m is the number of years when the net irrigation water requirement is greater than Q.
Step 6: Select different typical years and repeat steps 4 to 5 to obtain the relationship between the net irrigation water requirement and the probability of the irrigation water requirement under different crop planting structures.

Study Area
The Jingdian Irrigation District is situated in the middle of Gansu Province, China ( Figure 2). Its north side is bordered by the Tengger Desert, its east side is in close proximity to the Yellow River, and at its south side are the Changling mountains. According to raw National Meteorological data for Jingtai County, the annual average precipitation in the irrigation area is 186 mm, the annual variation range is 103.5-298.4 mm, and the annual water surface evaporation intensity is 2123 mm (E601 evaporating dish). The annual average temperature is 8.8 °C. Figure 3 shows the monthly average precipitation and temperature. The maximum and minimum temperatures in the Jingdian Irrigation District occur in July and January, respectively, and the maximum and minimum precipitation occur in August and December. The total area of the Jingdian Irrigation District is 1496 km 2 , and the irrigation area is 607.7 km 2 , which includes 201.1 km 2 in phase I and 406.6 km 2 in phase II. The main crop types that are found in the irrigation area include spring wheat, barley, corn, soybean, flax, potato, melon, medlar, and forest fruit. From the test results of soil samples, it was found that the soil texture of the root layer in the irrigation area is mainly sandy loam.

Data Source
Meteorological data for the Jingdian Irrigation District from 1960 to 2020, which includes rainfall (mm), average temperature (°C), sunshine hours (h), relative humidity (%), and wind speed (m/s), were obtained from the Jingtai County National Meteorological Station and China Meteorological Administration. The data for the crop planting area and planting structure in the irrigation area between 2000 and 2020 were obtained from the Jingtaichuan Electric Power Lifting Irrigation Water Resources Utilization Center, and the crop coefficient (Table 1) was obtained from [34].

Data Source
Meteorological data for the Jingdian Irrigation District from 1960 to 2020, which includes rainfall (mm), average temperature ( • C), sunshine hours (h), relative humidity (%), and wind speed (m/s), were obtained from the Jingtai County National Meteorological Station and China Meteorological Administration. The data for the crop planting area and planting structure in the irrigation area between 2000 and 2020 were obtained from the Jingtaichuan Electric Power Lifting Irrigation Water Resources Utilization Center, and the crop coefficient (Table 1) was obtained from [34].   The daily crop coefficient calculated based on the crop planting area in the Jingdian Irrigation District is shown in Figure 5. The crop growth period in the entire irrigation area is March to October. The crop coefficient for each typical year increases gradually between March and June and decreases gradually between July and October. From 2000 to 2020, the crop coefficient for the months between March and June demonstrated a yearto-year decrease, and the crop coefficient for each month between July and October demonstrated a year-to-year increase. This is mainly due to the growth period of spring wheat being March to June. As the planting proportion of spring wheat decreases, the crop coefficient between March and June decreases each year. The growth period for maize is from July to October. As the maize planting proportion increases, the crop coefficient between July and October also increases each typical year.  Figure 6 displays the calculation results for crop water demand. The crop water demand of each typical year demonstrates a single peak change, and the peak value concentration is from May to June. This is due to the development period for most crops in the irrigation area being concentrated between May and June and the crop coefficient during the development period being relatively large, meaning that the reference evapotranspiration is also large. In addition, as a result of the influence of climate factors, the peak value of reference evapotranspiration is also in June, meaning that the peak value of crop The daily crop coefficient calculated based on the crop planting area in the Jingdian Irrigation District is shown in Figure 5. The crop growth period in the entire irrigation area is March to October. The crop coefficient for each typical year increases gradually between March and June and decreases gradually between July and October. From 2000 to 2020, the crop coefficient for the months between March and June demonstrated a year-to-year decrease, and the crop coefficient for each month between July and October demonstrated a year-to-year increase. This is mainly due to the growth period of spring wheat being March to June. As the planting proportion of spring wheat decreases, the crop coefficient between March and June decreases each year. The growth period for maize is from July to October. As the maize planting proportion increases, the crop coefficient between July and October also increases each typical year.  The daily crop coefficient calculated based on the crop planting area in the Jingdian Irrigation District is shown in Figure 5. The crop growth period in the entire irrigation area is March to October. The crop coefficient for each typical year increases gradually between March and June and decreases gradually between July and October. From 2000 to 2020, the crop coefficient for the months between March and June demonstrated a yearto-year decrease, and the crop coefficient for each month between July and October demonstrated a year-to-year increase. This is mainly due to the growth period of spring wheat being March to June. As the planting proportion of spring wheat decreases, the crop coefficient between March and June decreases each year. The growth period for maize is from July to October. As the maize planting proportion increases, the crop coefficient between July and October also increases each typical year.  Figure 6 displays the calculation results for crop water demand. The crop water demand of each typical year demonstrates a single peak change, and the peak value concentration is from May to June. This is due to the development period for most crops in the irrigation area being concentrated between May and June and the crop coefficient during the development period being relatively large, meaning that the reference evapotranspiration is also large. In addition, as a result of the influence of climate factors, the peak value of reference evapotranspiration is also in June, meaning that the peak value of crop  Figure 6 displays the calculation results for crop water demand. The crop water demand of each typical year demonstrates a single peak change, and the peak value concentration is from May to June. This is due to the development period for most crops in the irrigation area being concentrated between May and June and the crop coefficient during the development period being relatively large, meaning that the reference evapotranspiration is also large. In addition, as a result of the influence of climate factors, the peak value of reference evapotranspiration is also in June, meaning that the peak value of crop water demand each year is mainly concentrated in the same month. From 2000 to 2020, the crop water demand between March and June mainly demonstrated an increasing trend each year, whereas the crop water demand between July and October demonstrated a yearly decreasing trend. This is due to the change of the planting proportions of spring wheat and corn, resulting in the crop coefficient and crop water demand changing accordingly.

Effects of Crop Planting Structure on Net Irrigation Water Requirement
water demand each year is mainly concentrated in the same month. From 2000 to 2020, the crop water demand between March and June mainly demonstrated an increasing trend each year, whereas the crop water demand between July and October demonstrated a yearly decreasing trend. This is due to the change of the planting proportions of spring wheat and corn, resulting in the crop coefficient and crop water demand changing accordingly.

Joint Distribution of Precipitation and Reference Evapotranspiration
From data relating to the monthly precipitation and reference evapotranspiration of Jingdian Irrigation District between 1960 and 2020, the marginal distribution parameters of precipitation and reference evapotranspiration can be fitted. Table 2 shows that the monthly precipitation for every month, with the exception of March, passed the K-S test at the significance level of 0.05, which indicates that the fitting effects of the three distribution functions are good (Figure 7). The AIC value is used for judging the optimal precipitation distribution function. The AIC of gamma distribution or most months is the minimum, so gamma distribution is the optimal distribution function of monthly precipitation. Table 3 shows that the reference evapotranspiration for every month passed the K-S test at the significance level of 0.05, which indicates that the fitting effects of the three distribution functions are good (Figure 8). The AIC of lognormal distribution for most months is the minimum, meaning that the optimal distribution function of monthly reference evapotranspiration is lognormal distribution.

Joint Distribution of Precipitation and Reference Evapotranspiration
From data relating to the monthly precipitation and reference evapotranspiration of Jingdian Irrigation District between 1960 and 2020, the marginal distribution parameters of precipitation and reference evapotranspiration can be fitted. Table 2 shows that the monthly precipitation for every month, with the exception of March, passed the K-S test at the significance level of 0.05, which indicates that the fitting effects of the three distribution functions are good (Figure 7). The AIC value is used for judging the optimal precipitation distribution function. The AIC of gamma distribution or most months is the minimum, so gamma distribution is the optimal distribution function of monthly precipitation. Table 3 shows that the reference evapotranspiration for every month passed the K-S test at the significance level of 0.05, which indicates that the fitting effects of the three distribution functions are good (Figure 8). The AIC of lognormal distribution for most months is the minimum, meaning that the optimal distribution function of monthly reference evapotranspiration is lognormal distribution. AIC values in bold indicate the optimal distribution function according to Akaike information criterion, and * means passed the K-S test at the 0.05 significance level.    For the further analysis of the influence that the coupling effect of precipitation and reference evapotranspiration has on the net irrigation water requirement, a joint distribution model of two related variables must be established. Table 4 displays the parameter fitting results or the joint distribution function of monthly precipitation and reference evapotranspiration based on the Frank copula function. From Figure 9, the probability density map demonstrates a "saddle" shape. When one of precipitation and reference evapotranspiration is large and the other is small, the occurrence probability is high. When both are simultaneously large or small, the occurrence probability is low. It can be seen from Table 4 and Figure 9 that the larger the θ is, the smoother the joint distribution probability density graph is, and the smaller the θ is, the steeper the graph is. This shows that with the decrease of θ, the opposite trend of precipitation and reference evapotranspiration is more obvious. For the further analysis of the influence that the coupling effect of precipitation and reference evapotranspiration has on the net irrigation water requirement, a joint distribution model of two related variables must be established. Table 4 displays the parameter fitting results or the joint distribution function of monthly precipitation and reference evapotranspiration based on the Frank copula function. From Figure 9, the probability density map demonstrates a "saddle" shape. When one of precipitation and reference evapotranspiration is large and the other is small, the occurrence probability is high. When both are simultaneously large or small, the occurrence probability is low. It can be seen from Table 4 and Figure 9 that the larger the θ is, the smoother the joint distribution probability density graph is, and the smaller the θ is, the steeper the graph is. This shows that with the decrease of θ, the opposite trend of precipitation and reference evapotranspiration is more obvious.

Uncertainty Analysis of Net Irrigation Water Requirement
The net irrigation water requirement under different probabilities for each typical year was calculated using the CFMC method as shown in Table 5. Table 5 displays the monthly net irrigation water requirement under different probabilities in each typical year. As the probability of irrigation increases, net irrigation water requirement also increases. From Figure 10, the monthly net irrigation water requirement for each typical year when the probability of irrigation is 75% can be seen. The monthly net irrigation water requirement demonstrates a single peak distribution, and the maximum value that is affected by crop water demand is between May and June. For each typical year, the net irrigation water requirement from March to June demonstrated an increasing trend, whereas the net irrigation water requirement from July to October demonstrated a decreasing trend. From 2000 to 2020, year-to-year, the net irrigation water requirement from

Uncertainty Analysis of Net Irrigation Water Requirement
The net irrigation water requirement under different probabilities for each typical year was calculated using the CFMC method as shown in Table 5. Table 5 displays the monthly net irrigation water requirement under different probabilities in each typical year. As the probability of irrigation increases, net irrigation water requirement also increases. From Figure 10, the monthly net irrigation water requirement for each typical year when the probability of irrigation is 75% can be seen. The monthly net irrigation water requirement demonstrates a single peak distribution, and the maximum value that is affected by crop water demand is between May and June. For each typical year, the net irrigation water requirement from March to June demonstrated an increasing trend, whereas the net irrigation water requirement from July to October demonstrated a decreasing trend. From 2000 to 2020, year-to-year, the net irrigation water requirement from March to June decreased, with the contrary from July to October. The changing trend of net irrigation water requirement is consistent with that of crop water demand, which indicates that crop water demand significantly affects the net irrigation water requirement.

Advantages of the CFMC Model
Although previous research on the joint distribution of precipitation and reference evapotranspiration has been conducted [30,41], these studies failed to quantify the impact that joint distribution has on net irrigation water requirement. The main reason is that

Advantages of the CFMC Model
Although previous research on the joint distribution of precipitation and reference evapotranspiration has been conducted [30,41], these studies failed to quantify the impact that joint distribution has on net irrigation water requirement. The main reason is that although a copula function can establish the coupling relationship between precipitation and reference evapotranspiration, the number of samples limits the statistical analysis of the net irrigation water requirement. However, determining the probability of irrigation water requirement is very important in the design of irrigation. If the probability of irrigation water requirement is determined by hydrological years, the problem of the uneven time distribution of the irrigation water requirement exists. With the aim of overcoming some of the shortcomings of previous research, and by using a Monte Carlo simulation, this paper addressed the problems of insufficient samples and uneven time distribution of the net irrigation water requirements. The CFMC method considers the influence of the joint distribution of precipitation and reference evapotranspiration on net irrigation water requirement but can also calculate the probability of the irrigation water requirement for a typical year of planting structures. The example of the Jingdian Irrigation District demonstrates the application value of the CFMC method, showing that it can provide a reference for the allocation of agricultural water resources and is certainly significant for the efficient utilization of water resources in the irrigation area.

Conclusions
Using a combination of a copula function and the Monte Carlo method, a CFMC method for calculating the probability of the net irrigation water requirement under a typical year of planting structures was proposed in this paper. The research results of the Jingdian Irrigation District demonstrated that: (1) The optimal distribution functions of precipitation and reference evapotranspiration in the irrigation area are gamma and lognormal distributions.
(2) The probability density diagram of the joint distribution function of precipitation and reference evapotranspiration demonstrates a "saddle" shape. When either precipitation or reference evapotranspiration is larger and the other is smaller, the occurrence probability is high, but when both are simultaneously large or small, the occurrence probability is low.
(3) The maximum monthly net irrigation water requirement is between May and June. From 2000 to 2020, the net irrigation water requirement between March and June demonstrated an increasing trend year-to-year, whereas the net irrigation water requirement between July and October demonstrated a decreasing yearly trend. As the probability of irrigation water requirement increases, the net irrigation water requirement in the irrigation area also increases.
(4) Climate uncertainty has a significant impact on the regional irrigation water requirement, and this study provides a quantitative analysis of agricultural water requirements under different probabilities. For the water resources management department of a basin, the results of this study can provide a reference for the assessment of the water shortage risk and economic benefits in the basin.

Institutional Review Board Statement: Not applicable.
Informed Consent Statement: Not applicable.

Data Availability Statement:
The data relevant to this paper can be found in the Supplementary Materials.

Conflicts of Interest:
The authors declare that they have no known competing financial interest.