Characterization of legacy organic carbon in a culturally eutrophic lake : the role of the historic carbon deposition in the time course and extent of lake recovery

The time course of lake recovery after a reduction in external loading of nutrients is often controlled by conditions in the sediment. Remediation of eutrophication is hindered by the presence of legacy organic carbon deposits, that exert a demand on the terminal electron acceptors of the lake and contribute to problems such as internal nutrient recycling, absence of sediment macrofauna, and flux of toxic metal species into the water column. Being able to quantify the timing of a lake’s response requires determination of the magnitude and lability, i.e., the susceptibility to biodegradation, of the organic carbon within the legacy deposit. This characterization is problematic for organic carbon in sediments because of the presence of different fractions of carbon, which vary from highly labile to refractory. The lability of carbon under varied conditions was tested with a bioassay approach. It was found that the majority of the organic material found in the sediments is conditionally-labile, where mineralization potential is dependent on prevailing conditions. High labilities were noted under oxygenated conditions and a favorable temperature of 30 °C. Lability decreased when oxygen was removed, and was further reduced when the temperature was dropped to the hypolimnetic average of 8° C . These results indicate that reversible preservation mechanisms exist in the sediment, and are able to protect otherwise labile material from being mineralized under in situ conditions. The concept of an active sediment layer, a region in the sediments in which diagenetic reactions occur (with nothing occurring below it), was examined through three lines of evidence. Initially, porewater profiles of oxygen, nitrate, sulfate/total sulfide, ETSA (Electron Transport System Activity- the activity of oxygen, nitrate, iron/manganese, and sulfate), and methane were considered. It was found through examination of the porewater profiles that the edge of diagenesis occurred around 15-20 cm. Secondly, historical and contemporary TOC profiles were compared to find the point at which the profiles were coincident, indicating the depth at which no change has occurred over the (13 year) interval between core collections. This analysis suggested that no diagenesis has occurred in Onondaga Lake sediment below a depth of 15 cm. Finally, the time to 99% mineralization, the t99, was viewed by using a literature estimate of the kinetic rate constant for diagenesis. A t99 of 34 years, or approximately 30 cm of sediment depth, resulted for the slowly decaying carbon fraction. Based on these three lines of evidence , an active sediment layer of 15-20 cm is proposed for Onondaga Lake, corresponding to a time since deposition of 15-20 years. While a large legacy deposit of conditionally-labile organic material remains in the sediments of Onondaga Lake, it becomes clear that preservation, mechanisms that act to shield labile organic carbon from being degraded, protects this material from being mineralized and exerting a demand on the terminal electron acceptors of the lake. This has major implications for management of the lake, as it defines the time course of lake recovery following a reduction in nutrient loading.

Retention of canopy trees as biological legacies for balancing woody biomass production and biodiversity in managed aspen forests of the Great Lakes Region

Green-tree retention under the conceptual framework of ecological forestry has the potential to provide both biomass feedstock for industry and maintain quality wildlife habitat. I examined the effects of retained canopy trees as biological legacies (“legacy trees”) in aspen (Populus spp.) forests on above-ground live woody biomass, understory plant floristic quality, and bird diversity. Additionally, I evaluated habitat quality for a high conservation priority species, the Golden-winged Warbler (Vermivora chrysoptera). I selected 27 aspen-dominated forest stands in northern Wisconsin with nine stands in each of three legacy tree retention treatments (conifer retention, hardwood retention, and clearcuts or no retention) across a chronosequence (4-36 years post-harvest). Conifer retention stands had greater legacy tree and all tree species biomass but lower regenerating tree biomass than clearcuts. Coniferous but not hardwood legacy trees appeared to suppress regenerating tree biomass. I evaluated the floristic quality of the understory plant assemblage by estimating the mean coefficient of conservatism (C). Mean C was lower in young stands than in middle-age or old stands; there was a marginally significant (p=0.058) interaction effect between legacy tree retention treatment and stand age. Late-seral plant species were positively associated with stand age and legacy tree diameter or age revealing an important relationship between legacy tree retention and stand development. Bird species richness was greatest in stands with hardwood retention particularly early in stand development. Six conservation priority bird species were indicators of legacy tree retention or clearcuts. Retention of legacy trees in aspen stands provided higher quality nest habitat for the Golden-winged Warbler than clearcuts based on high pairing success and nesting activity. Retention of hardwoods, particularly northern red oak (Quercus rubra), yielded the most consistent positive effects in this study with the highest bird species richness and the highest quality habitat for the Golden-winged Warbler. This treatment maintained stand biomass comparable to clearcuts and did not suppress regenerating tree biomass. In conclusion, legacy tree retention can enhance even-aged management techniques to produce a win-win scenario for the conservation of declining bird species and late-seral understory plants and for production of woody biomass feedstock from naturally regenerating aspen forests.

