Selecting parameters for calibration via sensitivity analysis: An individual-based model of mosquitofish population dynamics

A stochastic individual-based model (IBM) of mosquitofish population dynamics in experimental ponds was constructed in order to increase, virtually, the number of replicates of control populations in an ecotoxicology trial, and thus to increase the statistical power of the experiments. In this context, great importance had to be paid to model calibration as this conditions the use of the model as a reference for statistical comparisons. Accordingly, model calibration required that both mean behaviour and variability behaviour of the model were in accordance with real data. Currently, identifying parameter values from observed data is still an open issue for IBMs, especially when the parameter space is large. Our model included 41 parameters: 30 driving the model expectancy and 11 driving the model variability. Under these conditions, the use of “Latin hypercube” sampling would most probably have “missed” some important combinations of parameter values. Therefore, complete factorial design was preferred. Unfortunately, due to the constraints of the computational capacity, cost-acceptable “complete designs” were limited to no more than nine parameters, the calibration question becoming a parameter selection question. In this study, successive “complete designs” were conducted with different sets of parameters and different parameter values, in order to progressively narrow the parameter space. For each “complete design”, the selection of a maximum of nine parameters and their respective n values was carefully guided by sensitivity analysis. Sensitivity analysis was decisive in selecting parameters that were both influential and likely to have strong interactions. According to this strategy, the model of mosquitofish population dynamics was calibrated on real data from two different years of experiments, and validated on real data from another independent year. This model includes two categories of agents; fish and their living environment. Fish agents have four main processes: growth, survival, puberty and reproduction. The outputs of the model are the length frequency distribution of the population and the 16 scalar variables describing the fish populations. In this study, the length frequency distribution was parameterized by 10 scalars in order to be able to perform calibration. The recently suggested notion of “probabilistic distribution of the distributions” was also applied to our case study, and was shown to be very promising for comparing length frequency distributions (as such).     

Simulating the fate of Florida Snowy Plovers with sea-level rise: Exploring research and management priorities with a global uncertainty and sensitivity analysis perspective

Changes in coastal habitats due to sea-level rise provide an uncertain, yet significant threat to shoreline dependent birds. Rising sea levels can cause habitat fragmentation and loss which can result in considerable reduction in their foraging and nesting areas. Computational models and their algorithmic assumptions play an integral role in exploring potential mitigation responses to uncertain and potentially adverse ecological outcomes. The presence of uncertainty in metapopulation models is widely acknowledged but seldom considered in their development and evaluation, specifically the effects of uncertain model inputs on the model outputs. This paper was aimed to (1) quantify the contribution of each uncertain input factor to the uncertainty in the output of a metapopulation model which evaluated the effects of long-term sea-level rise on the population of Snowy Plovers (Charadrius alexandrinus) found in the Gulf Coast of Florida, and (2) determine the ranges of model inputs that produced a specific output for the purpose of formulating environmental management decisions. This was carried out by employing global sensitivity and uncertainty analysis (GSA) using two generic (model independent) methods, the qualitative screening Morris method and a quantitative variance-based Sobol’ method coupled with Monte Carlo filtering. The analyses were applied to three density dependence scenarios: assuming a ceiling-type density dependence, assuming a contest-type density dependence, and assuming that density dependence is uncertain as to being ceiling- or contest-dependent. The sources of uncertainty in the outputs depended strongly on the type of density dependence considered in the model. In general, uncertainty in the outputs highly depended on the uncertainty in stage matrix elements (fecundity, adult survival, and juvenile survival), dispersal rate from central areas with low current populations (the “Big Bend” area of Florida) to the northern, panhandle populations, the maximum growth rate, and density dependence type. Our results showed that increasing the maximum growth rate to a value of 1.2 or larger will increase the final average population of Snowy Plovers assuming a contest-type density dependence. Results suggest that studies that further quantify which density dependence relationship best describes Snowy Plover population dynamics should be conducted since this is the main driver of uncertainty in model outcomes. Furthermore, investigating the presence of Snowy Plovers in the Big Bend region may be important for providing connection between the panhandle and peninsula populations.     

Estimating the critical phosphorus loading of shallow lakes with the ecosystem model PCLake: Sensitivity, calibration and uncertainty

