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.

SourceResearch connectionFrom yearThrough year

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

At least 1,207 records · Page 67Linked to original sources

USGS2010

Impacts and predictions of coastal change during hurricanes

Beaches serve as a natural barrier between the ocean and inland communities, ecosystems, and resources. These dynamic environments move and change in response to winds, waves, and currents. During a powerful hurricane, changes to beaches can be large, and the results are sometimes catastrophic. Lives are lost, communities are destroyed, and millions of dollars are spent on rebuilding. There is a clear need to identify areas of our coastline that are likely to experience extreme and devastating erosion during a hurricane. It is also important to determine risk levels associated with development in areas where the land shifts and moves with each landfalling storm. The U.S. Geological Survey (USGS) provides scientific support for hurricane planning and response. Using observations of beach changes and models of waves and storm surge, we are predicting how the coast will respond to hurricanes and identifying areas vulnerable to extreme coastal changes.

Fact Sheet↗

USGS2023

Bayesian hierarchical modeling for probabilistic estimation of tsunami amplitude from far-field earthquake sources

Evaluation of tsunami disaster risk for a coastal region requires reliable estimation of tsunami hazard, for example, wave amplitude close to the shore. Observed tsunami data are scarce and have poor spatial coverage, and for this reason probabilistic tsunami hazard analysis (PTHA) traditionally relies on numerical simulation of “synthetic” tsunami generation and propagation toward the coast. Such an approach has been extensively studied in the past and it is widely recognized as an important disaster-risk mitigation tool. PTHA can not only provide less uncertain and spatially coherent hazard estimates in comparison with classical empirical data analysis which is restricted at the tide gauge stations, but also local inundation information. In this paper, we explore a purely statistical alternative to traditional PTHA for evaluation of tsunami amplitude hazard. Here, we use tide gauge measurements of tsunami amplitude along the western United States, specifically California and Oregon, and develop a spatial Bayesian hierarchical model (BHM) to assess tsunami hazard from far-field earthquake sources at various recurrence intervals. The configuration of our model incorporates latent Gaussian fields that utilize information on the distance between tide gauges as well as on the continental shelf width, that is, a covariate linked to potential dissipative effects on wave energy as the tsunami travels over shallow water. Through our BHM, we produce spatially continuous probabilistic maps of far-field tsunami hazard which can aid comprehensive tsunami disaster risk reduction and management.

California, Oregon, Washington↗

USGS2024

Isotopic evidence against North Pacific Deep Water formation during late Pliocene warmth

Several modelling and observational studies suggest deep water formation in the subpolar North Pacific as a possible alternative mode of thermohaline circulation that occurred in the warm Pliocene, a time when global atmospheric partial pressure of carbon dioxide was like the modern atmosphere (~400 ppm). We test this hypothesis by measuring the δ 13 C of the benthic foraminifer Cibicidoides wuellerstorfi collected from northernmost Pacific mid-Piacenzian Warm Period (3.264–3.025 Myr ago) sediments. The data reveal progressively more isotopically negative dissolved inorganic carbon along a northward Equator-to-pole transect, the opposite of the expected Pliocene Pacific meridional overturning circulation signal. C. wuellerstorfi δ 13 C is also often more positive at the deeper Ocean Drilling Program (ODP) site 887 compared with the shallower ODP site 883, suggesting ‘bottom-up’ ventilation of the deep Pacific Ocean. We then present alkenone sea surface temperature and export-productivity data from ODP site 883, which suggest that late Pliocene subarctic North Pacific carbonate sedimentation was, at least in part, probably due to higher coccolithophore export production, rather than North Pacific Deep Water formation as previously argued. Therefore, we suggest it is unlikely that North Pacific Deep Water formation occurred in the mid-Piacenzian Warm Period, although a shallower overturning cell cannot be ruled out.

Nature Geoscience↗

USGS2006

