Search USGSSearch

SEARCH · Search USGS

Results for “Ocean Modelling”

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 253 records · Page 14Linked to original sources

The restricted gemuk group: A triassic to lower cretaceous succession in southwestern Alaska

New data from an Upper Triassic to Lower Cretaceous deep marine succession-the herein reinstated and restricted Gemuk Group-provide a vital piece of the puzzle for unraveling southwestern Alaska's tectonic history. First defined by Cady et al. in 1955, the Gemuk Group soon became a regional catchall unit that ended up as part of at least four different terranes. In this paper we provide the first new data in nearly half a century from the Gemuk Group in the original type area in Taylor Mountains quadrangle and from contiguous rocks to the north in Sleetmute quadrangle. Discontinuous exposure, hints of complex structure, the reconnaissance level of our mapping, and spotty age constraints together permit definition of only a rough stratigraphy. The restricted Gemuk Group is at least 2250 m thick, and could easily be at least twice as thick. The age range of the restricted Gemuk Group is tightened on the basis of ten radiolarian ages, two new bivalve ages, one conodont age, two U-Pb zircon ages on tuff, and U-Pb ages of 110 detrital zircons from two sandstones. The Triassic part of the restricted Gemuk Group, which consists of intermediate pillow lavas interbedded with siltstone, chert, and rare limestone, produced radiolarians, bivalves, and conodonts of Carnian and Norian ages. The Jurassic part appears to be mostly siltstone and chert, and yielded radiolarians of Hettangian- Sinemurian, Pliensbachian-Toarcian, and Oxfordian ages. Two tuffs near the Jurassic-Cretaceous boundary record nearby arc volcanism: one at 146 Ma is interbedded with red and green siltstone, and a second at ca. 137 Ma is interbedded with graywacke turbidites. Graywacke appears to be the dominant rock type in the LowerCretaceous part of the restricted Gemuk Group. Detrital zircon analyses were performed on two sandstone samples using SHRIMP. One sandstone yielded a dominant age cluster of 133-180 Ma; the oldest grain is only 316 Ma. The second sample is dominated by zircons of 130-154 Ma; the oldest grain is 292 Ma. The youngest zircons are probably not much older than the sandstone itself. Point counts of restricted Gemuk Group sandstones yield average ratios of 24/29/47 for Q/F/L, 15/83/2 for Ls/Lv/Lm, and 41/48/11 for Qm/P/K. In the field, sandstones of the restricted Gemuk Group are not easily distinguished from sandstones of the overlying Upper Cretaceous turbidite-dominated Kuskokwim Group. Petrographically, however, the restricted Gemuk Group has modal K-feldspar, whereas the Kuskokwim Group generally does not (average Qm/P/K of 64/36/0). Some K-feldspar-bearing graywacke that was previously mapped as Kuskokwim Group (Cady et al., 1955) is here reassigned to the restricted Gemuk Group. Major- and trace element geochemistry of shales from the restricted Gemuk Group and the Kuskokwim Group show distinct differences. The chemical index of alteration (CIA) is distinctly higher forshales of the Kuskokwim Group than for those of the restricted Gemuk Group, suggesting more intense weathering during deposition of the Kuskokwim Group. The restricted Gemuk Group represents an estimated 90-100 m.y. of deep-water sedimentation, first accompanied by submarine volcanism and later by nearby explosive arc activity. Two hypotheses are presented for the tectonic setting. One model that needs additional testing is that the restricted Gemuk Group consists of imbricated oceanic plate stratigraphy. Based on available information, our preferred model is that it was deposited in a back-arc, intra-arc, or forearc basin that was subsequently deformed. The terrane affinity of the restricted Gemuk Group is uncertain. The rocks of this area were formerly assigned to the Hagemeister subterrane of the Togiak terrane-a Late Triassic to Early Cretaceous arc-but our data show this to be a poor match. None of the other possibilities (e.g., Nukluk and Tikchik subterranes of the Goodnews terrane) is viable; hence, the terrane subdivision and distribution in southwestern Alaska may need

Special Paper of the Geological Society of America

Eastern Aleutian volcanic arc digital model - version 1.0

A 3-dimensional model (Figure 1) of the interaction of oceanic and continental tectonic plates along the eastern portion of the Aleutian volcanic arc helps in the visualization of basic tectonic, geodetic, and geophysical data in this active plate boundary region. The model is constrained by topographic, bathymetric, and seismic data and by the principle of isostasy. Examination of free-air gravity anomalies over the region indicates where the flexural strength of the down-going oceanic slab disturbs local isostatic balance and where low-density sediments have accumulated in the trench and forearc regions.

