Geology ReportsSearch

SEARCH · Geology Reports

Results for “Offshore”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 577 records · Page 32Linked to original sources

Freshwater tidal forests and estuarine wetlands may confer early life growth advantages for delta-reared Chinook Salmon

Large river deltas are complex ecosystems that are believed to play a pivotal role in promoting the early marine growth and survival of threatened Chinook Salmon Oncorhynchus tshawytscha . We used a fish bioenergetics model to assess the functional role of multiple delta habitats across a gradient of salinities and vegetation types, where consumption and growth rate potential (GRP) were considered as proxies for habitat quality. We subsequently compared our model output to empirical, or realized , growth estimates from scale circuli. In terms of consumption, prey energy density (ED prey ) was 46–86% higher in tidal freshwater forest than in any other habitat type, while estimated consumption rates (expressed as proportion of maximum daily consumption; P cmax ) were positively correlated with FL. These size‐related differences in P cmax led to a noticeable increase along a freshwater‐to‐saline gradient from roughly 0.25 in tidal freshwater forest to 0.55 in the offshore subtidal zone, yet despite higher observed P cmax values in nearshore and offshore habitats, the tidal freshwater forest and emergent salt marsh demonstrated the highest modeled GRP values. Similarly, realized growth rates for fish caught in tidal freshwater forest were up to 0.5% higher per day than for fish caught in the offshore area, but habitat‐level differences were overshadowed by allometry and rearing origin. Scales from unmarked fish (assumed to be of wild origin) indicated that they grew, on average, 11% faster than did hatchery fish; however, these differences were subtle and were more obvious at fork lengths <100 mm. Our results suggest that tidal forests and emergent marshes may offer early life growth advantages for wild Chinook Salmon, but that wild and hatchery fish can compensate as they move seaward by opportunistically consuming greater quantities of low‐energy density prey, taking advantage of pulses of larval forage fish, or by spending time in multiple interconnected habitat types.

Transactions of the American Fisheries Society

Modeling three-dimensional flow over spur-and-groove morphology

Spur-and-groove (SAG) morphology characterizes the fore reef of many coral reefs worldwide. Although the existence and geometrical properties of SAG have been well documented, an understanding of the hydrodynamics over them is limited. Here, the three-dimensional flow patterns over SAG formations, and a sensitivity of those patterns to waves, currents, and SAG geometry were characterized using the physics-based Delft3D-FLOW and SWAN models. Shore-normal shoaling waves over SAG formations were shown to drive two circulation cells: a cell on the lower fore reef with offshore flow over the spurs and onshore flow over the grooves, except near the seabed where velocities were always onshore, and a cell on the upper fore reef with offshore surface velocities and onshore bottom currents, which result in depth-averaged onshore and offshore flow over the spurs and grooves, respectively. The mechanism driving this flow results from the net of the radiation stress gradients and pressure gradient, which is balanced by the Reynolds stress gradients and bottom friction that differ over the spur and over the groove. Waves were the primary driver of variations in modelled flow over SAG, with the flow strength increasing for increasing wave heights and periods. Spur height, SAG wavelength, and the water depth at peak spur height were the dominant influences on the hydrodynamics, with spur heights directly proportional to the strength of SAG circulation cells. SAG formations with shorter SAG wavelengths only presented one circulation cell on the shallower portion of the reef, as opposed to the two circulation cells for longer SAG wavelengths. SAG formations with peak spur heights occurring in shallower water had stronger circulation than those with peak spur heights occurring in deeper water. These hydrodynamic patterns also likely affect coral and reef development through sediment and nutrient fluxes.

Coral Reefs

Identification of fresh submarine groundwater off the coast of San Diego, USA, using electromagnetic methods

