Next Article in Journal
Rain Drop Measurement Techniques: A Review
Next Article in Special Issue
Shedding Light on Increasing Trends of Phosphorus Concentration in Upper Austrian Rivers
Previous Article in Journal
Extraction and Preference Ordering of Multireservoir Water Supply Rules in Dry Years
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Anthropogenic Nitrogen and Phosphorus Emissions and Related Grey Water Footprints Caused by EU-27′s Crop Production and Consumption

1
Twente Water Center, University of Twente, P.O. Box 217, 7500 AE Enschede, The Netherlands
2
Institute for Ecological Economics, Vienna University of Economics and Business, Welthandelsplatz 1/D4, 1020 Vienna, Austria
3
Energy Systems Research Group (GISE), Polytechnic Faculty, National University of Asunción, Campus San Lorenzo C.C. 2111, Paraguay
*
Author to whom correspondence should be addressed.
Water 2016, 8(1), 30; https://doi.org/10.3390/w8010030
Submission received: 19 October 2015 / Revised: 5 January 2016 / Accepted: 14 January 2016 / Published: 20 January 2016

Abstract

:
Water is a prerequisite for life on our planet. Due to climate change and pollution, water availability for agricultural production, industry and households is increasingly put at risk. With agriculture being the largest water user as well as polluter worldwide, we estimate anthropogenic nitrogen and phosphorus emissions to fresh water related to global crop production at a spatial resolution level of 5 by 5 arc min and calculate the grey water footprints (GWF) related to EU-27′s crop production. A multiregional input-output model is used to trace the the GWF embodied in the final consumption of crop products by the EU-27. The total GWF related to crop production in the EU-27 in 2007 was 1 × 1012 m3/year. Spain contributed about 40% to this total. Production of cereals (wheat, rice and other cereals) take the largest share, accounting for 30% of the GWF, followed by fruits (17%), vegetables (14%), and oil crops (13%). The total agricultural GWF of the EU-27 related to crop consumption was 1830 billion m3/year, which is 3700 m3/year per capita on average. Overall, the EU-27 was able to externalize about 41% of the GWF to the rest of the world through imports of crop products.

1. Introduction

The growing global population, crop production, sewage emissions and fossil fuel combustion have significantly altered the global nitrogen (N) and phosphorus (P) cycles. These changes of global nutrient cycles have both positive and negative effects [1]. On the one hand, the past few decades have seen an increase in crop yield that can be attributed to the increased use of N and P fertilizer in agriculture [2]. On the other hand, there is an imbalance between the applied nutrients in agriculture in the form of artificial fertilizer or manure and the nutrients taken up by the crops. Such imbalances—with crops not taking up the complete N and P provided—result in a large fraction of the mobilized N and P entering the freshwater system, causing degradation of the water quality and eutrophication of coastal marine ecosystems [2]. As a consequence, human-induced eutrophication of lakes, estuaries and coastal seas have been reported since the early twentieth century [3,4,5,6].
International trade of agricultural products has increased significantly over the last few decades. Taking into account the above described effects of agriculture, this increase results in a globalization of agricultural pollution [7,8,9]. Importing countries are shifting the environmental damage of agricultural production to regions that produce and export the products [7,8,9,10,11,12]. Hence, consumers do not have to bear the environmental cost and, as they may be far away from the regions that produce and supply the products, they are less aware of the environmental damages caused [7,8,9,10,11]. In order to increase awareness as well as to correctly assign responsibilities for environmental impacts to the responsible final consumer, it is necessary to quantify and account for the N and P emissions produced along the supply chain of the traded agricultural products and to relate them to final consumption. This need to measure the environmental pressure generated by intensive agricultural activities has led to the development of indicators that include the Ecological Footprint [13], the grey water footprint (GWF) [10] and the nitrogen footprint [14]. The water footprint (WF) is an indicator of both consumptive water use (water evaporated or incorporated into the product) and water required to assimilate pollutants. We can distinguish the WF into three colours: the green, blue and grey WF. The green WF refers to the consumption of rainwater, the blue WF refers to consumption of surface and ground water. The GWF is an indicator of water pollution. It is defined as the volume of freshwater that is required to assimilate the load of pollutants based on natural background concentration and existing ambient water quality standards [15]. The WF of a nation can be distinguished into: “the WF of national production” and “the WF national consumption”. The WF of national production refers to the total freshwater volume consumed or polluted within the territory of the nation as a result of activities within the agricultural sector, the industrial sector, and the domestic water supply sector. The WF of national consumption, on the other hand, is the total volume of freshwater that is consumed or polluted to produce the goods consumed by the nation′s population. We can further divide the WF of national consumption into two components: the internal and external WF of national consumption. The internal WF refers to the use of domestic water to produce goods consumed by the nation′s population. The external WF, on the other hand, is the volume of water used in other countries to produce goods consumed by the population in the nation under consideration.
Hoekstra and Mekonnen [12] have shown that 430 billion m3 per year of grey water have been traded globally over the period 1996–2005. During the same period, the global total external GWF (i.e., the part of the WF of national consumption falling outside the nation considered) was 218 billion m3 per year. In a recent study, O′Bannon et al. [7] showed that the global external GWF increased 136% during the 1986–2010 period.
The aim of this study was to estimate the quantity of nutrients (N and P) emitted and building on these estimates to calculate the GWF of the EU-27 related to both production and consumption of crop products. Further, the sustainability of the EU-27′s GWF of consumption is assessed by identifying priority basins and products. The analysis was carried out using the EXIOBASE v2.2 database (TNO, The Hague, The Netherlands) (base year 2007). EXIOBASE is a global, detailed multi-regional environmentally extended Supply and Use/Input Output (MR EE SUT/IOT) database [16]. The use of EXIOBASE allows tracking the environmental pressure caused by production processed through the supply chains to the final consumer.
The paper is structured as follows: In the following section, the method and data applied in the study is described. This is followed by a section presenting the results of our analysis. The final section provides a discussion of our main findings and the limitation of the study.

2. Materials and Methods

