Wind-Driven Sediment Resuspension in the World's Fourth Largest Lake Contributes Substantial Phosphorus Load to the 11th Largest Lake.

Capturing precipitation-based episodes is a longstanding issue for estimating tributary loads; however, wind-driven resuspension in Lake Huron creates similar uncertainties in its estimated load to Lake Erie. Recent suggestions that the phosphorus load from Lake Huron is underestimated because sampling frequencies miss contributions from resuspension events are speculative because they did not include direct load measurements, address all resuspension regions, or assess the potential bioavailability of the load. We address these shortcomings by evaluating Lake Huron's nearshore regions, characterizing the biological availability of the load, and providing direct comparisons of load estimates with and without the resuspended load. We show that total phosphorus concentrations in Lake Huron and the St. Clair River are higher during resuspension events and that bioavailability of that material is comparable to that reported elsewhere. New load estimates, based on continuous turbidity measurements converted to phosphorus through P-turbidity relationships, were almost 90% higher than traditional load estimates, providing empirical evidence for the significantly underestimated previous load. This confirmation is important because if the Lake Huron load is not decreased, reductions from other sources would be needed to meet the overall reduction targets set by the binational Great Lakes Water Quality Agreement.

Advancing estuarine ecological forecasts: seasonal hypoxia in Chesapeake Bay.

Ecological forecasts are quantitative tools that can guide ecosystem management. The coemergence of extensive environmental monitoring and quantitative frameworks allows for widespread development and continued improvement of ecological forecasting systems. We use a relatively simple estuarine hypoxia model to demonstrate advances in addressing some of the most critical challenges and opportunities of contemporary ecological forecasting, including predictive accuracy, uncertainty characterization, and management relevance. We explore the impacts of different combinations of forecast metrics, drivers, and driver time windows on predictive performance. We also incorporate multiple sets of state-variable observations from different sources and separately quantify model prediction error and measurement uncertainty through a flexible Bayesian hierarchical framework. Results illustrate the benefits of (1) adopting forecast metrics and drivers that strike an optimal balance between predictability and relevance to management, (2) incorporating multiple data sources in the calibration data set to separate and propagate different sources of uncertainty, and (3) using the model in scenario mode to probabilistically evaluate the effects of alternative management decisions on future ecosystem state. In the Chesapeake Bay, the subject of this case study, we find that average summer or total annual hypoxia metrics are more predictable than monthly metrics and that measurement error represents an important source of uncertainty. Application of the model in scenario mode suggests that absent watershed management actions over the past decades, long-term average hypoxia would have increased by 7% compared to 1985. Conversely, the model projects that if management goals currently in place to restore the Bay are met, long-term average hypoxia would eventually decrease by 32% with respect to the mid-1980s.

Evaluating management options to reduce Lake Erie algal blooms using an ensemble of watershed models.

Reducing harmful algal blooms in Lake Erie, situated between the United States and Canada, requires implementing best management practices to decrease nutrient loading from upstream sources. Bi-national water quality targets have been set for total and dissolved phosphorus loads, with the ultimate goal of reaching these targets in 9-out-of-10 years. Row crop agriculture dominates the land use in the Western Lake Erie Basin thus requiring efforts to mitigate nutrient loads from agricultural systems. To determine the types and extent of agricultural management practices needed to reach the water quality goals, we used five independently developed Soil and Water Assessment Tool models to evaluate the effects of 18 management scenarios over a 10-year period on nutrient export. Guidance from a stakeholder group was provided throughout the project, and resulted in improved data, development of realistic scenarios, and expanded outreach. Subsurface placement of phosphorus fertilizers, cover crops, riparian buffers, and wetlands were among the most effective management options. But, only in one realistic scenario did a majority (3/5) of the models predict that the total phosphorus loading target would be met in 9-out-of-10 years. Further, the dissolved phosphorus loading target was predicted to meet the 9-out-of-10-year goal by only one model and only in three scenarios. In all scenarios evaluated, the 9-out-of-10-year goal was not met based on the average of model predictions. Ensemble modeling revealed general agreement about the effects of several practices although some scenarios resulted in a wide range of uncertainty. Overall, our results demonstrate that there are multiple pathways to approach the established water quality goals, but greater adoption rates of practices than those tested here will likely be needed to attain the management targets.