Climate change has a pronounced effect on water resources in many semiarid climates, causing populated areas such as San Diego County (USA), to become more vulnerable to water shortages in the coming decades. To prepare for decreased water supply, San Diego County is adopting policies to decrease water use and to develop additional local sources of water. One new local source of freshwater is produced by a desalination facility that purifies brackish groundwater from the coastal San Diego Formation. This formation has been studied extensively onshore, but little is known about the geology or groundwater quality offshore in the adjacent continental shelf. Because most groundwater systems are interconnected and complex, further analysis is needed to identify offshore geology, possible sequestration of freshwater in the shelf, and potential pathways for saltwater intrusion. This comprehensive understanding is important because seawater intrusion may limit use of the San Diego Formation and longevity of desalination facilities. Controlled-source electromagnetic methods are uniquely suited to detecting offshore groundwater as they are sensitive to changes in pore fluids such as the transition from fresh to brackish groundwater. This paper describes results from surface-towed electromagnetic surveys that mapped the pore-fluid salinity and possible fluid pathways in the continental shelf off the coast of San Diego. The results indicate a considerable volume of fresh-to-brackish groundwater sequestered in the shelf, both in continuous lenses and isolated pockets, that appear influenced by fault systems and shallow stratigraphy.

Baja California, California

Reconnaissance survey for potential energy storage and carbon dioxide storage resources of petroleum reservoirs in western Europe

Energy producers and utilities use oil and gas reservoirs for gas storage to meet peak seasonal demand or to supplement intermittent energy production. These reservoirs are also suitable for the long-term storage of carbon dioxide (CO 2 ), a greenhouse gas. This study reports on a reconnaissance analysis of the potential magnitude of storage resources in 9424 known oil and gas reservoirs from 24 countries within highly industrialized western Europe. To standardize the storage resources of the oil and gas reservoirs, their volumetric capacity is expressed in terms of metric tons (mass) of CO 2 . Estimates of recoverable oil and gas at the surface are converted to subsurface volumes and then converted to the equivalent mass of CO 2 at reservoir conditions. The results indicate 36.7 gigatons of CO 2 could be stored, with oil reservoirs accounting for 32% of that total and natural gas reservoirs comprising the remaining 68%. About four-fifths of the reservoir storage resource is offshore, with about three-fourths of that offshore resource at water depths of 200 m or less. Most countries do not have the reservoir storage resources to store 15 years of CO 2 at 2017 emission levels. With few exceptions the bulk of the storage is offshore for countries that do have at least 15 years of storage. The expansion of natural gas storage for strategic purposes in abandoned onshore gas reservoirs is not expected to seriously impact CO 2 storage. The contribution of this analysis is the description of the spatial distribution of potential storage and physical accessibility.

Natural Resources Research

Origins of seawater intrusion in a coastal aquifer - A case study of the Pajaro Valley, California

Seawater may enter and contaminate stratified coastal aquifers through a number of different pathways. These pathways and their relative contribution are examined in the Pajaro Valley, California, a coastal area with extensive groundwater development. This study considers three pathways of possible intrusion of the primary confined aquifer: (1) onshore leakage from brackish sources, the estuary and sloughs, through the confining layer; (2) near-shore leakage from the ocean through the confining layer; and (3) offshore flow from the ocean through the submarine canyon outcrop of the aquifer. Groundwater flow and seawater intrusion are simulated using an areal, two-dimensional solute-transport computer model. This analysis indicates that leakage through confining layers is the principal mechanism of recharge to the aquifer. Although lateral flow through the offshore outcrop contaminates the aquifer, as a whole, at a higher rate, vertical leakage through the sea floor initially is the main pathway of seawater intrusion to the onshore portion of the aquifer. It is likely that leakage generally is the dominant mechanism of recharge and initial cause of seawater intrusion for poorly-confined, stratified coastal aquifers. This analysis suggests that a significant time interval follows the initial observation of seawater intrusion, during which remedial action can be taken to control lateral flow through the offshore outcrop, which ultimately will be the largest component of future intrusion in these aquifers.

California

Arctic continental shelf morphology related to sea-ice zonation, Beaufort Sea, Alaska

