Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Ocean Modeling”

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 685 records · Page 38Linked to original sources

Maximum entropy derived statistics of sound speed structure in a fine-grained sediment inferred from sparse broadband acoustic measurements on the New England continental shelf

Marginal probability distributions for parameters representing an effective sound-speed structure of a fine-grained sediment are inferred from a data ensemble maximum entropy method that utilizes a sparse spatially distributed set of received pressure time series resulting from multiple explosive sources in a shallow-water ocean environment possessing significant spatial variability of the seabed. A remote sensing seabed acoustics experiment undertaken in March 2017 off the New England Shelf was designed so that multiple independent analyses could infer the statistical properties of the seabed. The current analysis incorporates the measured horizontal variability from interpretations of a subbottom profiling survey of the experimental area. An idealized range- and azimuth-dependent parameterization of the seabed is derived from identification of horizons within the seabed that define multiple sediment layers. A sparse set of explosive charges were deployed on circular tracks with radii of about 2, 4, and 6.5 km with an acoustic array at the center to correlate a set of random measurements to physical acoustic processes that characterize the seabed. The mean values of a surface sound speed ratio and a linear sound speed gradient for the fine-grained sediment layer derived from 12 data samples processed in the 25–275-Hz band provide an estimate of the effective sound-speed structure in a 130-km $^2$ area. The inferred sediment sound speed values are evaluated by predicting measured time series data not used in the statistical inference, and are also compared to historical measurements. Finally, the low-frequency maximum entropy estimate of the sediment sound speed along with physical measurements derived from piston core measurements are utilized to estimate the sediment grain bulk modulus. predictions made by the viscous grain shearing model.

Maine, Vermont, New Hampshire, Massachusetts, Rhod↗

Characterizing operational signatures of reservoirs with the SWOT satellite by comparing natural lake and reservoir dynamics

Due to a lack of management operations data, hydrological models may represent reservoirs as natural lakes, leading to poor discharge predictions in regulated basins. To parse seasonal operational signatures, we compare the dynamics of natural lake and reservoir systems across North America using Surface Water and Ocean Topography (SWOT) satellite observations and derived discharge estimates. Overall, reservoirs and their adjacent river reaches exhibit significantly greater variability (in standard deviation) than their natural counterparts across almost all SWOT observed (e.g. water surface elevation) and inferred (e.g. discharge) variables. Natural lakes show strong same-day correlations between inflow and outflow discharge (median Spearman R = 0.8), whereas 76% of reservoirs exhibit maximum correlation when outflow is lagged, suggesting operations buffer seasonal flow variability. Our findings indicate operations not only affect reservoir dynamics themselves but also have upstream and downstream consequences, which, when integrated into models, will offer more realistic hydrologic conditions.

Environmental Research Letters↗

Dynamics of marsh-derived sediments in lagoon-type estuaries

Salt marshes are valuable ecosystems that must trap sediments and accrete in order to counteract the deleterious effect of sea‐level rise. Previous studies have shown that the capacity of marshes to build up vertically depends on both autogenous and exogenous processes including eco‐geomorphic feedbacks and sediment supply from in‐land and coastal ocean. There have been numerous efforts to quantify the role played by the sediments coming from marsh edge erosion on the resistance of salt marshes to sea‐level rise. However, the majority of existing studies investigating the interplay between lateral and vertical dynamics use simplified modelling approaches and they do not consider that marsh retreat can affect the regional scale hydrodynamics and sediment retention in back‐barrier basins. In this study, we evaluated the fate of the sediments originating from marsh lateral loss by using high‐resolution numerical model simulations of Jamaica Bay, a small lagoonal estuary located in New York City. Our findings show that up to 42% of the sediments released during marsh edge erosion deposits on the shallow areas of the basin and over the vegetated marsh platforms, contributing positively to the sediment budget of the remaining salt marshes. Furthermore, we demonstrate that with the present‐day sediment supply from the ocean the system cannot keep pace with sea‐level rise even accounting for the sediment liberated in the bay through marsh degradation. Our study highlights the relevance of multiple sediment sources for the maintenance of the marsh complex.