There is a vast body of knowledge that eutrophication of lakes may cause algal blooms. Among lakes, shallow lakes are peculiar systems in that they typically can be in one of two contrasting (equilibrium) states that are self-stabilizing: a ‘clear’ state with submerged macrophytes or a ‘turbid’ state dominated by phytoplankton. Eutrophication may cause a switch from the clear to the turbid state, if the P loading exceeds a critical value. The ecological processes governing this switch are covered by the ecosystem model PCLake, a dynamic model of nutrient cycling and the biota in shallow lakes. Here we present an extensive analysis of the model, using a three-step procedure. (1) A sensitivity analysis revealed the key parameters for the model output. (2) These parameters were calibrated on the combined data on total phosphorus, chlorophyll-a, macrophytes cover and Secchi depth in over 40 lakes. This was done by a Bayesian procedure, giving a weight to each parameter setting based on its likelihood. (3) These weights were used for an uncertainty analysis, applied to the switchpoints (critical phosphorus loading levels) calculated by the model. The model was most sensitive to changes in water depth, P and N loading, retention time and lake size as external input factors, and to zooplankton growth rate, settling rates and maximum growth rates of phytoplankton and macrophytes as process parameters. The results for the ‘best run’ showed an acceptable agreement between model and data and classified nearly all lakes to which the model was applied correctly as either ‘clear’ (macrophyte-dominated) or ‘turbid’ (phytoplankton-dominated). The critical loading levels for a standard lake showed about a factor two uncertainty due to the variation in the posterior parameter distribution. This study calculates in one coherent analysis uncertainties in critical phosphorus loading, a parameter that is of great importance to water quality managers.     

A global Bayesian sensitivity analysis of the 1d SimSphere soil–vegetation–atmospheric transfer (SVAT) model using Gaussian model emulation

Sensitivity analysis consists of an integral and important validatory check of a computer simulation model before the code is used in performing any kind of analysis operation. The present paper demonstrates the use of a relatively new method and tool for conducting global sensitivity analysis (GSA) for environmental models, providing simultaneously the first GSA study of the widely used 1d soil–vegetation–atmospheric transfer (SVAT) model named SimSphere. A software platform called the Gaussian emulation machine for sensitivity analysis (GEM SA), which has been developed for performing a GSA via Bayesian theory, is applied to SimSphere model in order to identify the most responsive model inputs to the simulation of key model outputs, detect their interactions and derive absolute sensitivity measures concerning the model structure. This study is also very timely in that, use of this particular SVAT model is currently being considered to be used in a scheme being developed for the operational retrieval of the soil surface moisture content by National Polar-orbiting Operational Environmental Satellite System (NPOESS), in a series of satellite platforms that are due to be launched in the next 12 years starting from 2016.The employed GSA method was found capable of identifying the most responsive SimSphere inputs and also of capturing their key interactions for each of the simulated target quantities on which the GSA was conducted. The most sensitive model inputs were the topography parameters (slope, aspect) as well as the fractional vegetation cover and soil surface moisture availability. The implications of these findings for the future use of SimSphere are discussed.     

Environ indicator sensitivity to flux uncertainty in a phosphorus model of Lake Sidney Lanier, USA

Effective environmental impact assessment and management requires improved understanding of the organization and transformation of ecosystems in which independent agents are linked through an intricate network of energy, matter, and informational interactions. While advances have been made, we still lack a complete understanding of the processes that create, constrain, and sustain ecosystems. Network environ analysis (NEA) provides one approach for building novel ecosystem insights, but it is model dependent. As ecological modeling is an imprecise art, often complicated by inadequate empirical data, the utility of NEA may be limited by model uncertainty. Here, we investigate the sensitivity of NEA indicators of ecosystem growth and development to flow and storage uncertainty in a phosphorus model of Lake Sidney Lanier, USA. The indicators are total system throughflow (TST), total system storage (TSS), total boundary input (Boundary), Finn cycling index (FCI), ratio of indirect-to-direct flows (Indirect/Direct), indirect flow index (IFI), network aggradation (AGG), network homogenization (HMG), and network amplification (AMP). Our results make two primary contributions. First, they demonstrate that five of the indicators – FCI, Indirect/Direct, IFI, AGG and HMG – are relatively robust to the flow and storage uncertainty in the Lake Lanier model. This stability lets us draw robust conclusions about the Lake Lanier ecosystem organization (e.g., phosphorus flux in the lake is dominated by internal processes) in spite of uncertainties in the model. Second, we show that the majority of the indicators co-vary and that most of their common variation could be mapped onto two latent factors, which we interpret as (1) system integration and (2) boundary influences.     

Experimental and numerical analysis of flow instabilities in rectangular shallow basins