Landsat-1 and NOAA satellite imagery for the winter 1972–1973, and a variety of ice and sea-floor data were used to study sea-ice zonation and dynamics and their relation to bottom morphology and geology on the Beaufort Sea continental shelf of arctic Alaska. In early winter the location of the boundary between undeformed fast ice and westward-drifting pack ice of the Pacific Gyre is controlled by major coastal promontories. Pronounced linear pressure- and shear-ridges, as well as hummock fields, form along this boundary and are stabilized by grounding, generally between the 10- and 20-m isobaths. Slippage along this boundary occurs intermittently at or seaward of the grounded ridges, forming new grounded ridges in a widening zone, the stamukhi zone, which by late winter extends out to the 40-m isobath. Between intermittent events along the stamukhi zone, pack-ice drift and slippage is continuous along the shelf edge, at average rates of 3–10 km/day. Whether slippage occurs along the stamukhi zone or along the shelf edge, it is restricted to a zone several hundred meters wide, and ice seaward of the slip face moves at uniform rates without discernible drag effects. A causal relationship is seen between the spatial distribution of major ice-ridge systems and offshore shoals downdrift of major coastal promontories. The shoals appear to have migrated shoreward under the influence of ice up to 400 m in the last 25 years. The sea floor seaward of these shoals within the stamukhi zone shows high ice-gouge density, large incision depths, and a high degree of disruption of internal sedimentary structures. The concentration of large ice ridges and our sea floor data in the stamukhi zone indicate that much of the available marine energy is expended here, while the inner shelf and coast, where the relatively undeformed fast ice grows, are sheltered. There is evidence that anomalies in the overall arctic shelf profile are related to sea-ice zonation, ice dynamics, and bottom processes. A proposed ice zonation, including zones of (1) bottom-fast ice, (2) floating fast ice, (3) stamukhi, and (4) seasonal pack ice, emphasizes ice interaction with the shelf surface and differs from previous zonation. Certain aspects of the results reported here are directly applicable to planned offshore developments in the Prudhoe Bay oil field. Properly placed artificial structures similar to offshore shoals should be able to withstand the forces of the ice, serve to modify the observed ice zonation, and might be used to make the environment less hostile to human activities.

Alaska

Stability of a very coarse-grained beach at Carmel, California

Monastery Beach at Carmel, California, is a pocket beach composed of very coarse to granular sediment. In profile, the beach has a well-defined berm crest; a steep foreshore; and a gently sloping, barless offshore covered by large, long-crested oscillation ripples. Carmel Submarine Canyon heads a few hundred meters offshore of the beach, and San Jose Creek, a small ephemeral steam, ponds onshore of the central part of the berm. Wave conditions vary greatly during a year because the beach lies open to the Pacific Ocean for azimuths between 270°–322°N whence come a variety of wave types. Even with a variable wave climate, Monastery Beach has maintained a swell profile for almost three years. Aperiodic beach surveys show that the beach responds little to seasonal changes in wave climate. Four survey lines maintained the same swell profile throughout the study period. The fifth line maintained a stable profile only across the foreshore; the berm was twice artificially breached during storms to prevent upstream flooding along San Jose Creek. In comparison, Carmel Beach, a nearby beach composed of medium sand, commonly alternates between swell and storm profiles. The increased stability of Monastery Beach relative to Carmel Beach is attributed to two factors: grain size differences and location within Carmel Bay. Rebuilding proceeded very slowly along the breached part of the berm at Monastery Beach. The probable cause of such a low recovery rate is that oscillation ripples trapped the sand that was carried offshore when San Jose Creek eroded the beach. The ripples, which are active under high-energy conditions, approach dormancy under low-energy conditions. Each ripple, therefore, acts like a reservoir, retaining sand during most swell conditions.

California

Tidal-cycle changes in oscillation ripples on the inner part of an estuarine sand flat