Open-File Report

Joint spatiotemporal models to predict seabird densities at sea

Introduction: Seabirds are abundant, conspicuous members of marine ecosystems worldwide. Synthesis of distribution data compiled over time is required to address regional management issues and understand ecosystem change. Major challenges when estimating seabird densities at sea arise from variability in dispersion of the birds, sampling effort over time and space, and differences in bird detection rates associated with survey vessel type. Methods: Using a novel approach for modeling seabirds at sea, we applied joint dynamic species distribution models (JDSDM) with a vector-autoregressive spatiotemporal framework to survey data collected over nearly five decades and archived in the North Pacific Pelagic Seabird Database. We produced monthly gridded density predictions and abundance estimates for 8 species groups (77% of all birds observed) within Cook Inlet, Alaska. JDSDMs included habitat covariates to inform density predictions in unsampled areas and accounted for changes in observed densities due to differing survey methods and decadal-scale variation in ocean conditions. Results: The best fit model provided a high level of explanatory power (86% of deviance explained). Abundance estimates were reasonably precise, and consistent with limited historical studies. Modeled densities identified seasonal variability in abundance with peak numbers of all species groups in July or August. Seabirds were largely absent from the study region in either fall (e.g., murrelets) or spring (e.g., puffins) months, or both periods (shearwaters). Discussion: Our results indicated that pelagic shearwaters ( Ardenna spp.) and tufted puffin ( Fratercula cirrhata ) have declined over the past four decades and these taxa warrant further investigation into underlying mechanisms explaining these trends. JDSDMs provide a useful tool to estimate seabird distribution and seasonal trends that will facilitate risk assessments and planning in areas affected by human activities such as oil and gas development, shipping, and offshore wind and renewable energy.

Alaska

Tsunami modelling with adaptively refined finite volume methods

Numerical modelling of transoceanic tsunami propagation, together with the detailed modelling of inundation of small-scale coastal regions, poses a number of algorithmic challenges. The depth-averaged shallow water equations can be used to reduce this to a time-dependent problem in two space dimensions, but even so it is crucial to use adaptive mesh refinement in order to efficiently handle the vast differences in spatial scales. This must be done in a 'wellbalanced' manner that accurately captures very small perturbations to the steady state of the ocean at rest. Inundation can be modelled by allowing cells to dynamically change from dry to wet, but this must also be done carefully near refinement boundaries. We discuss these issues in the context of Riemann-solver-based finite volume methods for tsunami modelling. Several examples are presented using the GeoClaw software, and sample codes are available to accompany the paper. The techniques discussed also apply to a variety of other geophysical flows. ?? 2011 Cambridge University Press.

Acta Numerica

A geologic interpretation of seismic-refraction results in northeastern California