Free surface flows in several shallow rectangular basins have been analyzed experimentally, numerically and theoretically. Different geometries, characterized by different widths and lengths, are considered as well as different hydraulic conditions. First, the results of a series of experimental tests are briefly depicted. They reveal that, under clearly identified hydraulic and geometrical conditions, the flow pattern is found to become non-symmetric, in spite of the symmetrical inflow conditions, outflow conditions and geometry of the basin. This non-symmetric motion results from the growth of small disturbances actually present in the experimental initial and boundary conditions. Second, numerical simulations are conducted based on a depth-averaged approach and a finite volume scheme. The simulation results reproduce the global pattern of the flow observed experimentally and succeed in predicting the stability or instability of a symmetric flow pattern for all tested configurations. Finally, an analytical study provides mathematical insights into the conditions under which the symmetric flow pattern becomes unstable and clarifies the governing physical processes.     

Uncertainty analysis of transient population dynamics

Two types of demographic analyses, perturbation analysis and uncertainty analysis, can be conducted to gain insights about matrix population models and guide population management. Perturbation analysis studies how the perturbation of demographic parameters (survival, growth, and reproduction parameters) may affect the population projection, while uncertainty analysis evaluates how much uncertainty there is in population dynamic predictions and where the uncertainty comes from. Previously, both perturbation analysis and uncertainty analysis were conducted on the long-term population growth rate. However, the population may not reach its equilibrium state, especially when there is management by harvesting or hunting. Recently, there has been an increased interest in short-term transient dynamics, which can differ from asymptotic long-term dynamics. There are currently techniques to conduct perturbation analyses of short-term transient dynamics, but no techniques have been proposed for uncertainty analysis of such dynamics. In this study, we introduced an uncertainty analysis technique, the general Fourier Amplitude Sensitivity Test (FAST), to study uncertainties in transient population dynamics. The general FAST is able to identify the amount of uncertainty in transient dynamics and contributions by different demographic parameters. We applied the general FAST to a mountain goat (Oreamnos americanus) matrix population model to give a clear illustration of how uncertainty analysis can be conducted for transient dynamics arising from matrix population models.     

Operating mechanism and set pair analysis model of a sustainable water resources system

There is no alternative to the world’s water resources, and their increasing scarcity is making it difficult to meet the world population’s water needs. This paper presents a sustainable water resources system (SWRS) and analyzes the operating mechanism that makes it possible to evaluate the status of such a system. A SWRS can be described as a complex coupling system that integrates water resources, social, economic and ecological systems into a whole. The SWRS’s operating mechanism is composed of dynamic, resistance and coordination components, and it interacts with and controls the system’s evolution process. The study introduces a new approach, set pair analysis theory, to measure the state of a SWRS, and an evaluation index system is established using the subsystems and operating mechanism of a SWRS. The evaluation index system is separated into three levels (goal level, criteria level and index level) and divides the index standard into five grades. An evaluation model of the SWRS based on set pair analysis theory is constructed, and an example of SWRS evaluation in Shanghai is presented. The connection degrees of the index in the three levels are calculated, and the connection degree of the goal index is calculated to be 0.342, which classifies the city’s SWRS condition as grade 2. The sustainable use of water resources in the region is determined to be at a relatively adequate level that meets the requirements of sustainable development.     

Evaluating the environmental flows of China's Wolonghu wetland and land use changes using a hydrological model, a water balance model, and remote sensing

Environmental flows are critical to sustaining a variety of plant and animal communities in wetlands. However, evaluation of environmental flows is hampered by the problem of hydrological and ecological data shortage, especially in many developing countries such as China. Based on a hydrological model, a water balance model and remote sensing data, we assessed the environmental flows of China's Wolonghu wetland with limited data. The hydrological model provides input data for the water balance model of the wetland, and the remote sensing data can be used to assess land use changes. Integration of these two models with the remote sensing data revealed both the environmental flows of the Wolonghu wetland and the relationships between these environmental flows and land use changes. The results demonstrate that environmental flows have direct and indirect influences on the wetland ecosystem and should be linked to sustainable wetland management.     

Uncertainty analysis of climate change mitigation options in the forestry sector using a generic carbon budget model

Industrialized countries agreed on a reduction of greenhouse gas emissions under the Kyoto Protocol. Many countries elected forest management activities and the resulting net balance of carbon emissions and removals of non-CO₂ greenhouse gases by forest management in their climate change mitigation measures. In this paper a generic dynamic forestry model (FORMICA) is presented. It has an empirical basis. Several modules trace C pools relevant for the Kyoto Protocol and beyond: biomass, litter, deadwood and soil, and harvested wood products. The model also accounts for the substitution of fossil fuels by wood products and bioenergy.     

