Calculation of the unit emergy value of water in an Italian watershed

Emergy is a thermodynamics-based entity that enables the implementation of a holistic environmental accounting system. It contributes to identify and measure all the inputs (energy and matter) supporting a given system, expressed in a common unit, namely solar emergy joule (sej). The emergy per unit product (called unit emergy value, UEV), is a measure of the environmental cost of a given resource. It is specific of the system/process and gives information on the dynamics, components and functioning of it. This paper presents the emergy evaluation of water resources within the watershed of the river Sieve, located in the Province of Florence (Italy). Along the river, an artificial basin has been created by means of a dam to preserve water quantity and quality, and to protect the Florentine area from dangerous floods and inundations. Different UEVs of water can be identified along the course of the river, especially upstream and downstream of the dam. These values quantify both the environmental and human efforts made to ensure and regulate the presence of water at different points of the river. The UEV of water flowing in the river increases from 1.35 × 10⁵ sej/g upstream, to 5.80 × 10⁵ sej/g downstream of the dam, depending mainly on man-made infrastructure. Along the watershed, three different systems of extraction, purification and distribution of water have been chosen on the basis of their dimension, type and location. UEVs of water distributed and the emergy investment necessary to implement different water management strategies are presented. The value of water purified and distributed decreases from 2.00 × 10⁶ sej/g for the smallest plant in the mountainous area, to 1.72 × 10⁶ sej/g for the largest plant, in the city of Florence, depending on production efficiency.     

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.     

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.     

Estimating the potential impact of vegetation on the water cycle requires accurate soil water parameter estimation

It is well known that vegetation dynamics at the catchment scale depends on the prevailing weather and soil moisture conditions. Soil moisture, however, is not equally distributed in space due to differences in topography, weather patterns, soil properties and the type and amount of vegetation cover. To elucidate the complex interaction between vegetation and soil moisture, the dynamic vegetation model LPJ-GUESS (Smith et al., 2001), which provides estimations of vegetation dynamics, but does not consider lateral water fluxes was coupled with the hydrological TOPMODEL (cf. Beven, 2001) in order to be able to evaluate the importance of these lateral fluxes. The new model LG-TM was calibrated and validated in two climatically different mountain catchments. The estimations of runoff were good, when monthly and weekly time scales were considered, although the low flow periods at winter time were somewhat underestimated. The uncertainty in the climate induced change vegetation carbon storage caused by the uncertainty in soil parameters was up to 3-5 kg C m⁻² (depending on elevation and catchment), compared to the total change in vegetation carbon storage of 5-9 kg C m⁻². Therefore accurate estimates of the parameters influencing the water holding capacity of the soil, for example depth and porosity, are necessary when estimating future changes in vegetation carbon storage. Similarly, changes in plant transpiration due to climatic changes could be almost double as high (88 mm m⁻²) in the not calibrated model compared to the new model version (ca 50 mm m⁻² transpiration change). The uncertainties in these soil properties were found to be more important than the lateral water exchange between grid cells, even in steep topography at least for the temporal and spatial resolution used here.     

Local calibration of coliforms parameters of water quality problem at Igapó I Lake, Londrina, Paraná, Brazil

This article presents results concerning the local calibration of the transport parameters (longitudinal and transversal diffusions and decay coefficient) for a two-dimensional problem of water quality at Igapó I Lake, located in Londrina, Paraná, Brazil, using fecal coliforms as an indicator of water quality. The simulation of fecal coliforms concentrations all over the water body is conducted by means of a structured discretization of the geometry of Igapó I Lake, together with the finite difference and finite element methods. By using the velocity field, modeled by the Navier-Stokes and Poisson equations, the flow of fecal coliforms is described by means of a transport model, which considers advective and diffusive processes, as well as a process of fecal coliforms decay. In the checkpoint, the longitudinal and transversal diffusion coefficients and the coliforms decay coefficient that best fitted the value of the fecal coliforms concentration were D_x = D_y = 0.001 m²/h and k = 0.5 d⁻¹ = 0.02083 h⁻¹. A qualitative and quantitative analysis of the numerical simulations conducted in function of the diffusion coefficients and of the coliforms decay parameter provided a better understanding of the local water quality at Igapó I Lake.     

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.     

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.     