In 1981, the U.S. Geological Survey conducted a seismic-refraction experiment in northeastern California designed to study the Klamath Mountains, Cascade Range, Modoc Plateau, and Basin and Range provinces. Key profiles include 135-km-long, north-south lines in the Klamath Mountains and Modoc Plateau provinces and a 260-km-long, east-west line crossing all of the provinces. The seismic-velocity models for the Klamath and Modoc lines are comparatively homogeneous laterally but are quite different from each other. The Klamath model is finely layered from the surface to at least 14-km depth, consisting of a series of high-velocity layers (6.1–6.7 km/s), ranging in thickness from 1 to 4 km, with alternating positive and negative velocity gradients. A layer with an unreversed velocity of 7.0 km/s extends from 14 km to an unknown depth. The Modoc model, in contrast, is relatively thickly layered and has lower velocities than does the Klamath model at all depths down to 25 km. An upper layer, 4.5 km thick, of low-velocity material (2.1–4.4 km/s) overlies a basement with a considerably higher velocity (6.2 km/s). Velocity increases slowly with depth, with a small velocity step (to 6.4 km/s) at 11 km and a 7.0-km/s layer beginning at 25-km depth. Moho is probably 38–45 km deep under the Modoc Plateau, but its depth is unknown under the Klamath Mountains. A combined velocity-density model for the east-west line consists of a western part similar in configuration to the Klamath velocity model, an eastern part similar to the Modoc velocity model, and laterally changing velocity-density structure in between, in the Cascade Range. Beneath its upper layer, the velocity model for the Modoc Plateau is similar to that determined by other researchers for the adjacent Sierra Nevada. The velocity model is unlike those for rift areas, to which the Modoc Plateau has been compared by some authors. We theorize that beneath a veneer of volcanic and sedimentary rocks (the upper layer), the Modoc Plateau is underlain by a basement of granitic and metamorphic rocks that, like rocks in the Sierra Nevada, are the roots of one or more magmatic arcs. The fine layering in the Klamath seismic-velocity model is consistent with the geologic structure of the Klamath Mountains, characterized by imbricate thrusting of oceanic rock layers of various compositions and ages. Independent modeling of aeromagnetic data indicates that the base of the Trinity ultramafic sheet, the second major rock layer down in the structural sequence, corresponds to a velocity step to 6.7 km/s at 7-km depth in our model. The 6.7-km/s layer beneath the Trinity ultramafic sheet apparently corresponds to rocks of the central metamorphic belt, which are mafic schists. Rock units structurally deeper than rocks of the central metamorphic belt can be correlated with velocity layers below the 6.7-km/s layer, but with less certainty. In the model for the east-west line, the region of laterally changing velocity structure beneath the Cascade Range includes a 10-km step down to the east in the top of the 7.0-km/s layer. This region of lateral velocity change we interpret to be a fault, fold, or intrusive contact (or some combination of the three) between the stack of oceanic rock layers that underlie the Klamath Mountains and the buried roots of magmatic arcs inferred to underlie the Modoc Plateau. Magmas forming the modern Cascade Range arc apparently rise through this region.

GSA Bulletin

A simple approach to modeling light attenuation in the Sacramento-San Joaquin Delta using commonly available data

The diffuse attenuation coefficient of photosynthetically active radiation (KdPAR) is commonly used to predict light attenuation in aquatic productivity models, but obtaining measurements of PAR to compute KdPAR is difficult. In situ calculations of KdPAR require multiple measurements of PAR through the water column, and these measurements are infeasible for real-time recording. Instead, predictive models using surface-water measurements may be used. Traditional KdPAR models are based on open-ocean habitats and rely on chlorophyll—as a proxy measurement for phytoplankton abundance—as the main predictive parameter. However, elevated suspended sediments and dissolved organic materials may also affect KdPAR values of inland water bodies and estuaries. In this study, we leverage KdPAR calculations derived from in situ light measurements collected along with surface-water-quality parameters across the Sacramento-San Joaquin River Delta in California, USA (the Delta). Sampling occurred between January of 2013 and May of 2014. We also explored regional and seasonal effects, but these did not clearly affect the model. Ultimately, the best-performing model included surface-level turbidity only (R2 = 0.91). The simplicity of the model facilitates use of KdPAR estimates for a variety of purposes throughout the Delta, including euphotic depth calculations, and as inputs to primary-productivity and habitat-suitability models. We demonstrate the model’s usability with two open-sources data sets (one spatially dense, and one temporally dense), and estimate KdPAR, euphotic depth, and primary productivity within the Delta. We provide calculations for each estimation, allowing users to easily adopt these models and apply them to their own data or with open-sourced data, which are abundant.

California

Sedimentary processes in modern and ancient oceanic arc settings: evidence from the Jurassic Talkeetna Formation of Alaska and the Mariana and Tonga Arcs, western Pacific

Sediment deposited around oceanic volcanic ares potentially provides the most complete record of the tectonic and geochemical evolution of active margins. The use of such tectonic and geochemical records requires an accurate understanding of sedimentary dynamics in an arc setting: processes of deposition and reworking that affect the degree to which sediments represent the contemporaneous volcanism at the time of their deposition. We review evidence from the modern Mariana and Tonga arcs and the ancient arc crustal section in the Lower Jurassic Talkeetna Formation of south-central Alaska, and introduce new data from the Mariana Arc, to produce a conceptual model of volcaniclastic sedimentation processes in oceanic arc settings. All three arcs are interpreted to have formed in tectonically erosive margin settings, resulting in long-term extension and subsidence. Debris aprons composed of turbidites and debris flow deposits occur in the immediate vicinity of arc volcanoes, forming relatively continuous mass-wasted volcaniclastic records in abundant accommodation space. There is little erosion or reworking of old volcanic materials near the arc volcanic front. Tectonically generated topography in the forearc effectively blocks sediment flow from the volcanic front to the trench; although some canyons deliver sediment to the trench slope, most volcaniclastic sedimentation is limited to the area immediately around volcanic centers. Arc sedimentary sections in erosive plate margins can provide comprehensive records of volcanism and tectonism spanning < 10 My. The chemical evolution of a limited section of an oceanic arc may be best reconstructed from sediments of the debris aprons for intervals up to ~ 20 My but no longer, because subduction erosion causes migration of the forearc basin crust and its sedimentary cover toward the trench, where there is little volcaniclastic sedimentation and where older sediments are dissected and reworked along the trench slope.