Ecosystem network analysis indicators are generally robust to parameter uncertainty in a phosphorus model of Lake Sidney Lanier, USA

Understanding how data uncertainty influences ecosystem analysis is critical as we move toward ecosystem-based management. Here, we investigate how 18 Ecological Network Analysis (ENA) indicators that characterize ecosystem growth, development, and condition are affected by uncertainty in an ecosystem model of Lake Sidney Lanier (USA). We applied ENA to 122 plausible parameterizations of the ecosystem developed by Borrett and Osidele (2007, Ecological Modelling 200, 371-387), and then used the coefficient of variation (CV) to compare system indicator variability. We considered Total System Throughput (TST) as a measure of the underlying model uncertainty and tested three hypotheses. First, we hypothesized that non-ratio indicators whose calculation includes the TST would be at least as variable as TST if not more variable. Second, we postulated that indicators calculated as ratios, with TST in the numerator and denominator would tend to be less variable than TST because its influence will cancel. Last, we expected the Average Mutual Information (AMI) to be less variable than TST because it is a bounded function. Our work shows that the 18 indicators grouped into four categories. The first group has significantly larger CVs than the CV for TST. In this group, model uncertainty is amplified rendering these three indicators less useful. The second group of four indicators shows no significant difference in variability with respect to TST. Finally, there are two groups whose CV values are significantly lower than that for TST. The least variable group includes the ratio-based indicators and Average Mutual Information. Due to their low variability, we conclude that these indicators are the most robust to the parameter uncertainty and most useful for ecosystem assessment and comparative ecosystem analysis. In summary, this work suggests that we can be as certain, or more certain, in most of the selected ENA indicators as we are in the parameters of the model analyzed.     

Development of a coupled reservoir operation and water diversion model: Balancing human and environmental flow requirements

Maintaining a natural flow regime helps preserve the health of riverine ecosystems. Conventional studies on reservoir operations have focused mainly on identifying optimal operational schemes for satisfying human water demands. To systematically reflect the ecological effects of both natural and human-induced hydrologic alterations, water diversions downstream of the reservoirs should be considered as well. This research focused on a coupled reservoir operation and water diversion (CROWD) model, created through the integration of a reservoir operation model and a water diversion model. The proposed model considers both human and environmental flow requirements, and represents a compromise that balances ecological protection (preservation of the natural flow regime of a river) and human needs (reduced water shortages). In the reservoir operation model, the reservoir space is divided into three zones and different operating rules are developed for directing reservoir operation when water levels are in different zones; in the water diversion model, different water users are assigned different supply priorities with the instream flows no more than the minimum environmental flows having the highest priority; and the two models are coupled by the water mass balance between the two hydraulic facilities. The non-dominated-sorting genetic algorithm II (NSGA-II) was used to determine the parameters of the developed CROWD model and the model was applied to support the joint operational management of the Tanghe Reservoir and the Liaoyang Diversion in the Tang river basin, China. The resulting reservoir operation and water diversion schemes indicate that the CROWD model is useful for optimizing the operation of reservoirs and water diversion schemes. Moreover, it helps to analyze tradeoffs between human and environmental water needs, resulting in solutions that reduce the risk of water shortages and minimize ecological integrity disturbances.     

Identifying key life-traits for the dynamics and gene flow in a weedy crop relative: Sensitivity analysis of the GeneSys simulation model for weed beet (Beta vulgaris ssp. vulgaris)