Spatially distributed modeling of the long-term carbon balance of a boreal landscape

Spatially and temporally distributed information on the sizes of biomass carbon (C) pools (BCPs) and soil C pools (SCPs) is vital for improving our understanding of biosphere-atmosphere C fluxes. Because the sizes of C pools result from the integrated effects of primary production, age-effects, changes in climate, atmospheric CO₂ concentration, N deposition, and disturbances, a modeling scheme that interactively considers these processes is important. We used the InTEC model, driven by various spatio-temporal datasets to simulate the long-term C-balance in a boreal landscape in eastern Canada. Our results suggested that in this boreal landscape, mature coniferous stands had stabilized their productivity and fluctuated as a weak C-sink or C-source depending on the interannual variations in hydrometeorological factors. Disturbed deciduous stands were larger C-sinks (NEP₂₀₀₄ = 150 gC m⁻² yr⁻¹) than undisturbed coniferous stands (e.g. NEP₂₀₀₄ = 8 gC m⁻² yr⁻¹). Wetlands had lower NPP but showed temporally consistent C accumulation patterns. The simulated spatio-temporal patterns of BCPs and SCPs were unique and reflected the integrated effects of climate, plant growth and atmospheric chemistry besides the inherent properties of the C pool themselves. The simulated BCPs and SCPs generally compared well with the biometric estimates (BCPs: r = 0.86, SCPs: r = 0.84). The largest BCP biases were found in recently disturbed stands and the largest SCP biases were seen in locations where moss necro-masses were abundant. Reconstructing C pools and C fluxes in the ecosystem in such a spatio-temporal manner could help reduce the uncertainties in our understanding of terrestrial C-cycle.     

Carbon, nitrogen, oxygen and sulfide budgets in the Black Sea: A biogeochemical model of the whole water column coupling the oxic and anoxic parts

Carbon, nitrogen, oxygen and sulfide budgets are derived for the Black Sea water column from a coupled physical-biogeochemical model. The model is applied in the deep part of the sea and simulates processes over the whole water column including the anoxic layer that extends from ?115 m to the bottom (?2000 m). The biogeochemical model involves a refined representation of the Black Sea foodweb from bacteria to gelatinous carnivores. It includes notably a series of biogeochemical processes typical for oxygen deficient conditions with, for instance, bacterial respiration using different types of oxidants (i.e denitrification, sulfate reduction), the lower efficiency of detritus degradation, the ANAMMOX (ANaerobic AMMonium OXidation) process and the occurrence of particular redox reactions. The model has been calibrated and validated against all available data gathered in the Black Sea TU Ocean Base and this exercise is described in Gregoire et al. (2008). In the present paper, we focus on the biogeochemical flows produced by the model and we compare model estimations with the measurements performed during the R.V. KNORR expedition conducted in the Black Sea from April to July 1988 (Murray and the Black Sea Knorr Expedition, 1991). Model estimations of hydrogen sulfide oxidation, metal sulfide precipitation, hydrogen sulfide formation in the sediments and water column, export flux to the anoxic layer and to the sediments, denitrification, primary and bacterial production are in the range of field observations.With a simulated Gross Primary Production (GPP) of 7.9 mol C m⁻² year⁻¹ and a Community Respiration (CR) of 6.3 mol C m⁻² year⁻¹, the system is net autotrophic with a Net Community Production (NCP) of 1.6 mol C m⁻² year⁻¹. This NCP corresponds to 20% of the GPP and is exported to the anoxic layer. In order to model Particulate Organic Matter (POM) fluxes to the bottom and hydrogen sulfide profiles in agreement with in situ observations, we have to consider that the degradation of POM in anoxic conditions is less efficient that in oxygenated waters as it has often been observed (see discussion in Hedges et al., 1999). The vertical POM profile produced by the model can be fitted to the classic power function describing the oceanic carbon rate (C_R=Z^−α) using an attenuation coefficient α of 0.36 which is the value proposed for another anoxic environment (i.e. the Mexico Margin) by Devol and Hartnett (2001). Due to the lower efficiency of detritus degradation in anoxic conditions and to the aggregation of particles that enhanced the sinking, an important part of the export to the anoxic layer (i.e. 33%, 0.52 mol C m⁻² year⁻¹) escapes remineralization in the water column and reaches the sediments. Therefore, sediments are active sites of sulfide production contributing to 26% of the total sulfide production.In the upper layer, the oxygen dynamics is mainly governed by photosynthesis and respiration processes as well as by air-sea exchanges. ?71% of the oxygen produced by phytoplankton (photosynthesis+nitrate reduction) is lost through respiration, ?21% by outgasing to the atmosphere, ?5% through nitrification and only ?2% in the oxidation of reduced components (e.g. Mn²⁺, Fe²⁺, H₂S).The model estimates the amount of nitrogen lost through denitrification at 307 mmol N m⁻² year⁻¹ that can be partitioned into a loss of ?55% through the use of nitrate for the oxidation of detritus in low oxygen conditions, ?40% in the ANAMMOX process and the remaining ?5% in the oxidation of reduced substances by nitrate.In agreement with data analysis performed on long time series collected since the 1960s (Konovalov and Murray, 2001), the sulfide and nitrogen budgets established for the anoxic layer are not balanced in response to the enhanced particle fluxes induced by eutrophication: the NH₄ and H₂S concentrations increase.     