Alaska

Ice ages and the thermal equilibrium of the earth

A model for climatic change, and particularly for the changes of the late Cenozoic, involves as its primary variables the albedo of the earth and the heat storage characteristics of the oceans. Geography exerts a strong influence. The model proposed does not require metaphysical causes or astronomical events other than known variations in the earth's orbit.

Journal of Research of the U.S. Geological Survey

Analysis of the Shallow Groundwater Flow System at Fire Island National Seashore, Suffolk County, New York

Fire Island National Seashore (FIIS) occupies 42 kilometers of the barrier island for which it is named that lies off the southern shore of Suffolk County, N.Y. Freshwater in the highly permeable, sandy aquifer underlying Fire Island is bounded laterally by marine surface waters and at depth by saline groundwater. Interspersed throughout FIIS are 17 pre-existing residential communities that in summer months greatly increase in population through the arrival of summer residents and vacationers; in addition, the National Park Service (NPS) has established several facilities on the island to accommodate visitors to FIIS. The 2.2 million people estimated by the NPS to visit Fire Island annually impact groundwater quality through the release of waste-derived contaminants, such as nutrients, pathogens, and organic compounds, into the environment. Waste-contaminated groundwater can move through the aquifer and threaten the ecological health of the adjacent back-barrier estuaries to which much of the groundwater ultimately discharges. In 2004, the U.S. Geological Survey (USGS), in cooperation with the NPS, began a 3-year investigation to (1) collect groundwater levels and water-quality (nutrient) samples, (2) develop a three-dimensional model of the shallow (water-table) aquifer system and adjacent marine surface waters, and (3) calculate nitrogen loads in simulated groundwater discharges from the aquifer to back-barrier estuaries and the ocean. The hydrogeology of the shallow aquifer system was characterized from the results of exploratory drilling, geophysical surveying, water-level monitoring, and water-quality sampling. The investigation focused on four areas-the communities of Kismet and Robbins Rest, the NPS Visitor Center at Watch Hill, and the undeveloped Otis Pike Fire Island High Dune Wilderness. Thirty-five observation wells were installed within FIIS to characterize subsurface hydrogeology and establish a water-table monitoring network in the four study areas. A variable-density model of the shallow aquifer system and adjacent marine surface waters was developed to simulate groundwater flow patterns and rates. Nitrogen loads from the shallow aquifer system were calculated from representative total nitrogen (TN) concentrations and simulated groundwater discharges to back-barrier estuaries and the ocean. The model simulates groundwater directions, velocities, and discharge rates under 2005 mean annual conditions. Groundwater budgets were developed for recharge areas of similar land use that contribute freshwater to back-barrier estuaries, the ocean, and subsea-discharge zones. Total freshwater discharge from the shallow aquifer system is about 43,500 cubic meters per day (m3/d) (79.8 percent) to back-barrier estuaries and about 10,200 m3/d (18.7 percent) to the ocean; about 836 m3/d (1.5 percent) may exit the system as subsea underflow. The total contribution of fresh groundwater to shoreline discharge zones amounts to about 53,700 m3/d (98.5 percent). The median age of freshwater discharged to back-barrier estuaries and the ocean was 3.4 years, and the 95th-percentile age was 20 years. The TN concentrations and loads under 2005 mean annual conditions for areas that contribute fresh groundwater to back-barrier estuaries and the ocean were calculated for the principal land uses on Fire Island. The overall TN load from the shallow aquifer system to shoreline discharge zones is about 16,200 kilograms per year (kg/yr) (82.2 percent) to back-barrier estuaries and about 3,500 kg/yr (17.8 percent) to the ocean. The overall TN load to marine surface waters amounts to about 19,700 kg/yr-roughly 6 percent of the annual TN load from shallow groundwater entering the South Shore Estuary Reserve (SSER) from the Suffolk County mainland, which is about 345,000 kg/yr. In contrast to the TN load from shallow groundwater for the SSER watershed, which annually yields about 353 kilograms per square kilometer (kg/km2), the overall TN loa