Uncertainty in critical source area predictions from watershed-scale hydrologic models.

Watershed-scale hydrologic models are frequently used to inform conservation and restoration efforts by identifying critical source areas (CSAs; alternatively 'hotspots'), defined as areas that export relatively greater quantities of nutrients and sediment. The CSAs can then be prioritized or 'targeted' for conservation and restoration to ensure efficient use of limited resources. However, CSA simulations from watershed-scale hydrologic models may be uncertain and it is critical that the extent and implications of this uncertainty be conveyed to stakeholders and decision makers. We used an ensemble of four independently developed Soil and Water Assessment Tool (SWAT) models and a SPAtially Referenced Regression On Watershed attributes (SPARROW) model to simulate CSA locations for flow, phosphorus, nitrogen, and sediment within the ~17,000-km2 Maumee River watershed at the HUC-12 scale. We then assessed uncertainty in CSA simulations determined as the variation in CSA locations across the models. Our application of an ensemble of models - differing with respect to inputs, structure, and parameterization - facilitated an improved accounting of CSA prediction uncertainty. We found that the models agreed on the location of a subset of CSAs, and that these locations may be targeted with relative confidence. However, models more often disagreed on CSA locations. On average, only 16%-46% of HUC-12 subwatersheds simulated as a CSA by one model were also simulated as a CSA by a different model. Our work shows that simulated CSA locations are highly uncertain and may vary substantially across models. Hence, while models may be useful in informing conservation and restoration planning, their application to identify CSA locations would benefit from comprehensive uncertainty analyses to avoid inefficient use of limited resources.

Lake Huron's Phosphorus Contributions to the St. Clair-Detroit River Great Lakes Connecting Channel.

The United States and Canada called for a 40% load reduction of total phosphorus from 2008 levels entering the western and central basins of Lake Erie to achieve a 6000 MTA target and help reduce its central basin hypoxia. The Detroit River is a significant source of total phosphorus to Lake Erie; it in turn has been reported to receive up to 58% of its load from Lake Huron when accounting for resuspended sediment loads previously unmonitored at the lake outlet. Key open questions are where does this additional load originate, what drives its variability, and how often does it occur. We used a hydrodynamic model, satellite images of resuspension events and ice cover, wave hindcasts, and continuous turbidity measurements at the outlet of Lake Huron to determine where in Lake Huron the undetected load originates and what drives its variability. We show that the additional sediment load, and likely phosphorus, is from wave-induced Lake Huron sediment resuspension, primarily within 30 km of the southeastern shore. When the flow is from southwest or down the center of the lake, the resuspended sediment is not detected at Canada's sampling station at the head of the St. Clair River.

Ensemble modeling informs hypoxia management in the northern Gulf of Mexico.

A large region of low-dissolved-oxygen bottom waters (hypoxia) forms nearly every summer in the northern Gulf of Mexico because of nutrient inputs from the Mississippi River Basin and water column stratification. Policymakers developed goals to reduce the area of hypoxic extent because of its ecological, economic, and commercial fisheries impacts. However, the goals remain elusive after 30 y of research and monitoring and 15 y of goal-setting and assessment because there has been little change in river nitrogen concentrations. An intergovernmental Task Force recently extended to 2035 the deadline for achieving the goal of a 5,000-km2 5-y average hypoxic zone and set an interim load target of a 20% reduction of the spring nitrogen loading from the Mississippi River by 2025 as part of their adaptive management process. The Task Force has asked modelers to reassess the loading reduction required to achieve the 2035 goal and to determine the effect of the 20% interim load reduction. Here, we address both questions using a probabilistic ensemble of four substantially different hypoxia models. Our results indicate that, under typical weather conditions, a 59% reduction in Mississippi River nitrogen load is required to reduce hypoxic area to 5,000 km2 The interim goal of a 20% load reduction is expected to produce an 18% reduction in hypoxic area over the long term. However, due to substantial interannual variability, a 25% load reduction is required before there is 95% certainty of observing any hypoxic area reduction between consecutive 5-y assessment periods.