Oscillation ripples form on subaqueous sand beds when wave-generated, near-bottom water motions are strong enough to move sand grains. The threshold of grain motion is the lower bound of the regime of oscillation ripples and the onset of sheet flow is the upper bound. Based on the relation between ripple spacing and orbital diameter, three types of symmetrical ripples occur within the ripple regime. In the lower part of the ripple regime (orbital ripples), spacing is proportional to orbital diameter; in the upper part (anorbital ripples) spacing is independent of orbital diameter. Between these regions occurs a transitional region (suborbital ripples). Oscillation ripples develop on a sandy tidal flat in Willapa Bay, Washington, as a result of waves traversing the area when it is submerged. Because wave energy is usually low within the bay, the ripples are primarily orbital in type. This means that their spacing should respond in a systematic way to changes in wave conditions. During the high-water parts of some tidal cycles, ripples near the beach decrease in spacing during the latter stage of the ebb tide while ripples farther offshore do not change. Observations made over several tidal cycles show that the zone of active ripples shifts on- or offshore in response to different wave conditions. Detailed bed profiles and current measurements taken during the high-water part of spring tides show the manner in which the oscillation ripples change with changes in orbital diameter. Changes in ripple spacing at the study site could be correlated with changes in orbital diameter in the manner suggested by the criterion for orbital ripples. However, there appeared to be a lag time between a decrease in orbital diameter and the corresponding decrease in ripple spacing. Absence of change during a tidal cycle could be attributed to orbital velocities below the threshold for grain motion that negated the effects of changes in orbital diameter. Because changes in sand-flat ripples depend both upon changes in orbital diameter and upon the magnitude of the orbital velocity, exposed ripples were not necessarily produced during the preceding high tide. In fact, some ripples may have been just produced, while others, farther offshore, may have been produced an unknown number of tides earlier. Therefore, when interpreting past wave conditions over tidal flats from low-tide ripples, one must remember that wave periods have to be short enough to produce velocities greater than the threshold velocity for the orbital diameters calculated from the observed ripple spacings.

Washington

Modes of cross-shore sediment transport on the shoreface of the Middle Atlantic Bight

The mechanisms responsible for onshore and offshore sediment fluxes across the shoreface zone seaward of the surf zone were examined in a 3-year field study. The study was conducted in the southern part of the Middle Atlantic Bight in the depth region 7–17 m using instrumented tripods supporting electromagnetic current meters, pressure sensors, suspended sediment concentration sensors, and sonar altimeters. The observations embraced fairweather, moderate energy, swell-dominated, and storm conditions. Cross-shore mean flows ranged from near zero during fairweather to > 20 cm s −1 during the storm; oscillatory flows were on the order of 10 cm s −1 during fairweather and 100 cm s −1 during the storm. Suspended sediment concentrations at about 10 cm above the bed were < 0.1 kg m −3 under fairweather conditions, 1–2 kg m −3 under moderate swell conditions, and > 5 kg m −3 during the storm. Three methods were applied to evaluate the relative importance of incident waves, long-period oscillations, mean flows and gravity in effecting shoreward or seaward sediment flux: (1) an energetics transport model was applied to instantaneous near-bottom velocity data, (2) higher moments of near-bottom flows were estimated and compared, and (3) suspended sediment fluxes were estimated directly from the instantaneous products of cross-shore velocity and suspended sediment concentration. The results show that measurable contributions were made by all four of the processes. Most significantly, mean flows were seen to dominate and cause offshore fluxes during the storm and to contribute significantly to onshore and offshore flux during fairweather and moderate energy. Incident waves were, in all cases, the major source of bed shear stress but also caused shoreward as well as seaward net sediment advection. Low-frequency effects involving wave groups and long-period waves made secondary contributions to cross-shore sediment flux. Contrary to expectations, low-frequency fluxes were just as often shoreward as seaward. Whereas cross-correlations between suspended sediment concentration and the instantaneous near-bottom current speed were high and in phase under storm conditions, they were weak and out of phase during fairweather conditions. This suggests that simple energetics models are probably inadequate for predicting fairweather transport of suspended sediment.

Middle Atlantic Bight

Bottom stress estimates and sand transport on northern California inner continental shelf

Measurements of velocities and light transmission in the bottom boundary layer on the continental shelf off northern California demonstrate the importance of storms in the transport of sediment along the coast and offshore in this region. Time-series estimates of bottom stress obtained from a combined wave-current bottom boundary layer model in which wave and current measurements from the Geoprobe tripod were used as input show high stress values of 10 dynes cm −2 during two distinct storm events in early February and early March, 1991. These stresses induce significant offshore sediment transport, achieving maximum values of about 0.5 g cm −1 s −1 . The net transport over the entire measurement period from 30 January 1991 to 13 March 1991 was along the coast toward the north and offshore. This transport pattern explains slow migration of low amplitude, broad crescentic dunes along and across this portion of the inner continental shelf.