Scientific Investigations Report

Productivity and linkages of the food web of the southern region of the western Antarctic Peninsula continental shelf

The productivity and linkages in the food web of the southern region of the west Antarctic Peninsula continental shelf were investigated using a multi-trophic level mass balance model. Data collected during the Southern Ocean Global Ocean Ecosystem Dynamics field program were combined with data from the literature on the abundance and diet composition of zooplankton, fish, seabirds and marine mammals to calculate energy flows in the food web and to infer the overall food web structure at the annual level. Sensitivity analyses investigated the effects of variability in growth and biomass of Antarctic krill ( Euphausia superba ) and in the biomass of Antarctic krill predators on the structure and energy fluxes in the food web. Scenario simulations provided insights into the potential responses of the food web to a reduced contribution of large phytoplankton (diatom) production to total primary production, and to reduced consumption of primary production by Antarctic krill and mesozooplankton coincident with increased consumption by microzooplankton and salps. Model-derived estimates of primary production were 187–207 g C m −2 y −1 , which are consistent with observed values (47–351 g C m −2 y −1 ). Simulations showed that Antarctic krill provide the majority of energy needed to sustain seabird and marine mammal production, thereby exerting a bottom-up control on higher trophic level predators. Energy transfer to top predators via mesozooplanton was a less efficient pathway, and salps were a production loss pathway because little of the primary production they consumed was passed to higher trophic levels. Increased predominance of small phytoplankton (nanoflagellates and cryptophytes) reduced the production of Antarctic krill and of its predators, including seabirds and seals.

Progress in Oceanography

The vegetation outlook (VegOut): a new method for predicting vegetation seasonal greenness

The vegetation outlook (VegOut) is a geospatial tool for predicting general vegetation condition patterns across large areas. VegOut predicts a standardized seasonal greenness (SSG) measure, which represents a general indicator of relative vegetation health. VegOut predicts SSG values at multiple time steps (two to six weeks into the future) based on the analysis of "historical patterns" (i.e., patterns at each 1 km grid cell and time of the year) of satellite, climate, and oceanic data over an 18-year period (1989 to 2006). The model underlying VegOut capitalizes on historical climate-vegetation interactions and ocean-climate teleconnections (such as El Niño and the Southern Oscillation, ENSO) expressed over the 18-year data record and also considers several environmental characteristics (e.g., land use/cover type and soils) that influence vegetation's response to weather conditions to produce 1 km maps that depict future general vegetation conditions. VegOut provides regionallevel vegetation monitoring capabilities with local-scale information (e.g., county to sub-county level) that can complement more traditional remote sensing-based approaches that monitor "current" vegetation conditions. In this paper, the VegOut approach is discussed and a case study over the central United States for selected periods of the 2008 growing season is presented to demonstrate the potential of this new tool for assessing and predicting vegetation conditions.

GIScience and Remote Sensing

Regional estimates of reef carbonate dynamics and productivity Using Landsat 7 ETM+, and potential impacts from ocean acidification

Using imagery at 30 m spatial resolution from the most recent Landsat satellite, the Landsat 7 Enhanced Thematic Mapper Plus (ETM+), we scale up reef metabolic productivity and calcification from local habitat-scale (10 -1 to 10 0 km 2 ) measurements to regional scales (10 3 to 10 4 km 2 ). Distribution and spatial extent of the North Florida Reef Tract (NFRT) habitats come from supervised classification of the Landsat imagery within independent Landsat-derived Millennium Coral Reef Map geomorphologic classes. This system minimizes the depth range and variability of benthic habitat characteristics found in the area of supervised classification and limits misclassification. Classification of Landsat imagery into 5 biotopes (sand, dense live cover, sparse live cover, seagrass, and sparse seagrass) by geomorphologic class is >73% accurate at regional scales. Based on recently published habitat-scale in situ metabolic measurements, gross production (P = 3.01 ?? 10 9 kg C yr -1 ), excess production (E = -5.70 ?? 10 8 kg C yr -1 ), and calcification (G = -1.68 ?? 10 6 kg CaCO 3 yr -1 ) are estimated over 2711 km 2 of the NFRT. Simple models suggest sensitivity of these values to ocean acidification, which will increase local dissolution of carbonate sediments. Similar approaches could be applied over large areas with poorly constrained bathymetry or water column properties and minimal metabolic sampling. This tool has potential applications for modeling and monitoring large-scale environmental impacts on reef productivity, such as the influence of ocean acidification on coral reef environments. ?? Inter-Research 2009.

