Next Article in Journal
Use of Pseudomonas protegens to Control Root Rot Disease Caused by Boeremia exigua var. exigua in Industrial Chicory (Cichorium intybus var. sativum Bisch.)
Next Article in Special Issue
Comprehensive Evaluation of Tomato Growth Status under Aerated Drip Irrigation Based on Critical Nitrogen Concentration and Nitrogen Nutrient Diagnosis
Previous Article in Journal
Physiological and Molecular Responses of Pyrus pyraster Seedlings to Salt Treatment Analyzed by miRNA and Cytochrome P450 Gene-Based Markers
Previous Article in Special Issue
Exogenous γ-Aminobutyric Acid Can Improve Seed Germination and Seedling Growth of Two Cotton Cultivars under Salt Stress
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Global Sensitivity Analysis of the Advanced ORYZA-N Model with Different Rice Types and Irrigation Regimes

1
State Key Laboratory of Efficient Utilization of Agricultural Water Resources, China Agricultural University, Beijing 100083, China
2
Institute of Environment and Sustainable Development in Agriculture, Chinese Academy of Agricultural Sciences, Beijing 100081, China
3
Centro de Ciência e Tecnologia do Ambiente e do Mar (MARETEC), Instituto Superior Técnico, Universidade de Lisboa, Av. Rovisco Pais nº. 1, 1049-001 Lisboa, Portugal
*
Authors to whom correspondence should be addressed.
Plants 2024, 13(2), 262; https://doi.org/10.3390/plants13020262
Submission received: 9 November 2023 / Revised: 5 December 2023 / Accepted: 6 December 2023 / Published: 17 January 2024
(This article belongs to the Special Issue Water and Nitrogen Management in Soil-Crop System II)

Abstract

:
Identifying important parameters in crop models is critical for model application. This study conducted a sensitivity analysis of 23 selected parameters of the advanced rice model ORYZA-N using the Extended FAST method. The sensitivity analysis was applied for three rice types (single-season rice in cold regions and double-season rice (early rice and late rice) in subtropical regions) and two irrigation regimes (traditional flood irrigation (TFI) and shallow–wet irrigation (SWI)). This study analyzed the parameter sensitivity of six crop growth outputs at four developmental stages and yields. Furthermore, we compared the variation in parameter sensitivity on model outputs between TFI and SWI scenarios for single-season rice, early rice, and late rice. Results indicated that parameters RGRLMX, FRPAR, and FLV0.5 significantly affected all model outputs and varied over developmental stages. Water stress in paddy fields caused by water-saving irrigation had more pronounced effects on single-season rice than on double-season rice.

1. Introduction

China is among the biggest consumers of nitrogen fertilizer worldwide, and its rice fields account for around 40% of global nitrogen fertilizer use [1]. It is well-known that water and nitrogen are the two major elements that affect rice growth and yield formation in field management. Traditional family-holder agriculture usually applies excessive irrigation and fertilization management approaches for paddy fields in pursuit of high yields, which brings many environmental problems such as soil consolidation, surface source pollution, and greenhouse effects [2,3,4]. Investigating reasonable water and nitrogen management practices in rice fields by means of field trials is often time-consuming and laborious, and the accuracy of measurement results is limited by the experimental conditions [5,6]. In addition, field trials are often influenced by different hydrological year patterns and agrometeorological conditions. Crop models are proposed on the basis of agricultural processes and mathematical methods, which can effectively complement the deficiencies in field experiments [7,8,9]. In recent years, crop models have been widely used in yield prediction, agro-ecological process simulation, field experimental management, and climate change reflections [10,11]. Accordingly, the reliability of predictions is crucial for farmers, researchers, and policy makers.
Crop models (generally developed based on multi-processes) have complicated and nonlinear properties involving numerous parameters in various modules (i.e., soil water, crop growth, management practices, meteorological environment, etc.). Many parameters and their value ranges vary greatly for crop varieties, soil environments, and agronomy management, thus limiting their direct determination through experiments [12,13]. Furthermore, the modules and parameters in crop models are numerous and exist in equifinality. It then becomes difficult to obtain suitable parameter estimation results by a trial-and-error method. Meanwhile, the parameter calibration results are easily influenced by the testers’ subjectivity, resulting in great uncertainty [14].
Parameter sensitivity analysis is usually applied to find out essential parameters in a crop model and greatly reduce the dimensionality of parameter space, which is significantly important for the improvement of model structure and parameter estimation efficiency. Sensitivity analysis (SA) is an effective approach to compute the impacts of unknown factors (parameters or driving variables) on the model’s output variables in a qualitative or quantitative manner [15,16,17]. SA methods are generally classified into two types: local and global SA techniques. Global SA is widely used in crop models, and the most representative methods are the Morris screening method, FAST (Fourier Amplitude Sensitivity Test), Sobol, and Extended FAST (EFAST) based on the variance method [18,19,20,21]. Extended FAST is an outstanding global SA method that has been adapted to many crop models, such as WOFOST [7], EPIC [22], AquaCrop [5], WARM [10], and APSIM [23].
Model parameters usually have different sensitivities under different agro-meteorological conditions [8,24]. Previous studies show that key parameters may be different between single-season rice and double-season rice (i.e., early rice and late rice) in cold regions and tropical/subtropical regions due to their various growth environments (e.g., cultivars, temperature, light, and radiation conditions) [11,25,26]. In the process of model parameter calibration, it is necessary to analyze the parameter sensitivity to different crop growth environments in order to produce more informative results.
The ORYZA model, which is extensively used worldwide, has demonstrated its reliability in simulating rice growth and paddy water balance across various experimental sites and cultivars [11,26,27,28]. Due to insufficient evidence and uncertainties in field experiments, the coupled effects of water and nitrogen in paddy fields have remained controversial over recent years [15,29]. In addition, the sensitivity of nitrogen-related parameters in crop models is still unknown. However, the validation of water and nitrogen coupling effects in paddy rice fields is crucial for simulating rice growth and development, formulating reasonable water and nitrogen management measures, and improving water and nitrogen use efficiency in paddy fields.
The ORYZA-N model (proposed by Gao et al. (2023) [30]) is an advanced soil–water–nitrogen–plant coupling model for paddy rice fields, which particularly strengthens the simulation capability of nitrogen balance and proves its accuracy and reliability. The stability and universal applicability of this advanced model still deserve to be further explored, especially for its use in different climatic zones [30,31,32]. In this study, we applied the Simlab software (ver. 2.2) to implement the Extended FAST method and achieve the following objectives: (1) to improve the application of ORYZA-N model to rice in typical climate regions (i.e., the single-season rice in cold regions and the double-season rice in subtropical regions); (2) based on the fully coupled water and nitrogen characteristics of the ORYZA-N model, to clarify the sensitive parameters on biomass and nitrogen related outputs at different development stages; and (3) to test the model stability and applicability in rice growth, yield formation, and nitrogen uptake simulation in typical climate regions. This study aims to simplify the input parameterization process for the ORYZA-N model and provide important calibration and evaluation suggestions for future model applications to rice cultivars in different climate regions.

2. Results and Discussion

2.1. Sensitivity Indices of Parameters for Model Crop-Related Outputs at Different Development Stages

2.1.1. Basic Vegetative Phase