Interpreting spatial heterogeneity of crop yield with a process model and remote sensing

A process-based crop growth model (Vegetation Interface Processes (VIP) model) is used to estimate crop yield with remote sensing over the North China Plain. Spatial pattern of the key parameter—maximum catalytic capacity of Rubisco (V_cmax) for assimilation is retrieved from Normalized Difference of Vegetation Index (NDVI) from Terra-MODIS and statistical yield records. The regional simulation shows that the agreements between the simulated winter wheat yields and census data at county-level are quite well with R² being 0.41-0.50 during 2001-2005. Spatial variability of photosynthetic capacity and yield in irrigated regions depend greatly on nitrogen input. Due to the heavy soil salinity, the photosynthetic capacity and yield in coastal region is less than 50 μmol C m⁻² s⁻¹ and 3000 kg ha⁻¹, respectively, which are much lower than that in non-salinized region, 84.5 μmol C m⁻² s⁻¹ and 5700 kg ha⁻¹. The predicted yield for irrigated wheat ranges from 4000 to 7800 kg ha⁻¹, which is significantly larger than that of rainfed, 1500-3000 kg ha⁻¹. According to the path coefficient analysis, nitrogen significantly affects yield, by which water exerts noticeably indirect influences on yield. The effect of water on yield is regulated, to a certain extent, by crop photosynthetic capacity and nitrogen application. It is believed that photosynthetic parameters retrieved from remote sensing are reliable for regional production prediction with a process-based model.     

National water footprint in an input–output framework—A case study of China 2002

Water footprint is an indicator showing the consumption impact on environment with water equivalent, which allows a detailed quantification of water use directly and indirectly in sectors for the domestic and abroad final consumption. This paper presents a framework of calculating national water footprint (NWF) with input–output method on China 2002. The results indicated that the NWF of China was 381 m³/cap yr in 2002. A new indicator termed as national water footprint intensity (NWFI) is also derived from NWF to evaluate the water consumption intensities of different sectors. Then the evaluation of virtual water trade in sectors is followed, the results of which give a detailed quantification of net virtual water import of sectors, verifying that China is a net virtual water exporter concerning the whole national sectors, which is different from the previous studies. Finally, it is suggested that the sectors with high NWFI and volume of net virtual water export should be regarded as the priority of Chinese water-saving strategy.     

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).     

Modeling estuarine-shelf exchanges in a deltaic estuary: Implications for coastal carbon budgets and hypoxia