Climate Change and Nutrient Loading in the Western Lake Erie Basin: Warming Can Counteract a Wetter Future.

In the past 20 years, Lake Erie has experienced a resurgence of harmful algal blooms and hypoxia driven by increased nutrient loading from its agriculturally dominated watersheds. The increase in phosphorus loading, specifically the dissolved reactive portion, has been attributed to a combination of changing climate and agricultural management. While many management practices and strategies have been identified to reduce phosphorus loads, the impacts of future climate remain uncertain. This is particularly the case for the Great Lakes region because many global climate models do not accurately represent the land-lake interactions that govern regional climate. For this study, we used midcentury (2046-2065) climate projections from one global model and four regional dynamically downscaled models as drivers for the Soil and Water Assessment Tool configured for the Maumee River watershed, the source of almost 50% of Lake Erie's Western Basin phosphorus load. Our findings suggest that future warming may lead to less nutrient runoff due to increased evapotranspiration and decreased snowfall, despite projected moderate increases in intensity and overall amount of precipitation. Results highlight the benefits of considering multiple environmental drivers in determining the fate of nutrients in the environment and demonstrate a need to improve approaches for climate change assessment using watershed models.

Evaluating the Impact of Legacy P and Agricultural Conservation Practices on Nutrient Loads from the Maumee River Watershed.

The recent resurgence of hypoxia and harmful algal blooms in Lake Erie, driven substantially by phosphorus loads from agriculture, have led the United States and Canada to begin developing plans to meet new phosphorus load targets. To provide insight into which agricultural management options could help reach these targets, we tested alternative agricultural-land-use and land-management scenarios on phosphorus loads to Lake Erie. These scenarios highlight certain constraints on phosphorus load reductions from changes in the Maumee River Watershed (MRW), which contributes roughly half of the phosphorus load to the lake's western basin. We evaluate the effects on phosphorus loads under nutrient management strategies, reduction of fertilizer applications, employing vegetative buffers, and implementing widespread cover crops and alternative cropping changes. Results indicate that even if fertilizer application ceased, it may take years to see desired decreases in phosphorus loads, especially if we experience greater spring precipitation or snowmelt. Scenarios also indicate that widespread conversions to perennial crops that may be used for biofuel production are capable of substantially reducing phosphorus loads. This work demonstrates that a combination of legacy phosphorus, land management, land use, and climate should all be considered when seeking phosphorus-loading solutions.

Engaging Stakeholders To Define Feasible and Desirable Agricultural Conservation in Western Lake Erie Watersheds.

Widespread adoption of agricultural conservation measures in Lake Erie's Maumee River watershed may be required to reduce phosphorus loading that drives harmful algal blooms and hypoxia. We engaged agricultural and conservation stakeholders through a survey and workshops to determine which conservation practices to evaluate. We investigated feasible and desirable conservation practices using the Soil and Water Assessment Tool calibrated for streamflow, sediment, and nutrient loading near the Maumee River outlet. We found subsurface placement of phosphorus applications to be the individual practice most influential on March-July dissolved reactive phosphorus (DRP) loading from row croplands. Perennial cover crops and vegetated filter strips were most effective for reducing seasonal total phosphorus (TP) loading. We found that practices effective for reducing TP and DRP load were not always mutually beneficial, culminating in trade-offs among multiple Lake Erie phosphorus management goals. Adoption of practices at levels considered feasible to stakeholders led to nearly reaching TP targets for western Lake Erie on average years; however, adoption of practices at a rate that goes beyond what is currently considered feasible will likely be required to reach the DRP target.

Record-setting algal bloom in Lake Erie caused by agricultural and meteorological trends consistent with expected future conditions.

In 2011, Lake Erie experienced the largest harmful algal bloom in its recorded history, with a peak intensity over three times greater than any previously observed bloom. Here we show that long-term trends in agricultural practices are consistent with increasing phosphorus loading to the western basin of the lake, and that these trends, coupled with meteorological conditions in spring 2011, produced record-breaking nutrient loads. An extended period of weak lake circulation then led to abnormally long residence times that incubated the bloom, and warm and quiescent conditions after bloom onset allowed algae to remain near the top of the water column and prevented flushing of nutrients from the system. We further find that all of these factors are consistent with expected future conditions. If a scientifically guided management plan to mitigate these impacts is not implemented, we can therefore expect this bloom to be a harbinger of future blooms in Lake Erie.