The benefits of genetically modified herbicide-tolerant (GMHT) sugar beet (Beta vulgaris) varieties stem from their presumed ability to improve weed control and reduce its cost, particularly targeting weed beet, a harmful annual weedy form of the genus Beta (i.e. B. vulgaris ssp. vulgaris) frequent in sugar beet fields. As weed beet is totally interfertile with sugar beet, it is thus likely to inherit the herbicide-tolerance transgene through pollen-mediated gene flow. Hence, the foreseeable advent of HT weed beet populations is a serious threat to the sustainability of GM sugar beet cropping systems. For studying and quantifying the long-term effects of cropping system components (crop succession and cultivation techniques) on weed beet population dynamics and gene flow, we developed a biophysical process-based model called GeneSys-Beet in a previous study. In the present paper, the model was employed to identify and rank the weed life-traits as function of their effect on weed beet densities and genotypes, using a global sensitivity analysis to model parameters. Monte Carlo simulations with simultaneous randomization of all life-trait parameters were carried out in three cropping systems contrasting for their risk for infestation by HT weed beets. Simulated weed plants and bolters (i.e. beet plants with flowering and seed-producing stems) were then analysed with regression models as a function of model parameters to rank processes and life-traits and quantify their effects. Key parameters were those determining the timing and success of growth, development, seed maturation and the physiological end of seed production. Timing parameters were usually more important than success parameters, showing for instance that optimal timing of weed management operations is more important than its exact efficacy. The ranking of life-traits though depended on the cropping system and, to a lesser extent, on the target variable (i.e. GM weeds vs. total weed population). For instance, post-emergence parameters were crucial in rotations with frequent sugar beet crops whereas pre-emergence parameters were most important when sugar beet was rare. In the rotations with frequent sugar beet and insufficient weed control, interactions between traits were small, indicating diverse populations with contrasted traits could prosper. Conversely, when sugar beet was rare and weed control optimal, traits had little impact individually, indicating that a small number of optimal combinations of traits would be successful. Based on the analysis of sugar beet parameters and genetic traits, advice for the future selection of sugar beet varieties was also given. In climatic conditions similar to those used here, the priority should be given to limiting the presence of hybrid seeds in seed lots rather than decreasing varietal sensitivity to vernalization.     

Scenario analysis of water pollution control in the typical peri-urban river using a coupled hydrodynamic-water quality model

The water quality pollution and ecological deterioration in peri-urban rivers are usually serious under rapid urbanization and economic growth. In the study, a typical peri-urban river, Nansha River, was selected as a case study to discuss the scheme of peri-urban river rehabilitation. Located in the north part of the Beijing central region, the Nansha River watershed has been designated as an ecologically friendly garden-style area with high-tech industry parks and upscale residential zones. However, the Nansha River is currently seriously contaminated by urban and rural pollutants from both nonpoint sources (NPS) and point sources (PS). In this study, the pollutant loads from point sources and nonpoint sources in the Nansha River watershed were first assessed. A coupled model, derived from the Environmental Fluid Dynamics Code and Water Quality Analysis Simulation Program, was developed to simulate the hydrodynamics and water quality in the Nansha River. According to the characteristics of the typical peri-urban river, three different PS and NPS control scenarios were designed and examined by modeling analyses. Based on the results of the scenario analysis, a river rehabilitation scheme was recommended for implementation.     

Using global sensitivity analysis of demographic models for ecological impact assessment

Population viability analysis (PVA) is widely used to assess population‐level impacts of environmental changes on species. When combined with sensitivity analysis, PVA yields insights into the effects of parameter and model structure uncertainty. This helps researchers prioritize efforts for further data collection so that model improvements are efficient and helps managers prioritize conservation and management actions. Usually, sensitivity is analyzed by varying one input parameter at a time and observing the influence that variation has over model outcomes. This approach does not account for interactions among parameters. Global sensitivity analysis (GSA) overcomes this limitation by varying several model inputs simultaneously. Then, regression techniques allow measuring the importance of input‐parameter uncertainties. In many conservation applications, the goal of demographic modeling is to assess how different scenarios of impact or management cause changes in a population. This is challenging because the uncertainty of input‐parameter values can be confounded with the effect of impacts and management actions. We developed a GSA method that separates model outcome uncertainty resulting from parameter uncertainty from that resulting from projected ecological impacts or simulated management actions, effectively separating the 2 main questions that sensitivity analysis asks. We applied this method to assess the effects of predicted sea‐level rise on Snowy Plover (Charadrius nivosus). A relatively small number of replicate models (approximately 100) resulted in consistent measures of variable importance when not trying to separate the effects of ecological impacts from parameter uncertainty. However, many more replicate models (approximately 500) were required to separate these effects. These differences are important to consider when using demographic models to estimate ecological impacts of management actions.     

Comparison of sensitivity analysis techniques: A case study with the rice model WARM