The export of wetland-derived materials to the coastal ocean (i.e., the “Outwelling” hypothesis) has received considerable attention over the past several decades. While a number of studies have shown that estuaries export appreciable amounts of nutrients and carbon, few studies have attempted to estimate the importance of estuarine sources for the coastal carbon budgets in river-dominated coastal ecosystems. A novel tidal prism model was developed to examine estuarine-shelf exchanges in the Barataria estuary, a deltaic estuary located in the north-central Gulf of Mexico. This estuary has been the site of a massive wetland loss, and it has been hypothesized that carbon export from the eroding coastal wetlands supports the development of a large hypoxic zone in the coastal Gulf of Mexico. The model results show that the Barataria estuary receives nitrogen through the tidal passes and releases carbon to the coastal ocean. The mean calculated tidal water discharge of 6930 m³ s⁻¹ is equivalent to about 43% of the lower Mississippi River discharge. The annual total organic carbon (TOC) export is 109 million kg, or 57 gC m² yr⁻¹ when prorated to the total water area of the estuary. This carbon export is equivalent to a loss of 0.5 m of wetland soil horizon over an area of 8.4 km², and accounts for about 34% of the observed annual wetland loss in the estuary between 1978 and 2000. Compared to the lower Mississippi River, the Barataria estuary appears to be a very small source of TOC for the northern Gulf of Mexico (2.7% of riverine TOC), and is unlikely to have a significant influence on the development of the Gulf's hypoxia.     

Rookery Bay and Naples Bay circulation simulations: Applications to tides and fresh water inflow regulation

A three-dimensional, baroclinic numerical circulation model is applied to the combined Rookery Bay and Naples Bay estuary complex. We focus on two aspects of the model simulations: the tidal circulation, for which we have observations of sea level and currents at several stations spanning Rookery Bay, and the salinity variations, for which we have observations at one station. After establishing model veracity through quantitative comparisons with observations, we use the model to address salinity changes that could arise due to regulatory fresh water flow rate modifications. Two different sets of experiments are performed. The first considers the salinity changes in Rookery Bay by increasing the fresh water flow rates into Henderson Creek through the addition of either 1.4 m³/s or 2.8 m³/s to the ambient flow rates. The second considers the salinity changes in Naples Bay by decreasing the fresh water flow rates through the subtraction of 9.9 m³/s from the Golden Gate Canals’ ambient flow rates. All of these scenarios result in significant, quantifiable salinity changes within the Rookery Bay and Naples Bay estuary complexes, demonstrating the environmental management utility of using high resolution, three-dimensional, numerical circulation models, as applied herein, for assessing the water properties of complex, multi-connected estuarine systems.     

Estimating soil carbon turnover using radiocarbon data: A case-study for European Russia

Turnover rates of soil carbon for 20 soil types typical for a 3.7 million km² area of European Russia were estimated based on ¹⁴C data. The rates are corrected for bomb radiocarbon which strongly affects the topsoil ¹⁴C balance. The approach is applied for carbon stored in the organic and mineral layers of the upper 1 m of the soil profile. The turnover rates of carbon in the upper 20 cm are relatively high for forest soils (0.16–0.78% year⁻¹), intermediate for tundra soils (0.25% year⁻¹), and low for grassland soils (0.02–0.08% year⁻¹) with the exception of southern Chernozems (0.32% year⁻¹). In the soil layer of 20–100 cm depth, the turnover rates were much lower for all soil types (0.01–0.06% year⁻¹) except for peat bog soils of the southern taiga (0.14% year⁻¹). Combined with a map of soil type distribution and a dataset of several hundred soil carbon profiles, the method provides annual fluxes for the slowest components of soil carbon assuming that the latter is in equilibrium with climate and vegetation cover. The estimated carbon flux from the soil is highest for forest soils (12–147 gC/(m² year)), intermediate for tundra soils (33 gC/(m² year)), and lowest for grassland soils (1–26 gC/(m² year)). The approach does not distinguish active and recalcitrant carbon fractions and this explains the low turnover rates in the top layer. Since changes in soil types will follow changes in climate and land cover, we suggest that pedogenesis is an important factor influencing the future dynamics of soil carbon fluxes. Up to now, both the effect of soil type changes and the clear evidence from ¹⁴C measurements that most soil organic carbon has a millennial time scale, are basically neglected in the global carbon cycle models used for projections of atmospheric CO₂ in 21st century and beyond.     

