Identifying source information after river chemical spill occurrences is critical for emergency responses. However, the inverse uncertainty characteristics of this kind of pollution source inversion problem have not yet been clearly elucidated. To fill this gap, stochastic analysis approaches, including a regional sensitivity analysis method, identifiability plot and perturbation methods, were employed to conduct an empirical investigation on generic inverse uncertainty characteristics under a well-accepted uncertainty analysis framework. Case studies based on field tracer experiments and synthetic numerical tracer experiments revealed several new rules. For example, the release load can be most easily inverted, and the source location is responsible for the largest uncertainty among the source parameters. The diffusion and convection processes are more sensitive than the dilution and pollutant attenuation processes to the optimization of objective functions in terms of structural uncertainty. The differences among the different objective functions are smaller for instantaneous release than for continuous release cases. Small monitoring errors affect the inversion results only slightly, which can be ignored in practice. Interestingly, the estimated values of the release location and time negatively deviate from the real values, and the extent is positively correlated with the relative size of the mixing zone to the objective river reach. These new findings improve decision making in emergency responses to sudden water pollution and guide the monitoring network design.
In this paper we describe and test a sub-model that integrates the cycling of carbon (C), nitrogen (N) and phosphorus (P) in the Soil Water Assessment Tool (SWAT) watershed model. The core of the sub-model is a multi-layer, one-pool soil organic carbon (SC) algorithm, in which the decomposition rate of SC and input rate to SC (through decomposition and humification of residues) depend on the current size of SC. The organic N and P fluxes are coupled to that of C and depend on the available mineral N and P, and the C:N and N:P ratios of the decomposing pools. Tillage explicitly affects the soil organic matter turnover rate through tool-specific coefficients. Unlike most models, the turnover of soil organic matter does not follow first order kinetics. Each soil layer has a specific maximum capacity to accumulate C or C saturation (Sx) that depends on texture and controls the turnover rate. It is shown in an analytical solution that Sx is a parameter with major influence in the model C dynamics. Testing with a 65-yr data set from the dryland wheat growing region in Oregon shows that the model adequately simulates the SC dynamics in the topsoil (top 0.3 m) for three different treatments. Three key model parameters, the optimal decomposition and humification rates and a factor controlling the effect of soil moisture and temperature on the decomposition rate, showed low uncertainty as determined by generalized likelihood uncertainty estimation. Nonetheless, the parameter set that provided accurate simulations in the topsoil tended to overestimate SC in the subsoil, suggesting that a mechanism that expresses at depth might not be represented in the current sub-model structure. The explicit integration of C, N, and P fluxes allows for a more cohesive simulation of nutrient cycling in the SWAT model. The sub-model has to be tested in forestland and rangeland in addition to agricultural land, and in diverse soils with extreme properties such high or low pH, an organic horizon, or volcanic soils. 相似文献