The sensitivity indices of 23 model parameters for crop-related outputs WAGT (a–c), WST (d–f), and LAI (g–i) at the end of the basic vegetative phase (DVS = 0.4) for early rice, late rice, and single-season rice under TFI and SWI regime were presented in Figure 1. At this stage, the parameter sensitivities of the same output variables for each rice type hardly differed under different irrigation regimes. However, the sensitivity effects of the 23 parameters varied a lot for different output variables of different rice types. The SA results indicated that the developmental rate in the juvenile phase (DVRJ), the sensitivity effects of the photosynthetically active sunlight energy radiation fraction parameter (FRPAR), and the maximum relative growth rate of leaf area (RGRLMX) were the most sensitive parameters affecting almost all crop-related outputs, with the exception of LAI of early rice, which was insensitive to the FRPAR parameter. The phenological development rate determines the rate at which the crop utilizes the average daily ambient temperature and photoperiod to generate assimilates [26]; thus, the parameter DVRJ was one of the most sensitive parameters at the beginning of rice growth. Due to the variation in early growth ambient temperatures, with single-season rice having the longest growth period, followed by early rice and late rice, the effect of the radiation coefficient FRPAR and its accumulation during the growing period could not be neglected for each rice type.
RGRLMX was the most sensitive parameter for early rice and late rice, but for single-season rice, it became the second most sensitive parameter, replaced by FRPAR at the top. FLV0.5 (at the DVS = 0.5, the fraction of shoot dry matter allocated to leaves) was found to be a sensitive parameter for WST in early and late rice, and LAI in single-season rice. Its effect on WAGT was negligible. For the output variable LAI, the sensitivity results were more similar for late rice and single-season rice, although they had different values for the SA results. However, there were significant differences between early rice and late rice, with LAI for early rice being sensitive only to DVRJ and RGRLMX. Generally, the results of this phenological phase were consistent with Tan et al. (2016) [8]. Due to the leaf morphology, higher latitude, and lower plant density [33,34,35], single-season rice in cold regions was more sensitive to FRPAR than double-season rice in subtropical regions, with the total-order sensitivity indices reaching over 0.4. It should be noted that in this model, the allocation coefficients of organs between two adjacent developmental stages were varied linearly. The sensitivity results of RGRLMX on LAI, WAGT, and WST were consistent with the results drawn from the ORYZA (v3) model [36], while results for other crop parameters were inconsistent. This was due to the different model structures of ORYZA-N and ORYZA (v3), with the former model providing a comprehensive improvement of the nitrogen module and other modules.

2.1.2. Photoperiod-Sensitive Phase

At this stage, the biomass of panicles (WSO) began to appear. The results of SA for TFI and SWI were not significantly different for various rice varieties (Figure 2). Output results for early rice and late rice were similar for all biomass components (i.e., WAGT, WST, LAI, and WSO), and all outputs were sensitive to FRPAR. Some other outputs had inconsistent sensitivity results. For example, WAGT, LAI, and WSO were sensitive to FLV0.5, but WST was insensitive; WAGT and WST were sensitive to the phenological development rate parameter during the juvenile phase and photoperiod-sensitive phase (i.e., DVRJ and DVRI) and RGRLMX, but their effects were negligible for LAI and WSO. FST0.75 and FLV0.75 (the coefficients of shoot dry matter allocated to stems and leaves at the DVS = 0.75) only influenced WSO because the fraction of dry matter allocated to the panicles after panicle formation could be obtained by minus the portion of leaves and stems from 1.
For single-season rice, FRPAR was the most sensitive parameter affecting all crop-related outputs, ranging from 0.43 to 0.54, because of the higher latitude and light radiation. FLV0.5 was the second most sensitive parameter to LAI and WAGT, the fourth most sensitive to WSO, and insensitive to WST. WSO of single rice was also influenced by FST0.75 and FLV0.75, agreeing with WSO of double-season rice. The sensitivity observed in all crop-related outputs to RGRLMX and FLV0.5 agreed with Tan et al. (2020) [36] and Yu et al. (2023) [31], while for other parameters, results were inconsistent due to the same reasoning as in the former phase. It was worth noting that the radiation parameter (FRPAR) became the most sensitive parameter, affecting all outputs from this stage. Radiation had a positive effect on rice growth and development, and high solar radiation interception by the rice canopy was closely related to high biomass production and grain yield [37,38]. In agreement with Tu et al. (2022) [39], radiation was one of the most influential factors in rice growth.

2.1.3. Panicle-Formation Phase

At the panicle-formation phase, different irrigation regimes slightly altered the sensitivity values of some parameters, such as DVRJ on all outputs of early and late rice, but did not significantly affect the determination of parameter sensitivity levels (Figure 3). FRPAR was the most sensitive parameter affecting almost every crop output, except LAI for double-season rice, which was second to FLV0.5. Most outputs were sensitive to FLV0.5, except WST and WSO of single-season rice and early rice (Figure 3). Because single-season rice and early rice had similar environmental growth temperatures during the early period, they exhibited similar sensitivity results.
Additionally, the difference in the level of parameter sensitivity resulting from the types of rice was present only in those parameters that had an equal effect on the variation of output, which showed a poor rule of variation. With the exception of the radiation parameters, WSO was strongly influenced by the dry matter partition coefficients of the organs within the stage. Since shallow–wet irrigation did not cause severe water stress, the sensitivity results of the different irrigation regimes were not significantly different until the current phase, and the crop-related output sensitivity was almost the same under both irrigation regimes. The sensitivity results were consistent with the calibration and validation procedures in our previous papers [11,30].

2.1.4. Grain-Filling Phase

All outputs were sensitive to parameter DVRJ during the vegetative growth period (before panicles form) but not sensitive during the reproductive growth period (after panicles form). Compared with previous stages, the effect of rice type on the sensitivity results for each section of biomass hardly changed and was mainly concentrated on single-season and double-season rice, whereas the sensitivity levels varied a little for double-season rice itself (Figure 4). The difference in sensitivity between single-season and double-season rice should be attributed to the varying ambient temperatures and climates of the planting areas [23]. FRPAR and FLV0.5 were the most sensitive and second most sensitive parameters influencing all outputs, with WST also affected by FST0.75. LAI differed from other outputs for both single-season and double-season rice, and was sensitive to multiple parameters, such as FRPAR, FLV0.5, DVRJ, the fraction of leaf death coefficient at the DVS = 1.0 and 1.6 (DRLV1.0 and DRLV1.6), the phenological development rate in reproductive phase DVRR, and FLV0.75 (Figure 4g–i). LAI at rice maturity was sensitive to leaf death rate (DRLV1.0, DRLV1.6) due to the fact that LAI was closely related to leaf biomass and decreased when the senescence death rate of leaves accelerated. This was consistent with the findings of Liu et al. (2019) [23]. At the same time, leaves that could photosynthesize were reduced, and assimilated dry matter mass was consequently decreased so that the panicle biomass was also affected. The sensitivity results were consistent with Liu et al. (2019) [23] but differed from Tan et al. (2016 and 2020) [8,36], where the effects of DRLV1.0 and DRLV1.6 were null. This was due to the difference in temperature and climate of the growing environment.
The effects of different irrigation regimes (i.e., TFI and SWI) on the sensitivity of model output variables were mainly observed in all model outputs for single-season rice. For double-season rice, the irrigation regime mainly affected LAI. LAI and WSO for single-season rice were sensitive to multiple parameters, and the interaction among the parameters had a significant effect on parameter sensitivity, indicating that for single-season rice, SWI could cause a degree of water stress by rapidly drying out the ponded water layer during rice maturity period, which was detrimental to leaf growth and rice yield formation (Figure 4i). In addition, the interaction effect between parameters increased, especially for LAI and WSO in single-season rice (Figure 4i,l). Since the ponded water layer in cold regions has a thermal insulation effect on crops, single-season rice in cold regions was more sensitive to irrigation regimes (i.e., the water stress caused by irrigation amount) than double-season rice [11,40]. The interactions among parameters were stronger for single-season rice under water-saving irrigation conditions (SWI) since the crop experienced greater water and nitrogen stress, similar to other modeling results of Liang et al. (2017) [15] and Wang et al. (2021) [41].

2.2. Sensitivity Indices of Parameters for Model Nitrogen-Related Outputs at Different Development Stages

2.2.1. Basic Vegetative Phase

At the basic vegetative stage, the different irrigation regimes barely affected the sensitivity levels of the nitrogen-related outputs ANCR and ANLV (Figure 5). The sensitivity results between ANCR and ANLV for different rice types were also in good agreement. Nevertheless, there was no pattern in the sensitivity of these outputs to the selected parameters. For the same output, the impact parameters under different climatic conditions were roughly the same, but the order tended to be different, which was consistent with the findings of Liu et al. (2019) [23]. For example, DVRJ was the most sensitive parameter affecting ANCR and ANLV in early rice, the second sensitive parameter in late rice, as well as the third/fourth sensitive parameter in single-season rice. RGRLMX was the most sensitive parameter for late rice, while FRPAR had the greatest effect on single-season rice (Figure 5). FLV0.5 should also be considered when calibrating parameter values for different rice types. Note that ANCR and ANLV of late rice were sensitive to the coefficient of solute uptake from soil by rice roots at DVS = 0 and DVS = 0.5 (KR0 and KR0.5). This was because nitrogen content within rice affects photosynthetic assimilation capacity, and a decrease in nitrogen content was accompanied by a decline in the ability to capture light [26,42]. The outputs ANCR and ANLV were more sensitive to KR0 and KR0.5 in late rice, where the ambient temperature was suitable during the seedling stage of rice to meet the demands of rapid rice growth. In general, the sensitivity levels of these influential parameters were the same, although some sensitivity indices varied slightly.