The sensitivity of carbon sequestration to harvesting and climate conditions in a temperate cypress forest: Observations and modeling

The role of disturbance and climate factors in determining the forest carbon balance was investigated at a Japanese cypress forest in central Japan with eddy flux measurements, tree-ring analyses, and a terrestrial biosphere model. The forest was established as a plantation after intermittent harvesting and replanting between 1959 and 1977, and acted as a strong carbon sink of approximately 500 g C m⁻² year⁻¹ for the measurement years between 2001 and 2007. A terrestrial biosphere model, BIOME-BGC, was validated using the eddy flux data at daily to interannual timescales, and the tree-ring width data at interannual to decadal timescales. According to the model simulation, during the observation period 270 ± 55 g C m⁻² year⁻¹ was additionally sequestered due to the indirect effects of the harvesting and planting, whereas the increase of CO₂ concentration and the change in climate increased the sink of 110 ± 40 and 30 ± 80 g C m⁻² year⁻¹, respectively. The model simulation shows that the forest is now recovering from harvesting, and that harvesting is a more important determinant of the current carbon sink than either interannual climate anomalies or increased atmospheric CO₂ concentration. We found that harvesting with long rotation length could be effective management activity in order to increase carbon sequestration, if the harvested timber is converted into products with long lifecycles.     

Resistance and re-organization of an ecosystem in response to biological invasion: Some hypotheses

Using a dynamic model of Lake Chozas developed by Marchi et al. (2011), we tested three hypotheses about recovery of the indigenous community and water quality after radical changes caused by introduction of an invasive allochthonous crayfish, Procambarus clarkii:

1.

Can the lake resist the pressure of an invasive species, like P. clarkii, by adaptation?

2.

Can the ecosystem recover when all the crayfish are removed and low phosphorus concentrations persist in inflow water?

3.

Does the simulated recovery of submerged vegetation occur at a total phosphorus concentration below 100 mg TP m⁻³, as estimated by Scheffer et al. (1993), Scheffer (1997), Jeppesen et al. (1998) and Zhang et al. (2003)?

We obtained the following answers:

1.

Lake Chozas can at least partly resist by adaptation. A combination of possible parameter changes could lead to a significant increase in eco-exergy.

2.

Removal of the phosphorus represented by crayfish (by harvesting) implies complete recovery of the lake and its eco-exergy, albeit not necessarily with the same organisms having the same properties.

3.

The expected hysteresis created by introduction and harvesting of crayfish is observed under the following conditions: phytoplankton dominance at total phosphorus ≥ about 200-250 mg TP m⁻³ and submerged vegetation returns at total phosphorus < 100 mg TP m⁻³.

     

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.     

Estimation of annual potential evapotranspiration at regional scale based on the effect of moisture on soil respiration

Potential evapotranspiration (PET) is an important component of water cycle. For traditional models derived from the principle of aerodynamics and the surface energy balance, its calculation always includes many parameters, such as net radiation, water vapor pressure, air temperature and wind speed. We found that it can be acquired in an easier way in specific regions. In this study, a new PET model (PETP model) derived from two empirical models of soil respiration was evaluated using the Penman-Monteith equation as a standard method. The results indicate that the PETP model estimation concur with the Penman-Monteith equation in sites where annual precipitation ranges from 717.71 mm to 1727.37 mm (R² = 0.68, p = 0.0002), but show large discrepancies in all sites (R² = 0.07, p = 0.1280). Then we applied our PETP model at the global scale to the regions with precipitation higher than 700 mm using 2.5° CMAP data to obtain the annual PET for 2006. As expected, the spatial pattern is satisfactory overall, with the highest PET values distributed in the lower latitudes or coastal regions, and with an average of 1292.60 ± 540.15 mm year⁻¹. This PETP model provides a convenient approach to estimate PET at regional scales.     

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?