New York↗

Influence of net freshwater supply on salinity in Florida Bay

An annual water budget for Florida Bay, the large, seasonally hypersaline estuary in the Everglades National Park, was constructed using physically based models and long‐term (31 years) data on salinity, hydrology, and climate. Effects of seasonal and interannual variations of the net freshwater supply (runoff plus rainfall minus evaporation) on salinity variation within the bay were also examined. Particular attention was paid to the effects of runoff, which are the focus of ambitious plans to restore and conserve the Florida Bay ecosystem. From 1965 to 1995 the annual runoff from the Everglades into the bay was less than one tenth of the annual direct rainfall onto the bay, while estimated annual evaporation slightly exceeded annual rainfall. The average net freshwater supply to the bay over a year was thus approximately zero, and interannual variations in salinity appeared to be affected primarily by interannual fluctuations in rainfall. At the annual scale, runoff apparently had little effect on the bay as a whole during this period. On a seasonal basis, variations in rainfall, evaporation, and runoff were not in phase, and the net freshwater supply to the bay varied between positive and negative values, contributing to a strong seasonal pattern in salinity, especially in regions of the bay relatively isolated from exchanges with the Gulf of Mexico and Atlantic Ocean. Changes in runoff could have a greater effect on salinity in the bay if the seasonal patterns of rainfall and evaporation and the timing of the runoff are considered. One model was also used to simulate spatial and temporal patterns of salinity responses expected to result from changes in net freshwater supply. Simulations in which runoff was increased by a factor of 2 (but with no change in spatial pattern) indicated that increased runoff will lower salinity values in eastern Florida Bay, increase the variability of salinity in the South Region, but have little effect on salinity in the Central and West Regions.

Florida↗

Utilizing chromophoric dissolved organic matter measurements to derive export and reactivity of dissolved organic carbon exported to the Arctic Ocean: A case study of the Yukon River, Alaska

[1] The quality and quantity of dissolved organic matter (DOM) exported by Arctic rivers is known to vary with hydrology and this exported material plays a fundamental role in the biogeochemical cycling of carbon at high latitudes. We highlight the potential of optical measurements to examine DOM quality across the hydrograph in Arctic rivers. Furthermore, we establish chromophoric DOM (CDOM) relationships to dissolved organic carbon (DOC) and lignin phenols in the Yukon River and model DOC and lignin loads from CDOM measurements, the former in excellent agreement with long‐term DOC monitoring data. Intensive sampling across the historically under‐sampled spring flush period highlights the importance of this time for total export of DOC and particularly lignin. Calculated riverine DOC loads to the Arctic Ocean show an increase from previous estimates, especially when new higher discharge data are incorporated. Increased DOC loads indicate decreased residence times for terrigenous DOM in the Arctic Ocean with important implications for the reactivity and export of this material to the Atlantic Ocean.

Geophysical Research Letters↗

Quantifying the influence of sea ice on ocean microseism using observations from the Bering Sea, Alaska

Microseism is potentially affected by all processes that alter ocean wave heights. Because strong sea ice prevents large ocean waves from forming, sea ice can therefore significantly affect microseism amplitudes. Here we show that this link between sea ice and microseism is not only a robust one but can be quantified. In particular, we show that 75–90% of the variability in microseism power in the Bering Sea can be predicted using a fairly crude model of microseism damping by sea ice. The success of this simple parameterization suggests that an even stronger link can be established between the mechanical strength of sea ice and microseism power, and that microseism can eventually be used to monitor the strength of sea ice, a quantity that is not as easily observed through other means.

Alaska↗

Tectonic evolution of Gorda Ridge inferred from sidescan sonar images