Proximate composition, energetic value, and relative abundance of prey fish from the inshore eastern Bering Sea: Implications for piscivorous predators

Changing ocean conditions and subsequent shifts in forage fish communities have been linked to numerical declines of some piscivorous marine birds and mammals in the North Pacific. However, limited information about fish communities is available for some regions, including nearshore waters of the eastern Bering Sea, where many piscivores reside. We determined proximate composition and energetic value of a suite of potential forage fish collected from an estuary on the Yukon-Kuskokwim Delta, Alaska, during 2002 and 2003. Across species, energy density ranged from 14.5 to 20.7 kJ g −1 dry mass and varied primarily as a function of lipid content. Total energy content was strongly influenced by body length and we provide species-specific predictive models of total energy based on this relationship; some models may be improved further by incorporating year and date effects. Based on observed energetic differences, we conclude that variation in fish size, quantity, and species composition of the prey community could have important consequences for piscivorous predators.

Polar Biology↗

USGS2005

Seismic precursory patterns before a cliff collapse and critical point phenomena

We analyse the statistical pattern of seismicity before a 1-2 103 m3 chalk cliff collapse on the Normandie ocean shore, Western France. We show that a power law acceleration of seismicity rate and energy in both 40 Hz-1.5 kHz and 2 Hz-10kHz frequency range, is defined on 3 orders of magnitude, within 2 hours from the collapse time. Simultaneously, the average size of the seismic events increases toward the time to failure. These in situ results are derived from the only station located within one rupture length distance from the rock fall rupture plane. They mimic the "critical point" like behavior recovered from physical and numerical experiments before brittle failures and tertiary creep failures. Our analysis of this first seismic monitoring data of a cliff collapse suggests that the thermodynamic phase transition models for failure may apply for cliff collapse. Copyright 2005 by the American Geophysical Union.

Geophysical Research Letters↗

USGS2020

Geoacoustic inversion for a New England mud patch sediment using the silt-suspension theory of marine mud

This article provides an application of the silt-suspension theory to a Bayesian-inference inversion for the geo-acoustic parameters in marine mud. The theory, with consequences that have been developed recently, postulates a suspension of water and clay mineral card-houses that supports moderately dilute concentrations of silt particles. The approach is an example of a physically based model inversion, in which parameters representing physical mud-layer properties are obtained by inversion and used to produce estimates of geoacoustic properties, including their frequency dependence. The acoustic data are from a combustive source signal propagated along a track, located over several meters of fine-grained mud in the New England Mud Patch, to a single hydrophone on a receiver array during the 2017 Seabed Characterization Experiment. Data extracted from a nearby piston core inform the physical modeling, with selections of inversion parameters guided by both sensitivity analyses and bounds from archival and core measurements. Results show the feasibility of this inversion approach. The estimates of mud density and sound speed are close to values obtained independently. The frequency dependence of attenuation is estimated over the full low-frequency source band and has an approximate power exponent of 1.72.

IEEE Journal of Oceanic Engineering↗

USGS1990

Integration of COCORP deep reflection and magnetic anomaly analysis in the southeastern United States: Implications for origin of the Brunswick and East Coast magnetic anomalies: Alternative interpretation and reply

Integration of magnetic anomaly analysis with COCORP deep reflection data from the southeastern United States provides three new constraints on the interpretation of the Brunswick and East Coast magnetic anomalies, as well as on the reflection data. These are as follows. (1) The source of the Brunswick anomaly lies within the deep crust. This anomaly is not caused by a Mesozoic rift basin, as proposed by some workers. (2) A simple, seaward-dipping, high- susceptibility slab model can explain both the Brunswick and East Coast magnetic anomalies. The along-strike change in character of the two anomalies results largely from a change in azimuth of the source body. (3) Beneath the southeastern United States, this source body dips south, lies immediately on the south flank of the prominent southward-dipping reflective zone revealed on COCORP surveys, and was previously associated with the Alleghanian suture between North America and Africa. These results imply that a dipping, highly magnetized zone in the upper plate of the Alleghanian suture is responsible for both the Brunswick and East Coast magnetic anomalies. The high- susceptibility material responsible for these anomalies might be mafic lower continental or oceanic crust thrust upward during Alleghanian continental collision, or mafic igneous material intruded into the upper plate of the suture zone during subsequent Mesozoic rifting, or both. The latter hypothesis implies that the Alleghanian suture acted, as a zone of weakness (a repository ?) which was reactivated to control the site of ultimate Atlantic rifting and possibly initial sea-floor spreading.

Geological Society of America Bulletin↗

USGS2025

Survival, travel time, and use of migration routes by juvenile steelhead in a modified river estuary

Greater understanding of the survival, travel time, and spatial distribution of juvenile salmonids among migration routes between their natal streams and the ocean is critical to the recovery of these threatened species. In the Sacramento–San Joaquin River Delta (Delta), a highly modified estuary in central California, USA, there is a critical need to evaluate how water management (e.g., water pumping) and environmental factors (e.g., water flow) impact these populations. While management actions can affect some environmental variables in the Delta, only recently have studies begun to uncover associations between these variables and key demographic parameters. In this study, we examine the effects of freshwater flows, water exports, tidal environment, and a temporary barrier on juvenile steelhead ( Oncorhynchus mykiss ) survival, travel times, and migration routing using a multiyear acoustic telemetry dataset and recent advancements in Bayesian multistate mark-recapture modeling. We found that no single covariate explained variation in juvenile steelhead population dynamics across the entire Delta, but that separate regions within the Delta showed association with specific environmental factors.

California↗

USGS2005

Authigenesis of trace metals in energetic tropical shelf environments

We evaluated authigenic changes of Fe, Mn, V, U, Mo, Cd and Re in suboxic, periodically remobilized, tropical shelf sediments from the Amazon continental shelf and the Gulf of Papua. The Cd/Al, Mo/Al, and U/Al ratios in Amazon shelf sediments were 82%, 37%, and 16% less than those in Amazon River suspended sediments, respectively. Very large depletions of U previously reported in this environment were not observed. The Cd/Al ratios in Gulf of Papua sediments were 76% lower than measurements made on several Papua New Guinea rivers, whereas U/Al ratios in the shelf sediments were enriched by approximately 20%. Other metal/Al ratios in the Papua New Guinea river suspended sediments and continental shelf sediments were not distinguishably different. Comparison of metal/Al ratios to grain size distributions in Gulf of Papua samples indicates that our observations cannot be attributed to differences in grain size between the river suspended sediments and continental shelf sediments. These two shelves constitute a source of dissolved Cd to the world ocean equal to 29–100% of the dissolved Cd input from rivers, but only 3% of the dissolved Mo input and 4% of the dissolved U input. Release of Cd, Mo, and U in tropical shelf sediments is likely a result of intense Fe and Mn oxide reduction in pore waters and resuspension of the sediments. Since we do not observe depletions of particulate Fe and Mn in the shelf sediments most of these dissolved metals must reoxidize in the overlying waters and reprecipitate. As Cd exhibits the largest losses on these tropical shelves, we examined the ability of newly formed Fe and Mn oxides to adsorb dissolved Cd using a geochemical diffuse double-layer surface complexation model and found the oxide surfaces are relatively ineffective at readsorbing Cd in seawater due to surface-site competition by Mg and Ca. If the remobilization and reoxidation of Fe and Mn occurs frequently enough before sediment is buried significant amounts of Cd may be removed from the oxide surfaces. Because a much greater percentage of Mn than Fe becomes remobilized in these shelf sediments, metals closely associated with Mn oxides (like Cd) are more likely to show losses during deposition.

Continental Shelf Research↗

USGS2005

Thermal state and complex geology of a heterogeneous salty crust of Jupiter's satellite, Europa