Public preferences of the Great Lakes Environment : a Lake Michigan pilot study

The Great Lakes watershed is home to over 40 million people, and the health of the Great Lakes ecosystem is vital to the overall economic, societal, and environmental health of the U.S. and Canada. However, environmental issues related to them are sometimes overlooked. Policymakers and the public face the challenges of balancing economic benefits with the need to conserve and/or replenish regional natural resources to ensure long term prosperity. From the literature review, nine critical stressors of ecological services were delineated, which include pollution and contamination, agricultural erosion, non-native species, degraded recreational resources, loss of wetlands habitat, climate change, risk of clean water shortage, vanishing sand dunes, and population overcrowding; this list was validated through a series of stakeholder discussions and focus groups in Grand Rapids. Focus groups were conducted in Grand Rapids to examine the awareness of, concern with, and willingness to expend resources on these stressors. Stressors that the respondents have direct contact with tend to be the most important. The focus group results show that concern related to pollution and contamination is much higher than for any of the other stressors. Low responses to climate change result in recommendations for outreach programs.

Great Lakes tributary phosphorus bioavailability

Information on phosphorus bioavailability can provide water quality managers with the support required to target point source and watershed loads contributing most significantly to water quality conditions. This study presents results from a limited sampling program focusing on the five largest sources of total phosphorus to the U.S. waters of the Great Lakes. The work provides validation of the utility of a bioavailability-based approach, confirming that the method is robust and repeatable. Chemical surrogates for bioavailability were shown to hold promise, however further research is needed to address site-to-site and seasonal variability before a universal relationship can be accepted. Recent changes in the relative contribution of P constituents to the total phosphorus analyte and differences in their bioavailability suggest that loading estimates of bioavailable P will need to address all three components (SRP, DOP and PP). A bioavailability approach, taking advantage of chemical surrogate methodologies is recommended as a means of guiding P management in the Great Lakes.

Critical evaluation of diameter increment modelling in the Great Lakes region

Simulations of forest stand dynamics in a modelling framework including Forest Vegetation Simulator (FVS) are diameter driven, thus the diameter or basal area increment model needs a special attention. This dissertation critically evaluates diameter or basal area increment models and modelling approaches in the context of the Great Lakes region of the United States and Canada. A set of related studies are presented that critically evaluate the sub-model for change in individual tree basal diameter used in the Forest Vegetation Simulator (FVS), a dominant forestry model in the Great Lakes region. Various historical implementations of the STEMS (Stand and Tree Evaluation and Modeling System) family of diameter increment models, including the current public release of the Lake States variant of FVS (LS-FVS), were tested for the 30 most common tree species using data from the Michigan Forest Inventory and Analysis (FIA) program. The results showed that current public release of the LS-FVS diameter increment model over-predicts 10-year diameter increment by 17% on average. Also the study affirms that a simple adjustment factor as a function of a single predictor, dbh (diameter at breast height) used in the past versions, provides an inadequate correction of model prediction bias. In order to re-engineer the basal diameter increment model, the historical, conceptual and philosophical differences among the individual tree increment model families and their modelling approaches were analyzed and discussed. Two underlying conceptual approaches toward diameter or basal area increment modelling have been often used: the potential-modifier (POTMOD) and composite (COMP) approaches, which are exemplified by the STEMS/TWIGS and Prognosis models, respectively. It is argued that both approaches essentially use a similar base function and neither is conceptually different from a biological perspective, even though they look different in their model forms. No matter what modelling approach is used, the base function is the foundation of an increment model. Two base functions – gamma and Box-Lucas – were identified as candidate base functions for forestry applications. The results of a comparative analysis of empirical fits showed that quality of fit is essentially similar, and both are sufficiently detailed and flexible for forestry applications. The choice of either base function in order to model diameter or basal area increment is dependent upon personal preference; however, the gamma base function may be preferred over the Box-Lucas, as it fits the periodic increment data in both a linear and nonlinear composite model form. Finally, the utility of site index as a predictor variable has been criticized, as it has been widely used in models for complex, mixed species forest stands though not well suited for this purpose. An alternative to site index in an increment model was explored, using site index and a combination of climate variables and Forest Ecosystem Classification (FEC) ecosites and data from the Province of Ontario, Canada. The results showed that a combination of climate and FEC ecosites variables can replace site index in the diameter increment model.