2.2.2. Photoperiod-Sensitive Phase

The growth stage of DVS = 0.4–0.65 was the rapid growth period of rice stems and leaves development, and the whole parameter sensitivity was different from the previous stage, but the parameters that had impacts on the outputs were approximately the same, namely FRPAR, FLV0.5, RGRLMX, DVRI, and DVRJ (Figure 6). At this stage, for the selected model outputs ANCR and ANLV, the sensitivity results were relatively similar for early rice and late rice, but the parameter values for the SA results were different (Figure 6a–d). For example, for double-season rice, FRPAR was the most sensitive parameter affecting the ANCR of early rice, while FLV0.5 was the most sensitive parameter for other outputs. The SA values of ANCR and ANLV for single-season rice were also different (Figure 6c,f).
In addition, the sensitivity results were approximately the same for different irrigation regimes, while the sensitivity parameters differed for different rice types (mainly between single-season and double-season rice). The nitrogen-related outputs (ANCR and ANLV) of single-season rice were sensitive to FLV0.5, FRPAR, and RGRLMX, which was consistent with previous ORYZA-N model simulations using field trial data [11]. It was noteworthy that ANCR and ANLV of single-season rice under the SWI condition were also sensitive to the solute (i.e., NH4-N and NO3-N) diffusion coefficient below the soil surface (Kdif2). Nitrogen uptake and utilization in rice under water-saving irrigation was higher than that under flood irrigation [43], and thus, the influence of nitrogen diffusion on crop N uptake was also expanded.

2.2.3. Panicle-Formation Phase

During panicle formation, the sensitivity of double-season rice under different irrigation regimes was affected by paddy water conditions, whereas the sensitivity for single-season rice was irregular. Under TFI conditions, ANCR was mainly affected by the variation of KR0.5, FLV0.5, and FRPAR, although the three parameters had different magnitudes of sensitivity. KR0.5 was the most sensitive parameter of ANCR for early rice, the second sensitive parameter for late rice, and the third sensitive parameter for single-season rice (Figure 7a–c).
ANLV of double-season rice was influenced by FLV0.5, FST0.75, and FRPAR, which should not be neglected, and the sensitivity indices of the most sensitive parameter, FLV0.5, were about 0.5. While the output of single-season rice mainly depended on the distribution rate, death rate, and photosynthetic effective utilization of leaves, the sensitivity results of ANLV were similar to LAI. Due to the difference in climate types between double- and single-season rice growing regions, the effect of paddy water conditions (caused by different irrigation regimes) was more noticeable for single-season rice.
The sensitivity pattern difference between ANCR and ANLV is not limited to single-season rice. Sensitivity patterns changed across all rice types. For early rice, the major flip occurred at parameters KR05, FLV05, and FLV075; meanwhile, the changed parameters were KR05, KDif2, FLV05, and FRPAR for single-season rice. SWI had a large effect on the model nitrogen outputs sensitivity to parameters, especially on single-season rice, and comparison of Figure 7c,f showed that ANCR and ANLV of single-season rice at SWI were affected by a combination of multiple parameters, the most influential of which were Kdif2 and FLV0.5, respectively. Similar to the previous stage, due to the water stress caused by SWI during this period, when conducting model calibration and validation under controlled irrigation regimes, the effect of the solute diffusion below the soil surface coefficient Kdif2 on ANCR must be taken into consideration.

2.2.4. Grain-Filling Phase

Different irrigation regimes had significant effects on N-related outputs for different rice types. At this phase, the sensitivity parameters, sensitivity change patterns, and sensitivity change magnitude were quite different on ANCR and ANLV of each rice type, compared with the previous periods. ANCR was sensitive to KR0.5, FLV0.5, and FRPAR in both single- and double-season rice, and the sensitivity values of the SA results were different (Figure 8). For example, KR0.5 was the most sensitive parameter for early rice (STi > 0.5), the second most sensitive parameter for late rice (STi was about 0.3), and the third most sensitive parameter for single-season rice (STi > 0.1), which was different from the sensitivity values at the panicle-formation stage (Figure 8). Under TFI, for early rice, ANCR was influenced by KR0.5 and FLV0.5, while for late rice and single-season rice, ANCR was affected by KR0.5, FLV0.5, and FRPAR.
The interaction among parameters under SWI was obvious; the ANCR and ANLV of all rice types were affected by a combination of multiple parameters, with the main sensitive parameters of double-season rice being KR0.5, FLV0.5, DVRJ, FRPAR, and RGRLMX (Figure 8a–c). This indicated that the interaction effects of nitrogen-related outputs were greater under water stress conditions, where priority should be provided to the soil water-related parameters [44]. The sensitivity results of ANCR and ANLV for single-season rice were both mainly influenced by Kdif2 and FRPAR. When focusing on the model output variable ANLV for early rice and single-season rice, the effects of DRLV1.0 and DRLV1.6 should be considered during parameter calibration. Overall, the parameter calibration and application of the nitrogen-related output should focus on the nitrogen parameters KR0.5 and Kdif2, the crop parameter FLV0.5, and the radiation parameter FRPAR.

2.3. Sensitivity Indices of Parameters for Yield at Maturity

Yields are usually the most focused on component of the crop models. As a result, for WRR14 (yield with 14% moisture at maturity), three rice types and two irrigation regimes were simulated, and the Si and STi were computed for the 23 parameters (Figure 9). The results were comparable to WSO at the maturity phase (as shown in Figure 4j–l), but the differences due to rice type (i.e., growth ambient environment) were more noticeable. It was important to note that the ORYZA-N model’s WRR14 was computed using the grains number growth rate, not just the WSO at ripe.
For WRR14, FRPAR was the most important parameter, followed by FLV0.5, FLV0.75, and DVRR. For both TFI and SWI of double-season rice, the effects of all other parameters were negligible. As for single-season rice, the sensitivity of parameters under TFI was consistent with that of double-season rice, and the effect of Kdif2 on the sensitivity results was not ignorable in SWI. The developmental rate of each stage (DVRJ, DVRI, DVRP, DVRR) had a varying impact on dry biomass and nitrogen content of each organ, which also indirectly influenced the yield formation. This was consistent with the results of Yu et al. (2023) [31] and Xu et al. (2018) [38]. Since the ponded water layer of single-season rice in cold regions has the role of maintaining soil temperature in the root zone and sustaining crop growth during the growing season, controlled irrigation regimes would affect the growth and yield formation of single-season rice in cold regions at some degree, and the interaction effect among the parameters was higher than other scenarios [11,39].

3. Materials and Methods