The complex geology of Europa is evidenced by many tectonic and cryomagmatic resurfacing structures, some of which are "painted" into a more visible expression by exogenic alteration processes acting on the principal endogenic cryopetrology. The surface materials emplaced and affected by this activity are mainly composed of water ice in some areas, but in other places there are other minerals involved. Non-ice minerals are visually recognized by their low albedo and reddish color either when first emplaced or, more likely, after alteration by Europan weathering processes, especially sublimation and alteration by ionizing radiation. While red chromophoric material could be due to endogenic production of solid sulfur allotropes or other compounds, most likely the red substance is an impurity produced by radiation alteration of hydrated sulfate salts or sulphuric acid of mainly internal origin. If the non-ice red materials or their precursors have a source in the satellite interior, and if they are not merely trace contaminants, then they can play an important role in the evolution of the icy crust, including structural differentiation and the internal dynamics. Here we assume that these substances are major components of Europa's cryo/hydrosphere, as some models have predicted they should be. If this is an accurate assumption, then these substances should not be neglected in physical, chemical, and biological models of Europa, even if major uncertainties remain as to the exact identity, abundance, and distribution of the non-ice materials. The physical chemical properties of the ice-associated materials will contribute to the physical state of the crust today and in the geological past. In order to model the influence of them on the thermal state and the geology, we have determined the thermal properties of the hydrated salts. Our new lab data reveal very low thermal conductivities for hydrated salts compared to water ice. Lower conductivities of salty ice would produce steeper thermal gradients than in pure ice. If there are salt-rich layers inside the crust, forming salt beds over the seafloor or a briny eutectic crust, for instance, the high thermal gradients may promote endogenic geological activity. On the seafloor, bedded salt accumulations may exhibit high thermochemical gradients. Metamorphic and magmatic processes and possible niches for thermophilic life at shallow suboceanic depths result from the calculated thermal profiles, even if the ocean is very cold. ?? 2004 Elsevier Inc. All rights reserved.

Icarus↗

USGS2021

Late Quaternary deglaciation of Prince William Sound, Alaska

To understand the timing of deglaciation of the northernmost marine-terminating glaciers of the Cordilleran Ice Sheet (CIS), we obtained 26 10 Be surface-exposure ages from glacially scoured bedrock surfaces in Prince William Sound (PWS), Alaska. We sampled six elevation transects between sea level and 620 m and spanning a distance of 14 to 70 km along ice flow paths. Most transect age–elevation patterns could not be explained by a simple model of thinning ice; the patterns provide evidence for lingering ice cover and possible inheritance. A reliable set of 20 ages ranges between 17.4 ± 2.0 and 11.6 ± 2.8 ka and indicates ice receded from northwestern PWS around 14.3 ± 1.6 ka, thinned at a rate of ~120–160 m/ka, and retreated from sea-level sites at 12.9 ± 1.1 ka at a rate of 20 m/yr. The retreat rate likely slowed as glaciers retreated into northern PWS. These results are consistent with the growing body of reported deglacial constraints on collapse of ice sheets along the Alaska margin indicating collapse of the CIS soon after 17 ka. These data are consistent with paleotemperature data indicating that a warming North Pacific Ocean caused catastrophic collapse of this part of the CIS.

Alaska↗

USGS1975

Structure and Paleozoic stratigraphy of a complex of thrust plates in the Fish Creek Reservoir area, south-central Idaho

