Emergy analysis applied to the estimation of the recovery of costs for water services under the European Water Framework Directive

In this paper, the European Union's Water Framework Directive 2000/60/EC (WFD) that is intended to foster protection of water resources is examined, focusing on the improvement of ecological and chemical quality of surface and groundwater. The WFD includes the concept of full cost recovery (FCR) in accordance with the Polluter-Pays Principle, as one of the tools of an adequate and sustainable water resource management system. The WFD defines three different costs associated with water: resource costs (RC), financial costs (FC), and environmental costs (ECs).The FCR of water is examined from a biophysical perspective using emergy evaluation to: (1) establish resource values of water from different sources, (2) establish the full economic costs associated with supplying water, and (3) the societal costs of water that is used incorrectly; from which the resource costs, financial costs, and environmental costs, respectively, can be computed. Financial costs are the costs associated with providing water including energy, materials, labor and infrastructure. The emergy based monetary values vary between 0.15 and 1.73 €/m³ depending on technology. The emergy based, global average resource value (from which resource costs can be computed) is derived from two aspects of water: its chemical potential and its geopotential. The chemical potential monetary value of different sources such as rain, groundwater, and surface water derived from global averages of emergy inputs varies from 0.03 to 0.18 €/m³, depending on source, and the geopotential values vary from 0.03 to 2.40 €/m³, depending on location in the watershed. The environmental costs of water were averaged for the county of Spain and were 1.42 €/m³.Time of year and spatial location within the watershed ultimately influence the resource costs (computed from emergy value of chemical potential and geopotential energy) of water. To demonstrate this spatial and temporal variability, a case study is presented using the Foix watershed in northeastern Spain. Throughout the year, the resource value of water varies from 0.21 to 3.17 €/m³, depending on location within the watershed. It is concluded that FCR would benefit from the evaluation of resource costs using spatially and temporally explicit emergy accounting.     

Updated evaluation of exergy and emergy driving the geobiosphere: A review and refinement of the emergy baseline

Crucial to the method of emergy synthesis are the main driving emergy flows of the geobiosphere to which all other flows are referenced. They form the baseline for the construction of tables of Unit Emergy Values (UEVs) to be used in emergy evaluations. We provide here an updated calculation of the geobiosphere emergy baseline and UEVs for tidal and geothermal flows. First, we recalculate the flows using more recent values that have resulted from satellite measurements and generally better measurement techniques. Second, we have recalculated these global flows according to their available energy content (exergy) in order to be consistent with Odum's (1996) definition of emergy. Finally, we have reinterpreted the interaction of geothermal energy with biosphere processes thus changing the relationship between geothermal energy and the emergy baseline. In this analysis we also acknowledge the significant uncertainties related to most estimates of global data. In all, these modifications to the methodology have resulted in changes in the transformities for tidal momentum and geothermal energy and a minor change in the emergy baseline from 15.8E24 seJ/J to 15.2E24 seJ/J. As in all fields of science basic constants and standards are not really constant but change according to new knowledge. This is especially true of earth and ecological sciences where a large uncertainty is also to be found. As a consequence, while these are the most updated values today, they may change as better understanding is gained and uncertainties are reduced.     

The solar transformity of petroleum fuels

Petroleum fuels are the primary energy basis for transportation and industry. They are almost always an important input to the economic and social activities of humanity. Emergy analyses require accurate estimates with specified uncertainty for the transformities of major energy and material inputs to economic and environmental systems. In this study, the oil refining processes in Italy and the United States were examined to estimate the transformity and specific emergy of petroleum derivatives. Based on our assumptions that petroleum derivatives are splits of a complex hydrocarbon mixture and that the emergy is split based on the fraction of energy in a product, we estimated that the transformity of petroleum derivatives is 65,826 sej/J ± 1.4% relative to the 9.26E+24 sej/year planetary baseline. Estimates of the specific emergies of the various liquid fuels from Italian and U.S. refineries are within 2% of one another and the relationship of particular values varies with the refinery design. Our average transformity is only 1.7% larger than the current estimate for petroleum fuels determined by back calculation, confirming the accuracy of this transformity in existing emergy analyses. The model uncertainty between using energy or mass to determine how emergy is split was less that 2% in the estimate of both the transformity and specific emergy of liquid fuels, but larger for solid and gaseous products. This study is a contribution to strengthen the emergy methodology, providing data that can be useful in the analysis of many human activities.     