In this study, we performed a parameter sensitivity analysis (SA) of the ORYZA-N model using the Extended FAST method. The ORYZA model employed the development stages (DVS) of rice growth to determine the physiological phase of the rice plant, which is characterized by the emergency of the different organs and their appearance [26]. The whole rice growth process was separated into four phenological phases: the basic vegetative phase (DVS = 0–0.4), the photoperiod-sensitive phase (DVS = 0.4–0.65), the panicle formation phase (DVS = 0.65–1.0), and the grain-filling phase (DVS = 1.0–2.0) [26]. To assess the sensitivity of model output variables to changes in the selected input parameters, we used the Simlab software [45] (https://ec.europa.eu/jrc/en/samo/simlab (accessed on 10 April 2022)) to calculate the first-order sensitivity indices and the total-order sensitivity indices of each parameter at the end of every phenological phase.

3.1. ORYZA-N Model

The ORYZA series model is a dynamic and mechanistic eco-physiological model that was jointly developed by International Rice Research Institution and Wageningen University. The model is designed to simulate rice growth, paddy water balance, and field management and is widely used in various regions around the world and for different cultivars [26,28,46,47]. The crop nitrogen uptake method in ORYZA2000 model is a simple bookkeeping algorithm with a daily given rate of crop nitrogen uptake, which is insufficient to quantify the nitrogen balance of the soil-crop system. To simulate and quantify the water and nitrogen balance items in paddy rice effectively, the ORYZA-N model was improved based on the ORYZA2000 version. Three new modules (i.e., a soil temperature module, an ammonium nitrogen (NH4-N) and nitrate nitrogen (NO3-N) transport module, and a soil organic carbon and nitrogen turnover module) were developed and incorporated into the ORYZA-N model. The soil water module was also improved by enhancing the calculation procedures for root water uptake and soil water drainage. ORYZA-N model was calibrated and validated using two years of field experiment data, and the results showed its reliability in the aspect of simulating the agro-hydrological processes, crop growth, and especially water–nitrogen transformation in paddy rice fields [30]. More detailed information on the ORYZA-N model can be found in Gao et al. (2023) [30].

3.2. Extended FAST Method

This study applied the Extended FAST method, a variance-based global SA approach, to evaluate the impact of diverse parameters from different modules on the model outputs, including dry biomass, crop nitrogen content, leaf area index (LAI), and yield, and calculate higher-order interactions among the parameters [17,48]. Extended FAST was developed on the basis of the FAST algorithm (Fourier Amplitude Sensitivity Test [49]) and the Sobol algorithm [50]. It combines the computing efficiency of the FAST method with Sobol’s ability to calculate the total effect [5,51]. The Extended FAST method used in this study has the advantages of high efficiency and accuracy and is capable of analyzing the interaction effects between parameters. Extended FAST has gained popularity in various fields, such as hydrological, ecological, and crop modeling, as well as climate change studies, owing to the advantages mentioned above [5,10,12].
The total variance V(Y) of model outputs was calculated as follows [52]:
V ( Y ) = i = 1 n V i + i < j n n V i j + + V 12 n
where V(Y) indicated the total variance of model output Y induced by ORYZA-N model, Vi ( V i = V [ E ( Y | x i ) ] ) represented the first order for each factor, and Vij ( V i j = V [ E ( Y | x i ,   x j ) ] V i V j ) to V12…n represented the interactions among n factors. The variance of the conditional expectation V [ E ( Y | x i ) ] , also known as the main effect, was used to determine the significance of xi on the variance of Y, indicating the sensitivity of Y to xi. Extended FAST calculated two indices: the first-order sensitivity indices (Si) and total-order sensitivity indices (STi). The Si corresponding to a single factor (xi) was calculated as follows:
S i = V i / V ( Y ) = V [ ( E ( Y / x i ) ] / V ( Y )
The STi that correspond to a single factor (index i) and the interaction of more factors that involve the index i and at least ji from 1 to n was calculated by Equation (3):
S T i = S i + j i S i j + + S 1 i n = V ( Y ) V i V ( Y )
Higher values denoted a greater contribution of the variable on the outcomes and vice versa.

3.3. Experiment and Data

Field experiment for single-season rice was carried out in 2019 at the Qing’an Irrigation Experimental Station (46°58′ N, 127°39′ E), located in the typical cold region (Songnen Plain) of northeast China. The region is characterized by a middle temperate and sub-humid continental monsoon climate, with an average annual temperature of 2.5 °C and an average annual precipitation of 635 mm. Field experimental data for double-season rice of 2012 was collected from the Jiangxi Irrigation Experiment Station in Nanchang, China (located at 28°26′ N, 115°28′ E), with a warm and humid subtropical monsoon climate. The average annual temperature and precipitation of this region were 18.1 °C and 1636 mm, respectively.
The main soil texture of both study sites is silty loam, with the main soil physical characteristics listed in Table 1. The experimental data and the phenology investigation dates used in this work were obtained in 2019 and 2012 for single-season rice in a cold region and double-season rice in a subtropical region, respectively. According to the rainfall pattern, the selected experimental years for both single- and double-season rice were wet years so that differences in the irrigation regime sensitivity caused by the climatic conditions of the experimental years (mainly rainfall) could be neglected. The dates associated with phenology, such as sowing, emergency, transplanting, panicle initiation, blooming, and maturity, were observed and recorded (Table 2). In addition, daily meteorological data, including minimum and maximum temperature, sunshine hours, relative humidity, and wind speed during the experimental years, were obtained from the nearby Tieli and Jiangxi weather stations. These data were used to calculate the reference evapotranspiration (ETo) according to the FAO Penman–Monteith method [53]. Nitrogen fertilizer application schedules for the single-season rice, early rice, and late rice during the selected experimental years are provided in Table 3 and were defined following local farmers’ practices. More detailed information regarding observations in single-season rice and double-season rice can be found in Gao et al. (2021 and 2023) [11,30] and Tan et al. (2016) [8], respectively.

3.4. Model Parameters and Sensitivity Analysis

In this work, 23 parameters in soil water, crop, and nitrogen modules were considered for SA (Table 4). The output state variables selected for the ORYZA-N model were dry aboveground green biomass (WAGT, kg ha−1), the leaf area index (LAI, m2 m−2), the dry biomass of stems (WST, kg ha−1) and panicles (WSO, kg ha−1), the nitrogen content in dry aboveground green biomass (ANCR, kg N ha−1) and dry biomass of leaves (ANLV, kg N ha−1). These outputs were analyzed at the end of the four phenological phases (i.e., DVS = 0.4, 0.65, 1.0, and 2.0).
First, the parameter values of model “warm-up” were determined based on the recommended/default parameter values of single-season rice in cold regions and double-season rice in subtropical regions from Gao et al. (2021 and 2023) [11,30] and Tan et al. (2016 and 2020) [8,36], respectively. Then, based on these parameter values and field experiment data, a rational calibration of the crop growth process was carried out for different rice types. The calibration was conducted using a trial-and-error approach, particularly aiming to provide the basic reference parameter values for SA (the comparison results are shown in Supplementary Figure S1). Both simulated and measured values showed good agreement; thus, the parameter values used in these cases could be applied as base values for SA. Next, 23 parameters related to water stress, crop growth, and nitrogen uptake/transformation were selected according to the calibration process and set to the base values obtained in the previous step. Detailed SA was finally conducted to identify the sensitive parameters for three rice types at four development stages.
Table 4 lists the base values and variance range of the selected parameters in the ORYZA-N model. The variation range of selected parameters was set as ±30% of base values as recommended by Tan et al. (2017) [29], while for four phenological parameters, the variance range of three rice types was set as ±10% because of the low ambient temperatures in the early stages of crop growth and a short suitable sowing period.
To evaluate the influence of different irrigation regimes on the parameter sensitivity results, traditional flood irrigation (TFI) and shallow–wet irrigation (SWI) regimes were considered in simulations of single-season rice, early rice, and late rice. TFI and SWI regimes were the most popular irrigation regimes in both experimental areas. The details about irrigation management at each growth stage are presented in Table 5.
In this study, the sensitivity indices Si and STi were calculated by the Simlab (ver. 2.2) software using the Extended FAST method, as described by Tarantola [54]. Based on the knowledge of the convergence of sensitivity indices, a parameter sample size of N = 15,000 was set for the SA. This method included the following four steps: (1) Definition of parameters and their distributions. There were 6 soil water-related parameters, 12 crop-related parameters, and 5 nitrogen-related parameters. The default values and ranges are listed in Table 4. (2) Generation of sets of targeted parameters by terms of Monte Carlo distributions sampling. (3) Evaluation of model outputs for each set of selected input parameters. The ORYZA-N model was called in this process, and the results of the runs with different parameter combinations generated by Simlab were recorded through the Python code. (4) Computation of Si and STi from the outputs resulting from variations of each input parameter produced by Simlab.

4. Conclusions

In this study, the Extended FAST sensitivity analysis (SA) method was used to analyze the parameter sensitivities for the ORYZA-N model in different rice growth phases. The SA was carried out for three rice types (single-season rice, early rice, and late rice) and two irrigation regimes (traditional flood irrigation (TFI) and shallow–wet irrigation (SWI)) in two typical climatic conditions (subtropical and cold temperate regions). Twenty-three parameters related to the soil-water module, crop module, and nitrogen module were selected from the ORYZA-N model to identify the sensitivity effects on six model outputs of crop growth (aboveground biomass, WAGT; dry biomass of stems and panicles, WST and WSO; leaf area index, LAI; nitrogen content in dry aboveground biomass and dry biomass of leaves, ANCR and ANLV).
In the ORYZA-N, there were only a few parameters that had significant effects on the model outputs, which could simplify the calibration procedure and provide guidance during the model application. Different rice types showed variations in the sensitivity degrees. The results showed that RGRLMX, FLV0.5, and FRPAR were the most sensitive parameters for both crop-related and nitrogen-related outputs. The sensitivity ranking of the above parameters varied among developmental stages, which should be particularly paid attention. For the three rice types, the parameter sensitivities to most outputs were similar for early rice in sub-tropical regions and single-season rice in cold regions but different from late rice. Throughout the growing periods, DVRJ, RGRLMX, and DVRI were the most sensitive parameters for biomass and nitrogen uptake during the vegetative growth period. During the reproductive growth period, FRPAR and FLV0.5 had a more prominent effect on all outputs. ANCR was sensitive to KR0.5, while ANLV was closely correlated to FLV0.75 and leaf death rate coefficients. For double-season rice (early rice and late rice), the sensitivity performance was similar under TFI and SWI regimes; meanwhile, we found that the interaction effects were stronger for water deficit conditions. For single-season rice, all outputs were sensitive to FRPAR, with stronger and more sensitive interactions than for double-season rice under water deficit conditions. For all rice types and irrigation regimes, FRPAR was the most sensitive parameter for yield (i.e., the output WRR14), followed by FLV0.5. The effect of the parameter Kdif2 on nitrogen and yield outputs could not be negligible under SWI for single-season rice.
Results also indicated that ORYZA-N has good computational and structural stability. The fully coupled water–nitrogen structure also verified the significance of the nitrogen parameter simulation process for rice modeling, which was valuable to the improvement of nitrogen simulation for other crop models. The base parameter set in this study could also complement the parameter database for rice modeling and provide useful information for the calibration and application of the advanced ORYZA-N model in different climatic regions.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/2223-7747/13/2/262/s1. Figure S1: Model calibration comparisons of simulated and measured crop growth indicators for single-season rice (a,d), early rice (b,e), and late rice (c,f) involving dry aboveground green biomass (WAGT), stem biomass (WST), green leaves (WLVG) and panicles (WSO), and leaf area index (LAI). (This figure is adapted from Gao et al., 2023 [30]).

Author Contributions

Conceptualization, Y.G. and X.X.; methodology, J.T.; software, Y.G.; validation, Y.G. and Q.H.; formal analysis and investigation, C.S.; resources and project administration, G.H.; data curation, C.S. and J.T.; writing—original draft preparation, Y.G. and A.R.O.; writing—review and editing, T.B.R., J.T. and X.X.; supervision, X.X.; funding acquisition, X.X., C.S. and T.B.R. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the 14th Five-year National Key Research and Development Program of the Chinese Ministry of Science and Technology: 2023YFD2300301, National Natural Science Foundation of China: 52022108, National Natural Science Foundation of China: 52379053, 13th Five-year National Key Research and Development Program of the Chinese Ministry of Science and Technology: 2017YFD0300403, Fundação para a Ciência e Tecnologia: CEECIND/01152/2017.

Data Availability Statement

Data are contained within the article and Supplementary Materials.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. FAO. FAOSTAT (Statistical Database of the Food and Agricultural Organization of the United Nations). Available online: http://faostat.fao.org/site/405/default.aspx (accessed on 5 February 2023).
  2. Ahmad, Z.U.; Sanin, M.; Lian, Q.Y.; Zappi, M.; Gang, D.D. Nonpoint Source Pollution. Water Environ. Res. 2017, 89, 1580–1602. [Google Scholar] [CrossRef]
  3. Guo, Y.; Chen, Y.; Searchinger, T.D.; Zhou, M.; Pan, D.; Yang, J.; Wu, L.; Cui, Z.; Zhang, W.; Zhang, F.; et al. Air quality, nitrogen use efficiency and food security in China are improved by cost-effective agricultural nitrogen management. Nat. Food 2020, 1, 648–658. [Google Scholar] [CrossRef] [PubMed]
  4. Li, Y.; Huang, G.; Chen, Z.; Xiong, Y.; Huang, Q.; Xu, X.; Huo, Z. Effects of irrigation and fertilization on grain yield, water and nitrogen dynamics and their use efficiency of spring wheat farmland in an arid agricultural watershed of Northwest China. Agric. Water Manag. 2022, 260, 107277. [Google Scholar] [CrossRef]
  5. Jin, X.; Li, Z.; Nie, C.; Xu, X.; Feng, H.; Guo, W.; Wang, J. Parameter sensitivity analysis of the AquaCrop model based on extended fourier amplitude sensitivity under different agro-meteorological conditions and application. Field Crops Res. 2018, 226, 1–15. [Google Scholar] [CrossRef]
  6. Ren, D.; Xu, X.; Engel, B.; Huang, G. Growth responses of crops and natural vegetation to irrigation and water table changes in an agro-ecosystem of Hetao, upper Yellow River basin: Scenario analysis on maize, sunflower, watermelon and tamarisk. Agric. Water Manag. 2018, 199, 93–104. [Google Scholar] [CrossRef]
  7. Gilardelli, C.; Confalonieri, R.; Cappelli, G.A.; Bellocchi, G. Sensitivity of WOFOST-based modelling solutions to crop parameters under climate change. Ecol. Model. 2018, 368, 1–14. [Google Scholar] [CrossRef]
  8. Tan, J.; Cui, Y.; Luo, Y. Global sensitivity analysis of outputs over rice-growth process in ORYZA model. Environ. Model. Softw. 2016, 83, 36–46. [Google Scholar] [CrossRef]
  9. Steduto, P.; Hsiao, T.C.; Raes, D.; Fereres, E. AquaCrop—The FAO Crop Model to Simulate Yield Response to Water: I. Concepts and Underlying Principles. Agron. J. 2009, 101, 426–437. [Google Scholar] [CrossRef]
  10. Confalonieri, R.; Bellocchi, G.; Bregaglio, S.; Donatelli, M.; Acutis, M. Comparison of sensitivity analysis techniques: A case study with the rice model WARM. Ecol. Model. 2010, 221, 1897–1906. [Google Scholar] [CrossRef]
  11. Gao, Y.; Xu, X.; Sun, C.; Ding, S.; Huo, Z.; Huang, G. Parameterization and modeling of paddy rice (Oryza sativa L. ssp. japonica) growth and water use in cold regions: Yield and water-saving analysis. Agric. Water Manag. 2021, 250, 106864. [Google Scholar] [CrossRef]
  12. Wang, J.; Li, X.; Lu, L.; Fang, F. Parameter sensitivity analysis of crop growth models based on the extended Fourier Amplitude Sensitivity Test method. Environ. Model. Softw. 2013, 48, 171–182. [Google Scholar] [CrossRef]
  13. Vazquez-Cruz, M.A.; Guzman-Cruz, R.; Lopez-Cruz, I.L.; Cornejo-Perez, O.; Torres-Pacheco, I.; Guevara-Gonzalez, R.G. Global sensitivity analysis by means of EFAST and Sobol’ methods and calibration of reduced state-variable TOMGRO model using genetic algorithms. Comput. Electron. Agric. 2014, 100, 1–12. [Google Scholar] [CrossRef]
  14. He, D.; Wang, E.; Wang, J.; Robertson, M.J. Data requirement for effective calibration of process-based crop models. Agric. For. Meteorol. 2017, 234–235, 136–148. [Google Scholar] [CrossRef]
  15. Liang, H.; Qi, Z.; DeJonge, K.C.; Hu, K.; Li, B. Global sensitivity and uncertainty analysis of nitrate leaching and crop yield simulation under different water and nitrogen management practices. Comput. Electron. Agric. 2017, 142, 201–210. [Google Scholar] [CrossRef]
  16. Saltelli, A.; Tarantola, S.; Campolongo, F.; Ratto, M. Sensitivity Analysis in Practice; Joint Research Centre of the European Commission: Ispra, Italy, 2004. [Google Scholar]
  17. Saltelli, A.; Tarantola, S.; Chan, K.P.S. A Quantitative Model-Independent Method for Global Sensitivity Analysis of Model Output. Technometrics 1999, 41, 39–56. [Google Scholar] [CrossRef]
  18. Ge, Q.; Menendez, M. Extending Morris method for qualitative global sensitivity analysis of models with dependent inputs. Reliab. Eng. Syst. Saf. 2017, 162, 28–39. [Google Scholar] [CrossRef]
  19. Kucherenko, S.; Klymenko, O.V.; Shah, N. Sobol’ indices for problems defined in non-rectangular domains. Reliab. Eng. Syst. Saf. 2017, 167, 218–231. [Google Scholar] [CrossRef]
  20. Morris, M.D. Factorial Sampling Plans for Preliminary Computational Experiments. Technometrics 1991, 33, 161–174. [Google Scholar] [CrossRef]
  21. Sobol, I.M. Sensitivity Estimates for Nonlinear Mathematical Models. Math. Model. Comput. Exp. 1993, 4, 407–414. [Google Scholar]
  22. Zhang, X.; Guo, H.; Wang, R.; Lin, D.; Gao, Y.; Lian, F.; Wang, J.A. Identification of the Most Sensitive Parameters of Winter Wheat on a Global Scale for Use in the EPIC Model. Agron. J. 2017, 109, 58–70. [Google Scholar] [CrossRef]
  23. Liu, J.; Liu, Z.; Zhu, A.-X.; Shen, F.; Lei, Q.; Duan, Z. Global sensitivity analysis of the APSIM-Oryza rice growth model under different environmental conditions. Sci. Total Environ. 2019, 651, 953–968. [Google Scholar] [CrossRef] [PubMed]
  24. DeJonge, K.C.; Ascough, J.C.; Ahmadi, M.; Andales, A.A.; Arabi, M. Global sensitivity and uncertainty analysis of a dynamic agroecosystem model under different irrigation treatments. Ecol. Model. 2012, 231, 113–125. [Google Scholar] [CrossRef]
  25. Soundharajan, B.; Sudheer, K.P. Sensitivity analysis and auto-calibration of ORYZA2000 using simulation-optimization framework. Paddy Water Environ. 2011, 11, 59–71. [Google Scholar] [CrossRef]
  26. Bouman, B.A.M.; Kropff, M.J.; Tuong, T.P.; Wopereis, M.C.S.; Ten Berge, H.F.M.; van Laar, H.H. ORYZA2000: Modeling Lowland Rice; International Rice Research Institute: Los Baños, Philippines; Wageningen University and Research Centre: Wageningen, The Netherlands, 2001. [Google Scholar]
  27. Bouman, B.A.M.; van Laar, H.H. Description and evaluation of the rice growth model ORYZA2000 under nitrogen-limited conditions. Agric. Syst. 2006, 87, 249–273. [Google Scholar] [CrossRef]
  28. Wang, X.; Li, T.; Yang, X.; Zhang, T.; Liu, Z.; Guo, E.; Liu, Z.; Qu, H.; Chen, X.; Wang, L.; et al. Rice yield potential, gaps and constraints during the past three decades in a climate-changing Northeast China. Agric. For. Meteorol. 2018, 259, 173–183. [Google Scholar] [CrossRef]
  29. Tan, J.; Cui, Y.; Luo, Y. Assessment of uncertainty and sensitivity analyses for ORYZA model under different ranges of parameter variation. Eur. J. Agron. 2017, 91, 54–62. [Google Scholar] [CrossRef]
  30. Gao, Y.; Sun, C.; Ramos, T.B.; Huo, Z.; Huang, G.; Xu, X. Modeling nitrogen dynamics and biomass production in rice paddy fields of cold regions with the ORYZA-N model. Ecol. Model. 2023, 475, 110184. [Google Scholar] [CrossRef]
  31. Yu, Q.; Cui, Y.; Liu, L. Assessment of the parameter sensitivity for the ORYZA model at the regional scale—A case study in the Yangtze River Basin. Environ. Model. Softw. 2023, 159, 105575. [Google Scholar] [CrossRef]
  32. Liang, H.; Xu, J.; Hou, H.; Qi, Z.; Yang, S.; Li, Y.; Hu, K. Modeling CH4 and N2O emissions for continuous and noncontinuous flooding rice systems. Agric. Syst. 2022, 203, 103528. [Google Scholar] [CrossRef]
  33. Rocha, A.V.; Appel, R.; Bret-Harte, M.S.; Euskirchen, E.S.; Salmon, V.; Shaver, G. Solar position confounds the relationship between ecosystem function and vegetation indices derived from solar and photosynthetically active radiation fluxes. Agric. For. Meteorol. 2021, 298–299, 108291. [Google Scholar] [CrossRef]
  34. Yan, Z.; Zhou, Q.; Teng, M.; Ji, H.; Zhang, J.; He, W.; Ye, Y.; Wang, B.; Wang, P. High planting density and leaf area index of masson pine forest reduce crown transmittance of photosynthetically active radiation. Glob. Ecol. Conserv. 2019, 20, e00759. [Google Scholar] [CrossRef]
  35. Zhu, X.; He, H.; Liu, M.; Yu, G.; Sun, X.; Gao, Y. Spatio-temporal variation of photosynthetically active radiation in China in recent 50 years. J. Geogr. Sci. 2010, 20, 803–817. [Google Scholar] [CrossRef]
  36. Tan, J.; Cui, Y.; Wang, W. Global sensitivity analysis for ORYZA(v3) model under different rice growing environments in China. Trans. Chin. Soc. Agric. Eng. 2020, 36, 153–163. [Google Scholar] [CrossRef]
  37. Huang, M.; Jiang, L.; Zou, Y.; Zhang, W. On-farm assessment of effect of low temperature at seedling stage on early-season rice quality. Field Crops Res. 2013, 141, 63–68. [Google Scholar] [CrossRef]
  38. Xu, L.; Zhan, X.; Yu, T.; Nie, L.; Huang, J.; Cui, K.; Wang, F.; Li, Y.; Peng, S. Yield performance of direct-seeded, double-season rice using varieties with short growth durations in central China. Field Crops Res. 2018, 227, 49–55. [Google Scholar] [CrossRef]
  39. Tu, D.; Wu, W.; Xi, M.; Zhou, Y.; Xu, Y.; Chen, J.; Shao, C.; Zhang, Y.; Zhao, Q. Effect of Temperature and Radiation on Indica Rice Yield and Quality in Middle Rice Cropping System. Plants 2022, 11, 2697. [Google Scholar] [CrossRef] [PubMed]
  40. Liu, L.; Ouyang, W.; Liu, H.; Zhu, J.; Ma, Y.; Wu, Q.; Chen, J.; Zhang, D. Potential of paddy drainage optimization to water and food security in China. Resoiurce Conserv. Recycl. 2021, 171, 105624. [Google Scholar] [CrossRef]
  41. Wang, Y.; Guo, F.; Shen, H.; Xing, X.; Ma, X. Global Sensitivity Analysis and Evaluation of the DSSAT Model for Summer Maize (Zea mays L.) Under Irrigation and Fertilizer Stress. Int. J. Plant Prod. 2021, 15, 523–539. [Google Scholar] [CrossRef]
  42. Hao, S.; Ryu, D.; Western, A.W.; Perry, E.; Bogena, H.; Franssen, H.J.H. Global sensitivity analysis of APSIM-wheat yield predictions to model parameters and inputs. Ecol. Model. 2024, 487, 110551. [Google Scholar] [CrossRef]
  43. Cui, Y.; Li, Y.; Lu, G.; Sha, Z. Nitrogen movement and transformation with different water supply for paddy rice. Adv. Water Sci. 2004, 15, 280–285. [Google Scholar]
  44. Sun, H.; Zhu, Y.; Yang, J.; Wang, X. Global sensitivity analysis for an integrated model for simulation of nitrogen dynamics under the irrigation with treated wastewater. Environ. Sci. Pollut. Res. 2015, 22, 16664–16675. [Google Scholar] [CrossRef]
  45. Simlab. SIMLAB: Sensitivity Analysis Software—Joint Research Centre of the European Commission. Available online: https://ec.europa.eu/jrc/en/samo/simlab (accessed on 10 April 2022).
  46. Amiri, E.; Rezaei, M.; Rezaei, E.E.; Bannayan, M. Evaluation of Ceres-Rice, Aquacrop and Oryza2000 Models in Simulation of Rice Yield Response to Different Irrigation and Nitrogen Management Strategies. J. Plant Nutr. 2014, 37, 1749–1769. [Google Scholar] [CrossRef]
  47. Zhang, J.; Feng, L.; Zou, H.; Liu, D.L. Using ORYZA2000 to model cold rice yield response to climate change in the Heilongjiang province, China. Crop J. 2015, 3, 317–327. [Google Scholar] [CrossRef]
  48. Saltelli, A.; Ratto, M.; Andres, T.; Campolongo, F.; Cariboni, J.; Gatelli, D.; Saisana, M.; Tarantola, S. Global Sensitivity Analysis: The Primer; John Wiley & Sons: Chichester, UK, 2008. [Google Scholar]
  49. Cukier, R.I.; Levine, H.B.; Shuler, K.E. Nonlinear sensitivity analysis of multiparameter model systems. J. Comput. Phys. 1978, 26, 1–42. [Google Scholar] [CrossRef]
  50. Sobol, I.M. On the Distribution of Points in a Cube and the Appeoximate Evaluation of Integral. USSR Comput. Math. Math. Phys. 1967, 7, 86–112. [Google Scholar] [CrossRef]
  51. Chen, E.J.; Chen, Y.-Y.; Wei, L.-C.; Luo, H.-B.; Gong, W. Uncertainty and Sensitivity Analysis of Significant Parameters for Superlarge Diameter Shield Excavation. Adv. Civ. Eng. 2021, 2021, 8819393. [Google Scholar] [CrossRef]
  52. Makowski, D.; Naud, C.; Jeuffroy, M.-H.; Barbottin, A.; Monod, H. Global sensitivity analysis for calculating the contribution of genetic parameters to the variance of crop model prediction. Reliab. Eng. Syst. Saf. 2006, 91, 1142–1147. [Google Scholar] [CrossRef]
  53. Allen, R.G.; Pereira, L.S.; Raes, D.; Smith, M. Crop Evapotranspiration—Guidelines for Computing Crop Water Requirements—FAO Irrigation and Drainage Paper 56; FAO: Rome, Italy, 1998. [Google Scholar]
  54. Tarantola, A. Inverse Problem Theory and Methods for Model Parameter Estimation; SIAM: Society for Industrial and Applied Mathematics: Philadelphia, PA, USA, 2005. [Google Scholar]
Figure 1. First-order sensitivity indices (Si) and interaction effects of the ORYZA-N model at DVS = 0.4: crop-related outputs WAGT (ac), WST (df), and LAI (gi) for early rice, late rice, and single-season rice under traditional flood irrigation and shallow–wet irrigation conditions.
Figure 1. First-order sensitivity indices (Si) and interaction effects of the ORYZA-N model at DVS = 0.4: crop-related outputs WAGT (ac), WST (df), and LAI (gi) for early rice, late rice, and single-season rice under traditional flood irrigation and shallow–wet irrigation conditions.
Plants 13 00262 g001
Figure 2. First-order sensitivity indices (Si) and interaction effects of the ORYZA-N model at DVS = 0.65: crop-related outputs WAGT (ac), WST (df), LAI (gi), and WSO (jl) for early rice, late rice, and single-season rice under traditional flood irrigation and shallow–wet irrigation conditions.
Figure 2. First-order sensitivity indices (Si) and interaction effects of the ORYZA-N model at DVS = 0.65: crop-related outputs WAGT (ac), WST (df), LAI (gi), and WSO (jl) for early rice, late rice, and single-season rice under traditional flood irrigation and shallow–wet irrigation conditions.
Plants 13 00262 g002
Figure 3. First-order sensitivity indices (Si) and interaction effects of the ORYZA-N model at DVS = 1.0: crop-related outputs WAGT (ac), WST (df), LAI (gi), and WSO (jl) for early rice, late rice, and single-season rice under traditional flood irrigation and shallow–wet irrigation conditions.
Figure 3. First-order sensitivity indices (Si) and interaction effects of the ORYZA-N model at DVS = 1.0: crop-related outputs WAGT (ac), WST (df), LAI (gi), and WSO (jl) for early rice, late rice, and single-season rice under traditional flood irrigation and shallow–wet irrigation conditions.
Plants 13 00262 g003
Figure 4. First-order sensitivity indices (Si) and interaction effects of the ORYZA-N model at DVS = 2.0: crop-related outputs WAGT (ac), WST (df), LAI (gi), and WSO (jl) for early rice, late rice, and single-season rice under traditional flood irrigation and shallow–wet irrigation conditions.
Figure 4. First-order sensitivity indices (Si) and interaction effects of the ORYZA-N model at DVS = 2.0: crop-related outputs WAGT (ac), WST (df), LAI (gi), and WSO (jl) for early rice, late rice, and single-season rice under traditional flood irrigation and shallow–wet irrigation conditions.
Plants 13 00262 g004
Figure 5. First-order sensitivity indices (Si) and interaction effects of the ORYZA-N model at DVS = 0.4: nitrogen-related outputs ANCR (ac) and ANLV (df) for early rice, late rice, and single-season rice under traditional flood irrigation and shallow–wet irrigation conditions.
Figure 5. First-order sensitivity indices (Si) and interaction effects of the ORYZA-N model at DVS = 0.4: nitrogen-related outputs ANCR (ac) and ANLV (df) for early rice, late rice, and single-season rice under traditional flood irrigation and shallow–wet irrigation conditions.
Plants 13 00262 g005
Figure 6. First-order sensitivity indices (Si) and interaction effects of the ORYZA-N model at DVS = 0.65: nitrogen-related outputs ANCR (ac) and ANLV (df) for early rice, late rice, and single-season rice under traditional flood irrigation and shallow–wet irrigation conditions.
Figure 6. First-order sensitivity indices (Si) and interaction effects of the ORYZA-N model at DVS = 0.65: nitrogen-related outputs ANCR (ac) and ANLV (df) for early rice, late rice, and single-season rice under traditional flood irrigation and shallow–wet irrigation conditions.
Plants 13 00262 g006
Figure 7. First-order sensitivity indices (Si) and interaction effects of the ORYZA-N model at DVS = 1.0: nitrogen-related outputs ANCR (ac) and ANLV (df) for early rice, late rice, and single-season rice under traditional flood irrigation and shallow–wet irrigation conditions.
Figure 7. First-order sensitivity indices (Si) and interaction effects of the ORYZA-N model at DVS = 1.0: nitrogen-related outputs ANCR (ac) and ANLV (df) for early rice, late rice, and single-season rice under traditional flood irrigation and shallow–wet irrigation conditions.
Plants 13 00262 g007
Figure 8. First-order sensitivity indices (Si) and interaction effects of the ORYZA-N model at DVS = 2.0. Nitrogen-related outputs ANCR (ac) and ANLV (df) for early rice, late rice, and single-season rice under traditional flood irrigation and shallow–wet irrigation conditions.
Figure 8. First-order sensitivity indices (Si) and interaction effects of the ORYZA-N model at DVS = 2.0. Nitrogen-related outputs ANCR (ac) and ANLV (df) for early rice, late rice, and single-season rice under traditional flood irrigation and shallow–wet irrigation conditions.
Plants 13 00262 g008
Figure 9. First-order (Si) and total-order sensitivity indices (STi) results of the ORYZA-N model for crop yields of early (ER), late (LR), and single-season rice (SR) under traditional flood irrigation (TFI) and shallow–wet irrigation (SWI) conditions.
Figure 9. First-order (Si) and total-order sensitivity indices (STi) results of the ORYZA-N model for crop yields of early (ER), late (LR), and single-season rice (SR) under traditional flood irrigation (TFI) and shallow–wet irrigation (SWI) conditions.
Plants 13 00262 g009
Table 1. Soil physical properties for different soil layers of the experimental area.
Table 1. Soil physical properties for different soil layers of the experimental area.
Soil Horizon (cm)Soil Particle-Size Distribution (%)Bulk Density
(g cm−3)
Soil Organic Matter (%)
Sand
(2.0–0.05 mm)
Silt
(0.05–0.002 mm)
Clay (<0.002 mm)
Single-season rice in Heilongjiang province, China
0–1020.067.212.81.472.3
10–259.176. 914.01.483.3
25–4011.374.913.81.532.5
40–8017.871.710.51.294.0
Double-season rice in Jiangxi province, China
0–208.170.021.91.3432.1
20–5011.367.221.51.6171.4
50–1006.471.223.41.6750.9
Table 2. The observed phenology dates of the single-season rice, early rice, and late rice.
Table 2. The observed phenology dates of the single-season rice, early rice, and late rice.
PhenologySowingTransplantingPanicle InitiationFloweringMaturity
Dates for single-season rice15 April 201916 May 201910 July 20192 August 201921 September 2019
Dates for early rice22 March 201221 April 201223 May 201214 June 201210 July 2012
Dates for late rice21 June 201219 July 201222 August 201214 September 201219 October 2012
Note: single-season rice was observed in 2019 at Qing’an Irrigation Experiment Station, Heilongjiang province, China; early rice and late rice were observed in 2012 at Jiangxi Irrigation Experiment Station, Jiangxi province, China.
Table 3. Nitrogen fertilizer application schedule for the single-season rice, early rice, and late rice during the experimental years.
Table 3. Nitrogen fertilizer application schedule for the single-season rice, early rice, and late rice during the experimental years.
TreatmentSingle-Season Rice (2019)Early Rice (2012)Late Rice (2012)
Date
(Month–Day)
Net Nitrogen Amount
(kg N ha−1)
Date (Month–Day)Net Nitrogen Amount
(kg N ha−1)
Date (Month–Day)Net Nitrogen Amount
(kg N ha−1)
Base fertilizerApril 2757.6April 1990July 1790
Tillering fertilizerMay 2933.5April 2754July 2954
Panicle fertilizerJune 1546.4 May 2636August 2236
Total 137.5 180 180
Table 4. Definition, base values, and variance range of the selected parameters in the ORYZA-N model.
Table 4. Definition, base values, and variance range of the selected parameters in the ORYZA-N model.
ParametersDefinitionUnitDefault Values for Single-Season Rice (Variation Ranges)Default Values for Early Rice (Variation Ranges)Default Values for Late Rice
(Variation Ranges)
ULLSUpper limit of leaf rollingkPa74.13 (±30%)74.13 (±30%)74.13 (±30%)
LLLSLower limit of leaf rollingkPa794.33 (±30%)794.33 (±30%)794.33 (±30%)
ULDLUpper limit of drought-induced dead leaveskPa630.95 (±30%)630.95 (±30%)630.95 (±30%)
LLDLLower limit of drought-induced dead leaveskPa1584.89 (±30%)1584.89 (±30%)1584.89 (±30%)
ULLEUpper limit of leaf expansionkPa1.45 (±30%)1.45 (±30%)1.45 (±30%)
LLLELower limit of leaf expansionkPa1404 (±30%)1404 (±30%)1404 (±30%)
RGRLMXMaximum relative growth rate of leaf area(°Cd)−10.0088 (±30%)0.00877 (±30%)0.00904 (±30%)
FRPARFraction of photosynthetically active sunlight energy-0.45 (±30%)0.5 (±30%)0.62 (±30%)
DVRJDevelopment rate during the juvenile phase(°Cd)−10.000890 (±10%)0.000995 (±10%)0.001081 (±10%)
DVRIDevelopment rate during the photoperiod-sensitive phase(°Cd)−10.000758 (±10%)0.000758 (±10%)0.000758 (±10%)
DVRPDevelopment rate during the panicle phase(°Cd)−10.001055 (±10%)0.000981 (±10%)0.000853 (±10%)
DVRRDevelopment rate in reproductive phase (post-anthesis)(°Cd)−10.002000 (±10%)0.002139 (±10%)0.001728 (±10%)
FLV0.5Fraction of shoot dry matter partitioned to the leaves at DVS = 0.5-0.5 (±30%)0.6 (±30%)0.6 (±30%)
FLV0.75Fraction of shoot dry matter partitioned to the leaves at DVS = 0.75-0.3 (±30%)0.3 (±30%)0.3 (±30%)
FST0.75Fraction of shoot dry matter partitioned to the stems at DVS = 0.75-0.3 (±30%)0.7 (±30%)0.7 (±30%)
FST1.0Fraction of shoot dry matter partitioned to the stems at DVS = 1.0-0.6 (±30%)0.4 (±30%)0.4 (±30%)
DRLV1.0Fraction of leaf death coefficient at the DVS = 1.0-0.024 (±30%)0.015 (±30%)0.003 (±30%)
DRLV1.6Fraction of leaf death coefficient at the DVS = 1.6-0.020 (±30%)0.025 (±30%)0.005 (±30%)
Kdif1Solute diffusion coefficient of the first soil layer, which is related to solute type and conduction methodmm2
d−1
100 (±30%)100 (±30%)100 (±30%)
Kdif2Solute diffusion coefficient of the next 9 soil layermm2
d−1
600 (±30%)600 (±30%)600 (±30%)
KR0.0Relative uptake of solutes by roots from soil at DVS = 0.0-5 (±30%)5 (±30%)5 (±30%)
KR0.5Relative uptake of solutes by roots from soil at DVS = 0.5-3 (±30%)3 (±30%)3 (±30%)
KR2.0Relative uptake of solutes by roots from soil at DVS = 2.0-3 (±30%)3 (±30%)3 (±30%)
Note: single-season rice was observed in 2019 at Qing’an Irrigation Experiment Station, Heilongjiang province, China; early rice and late rice were observed in 2012 at Jiangxi Irrigation Experiment Station, Jiangxi province, China.
Table 5. Irrigation controlled criteria for two different irrigation regimes in various rice growing stages.
Table 5. Irrigation controlled criteria for two different irrigation regimes in various rice growing stages.
Irrigation
Regimes
SeedingEarly
Tillering
Middle
Tillering
Later
Tillering
Jointing
and Booting
BloomingGroutingYellow Ripe
Lower Limit Criterion—Upper Limit Criterion
Traditional flood irrigation (TFI)30 mm–
50 mm
30 mm–
50 mm
30 mm–
50 mm
Drainage30 mm–
50 mm
30 mm–
50 mm
30 mm–
50 mm
Naturally drying
Shallow–wet irrigation
(SWI)
10 mm–
30 mm
80%θs
15 mm
80%θs
15 mm
Drainage10 mm–
30 mm
10 mm–
30 mm
80%θs
15 mm
Naturally drying
Note: θs is the volumetric saturated soil water content in root layers; the lower limit is expressed as the ponded water depth or the percentage of θs, and the irrigation is applied when reaching the lower limit criterion; and the upper limit is the desired ponded water depth after irrigation.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Gao, Y.; Sun, C.; Ramos, T.B.; Tan, J.; Oliveira, A.R.; Huang, Q.; Huang, G.; Xu, X. Global Sensitivity Analysis of the Advanced ORYZA-N Model with Different Rice Types and Irrigation Regimes. Plants 2024, 13, 262. https://doi.org/10.3390/plants13020262

AMA Style

Gao Y, Sun C, Ramos TB, Tan J, Oliveira AR, Huang Q, Huang G, Xu X. Global Sensitivity Analysis of the Advanced ORYZA-N Model with Different Rice Types and Irrigation Regimes. Plants. 2024; 13(2):262. https://doi.org/10.3390/plants13020262

Chicago/Turabian Style

Gao, Ya, Chen Sun, Tiago B. Ramos, Junwei Tan, Ana R. Oliveira, Quanzhong Huang, Guanhua Huang, and Xu Xu. 2024. "Global Sensitivity Analysis of the Advanced ORYZA-N Model with Different Rice Types and Irrigation Regimes" Plants 13, no. 2: 262. https://doi.org/10.3390/plants13020262

APA Style

Gao, Y., Sun, C., Ramos, T. B., Tan, J., Oliveira, A. R., Huang, Q., Huang, G., & Xu, X. (2024). Global Sensitivity Analysis of the Advanced ORYZA-N Model with Different Rice Types and Irrigation Regimes. Plants, 13(2), 262. https://doi.org/10.3390/plants13020262

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop