Ecology of larval fishes and large zooplankton in the Keweenaw Current region of Lake Superior, with special focus on lake herring, Coregonus artedi

I assessed the influence of the Keweenaw Current and spring thermal bar on the distribution of larval fishes and large zooplankton in Lake Superior. In 1998 and 1999, samples were collected from inshore (0.2 – 3.0 km from shore) and offshore (5.0 – 9.0 km from shore) locations on three transects off the western coast of the Keweenaw Peninsula, Michigan. For larval fishes, density and size distribution patterns of lake herring (Coregonus artedi), rainbow smelt (Osmerus mordax), burbot (Lota lota), deepwater sculpin (Myoxocephalus thompsoni), and spoonhead sculpin (Cottus ricei) suggest a seasonal inshore to offshore movement. For zooplankton, seasonal warming appeared to be the major factor that limited planktonic catches of the primarily benthic Mysisrelicta and Diporeia spp., while simultaneously stimulated growth and reproduction of the cladocerans Daphnia spp., Holopedium gibberum, and Bythotrephes cederstroemi. In contrast, calanoid copepods as a group were abundant throughout the entire sampling season. The greatest abundances of zooplankton were generally encountered offshore, even for the cladocerans, which apparently expanded from inshore to offshore locations with seasonal warming. In 2000, sampling efforts focused on lake herring. Samples were collected from surface waters at 0.1 – 17.0 km from shore on two transects. Lake herring larvae were also reared in the laboratory from eggs in order to validate the use of otolith microstructure for aging. Increment deposition was not statistically different from a daily rate starting from 28 days after hatching, near the time of yolk-sac absorption, but larvae with lower growth rates could not be aged as accurately. In Lake Superior, lake herring tended to be slightly more abundant, larger, and older at inshore locations, but a dense patch of younger larvae was also encountered 7 – 13 km from shore. The distribution iiipatterns suggest that larvae were transported by prevailing currents into the study region, possibly from the more productive spawning regions in western Lake Superior. Growth rates were suppressed at offshore locations where temperatures were less than 8°C. These results indicate that lake herring larvae may be transported far distances from spawning concentrations by longshore currents, and water temperatures may largely control their growth.

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.

Study of particle movement in the nearshore region of Lake Superior with radioisotope tracers

Biogeochemical processes in the coastal region, including the coastal area of the Great Lakes, are of great importance due to the complex physical, chemical and biological characteristics that differ from those on either the adjoining land or open water systems. Particle-reactive radioisotopes, both naturally occurring (210Pb, 210Po and 7Be) and man-made (137Cs), have proven to be useful tracers for these processes in many systems. However, a systematic isotope study on the northwest coast of the Keweenaw Peninsula in Lake Superior has not yet been performed. In this dissertation research, field sampling, laboratory measurements and numerical modeling were conducted to understand the biogeochemistry of the radioisotope tracers and some particulate-related coastal processes. In the first part of the dissertation, radioisotope activities of 210Po and 210Pb in a variability of samples (dissolved, suspended particle, sediment trap materials, surficial sediment) were measured. A completed picture of the distribution and disequilibrium of this pair of isotopes was drawn. The application of a simple box model utilizing these field observations reveals short isotope residence times in the water column and a significant contribution of sediment resuspension (for both particles and isotopes). The results imply a highly dynamic coastal region. In the second part of this dissertation, this conclusion is examined further. Based on intensive sediment coring, the spatial distribution of isotope inventories (mainly 210Pb, 137Cs and 7Be) in the nearshore region was determined. Isotope-based focusing factors categorized most of the sampling sites as non- or temporary depositional zones. A twodimensional steady-state box-in-series model was developed and applied to individual transects with the 210Pb inventories as model input. The modeling framework included both water column and upper sediments down to the depth of unsupported 210Pb penetration. The model was used to predict isotope residence times and cross-margin fluxes of sediments and isotopes at different locations along each transect. The time scale for sediment focusing from the nearshore to offshore regions of the transect was on the order of 10 years. The possibility of sediment longshore movement was indicated by high inventory ratios of 137Cs: 210Pb. Local deposition of fine particles, including fresh organic carbon, may explain the observed distribution of benthic organisms such as Diporeia. In the last part of this dissertation, isotope tracers, 210Pb and 210Po, were coupled into a hydrodynamic model for Lake Superior. The model was modified from an existing 2-D finite difference physical-biological model which has previously been successfully applied on Lake Superior. Using the field results from part one of this dissertation as initial conditions, the model was used to predict the isotope distribution in the water column; reasonable results were achieved. The modeling experiments demonstrated the potential for using a hydrodynamic model to study radioisotope biogeochemistry in the lake, although further refinements are necessary.

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.

Paleomagnetic Study of the Portage Lake Volcanics exposed in the Quincy Mine