Applications of LiDAR remote sensing of forest structure in the Upper Great Lakes region, USA

We used active remote sensing technology to characterize forest structure in a northern temperate forest on a landscape- and local-level in the Upper Peninsula of Michigan. Specifically, we used a form of active remote sensing called light detection and ranging (e.g., LiDAR) to aid in the depiction of current forest structural stages and total canopy gap area estimation. On a landscape-level, LiDAR data are shown not only to be a useful tool in characterizing forest structure, in both coniferous and deciduous forest cover types, but also as an effective basis for data-driven surrogates for classification of forest structure. On a local-level, LiDAR data are shown to be a benchmark reference point to evaluate field-based canopy gap area estimations, due to the highly accurate nature of such remotely sensed data. The application of LiDAR remote sensed data can help facilitate current and future sustainable forest management.

DEVELOPMENT OF A CONCEPTUAL MODEL AND PARAMETERIZATION OF A MATHEMATICAL MODEL FOR SEDIMENT PHOSPHORUS DIAGENESIS, SORPTION AND RELEASE FROM ONONDAGA LAKE SYRACUSE, NEW YORK

Eutrophication is a persistent problem in many fresh water lakes. Delay in lake recovery following reductions in external loading of phosphorus, the limiting nutrient in fresh water ecosystems, is often observed. Models have been created to assist with lake remediation efforts, however, the application of management tools to sediment diagenesis is often neglected due to conceptual and mathematical complexity. SED2K (Chapra et al. 2012) is proposed as a "middle way", offering engineering rigor while being accessible to users. An objective of this research is to further support the development and application SED2K for sediment phosphorus diagenesis and release to the water column of Onondaga Lake. Application of SED2K has been made to eutrophic Lake Alice in Minnesota. The more homogenous sediment characteristics of Lake Alice, compared with the industrially polluted sediment layers of Onondaga Lake, allowed for an invariant rate coefficient to be applied to describe first order decay kinetics of phosphorus. When a similar approach was attempted on Onondaga Lake an invariant rate coefficient failed to simulate the sediment phosphorus profile. Therefore, labile P was accounted for by progressive preservation after burial and a rate coefficient which gradual decreased with depth was applied. In this study, profile sediment samples were chemically extracted into five operationally-defined fractions: CaCO3-P, Fe/Al-P, Biogenic-P, Ca Mineral-P and Residual-P. Chemical fractionation data, from this study, showed that preservation is not the only mechanism by which phosphorus may be maintained in a non-reactive state in the profile. Sorption has been shown to contribute substantially to P burial within the profile. A new kinetic approach involving partitioning of P into process based fractions is applied here. Results from this approach indicate that labile P (Ca Mineral and Organic P) is contributing to internal P loading to Onondaga Lake, through diagenesis and diffusion to the water column, while the sorbed P fraction (Fe/Al-P and CaCO3-P) is remaining consistent. Sediment profile concentrations of labile and total phosphorus at time of deposition were also modeled and compared with current labile and total phosphorus, to quantify the extent to which remaining phosphorus which will continue to contribute to internal P loading and influence the trophic status of Onondaga Lake. Results presented here also allowed for estimation of the depth of the active sediment layer and the attendant response time as well as the sediment burden of labile P and associated efflux.