Marine Ecology Progress Series

Baseline monitoring of the western Arctic Ocean estimates 20% of the Canadian Basin surface waters are undersaturated with respect to aragonite

Marine surface waters are being acidified due to uptake of anthropogenic carbon dioxide, resulting in surface ocean areas of undersaturation with respect to carbonate minerals, including aragonite. In the Arctic Ocean, acidification is expected to occur at an accelerated rate with respect to the global oceans, but a paucity of baseline data has limited our understanding of the extent of Arctic undersaturation and of regional variations in rates and causes. The lack of data has also hindered refinement of models aimed at projecting future trends of ocean acidification. Here, based on more than 34,000 data records collected in 2010 and 2011, we establish a baseline of inorganic carbon data (pH, total alkalinity, dissolved inorganic carbon, partial pressure of carbon dioxide, and aragonite saturation index) for the western Arctic Ocean. This data set documents aragonite undersaturation in ~20% of the surface waters of the combined Canada and Makarov basins, an area characterized by recent acceleration of sea ice loss. Conservative tracer studies using stable oxygen isotopic data from 307 sites show that while the entire surface of this area receives abundant freshwater from meteoric sources, freshwater from sea ice melt is most closely linked to the areas of carbonate mineral undersaturation. These data link the Arctic Ocean’s largest area of aragonite undersaturation to sea ice melt and atmospheric CO 2 absorption in areas of low buffering capacity. Some relatively supersaturated areas can be linked to localized biological activity. Collectively, these observations can be used to project trends of ocean acidification in higher latitude marine surface waters where inorganic carbon chemistry is largely influenced by sea ice meltwater.

PLoS ONE

Ten years of studies on Maryland's inner Continental Margin and coastal bays

During the past ten years of the Association of American State Geologists-Mineral Management Service Continental Margins Program, the Maryland Geological Survey investigated the sedimentological, paleontological, stratigraphical and geophysical character of Maryland's inner continental shelf. Based on seismic records and sedimentological analyses completed during the first four years, a late Quaternary stratigraphic model was developed. Five distinct stratigraphic units were identified and described on the Maryland inner shelf. These units represent late Pleistocene interglacial deposits, the oldest of which corresponds to pre-Illinoian (oxygen-isotope stages 7 and/or 9) transgressive shelf sands. Overlying the Q1 unit, the Q2 unit is a 6-meter thick mud sequence of oxygen-isotope stage 5 (128-75 ka) age. Units Q3 and Q4 representing fluvial and leading edge estuarine deposits (oxygen-isotope stages 4, 3 and 2) filled numerous paleochannels that were incised into units Q2 and Q1. Modern trailing-edge transgressive shelf shoals (Unit Q5) discontinuously cap the sequence. The 5th and 6th years studies reported on the economic minerals of surficial and cored sediments. Vibracores collected off the Maryland's shelf during previous studies were analyzed for mineral types and abundances, weight percent of general size fractions, and heavy mineral (HM) content. Mineralogic maturity indices were compiled to correlate the THM and economic heavy minerals (EHM) abundances with position offshore, sediment type, and the indices themselves. For the 7th year, the Maryland Geological Survey re-examined geophysical records and lithological data originally collected by the Army Corps of Engineers to locate and assess beach fill borrow areas for the Ocean City Beach Replenishment Project. Data from 163 vibracores and over 300 kilometers of high-resolution seismic profile records collected off Ocean City, Maryland, supported the stratigraphic model developed by MGS during the first four years of the AASG-MMS program. The Maryland coastal bays became the focus of study during the 8th and 9th years during which seismic records, cores and surficial sediment were collected in Isle of Wight and Assawoman Bays. Shallow pretransgression surface was mapped, relating the existing streams to offshore paleochannels. The tenth year study focused on developing a repository for vibracores collected on Maryland's inner continental shelf.During the past ten years of the Association of American State Geologists-Mineral Management Service Continental Margins Program, the Maryland Geological Survey investigated the sedimentological, paleontological, stratigraphical and geophysical character of Maryland's inner continental shelf. Based on seismic records and sedimentological analysis completed during the first four years, a late quaternary stratigraphic model was developed. Five distinct stratigraphic units were identified and described on the Maryland inner shelf. These units represent late Pleistocene interglacial deposits, the oldest of which corresponds to pre-Illinoian transgressive shelf sands.

Marine Georesources and Geotechnology

Wave-current interaction in Willapa Bay

This paper describes the importance of wave-current interaction in an inlet-estuary system. The three-dimensional, fully coupled, Coupled Ocean-Atmosphere-Wave-Sediment Transport (COAWST) modeling system was applied in Willapa Bay (Washington State) from 22 to 29 October 1998 that included a large storm event. To represent the interaction between waves and currents, the vortex-force method was used. Model results were compared with water elevations, currents, and wave measurements obtained by the U.S. Army Corp of Engineers. In general, a good agreement between field data and computed results was achieved, although some discrepancies were also observed in regard to wave peak directions in the most upstream station. Several numerical experiments that considered different forcing terms were run in order to identify the effects of each wind, tide, and wave-current interaction process. Comparison of the horizontal momentum balances results identified that wave-breaking-induced acceleration is one of the leading terms in the inlet area. The enhancement of the apparent bed roughness caused by waves also affected the values and distribution of the bottom shear stress. The pressure gradient showed significant changes with respect to the pure tidal case. During storm conditions the momentum balance in the inlet shares the characteristics of tidal-dominated and wave-dominated surf zone environments. The changes in the momentum balance caused by waves were manifested both in water level and current variations. The most relevant effect on hydrodynamics was a wave-induced setup in the inner part of the estuary.

Washington

The NASA hydrological forecast system for food and water security applications

Many regions in Africa and the Middle East are vulnerable to drought and to water and food insecurity, motivating agency efforts such as the U.S. Agency for International Development’s (USAID) Famine Early Warning System Network (FEWS NET) to provide early warning of drought events in the region. Each year these warnings guide life-saving assistance that reaches millions of people. A new NASA multi-model, remote sensing-based hydrological forecasting and analysis system, NHyFAS, has been developed to support such efforts by improving the FEWS NET’s current early warning capabilities. NHyFAS derives its skill from two sources: (i) accurate initial conditions, as produced by an offline land modeling system through the application and/or assimilation of various satellite data (precipitation, soil moisture, and terrestrial water storage); and (ii) meteorological forcing data during the forecast period as produced by a state-of-the-art ocean-land-atmosphere forecast system. The land modeling framework used is the Land Information System (LIS), which employs a suite of land surface models, allowing multi-model ensembles and multiple data assimilation strategies to better estimate land surface conditions. An evaluation of NHyFAS shows that its one-to-five month forecasts successfully capture known historic drought events. The system also benefits from strong collaboration with end-user partners in Africa and the Middle East, who provide insights on strategies to formulate and communicate early warning indicators to water and food security communities. The additional lead time provided by this system will increase the speed, accuracy and efficacy of humanitarian disaster relief, helping to save lives and livelihoods.

Bulletin of the American Meteorological Society

Stratiform barite deposits in the Roberts Mountains allochthon, Nevada: A review of potential analogs in modern sea-floor environments