Carrying capacity of an uninhabited island off the southwestern coast of Korea

An emergy evaluation was carried out to assess the carrying capacity of a small, uninhabited island (Woosedo) off the southwestern coast of Korea. The sea area within 1 km from the high tide level of the island was included in the evaluation. The total environmental emergy input to Woosedo was 1.66E19 sej/yr, with the most emergy contribution from the tidal energy. The land and marine ecosystems of Woosedo contributed 4.97 million Em$ (7600 Em$/ha/yr) to the Korean economy annually. If Woosedo was developed to the national average at the emergy investment ratio of 2.86, its carrying capacity was estimated at 1034 people at the current living standard of Korea. With this population, the island system would not be sustainable with a very low emergy sustainability index of 0.36. At the same living standard used in the developed scenario, the carrying capacity of the island would be 370 people for a sustainable development scenario and 270 people if the renewable emergy were the only source to support the population. The emergy contribution of the marine ecosystem of the island was the major source of support in determining the level of carrying capacity of the island.     

Uncertainty characterization for emergy values

While statistical estimation of uncertainty has not typically accompanied published emergy values, as with any other quantitative model, uncertainty is embedded in these values, and lack of uncertainty characterization makes their accuracy not only opaque, it also prevents the use of emergy values in statistical tests of hypotheses. This paper first attempts to describe sources of uncertainty in unit emergy values (UEVs) and presents a framework for estimating this uncertainty with analytical and stochastic models, with model choices dependent upon on how the UEV is calculated and what kind of uncertainties are quantified. The analytical model can incorporate a broader spectrum of uncertainty types than the stochastic model, including model and scenario uncertainty, which may be significant in emergy models, but is only appropriate for the most basic of emergy calculations. Although less comprehensive in its incorporation of uncertainty, the proposed stochastic method is suitable for all types of UEVs. The distributions of unit emergy values approximate the lognormal distribution with variations depending on the types of uncertainty quantified as well as the way the UEVs are calculated. While both methods of estimating uncertainty in UEVs have their limitations in their presented stage of development, this paper provides methods for incorporating uncertainty into emergy, and demonstrates how this can be depicted and propagated so that it can be used in future emergy analyses and permit emergy to be more readily incorporated into other methods of environmental assessment, such as LCA.     

Assessment of a large watershed in Brazil using Emergy Evaluation and Geographical Information System

Humanity's future depends on the preservation of natural ecosystems that supply resources and absorb pollutants. Rural and urban productions are currently based on chemical products made from petroleum, which are responsible for high negative impacts on the Biosphere. In order to prevent those impacts, efficient public policies seeking for sustainable development are necessary. Aiming to assess the load on the environment (considering the gratuitous contributions of natural systems—a donor's perspective) due to human-dominated process, a scientific tool called Emergy Evaluation has been applied in different production systems, including crops and farms. However, there is still a lack of emergy studies in the context of watersheds, probably due to the difficulty of collecting raw data. The present work aims to carry out an assessment of Mogi-Guaçu and Pardo watershed, through the combined use of Emergy Evaluation and Geographical Information System. The agricultural and natural land uses were considered, while urban areas were excluded. Emergy flows (expressed in seJ ha⁻¹ yr⁻¹) obtained for all agricultural and natural land uses were expanded for the whole watershed and the emergy indices were calculated. The results show that the watershed has: low renewability (%R = 32%); low capture of natural resources through high external economic investment (EYR = 1.86); low dependence on natural resources (EIR = 1.16); and moderate load on the environment (ELR = 2.08). Considering a scenario where sugar-cane crops, orchards and pasture areas are converted from conventional to organic management, watershed's emergy performance improved, reaching a new renewability of 38%, but it is still not enough to be considered sustainable.     

Estimating Plecoglossus altivelis altivelis migration using a mass balance model expressed by hydrological distribution parameters in a major limpid river basin in Japan

