The Water Implications of Greenhouse Gas Mitigation : Effects on Land Use , Land Use Change , and Forestry

This study addresses the water quantity and quality implications of greenhouse gas mitigation efforts in agriculture and forestry. This is done both through a literature review and a case study. The case study is set in the Missouri River Basin (MRB) and involves integration of a water hydrology model and a land use model with an econometric model estimated to make the link. The hydrology model (Soil and Water Assessment Tool, SWAT) is used to generate a multiyear, multilocation dataset that gives estimated water quantity and quality measures dependent on land use. In turn, those data are used in estimating a quantile regression model linking water quantity and quality with climate and land use. Additionally, a land use model (Forest and Agricultural Sector Optimization Model with Greenhouse Gases, FASOMGHG) is used to simulate the extent of mitigation strategy adoption and land use implications under alternative carbon prices. Then, the land use results and climate change forecasts are input to the econometric model and water quantity/quality projections developed. The econometric results show that land use patterns have significant influences on water quantity. Specifically, an increase in grassland significantly decreases water quantity, with forestry having mixed effects. At relatively high quantiles, land use changes from cropped land to grassland reduce water yield, while switching from cropping or grassland to forest yields more water. It also shows that an increase in cropped land use significantly degrades water quality at the 50% quantile and moving from cropped land to either forest or pasture slightly improves water quality at the 50% quantile but significantly worsens water quality at the 90% quantile. In turn, a simulation exercise shows that water quantity slightly increases under mitigation activity stimulated by lower carbon prices but significantly decreases under higher carbon prices. For water quality, when carbon prices are low, water quality is degraded under most mitigation alternatives but quality improves under higher carbon prices.


Introduction
The Intergovernmental Panel on Climate Change (IPCC) [1] indicates greenhouse gas emissions (GHGEs) are a main driver of climate change, and the agriculture and forestry (AF) sector can alter operations to reduce net emissions [2][3][4][5][6][7][8].The logic for employing AF mitigation involves, among other things, the size of sectoral GHGEs.Estimates indicate AF is the source of 30% of global GHGEs predominately through emissions of methane, nitrous oxide, and carbon dioxide [9].Additionally, an estimated 25% of historical carbon releases have come from the sectors mainly in the form of lost sequestration [10].In terms of mitigation, AF can pursue several major types of mitigation actions.First, AF may manipulate enterprise management to reduce cropping-, livestock-, or forest-based emissions.Secondly, AF may enhance sequestration by creating or expanding land-based sinks, retaining more of the carbon that cycling in and out of AF each year.This is done by reducing tillage intensity, altering land use towards less disturbed regimes, reducing deforestation, enhancing forest management, and pursuing afforestation [11,12].Thirdly, AF may produce substitute products for emission-intensive products like fossil fuels, steel and concrete, displacing emissions involved with their manufacture and use [13].Finally, AF may develop and utilize technical advances that increase yields while not increasing GHGEs, thus reducing emissions per unit produced and allowing less land to be used to produce a given amount and perhaps less input use [14].
The abovementioned AF mitigation measures have potential co-benefits or adverse side-effects in several dimensions such as institutional, social, economic, or environmental sides [8,15,16].Specifically, on the environmental side, AF mitigation measures can influence land degradation/restoration, biodiversity, soil quality, and water resources.Here, we focus on water impacts.On balance, a comprehensive study of water quantity and quality effects of a wide array of mitigation possibilities is not available.We attempt to fill that gap.
A number of studies have called for comprehensive examinations of the water effects of mitigation activities.Gupta et al. [17] indicated that many empirical studies of agricultural impacts on water did not account for the fundamental principles of soil water storage, water infiltration, and surface runoff.Medhi et al. [18] argued that combined climate and land use change scenarios should be considered to analyze the impacts of mitigation alternatives on water.Jackson et al. [15] indicate that water quantity effects usually occur through direct alterations in irrigation water use plus alterations in water run off or groundwater infiltration.In general, water quality effects occur when agriculture, forestry, and other land use mitigation strategies alter erosion rates, input usage, and animal manure quantity, in turn altering runoff and infiltration of sedimentation, manure, and chemicals.Others have shown water quality can be improved by use of AF GHGE-mitigating practices like adoption of no till or buffer strips [19,20].
This study does a comprehensive water quality and quantity evaluation of a set of mitigation possibilities in a case study setting.To do this, we employ a hydrological model to simulate water processes coupled with an econometric synthesis and then a land use, economic sector model to incorporate land use change influences.The study is done in the context of the Missouri River Basin (MRB).

Literature Based Findings on Mitigation and Water
The following discusses the water quantity and quality implications of AF mitigation strategies under four broad categories: AF management, land use change, bioenergy, and technological progress.