Gorda Ridge is the southern segment of the Juan de Fuca Ridge complex, in the north-east Pacific. Along-strike spreading-rate variation on Gorda Ridge and deformation of Gorda Plate are evidence for compression between the Pacific and Gorda Plates. GLORIA sidescan sonographs allow the spreading fabric associated with Gorda Ridge to be mapped in detail. Between 5 and 2 Ma, a pair of propagating rifts re-orientated the northern segment of Gorda Ridge by about 10?? clockwise, accommodating a clockwise shift in Pacific-Juan de Fuca plate motion that occurred around 5 Ma. Deformation of Gorda Plate, associated with southward decreasing spreading rates along southern Gorda Ridge, is accommodated by a combination of clockwise rotation of Gorda Plate crust, coupled with left-lateral motion on the original normal faults of the ocean crust. Segments of Gorda Plate which have rotated by different amounts are separated by narrow deformation zones across which sharp changes in ocean fabric trend are seen. Although minor lateral movement may occur on these NW to WNW structures, no major right-lateral movement, as predicted by previous models, is observed. ?? 1988 Kluwer Academic Publishers.

Marine Geophysical Research↗

Future coral reef habitat marginality: Temporal and spatial effects of climate change in the Pacific basin

Marginal reef habitats are regarded as regions where coral reefs and coral communities reflect the effects of steady-state or long-term average environmental limitations. We used classifications based on this concept with predicted time-variant conditions of future climate to develop a scenario for the evolution of future marginality. Model results based on a conservative scenario of atmospheric CO2 increase were used to examine changes in sea surface temperature and aragonite saturation state over the Pacific Ocean basin until 2069. Results of the projections indicated that essentially all reef locations are likely to become marginal with respect to aragonite saturation state. Significant areas, including some with the highest biodiversity, are expected to experience high-temperature regimes that may be marginal, and additional areas will enter the borderline high temperature range that have experienced significant ENSO-related bleaching in the recent past. The positive effects of warming in areas that are presently marginal in terms of low temperature were limited. Conditions of the late 21st century do not lie outside the ranges in which present-day marginal reef systems occur. Adaptive and acclimative capabilities of organisms and communities will be critical in determining the future of coral reef ecosystems.

Conference Paper↗

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↗

Mysterious tsunami in the Caribbean Sea following the 2010 Haiti earthquake possibly generated by dynamically triggered early aftershocks

Dynamically triggered offshore aftershocks, caused by passing seismic waves from main shocks located on land, are currently not considered in tsunami warnings. The M7.0 2010 Haiti earthquake epicenter was located on land 27 km north of the Caribbean Sea and its focal mechanism was oblique strike-slip. Nevertheless, a tsunami recorded on a Caribbean Deep-Ocean Assessment and Reporting of Tsunami (DART) buoy and a tide gauge produced runup heights of 1–3 m along Haiti southeast coast. Earthquake finite-fault model inversions of the DART waveform suggest that a reverse fault doublet with magnitudes of M6.8 and M6.5 located 85 km southwest of the epicenter may have excited the tsunami. This doublet collocates with dynamically triggered aftershocks, derived from back-projection analysis, that occurred 20-60 s after the main shock of the Haiti earthquake. The aftershocks are within a region of maximum dynamic strain predicted by the main shock, on a possibly tectonically active submarine ridge southwest of Haiti's Southern Peninsula. The agreement between the tsunami finite-fault source models and the seismic and tectonic evidence suggests that earthquakes on land, even strike-slip faults, can generate tsunamis by dynamically triggering offshore aftershocks.

Earth and Planetary Science Letters↗

Generation and propagation of nonlinear internal waves in Massachusetts Bay

During the summer, nonlinear internal waves (NLIWs) are commonly observed propagating in Massachusetts Bay. The topography of the area is unique in the sense that the generation area (over Stellwagen Bank) is only 25 km away from the shoaling area, and thus it represents an excellent natural laboratory to study the life cycle of NLIWs. To assist in the interpretation of the data collected during the 1998 Massachusetts Bay Internal Wave Experiment (MBIWE98), a fully nonlinear and nonhydrostatic model covering the generation/shoaling region was developed, to investigate the response of the system to the range of background and driving conditions observed. Simplified models were also used to elucidate the role of nonlinearity and dispersion in shaping the NLIW field. This paper concentrates on the generation process and the subsequent evolution in the basin. The model was found to reproduce well the range of propagation characteristics observed (arrival time, propagation speed, amplitude), and provided a coherent framework to interpret the observations. Comparison with a fully nonlinear hydrostatic model shows that during the generation and initial evolution of the waves as they move away from Stellwagen Bank, dispersive effects play a negligible role. Thus the problem can be well understood considering the geometry of the characteristics along which the Riemann invariants of the hydrostatic problem propagate. Dispersion plays a role only during the evolution of the undular bore in the middle of Stellwagen Basin. The consequences for modeling NLIWs within hydrostatic models are briefly discussed at the end.