We examined the influence of several hydrological and meteorological parameters on the migratory movements of ayu Plecoglossus altivelis altivelis in central Japan. When comprehensively evaluating rivers and ayu behaviour on a catchment scale, the subjects of analysis typically include human activities and hydrological and meteorological phenomena. However, limiting analyses to such factors may be too restrictive when human activities are being conducted. Accordingly, we incorporated a biological viewpoint into the evaluation method, analysing hydrological data (river discharge, river water temperature, sea water temperature) to determine watershed characteristics and examining the relationship between these characteristics and the habitat conditions of ayu. Then we constructed a numerical model for ayu migratory runs that incorporated ayu ecology and watershed characteristics. Analyses of ayu movements from a lower estuarine dam demonstrated that downstream displacements were associated with high water flows of more than 200 m³ s⁻¹ at the beginning of summer. We conclude that it is important to consider the effects of environmental parameters on the movements of different fish species to understand the causes of spatial variation in fish distribution in lowland rivers.     

Emergy analysis of the urban metabolism of Beijing

Cities can be modeled as if they were superorganisms with a range of metabolic processes. Research on this urban metabolism can contribute to solving urban environmental problems by revealing details of the metabolic throughput of the system. A key issue is how to find a common basis for measuring the environmental and economic values. By providing a single unified unit, emergy theory integrates the natural and socioeconomic systems and thoroughly evaluates a system's metabolism. We analyzed Beijing's urban metabolic system using emergy synthesis to evaluate its environmental resources, economy, and environmental and economic relations with the regions outside the city during 14 years of development. We compared Beijing's emergy indices with those of five other Chinese cities and of China as a whole to assess Beijing's relative development status. These indices are the emergy self-support ratio (metabolic dependence), the environmental load ratio (metabolic loading), empower density (metabolic pressure), emergy used per person (metabolic intensity per capita), and the monetary equivalent of emergy (emdollars; metabolic intensity). Based on our emergy analysis, Beijing's socioeconomic system is not self-sufficient, and depends greatly on external environmental resources. Its GDP is supported by a high percentage of emergy purchased from outside the city. During the study period, Beijing's urban system showed an increasing dependence on external resources for its economic development. Beijing's loading and pressure on the ecological environment is continuously increasing, accompanied by continuously increasing human emergy consumption. In the future, it will become increasingly necessary to improve Beijing's metabolic efficiency.     

Methods for estimating the uncertainty in emergy table-form models

Emergy studies have suffered criticism due to the lack of uncertainty analysis and this shortcoming may have directly hindered the wider application and acceptance of this methodology. Recently, to fill this gap, the sources of uncertainty in emergy analysis were described and analytical and stochastic methods were put forward to estimate the uncertainty in unit emergy values (UEVs). However, the most common method used to determine UEVs is the emergy table-form model, and only a stochastic method (i.e., the Monte Carlo method) was provided to estimate the uncertainty of values calculated in this way. To simplify the determination of uncertainties in emergy analysis using table-form calculations, we introduced two analytical methods provided by the Guide to the Expression of Uncertainty in Measurement (GUM), i.e., the Variance method and the Taylor method, to estimate the uncertainty of emergy table-form calculations for two different types of data, and compared them with the stochastic method in two case studies. The results showed that, when replicate data are available at the system level, i.e., the same data on inputs and output are measured repeatedly in several independent systems, the Variance method is the simplest and most reliable method for determining the uncertainty of the model output, since it considers the underlying covariance of the inputs and requires no assumptions about the probability distributions of the inputs. However, when replicate data are only available at the subsystem level, i.e., repeat samples are measured on subsystems without specific correspondence between an output and a certain suite of inputs, the Taylor method will be a better option for calculating uncertainty, since it requires less information and is easier to understand and perform than the Monte Carlo method.     

烟台市生态经济系统的能值分析

利用能值分析方法分析了2013年烟台市生态经济系统的能值流动状况,并提出了相应的能值利用调控对策。2013年烟台市生态经济系统的能值总用量为2.630×1023 sej,人均能值用量为4.044×10^16 sej,能值自给率为22.6%,能值密度为19.14×10^12sej/m^2,能值货币比率为2.848×10^12sej/,电力能值占能值总用量的比例为5.6%,由可更新资源及其产品支持的人口承载力为16.19×10^4人,由可更新资源及其产品、进口产品及技术共同支持的人口承载量为518.6×10^4人,能值可持续指标值为4.986。将上述计算结果与其它国家和地区进行比较表明：虽然目前烟台市经济较发达,能值利用效率和人民生活水平较高,生态经济系统的能值使用总体上符合区域可持续发展的要求,但系统能值自给率较低,生态环境比较脆弱,经济增长对外部输入的能源和资源依赖性较强。据此,提出了调整产业结构、转变经济增长方式、加强电力资源开发、发展科学教育事业等能值利用的调控对策。     