INTEGRATED ASSESSMENT OF ANTHROPOGENIC, CLIMATE, AND POLICY INDUCED CHANGES OF PHOSPHORUS EXPORT IN THE UNITED STATES LAURENTIAN GREAT LAKES WATERSHEDS

Anthropogenic activities have increased phosphorus (P) loading in tributaries to the Laurentian Great Lakes resulting in eutrophication in small bays to most notably, Lake Erie. Changes to surface water quality from P loading have resulted in billions of dollars in damage and threaten the health of the world’s largest freshwater resource. To understand the factors affecting P delivery with projected increasing urban lands and biofuels expansion, two spatially explicit models were coupled. The coupled models predict that the majority of the basin will experience a significant increase in urban area P sources while the agriculture intensity and forest sources of P will decrease. Changes in P loading across the basin will be highly variable spatially. Additionally, the impacts of climate change on high precipitation events across the Great Lakes were examined. Using historical regression relationships on phosphorus concentrations, key Great Lakes tributaries were found to have future changes including decreasing total loads and increases to high-flow loading events. The urbanized Cuyahoga watersheds exhibits the most vulnerability to these climate-induced changes with increases in total loading and storm loading , while the forested Au Sable watershed exhibits greater resilience. Finally, the monitoring network currently in place for sampling the amount of phosphorus entering the U.S. Great Lakes was examined with a focus on the challenges to monitoring. Based on these interviews, the research identified three issues that policy makers interested in maintaining an effective phosphorus monitoring network in the Great Lakes should consider: first, that the policy objectives driving different monitoring programs vary, which results in different patterns of sampling design and frequency; second, that these differences complicate efforts to encourage collaboration; and third, that methods of funding sampling programs vary from agency to agency, further complicating efforts to generate sufficient long-term data to improve our understanding of phosphorus into the Great Lakes. The dissertation combines these three areas of research to present the potential future impacts of P loading in the Great Lakes as anthropogenic activities, climate and monitoring changes. These manuscripts report new experimental data for future sources, loading and climate impacts on phosphorus.

Carbon in the Peatlands in the Great Lakes Region

Peatlands cover only ~3% of the global land area, but store ~30% of the worlds' soil carbon. There are many different peat types that store different amounts of carbon. Most inventories of carbon storage in northern peatlands have been conducted in the expansive Sphagnum dominated peatlands. Although, northern white cedar peatlands (NW cedar, Thuja occidentalis L.) are also one of the most common peatland types in the Great Lakes Region, occupying more than 2 million hectares. NW cedar swamps are understudied, due in part to the difficulties in collection methods. General lack of rapid and consistent sampling methods has also contributed in a lack of carbon stock quantification for many peatlands. The main objective of this thesis is to quantify: 1) to evaluate peat sampling methods 2) the amount of C-stored and the rates of long-term carbon accumulation in NW cedar peatlands. We sampled 38 peatlands separated into four categories (black ash, NW cedar swamp, sedge, and Sphagnum) during the summers of 2011/2012 across northern MN and the Upper Peninsula of MI. Basal dates of peat indicate that cedar peatlands were between 1970-7790 years old. Cedar peatlands are generally shallower than Sphagnum peat, but due to their higher bulk density, hold similar amounts of carbon with our sites averaging ~800 MgC ha-1. We estimate that NW cedar peatlands store over 1.7 Gt of carbon in the Great Lakes Region. Each of the six methods evaluated had a different level of accuracy and requires varying levels of effort and resources. The depth only method and intermittent sampling method were the most accurate methods of peatland sampling.

Evaluation of sugar maple dieback in the Upper Great Lakes region and development of a forest health youth education program