Massachusetts↗

One-dimensional wave bottom boundary layer model comparison: Specific eddy viscosity and turbulence closure models

Six one-dimensional-vertical wave bottom boundary layer models are analyzed based on different methods for estimating the turbulent eddy viscosity: Laminar, linear, parabolic, k —one equation turbulence closure, k−ε —two equation turbulence closure, and k−ω —two equation turbulence closure. Resultant velocity profiles, bed shear stresses, and turbulent kinetic energy are compared to laboratory data of oscillatory flow over smooth and rough beds. Bed shear stress estimates for the smooth bed case were most closely predicted by the k−ω model. Normalized errors between model predictions and measurements of velocity profiles over the entire computational domain collected at 15° intervals for one-half a wave cycle show that overall the linear model was most accurate. The least accurate were the laminar and k−ε models. Normalized errors between model predictions and turbulence kinetic energy profiles showed that the k−ω model was most accurate. Based on these findings, when the smallest overall velocity profile prediction error is required, the processing requirements and error analysis suggest that the linear eddy viscosity model is adequate. However, if accurate estimates of bed shear stress and TKE are required then, of the models tested, the k−ω model should be used.

Journal of Waterway, Port, Coastal and Ocean Engin↗

Simulated hydrologic responses to proposed wastewater-returnflow scenarios in Falmouth, Massachusetts

The Cape Cod aquifer is the sole source of drinking water for communities on Cape Cod, Massachusetts, including the Town of Falmouth, where the aquifer is currently threatened by contamination from septic-system-derived nitrogen. To address this problem, the Town is proposing to sewer areas of Falmouth, treat the wastewater at the Town’s Main Wastewater Treatment Facility (a nitrogen removing/tertiary treatment facility), and discharge the treated wastewater to an ocean outfall pipe in Nantucket Sound. The U.S. Geological Survey, in cooperation with the Town of Falmouth, updated a three-dimensional steady-state groundwater flow model to represent current (defined as 2019–23) average hydrologic conditions and to simulate the long-term average freshwater hydrologic response to two wastewater-return-flow scenarios. Scenario 1 involves the sewering of all properties south of Route 28 in Falmouth, which approximates the Town’s possible sewer expansion over the next 20–30 years. Scenario 2 involves sewering of all properties in Falmouth to demonstrate the maximum potential effect of sewering on the aquifer. Overall, the simulated hydrologic response of water-table altitudes and streamflow in both scenarios was relatively small compared to fluctuations from natural recharge. In scenario 1, the water-table altitude decreased by about 0.1 feet south of Route 28, where the conversion to municipal sewers removed wastewater-return flow from onsite septic systems. The water-table altitude decreased by about 0.1–0.2 feet over a larger area in Falmouth under town-wide sewering in scenario 2. The greatest decrease in water-table altitude in both scenarios occurred near the Main Wastewater Treatment Facility, with a decrease of about 1.1 feet in scenario 1 and about 1.3 feet in scenario 2. Simulated decreases in streamflow also were estimated for six selected streams in Falmouth and Mashpee. In both scenarios, the largest simulated decreases in streamflow were at the Coonamessett River, which is the closest stream to the Main Wastewater Treatment Facility. In scenario 1, the average annual decrease in flow at the Coonamessett River was 0.1 cubic feet per second, a 1.1 percent decrease from current (2019–23) conditions. In scenario 2, streamflow at the Coonamessett River decreased by 0.6 cubic feet per second, a 5.4 percent decrease from current (2019–23) conditions.

Massachusetts↗

A delayed action oscillator shared by biennial, interannual, and decadal signals in the Pacific Basin