Assessing spawning ground hydraulic suitability for Chinese sturgeon (Acipenser sinensis) from horizontal mean vorticity in Yangtze River

The Chinese sturgeon, Acipenser sinensis, is an anadromous protected species. Its migration pattern has been blocked since the construction of Gezhouba dam and the length of the natural spawning site reduced to less than 7 km. However, the fish eventually established an alternative spawning ground in the narrow downstream reach of Gezhouba dam. In this article, we applied Delft3D-Flow model to simulate the hydraulic suitability of the spawning ground downstream Gezhouba dam. Horizontal mean vorticity was used to assess the hydraulic environment of spawning ground. The flow field state was determined through model simulation and field-measured data used to validate the model. The computational method was improved by calculating absolute horizontal mean vorticity from estimates the literature. The final vorticity was determined from the simulation output and its distribution pattern retrieved. The horizontal mean vorticity range was 0.71–4.61 10⁻³ s⁻¹ for the entire spawning grounds, with egg mass field upper limit of 1.0 × 10⁻³ s⁻¹. Vorticity strength selection of Chinese sturgeon spawning can enhance our understanding of egg fertilization rate, hence the protection of fertilized eggs. Furthermore, the results of the study would add to existing scientific database for spawning ground hydraulic environment protection, especially in the ecological regulation drive of the Three Gorges Reservoir.     

Recycling flows in emergy evaluation: A mathematical paradox?

This paper is a contribution to the emergy evaluation of systems involving recycling or reuse of waste. If waste exergy (its residual usefulness) is not negligible, wastes could serve as input to another process or be recycled. In cases of continuous waste recycle or reuse, what then is the role of emergy? Emergy is carried by matter and its value is shown to be the product of specific energy with mass flow rate and its transformity. This transformity (τ) given as the ratio of the total emergy input and the useful available energy in the product (exergy) is commonly calculated over a specific period of time (usually yearly) which makes transformity a time dependent factor. Assuming a process in which a part of the non-renewable input is an output (waste) from a previous system, for the waste to be reused, an emergy investment is needed. The transformity of the reused or recycled material should be calculated based on the pathway of the reused material at a certain time (T) which results in a specific transformity value (τ). In case of a second recycle of the same material that had undergone the previous recycle, the material pathway has a new time (T + T₁) which results in a transformity value (τ₁). Recycling flows as in the case of feedback is a dynamic process and as such the process introduces its own time period depending on its pathway which has to be considered in emergy evaluations. Through the inspiration of previous emergy studies, authors have tried to develop formulae which could be used in such cases of continuous recycling of material in this paper. The developed approach is then applied to a case study to give the reader a better understanding of the concept. As a result, a ‘factor’ is introduced which could be included on emergy evaluation tables to account for subsequent transformity changes in multiple recycling. This factor can be used to solve the difficulties in evaluating aggregated systems, serve as a correction factor to up-level such models keeping the correct evaluation and also solve problems of memory loss in emergy evaluation. The discussion deals with the questions; is it a pure mathematical paradox in the rules of emergy? Is it consistent with previous work? What were the previous solutions to avoid the cumulative problem in a reuse? What are the consequences?     

Economic and environmental performance of electricity production in Finland: A multicriteria assessment framework

Meeting environmental, economic, and societal targets in energy policy is complex and requires a multicriteria assessment framework capable of exploring trade-offs among alternative energy options. In this study, we integrated economic analysis and biophysical accounting methods to investigate the performance of electricity production in Finland at plant and national level. Economic and environmental costs of electricity generation technologies were assessed by evaluating economic features (direct monetary production cost), direct and indirect use of fossil fuels (GER cost), environmental impact (CO₂ emissions), and global environmental support (emergy cost). Three scenarios for Finland's energy future in 2025 and 2050 were also drawn and compared with the reference year 2008. Accounting for an emission permit of 25 €/t CO₂, the production costs calculated for CHP, gas, coal, and peat power plants resulted in 42, 67, 68, and 74 €/MWh, respectively. For wind and nuclear power a production cost of 63 and 35 €/MWh were calculated. The sensitivity analysis confirmed wind power's competitiveness when the price of emission permits overcomes 20 €/t CO₂. Hydro, wind, and nuclear power were characterized by a minor dependence on fossil fuels, showing a GER cost of 0.04, 0.13, and 0.26 J/J_e, and a value of direct and indirect CO₂ emissions of 0.01, 0.04, and 0.07 t CO₂/MWh. Instead, peat, coal, gas, and CHP plants showed a GER cost of 4.18, 4.00, 2.78, and 2.33 J/J_e. At national level, a major economic and environmental load was given by CHP and nuclear power while hydro power showed a minor load in spite of its large production. The scenario analysis raised technological and environmental concerns due to the massive increase of nuclear power and wood biomass exploitation. In conclusion, we addressed the need to further develop an energy policy for Finland's energy future based on a diversified energy mix oriented to the sustainable exploitation of local, renewable, and environmentally friendly energy sources.     