California

Nearshore bars and the break-point hypothesis

The set of hypotheses calling for bar formation at the break point was tested with field data. During two different experiments, waves were measured across the surf zone coincident with the development of a nearshore bar. We use a criterion, based on the wave height to depth ratio, to determine the offshore limit of the inner surf zone. During the first experiment, the bar became better developed and migrated offshore while remaining well within the inner surf zone. During the second experiment, the surf zone was narrower and we cannot rule out the possibility of break point processes contributing to bar development. We conclude that bars are not necessarily coupled with the break point and can become better developed and migrate offshore while being in the inner surf zone landward from initial wave breaking in the outer surf zone.

Coastal Engineering

Interpreting the spatio-temporal patterns of sea turtle strandings: Going with the flow

Knowledge of the spatial and temporal distribution of specific mortality sources is crucial for management of species that are vulnerable to human interactions. Beachcast carcasses represent an unknown fraction of at-sea mortalities. While a variety of physical (e.g., water temperature) and biological (e.g., decomposition) factors as well as the distribution of animals and their mortality sources likely affect the probability of carcass stranding, physical oceanography plays a major role in where and when carcasses strand. Here, we evaluate the influence of nearshore physical oceanographic and wind regimes on sea turtle strandings to decipher seasonal trends and make qualitative predictions about stranding patterns along oceanfront beaches. We use results from oceanic drift-bottle experiments to check our predictions and provide an upper limit on stranding proportions. We compare predicted current regimes from a 3D physical oceanographic model to spatial and temporal locations of both sea turtle carcass strandings and drift bottle landfalls. Drift bottle return rates suggest an upper limit for the proportion of sea turtle carcasses that strand (about 20%). In the South Atlantic Bight, seasonal development of along-shelf flow coincides with increased numbers of strandings of both turtles and drift bottles in late spring and early summer. The model also predicts net offshore flow of surface waters during winter - the season with the fewest relative strandings. The drift bottle data provide a reasonable upper bound on how likely carcasses are to reach land from points offshore and bound the general timeframe for stranding post-mortem (< two weeks). Our findings suggest that marine turtle strandings follow a seasonal regime predictable from physical oceanography and mimicked by drift bottle experiments. Managers can use these findings to reevaluate incidental strandings limits and fishery takes for both nearshore and offshore mortality sources. ?? 2005 Elsevier Ltd. All rights reserved.

Biological Conservation

Empirical parameterization of setup, swash, and runup

Using shoreline water-level time series collected during 10 dynamically diverse field experiments, an empirical parameterization for extreme runup, defined by the 2% exceedence value, has been developed for use on natural beaches over a wide range of conditions. Runup, the height of discrete water-level maxima, depends on two dynamically different processes; time-averaged wave setup and total swash excursion, each of which is parameterized separately. Setup at the shoreline was best parameterized using a dimensional form of the more common Iribarren-based setup expression that includes foreshore beach slope, offshore wave height, and deep-water wavelength. Significant swash can be decomposed into the incident and infragravity frequency bands. Incident swash is also best parameterized using a dimensional form of the Iribarren-based expression. Infragravity swash is best modeled dimensionally using offshore wave height and wavelength and shows no statistically significant linear dependence on either foreshore or surf-zone slope. On infragravity-dominated dissipative beaches, the magnitudes of both setup and swash, modeling both incident and infragravity frequency components together, are dependent only on offshore wave height and wavelength. Statistics of predicted runup averaged over all sites indicate a - 17 cm bias and an rms error of 38 cm: the mean observed runup elevation for all experiments was 144 cm. On intermediate and reflective beaches with complex foreshore topography, the use of an alongshore-averaged beach slope in practical applications of the runup parameterization may result in a relative runup error equal to 51% of the fractional variability between the measured and the averaged slope.

Coastal Engineering

Prediction and assimilation of surf-zone processes using a Bayesian network: Part II: Inverse models