Biennial, interannual, and decadal signals in the Pacific basin are observed to share patterns and evolution in covarying sea surface temperature (SST), 18??C isotherm depth (Z18), zonal surface wind (ZSW), and wind stress curl (WSC) anomalies from 1955 to 1999. Each signal has warm SST anomalies propagating slowly eastward along the equator, generating westerly ZSW anomalies in their wake. These westerly ZSW anomalies produce cyclonic WSC anomalies off the equator which pump baroclinic Rossby waves in the western/central tropical North Pacific Ocean. These Rossby waves propagate westward, taking ???6, ???12, and ???36 months to reach the western boundary near ???7??N, ???12??N, and ???18??N on biennial, interannual, and decadal period scales, respectively. There, they reflect as equatorial coupled waves, propagating slowly eastward in covarying SST, Z18, and ZSW anomalies, taking ???6, ???12, and ???24 months to reach the central/eastern equatorial ocean. These equatorial coupled waves produce a delayed-negative feedback to the warm SST anomalies there. The decrease in Rossby wave phase speed with latitude, the increase in meridional scale of equatorial SST anomalies with period scale, and the associated increase in latitude of Rossby wave forcing are consistent with the delayed action oscillator (DAO) model used to explain El Nin??o. However, this is not true of the western-boundary reflection of Rossby waves into slow equatorial coupled waves. This requires modification of the extant DAO model. We construct a modified DAO model, demonstrating how the various mechanisms and the size and sources of their delays yield the resulting frequency of each signal.

Journal of Geophysical Research C: Oceans↗

Simulated effects of groundwater withdrawals from aquifers in Ocean County and vicinity, New Jersey

Rapid population growth since the 1930s in Ocean County and vicinity, New Jersey, has placed increasing demands upon the area’s freshwater resources. To examine effects of groundwater withdrawals, a three-dimensional groundwater-flow model was developed to simulate the groundwater-flow systems of five area aquifers: the unconfined Kirkwood-Cohansey aquifer system and Vincentown aquifer, and three confined aquifers— the Rio Grande water-bearing zone, the Atlantic City 800-foot sand, and the Piney Point aquifer. The influence of withdrawals is evaluated by using transient groundwater-flow model simulations that incorporate three withdrawal schemes. These are (1) no-withdrawal conditions; (2) 2000–03 withdrawal conditions, using reported monthly withdrawals at all production wells from January 2000 through December 2003; and (3) maximum-allocation withdrawal conditions using the maximum withdrawal allowed by New Jersey Department of Environmental Protection permits at each well. Particle tracking analysis, using results from model simulations, delineated particle flow paths from production wells to the point of recharge, and estimated particle travel times. Compared with no-withdrawal conditions, 2000–03 withdrawal conditions reduced the amount of groundwater flow out of the Kirkwood-Cohansey aquifer system into streams, increased the net flow of water into other layers, reduced net flow into or out of storage, and reduced flow from the Kirkwood-Cohansey aquifer system to constant head cells. Freshwater discharging to the Barnegat Bay-Little Egg Harbor estuary from streams and groundwater is essential to maintaining the ecology of the bay. Examination of selected stress periods indicates that simulated base flow in streams flowing into the Barnegat Bay-Little Egg Harbor estuary is reduced by as much as 49 cubic feet per second for 2000 to 2003 withdrawal conditions when compared with no-withdrawal conditions. In the three confined aquifers, water levels during periods of low recharge and high withdrawals, and high recharge and low withdrawals, were examined to determine seasonal effects on the confined flow systems. The simulated potentiometric surface of the Rio Grande water-bearing zone and the Atlantic City 800-foot sand during selected stress periods indicates substantial declines from no-withdrawal conditions to 2000–03 conditions as a result of groundwater withdrawals. Cones of depression in Toms River Township, Seaside Heights and Seaside Park Boroughs, and Barnegat Light Borough developed in the potentiometric surface of the Piney Point aquifer in response to withdrawals. Maximum-allocation withdrawals decreased flow out of the Kirkwood-Cohansey aquifer system to constant head cells, increased flow out of the aquifer system to adjacent and lower layers, and reduced groundwater discharge to streams when compared with 2000–03 withdrawal conditions. Increases in withdrawals from the Rio Grande water-bearing zone, the Atlantic City 800-foot sand, and the Piney Point aquifer result in an increase in simulated net groundwater flow into these aquifers. Base-flow reduction from 2000–03 conditions to maximum-allocation conditions of 25 to 29 cubic feet per second in all streams draining to the Barnegat Bay-Little Egg Harbor also is indicated. Potentiometric surfaces of the Rio Grande water-bearing zone, Atlantic City 800-foot sand, and the Piney Point aquifer during two stress periods of simulated maximum-allocation withdrawal conditions indicated the expansion of several cones of depression developed during 2000–03 withdrawals. Simulation of average 2000–03 withdrawal conditions indicated the extent to which the groundwater-flow system is susceptible to potential saltwater intrusion into near-shore wells. Travel time from recharge to discharge location ranged from 11 to approximately 50,700 years in near-shore Kirkwood-Cohansey aquifer system wells. Those in Seaside Heights Borough, in Island Beach State Park (Berkeley Township), and in Ship Bottom Borough have particle travel times from 140 to 12,000 years and flow paths that originated under Barnegat Bay or the Atlantic Ocean from the simulation of average maximum-allocation withdrawal conditions. Travel time along flow paths to wells screened in the Rio Grande water-bearing zone and the Atlantic City 800-foot sand from recharge to discharge point ranged from nearly 530 years to greater than 3.73 million years from the simulation of average 2000–03 withdrawal conditions. Particle tracking indicated that most wells screened in these aquifers derived a large part of their recharge from the Oswego River Basin, with a small portion of flow originating either beneath Barnegat Bay or to the east beneath the Atlantic Ocean. Travel time along flow paths that start beneath either Barnegat Bay or the Atlantic Ocean ranged from 2,300 to approximately 134,000 years from the simulation of average maximum-allocation withdrawal conditions."