Nitrogen and phosphorus emissions to fresh water. Nitrogen loss through leaching (OUTlea, kg/ha) was calculated using the regression model developed by De Willigen [17]:
O U T l e a = ( 0.0463 + 0.0037 × ( P C × L ) ) × ( N f + α × N o m N u )
where P is the annual precipitation (mm/year), C the clay content (%), L the layer thickness or rooting depth (m), Nf the mineral and manure N (kg N/ha), α the decomposition rate (=1.6% per year), Nom the amount of N in soil organic matter (kg N/ha), and Nu the N uptake by harvested crop and crop residues removed from the field (kg of N/ha per year). Nom is calculated by multiplying bulk density (BD, kg/m3), total N in soil (TOTN, in kg/kg), and layer thickness L. The data on C, BD, and TOTN were obtained from the derived soil properties on a 5 × 5 arc-min global grid (version 1.2) from ISRICWISE [18]. The precipitation data for the period 2000–2010 were obtained from the Climate Research Unit of the University of East Anglia (version CRU TS 2.0) [19]. The 30-arc-min precipitation data were converted into 5-arc-min resolution data by assigning the 30-arc-min values to each of the 36 5-arc-min grid cells contained in the 30 by 30 arc-min grid cell. The rooting depth (L, m) for individual crops were obtained from Allen et al. [20].
The first part of the regression equation (0.0463 + 0.0037 × (P/C × L)) determines the fraction of the total mobile N that will leach, and the second part (Nf + α × Nom − Nu) determines the available mobile N [21]. To keep the fractional losses less than one, following De Willigen [17], an upper limit of 250 mm/(% m) was set for the combined parameter (P/C × L). Since perennials crops can take up N throughout the year, for these crops, the amount of N in soil organic matter was multiplied by 0.5 in order to further limit overestimation of N leaching [22]. We believe this assumption is reasonable given the fact that other authors have set an upper limit on the total leaching arguing the model might overestimate the leaching volume. For example, Liu et al. [23] have set an upper limit on the total leached amount not to exceed 20% of the total mineral fertilizer application and 50% of the N from manure arguing that the model may give unreasonably high leaching values.
We estimated the human-induced (i.e., due to artificial fertilizer and manure application) leaching quantity (Anthroplea) by multiplying the total leached volume (OUTlea) by the fraction of N input in the form of fertilizer and manure (Nf) to the total N input including the amount of N in the soil organic content:
A n t h r o p l e a = O U T l e a × N f ( N f + α × N o m )
Modeling the emission of P from diffuse sources at a global scale is more complicated than estimating N. Therefore, we have adopted the approach of Bouwman et al. [24], where the amount of P emitted to the water system is assumed to be 12.5% of the P input from the fertilizer and manure application. This value was obtained from the increase of total P river export between 1970 and 2000 as modeled by the global NEWS model [25] and corrected for the global average P retention of 20% in river systems [26].
IFA et al. [27] were the primary source of artificial fertilizer application rate per crop per country. This dataset was further complemented with FAO [28] and Heffer [29]. These application rates were adjusted to fit FAO [30] country average artificial fertilizer consumption per year for the period 2007. The manure input was calculated at grid cell level by multiplying livestock density (taken from FAO [31]) by the animal-specific excretion rates (taken from Sheldrick et al. [32]). The N and P removal with harvested crops was estimated by multiplying the crop yield by the crop-specific N and P content. Similarly, the N and P removal with crop residues was calculated by multiplying the yield of crop residue by the N and P content of the crop residue and a residual removal factor. The crop yield at 5 arc min per crop was taken from Mekonnen and Hoekstra [33]. The crop-specific N and P contents of major crop and crop residues were taken from IPNI [34]. For other crops and crop groupings, values from FAO [22] and Roy et al. [35] were used. For nuts and spices, we have adopted the values of fruits and vegetables respectively from FAO [22]. The crop residues removal factors were obtained from different sources [22,36,37,38,39,40].
Calculating the grey water footprint. The GWFs of N and P are first calculated separately, and then the overall GWF of N and P is determined by the nutrient that is most critical, i.e., the nutrient with the larger GWF. Following Hoekstra et al. [15], the GWF (WFgrey, m3/year) is calculated by dividing the N and P nutrients load (L, kg/year) by the difference between the ambient water quality standard for N and P nutrients (the maximum acceptable concentration (cmax, mg/L) and their natural concentration in the receiving water body (cnat, in mg/L)):
W F g r e y = L ( C max C n a t )
The natural concentration in a receiving water body is the concentration of the water systems that are in pristine condition before human disturbances occur in the catchment. We have adopted the maximum acceptable and natural concentration values for total N and total P from Liu et al. [41]. The maximum acceptable concentration values are 3.1 mg N/L and 0.95 mg P/L for total N and total P, respectively, and the natural concentration value of 1.5 mg N/L and 0.52 mg P/L for total N and total P, respectively. Although different ecosystems respond differently to nutrient loads, requiring different maximum allowable concentration values, we have taken a single value for the whole world. Besides, natural concentrations vary from basin to basin. However, obtaining or estimating basin-specific values for both maximum allowable and natural concentration is an elaborate task and impossible for a global study like this one.
To calculate the GWF related to the EU-27 consumption of goods and services, we extend the monetary MRIO model EXIOBASE 2.2 by the physical extensions consisting of GWF values related to the production of 146 different crops in the 48 different countries and regions of the EXIOBASE model. By allocating the GWF to the sectors producing the specific crops, the GWF is “linked” to the trade flows of the crops and ends with them in the countries of final consumption.
The GWF in a river basin is environmentally unsustainable, creating an environmental hotspot when the water pollution level exceeds the waste assimilation capacity. Those basins where the water pollution level creates environmental hotspot are identified as “priority river basins.” The priority river basins and products related to EU-27′s WF of consumption were identified on the basis of the size of the GWF of EU-27′s consumption and the water pollution level (WPL) obtained from Liu et al. [41]. Water pollution levels (WPL) can be calculated as the ratio of the total GWF in a catchment to the actual runoff from that catchment. A water pollution level of one means that the waste assimilation capacity has been fully consumed, as the complete runoff is “needed” to assimilate the emissions. When the water pollution level exceeds one, ambient water quality standards are violated. If the EU-27′s consumption causes a GWF in a basin with WPL higher than one, the products contributing a large share to this GWF are identified as priority products as per the definition the basin is defined as priority river basin (see above).

3. Results