Assessing geobiosphere work of generating global reserves of coal, crude oil, and natural gas

A teacher of ours used to say, “Like ice in a fire, something for nothing you will never acquire”, which is a poetic equivalent of “there is no such a thing as a free lunch”. Human economies are dependent on high quality fossil fuels and will likely continue depending on them for some time to come. Value of a resource is not only what one pays for it, or what can be extracted from it, but also value can be attributed to the “effort” required in its production. In this analysis we apply the emergy synthesis method to evaluate the work invested by the geobiosphere to generate the global storages of fossil energy resources. The upgrading of raw resources to secondary fuels is also evaluated. The analysis relies on published estimates of historic, global net primary production (NPP) on land and oceans, published preservation and conversion factors of organic matter, and assessments of the present total global storages of coal, petroleum, and natural gas. Results show that the production of coal resources over geologic time required between 6.63E4 (±0.51E4) seJ/J and 9.71E4 (±0.79E4) seJ/J, while, oil and natural gas resources required about 1.48E5 (±0.07 E5) seJ/J and 1.70E5 (±0.06E5) seJ/J, respectively. These values are between 1.5 and 2.5 times larger than previous estimates and acknowledge a far greater power of fossil fuels in driving and shaping modern society.     

Formulation and testing of a novel river nitrification model

The nitrification process in many river water quality models has been approximated by a simple first order dependency on the water column ammonia concentration, while the benthic contribution has routinely been neglected. In this study a mathematical framework was developed for sediment bed nitrification based on mass transfer theory and Monod bacterial growth kinetics. The model describes ammonia transport across the boundary layer and consumption within the biofilm to quantify the overall nitrification flux. Model results suggest that nitrification is usually controlled by the boundary layer thickness, and we estimated a nitrification velocity range between 0.14 and 0.97 m d⁻¹, assuming typical boundary thicknesses of 0.1–1.0 mm. These ranges compared favorably with reported literature values, including our own measurements. The model was applied to several river systems of different depths where nitrification rates and river depths were available. Assuming that nitrification is exclusively a benthic process, the average velocity of all the rivers evaluated was 0.85 m d⁻¹ (r² = 0.72).     

The solar transformity of oil and petroleum natural gas

This paper presents an emergy evaluation of the biogeochemical process of petroleum formation. Unlike the previous calculation, in which the transformity of crude oil was back calculated from the relative efficiency of electricity production and factors relating coal to transportation fuels and transportation fuels to crude oil, we analyzed the geochemical process of petroleum formation (naftogenesis) to determine the transformities of oil and natural gas. We assumed that the process of oil and gas production is a steady state process in which all the emergy required is captured in the initial input. For such a system, we can use the mass concentration of the initial input to determine the specific emergy and transformity of the products. We used the maximum photosynthetic yield in Joules of phytoplankton organic matter per Joule of sunlight as the starting point. From this initial assumption, we traced the energy transformations in the oil and gas formation process through photosynthesis, death and decay of the phytoplankton, and diagenesis to kerogen production and from kerogen through catagenesis to petroleum formation. Our results show that both methods converge to similar values for oil (∼54,200 solar emJoules per Joule (sej/J)) and petroleum natural gas (43,500 sej/J) increasing our confidence in the results of past emergy analyses and providing a firm basis for the calculation of transformities for oil and gas derivatives.     

Modeling the carbon dynamics of the La Grande hydroelectric complex in northern Quebec