The United States is a net importer of barite, a critical mineral for the oil and gas industry; more than 80 percent of current domestic consumption of barite is imported from China. Nearly all of the domestic production of barite comes from stratiform deposits in Nevada. The 'modern analogs' approach adopted in this review can contribute to improving deposit models and the long-term resource picture in the United States. Massive barite deposits in Nevada are interlayered with deep-water siliceous sedimentary strata of Paleozoic age within the Roberts Mountains Allochthon. Although the barite deposits formed along the long-lived, tectonically active margin of western North America before the Antler orogeny, uncertainty still exists regarding key aspects of their genesis, especially with respect to tectonic setting and depositional processes. Proposed tectonic settings include a continental slope adjacent to an ocean basin, and a rifted basin formed on continental crust. A margin dominated by other stress configurations and strike-slip faulting may also have been present during the Paleozoic. Hypotheses for the genesis of sediment-hosted, stratiform barite deposits in Nevada can be grouped into two categories: (1) a synsedimentary hydrothermal model and (2) an ocean-circulation/productivity-zone model. Both models include a stage involving bacterial reduction of seawater sulfate. Disseminated, diagenetic, and hydrothermal barite deposits are widespread features of the modern ocean. In the Pacific Ocean, disseminated barite deposits (max 9 weight percent BaSO4) are forming below high-productivity zones (for example, in equatorial belts) and on the flanks of ocean ridges (as fallout of 'black smoker' particles). Massive barite deposits of diagenetic (in the subsea floor and at cold seeps) and hydrothermal origins are present in several sediment-covered tectonic settings, including ocean ridges (Escanaba Trough and Guaymas Basin), oceanic transform faults (Blanco Fracture Zone), marginal basins (Sea of Okhotsk), convergent margins (Peru, Oregon, Alaska), and transform margins (California Continental Borderland). Both hydrothermal and diagenetic barite deposits may be present in some environments (for example, the California California Borderland). On the basis of a consideration of tectonic settings and a comparison of deposit attributes (associated rock types, size, structure), mineralogy (BaSO4 content, SiO2 content), and geochemistry (S- and Sr-isotopic ratios) for modern and ancient massive barite deposits, cold seeps along transform margins (or, possibly, marginal basins) represent the most promising present-day metallogenetic analogs for stratiform barite. Hydrothermal systems can also produce high-grade barite, but the ubiquity of associated sulfide mineralization on the modern sea floor (and the paucity of sulfides in Nevada deposits) is problematic. On the basis of data from modern barite deposits, the presence of vent-specific faunas (tubeworms) and the variation in d34S values for barite (related to bacterial reduction of seawater sulfate) may not permit discrimination between a diagenetic or hydrothermal origin for ancient barite deposits.

Nevada

Calibration of carbonate-water triple oxygen isotope fractionation: Seeing through diagenesis in ancient carbonates

High precision triple oxygen isotope measurements of carbonates can better constrain temperatures and oxygen isotope compositions of seawater through geologic time than 18 O/ 16 O measurements alone, but lack of a definitive calibration has hindered progress. In this study, we fluorinated both carbonate and water samples to measure quantitatively the triple oxygen isotope composition of each phase. We compared the oxygen isotope fractionation between carbonate and water for different carbonate materials: calcite synthesized with and without carbonic anhydrase, abiogenic calcite from Devils Hole, and extant biogenic calcite and aragonite of marine origin. We found similar 1000 lnα 18 O cc-wt values for all materials and combined the results with the high temperature experimental data of O'Neil et al. (1969) , resulting in the following fractionation equation ( T in Kelvins) 1000 ln &#x3B1; 18 O cc - w t = 2.84 ( &#xB1; 0.02 ) &#xD7; 10 6 T 2 - 2.96 ( &#xB1; 0.19 ) "> 1000lnα18Occ-wt=2.84(±0.02)×106T2-2.96(±0.19) . The calcite triple oxygen isotope values yielded a θ - T relationship of θ cc-wt = –1.39(±0.01)/ T + 0.5305 whereas the aragonite triple oxygen isotope values yielded a θ - T relationship of θ ara-wt = –1.53(±0.02)/ T + 0.5305. The calcite-water triple oxygen isotope equilibrium fractionation equation for natural samples is &#x394; 17 &#x2032; O cc - &#x394; 17 &#x2032; O wt = 2.84 ( &#xB1; 0.02 ) &#xD7; 10 6 T 2 - 2.96 ( &#xB1; 0.19 ) - 1.39 ( &#xB1; 0.01 ) T + 0.5305 - &#x3BB; "> Δ17′Occ-Δ17′Owt=2.84(±0.02)×106T2-2.96(±0.19)-1.39(±0.01)T+0.5305-λ . The combined 1000 lnα 18 O and 1000 lnα 17 O relationships can be used to assess equilibrium in ancient samples and to evaluate potential secular changes in the δ 18 O value of seawater. Most of the Phanerozoic samples analyzed in this study, which were determined to be pristine in previous studies, have undergone some level of diagenesis. Two samples appear to preserve their original oxygen isotope compositions and suggest a cool ocean with a δ 18 O value similar to the modern ocean. Using a fluid-rock interaction model, we can “see through” the diagenetic process and estimate the triple oxygen isotope composition of the carbonate prior to alteration. In doing so, we show that for the time intervals and sample locations measured in this study, Phanerozoic oceans had a comparable range of oxygen isotope compositions and temperatures as modern seawater.

Geochemica et Cosmochimica Acta