Acer saccharum Marsh., is one of the most valuable trees in the northern hardwood forests. Severe dieback was recently reported by area foresters in the western Upper Great Lakes Region. Sugar Maple has had a history of dieback over the last 100 years throughout its range and different variables have been identified as being the predisposing and inciting factors in different regions at different times. Some of the most common factors attributed to previous maple dieback episodes were insect defoliation outbreaks, inadequate precipitation, poor soils, atmospheric deposition, fungal pathogens, poor management, or a combination of these. The current sugar maple dieback was evaluated to determine the etiology, severity, and change in dieback on both industry and public lands. A network of 120 sugar maple health evaluation plots was established in the Upper Peninsula, Michigan, northern Wisconsin, and eastern Minnesota and evaluated annually from 2009-2012. Mean sugar maple crown dieback between 2009-2012 was 12.4% (ranging from 0.8-75.5%) across the region. Overall, during the sampling period, mean dieback decreased by 5% but individual plots and trees continued to decline. Relationships were examined between sugar maple dieback and growth, habitat conditions, ownership, climate, soil, foliage nutrients, and the maple pathogen sapstreak. The only statistically significant factor was found to be a high level of forest floor impacts due to exotic earthworm activity. Sugar maple on soils with lower pH had less earthworm impacts, less dieback, and higher growth rates than those on soils more favorable to earthworms. Nutritional status of foliage and soil was correlated with dieback and growth suggesting perturbation of nutrient cycling may be predisposing or contributing to dieback. The previous winter's snowfall totals, length of stay on the ground, and number of days with freezing temperatures had a significant positive relationship to sugar maple growth rates. Sapstreak disease, Ceratocystis virescens, may be contributing to dieback in some stands but was not related to the amount of dieback in the region. The ultimate goal of this research is to help forest managers in the Great Lakes Region prevent, anticipate, reduce, and/or salvage stands with dieback and loss in the future. An improved understanding of the complex etiology associated with sugar maple dieback in the Upper Great Lakes Region is necessary to make appropriate silvicultural decisions. Forest Health education helps increase awareness and proactive forest management in the face of changing forest ecosystems. Lessons are included to assist educators in incorporating forest health into standard biological disciplines at the secondary school curricula.

SYSTEM DYNAMICS MODELING AS A QUANTITATIVE-QUALITATIVE FRAMEWORK FOR SUSTAINABLE WATER RESOURCES MANAGEMENT: INSIGHTS FOR WATER QUALITY POLICY IN THE GREAT LAKES REGION

Early water resources modeling efforts were aimed mostly at representing hydrologic processes, but the need for interdisciplinary studies has led to increasing complexity and integration of environmental, social, and economic functions. The gradual shift from merely employing engineering-based simulation models to applying more holistic frameworks is an indicator of promising changes in the traditional paradigm for the application of water resources models, supporting more sustainable management decisions. This dissertation contributes to application of a quantitative-qualitative framework for sustainable water resources management using system dynamics simulation, as well as environmental systems analysis techniques to provide insights for water quality management in the Great Lakes basin. The traditional linear thinking paradigm lacks the mental and organizational framework for sustainable development trajectories, and may lead to quick-fix solutions that fail to address key drivers of water resources problems. To facilitate holistic analysis of water resources systems, systems thinking seeks to understand interactions among the subsystems. System dynamics provides a suitable framework for operationalizing systems thinking and its application to water resources problems by offering useful qualitative tools such as causal loop diagrams (CLD), stock-and-flow diagrams (SFD), and system archetypes. The approach provides a high-level quantitative-qualitative modeling framework for "big-picture" understanding of water resources systems, stakeholder participation, policy analysis, and strategic decision making. While quantitative modeling using extensive computer simulations and optimization is still very important and needed for policy screening, qualitative system dynamics models can improve understanding of general trends and the root causes of problems, and thus promote sustainable water resources decision making. Within the system dynamics framework, a growth and underinvestment (G&U) system archetype governing Lake Allegan's eutrophication problem was hypothesized to explain the system's problematic behavior and identify policy leverage points for mitigation. A system dynamics simulation model was developed to characterize the lake's recovery from its hypereutrophic state and assess a number of proposed total maximum daily load (TMDL) reduction policies, including phosphorus load reductions from point sources (PS) and non-point sources (NPS). It was shown that, for a TMDL plan to be effective, it should be considered a component of a continuous sustainability process, which considers the functionality of dynamic feedback relationships between socio-economic growth, land use change, and environmental conditions. Furthermore, a high-level simulation-optimization framework was developed to guide watershed scale BMP implementation in the Kalamazoo watershed. Agricultural BMPs should be given priority in the watershed in order to facilitate cost-efficient attainment of the Lake Allegan's TP concentration target. However, without adequate support policies, agricultural BMP implementation may adversely affect the agricultural producers. Results from a case study of the Maumee River basin show that coordinated BMP implementation across upstream and downstream watersheds can significantly improve cost efficiency of TP load abatement.