Assessing biophysical controls on Gulf of Mexico hypoxia through probabilistic modeling.

A mechanistic model was developed to predict midsummer bottom-water dissolved oxygen (BWDO) concentration and hypoxic area on the Louisiana shelf of the northern Gulf of Mexico, USA (1985-2011). Because of its parsimonious formulation, the model possesses many of the benefits of simpler, more empirical models, in that it is computationally efficient and can rigorously account for uncertainty through Bayesian inference. At the same time, the model incorporates important biophysical processes such that its parameterization can be informed by field-measured biological and physical rates. The model is used to explore how freshwater flow, nutrient load, benthic oxygen demand, and wind velocity affect hypoxia on the western and eastern sections of the shelf, delineated by the Atchafalaya River outfall. The model explains over 70% of the variability in BWDO on both shelf sections, and outperforms linear regression models developed from the same input variables. Model results suggest that physical factors (i.e., wind and flow) control a larger portion of the year-to-year variability in hypoxia than previously thought, especially on the western shelf, though seasonal nutrient loads remain an important driver of hypoxia, as well. Unlike several previous Gulf hypoxia modeling studies, results do not indicate a temporal shift in the system's propensity for hypoxia formation (i.e., no regime change). Results do indicate that benthic oxygen demand is a substantial BWDO sink, and a better understanding of the long-term dynamics of this sink is required to better predict how the size of the hypoxic zone will respond to proposed reductions in nutrient loading.

A scenario and forecast model for Gulf of Mexico hypoxic area and volume.

For almost three decades, the relative size of the hypoxic region on the Louisiana-Texas continental shelf has drawn scientific and policy attention. During that time, both simple and complex models have been used to explore hypoxia dynamics and to provide management guidance relating the size of the hypoxic zone to key drivers. Throughout much of that development, analyses had to accommodate an apparent change in hypoxic sensitivity to loads and often cull observations due to anomalous meteorological conditions. Here, we describe an adaptation of our earlier, simple biophysical model, calibrated to revised hypoxic area estimates and new hypoxic volume estimates through Bayesian estimation. This application eliminates the need to cull observations and provides revised hypoxic extent estimates with uncertainties corresponding to different nutrient loading reduction scenarios. We compare guidance from this model application, suggesting an approximately 62% nutrient loading reduction is required to reduce Gulf hypoxia to the Action Plan goal of 5000 km(2), to that of previous applications. In addition, we describe for the first time, the corresponding response of hypoxic volume. We also analyze model results to test for increasing system sensitivity to hypoxia formation, but find no strong evidence of such change.

Spatial and temporal trends in Lake Erie hypoxia, 1987-2007.

Hypoxic conditions, defined as dissolved oxygen (DO) concentrations below 2 mg/L, are a regular summertime occurrence in Lake Erie, but the spatial extent has been poorly understood due to sparse sampling. We use geostatistical kriging and conditional realizations to provide quantitative estimates of the extent of hypoxia in the central basin of Lake Erie for August and September of 1987 to 2007, along with their associated uncertainties. The applied geostatistical approach combines the limited in situ DO measurements with auxiliary data selected using the Bayesian Information Criterion. Bathymetry and longitude are found to be highly significant in explaining the spatial distribution of DO, while satellite observations of sea surface temperature and satellite chlorophyll are not. The hypoxic extent was generally lowest in the mid-1990s, with the late 1980s (1987, 1988) and the 2000s (2003, 2005) experiencing the largest hypoxic zones. A simple exponential relationship based on the squared average measured bottom DO explains 97% of the estimated variability in the hypoxic extent. The change in the observed maximum extent between August and September is found to be sensitive to the corresponding variability in the hypolimnion thickness.

Retrospective analysis of midsummer hypoxic area and volume in the northern Gulf of Mexico, 1985-2011.