A mechanistic semi-empirical carbon cycle model of the La Grande reservoir complex in northern Quebec, Canada was conceived in order to investigate the climate impact of such a large alteration of the continental water cycle. The model includes inputs from the drainage basin, organic matter release from flooded soils, CO₂ emissions across the water-atmosphere interface and sedimentation. Most input data stems from previous research by our group on those ecosystems. The model includes the seven reservoirs of the La Grande complex and was run for periods of 50 and 100 years. Terrigeneous dissolved, particulate and suspended soil carbon fluxes and concentrations were computed. Over 100 years, 31.3 × 10¹² g C are released from flooded soils, equivalent to 28-29% of inputs from the drainage basin. 40-74% of dissolved organic carbon is mineralized. CO₂ fluxes over 100 years are 50.5-79.8 × 10¹² g C, 46.4-67.9 × 10¹² g C more than in the absence of reservoirs. The increase in mineralization of organic matter and in CO₂ emissions is a result of the increase in cumulated water residence time due to the creation of the reservoirs. Changes in other carbon sinks and sources likely offset a part of this additional carbon flux to the atmosphere. In the first years following flooding of the reservoir, organic carbon release from flooded soils exceeds CO₂ emissions, implying the downstream export of large quantities of eroded soil organic carbon. After this initial period, CO₂ emissions are fuelled by organic carbon originating from the drainage basin.     

Modelling llama population development under environmental and market constraints in the Bolivian highlands

Changes in the size of animal populations over time are mainly determined by demographic and environmental factors. Livestock population dynamics are additionally influenced by harvesting decisions taken by herders. In Bolivia, not much is known about current llama husbandry and the main influencing factors determining population sizes. We collected data on demography, environmental factors and market values affecting the current and future llama population in three different regions in Bolivia. We generated a population model and assessed the future development of the llama population including environmental factors (rangeland carrying capacity, disturbance phenomena), herd structure and dynamics, and economic market demands. We calibrated and validated the llama model on the basis of 20-year data sets of the regions of Oruro, Potosi and La Paz, respectively. Model calibration by means of the Gauss-Marquardt-Levenberg algorithm yielded a model efficiency of 0.94. For model validation, however, the simulation slightly overestimated the observed llama population yielding model efficiencies of 0.91 and 0.87 for Potosi and La Paz, respectively. Model outcomes were most sensitive to death and birth rates of juveniles and death rate of females compared to environmental or other demographic factors. Population trajectories approached an overall carrying capacity for Oruro, Potosi and La Paz of 8.8 × 10⁵, 9.1 × 10⁵, and 9.0 × 10⁵ llama individuals after 100 years of simulation. Hence, detailed monitoring of demographic, environmental, and economic factors can improve predictions of llama population development over time. Further management should focus on improving birth rates and lowering female mortality through providing supplemental food and shelters against the harsh environmental conditions of the Andean highlands.     

Emergy and emergy algebra explained by means of ingenuous set theory

Emergy is an important concept that has originated several effects in ecology, systems ecology and sustainability science. Its communication, however, has always presented several problems, since it does not follow the same rules of conservation as other energy-based approaches. Attempts have been made to clarify emergy by means of more formal/mathematical approaches, but the problem persists. In this paper, we have introduced a view of emergy and of its algebra based on ingenuous set theory. By means of this simple tool, emergy can be defined as the set of solar exergy that is directly and indirectly necessary to make a product. The operation that correctly sums the emergy “carried” by the inputs to a process is the union. This definition and the operation of union are able to account for all the rules of emergy algebra.     

Control division of agricultural non-point source pollution at medium-sized watershed scale in Southeast China

This paper presents the study carried out for controlling agricultural non-point source pollution (NSP) in a medium-sized watershed covering 1.47 × 10⁴ km² in Southeast China using quantitative analysis coupled with geographic information system (GIS), universal soil loss equation (USLE), soil conservation service-curve number (SCS-CN), nutrient loss equations, and annualized agricultural nonpoint source model (AnnAGNPS). Based on the quantitative results derived from GIS and environmental models, five control division units were generated for NSP control in Jiulong River watershed, namely, controlling unit for soil losses, controlling unit for livestock breeding and soil losses, controlling unit for excessive fertilizer use and livestock breeding, controlling unit for soil losses and fertilizer use, and controlling unit for excessive fertilizer use and soil losses. This study proved that integrating GIS with environmental models can be adopted to efficiently evaluate major sources and contributors of NSP, and identify the critical source areas of NSP, which enables adjusting measures to local conditions by further control division units developed through such study for control and management of water quality degradation induced by NSP in the Jiulong River watershed.