AF Management
Crop management alternatives for mitigation involve manipulation of tillage, mix of crops grown, irrigation practices, cover crop usage, and fertilization among other possibilities (for an extensive list, see [3,4,21]).Water quantity is affected through: (a) changing tillage practices to retain more residue, which can increase sequestration while also reducing runoff and possibly aquifer infiltration [22]; (b) altering crop mix, like shifting from rice to row crops, which can reduce water diversion, consumption, and evaporation [23][24][25]; and (c) improving irrigation efficiency, which reduces water diversion and evaporative losses plus reduces runoff [26,27].Also, recharge of groundwater may be increased by additional use of deep tillage [28] but can be reduced because of less intensive conservation tillage [29].Additionally, when decreasing input intensity, many studies have found that this combined with climate change is likely to decrease water yield but increase water quality [18,30,31].
Mitigation strategies also alter water quality particularly when runoff/infiltration quantity and chemical content are reduced by tillage shifts, reduced fertilization, or conservation practice usage [19,[32][33][34][35][36].For example, Honisch et al. [19] monitored water quality under sustainable farming practices in Bavaria, Germany, finding significant reductions of N loads and phosphate loads after 4 years.
Animal-related mitigation approaches include manure management through digesters and land application, animal stocking intensity/herd size reduction, species choice, grazing land management, and feeding practice alteration.Water quality is affected by use of these practices through changes in runoff of manure, sediments, and nutrients plus indirect cropping effects of altered feed demands.Water quantity is affected by manure handling choice, where systems involving washing increase quantity used and runoff.Fast, mechanical removal of manure solids not only reduces water use but also improves water quality [37,38].Land-based application of manure as a fertilizer source can increase nutrient runoff and degrade water quality [39,40].Manipulation of animal stocking intensity/herd size alteration alters manure loads and land runoff potential (by altering vegetative cover), in turn altering water infiltration/runoff and nutrient/sediment runoff.Switches to non ruminant species as opposed to ruminant animal species reduces enteric fermentation while altering feed consumption and the amount of manure produced per unit of animal protein.This in turn affects water quantity and water quality via the means discussed above [41].Alterations in management of extensive grazing lands can be done via stocking rate reduction, changes in fertilization, improved fire management, reduced brush incidence, and altered plant species mix, all of which in turn alter water infiltration and runoff plus associated sediments and chemicals reaching ground and surface waters.
Forest management mitigation strategies include afforestation, rotation length extension, improved silviculture, and reduced fire incidence, as discussed in Murray et al. [12], Nabuurs et al. [42], and Wall [43].In terms of water quality, afforestation of crop, pasture, or other land, particularly in the cropland case, could reduce the usage of chemical inputs and soil disturbance, thus reducing chemical runoff and sedimentation and in turn altering water quality [44][45][46].However, this generally increases the amount of water consumed, as trees use more than grass and crops in many cases [15].Altering rotation length would reduce sedimentation from soils disturbed by harvesting and alter water quantity through reductions in runoff but increase in vegetative use.Fire management by thinning, removal of dead materials, and controlled burns increase sediment runoff due to soil disturbance but then reduces runoff of ash, nutrients, and sediments when fires occur.

Land Use Change
Mitigation actions being undertaken to reduce GHGEs would change land use systems [47][48][49][50] and in turn alter water resources.Land use change in an agricultural context involves transformations between cropping, pasture/range, forest, wetlands, or urban spaces.For mitigation, this generally involves movements to less intensive usages, with cropland moving to pasture/range, forestry or wetlands plus pasture moving to forest or wetlands.
Land use change affects water use and water quality with the effect dependent on the original and final land uses along with irrigation status, vegetation mix, and production practices/input use.A number of studies have addressed the water quantity implications of land use change.Leterme and Mallants [51] show that groundwater recharge will increase when cropland is used for maize but will decrease if land is converted to pasture or forest, implying the more complete the land cover, the less the infiltration.Bhardwaj et al. [52] suggest that cropland movement into energy crops will increase water use.Water runoff is also altered when croplands are moved to forests, with runoff predominantly decreasing, indicating that larger plants use more water [15].Frankenberger [53] examined how runoff is affected by land use in Indiana, finding on fields planted to corn or soybeans that with a 4-in rain, runoff volume fell from 97% to between 12% to 30% when the land was altered to forest, pasture, or turf grass.This is likely due to a combination of these land uses slowing runoff, possibly stimulating more infiltration and increased vegetative water use.
Water quality is also affected by land use [54][55][56][57][58][59].Erosion is a major force and is altered by mitigation practices that alter tillage intensity or soil disturbance [8,[60][61][62][63].For example, erosion associated with the loss of soil carbon can be reduced by windbreaks, grassland replacement, and afforestation [8,62].Furthermore, water quality is impacted by rates of chemical or animal manure fertilization [64][65][66].It is also altered by practices that alter water infiltration, like buffer strips [19,67,68].The impacts on water quality are varied.For example, Pattanayak et al. [69] and Townsend et al. [70] showed afforestation leading to increased water quality, while Jackson et al. [15] reviewed cases exhibiting the opposite.Thus, impacts on water quantity and quality differ by situation and prior land use/management.

Bioenergy
Bioenergy is a mitigation alternative with biogenic feedstocks displacing fossil fuel use [2,13].The same would be true for water, where the direct effects involve water used in processing and production.The Iowa Department of Natural Resources [71] showed quantity can be an issue, with aquifers in Iowa being drawn down by 17.1 ft over 10 years when ethanol facilities were installed.
Quality is also affected as bioenergy feedstock sources include conventional crops and their residues, energy crops, and animal wastes.Raising or recovering such feedstocks can alter runoff of chemical fertilizer, pesticides, and erosion along with vegetative water use and runoff volume [72].

Technological Progress
Technological progress is another mitigation strategy and can also influence water quantity and quality.This is particularly true when: (1) yield increases are achieved with approximately the same inputs as equal amounts that can be produced with smaller acreage and thus less inputs per unit product; and/or when (2) the same product can be produced with less inputs reducing chemical runoff.For example, recently US corn yields have risen without increasing nitrogen fertilization per acre due to such things as improvements in genetics, pest resistance, drought tolerance, and nutrient efficiency/application practices [14].
Technology can also have quantity impacts when the amount of water used per acre or the land area of crops changes, although this can be either positive or negative.Quality impacts occur when the technological developments induce less input use or crop mix/management shifts that lower chemical input usage and erosion.Baker et al. [14] examined technology effects in a US setting, showing that a continuation of recent trends leads to both reduced GHG emissions and water quality improvements.However, the interrelationship is complex, with many possibilities having different water implications.

Empirical Investigation on Mitigation and Water
This study aims to examine the water implications of mitigation strategy use via a four-step process.First, we run a hydrological model (Soil and Water Assessment Tool, SWAT) [73,74] across a broad study region generating results on water quantity and quality.Second, we take the SWAT results and estimate summary equations over them that indicate how land use alters water quantity and quality attributes.Third, we use a land use, economic sector model (the Forest and Agricultural Sector Optimization Model with Greenhouse Gases, hereafter abbreviated as FASOMGHG) [75,76] to simulate land use changes under different mitigation alternatives.Fourth, we evaluate the SWAT summary functions with the FASOMGHG land use changes to examine the composite effects on water.Details on the study area, SWAT use, and climate data follow.

Study Area
The MRB (Figure 1) is the largest river basin in the continental United States and includes all of Nebraska and parts of Colorado, Iowa, Kansas, Minnesota, Missouri, Montana, North Dakota, South Dakota, and Wyoming, covering a total of 411 counties.The United States Department of Agriculture (USDA) [77] estimates that 29% of the basin land is cropped, 49% is for grazing, 9% is forest, 10% is a mix of permanent pasture, hay land, water, wetland, horticulture, and barren land, and the remaining 3% is urban areas.Most of the grazing land is located in the western and central parts of the basin, while most of the forestland is located in west and central Missouri.
The 2007 Census of Agriculture reported that the MRB produced about $49 billion in agricultural sales, 45% of which is from crops and 55% from livestock.The principal crops grown in the basin are corn, soybeans (mainly in the eastern portion of the basin), and wheat (mainly in the western portion of the basin).Cow calf and feedlot production are the primary livestock enterprises.The USDA [77] indicates that 65% of the cropland was irrigated and 76% of both cropland and pasture were treated with pesticides.
Sustainability 2018, 10, x FOR PEER REVIEW 5 of 22 and the remaining 3% is urban areas.Most of the grazing land is located in the western and central parts of the basin, while most of the forestland is located in west and central Missouri.
The 2007 Census of Agriculture reported that the MRB produced about $49 billion in agricultural sales, 45% of which is from crops and 55% from livestock.The principal crops grown in the basin are corn, soybeans (mainly in the eastern portion of the basin), and wheat (mainly in the western portion of the basin).Cow calf and feedlot production are the primary livestock enterprises.The USDA [77] indicates that 65% of the cropland was irrigated and 76% of both cropland and pasture were treated with pesticides.

SWAT Generated Data
To do the estimation phase in this analysis, water information related to land use in MRB is needed.We generated such data using the SWAT model, as it considers many hydrologic processes plus climate, land use, and management practices.SWAT is a comprehensive semi-distributed and physically based model that simulates water quantity and quality within catchments in hydrologic response units (HURs) at a daily time-step [78].
The SWAT input and output data contains records on water runoff, water quality, land use, climate, irrigation water use, and land use change.SWAT input data covers hydrography, terrain, land use, soil, soil drainage-tile, weather, and management practices [74].SWAT output contains simulated water runoff and water quality indicators on a monthly basis for the 13,437 sub-basins in the MRB.SWAT was run for the 1990-2010 period and is well-calibrated and validated in this study area [79,80].
One estimation focus is water quantity and quality effects.For water quality, following Cude [81], we selected two quality indicators: total nitrogen, including ammonia, nitrate, and nitrite in surface runoff, and total phosphorus in surface runoff, which were rated from 10 to 100.Cude [81] in turn uses a subindex (SI) to convert the water quality indicators into a relative quality rating.The formulae for the total nitrogen and total phosphorus subindices are  Total Nitrogen (N)  Total Phosphorus (P)

SWAT Generated Data
To do the estimation phase in this analysis, water information related to land use in MRB is needed.We generated such data using the SWAT model, as it considers many hydrologic processes plus climate, land use, and management practices.SWAT is a comprehensive semi-distributed and physically based model that simulates water quantity and quality within catchments in hydrologic response units (HURs) at a daily time-step [78].
The SWAT input and output data contains records on water runoff, water quality, land use, climate, irrigation water use, and land use change.SWAT input data covers hydrography, terrain, land use, soil, soil drainage-tile, weather, and management practices [74].SWAT output contains simulated water runoff and water quality indicators on a monthly basis for the 13,437 sub-basins in the MRB.SWAT was run for the 1990-2010 period and is well-calibrated and validated in this study area [79,80].
One estimation focus is water quantity and quality effects.For water quality, following Cude [81], we selected two quality indicators: total nitrogen, including ammonia, nitrate, and nitrite in surface runoff, and total phosphorus in surface runoff, which were rated from 10 to 100.Cude [81] in turn uses a subindex (SI) to convert the water quality indicators into a relative quality rating.The formulae for the total nitrogen and total phosphorus subindices are • Total Phosphorus (P) where SI N and SI P are the subindices for the nitrogen-and phosphorus-related measures.A single water quality index (WQI) is then formed from these subindices by using the unweighted harmonic square mean, as in Swamee and Tyagi [82]: where I is the water quality index, M is the number of subindices, and S i is the ith subindex.
The SWAT output describes conditions by sub-basin.For the estimation, we aggregated the results to the county level.In doing this, when sub-basins were distributed across several counties, we did a weighted average based on the proportion of each sub-basin that falls within each county.
Land use is categorized in SWAT into the following categories: (1) continuous dryland crops, (2) irrigated crops, (3) crop rotations, (4) alfalfa and hay, (5) evergreen forest and deciduous forest, (6) pasture/range, and (7) urban area.In the SWAT simulations, land use is assumed to be time invariant, and Figure 2 portrays the incidence of each land use across the MRB.Although land use does not change over time, the water runoff data was calibrated to real values.In our analysis, we will mainly looked at how land use changes between cropland, grassland, and forest land affect water quality and quantity, and thus we aggregated the cropped lands (categories one to four above) into one broad category: cropped land.where N SI and P SI are the subindices for the nitrogen-and phosphorus-related measures.A single water quality index (WQI) is then formed from these subindices by using the unweighted harmonic square mean, as in Swamee and Tyagi [82]: where I is the water quality index, M is the number of subindices, and i S is the ith subindex.
The SWAT output describes conditions by sub-basin.For the estimation, we aggregated the results to the county level.In doing this, when sub-basins were distributed across several counties, we did a weighted average based on the proportion of each sub-basin that falls within each county.
Land use is categorized in SWAT into the following categories: (1) continuous dryland crops, (2) irrigated crops, (3) crop rotations, (4) alfalfa and hay, (5) evergreen forest and deciduous forest, (6) pasture/range, and (7) urban area.In the SWAT simulations, land use is assumed to be time invariant, and Figure 2 portrays the incidence of each land use across the MRB.Although land use does not change over time, the water runoff data was calibrated to real values.In our analysis, we will mainly looked at how land use changes between cropland, grassland, and forest land affect water quality and quantity, and thus we aggregated the cropped lands (categories one to four above) into one broad category: cropped land.

Climatic Data
Monthly climate data summarizing temperature and precipitation averages and extremes were drawn from the National Oceanic and Atmospheric Administration (NOAA) Satellite and Information Service, National Climatic Data Center (NCDC) for the period 1990-2010 [83].Those data include: (1) number of days with precipitation greater than or equal to 0.1, 0.5, and 1.0 in, respectively; (2) number of days with minimum temperature less than or equal to 0.0 and 32.0 • F, respectively; (3) number of days with maximum temperature greater than or equal to 90.0 • F; (4) total precipitation measured in millimeters; and (5) monthly mean temperature measured in • F. The NOAA data often contains multiple weather stations in a county and these were averaged across all stations in that county to form the dataset.We also considered the influences of the El Niño Southern Oscillation phase, which is categorized into three phases-warm phase (El Niño), cold phase (La Niña), and neutral phase-as reported by the Japan Meteorological Agency, Tokyo Climate Center [84].In the period 1990-2010, there were five warm phase years (1991, 1997, 2002, 2006, and 2009) and four cold phase years (1998, 1999, 2007, and 2010).

Quantile Regression over Panel Data
Quantile regression over panel data [85] was used to summarize SWAT results on how land use affects water runoff quantity and surface water quality.Quantile regression lets us estimate the median and other break points.We estimated 10%, 25%, 50% (median), 75%, and 90% quantiles.The resultant estimation for say the 90% quantile gives a level of water quantity that is exceeded only 10% of the time, with 90% of the observations being equal or smaller.
Quantile regression for panel data following Koenker [85] involves estimating the function where y it denotes water quantity or quality measured in region i at time t and X T it is a vector of independent variables that influence the water quantity and quality (including land use, climate factors, and water quantity (for the quality estimation).β τ are parameters to be estimated, c i is the regional fixed effect, and u it is the error term.In terms of the quantiles given, we are estimating quantile τ, and the conditional τth quantile of y it is Following Koenker [85], we can estimate Model (5) for 10 quantiles simultaneously by solving the following minimization problem: where w k is the probability weight controlling the relative influence of the associated quantiles τ k and ρ τk denotes the piecewise linear quantile loss function.Koenker [85] named Problem (6) the penalized quantile regression with fixed effects approach, and we can obtain the fixed effect estimators while λ→0.

FASOMGHG Model
A national agricultural land use, sector model was used to generate information on mitigation strategies adopted and their land use consequences under alternative mitigation incentives in the form of carbon prices.The model used was FASOMGHG [75,76,86,87], which is a dynamic, nonlinear, and price endogenous programming model of the US forest and agricultural sectors that simulates forest and agricultural land allocation and management over time in a perfectly competitive set of markets.This model represents: forest, crop, and pasture/range land use; agricultural crop and livestock production; livestock feeding; agricultural processing; log production; forest processing; carbon sequestration; CO 2 /non-CO 2 GHG emissions; wood product markets; agricultural markets; and GHG payments [75].FASOMGHG simulates intertemporal factor and commodity market equilibria which are imposed by the first order conditions resulting from maximizing intertemporal economic welfare.More on the mathematical structure of the model is given in Appendix A and in Adams et al. [75].
In terms of this project, FASOMGHG was run under a range of prices paid per metric ton of GHGEs reduced, sequestered, or offset.In turn, FASOMGHG chose a varying portfolio of agricultural GHG mitigation possibilities composed of manipulations in land use, AF management, and biofuel production across a variety of explicitly modeled sequestration, emission reduction, and biofuels-related possibilities, as done previously by McCarl and Schneider [2], Lee [88], Lee et al. [89], and Murray et al. [12].

Estimation Results
As stated above, summary function regressions were estimated over the SWAT county-level results arising from a 1990-2010 simulation.Table 1 reports descriptive statistics on those SWAT water results, the SWAT input data, and the regional climate data.Table 2 reports quantiles in that data for the water quantity yield and water quality.The average water quantity is around 7.2 mm per sub-basin per month, which is close to the value reported at the 75% quantile, meaning that 75% of the SWAT-generated MRB water quantities are below average.On the other hand, the average water quality index is around 19 and we find that 88% of the observations exhibit a lower quality index than that with about half having the worst quality index value (WQI = 10).This shows the importance of examining the whole distribution and is why we used quantile regression.To begin examination of the interrelationship between land use and water quality, we partitioned the data into different classes based on the water quality index, as shown in Table 3. There, we see that the counties exhibiting the worst water quality have 53.8% of lands being cropped, while those with better water quality have 31.8% of the land cropped.On the other hand, grass coverage in areas exhibiting better water quality averages 40.2%, while it averages 15.7% in areas with the worst water quality.This indicates that land with more grass coverage as opposed to crops tends to have higher water quality, which is not surprising since nitrogen and phosphorus typically are applied on croplands.We also observe from Table 3 that forest coverage appears to be associated with a slight increase in water quality at the low quantiles but not at the highest forest coverage shares.Furthermore, Table 3 indicates urban coverage worsens water quality.These relationships will be further examined in the econometric analysis.Source: the SWAT model (2013).Note: 1. Water quality is an index conducted by two indicators, total nitrogen and total phosphorus, and scaled from 10 (the worst water quality) to 100 (the best water quality).

Quantile Regression Results for Water Quantity
The quantile regressions are estimated using the R package rqpd [85,90].Results relevant to the water quantity influence of land use and climate are reported in Table 4.
Land use patterns have significant influences on water quantity across all quantiles.The proportion of land used for crops and forests increases water quantity, but the incremental rate falls as the share of cropland increases, which is in conflict with Molina-Navarro et al. [91].We also find grass coverage has negative influences on water quantity, with positive effects of squared terms at all quantiles.Urban land positively affects water quantity but only at lower quantiles.We further calculated land use change effects on water yield.Figures 3-5 depict the effects of land reallocations from cropped usage to grassland, from cropped usage to forest, and from grassland to forest, respectively, where other land shares and climate are held constant.The water yield baseline for computing land switch effects is the total water yield under the land shares of 43% cropped land, 28% grassland, and 8% forest-the mean values reported in Table 1.
yield baseline for computing land switch effects is the total water yield under the land shares of 43% cropped land, 28% grassland, and 8% forest-the mean values reported in Table 1.As land use changes from cropped land to grassland (Figure 3), water quantity reduces at all quantiles, with lower effects at higher quantiles as also found in [51,53].Namely, when the yield baseline for computing land switch effects is the total water yield under the land shares of 43% cropped land, 28% grassland, and 8% forest-the mean values reported in Table 1.As land use changes from cropped land to grassland (Figure 3), water quantity reduces at all quantiles, with lower effects at higher quantiles as also found in [51,53].Namely, when the As land use changes from cropped land to grassland (Figure 3), water quantity reduces at all quantiles, with lower effects at higher quantiles as also found in [51,53].Namely, when the proportion of cropped land falls from 43% to 32% with a corresponding increase in grassland, then water yield decreases by −1.5 mm at the 75% and 90% quantiles, −1.3 mm at the 50% quantile, −0.4 mm quantile at the 25% quantile, and −0.2 mm at the 10% quantile.It implies that grasslands reduce runoff particularly at relatively high quantiles of water yield.This is likely due to more complete land coverage slowing runoff and increased infiltration.
Sustainability 2018, 10, x FOR PEER REVIEW 12 of 22 proportion of cropped land falls from 43% to 32% with a corresponding increase in grassland, then water yield decreases by −1.5 mm at the 75% and 90% quantiles, −1.3 mm at the 50% quantile, −0.4 mm quantile at the 25% quantile, and −0.2 mm at the 10% quantile.It implies that grasslands reduce runoff particularly at relatively high quantiles of water yield.This is likely due to more complete land coverage slowing runoff and increased infiltration.
Figure 5.The total effects on water quantity when switching land from grass to forest.
On the other hand, switching land use from cropping to forest (Figure 4) significantly increases water quantity.The result is different from the findings in [15,53].Moreover, the land switch effects on water yield from cropping to forest takes the form of an inverse U at all quantiles, implying that adding forest land increases water yield up to a critical point then decreases it.Switching from cropping to forest yields more water at the higher quantiles.We also find the same pattern of water yield when grass land switches to forest (Figure 5).
Across these results, afforestation from either cropped land or grassland has greater water quantity benefits than does the switch from cropped land to grassland.
On the climate side, the results show, unsurprisingly, that water quantity is principally influenced by precipitation, with increases as conditions get wetter.Also, a measure of precipitation extremes (the number of days with greater than 1 in of precipitation) is also associated with increased quantity, a result that conflicts with the findings in Zhang et al. [92].Average temperature, on the other hand, does not have significant effects across most quantile equations but extreme temperature events do.The number of cold days (≤32 °F ) has a positive effect on water quantity, while the number of hot days (≥90 °F ) has a negative impact.
On El Niño and La Niña events, we find El Niño reduces MRB water quantity while La Niña increases it at all quantiles.This corresponds with NOAA findings that MRB winters are warmer and drier when El Niño occurs, while La Niña leads to cooler and wetter winters [93].The results also show that both El Niño and La Niña have greater effects on the higher quantiles.

Quantile Regression Results for Water Quality
Water quality quantile regressions are given in Table 5, although we had to drop the 10% quantile since almost all of observations in that quantile exhibit the lowest water quality value of 10.On the other hand, switching land use from cropping to forest (Figure 4) significantly increases water quantity.The result is different from the findings in [15,53].Moreover, the land switch effects on water yield from cropping to forest takes the form of an inverse U at all quantiles, implying that adding forest land increases water yield up to a critical point then decreases it.Switching from cropping to forest yields more water at the higher quantiles.We also find the same pattern of water yield when grass land switches to forest (Figure 5).
Across these results, afforestation from either cropped land or grassland has greater water quantity benefits than does the switch from cropped land to grassland.
On the climate side, the results show, unsurprisingly, that water quantity is principally influenced by precipitation, with increases as conditions get wetter.Also, a measure of precipitation extremes (the number of days with greater than 1 in of precipitation) is also associated with increased quantity, a result that conflicts with the findings in Zhang et al. [92].Average temperature, on the other hand, does not have significant effects across most quantile equations but extreme temperature events do.The number of cold days (≤32 • F) has a positive effect on water quantity, while the number of hot days (≥90 • F) has a negative impact.
On El Niño and La Niña events, we find El Niño reduces MRB water quantity while La Niña increases it at all quantiles.This corresponds with NOAA findings that MRB winters are warmer and drier when El Niño occurs, while La Niña leads to cooler and wetter winters [93].The results also show that both El Niño and La Niña have greater effects on the higher quantiles.

Quantile Regression Results for Water Quality
Water quality quantile regressions are given in Table 5, although we had to drop the 10% quantile since almost all of observations in that quantile exhibit the lowest water quality value of 10.Additionally, most of the estimated coefficients are insignificant in the 25% and 75% quantiles, so we mainly discuss the land use impacts on the 50% and 90% quantiles.Previous studies have found that agricultural land use has impacts on water quality [94][95][96], while Sliva and Williams [97] found the opposite.In our analysis, the results show that moving from cropped land to either forest or pasture will slightly improve water quality at the 50% quantile but worsen water quality at the 90% quantile.An increase in cropped land use significantly degrades water quality at the 50% quantile.Previous studies have found that water quality is strongly influenced by agricultural land use and land use change [44,59,91,95,96,98].For example, Lee et al. [59] found that the proportions of land use of urban, agricultural, and forest are correlated with water quality in South Korea.Gunawardhana et al. [98] found that converting forested land to agricultural land degrades water quality in the Uma Oya catchment area in Sri Lanka.Miserendino et al. [99] and Tu [100] found that water quality increases as forest land increases, perhaps because of lower inorganic ions in forested land use [56,58,97,101].Urban area does not have significant impacts on water quality, which conflicts with the results from the study of DeFries and Eshleman [57], Miserendino et al. [99], and Ahearn et al. [102].
The effect of water quantity on water quality is positive at the 90% quantile but negative at the 75% quantile.The quantity effect on water quality at the 75% quantile is negative and similar to results at the 50% quantile, and it implies that increasing water quantity will improve the higher part of the water quality distribution (WQI > 90% quantile) but worsen lower water quality below that (WQI < 75% quantile), perhaps due to runoff.Our result is consistent with Molina-Navarro et al. [91] at higher quantiles in showing that water quality is closely linked with water quantity.
Next, we discuss the impacts of climatic factors.Notice that the impacts depend on current conditions since the estimated functions include squared terms for some climate variables.The result shows that water quality will worsen with lower precipitation levels and is improved as precipitation increases.Extreme precipitation has positive impacts on water quality probably since precipitation has dilution and flushing effects.On the other hand, extreme temperature degrades water quality.This is likely because extreme lower temperatures freeze water and slow down flows while extreme high temperatures cause higher evaporation or transpiration, in turn reducing water flows and resultant dilution.
The marginal effects from El Niño and La Niña have opposite results between the 50% and 90% quantiles, with the absolute magnitude at the 90% quantile being larger than that at the 50% quantile.Compared with the neutral phase, El Niño and La Niña significantly improve water quality at the 90% quantile, while both El Niño and La Niña have significant negative effects on water quality at the 50% quantile.

Analysis of Water Implications of Mitigation Strategy Choice
We now turn to examine the effects of mitigation efforts on water quantity and quality.We do this by using FASOMGHG model under different carbon prices to give us changes in: (1) afforestation; (2) crop fertilization; (3) crop tillage intensification; (4) crop management; (5) livestock enteric and manure management; (6) bioenergy feedstock incidence; and (7) forest management.The afforestation mitigation strategy converts crop and/or grasslands to forest, and the forest management strategy includes lengthened timber rotations, increased forest management intensity, and avoided deforestation.In simulating the effects at different carbon prices, we ran FASOMGHG one strategy at a time and later allowed the whole portfolio of strategies.
The runs were done under several hypothetical carbon price scenarios which contained rising prices.Specifically, the prices were initially set at $5, $10, $30, and $50 per metric ton CO 2 equivalent but then escalated at 5% per year.Additionally, we ran a baseline scenario with no carbon price.The considered shifts in land use were those under the carbon price scenario relative to those under the baseline scenario.
The version of FASOMGHG used has the conterminous United States divided into 11 regions.The Missouri River Basin consists of parts of the FASOMGHG regions of the Corn Belt, Northern Plains, and Rocky Mountains, and based on land use shares, we used 10.90% of the Corn Belt results, 45.72% of the Northern Plains, and 43.38%, of the Rocky Mountains.
FASOMGHG generated many results, but here, we will only focus on those with implications for water.The FASOMGHG results for cropped land, grassland, and forested land plus the regression evaluations for the 50% quantile are reported in Table 6.All the values use simulated land use changes during the five-year period beginning with 2025.The first carbon price scenario-$5 at a 5% increase-shows that in all of the single strategy mitigation runs, excepting the bioenergy case, there is a small increase in cropped land with an accompanying decrease in both grassland and forest land.When carbon prices are only paid for crop fertilization reduction, the results show cropped land use increases by 0.83%, while grassland and forests decrease by 0.78% and 0.16%, causing an increase in water quantity and a decrease in water quality.The opposite effects occur for bioenergy, forest management, and the case where all strategies are allowed with water quantity decreasing but quality increasing.Similar results are found under the $10 at a 5% increase carbon price scenario.There, all the single strategy mitigation cases expand cropping but reduce grassland and forests.

Effects under the Higher Carbon Price Scenarios
In the higher carbon price scenarios (those starting with $30 and $50 prices), all the single strategy mitigation results and the simultaneous case decrease cropping and forests but increase grassland.In these scenarios, water quantity is decreased while water quality increases across the board.When forest management is the only available strategy, this causes the largest water quantity decrease and water quality increase across the individual strategies.The all strategies case causes yet a larger quantity reduction (1.05 mm per acre in the $30 case) and the largest increase in water quality.This suggests that stronger AF mitigation efforts will result in reduced water quantity but increased water quality.We also see diminishing effects on quantity and quality as prices go higher.

Conclusions
This paper reports water quality and quantity implications of agricultural and forestry climate change mitigation strategies.To do this, we did an empirical study in the Missouri River Basin investigating the effects of differing levels of mitigation incentives and involving four phases.
First, we ran SWAT for 21 years, observing water quantity and quality effects on a monthly basis for the 13,437 hydrologic unit watershed levels.Then, we took those results and aggregated them to counties, in turn, doing a quantile regression estimation over panel data explaining the water quantity and quality effects of altered land use shares and climate.The result shows that an increase in grassland significantly decreases water quantity but increases quality, while changes in cropped land share have the opposite effect.We also found an increase in forested land share had mixed effects.
The third empirical phase examined the land use effects of an agricultural and forestry mitigation strategy adopted under alternative carbon prices.The fourth phase integrated the quantile regressions and the land use effects to project water quantity and quality effects of various carbon prices and mitigation strategies.The consequent results showed crop fertilization, afforestation, and livestock strategies generally led to increased water quantity and decreased quality under lower carbon price scenarios but led to the opposite result under higher carbon price scenarios.The results also showed that water quality is degraded under most mitigation alternatives except for the production of bioenergy feedstocks and the alteration of forest management.In these cases, these effects occurred when carbon price was low but was reversed with higher carbon prices.
Collectively, the results imply that AF mitigation has varied effects on water quantity and water quality, depending on the carbon price or the extent of mitigation activity.Generally, water quantity is slightly increased under lower carbon prices but then is reduced when carbon prices go higher.Water quality results reveal mixed effects at lower prices but an improvement under higher carbon prices.We also see the magnitudes of improvement shrink as the carbon price goes higher.

Figure 2 .
Figure 2. The proportion of each land use type in the Missouri River basin.

Figure 2 .
Figure 2. The proportion of each land use type in the Missouri River basin.

Figure 3 .
Figure 3.The total effects on water quantity when switching land from cropping to grassland.

Figure 4 .
Figure 4.The total effects on water quantity when switching land from cropping to forest.

Figure 3 .
Figure 3.The total effects on water quantity when switching land from cropping to grassland.

Figure 3 .
Figure 3.The total effects on water quantity when switching land from cropping to grassland.

Figure 4 .
Figure 4.The total effects on water quantity when switching land from cropping to forest.

Figure 4 .
Figure 4.The total effects on water quantity when switching land from cropping to forest.

Figure 5 .
Figure 5.The total effects on water quantity when switching land from grass to forest.

Table 1 .
Summary statistics over the dataset generated using SWAT.

Table 2 .
Quantile statistics over the water quantity and water quality dataset generated by SWAT.: Calculated over MRB simulation results from the SWAT model.Note: 1.The water quality index is described above and summarizes the results of indicators of total nitrogen and total phosphorus in surface runoff, and is scaled from 10 (the worst water quality) to 100 (the best water quality). Source

Table 3 .
Summary Statistics in the Missouri River basin based on water quality index.

Table 4 .
Estimation results of the effects on water quantity (per county per month).

Table 5 .
Estimation results of the effects on water quality index (per county per month).

Table 6 .
Effects on water when carbon prices are applied to alternative sets of mitigation strategies.