Pesky Pests of the Great Lakes State: Is Public Participation Influenced by Geographic Differences?

In Michigan, environmental issues, such as invasive species, are not geographically constrained, affecting citizens throughout the state. Regulations and management plans organized by scientists and officials are intended to address issues statewide, but these policies may not adequately tackle the threat from invasive species as it impacts different parts of the state at different times. Participation and contributions from citizens can offer insight into the impacts and changes non-native species have on the local ecosystem. However, chances to participate and contribute may be influenced by geographic location in the state. To understand if this was the case, this research studied publicly available documents and completed participant observations and semistructured interviews with participants, leaders, and officials included in invasive species management. Between the two study locations, Metro Detroit and the Western Upper Peninsula of Michigan, locational differences had some impact on opportunities to contribute to invasive species management. Population and the differences in the type of advertising used to alert citizens about events influenced access to participation opportunities. This research also revealed that this public policy issue lacks public involvement and contributions. Between the two locations, more involvement opportunities and organizations were present in Metro Detroit. However, it was the organizations themselves and their limited political involvement, and not geographic location, which had a greater impact on citizens' lack of participation in invasive species management.

Development of measurement and modeling techniques to quantify atmospheric deposition of persistent, bioaccumulative and toxic chemicals in the Great Lakes

Measurement and modeling techniques were developed to improve over-water gaseous air-water exchange measurements for persistent bioaccumulative and toxic chemicals (PBTs). Analytical methods were applied to atmospheric measurements of hexachlorobenzene (HCB), polychlorinated biphenyls (PCBs), and polybrominated diphenyl ethers (PBDEs). Additionally, the sampling and analytical methods are well suited to study semivolatile organic compounds (SOCs) in air with applications related to secondary organic aerosol formation, urban, and indoor air quality. A novel gas-phase cleanup method is described for use with thermal desorption methods for analysis of atmospheric SOCs using multicapillary denuders. The cleanup selectively removed hydrogen-bonding chemicals from samples, including much of the background matrix of oxidized organic compounds in ambient air, and thereby improved precision and method detection limits for nonpolar analytes. A model is presented that predicts gas collection efficiency and particle collection artifact for SOCs in multicapillary denuders using polydimethylsiloxane (PDMS) sorbent. An approach is presented to estimate the equilibrium PDMS-gas partition coefficient (Kpdms) from an Abraham solvation parameter model for any SOC. A high flow rate (300 L min-1) multicapillary denuder was designed for measurement of trace atmospheric SOCs. Overall method precision and detection limits were determined using field duplicates and compared to the conventional high-volume sampler method. The high-flow denuder is an alternative to high-volume or passive samplers when separation of gas and particle-associated SOCs upstream of a filter and short sample collection time are advantageous. A Lagrangian internal boundary layer transport exchange (IBLTE) Model is described. The model predicts the near-surface variation in several quantities with fetch in coastal, offshore flow: 1) modification in potential temperature and gas mixing ratio, 2) surface fluxes of sensible heat, water vapor, and trace gases using the NOAA COARE Bulk Algorithm and Gas Transfer Model, 3) vertical gradients in potential temperature and mixing ratio. The model was applied to interpret micrometeorological measurements of air-water exchange flux of HCB and several PCB congeners in Lake Superior. The IBLTE Model can be applied to any scalar, including water vapor, carbon dioxide, dimethyl sulfide, and other scalar quantities of interest with respect to hydrology, climate, and ecosystem science.