In the first part of the results (Section 3.1), we will present the nutrient emission values that are used as input to the EXIOBASE database. Section 3.2 presents the GWF related to the production of crops within the EU-27. The GWF from the production perspective is the total GWF within the territory of the 27 EU countries related to the production of crops. Part of the GWF within the EU-27 can be used to produce goods consumed within the region, and the rest used for producing goods for export to the rest of the world. In Section 3.3, we present the GWF as the result of EU-27 consumption of crop products. For this section, EXIOBASE is used to calculate the GWF embodied in the final consumption of crop products by the EU-27. Part of the GWF related to EU-27 consumption is located within the EU-27 (internal WF) and part of it is located outside the boarder of EU-27 (external WF).

3.1. Anthropogenic Nitrogen and Phosphorus Emissions to Fresh Water Related to EU-27′s Crop Production

Nitrogen and phosphorus emissions to fresh water at high spatial resolution are presented in Figure 1. The results show that the global total N and P emissions from agricultural production in 2007 were 31 and 2.9 million tonnes per year, respectively. However, as Figure 1 illustrates, there are marked spatial variations in the anthropogenic N and P emissions to fresh water across the globe. In Southeastern China and Northeastern India, N loads can reach above 100 kg N/ha. Other places with large N loads include Western Europe, Central US, and Central America. In Europe, mainly in Germany, Netherlands, Belgium, the UK, Denmark, and Poland, N loads of up to 75 kg N/ha were observed. On the other hand, anthropogenic P load is largest in Northeastern Spain in the Ebro River basin where it reaches as much 50 kg P/ha. In Northeastern China and Northern India, the anthropogenic P load of a little over 5 kg/ha was observed.
The anthropogenic N and P emissions in 2007 within the EU-27 were 0.9 and 0.4 million tonnes per year, respectively. The contribution of major countries to the total EU-27 N and P loads are presented in Table 1. Ireland contributed the largest share (23%) of anthropogenic N load, followed by Germany (19%), UK (14%), Netherlands and Poland, each with 10%. About 50% of the anthropogenic P load in the EU-27 was contributed by Spain, followed by France (11%), and Poland (6%).
In terms of per unit of harvested area, the largest anthropogenic emission occurred in Ireland with 690 kg N/ha and 26 kg P/ha, followed by Netherlands (157 kg N/ha and 17 kg P/ha), Belgium (41 kg N/ha and 9 kg P/ha), and UK (33 kg N/ha and 5 kg P/ha). An excessive nutrient input from manure is mainly responsible for the larger emission in these countries.
Figure 1. Map of global anthropogenic loads of (a) nitrogen and (b) phosphorus due to global production of crops shown in kg/ha on a 5 by 5 arc-min grid. Period: 2007.
Figure 1. Map of global anthropogenic loads of (a) nitrogen and (b) phosphorus due to global production of crops shown in kg/ha on a 5 by 5 arc-min grid. Period: 2007.
Water 08 00030 g001
Table 1. Anthropogenic nitrogen and phosphorus emissions to water for selected European countries and EU-27 total.
Table 1. Anthropogenic nitrogen and phosphorus emissions to water for selected European countries and EU-27 total.
ProductAnthropogenic Nitrogen EmissionsAnthropogenic Phosphorus Emissions
k Tonne/Yearkg/hak Tonne/Yearkg/ha
Belgium254169
Bulgaria11462
France705453
Germany17619222
Greece13462
Ireland214690826
Italy193213
Netherlands901571017
Poland918262
Portugal6554
Spain58521217
UK13933225
Other countries211382
Total EU-27933114255

3.2. The Grey Water Footprint of EU-27 Related to Crop Production

Figure 2 illustrates the GWFs due to N and P emissions related to the crop production in the EU-27 in 2007. The overall GWF is the larger of the two values; the total overall GWF related to agricultural production in the EU-27 was 1065 billion m3 per year (Table 2). About 78% of this GWF is due to the anthropogenic P emissions, while the remaining 22% is due to anthropogenic N load from crop fields mainly cereals, fruits, and vegetables. Cereals (wheat, rice and other cereals) take the largest share, accounting for 30% of the GWF. The other crops with major contribution to the total GWF are fruits (17%), vegetables (14%) and oil crops (13%). The large GWF by these crops can be explained by their large extent of the harvested area. Cereals (wheat, rice and other cereals) account for 65% of the harvested area, followed by oil crops with 18%, and vegetables, fruits and nuts together with 11%. The other reasons for difference in the contribution of different crops to the overall GWF of EU-27 are the difference in the artificial fertilizer and manure application rate, and the crops nutrient removal rate.
The contribution of different countries to the overall GWF of the EU-27 related to crop production is presented in Figure 2. About 40% of the overall GWF was contributed by Spain, mainly due to anthropogenic P loads from cropland. The other major countries contributing a large share to the EU-27′s GWF of production are Ireland (11%), Germany (9%), France (9%), and the UK (7%). In Ireland, Germany and the UK, the larger GWF is due to N loads, while in France anthropogenic load was the dominant one. The large GWF in Ireland is mainly due to large manure application rates, while in Germany, France, and the UK it is mainly due to the large agricultural harvested area combined with high nutrient (fertilizer and manure) application rate.
Table 2. Grey water footprint of production of the EU-27 (million m3 per year). Period: 2007.
Table 2. Grey water footprint of production of the EU-27 (million m3 per year). Period: 2007.
ProductNitrogenPhosphorusOverall
Other cereals105,093190,130190,130
Fruits22,631175,891175,891
Vegetables30,962144,285144,285
Oil crops139,67495,105139,674
Wheat106,408115,376115,376
Pulses11,77585,15185,151
fodder crops62,58163,12163,121
Roots and tubers51,32432,72051,324
Sugar crops41,00327,68141,003
Nuts958031,33031,330
Paddy rice1512,86912,869
Other crops184911,77511,775
Fibres4329852985
Total582,938988,4191,064,914
Figure 2. Contribution of different countries to the overall GWF of EU-27 related to production. Period: 2007.
Figure 2. Contribution of different countries to the overall GWF of EU-27 related to production. Period: 2007.
Water 08 00030 g002

3.3. The Grey Water Footprint Related to EU-27 Consumption