New Jersey↗

Geoelectrical structure of the central zone of Piton de la Fournaise volcano (Reunion)

A study of the geoelectrical structure of the central part of Piton de la Fournaise volcano (Reunion, Indian Ocean) was made using direct current electrical (DC) and transient electromagnetic soundings (TEM). Piton de la Fournaise is a highly active oceanic basaltic shield and has been active for more than half a million years. Joint interpretation of the DC and TEM data allows us to obtain reliable 1D models of the resistivity distribution. The depth of investigation is of the order of 1.5 km but varies with the resistivity pattern encountered at each sounding. Two-dimensional resistivity cross sections were constructed by interpolation between the soundings of the 1D interpreted models. Conductors with resistivities less than 100 ohm-m are present at depth beneath all of the soundings and are located high in the volcanic edifice at elevations between 2000 and 1200 m. The deepest conductor has a resistivity less than 20 ohm-m for soundings located inside the Enclos and less than 60-100 ohm-m for soundings outside the Enclos. From the resistivity distributions, two zones are distinguished: (a) the central zone of the Enclos; and (b) the outer zone beyond the Enclos. Beneath the highly active summit area, the conductor rises to within a few hundred meters of the surface. This bulge coincides with a 2000-mV self-potential anomaly. Low-resistivity zones are inferred to show the presence of a hydrothermal system where alteration by steam and hot water has lowered the resistivity of the rocks. Farther from the summit, but inside the Enclos the depth to the conductive layers increases to approximately 1 km and is inferred to be a deepening of the hydrothermally altered zone. Outside of the Enclos, the nature of the deep, conductive layers is not established. The observed resistivities suggest the presence of hydrated minerals, which could be found in landslide breccias, in hydrothermally altered zones, or in thick pyroclastic layers. Such formations often create perched water tables. The known occurrence of large eastward-moving landslides in the evolution of Piton de la Fournaise strongly suggests that large volumes of breccias should exist in the interior of the volcano; however, extensive breccia deposits are not observed at the bottom of the deep valleys that incise the volcano to elevations lower than those determined for the top of the conductors. The presence of the center of Piton de la Fournaise beneath the Plaine des Sables area during earlier volcanic stages (ca. 0.5 to 0.150 Ma) may have resulted in broad hydrothermal alteration of this zone. However, this interpretation cannot account for the low resistivities in peripheral zones. It is not presently possible to discriminate between these general interpretations. In addition, the nature of the deep conductors may be different in each zone. Whatever the geologic nature of these conductive layers, their presence indicates a major change of lithology at depth, unexpected for a shield volcano such as Piton de la Fournaise.