Permian, Pennsylvanian, Mississippian, Devonian, and Silurian marine rocks of diverse facies are brought together in a complex of six thrust sheets in the Fish Creek Reservoir area on the north edge of the Snake River Plain, Idaho. The lowest structural element, the parautochthon, is made of more than 450 m (1,500 ft) of folded and faulted Devonian miogeosynclinal carbonate rocks present in a 6.5-km 2 (2.5-mi 2 ) window. Along the east margin of the window, a sliver of continental margin transitional carbonate rocks of Early Devonian and Late Silurian age assigned to the Roberts Mountains Formation is thrust over the miogeosynclinal rocks. The window of middle Paleozoic rocks is overridden along the Fish Creek thrust fault by the flysch facies of the Copper Basin Formation, a turbidite-submarine-fan sequence more than 1,000 m (3,300 ft) thick, of Mississippian age. About 4.8 km (3 mi) southwest of the window, about 100 m (300 ft) of deepwater siliceous oceanic facies clastic rocks are exposed, which are assigned with question to the Milligen(?) Formation of Devonian age. These clastic rocks are interpreted to be thrust over the Copper Basin Formation. The highest structural elements are sequences more than 610 m (2,000 ft) thick of interbedded sandy and conglomeratic limestones, quartzites, and conglomerates and interbedded siltstones and argillites of the Wood River Formation of Middle Pennsylvanian to Early Permian age. The Wood River Formation is in thrust contact with the Milligen(?) Formation in the southwest part of the mapped area and with Copper Basin Formation along the west side of Fish Creek Reservoir. All the thrust sheets have moved eastward. The minimum distance moved is estimated from sedimentation models and facies reconstructions to range from perhaps several kilometres for the allochthon of the Roberts Mountains Formation to 48 km (30 mi) for the Milligen(?) Formation allochthon. The principal period of thrusting was post-Early Permian (post-Wood River Formation) and preEocene (pre-Challis Volcanics) and is of probable Sevier age. Middle Paleozoic rocks of the Milligen and Roberts Mountains Formations, however, also may have been involved in an earlier period of thrusting of latest Devonian to earliest Mississippian age related to the Antler orogeny. The thrust sheets were deformed into a northwest-trending dome in late Mesozoic time and were broken by basin-range faults during the Tertiary.

Idaho↗

USGS1998

Ground Water in the Southern Lihue Basin, Kauai, Hawaii

A multi-phased study of ground-water resources, including well drilling, aquifer tests, analysis of ground-water discharge, and numerical ground-water modeling, indicates that the rocks of the southern Lihue Basin, Kauai, have permeabilities that are much lower than in most other areas of ground-water development in the Hawaiian islands. The regional hydraulic conductivity of the Koloa Volcanics, which dominates fresh ground-water flow in the basin, is about 0.275 foot per day. The Waimea Canyon Basalt which surrounds the basin and underlies the Koloa Volcanics within the basin is intruded by dikes that reduce the bulk hydraulic conductivity of the rocks to about 1.11 feet per day. The low permeabilities result in steeper head gradients compared with other areas in the Hawaiian islands, and a higher proportion of ground-water discharging to streams than to the ocean. Water levels rise from near sea level at the coast to several hundreds of feet above sea level at the center of the basin a few miles inland. The high inland water levels are part of a completely saturated ground-water system. Because of the low regional hydraulic conductivity and high influx of water from recharge in the southern Lihue Basin, the rocks become saturated nearly to the surface and a variably saturated/unsaturated (perched) condition is not likely to exist. Streams incising the upper part of the aquifer drain ground water and keep the water levels just below the surface in most places. Streams thus play an important role in shaping the water table in the southern Lihue Basin. At least 62 percent of the ground water discharging from the aquifer in the southern Lihue Basin seeps to streams; the remainder seeps directly to the ocean or is withdrawn by wells.

Water-Resources Investigations Report↗

USGS1997

Review of electric and magnetic fields accompanying seismic and volcanic activity