The total agricultural GWF of the EU-27 related to consumption in 2007 was 1834 billion m3/year, which is 3700 m3/year per capita on average (Table 3). Looking at some selected European countries, the difference in the per-capita WF is quite large, ranging from 8900 m3/cap/year for Spain to 1600 m3/cap/year for Estonia.
The extent to which countries externalize their WF varies considerably. For the EU-27 as a whole, the share of the external WF was 41%. The share of the external GWF for the selected countries varies from 26% for Spain to 84% for Estonia. The extent to which single EU-27 countries externalize their GWF to other EU-27 countries or to the rest of the world also varies. For Spain, 42% of its GWF was external but only 4% was externalized to the EU-27, while 22% was to the rest of the world. Estonia, on the other hand, shifted 26% of its grey WF to the EU-27 and 57% to the rest of the world. The large share of the external WF shows that the EU-27 and the other countries displaced a large portion of the water pollution to the rest of the world.
Table 3. Grey water footprint of the EU-27′s and some selected countries related to consumption of agricultural products. Period: 2007.
Table 3. Grey water footprint of the EU-27′s and some selected countries related to consumption of agricultural products. Period: 2007.
CountryTotal GWF (Gm3/Year)Per capita GWF (m3/Cap/Year)
InternalExternal GWFTotal
Within EU-27Rest of the World
EU-27108075418343694
Estonia0.30.61.22.11608
Germany76741783283924
Netherlands171746804839
Poland461120772023
Spain29415903998931
UK72641292654322
The origin of the global GWF of the EU-27 related to agricultural products consumption is shown in Figure 3. The figure shows that, although the external GWF of the EU-27 is spread across the rest of the world, its extent is concentrated in only a few places. Large external GWFs of the EU-27 can be observed in the Mississippi river basin (USA), Parana River basin (shared by Argentina, Bolivia, Brazil and Paraguay), Indus River basin, Southeastern China in the Huai He River basin, Yellow River basin and Yangtze River basin, and the Nile Delta in Egypt. The large GWF (both internal and external) is the function of the volume of the different crops supplied to the EU-27 and the GWF per unit of product (m3/tonne) of the consumed crops. Figure 3 shows a number of darker areas within the EU-27, which can be explained by the fact that about 59% of the GWF of consumption of the EU-27 is internal (within EU-27). The other darker area is in Eastern China, which is due to both the large volume of oil crops and cereals imported from China, and an excessive fertilizer application rate in China.
The contribution of different crops to the average GWF of consumption of the EU is presented in Figure 4. One quarter of the GWF relates to the consumption of oil crops (including cotton, soybean, oil palm, sunflower, rapeseed and others). The difference in the GWF among the different crops is mainly due to the rate of application of artificial fertilizers and manure, plus the crops′ nutrient removal rate. For example, cereal crops have larger N removal rates compared to oil seeds (e.g., 55 kg N/tonne of wheat compared to 10 kg N/tonne of oil crops). In addition, the fraction of N leaching from the agricultural field depends on the precipitation and the clay content of the soil as shown in Equation 1. The leaching fraction will be larger for areas with low clay soil and high precipitation. The large GWF of oil crops comes mainly from China, which has relatively larger fertilizer application rate and anthropogenic N and P emissions. The other products contributing a large share to the EU′s GWF of consumption are other cereals (14%), fruits (12%), wheat (11%), and vegetables (10%). Eastern China and developed Western countries (EU and USA) are known to have an excessive nutrient application rate, while the fertilizer application rate in developing countries such as Africa is not sufficient to supply the crop nutrient requirement [2]. This is reflected in the high N and P loads in Figure 1 and large GWF in Figure 3, in Eastern China and Western Europe.
Figure 3. The global GWF of the EU-27 related to consumption of agricultural products shown in mm/year on a 5 by 5 arc-min grid. Data per grid cell have been calculated as the WF within a grid cell (in cubic meters per year) divided by the area of the grid cell (in 103 m2). Period: 2007.
Figure 3. The global GWF of the EU-27 related to consumption of agricultural products shown in mm/year on a 5 by 5 arc-min grid. Data per grid cell have been calculated as the WF within a grid cell (in cubic meters per year) divided by the area of the grid cell (in 103 m2). Period: 2007.
Water 08 00030 g003
Figure 4. The contribution of different crops to the overall GWF of the EU consumption of agricultural products. Period: 2007.
Figure 4. The contribution of different crops to the overall GWF of the EU consumption of agricultural products. Period: 2007.
Water 08 00030 g004

3.4. Priority Basins and Products

