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.

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.