A detailed paleomagnetic and rock-magnetic investigation was conducted on thirty six basaltic flows of the ~1095 Ma Portage Lake Volcanics. The flows were sampled along the East Adit of the Quincy Mine (Hancock, MI). Thirty two flows yielded well-defined primary magnetization directions carried by magnetite. A secondary magnetization component carried by hematite was also found in twenty nine flows. After correction for serial correlation between the flows, nineteen independent mean directions were calculated. The corresponding paleomagnetic pole is located at 25.5 °N, 182.1 °W (A95 = 3.5°). The new pole overlaps with the pole from the ~1087 Ma Lake Shore Traps suggesting a standstill of the North American plate during that time period. The low angular dispersion of virtual geomagnetic poles (S = 7.9°) suggests that the flows were erupted within a short time period, or that the strength of geomagnetic secular variation was lower than that of the recent field.

A Balance Shifted: The Unionization of Michigan’s Lake Copper District, 1935-1955

Three decades after the unsuccessful 1913-1914 strike at the Lake District copper mines of Michigan, workers organized as Local 584 of the International Union of Mine, Mill, and Smelter Workers (Mine Mill) signed a union contract with Calumet & Hecla Consolidated Copper Company. C & H was the last and most significant of the region���s three major copper mining companies to unionize during the three-year period from 1939 to 1942. This paper tells the untold history of the successful union drives in the Lake District’s copper mines, starting with Copper Range Company in 1939 and encompassing the subsequent unionizations of Quincy Mining Company and finally C & H. The paper develops thematic connections between the 1913-1914, including Mine Mill’s lineage to the Western Federation of Miners, parallel ethnic dimensions, and, most significantly, the contrasting role of state authority between the two time periods. The paper carries the Lake District’s labor history forward to 1955 to include United Steelworkers’ successful challenge to Mine Mill in 1950 and the strike of 1955. This history also incorporates source material from the papers of highly influential union organizer and representative Eugene Saari, material which to date has not been integrated into the labor history of the region. This paper has not yet been submitted.

A STUDY OF SO2 EMISSIONS AND GROUND SURFACE DISPLACEMENTS AT LASTARRIA VOLCANO, ANTOFAGASTA REGION, NORTHERN CHILE

Lastarria volcano (Chile) is located at the North-West margin of the `Lazufre' ground inflation signal (37x45 km²), constantly uplifting at a rate of ~2.5 cm/year since 1996 (Pritchard and Simons 2002; Froger et al. 2007). The Lastarria volcano has the double interest to be superimposed on a second, smaller-scale inflation signal and to be the only degassing area of the Lazufre signal. In this project, we compared daily SO2 burdens recorded by AURA's OMI mission for 2005-2010 with Ground Surface Displacements (GSD) calculated from the Advanced Synthetic Aperture Radar (ASAR) images for 2003-2010. We found a constant maximum displacement rate of 2.44 cm/year for the period 2003-2007 and 0.80- 0.95 cm/year for the period 2007-2010. Total SO2 emitted is 67.0 kT for the period 2005-2010, but detection of weak SO2 degassing signals in the Andes remains challenging owing to increased noise in the South Atlantic radiation Anomaly region.

PROBABILISTIC MODELING OF RAINFALL INDUCED LANDSLIDE HAZARD ASSESSMENT IN SAN JUAN LA LAGUNA, SOLOLÁ, GUATEMALA

The municipality of San Juan La Laguna, Guatemala is home to approximately 5,200 people and located on the western side of the Lake Atitlán caldera. Steep slopes surround all but the eastern side of San Juan. The Lake Atitlán watershed is susceptible to many natural hazards, but most predictable are the landslides that can occur annually with each rainy season, especially during high-intensity events. Hurricane Stan hit Guatemala in October 2005; the resulting flooding and landslides devastated the Atitlán region. Locations of landslide and non-landslide points were obtained from field observations and orthophotos taken following Hurricane Stan. This study used data from multiple attributes, at every landslide and non-landslide point, and applied different multivariate analyses to optimize a model for landslides prediction during high-intensity precipitation events like Hurricane Stan. The attributes considered in this study are: geology, geomorphology, distance to faults and streams, land use, slope, aspect, curvature, plan curvature, profile curvature and topographic wetness index. The attributes were pre-evaluated for their ability to predict landslides using four different attribute evaluators, all available in the open source data mining software Weka: filtered subset, information gain, gain ratio and chi-squared. Three multivariate algorithms (decision tree J48, logistic regression and BayesNet) were optimized for landslide prediction using different attributes. The following statistical parameters were used to evaluate model accuracy: precision, recall, F measure and area under the receiver operating characteristic (ROC) curve. The algorithm BayesNet yielded the most accurate model and was used to build a probability map of landslide initiation points. The probability map developed in this study was also compared to the results of a bivariate landslide susceptibility analysis conducted for the watershed, encompassing Lake Atitlán and San Juan. Landslides from Tropical Storm Agatha 2010 were used to independently validate this study’s multivariate model and the bivariate model. The ultimate aim of this study is to share the methodology and results with municipal contacts from the author's time as a U.S. Peace Corps volunteer, to facilitate more effective future landslide hazard planning and mitigation.

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.