A Bayesian network model has been developed to simulate a relatively simple problem of wave propagation in the surf zone (detailed in Part I). Here, we demonstrate that this Bayesian model can provide both inverse modeling and data-assimilation solutions for predicting offshore wave heights and depth estimates given limited wave-height and depth information from an onshore location. The inverse method is extended to allow data assimilation using observational inputs that are not compatible with deterministic solutions of the problem. These inputs include sand bar positions (instead of bathymetry) and estimates of the intensity of wave breaking (instead of wave-height observations). Our results indicate that wave breaking information is essential to reduce prediction errors. In many practical situations, this information could be provided from a shore-based observer or from remote-sensing systems. We show that various combinations of the assimilated inputs significantly reduce the uncertainty in the estimates of water depths and wave heights in the model domain. Application of the Bayesian network model to new field data demonstrated significant predictive skill (R 2 = 0.7) for the inverse estimate of a month-long time series of offshore wave heights. The Bayesian inverse results include uncertainty estimates that were shown to be most accurate when given uncertainty in the inputs (e.g., depth and tuning parameters). Furthermore, the inverse modeling was extended to directly estimate tuning parameters associated with the underlying wave-process model. The inverse estimates of the model parameters not only showed an offshore wave height dependence consistent with results of previous studies but the uncertainty estimates of the tuning parameters also explain previously reported variations in the model parameters.

Coastal Engineering

Currents, waves and sediment transport around the headland of Pt. Dume, California

Sediment transport past rocky headlands has received less attention compared to transport along beaches. Here we explore, in a field-based study, possible pathways for sediment movement adjacent to Point Dume, a headland in Santa Monica Bay, California. This prominent shoreline feature is a nearly symmetrical, triangular-shaped promontory interior to the Santa Monica Littoral Cell. We collected current, wave, and turbidity data for 74 days during which several wave events occurred, including one associated with a remote hurricane and another generated by the first winter storm of 2014. We also acquired sediment samples to quantify seabed grain-size distributions. Near-bottom currents towards the headland dominated on both of its sides and wave-driven longshore currents in the surf zone were faster on the exposed side. Bed shear stresseswere generated mostly by waves with minor contributions from currents, but both wave-driven and other currents contributed to sediment flux. On the wave-exposed west side of the headland, suspended sediment concentrations correlated with bed stress suggesting local resuspension whereas turbidity levels on the sheltered east side of the headland are more easily explained by advective delivery. Most of the suspended sediment appears to be exported offshore due to flow separation at the apex of the headland but may not move far given that sediment fluxes at moorings offshore of the apex were small. Further, wave-driven sediment flux in the surf zone is unlikely to pass the headland due to the discontinuity in wave forcing that causes longshore transport in different directions on each side of the headland. It is thus unlikely that sand is transported past the headland (specifically in a westerly direction), although some transport of finer fractions may occur offshore in deep water. These findings of minimal sediment flux past Point Dume are consistent with its role as a littoral cell boundary, although more complex multi-stage processes and unusual events may account for some transport at times.

California

Seismic displacement of gently-sloping coastal and marine sediment under multidirectional earthquake loading

Gentle sediment-laden slopes are typical of the onshore coastal zone and offshore continental shelf and slope. Coastal sediment are commonly young weakly consolidated materials that are well stratified, have low strength, and can mobilize shear displacements at low levels of stress. Seismically-driven plastic displacements of these sediment pose a hazard to coastal cities, buried onshore utilities, and offshore infrastructure like harbor protection and outfalls. One-dimensional rigid downslope-directed Newmark sliding block analyses have been used to predict earthquake deformations generally on steeper slopes that are modeled as frictional materials. This study probes the effect of multidirectional earthquake motions on inertial displacements of gently sloping ground of the coastal and offshore condition where soft-compliant soil is expected. Toward that objective, this investigation seeks to understand the effect on Newmark-type displacements of [1] multidirectional earthquake shaking and [2] soil compliance. In order to model multidirectional effects, the earthquake motions are rotated into the local slope strike- and dip-components. On gently sloping ground, including the strike component of motion always results in a larger and more accurate shear stress vector. Strike motions are found to contribute to downslope deformations on any declivity. Compliant response of the soil mass also influences the plastic displacements. The magnitude of seismic displacements can be estimated with a simplified model using only the estimated soil yield-acceleration ( k y ) and the peak ground velocity ( V max ) of the earthquake motions. Compliance effects can be effectively mapped using the concept of Plastic Displacement Response Spectra (PDRS).