Robust estimates of hypoxic extent (both area and volume) are important for assessing the impacts of low dissolved oxygen on aquatic ecosystems at large spatial scales. Such estimates are also important for calibrating models linking hypoxia to causal factors, such as nutrient loading and stratification, and for informing management decisions. In this study, we develop a rigorous geostatistical modeling framework to estimate the hypoxic extent in the northern Gulf of Mexico from data collected during midsummer, quasi-synoptic monitoring cruises (1985-2011). Instead of a traditional interpolation-based approach, we use a simulation-based approach that yields more robust extent estimates and quantified uncertainty. The modeling framework also makes use of covariate information (i.e., trend variables such as depth and spatial position), to reduce estimation uncertainty. Furthermore, adjustments are made to account for observational bias resulting from the use of different sampling instruments in different years. Our results suggest an increasing trend in hypoxic layer thickness (p = 0.05) from 1985 to 2011, but less than significant increases in volume (p = 0.12) and area (p = 0.42). The uncertainties in the extent estimates vary with sampling network coverage and instrument type, and generally decrease over the study period.

Advancing freshwater ecological forecasts: Harmful algal blooms in Lake Erie.

Ecological models help provide forecasts of ecosystem responses to natural and anthropogenic stresses. However, their ability to create reliable predictions requires forecasts with track records sufficiently long to build confidence, skill assessments, and treating uncertainty quantitatively. We use Lake Erie harmful algal blooms as a case study to help formalize ecological forecasting. Key challenges for models include uncertainty in the deterministic structure of the load-bloom relationship and the need to assess alternative drivers (e.g., biologically available phosphorus load, spring load, longer term cumulative load) with a larger dataset. We enhanced a Bayesian model considering new information and an expanded data set, test it through cross validation and blind forecasts, quantify and discuss its uncertainties, and apply it for assessing historical and future scenarios. Allowing a segmented relationship between bloom size and spring load indicates that loading above 0.15 Gg/month will have a substantially higher marginal impact on bloom size. The new model explains 84 % of interannual variability (9.09 Gg RMSE) when calibrated to the 19-year data set and 66 % of variability in cross validation (12.58 Gg RMSE). Blind forecasts explain 84 % of HAB variability between 2014 and 2020, which is substantially better than the actual forecast track record (R2 = 0.32) over this same period. Because of internal phosphorus recycling, represented by the long-term cumulative load, it could take over a decade for HABs to fully respond to loading reductions, depending on the pace of those reductions. Thus, the desired speed and endpoint of the lake's recovery should be considered when updating and adaptively managing load reduction targets. Results are discussed in the context of ecological forecasting best pactices: incorporate new knowledge and data in model construction; account for multiple sources of uncertainty; evaluate predictive skill through validation and hindcasting; and answer management questions related to both short-term forecasts and long-term scenarios.

Evaluating causes of trends in long-term dissolved reactive phosphorus loads to Lake Erie.

Renewed harmful algal blooms and hypoxia in Lake Erie have drawn significant attention to phosphorus loads, particularly increased dissolved reactive phosphorus (DRP) from highly agricultural watersheds. We use the Soil and Water Assessment Tool (SWAT) to model DRP in the agriculture-dominated Sandusky watershed for 1970-2010 to explore potential reasons for the recent increased DRP load from Lake Erie watersheds. We demonstrate that recent increased storm events, interacting with changes in fertilizer application timing and rate, as well as management practices that increase soil stratification and phosphorus accumulation at the soil surface, appear to drive the increasing DRP trend after the mid-1990s. This study is the first long-term, detailed analysis of DRP load estimation using SWAT.

Quantifying the impacts of stratification and nutrient loading on hypoxia in the northern Gulf of Mexico.

Stratification and nutrient loading are two primary factors leading to hypoxia in coastal systems. However, where these factors are temporally correlated, it can be difficult to isolate and quantify their individual impacts. This study provides a novel solution to this problem by determining the effect of stratification based on its spatial relationship with bottom-water dissolved oxygen (BWDO) concentration using a geostatistical regression. Ten years (1998-2007) of midsummer Gulf of Mexico BWDO measurements are modeled using stratification metrics along with trends based on spatial coordinates and bathymetry, which together explain 27-61% of the spatial variability in BWDO for individual years. Because stratification effects explain only a portion of the year-to-year variability in mean BWDO; the remaining variability is explained by other factors, with May nitrate plus nitrite river concentration the most important. Overall, 82% of the year-to-year variability in mean BWDO is explained. The results suggest that while both stratification and nutrients play important roles in determining the annual extent of midsummer hypoxia, reducing nutrient inputs alone will substantially reduce the average extent.