A detailed analysis of the GWF of the EU-27 related to consumption of agricultural products shows that about 87% of the total GWF occurs in basins where the WPL is beyond 1.0, i.e., the ambient water quality standard are violated (Table 4). The hotspots with respect to the GWF of the EU-27′s consumption of agricultural products are identified based on the WPL in the respective river basins where these GWFs are located. Table 4 presents the top 20 river basins ranked by size of the GWF of the EU related to the consumption of agricultural products. The table also shows the crop products with significant contribution to the GWF within the basin. In all 20 river basins, the WPL exceeds 1.0, indicating that in all basins the water quality standards are violated.
Table 4. The water pollution level and contribution of major crops in the priority basins supporting the largest share of the total GWF related to EU′s consumption.
Table 4. The water pollution level and contribution of major crops in the priority basins supporting the largest share of the total GWF related to EU′s consumption.
BasinPercentage the Total GWF EU-27′s Consumption Located in the BasinWater Pollution Level for the Nutrient with Larger GWFProducts with Significant Contribution to the GWF in the Basin (% Contribution)
Yangtze (Chang Jiang)8.5%2.50oil crops-39%, paddy rice-14%, roots and tubers-12%, vegetables-10%, other cereals-8%, wheat-7%
Ebro6.9%1.08fruits-28%, vegetables-18%, other cereals-16%, pulses-14%, oil crops-6%, nuts-5%
Douro4.2%1.53fruits-24%, other cereals-20%, vegetables-17%, pulses-12%, nuts-5%, oil crops-5%
Huai He3.0%7.08oil crops-38%, wheat-19%, other cereals-16%, vegetables-7%, roots and tubers-7%, paddy rice-6%
Guadalquivir2.9%8.31oil crops-24%, fruits-21%, vegetables-15%, pulses-11%, other cereals-9%
Danube2.8%1.09other cereals-29%, Wheat-17%, oil crops-12%, fruits-10%, vegetables-9%, fodder crops-8%
Guadiana2.6%1.70fruits-24%, vegetables-16%, other cereals-15%, pulses-12%, oil crops-11%, nuts-6%
Tejo2.4%1.00fruits-26%, vegetables-17%, other cereals-16%, pulses-13%, oil crops-9%, nuts-6%
Mississippi2.1%1.33oil crops-97%
Xi Jiang2.0%2.49oil crops-29%, sugar crops-21%, vegetables-14%, paddy rice-13%, other cereals-9%, roots and tubers-7%
Rhine1.5%2.59oil crops-39%, other cereals-23%, sugar crops-11%, wheat-9%, roots and tubers-6%
Wisla1.5%3.33fruits-21%, oil crops-18%, wheat-15%, roots and tubers-14%, sugar crops-14%, other cereals-10%, vegetables-8%
Elbe1.4%2.72oil crops-58%, other cereals-16%, sugar crops-9%, wheat-6%, vegetables-5%
Oder1.4%2.14oil crops-37%, sugar crops-16%, fruits-13%, wheat-13%, roots and tubers-7%, other cereals-7%, vegetables-5%
Huang He (Yellow)1.2%2.87other cereals-25%, wheat-24%, oil crops-23%, vegetables-6%, fruits-7%, roots and tubers-7%
Amur1.1%1.18oil crops-34%, other cereals-33%, wheat-12%
Loire1.1%1.58other cereals-44%, oil crops-22%, pulses-19%, Fruits-7%
Indus1.1%2.99wheat-50%, oil crops-35%, sugar crops-6%, pulses-6%
Uruguay1.0%1.02oil crops-98%
Yongding He0.9%8.98other cereals-29%, wheat-26%, oil crops-21%, fruits-9%, vegetables-7%
The Yangtze River basin (China) and Ebro River basin (Spain) are two of the most important priority basins, contributing about 8.5% and 6.9% to the overall GWF of the EU-27′s consumption of agricultural products, respectively. In the Yangtze River basin, oil crops production accounts for a significant share of the GWF (39%), followed by rice (14%) and root & tubers (12%). In the Ebro River basin, fruits are the dominant crops that account for the largest share of the GWF in the basin (28%), followed by vegetables (18%) and other cereals (16%).
The Douro (shared by Spain and Portugal) and Huai He (China) are the other priority basins, where 4.2% and 3% of EU′s grey WF is located, respectively. The major crops accounting for a large share of the GWF in the Douro basin are fruits production (24%), other cereals (20%), and vegetables (17%). In the Huai He basins, oil crops (mainly cotton production) contribute about 38% to the GWF, followed by wheat (19%) and other cereals (16%). Next in the list of the priority basins are the Guadalquivir and the Danube, where 2.9% and 2.8% of the GWF of EU′s consumption are located, respectively. The GWF in the Guadalquivir river basin is mainly due to the production of oil crops, fruits, and vegetables. In the Danube basin, the major crops responsible for the large GWF are mainly other cereals, wheat, and oil crops production.
The results of this basin prioritization provide the EU-27 with a clear focus on two priority regions for crop production and sourcing: within the EU-27 itself and Eastern Asia (mainly China). Among the 20 priority basins, 11 of them are within the EU-27 and 6 mainly in China.
There is diversity in the priority crops among the basins. However, five crops and crop categories stand out as the major priority crops: oil crops, other cereals, fruits, wheat and vegetables. These five products also account for about 72% of the overall GWF of the EU-27′s consumption (Figure 4). Oil crops (including cotton) are the most important crops contributing the largest share of the EU-27 consumers GWF. They are listed as the priority crops in all the priority basins.

4. Discussion and Conclusions