The considerable complexity often included in biophysical models leads to the need of specifying a large number of parameters and inputs, which are available with various levels of uncertainty. Also, models may behave counter-intuitively, particularly when there are nonlinearities in multiple input-output relationships. Quantitative knowledge of the sensitivity of models to changes in their parameters is hence a prerequisite for operational use of models. This can be achieved using sensitivity analysis (SA) via methods which differ for specific characteristics, including computational resources required to perform the analysis. Running SA on biophysical models across several contexts requires flexible and computationally efficient SA approaches, which must be able to account also for possible interactions among parameters. A number of SA experiments were performed on a crop model for the simulation of rice growth (Water Accounting Rice Model, WARM) in Northern Italy. SAs were carried out using the Morris method, three regression-based methods (Latin hypercube sampling, random and Quasi-Random, LpTau), and two methods based on variance decomposition: Extended Fourier Amplitude Sensitivity Test (E-FAST) and Sobol’, with the latter adopted as benchmark. Aboveground biomass at physiological maturity was selected as reference output to facilitate the comparison of alternative SA methods. Rankings of crop parameters (from the most to the least relevant) were generated according to sensitivity experiments using different SA methods and alternate parameterizations for each method, and calculating the top-down coefficient of concordance (TDCC) as measure of agreement between rankings. With few exceptions, significant TDCC values were obtained both for different parameterizations within each method and for the comparison of each method to the Sobol’ one. The substantial stability observed in the rankings seem to indicate that, for a crop model of average complexity such as WARM, resource intensive SA methods could not be needed to identify most relevant parameters. In fact, the simplest among the SA methods used (i.e., Morris method) produced results comparable to those obtained by methods more computationally expensive.     

A population model for the peach fruit moth, Carposina sasakii Matsumura (Lepidoptera: Carposinidae), in a Korean orchard system

A population model for the peach fruit moth, Carposina sasakii Matsumura, was constructed to understand the population dynamics of this pest species and to develop an effective management strategy for various orchard (apple, peach, apple + peach) systems. The model was structured by the five developmental stages of C. sasakii: egg, larva, pupa, larval-cocoon (overwintering larva), and adult. The model consisted of a series of component models: (1) a bimodal spring adult emergence model, (2) an adult oviposition model, (3) stage emergence models of eggs, larvae, and pupae, (4) a larval survival rate model in fruits, (5) a larval-cocoon formation model, and (6) an insecticide effect model. Simulations using the model described the typical patterns of C. sasakii adult abundance in various orchard systems well, and was specific to the composition of host plants: three adult abundance peaks (first peak, mid-season peak, and last peak) a year with decreased peaks after the first peak in monoculture orchards of late apple, two adult peaks a year with a much higher last peak in monoculture orchards of early peach, and three adult peaks a year with much higher later peaks in mixed orchards of late apple and early peach. The average deviation between model outputs and actual records for first and second adult peak dates was 2.8 and 3.9 d, respectively, in simulations without an insecticide effect. The deviation decreased when insecticide effects were incorporated into the model. We also performed a sensitivity analysis of our model, and suggest possible applications of the model.     

Ecological network analysis for water use systems—A case study of the Yellow River Basin

Ecological network analysis (ENA) is introduced in this paper as a promising approach to study water use systems. Information indices from ENA involving total system throughput (TST), ascendency and overhead are calculated here. Two related aspects including organization inherent in system structures and synthesized water use intensity related with sustainable development of water use systems are analyzed. The indices of ascendency and overhead are applied for analyzing and characterizing water use network organization. For comparison of sustainability of water use systems from integrated aspects of environment, society and economy and based on TST, a new indicator termed as total system throughput intensity (TSTI) is constructed incorporating parameters of land, precipitation, population, GDP and environmental flow, which can be used as a measure of sustainability in terms of synthesized water use intensity. The Yellow River Basin in China during 1998–2006 is chosen as the case study and divided into subsystems according to the six river sections as from source to Lanzhou (S1-L1), Lanzhou to Toudaoguai (L1-T), Toudaoguai to Longmen (T-L2), Longmen to Sanmenxia (L2-S2), Sanmenxia to Huayuankou (S2-H) and Huayuankou to the mouth of Bo Sea (H-B). The results show that (i) the organization levels of L1-T and H-B are better than those of S1-L1 and T-L2, with those of L2-S2 and S2-H the worst; (ii) the synthesized water use intensity has been improving, of which T-L2, L2-S2 and S2-H are at the highest levels while H-B the lowest. In addition, the comparison between TSTI and other metrics and the relationship between ascendency and TSTI are discussed, from which the importance of TSTI is reflected and the optimization criterions for sustainable development of six subsystems are derived. It can be concluded that the application of ENA in water use systems can provide new angles for water resource management to address the challenges of assessing and optimizing options to obtain more sustainable water use.