Incidental oligotrophication of North American Great Lakes.

Phytoplankton production is an important factor in determining both ecosystem stability and the provision of ecosystem goods and services. The expansive and economically important North American Great Lakes are subjected to multiple stressors and understanding their responses to those stresses is important for understanding system-wide ecological controls. Here we show gradual increases in spring silica concentration (an indicator of decreasing growth of the dominant diatoms) in all basins of Lakes Michigan and Huron (USA and Canadian waters) between 1983 and 2008. These changes indicate the lakes have undergone gradual oligotrophication coincident with and anticipated by nutrient management implementation. Slow declines in seasonal drawdown of silica (proxy for seasonal phytoplankton production) also occurred, until recent years, when lake-wide responses were punctuated by abrupt decreases, putting them in the range of oligotrophic Lake Superior. The timing of these dramatic production drops is coincident with expansion of populations of invasive dreissenid mussels, particularly quagga mussels, in each basin. The combined effect of nutrient mitigation and invasive species expansion demonstrates the challenges facing large-scale ecosystems and suggest the need for new management regimes for large ecosystems.

Simulating internal watershed processes using multiple SWAT models.

The need for effective water quality models to help guide management and policy, and extend monitoring information, is at the forefront of recent discussions related to watershed management. These models are often calibrated and validated at the basin outlet, which ensures that models are capable of evaluating basin scale hydrology and water quality. However, there is a need to understand where these models succeed or fail with respect to internal process representation, as these watershed-scale models are used to inform management practices and mitigation strategies upstream. We evaluated an ensemble of models-each calibrated to in-stream observations at the basin outlet-against discharge and nutrient observations at the farm field scale to determine the extent to which these models capture field-scale dynamics. While all models performed well at the watershed outlet, upstream performance varied. Models tended to over-predict discharge through surface runoff and subsurface drainage, while under-predicting phosphorus loading through subsurface drainage and nitrogen loading through surface runoff. Our study suggests that while models may be applied to predict impacts of management at the basin scale, care should be taken in applying the models to evaluate field-scale management and processes in the absence of data that can be incorporated at that scale, even with the use of multiple models.

Bias correction of climate model outputs influences watershed model nutrient load predictions.

Waterbodies around the world experience problems associated with elevated phosphorus (P) and nitrogen (N) loads. While vital for ecosystem functioning, when present in excess amounts these nutrients can impair water quality and create symptoms of eutrophication, including harmful algal blooms. Under a changing climate, nutrient loads are likely to change. While climate models can serve as inputs to watershed models, the climate models often do not adequately represent the distribution of observed data, generating uncertainties that can be addressed to some degree with bias correction. However, the impacts of bias correction on nutrient models are not well understood. This study compares 4 univariate and 3 multivariate bias correction methods, which correct precipitation and temperature variables from 4 climate models in the historical (1980-1999) and mid-century future (2046-2065) time periods. These variables served as inputs to a calibrated Soil and Water Assessment Tool (SWAT) model of Lake Erie's Maumee River watershed. We compared the performance of SWAT outputs driven with climate model outputs that were bias-corrected (BC) and not bias-corrected (no-BC) for dissolved reactive P, total P, and total N. Results based on graphical comparisons and goodness of fit metrics showed that the choice of BC method impacts both the direction of change and magnitude of nutrient loads and hydrological processes. While the Delta method performed best, it should be used with caution since it considers historical variable relationships as the basis for predictions, which may not hold true under future climate. Quantile Delta Mapping (QDM) and Multivariate Bias Correction N-dimensional probability density function transform (MBCn) BC methods also performed well and work well for non-stationary climate scenarios. Furthermore, results suggest that February-July cumulative load in the Maumee basin is likely to decrease in the mid-century as runoff and snowfall decrease, and evapotranspiration increases with warming temperatures.