Local and global sensitivity analysis for a prediction model of nitrogen loss in Southern China’s paddy fields via HYDRUS-1D

Abstract Nitrogen loss in paddy fields has been widely recognized as a significant contributor to nonpoint source pollution. Predicting this process through modeling is crucial, yet model parameters always carry uncertainty. Clarifying the time-dependent importance of the model parameters can help t...

Full description

Saved in:
Bibliographic Details
Main Authors: Shuhuai Wang, Juxiu Tong, Chen Huang
Format: Article
Language:English
Published: Nature Portfolio 2025-07-01
Series:Scientific Reports
Subjects:
Online Access:https://doi.org/10.1038/s41598-025-09858-3
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary:Abstract Nitrogen loss in paddy fields has been widely recognized as a significant contributor to nonpoint source pollution. Predicting this process through modeling is crucial, yet model parameters always carry uncertainty. Clarifying the time-dependent importance of the model parameters can help to better know the process effect such as precipitation or chemical reaction on nitrogen loss. Therefore, to rank the parameter importance, a global sensitivity analysis (GSA) named Sobol method is applied to a nitrogen loss model in paddy fields based on the soil mixing layer theory via modifying HYDRUS-1D model. To reduce the computational cost, local sensitivity analysis (LSA) is applied to the prediction model firstly, and three important model parameters, including soil mixing layer depth (d mix ), soil detachability coefficient (α) and precipitation intensity (p), are selected. Then, the Sobol method is applied to the prediction model to analyze the sensitivities of these three parameters. It is novel but reasonable that the Sobol sensitivity indices (including the first, second and total order indices, FOI, SOI and TOI) of d mix , α and p vary with time. The study results indicate that the importance of parameter varies with time during a rainfall. In surface runoff, α is most important at early times, while p becomes most important at later times for predicted urea and NO3 −-N concentrations. α is always the most sensitive parameter for the predicted NH4 +-N concentration in surface runoff. In soil, the GSA results are opposite for α and p. Generally, d mix is less important than α and p, and the interaction between each two parameters reflected by the SOIs has limited importance. d mix presented sensitivity in LSA but insensitivity in GSA. Sensitivity of α showed similar results in LSA (elasticity) and GSA (TOI and FOI), which decreases in surface runoff and increases in soil. All elasticities of p increase at first and decrease gradually later, while the GSA results of p vary oppositely in surface runoff and soil.
ISSN:2045-2322