New observations of magnetic, electric and electromagnetic field variations, possibly related to recent volcanic and seismic events, have been obtained on Mt. Unzen in Japan, Reunion Island in Indian Ocean, the Long Valley volcanic caldera in California, and for faults in China and Russia, California and several other locations. For volcanic events, contributions from different physical processes can be identified during the various eruption stages. Slow processes (weeks to months) include near-surface thermal demagnetization effects, piezomagnetic effects, and effects from rotation/displacement of magnetized material. Rapid processes (seconds to days) include piezomagnetic effects from instantaneous stress redistribution with explosive eruptions and electrokinetic effects from rupture of high pressure fluid compartments commonly encountered in volcanic regions. For seismic events, the observed coseismic offsets are instantaneous, provided care has been taken to ensure sensors are insensitive to seismic shaking and are in regions of low magnetic field gradient. Simple piezomagnetic dislocation models based on geodetically and seismically determined fault parameters generally match the observed signals in size and sign. Electrokinetic effects resulting from rupture of fluid filled compartments at hydrostatic to lithostatic pore pressures can generate transient signals in the frequency band 100 Hz to 0.01 Hz. However, large-scale fluid driven processes are not evident in near-field measurements in the epicentral region minutes to weeks before large earthquakes. The subset of ionospheric disturbances generated by trapped atmospheric pressure waves (also termed gravity waves and/or acoustic waves, traveling ionospheric disturbances or TID's) that are excited by earthquakes and volcanic eruptions are common and propagate to great distances. These are known and expected consequences of earthquakes, volcanic explosions (and other atmospheric disturbances), that must be identified and their effects removed from VLF/ULF electromagnetic field records before associating new observations of ionospheric disturbances with earthquake activity.

Surveys in Geophysics↗

USGS2020

Marine latitudinal diversity gradients, niche conservatism and out of the tropics and Arctic: Climatic sensitivity of small organisms

Aim The latitudinal diversity gradient (LDG) is a consequence of evolutionary and ecological mechanisms acting over long history, and thus is best investigated with organisms that have rich fossil records. However, combined neontological‐palaeontological investigations are mostly limited to large, shelled invertebrates, which keeps our mechanistic understanding of LDGs in its infancy. This paper aims to describe the modern meiobenthic ostracod LDG and to explore the possible controlling factors and the evolutionary mechanisms of this large‐scale biodiversity pattern. Location Present‐day Western North Atlantic. Taxon Ostracoda. Methods We compiled ostracod census data from shallow‐marine environments of the western North Atlantic Ocean. Using these data, we documented the marine LDG with multiple metrics of alpha, beta (nestedness and turnover) and gamma diversity, and we tested whether macroecological patterns could be governed by different environmental factors, including temperature, salinity, dissolved oxygen, pH and primary productivity. We also explored the geologic age distribution of ostracod genera to investigate the evolutionary mechanisms underpinning the LDG. Results Our results show that temperature and climatic niche conservatism are important in setting LDGs of these small, poorly dispersing organisms. We also found evidence for some dispersal‐driven spatial dynamics in the ostracod LDG. Compared to patterns observed in marine bivalves, however, dispersal dynamics were weaker and they were bi‐directional, rather than following the ‘out‐of‐the‐tropics’ model. Main conclusions Our detailed analyses revealed that meiobenthic organisms, which comprise two‐thirds of marine diversity, do not always follow the same rules as larger, better‐studied organisms. Our findings suggest that the understudied majority of biodiversity may be more sensitive to climate than well‐studied, large organisms. This implies that the impacts of ongoing Anthropocene climatic change on marine ecosystems may be much more serious than presently thought.

Journal of Biogeography↗

USGS2020

Sediment budget estimates for a highly impacted embayment with extensive wetland loss

External sediment supply is an important control on wetland morphology and vulnerability to storms, sea-level rise, and land use change. Constraining sediment supply and net budgets is difficult due to multiple timescales of variability in hydrodynamic forcing and suspended-sediment concentrations, as well as the fundamental limitations of measurement and modeling technologies. We used two independent observational campaigns and one hydrodynamic modeling effort to estimate the sediment supply to Jamaica Bay, New York, USA, an urbanized embayment with a history of extensive wetland loss. We found that all three estimates indicate a net import to the system, ranging from 36 x 106 – 74 x 106 kg/y, with a mean estimate of 55,000 t/y +/- 31,000 t/y, which compares well with a prior estimate derived from radionuclide tracers. Net sediment import is controlled by flood-ebb asymmetry in bed shear stress, which results in higher suspended sediment concentrations on flood tide relative to ebb. This indicates a seaward source of sediment, likely offshore marine deposits and potentially sediment from the adjacent Hudson River-Estuary that is resuspended by waves in the coastal ocean. Despite the net sediment import, a simple sediment budget suggests that the rate of supply is not sufficient to maintain the present geomorphic planform of the system relative to sea-level rise. The convergent estimates from independent methods provide reasonable guidance as context for sediment-based restoration efforts.

New York↗

USGS2018

Flow-mediated effects on travel time, routing, and survival of juvenile Chinook salmon in a spatially complex, tidally forced river delta

We evaluated the interacting influences of river flows and tides on travel time, routing, and survival of juvenile late-fall Chinook salmon ( Oncorhynchus tshawytscha ) migrating through the Sacramento–San Joaquin River Delta. To quantify these effects, we jointly modeled the travel time, survival, and migration routing in relation to individual time-varying covariates of acoustic-tagged salmon within a Bayesian framework. We used observed arrival times for detected individuals and imputed arrival times for undetected individuals to assign covariate values in each reach. We found travel time was inversely related to river inflow in all reaches, yet survival was positively related to inflow only in reaches that transitioned from bidirectional tidal flows to unidirectional flow with increasing inflows. We also found that the probability of fish entering the interior Delta, a low-survival reach, declined as inflow increased. Our study illustrates how river inflows interact with tides to influence fish survival during the critical transition between freshwater and ocean environments. Furthermore, our analytical framework introduces new techniques to integrate formally over missing covariate values to quantify effects of time-varying covariates.

California↗

USGS2000

Parent brine of the castile evaporites (Upper Permian), Texas and New Mexico

The Upper Permian (lower Ochoan) Castile Formation is a major evaporite sequence (∼10,000 km 3 ) of calcite, anhydrite, and halite in west Texas and southeastern New Mexico. Traditionally the Castile brine has been considered to have been derived from seawater. This tradition has recently been challenged by two versions of the closed-basin drawdown model. They call for deposition from a mixed brine, in part marine and in large part nonmarine. They propose drawdown of as much as 500 m to form a major sink for ground water issuing from the surrounding Capitan reef complex. A large fraction of the solute in the brine body is inferred to have been recycled from older Permian evaporites on the surrounding shelf. Strontium-isotope analyses show no evidence that meteoric ground water was contributed to the Castile brine. From a stratigraphic, geographic, and lithologic array of 65 samples of anhydrite, gypsum, and calcite, 59 have an 87 Sr/ 86 Sr ratio of 0.706923 (Δsw of -225.0), a ratio that is the same as that of strontium in early Ochoan ocean water. If considerable (>15%) influx of meteoric water had occurred, enough continental strontium would have been introduced to have resulted in higher ratios. Low bromide values (20-40 ppm) in Castile halite, which have been used to argue for meteoric influx and for recycled salt, probably resulted from diagenesis. During shallow burial by halite, centimeter-size, bottom-grown crystals of gypsum were altered to nodular anhydrite. The rising water of dehydration caused the halite to recrystallize. During the recrystallization, some bromide was expelled. Despite the large volume of water that evaporated annually from its surface (∼52 km 3 /yr, assuming an evaporation rate of 2 m/yr), the Castile brine body never completely desiccated. The surrounding shelf was flat, hot, and generally dry. It probably could not have supplied a significant volume of meteoric spring water to the basin over tens of thousands of years. More likely, during the entire history of the evaporite sequence, influx was dominantly marine. Marine ground water flowed through the Capitan Formation into the evaporite basin along its southern and possibly western margin probably with a rate of flow that was usually fast enough to prevent major drawdown of the brine surface.

Journal of Sedimentary Research↗