In the current study, we estimate the global N and P emission related to crops production. The study also provides important information on the level of nutrient emission to water and the GWF of the EU-27, both related to the production and consumption of crops.
Our estimate of global total N emission of 31 million tonnes of N per year is 20% lower than the estimate by Bouwman et al. [24] (39 million tonne N/year) and 35% above the estimate by Liu et al. [23] (23 million tonne N/year). Our estimate of global P emission of 2.9 million tonne P per year is 28% lower than the estimate by Bouwman et al. [24] (4 million tonne P/year). It should be noted that the estimate by Bouwman et al. [24] includes nutrients flows from grassland, which is not included in the current study. Besides, the two studies refer to two different time periods—while the values in the current study refer to the year 2007, the Bouwman et al. [24] study refers to the year 2000.
We further compared our estimate of total N and P emission in from the EU-27′s croplands to freshwater with results from earlier studies. Our estimate of total N emission from EU-27′s croplands is much lower than the other studies (Table 5): It varies from as low as 30% of the estimate by Leip et al. [42] to almost half of the estimate by Westhoek et al. [43]. Further, differences in the modeling approach, the large difference between our estimate and the earlier studies, can be explained by the following facts: (a) our estimate refers to cropland, while the estimate from the other studies include nutrients flows from grassland; (b) the other studies have considered inputs through N deposition and fixation, while we consider only fertilizer and manure inputs in croplands; and (c) the study by Leip et al. [44] and Westhoek et al. [43] have considered N flows from other countries to EU through crop imports. On the other hand, our estimate of the total P emission is 87% larger than the estimate by Grizzetti et al. [45]. The estimate by Grizzetti et al. [45] refers to nutrient flows to European seas after accounting for nutrient retention in the river networks, while our estimate accounts for all the emission that goes to the freshwater systems. Although the volume of the estimated nutrients emission vary greatly among the different studies, the findings of all of the studies, including ours, underline the need to control nutrient emission in EU-27 and around the world to manage the undesired consequences on the aquatic ecosystems.
Table 5. Comparison of the Estimated Overall Diffuse N Load from the EU-27′S Croplands to Fresh Water with the Results from Previous Studies.
Table 5. Comparison of the Estimated Overall Diffuse N Load from the EU-27′S Croplands to Fresh Water with the Results from Previous Studies.
StudyN Leaching and Runoff to Fresh Water from Diffuse Sources (Million Tonne N/Year)P Leaching and Runoff to Fresh Water from Diffuse Sources (Million Tonne P/Year)Study Period
Leip et al. [42]5.98--2002
Liep et al. [44]5.0--2004
Westhoek et al. [43]3.3--2004
Grizzetti et al. [45]4.770.232000
Current study1.68 *0.432007
* Note: This refers to the total N leaching and runoff as estimated by Equation (1), which includes both natural and anthropogenic emissions.
The nutrient emission estimates are derived based on a large number of inputs that may include uncertainties and possible errors. Possible sources of uncertainties include: the fertilizer and manure application rate, the crops nutrient uptake rate, the maximum allowable nutrient concentration, and the natural nutrient concentration. For the GWF in particular, the maximum allowable concentration and the natural nutrient concentrations values taken will have large impacts on the final result. For example, Liu et al. [41] have shown that a slight increase/decrease results in a large change in the final values of the WPL. In a recent study, Mekonnen and Hoekstra [46] showed that the uncertainty of the N related global GWF due to uncertainty in input data and parameters ranges from −33% to +60%. Such large uncertainty ranges highlights the need for further detailed assessment of the nutrient emission and GWF based on local data and better modeling tools. In general, depending on the elements included and the data on the maximum allowable concentration and natural concertation used, the final volume of the GWF differs substantially. To address this and other criticism directed at the GWF concept, Vanham and Bidoglio [47] suggest further standardization of the GWF method.
The other potential uncertainties are related to the models used in estimating both the N and P emission per grid cell. The regression model of De Willigen [17] used here in estimating N emission is based on experimental data, and could therefore give unreasonably high leaching values when used for input parameters (precipitation, clay content, and layer thickness) beyond those used in deriving the regression equation. Liu et al. [23], for example, have set an upper limit on the total leached amount not to exceed 20% of the total mineral fertilizer application and 50% of the N from manure, arguing that the model may give unreasonably high leaching values for extreme precipitation. On the other hand, the model accounts for the effect of precipitation neglecting the effect of irrigation, potentially underestimating the N leaching from irrigated fields. The P emission is also estimated using a single factor derived from extensive modeling exercise applied at a global scale. Since it does not account for the local conditions (e.g., precipitation and soil P content) that might influence P runoff from agricultural field, the estimates at grid level are subject to large uncertainties and errors. Given the fact that 78% of the GWF is related to anthropogenic P emission, further refinement of the approach in modeling P emission will improve the final result. Given the uncertainties in the input variables and the modeling, the results in this paper should be taken as indicative of the nutrient loads and the grey WF magnitudes rather than as an absolute value.
The study also does not assess the final effect of eutrophication of lakes, estuaries and coastal seas associated with the emitted nutrients, as this is beyond the scope of the present paper. In future studies, the assessment of such effects on fish, the threat on biodiversity, on the tourism and recreation activities would add value to the analysis results.
The increases in international trade of agricultural products in the last few decades have also increased the globalization of virtual water flows as well as agricultural pollution. Consumers in one part of the world are shifting the environmental damage brought about by fertilizer application to regions that produce and export the products often without being aware of the environmental damages. The methodology of tracing global virtual water flows through an environmentally extended input-output framework such as the EXIOBASE has proven to be a valuable tool for providing insight into the environmental pressure put by consumers in one part of the word on regions that produce and export the product. By combining the data on direct and indirect water consumption per country with data on water pollution, it is possible to trace water pollution induced by final consumption.
This study is a major contribution to understanding the sustainability of the EU′s consumption with respect to water pollution both within the EU and the rest of the world. By tracing the global GWF of the EU, the study made the link between the crops consumed in the EU to the water pollution problem in other places. It further assessed the geographically specific extent to which the EU has externalized its environmental damage by identifying priority basins and products. Prioritizing or identifying the basins with serious water pollution and the associated products responsible for the large GWF can help the EU policy makers and governments to formulate policies and strategies on how to share responsibility for reducing GWFs where it is most needed.
As a consequence, analyses like the one presented in this paper may further encourage EU policy makers to develop GWF reduction targets that account for the assimilation capacity of river basins. These targets can be translated into maximum acceptable loads per river basin to be downscaled to the specific crops production. In many regions of the world, crops receive excessive amounts of nutrient. In these regions, the excessive application of nutrients can be reduced without affecting agricultural productivities [2]. Reducing the application of N and P in agriculture provides both environmental and economic benefits. The EU could implement a number of regulations including the certification of imported crop products in order to promote environmentally sustainable crop production in other countries.

Acknowledgments

The EXIOBASE database used in this research was developed through EU′s Sixth and Seventh Framework Programmes. This research was partly supported by the European Community FP7 project CREEA (Compiling and Refining Environmental and Economic Accounts), grant agreement No. 265134.

Author Contributions