Bulletin of Volcanology↗

Reconciliation of geochronology and paleozoogeography for Quaternary marine terraces, San Luis Obispo Bay area, California, USA

In the San Luis Obispo Bay area of central California, interpretations of marine terrace ages have been hampered by inconsistent results from geochronological indicators (U-series ages of corals and correlations using amino acid racemization of mollusks) and seemingly contradictory paleozoogeographic aspects of fossil faunal assemblages. New U-series ages of corals, amino acid data from mollusks, and detailed analyses of fossil mollusk assemblages are presented that reconcile many of the apparently discordant results from previous studies. The two lowest-elevation terraces are Q1 (lower) and Q2 (upper); both are thought to date from high-sea stands of marine isotope stage (MIS) 5, the last interglacial complex. A combination of U-series dating and amino acid racemization results indicates that the Q1 terrace probably dates to the ∼80 ka high-sea stand of MIS 5a. U-series analyses of corals from the Q2 terrace show open-system histories, but consideration of two alternative open-system histories indicates that this terrace likely hosts corals dating to the high-sea stands of MIS 5c (∼105–95 ka) and MIS 5e (∼130–115 ka). Amino acid ratios in the fossil bivalve Leukoma staminea (Conrad, 1837) support the age differences between the two terraces and the open-system models suggested by the U-series data. Alkenone studies along much of the coast of California have shown that sea surface temperatures (SSTs) of the eastern Pacific Ocean were warmer than present during MIS 5e. During MIS 5c, on the other hand, SSTs were either similar to the present or cooler, but in either case, much cooler than during MIS 5e. The paleozoogeographic aspects of the molluscan faunas from the Q2 terrace support an interpretation that these deposits contain a mixture of fossils from both the MIS 5c and 5e high-sea stands. Fossils from the Q2 terrace include southern or southward-ranging species (likely dating to MIS 5e) and northern or northward-ranging species (likely dating to MIS 5c). A mixture of MIS 5e and MIS 5c fossils on the Q2 terrace is consistent with similar data reported for other localities in California. This phenomenon is explained by glacial isostatic adjustment (GIA) processes, which brought about a higher-than-present sea level during MIS 5c. GIA processes, combined with a low uplift rate, allowed mixing of MIS 5e fossils during the MIS 5c high-sea stand through capture of an MIS 5e terrace that had experienced minimal uplift in ∼20,000 years. Terrace reoccupation and fossil mixing such as that in the San Luis Obispo Bay area can be expected along any part of the Pacific Coast of North America that was affected by GIA processes and where uplift rates are low.

California↗

The effects of mixotrophy on the stability and dynamics of a simple planktonic food web

Recognition of the microbial loop as an important part of aquatic ecosystems disrupted the notion of simple linear food chains. However, current research suggests that even the microbial loop paradigm is a gross simplification of microbial interactions due to the presence of mixotrophs&mdash;organisms that both photosynthesize and graze. We present a simple food web model with four trophic species, three of them arranged in a food chain (nutrients&ndash;autotrophs&ndash;herbivores) and the fourth as a mixotroph with links to both the nutrients and the autotrophs. This model is used to study the general implications of inclusion of the mixotrophic link in microbial food webs and the specific predictions for a parameterization that describes open ocean mixed layer plankton dynamics. The analysis indicates that the system parameters reside in a region of the parameter space where the dynamics converge to a stable equilibrium rather than displaying periodic or chaotic solutions. However, convergence requires weeks to months, suggesting that the system would never reach equilibrium in the ocean due to alteration of the physical forcing regime. Most importantly, the mixotrophic grazing link seems to stabilize the system in this region of the parameter space, particularly when nutrient recycling feedback loops are included.

Theoretical Population Biology↗