Engineering Geology

Influence of sea level rise on iron diagenesis in an east Florida subterranean estuary

Subterranean estuary occupies the transition zone between hypoxic fresh groundwater and oxic seawater, and between terrestrial and marine sediment deposits. Consequently, we hypothesize, in a subterranean estuary, biogeochemical reactions of Fe respond to submarine groundwater discharge (SGD) and sea level rise. Porewater and sediment samples were collected across a 30-m wide freshwater discharge zone of the Indian River Lagoon (Florida, USA) subterranean estuary, and at a site 250 m offshore. Porewater Fe concentrations range from 0.5 μM at the shoreline and 250 m offshore to about 286 μM at the freshwater–saltwater boundary. Sediment sulfur and porewater sulfide maxima occur in near-surface OC-rich black sediments of marine origin, and dissolved Fe maxima occur in underlying OC-poor orange sediments of terrestrial origin. Freshwater SGD flow rates decrease offshore from around 1 to 0.1 cm/day, while bioirrigation exchange deepens with distance from about 10 cm at the shoreline to about 40 cm at the freshwater–saltwater boundary. DOC concentrations increase from around 75 μM at the shoreline to as much as 700 μM at the freshwater–saltwater boundary as a result of labile marine carbon inputs from marine SGD. This labile DOC reduces Fe-oxides, which in conjunction with slow discharge of SGD at the boundary, allows dissolved Fe to accumulate. Upward advection of fresh SGD carries dissolved Fe from the Fe-oxide reduction zone to the sulfate reduction zone, where dissolved Fe precipitates as Fe-sulfides. Saturation models of Fe-sulfides indicate some fractions of these Fe-sulfides get dissolved near the sediment–water interface, where bioirrigation exchanges oxic surface water. The estimated dissolved Fe flux is approximately 0.84 μM Fe/day per meter of shoreline to lagoon surface waters. Accelerated sea level rise predictions are thus likely to increase the Fe flux to surface waters and local primary productivity, particularly along coastlines where groundwater discharges through sediments.

Geochimica et Cosmochimica Acta

Evaluating the importance of abiotic and biotic drivers on Bythotrephes biomass in Lakes Superior and Michigan

The ability of planktivorous fishes to exert top-down control on Bythotrephes potentially has far-reaching impacts on aquatic food-webs, given previously described effects of Bythotrephes on zooplankton communities. We estimated consumption of Bythotrephes by planktivorous and benthivorous fishes, using bioenergetics and daily ration models at nearshore (18 m), intermediate (46 m), and offshore (110 m) depths along one western Lake Superior transect (April, and September-November) and two northern Lake Michigan transects (April, July, September). In Lake Superior, consumption (primarily by cisco Coregonus artedi) exceeded Bythotrephes production at all offshore sites in September-November (up to 396% of production consumed) and at the intermediate site in November (842%) with no evidence of consumption nearshore. By comparing Bythotrephes biomass following months of excessive consumption, we conservatively concluded that top-down control was evident only at the offshore site during September-October. In Lake Michigan, consumption by fishes (primarily alewife Alosa pseudoharengus) exceeded production at nearshore sites (up to 178%), but not in deeper sites (< 15%). Evidence for top-down control in the nearshore was not supported, however, as Bythotrephes never subsequently declined. Using generalized additive models, temperature, and not fish consumption, not zooplankton prey density, best explained variability in Bythotrephes biomass. The non-linear pattern revealed Bythotrephes to increase with temperature up to 16 &deg;C, and then decline between 16 and 23 &deg;C. We discuss how temperature likely has direct negative impacts on Bythotrephes when temperatures near 23 &deg;C, but speculate that predation also contributes to declining biomass when temperatures exceed 16 &deg;C.

Wisconsin and Michigan