Mesfin Mekonnen designed and performed research. Stephan Lutter and Aldo Martinez carried out the input-output modeling. Mesfin Mekonnen wrote the paper, and Stephan Lutter provided comments. All authors approved the final manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Galloway, J.N.; Cowling, E.B. Reactive nitrogen and the world: 200 Years of Change. AMBIO J. Hum. Environ. 2002, 31, 64–71. [Google Scholar] [CrossRef]
  2. Vitousek, P.M.; Naylor, R.; Crews, T.; David, M.B.; Drinkwater, L.E.; Holland, E.; Johnes, P.J.; Katzenberger, J.; Martinelli, L.A.; Matson, P.A.; et al. Nutrient imbalances in agricultural development. Science 2009, 324, 1519–1520. [Google Scholar] [CrossRef] [PubMed]
  3. Smith, V. Eutrophication of freshwater and coastal marine ecosystems a global problem. Environ. Sci. Pollut. Res. 2003, 10, 126–139. [Google Scholar] [CrossRef]
  4. Landsberg, J.H. The effects of harmful algal blooms on aquatic organisms. Rev. Fish. Sci. 2002, 10, 113–390. [Google Scholar] [CrossRef]
  5. Lewitus, A.J.; Horner, R.A.; Caron, D.A.; Garcia-Mendoza, E.; Hickey, B.M.; Hunter, M.; Huppert, D.D.; Kudela, R.M.; Langlois, G.W.; Largier, J.L.; et al. Harmful algal blooms along the north american west coast region: History, Trends, Causes, and Impacts. Harmful Algae 2012, 19, 133–159. [Google Scholar] [CrossRef] [Green Version]
  6. Anderson, D.; Glibert, P.; Burkholder, J. Harmful algal blooms and eutrophication: Nutrient Sources, Composition, and Consequences. Estuaries 2002, 25, 704–726. [Google Scholar] [CrossRef]
  7. O′Bannon, C.; Carr, J.; Seekell, D.A.; D′Odorico, P. Globalization of agricultural pollution due to international trade. Hydrol. Earth Syst. Sci. 2014, 18, 503–510. [Google Scholar] [CrossRef]
  8. Galloway, J.N.; Townsend, A.R.; Erisman, J.W.; Bekunda, M.; Cai, Z.; Freney, J.R.; Martinelli, L.A.; Seitzinger, S.P.; Sutton, M.A. Transformation of the nitrogen cycle: Recent Trends, Questions, and Potential Solutions. Science 2008, 320, 889–892. [Google Scholar] [CrossRef] [PubMed]
  9. Naylor, R.; Steinfeld, H.; Falcon, W.; Galloway, J.; Smil, V.; Bradford, E.; Alder, J.; Mooney, H. Losing the links between livestock and land. Science 2005, 310, 1621–1622. [Google Scholar] [CrossRef] [PubMed]
  10. Hoekstra, A.Y.; Chapagain, A.K. Globalization of Water: Sharing the Planet′s Freshwater Resources; Blackwell: Oxford, UK, 2008. [Google Scholar]
  11. Steen-Olsen, K.; Weinzettel, J.; Cranston, G.; Ercin, A.E.; Hertwich, E.G. Carbon, land, and water footprint accounts for the european union: Consumption, Production, and Displacements through International Trade. Environ. Sci. Technol. 2012, 46, 10883–10891. [Google Scholar] [CrossRef] [PubMed]
  12. Hoekstra, A.Y.; Mekonnen, M.M. The water footprint of humanity. Proc. Natl. Acad. Sci. USA 2012, 109, 3232–3237. [Google Scholar] [CrossRef] [PubMed]
  13. Wackernagel, M.; Onisto, L.; Bello, P.; Callejas Linares, A.; Susana López Falfán, I.; Méndez Garcı́a, J.; Isabel Suárez Guerrero, A.; Guadalupe Suárez Guerrero, M. National natural capital accounting with the ecological footprint concept. Ecol. Econ. 1999, 29, 375–390. [Google Scholar] [CrossRef]
  14. Leach, A.M.; Galloway, J.N.; Bleeker, A.; Erisman, J.W.; Kohn, R.; Kitzes, J. A nitrogen footprint model to help consumers understand their role in nitrogen losses to the environment. Environ. Dev. 2012, 1, 40–66. [Google Scholar] [CrossRef]
  15. Hoekstra, A.Y.; Chapagain, A.K.; Aldaya, M.M.; Mekonnen, M.M. The Water Footprint Assessment Manual: Setting the Global Standard; Earthscan: London, UK, 2011. [Google Scholar]
  16. Tukker, A.; Bulavskaya, T.; Giljum, S.; de Koning, A.; Lutter, S.; Simas, M.; Stadler, K.; Wood, R. The Global Resource Footprint of Nations. Carbon, Water, Land and Materials Embodied in Trade and Final Consumption Calculated with Exiobase 2.1.; The Netherlands Organisation for Applied Scientific Research/Leiden University/Vienna University of Economics and Business/Norwegian University of Science and Technology: Leiden/Delft, The Netherlands; Vienna, Austria; Trondheim, Norway, 2014. [Google Scholar]
  17. De Willigen, P. An Analysis of the Calculation of Leaching and Denitrification Losses as Practised in the Nutmon Approach; Plant Research International: Wageningen, The Netherlands, 2000; Volume 18. [Google Scholar]
  18. Batjes, N.H. Isric-Wise Derived Soil Properties on a 5 by 5 Arc-Min Global Grid (ver 1.2b); ISRIC—World Soil Information: Wageningen, The Netherlands, 2012. [Google Scholar]
  19. Mitchell, T.D.; Jones, P.D. An improved method of constructing a database of monthly climate observations and associated high-resolution grids. Int. J. Climatol. 2005, 25, 693–712. [Google Scholar] [CrossRef]
  20. Allen, R.G.; Pereira, L.S.; Raes, D.; Smith, M. Crop. Evapotranspiration: Guidline for Computing Crop Water Requirement; Food and Agriculture Organization: Rome, Italy, 1998. [Google Scholar]
  21. Smaling, E.M.A.; Roscoe, R.; Lesschen, J.P.; Bouwman, A.F.; Comunello, E. From forest to waste: Assessment of the Brazilian Soybean Chain, Using Nitrogen as a Marker. Agric. Ecosyst. Environ. 2008, 128, 185–197. [Google Scholar] [CrossRef]
  22. FAO. Scaling Soil Nutrient Balances; Food and Agriculture Organization: Rome, Italy, 2004; Volume 15. [Google Scholar]
  23. Liu, J.; You, L.; Amini, M.; Obersteiner, M.; Herrero, M.; Zehnder, A.J.B.; Yang, H. A high-resolution assessment on global nitrogen flows in cropland. Proc. Natl. Acad. Sci. USA 2010, 107, 8035–8040. [Google Scholar] [CrossRef] [PubMed]
  24. Bouwman, L.; Goldewijk, K.K.; van der Hoek, K.W.; Beusen, A.H.W.; van Vuuren, D.P.; Willems, J.; Rufino, M.C.; Stehfest, E. Exploring global changes in nitrogen and phosphorus cycles in agriculture induced by livestock production over the 1900–2050 period. Proc. Natl. Acad. Sci. USA 2013, 110, 20882–20887. [Google Scholar] [CrossRef] [PubMed]
  25. Seitzinger, S.P.; Mayorga, E.; Bouwman, A.F.; Kroeze, C.; Beusen, A.H.W.; Billen, G.; van Drecht, G.; Dumont, E.; Fekete, B.M.; Garnier, J.; et al. Global river nutrient export: A Scenario Analysis of Past and Future Trends. Glob. Biogeochem. Cycle 2010, 24, GB0A08. [Google Scholar] [CrossRef]
  26. Harrison, J.A.; Seitzinger, S.P.; Bouwman, A.F.; Caraco, N.F.; Beusen, A.H.W.; Vörösmarty, C.J. Dissolved inorganic phosphorus export to the coastal zone: Results from a Spatially Explicit, Global Model. Glob. Biogeochem. Cycle 2005, 19, GB4S03. [Google Scholar] [CrossRef]
  27. IFA; IFDC; IPI; PPI; FAO. Fertilizer Use by Crops; Food and Agriculture Organization: Rome, Italy, 2002. [Google Scholar]
  28. FAO. Fertistat on-Line Database; Food and Agriculture Organization: Rome, Italy, 2012. [Google Scholar]
  29. Heffer, P. Assessment of Fertilizer Use by Crop at the Global Level 2006/07–2007/08; International Fertilizer Industry Association: Paris, France, 2009. [Google Scholar]
  30. FAO. Faostat on-Line Database; Food and Agriculture Organization: Rome, Italy, 2012. [Google Scholar]
  31. FAO. Global Livestock Densities. Available online: http://www.fao.org/ag/againfo/resources/en/glw/GLW_dens.html (accessed on 20 January 2016).
  32. Sheldrick, W.; Keith Syers, J.; Lingard, J. Contribution of livestock excreta to nutrient balances. Nutr. Cycl. Agroecosyst. 2003, 66, 119–131. [Google Scholar] [CrossRef]
  33. Mekonnen, M.M.; Hoekstra, A.Y. The green, blue and grey water footprint of crops and derived crop products. Hydrol. Earth Syst. Sci. 2011, 15, 1577–1600. [Google Scholar] [CrossRef]
  34. IPNI. Ipni Estimates of Nutrient Uptake and Removal. International Plant Nutrition Institute. Available online: http://www.ipni.net/article/IPNI-3296 (accessed on 12 July 2012).
  35. Roy, R.N.; Finck, A.; Blair, G.J.; Tandon, H.L.S. Plant. Nutrition for Food Security a Guide for Integrated Nutrient Management; Food and Agriculture Organization: Rome, Italy, 2006; Volume 16. [Google Scholar]
  36. Ravindranath, N.H.; Somashekar, H.I.; Nagaraja, M.S.; Sudha, P.; Sangeetha, G.; Bhattacharya, S.C.; Abdul Salam, P. Assessment of sustainable non-plantation biomass resources potential for energy in india. Biomass Bioenergy 2005, 29, 178–190. [Google Scholar] [CrossRef]
  37. Graham, R.L.; Nelson, R.; Sheehan, J.; Perlack, R.D.; Wright, L.L. Current and potential U.S. Corn stover supplies. Agron. J. 2007, 99, 1–11. [Google Scholar] [CrossRef]
  38. Nelson, R.G. Resource assessment and removal analysis for corn stover and wheat straw in the eastern and midwestern united states—Rainfall and wind-induced soil erosion methodology. Biomass Bioenergy 2002, 22, 349–363. [Google Scholar] [CrossRef]
  39. Perlack, R.D.; Wright, L.L.; Turhollow, A.F.; Graham, R.L.; Stokes, B.J.; Erbach, D.C. Biomass as Feedstock for a Bioenergy and Bioproducts Industry: The Technical Feasibility of a Billion-Ton Annual Supply; DOE/GO-102005-2135ORNL/TM-2005/66; Oak Ridge National Laboratory, US Department of Energy: Oak Ridge, TN, USA, 2005. [Google Scholar]
  40. Krausmann, F.; Erb, K.H.; Gingrich, S.; Lauk, C.; Haberl, H. Global patterns of socioeconomic biomass flows in the year 2000: A Comprehensive Assessment of Supply, Consumption And Constraints. Ecol. Econ. 2008, 65, 471–487. [Google Scholar] [CrossRef]
  41. Liu, C.; Kroeze, C.; Hoekstra, A.Y.; Gerbens-Leenes, W. Past and future trends in grey water footprints of anthropogenic nitrogen and phosphorus inputs to major world rivers. Ecol. Indic. 2012, 18, 42–49. [Google Scholar] [CrossRef]
  42. Leip, A.; Achermann, B.; Billen, G.; Bleeker, A.; Bouwman, A.F.; de Vries, W.; Dragosits, U.; Döring, U.; Fernall, D.; Geupel, M.; et al. Integrating nitrogen fluxes at the european scale. In The European Nitrogen Assessment; Sutton, M.A., Howard, C.M., Erisman, J.W., Billen, G., Bleeker, A., Grennfelt, P., van Grinsven, H., Grizzetti, B., Eds.; Cambridge University Press: Cambridge, UK, 2011; pp. 345–376. [Google Scholar]
  43. Westhoek, H.; Lesschen, J.P.; Rood, T.; Wagner, S.; de Marco, A.; Murphy-Bokern, D.; Leip, A.; van Grinsven, H.; Sutton, M.A.; Oenema, O. Food choices, health and environment: Effects of Cutting Europe′s Meat and Dairy Intake. Glob. Environ. Chang. 2014, 26, 196–205. [Google Scholar] [CrossRef]
  44. Leip, A.; Weiss, F.; Lesschen, J.P.; Westhoek, H. The nitrogen footprint of food products in the european union. J. Agric. Sci. 2014, 152, 20–33. [Google Scholar] [CrossRef]
  45. Grizzetti, B.; Bouraoui, F.; Aloe, A. Changes of nitrogen and phosphorus loads to european seas. Glob. Change Biol. 2012, 18, 769–782. [Google Scholar] [CrossRef]
  46. Mekonnen, M.M.; Hoekstra, A.Y. Global gray water footprint and water pollution levels related to anthropogenic nitrogen loads to fresh water. Environ. Sci. Technol. 2015, 49, 12860–12868. [Google Scholar] [CrossRef] [PubMed]
  47. Vanham, D.; Bidoglio, G. A review on the indicator water footprint for the eu28. Ecol. Indic. 2013, 26, 61–75. [Google Scholar] [CrossRef]

Share and Cite

MDPI and ACS Style

Mekonnen, M.M.; Lutter, S.; Martinez, A. Anthropogenic Nitrogen and Phosphorus Emissions and Related Grey Water Footprints Caused by EU-27′s Crop Production and Consumption. Water 2016, 8, 30. https://doi.org/10.3390/w8010030

AMA Style

Mekonnen MM, Lutter S, Martinez A. Anthropogenic Nitrogen and Phosphorus Emissions and Related Grey Water Footprints Caused by EU-27′s Crop Production and Consumption. Water. 2016; 8(1):30. https://doi.org/10.3390/w8010030

Chicago/Turabian Style

Mekonnen, Mesfin M., Stephan Lutter, and Aldo Martinez. 2016. "Anthropogenic Nitrogen and Phosphorus Emissions and Related Grey Water Footprints Caused by EU-27′s Crop Production and Consumption" Water 8, no. 1: 30. https://doi.org/10.3390/w8010030

APA Style

Mekonnen, M. M., Lutter, S., & Martinez, A. (2016). Anthropogenic Nitrogen and Phosphorus Emissions and Related Grey Water Footprints Caused by EU-27′s Crop Production and Consumption. Water, 8(1), 30. https://doi.org/10.3390/w8010030

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop