David Schwab - Academia.edu (original) (raw)

Papers by David Schwab

Research paper thumbnail of A hybrid Lagrangian–Eulerian particle model for reacting pollutant dispersion in non-homogeneous non-isotropic turbulence

Physica A: Statistical Mechanics and its Applications, 2009

Current numerical methods for simulating biophysical processes in aquatic environments are typica... more Current numerical methods for simulating biophysical processes in aquatic environments are typically constructed in a grid-based Eulerian framework or as an individual-based model in a particle-based Lagrangian framework. Often, the biogeochemical processes and physical (hydrodynamic) processes occur at different time and space scales, and changes in biological processes do not affect the hydrodynamic conditions. Therefore, it is possible to develop an alternative strategy to grid-based approaches for linking hydrodynamic and biogeochemical models that can significantly improve computational efficiency for this type of linked biophysical model. In this work, we utilize a new technique that links hydrodynamic effects and biological processes through a property-carrying particle model (PCPM) in a Lagrangian/Eulerian framework. The model is tested in idealized cases and its utility is demonstrated in a practical application to Sandusky Bay. Results show the integration of Lagrangian and Eulerian approaches allows for a natural coupling of mass transport (represented by particle movements and random walk) and biological processes in water columns which is described by a nutrient-phytoplankton-zooplankton-detritus (NPZD) biological model. This method is far more efficient than traditional tracer-based Eulerian biophysical models for 3-D simulation, particularly for a large domain and/or ensemble simulations.

Research paper thumbnail of Particle Transport, Nutrient Cycling, and Algal Community Structure Associated with a Major Winter-Spring Sediment Resuspension Event in Southern Lake Michigan

Journal of Great Lakes Research, 2002

Over the past decade, intermittent satellite imagery revealed the presence of an extensive plume ... more Over the past decade, intermittent satellite imagery revealed the presence of an extensive plume of resuspended sediments in late winter-early spring with a clear offshore projection coinciding with the region of maximum sediment accumulation in the lake. The large scale of the plume implied that this process was important in sediment, and associated constituent, cycling and transport, but it had never been sampled due to severe conditions. The onset of the 1996 event coincided with a major March storm. Within a few days the plume was approximately 10 km wide and over 300 km in length, implying that the source of the reflective materials was widely distributed. An estimate of the total mass of resuspended sediment, 12 days after the storm, was similar to the annual external load of (sand-free) particulate material to the southern basin. The high turbidity plume persisted for over a month, progressing northward along the eastern shore with a major offshore transport feature. Sediment traps within this feature recorded a major mass flux event. The plume was sampled on two occasions and was found to contain 5 to 10 times as much suspended matter as open-lake locations outside the visible plume. Total particulate phosphorus was high within the plume making this episodic process important in sedimentwater exchange. The diatom community structure within the plume was significantly different from outside the plume and was characteristic of more eutrophic waters. Abundance of non-diatom phytoplankton and microbial food web organisms were highest at the plume edge. The episodic nature of this process makes it difficult to sample, but the scale makes it important in designing monitoring programs and massbalance modeling efforts.

Research paper thumbnail of Modeling wind-driven circulation during the March 1998 sediment resuspension event in Lake Michigan

Journal of Geophysical Research: Oceans, 2003

A three-dimensional primitive equation numerical ocean model was applied to Lake Michigan to simu... more A three-dimensional primitive equation numerical ocean model was applied to Lake Michigan to simulate hydrodynamic conditions during the March 1998 sediment resuspension event in southern Lake Michigan caused by a storm with winds up to 20 m/s. The hydrodynamic model is driven with surface winds derived from observed meteorological conditions at 18 land stations and a meteorological buoy and also with surface winds calculated using a mesoscale meteorological model. Current observations from 11 subsurface moorings showed that the model driven with observed winds was able to qualitatively simulate wind-driven currents but underestimated current speeds during the most significant wind event. In addition, a pronounced offshore flow in the area of observations was also underestimated. Hydrodynamic model results using the meteorological model winds as the forcing function showed significant improvement over model results which were based on observed winds proving the importance of mesoscale winds for current modeling in large lakes.

Research paper thumbnail of A model study of the coupled biological and physical dynamics in Lake Michigan

Ecological Modelling, 2002

A coupled physical and biological model was developed for Lake Michigan. The physical model was t... more A coupled physical and biological model was developed for Lake Michigan. The physical model was the Princeton ocean model (POM) driven directly by observed winds and net surface heat flux. The biological model was an eight-component, phosphorus-limited, lower trophic level food web model, which included phosphate and silicate for nutrients, diatoms and non-diatoms for dominant phytoplankton species, copepods and protozoa for dominant zooplankton species, bacteria and detritus. Driven by observed meteorological forcings, a 1-D modeling experiment showed a controlling of physical processes on the seasonal variation of biological variables in Lake Michigan: diatoms grew significantly in the subsurface region in early summer as stratification developed and then decayed rapidly in the surface mixed layer when silicate supplied from the deep stratified region was reduced as a result of the formation of the thermocline. The non-diatoms subsequently grew in mid and late summer under a limited-phosphate environment and then declined in the fall and winter as a result of the nutrient consumption in the upper eutrophic layer, limitation of nutrients supplied from the deep region and meteorological cooling and wind mixing. The flux estimates suggested that the microbial loop had a significant contribution in the growth of microzooplankton and hence, to the lower-trophic level food web system. The model results agreed with observations, suggesting that the

Research paper thumbnail of Influences of suspended sediments on the ecosystem in Lake Michigan: a 3-D coupled bio-physical modeling experiment

Ecological Modelling, 2002

The influence of suspended sediments on the Lake Michigan ecosystem was examined using a 3-dimens... more The influence of suspended sediments on the Lake Michigan ecosystem was examined using a 3-dimensional (3-D) coupled biological and physical model developed by Chen et al. (part I). The model was driven by the realistic meteorological forces observed in March 1998, with daily inputs of suspended sediment concentration that were derived from temporally and spatially interpolated satellite imagery. The model results show the significant impact of a seasonally recurring coastal resuspension plume on the spatial and temporal variation of the nutrients and plankton in southern Lake Michigan. The plume-released nutrients played an essential role in maintaining the nutrient level in the lake. Although the growth of phytoplankton in the plume depended on the availability of nutrients and light, the offshore decrease in phytoplankton biomass still satisfied the Sverdrup's relationship. Cross-shore fluxes of nutrients and phytoplankton were controlled by episodic wind events with a period of : 5-7 days: offshore during southward winds and onshore during northward winds. The flux estimates for biological variables suggest that the microbial food web is a key contributor to secondary production in southern Lake Michigan and the lower trophic level food web system could be dynamically divided into two decoupled loops: (1) detritus-bacteria-microzooplankton-large zooplankton; and (2) nutrient-phytoplankton-detritus. The model-predicted spatial distributions of nutrients and phytoplankton were in reasonable agreement with observations taken during the 1998 EEGLE interdisciplinary cruises, suggesting that the model was sufficiently robust to capture the basic characteristics of the Lake Michigan ecosystem during the plume event.

Research paper thumbnail of Impact of inertial period internal waves on fixed-depth primary production estimates

Journal of Plankton Research, 1988

ABSTRACT Fixed-depth subthermocline phytoplankton production estimates were compared to variable-... more ABSTRACT Fixed-depth subthermocline phytoplankton production estimates were compared to variable-depth estimates calculated along the simulated trajectory of an internal wave. The mean of 5000 Monte Carlo simulations, varying wave phase and amplitude, was not significantly different than the fixed-depth estimate for average Lake Michigan internal waves. Differences were significant, however, for wave amplitudes >5 m but only at some depths. Differences between the two estimates were related to differences in irradiance received and the portion of the photosynthesisirradiance ( P-I ) curve controlling production. Oscillating communities always receive more irradiance than fixed-depth communities and the magnitude of this increase is related to the extinction coefficient and the amplitude of the internal waves. Production was also estimated along an individual isotherm trajectory (isotherm-derived) and compared to fixed-depth production. Larger differences between these isotherm-derived and fixed-depth estimates were noted in some cases and were related to differences in the mean isotherm depth and the sampled (fixed) depth. If one accounts for the trajectory of the sampled community, fixed-depth estimates are reliable; however, if the trajectory is unknown or unaccounted for, any individual fixed-depth production estimate may not adequately measure in situ production.

Research paper thumbnail of A hybrid Lagrangian–Eulerian particle model for reacting pollutant dispersion in non-homogeneous non-isotropic turbulence

Physica A: Statistical Mechanics and its Applications, 2009

Current numerical methods for simulating biophysical processes in aquatic environments are typica... more Current numerical methods for simulating biophysical processes in aquatic environments are typically constructed in a grid-based Eulerian framework or as an individual-based model in a particle-based Lagrangian framework. Often, the biogeochemical processes and physical (hydrodynamic) processes occur at different time and space scales, and changes in biological processes do not affect the hydrodynamic conditions. Therefore, it is possible to develop an alternative strategy to grid-based approaches for linking hydrodynamic and biogeochemical models that can significantly improve computational efficiency for this type of linked biophysical model. In this work, we utilize a new technique that links hydrodynamic effects and biological processes through a property-carrying particle model (PCPM) in a Lagrangian/Eulerian framework. The model is tested in idealized cases and its utility is demonstrated in a practical application to Sandusky Bay. Results show the integration of Lagrangian and Eulerian approaches allows for a natural coupling of mass transport (represented by particle movements and random walk) and biological processes in water columns which is described by a nutrient-phytoplankton-zooplankton-detritus (NPZD) biological model. This method is far more efficient than traditional tracer-based Eulerian biophysical models for 3-D simulation, particularly for a large domain and/or ensemble simulations.

Research paper thumbnail of Particle Transport, Nutrient Cycling, and Algal Community Structure Associated with a Major Winter-Spring Sediment Resuspension Event in Southern Lake Michigan

Journal of Great Lakes Research, 2002

Over the past decade, intermittent satellite imagery revealed the presence of an extensive plume ... more Over the past decade, intermittent satellite imagery revealed the presence of an extensive plume of resuspended sediments in late winter-early spring with a clear offshore projection coinciding with the region of maximum sediment accumulation in the lake. The large scale of the plume implied that this process was important in sediment, and associated constituent, cycling and transport, but it had never been sampled due to severe conditions. The onset of the 1996 event coincided with a major March storm. Within a few days the plume was approximately 10 km wide and over 300 km in length, implying that the source of the reflective materials was widely distributed. An estimate of the total mass of resuspended sediment, 12 days after the storm, was similar to the annual external load of (sand-free) particulate material to the southern basin. The high turbidity plume persisted for over a month, progressing northward along the eastern shore with a major offshore transport feature. Sediment traps within this feature recorded a major mass flux event. The plume was sampled on two occasions and was found to contain 5 to 10 times as much suspended matter as open-lake locations outside the visible plume. Total particulate phosphorus was high within the plume making this episodic process important in sedimentwater exchange. The diatom community structure within the plume was significantly different from outside the plume and was characteristic of more eutrophic waters. Abundance of non-diatom phytoplankton and microbial food web organisms were highest at the plume edge. The episodic nature of this process makes it difficult to sample, but the scale makes it important in designing monitoring programs and massbalance modeling efforts.

Research paper thumbnail of Modeling wind-driven circulation during the March 1998 sediment resuspension event in Lake Michigan

Journal of Geophysical Research: Oceans, 2003

A three-dimensional primitive equation numerical ocean model was applied to Lake Michigan to simu... more A three-dimensional primitive equation numerical ocean model was applied to Lake Michigan to simulate hydrodynamic conditions during the March 1998 sediment resuspension event in southern Lake Michigan caused by a storm with winds up to 20 m/s. The hydrodynamic model is driven with surface winds derived from observed meteorological conditions at 18 land stations and a meteorological buoy and also with surface winds calculated using a mesoscale meteorological model. Current observations from 11 subsurface moorings showed that the model driven with observed winds was able to qualitatively simulate wind-driven currents but underestimated current speeds during the most significant wind event. In addition, a pronounced offshore flow in the area of observations was also underestimated. Hydrodynamic model results using the meteorological model winds as the forcing function showed significant improvement over model results which were based on observed winds proving the importance of mesoscale winds for current modeling in large lakes.

Research paper thumbnail of A model study of the coupled biological and physical dynamics in Lake Michigan

Ecological Modelling, 2002

A coupled physical and biological model was developed for Lake Michigan. The physical model was t... more A coupled physical and biological model was developed for Lake Michigan. The physical model was the Princeton ocean model (POM) driven directly by observed winds and net surface heat flux. The biological model was an eight-component, phosphorus-limited, lower trophic level food web model, which included phosphate and silicate for nutrients, diatoms and non-diatoms for dominant phytoplankton species, copepods and protozoa for dominant zooplankton species, bacteria and detritus. Driven by observed meteorological forcings, a 1-D modeling experiment showed a controlling of physical processes on the seasonal variation of biological variables in Lake Michigan: diatoms grew significantly in the subsurface region in early summer as stratification developed and then decayed rapidly in the surface mixed layer when silicate supplied from the deep stratified region was reduced as a result of the formation of the thermocline. The non-diatoms subsequently grew in mid and late summer under a limited-phosphate environment and then declined in the fall and winter as a result of the nutrient consumption in the upper eutrophic layer, limitation of nutrients supplied from the deep region and meteorological cooling and wind mixing. The flux estimates suggested that the microbial loop had a significant contribution in the growth of microzooplankton and hence, to the lower-trophic level food web system. The model results agreed with observations, suggesting that the

Research paper thumbnail of Influences of suspended sediments on the ecosystem in Lake Michigan: a 3-D coupled bio-physical modeling experiment

Ecological Modelling, 2002

The influence of suspended sediments on the Lake Michigan ecosystem was examined using a 3-dimens... more The influence of suspended sediments on the Lake Michigan ecosystem was examined using a 3-dimensional (3-D) coupled biological and physical model developed by Chen et al. (part I). The model was driven by the realistic meteorological forces observed in March 1998, with daily inputs of suspended sediment concentration that were derived from temporally and spatially interpolated satellite imagery. The model results show the significant impact of a seasonally recurring coastal resuspension plume on the spatial and temporal variation of the nutrients and plankton in southern Lake Michigan. The plume-released nutrients played an essential role in maintaining the nutrient level in the lake. Although the growth of phytoplankton in the plume depended on the availability of nutrients and light, the offshore decrease in phytoplankton biomass still satisfied the Sverdrup's relationship. Cross-shore fluxes of nutrients and phytoplankton were controlled by episodic wind events with a period of : 5-7 days: offshore during southward winds and onshore during northward winds. The flux estimates for biological variables suggest that the microbial food web is a key contributor to secondary production in southern Lake Michigan and the lower trophic level food web system could be dynamically divided into two decoupled loops: (1) detritus-bacteria-microzooplankton-large zooplankton; and (2) nutrient-phytoplankton-detritus. The model-predicted spatial distributions of nutrients and phytoplankton were in reasonable agreement with observations taken during the 1998 EEGLE interdisciplinary cruises, suggesting that the model was sufficiently robust to capture the basic characteristics of the Lake Michigan ecosystem during the plume event.

Research paper thumbnail of Impact of inertial period internal waves on fixed-depth primary production estimates

Journal of Plankton Research, 1988

ABSTRACT Fixed-depth subthermocline phytoplankton production estimates were compared to variable-... more ABSTRACT Fixed-depth subthermocline phytoplankton production estimates were compared to variable-depth estimates calculated along the simulated trajectory of an internal wave. The mean of 5000 Monte Carlo simulations, varying wave phase and amplitude, was not significantly different than the fixed-depth estimate for average Lake Michigan internal waves. Differences were significant, however, for wave amplitudes >5 m but only at some depths. Differences between the two estimates were related to differences in irradiance received and the portion of the photosynthesisirradiance ( P-I ) curve controlling production. Oscillating communities always receive more irradiance than fixed-depth communities and the magnitude of this increase is related to the extinction coefficient and the amplitude of the internal waves. Production was also estimated along an individual isotherm trajectory (isotherm-derived) and compared to fixed-depth production. Larger differences between these isotherm-derived and fixed-depth estimates were noted in some cases and were related to differences in the mean isotherm depth and the sampled (fixed) depth. If one accounts for the trajectory of the sampled community, fixed-depth estimates are reliable; however, if the trajectory is unknown or unaccounted for, any individual fixed-depth production estimate may not adequately measure in situ production.