USGS Science⌕ Search

SEARCH · USGS Science

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 847 records · Page 47Linked to original sources
USGS2011

Hydrogeology and simulation of groundwater flow in fractured rock in the Newark basin, Rockland County, New York

Groundwater in the Newark basin aquifer flows primarily through discrete water-bearing zones parallel to the strike and dip of bedding, whereas flow perpendicular to the strike is restricted, thereby imparting anisotropy to the groundwater flow field. The finite-element model SUTRA was used to represent bedrock structure in the aquifer by spatially varying the orientation of the hydraulic conductivity tensor to reflect variations in the strike and dip of the bedding. Directions of maximum and medium hydraulic conductivity were oriented parallel to the bedding, and the direction of minimum hydraulic conductivity was oriented perpendicular to the bedding. Groundwater flow models were prepared to simulate local flow in the vicinity of the Spring Valley well field and regional flow through the Newark basin aquifer. The Newark basin contains sedimentary rocks deposited as alluvium during the Late Triassic and is one of a series of basins that developed when Mesozoic rifting of the super continent Pangea created the Atlantic Ocean. The westward-dipping basin is filled with interbedded facies of coarse-grained to fine-grained rocks that were intruded by diabase associated with Jurassic volcanism. The Newark basin aquifer is bounded to the north and east by the Palisades sill and to the west by the Ramapo Fault. Although the general dip of bedding is toward the fault, mapping of conglomerate beds indicates the rocks are folded into broad anticlines and synclines. An alternative, more uniform pattern of regional structure, based on interpolated strike and dip measurements from a number of sources, has also been proposed. Two groundwater flow models (A for the former type of bedrock structure and B for the latter type) were developed to represent these alternative depictions of bedrock structure. Transient simulations were calibrated to reproduce measured water-level recoveries in a 9.3 mi&sup2 area surrounding the Spring Valley well field during a 5-day aquifer test in 1992. The models represented a 330-ft thick rock mass divided vertically into 10 equally spaced layers and were calibrated through nonlinear regression. Results of model B best matched the observed water-level recoveries with an estimated hydraulic conductivity of 9.5 ft/day, specific storage of 7.6 x 10 -6 ft -1, and K max : K min anisotropy ratio (hydraulic conductivity parallel to bedding: perpendicular to bedding) of 72:1. Model error was 50 percent greater in model A because the assumed structure did not match the actual strike of bedding in this area. Steady-state simulations of regional flow through the 85.4-mi2 modeled extent of the Newark basin aquifer represented both the alluvial aquifer beneath the Mawah River and the fractured bedrock. The rock mass was divided into two aquifer units: an upper 500-ft thick unit divided into 10 equally spaced layers through which most groundwater is assumed to flow and a lower unit divided into 7 layers with increasing thickness. Models were calibrated through nonlinear regression to average water levels measured in 140 wells from August 2005 through April 2007. Water levels simulated using the two models were similar and generally matched those observed, and the average recharge rate estimated using both models was 19 inches/year for the simulated period. Estimated transmissivity parallel to the strike of bedding (1,100 ft&sup2/d) was uniform in two transmissivity (T) zones in model A, but in model B the transmissivity of a high T zone (1,600 ft&sup2/d), delineated on the basis of aquifer test data, was slightly greater than in a low T zone (1,300 ft&sup2/d). The K max : K min anisotropy was estimated to be 58:1 in model A and 410:1 in model B, so the proportion of flow perpendicular to bedding is less in model B than in model A. Distributions of groundwater age simulated with models A and B are similar and indicate that most shallow ground-water (225 ft below the bedrock surface) is 5 t

Scientific Investigations Report↗
USGS2003

Evolution of CAM and C 4 carbon-concentrating mechanisms

Mechanisms for concentrating carbon around the Rubisco enzyme, which drives the carbon-reducing steps in photosynthesis, are widespread in plants; in vascular plants they are known as crassulacean acid metabolism (CAM) and C 4 photosynthesis. CAM is common in desert succulents, tropical epiphytes, and aquatic plants and is characterized by nighttime fixation of CO2. The proximal selective factor driving the evolution of this CO2-concentrating pathway is low daytime CO2, which results from the unusual reverse stomatal behavior of terrestrial CAM species or from patterns of ambient CO2 availability for aquatic CAM species. In terrestrials the ultimate selective factor is water stress that has selected for increased water use efficiency. In aquatics the ultimate selective factor is diel fluctuations in CO2 availability for palustrine species and extreme oligotrophic conditions for lacustrine species. C4 photosynthesis is based on similar biochemistry but carboxylation steps are spatially separated in the leaf rather than temporally as in CAM. This biochemical pathway is most commonly associated with a specialized leaf anatomy known as Kranz anatomy; however, there are exceptions. The ultimate selective factor driving the evolution of this pathway is excessively high photorespiration that inhibits normal C3 photosynthesis under high light and high temperature in both terrestrial and aquatic habitats. CAM is an ancient pathway that likely has been present since the Paleozoic era in aquatic species from shallow-water palustrine habitats. While atmospheric CO2 levels have undoubtedly affected the evolution of terrestrial plant carbon-concentrating mechanisms, there is reason to believe that past atmospheric changes have not played as important a selective role in the aquatic milieu since palustrine habitats today are not generally carbon sinks, and the selective factors driving aquatic CAM are autogenic. Terrestrial CAM, in contrast, is of increasing selective value under extreme water deficits, and undoubtedly, high Mesozoic CO2 levels reduced the amount of landscape perceived by plants as water limited. Late Tertiary and Quaternary reductions in atmospheric CO2, coupled with increasing seasonality, were probably times of substantial species radiation and ecological expansion for CAM plants. C4 photosynthesis occurs in only about half as many families as CAM, and three-fourths of C4 species are either grasses or sedges. Molecular phylogenies indicate C4 is a more recent innovation than CAM and that it originated in the mid-Tertiary, 20–30 Ma, although some data support an earlier origin. While the timing of the origin of C4 remains controversial, the nearly explosive increase in C4 species is clearly documented in the late Miocene, 4–7 Ma. Increasing seasonality has been widely suggested as an important climatic stimulus for this C4 expansion. Alternatively, based on models of photosynthetic quantum yield at different temperatures and CO2 concentration, it has been hypothesized that the late Miocene C4 expansion resulted from declining atmospheric CO2 levels. This model is most appropriate for explaining the transition from C3 grasslands to C4 grasslands but by itself may not be sufficient to explain the more likely scenario of a late Miocene transition from C3 woodland/ savanna to C4 grasslands. A largely unexplored hypothesis is that climatic changes in late Miocene altered disturbance regimes, in particular the incidence of fires, which today are often associated with maintenance of C4 grasslands. Oceanic charcoal sediments that appear to represent Aeolian deposits from continental wildfires follow a strikingly similar pattern of explosive increase in late Miocene. Climate, CO2, and disturbance are not mutually exclusive explanations and probably all acted in concert to promote the expansion of C4 grasslands. More recently, late Quaternary changes in CO2 may have been responsible for driving major changes in the landscape distribution of C4 species. The theory is sound; however, many of the studies cited in support of this model are open to alternative interpretations, and none has eliminated climatic factors as important selective agents. CAM and C4 evolution required coupling of biochemical pathways with structural changes in photosynthetic tissues, succulence in CAM and Kranz in C4. This was apparently accomplished by piecemeal evolution beginning with mechanisms for recapturing respiratory CO2, although this need not have been so in aquatic CAM species. It has been proposed that the extreme rarity of both pathways in the same plant results from biochemical and structural incompatibilities (Sage 2002). Equally important is the fact that the selective environments are quite different, with CAM evolution thriving on stressful sites inhospitable to C3 species whereas C4 evolution has selected for rapid growth capable of outcompeting associated C3 plants.

International Journal of Plant Sciences↗
USGS2000

Geophysical evidence for the evolution of the California Inner Continental Borderland as a metamorphic core complex

We use new seismic and gravity data collected during the 1994 Los Angeles Region Seismic Experiment (LARSE) to discuss the origin of the California Inner Continental Borderland (ICB) as an extended terrain possibly in a metamorphic core complex mode. The data provide detailed crustal structure of the Borderland and its transition to mainland southern California. Using tomographic inversion as well as traditional forward ray tracing to model the wide-angle seismic data, we find little or no sediments, low (≤6.6 km/s) P wave velocity extending down to the crust-mantle boundary, and a thin crust (19 to 23 km thick). Coincident multichannel seismic reflection data show a reflective lower crust under Catalina Ridge. Contrary to other parts of coastal California, we do not find evidence for an underplated fossil oceanic layer at the base of the crust. Coincident gravity data suggest an abrupt increase in crustal thickness under the shelf edge, which represents the transition to the western Transverse Ranges. On the shelf the Palos Verdes Fault merges downward into a landward dipping surface which separates "basement" from low-velocity sediments, but interpretation of this surface as a detachment fault is inconclusive. The seismic velocity structure is interpreted to represent Catalina Schist rocks extending from top to bottom of the crust. This interpretation is compatible with a model for the origin of the ICB as an autochthonous formerly hot highly extended region that was filled with the exhumed metamorphic rocks. The basin and ridge topography and the protracted volcanism probably represent continued extension as a wide rift until ∼13 m.y. ago. Subduction of the young and hot Monterey and Arguello microplates under the Continental Borderland, followed by rotation and translation of the western Transverse Ranges, may have provided the necessary thermomechanical conditions for this extension and crustal inflow.

California↗
USGS2023

Updates to the regional groundwater-flow model of the New Jersey Coastal Plain, 1980–2013

A 21-layer three-dimensional transient groundwater-flow model of the New Jersey Coastal Plain was developed and calibrated by the U.S. Geological Survey (USGS) in cooperation with the New Jersey Department of Environmental Protection to simulate groundwater-flow conditions during 1980–2013, incorporating average annual groundwater withdrawals and average annual groundwater recharge. This model is the third version of the New Jersey Coastal Plain regional groundwater-flow model that was initially developed as part of the USGS Regional Aquifer System Analysis (RASA) program. The model simulates groundwater flow in 11 aquifers and 10 intervening confining units of the New Jersey Coastal Plain to provide a regional overview of groundwater conditions. Averaged groundwater withdrawal data for 1980 to 2013 were used in the model. The 11 aquifers in New Jersey are, from shallowest to deepest, the Holly Beach water-bearing zone and the confined Cohansey aquifer in Cape May County; the Rio Grande water-bearing zone; the Atlantic City 800-foot sand; the Piney Point, Vincentown, and Wenonah-Mount Laurel aquifers; the Englishtown aquifer system; and the upper, middle, and lower aquifers of the Potomac-Raritan-Magothy (PRM) aquifer system. The model was developed with the MODFLOW–2005 numerical code and the UCODE parameter estimation technique and calibrated using water-level and base-flow observations. A total of 3,453 water-level observations from 392 wells in New Jersey and 48 wells in Delaware from 1983 to 2013 were used in model calibration, which includes historical water-level trends for 29 wells in New Jersey during 1980–2013 presented in time-series hydrographs. In addition, derived observations also were included by calculating the vertical gradient at 33 pairs of nested observation wells in New Jersey, for a total of 210 observations. Changes in water levels over time were calculated for 134 wells in New Jersey and four wells in Delaware where water levels had varied substantially (approximately 10 ft) over the 30-year span of synoptic water-level measurements, for a total of 767 observations. A total of 1,485 base-flow observations in 47 surface-water basins in New Jersey from 1980 to 2013 were used in model calibration. Updates to the groundwater-flow model include the conversion to a fully three-dimensional model from the previous quasi-three-dimensional model. The new model will allow for potential future uses such as particle tracking or simulation of variable-density groundwater flow that could not be accomplished with earlier versions of the model. Spatially and temporally variable recharge estimated by using a soil-water balance model resulted in a spatially and temporally finer discretization. The Rio Grande water-bearing zone was added to the model as an aquifer layer to refine estimates of simulated flow in Atlantic and Cape May Counties, New Jersey. Hydrogeologic parameters were updated to include the confining units in New Jersey and corresponding hydrogeologic units in Delaware and eastern Maryland. The simulated water levels for the New Jersey Coastal Plain aquifers were compared to water-level measurements made during 1980–2013. The average residual for 4,243 water-level observations for New Jersey (simulated water levels minus measured water levels) is 1.5 feet. The simulated water-level contours for the confined aquifers for 2013 were compared to potentiometric surfaces produced from water levels measured during 2013. Simulated water levels generally matched the 2013 potentiometric surfaces of the confined aquifers in the areas of large withdrawals. Hydrographs of wells in the confined Coastal Plain aquifers of New Jersey show that simulated water levels generally match the magnitude and seasonal variation of the observed water levels. Hydrographs of base flow for the 47 streamgaging stations in New Jersey indicate that most of the simulated and estimated data match reasonably well. Groundwater withdrawals are an important resource for water supply, agricultural, industrial, and commercial needs in the New Jersey Coastal Plain. Groundwater withdrawals from the New Jersey Coastal Plain aquifers have resulted in persistent, regionally extensive cones of depression in the Englishtown aquifer system and Wenonah-Mount Laurel aquifer in Ocean and Monmouth Counties; Wenonah-Mount Laurel and upper, middle, and lower PRM aquifers in Camden County; and Atlantic City 800-foot sand in Atlantic County. Because hydrologic stresses and water-management needs change with time, periodic updates to the groundwater-flow model are required to provide current information about hydrologic conditions in the New Jersey Coastal Plain and to maintain its usefulness as a tool to manage water resources and develop water-resource strategies. The current updates will support the continued application of this model as a tool for evaluating the regional effects of changes in groundwater withdrawals and of current and potential future water-management strategies on groundwater levels in the New Jersey Coastal Plain.

New Jersey Coastal Plain↗
USGS2012

Statistical analyses to support guidelines for marine avian sampling. Final report

Interest in development of offshore renewable energy facilities has led to a need for high-quality, statistically robust information on marine wildlife distributions. A practical approach is described to estimate the amount of sampling effort required to have sufficient statistical power to identify species-specific “hotspots” and “coldspots” of marine bird abundance and occurrence in an offshore environment divided into discrete spatial units (e.g., lease blocks), where “hotspots” and “coldspots” are defined relative to a reference (e.g., regional) mean abundance and/or occurrence probability for each species of interest. For example, a location with average abundance or occurrence that is three times larger the mean (3x effect size) could be defined as a “hotspot,” and a location that is three times smaller than the mean (1/3x effect size) as a “coldspot.” The choice of the effect size used to define hot and coldspots will generally depend on a combination of ecological and regulatory considerations. A method is also developed for testing the statistical significance of possible hotspots and coldspots. Both methods are illustrated with historical seabird survey data from the USGS Avian Compendium Database. Our approach consists of five main components: 1. A review of the primary scientific literature on statistical modeling of animal group size and avian count data to develop a candidate set of statistical distributions that have been used or may be useful to model seabird counts. 2. Statistical power curves for one-sample, one-tailed Monte Carlo significance tests of differences of observed small-sample means from a specified reference distribution. These curves show the power to detect "hotspots" or "coldspots" of occurrence and abundance at a range of effect sizes, given assumptions which we discuss. 3. A model selection procedure, based on maximum likelihood fits of models in the candidate set, to determine an appropriate statistical distribution to describe counts of a given species in a particular region and season. 4. Using a large database of historical at-sea seabird survey data, we applied this technique to identify appropriate statistical distributions for modeling a variety of species, allowing the distribution to vary by season. For each species and season, we used the selected distribution to calculate and map retrospective statistical power to detect hotspots and coldspots, and map pvalues from Monte Carlo significance tests of hotspots and coldspots, in discrete lease blocks designated by the U.S. Department of Interior, Bureau of Ocean Energy Management (BOEM). 5. Because our definition of hotspots and coldspots does not explicitly include variability over time, we examine the relationship between the temporal scale of sampling and the proportion of variance captured in time series of key environmental correlates of marine bird abundance, as well as available marine bird abundance time series, and use these analyses to develop recommendations for the temporal distribution of sampling to adequately represent both shortterm and long-term variability. We conclude by presenting a schematic “decision tree” showing how this power analysis approach would fit in a general framework for avian survey design, and discuss implications of model assumptions and results. We discuss avenues for future development of this work, and recommendations for practical implementation in the context of siting and wildlife assessment for offshore renewable energy development projects.

NOAA Technical Memorandum↗
USGS2013

Tidally influenced alongshore circulation at an inlet-adjacent shoreline

The contribution of tidal forcing to alongshore circulation inside the surfzone is investigated at a 7 km long sandy beach adjacent to a large tidal inlet. Ocean Beach in San Francisco, CA (USA) is onshore of a ∼150 km2 ebb-tidal delta and directly south of the Golden Gate, the sole entrance to San Francisco Bay. Using a coupled flow-wave numerical model, we find that the tides modulate, and in some cases can reverse the direction of, surfzone alongshore flows through two separate mechanisms. First, tidal flow through the inlet results in a barotropic tidal pressure gradient that, when integrated across the surfzone, represents an important contribution to the surfzone alongshore force balance. Even during energetic wave conditions, the tidal pressure gradient can account for more than 30% of the total alongshore pressure gradient (wave and tidal components) and up to 55% during small waves. The wave driven component of the alongshore pressure gradient results from alongshore wave height and corresponding setup gradients induced by refraction over the ebb-tidal delta. Second, wave refraction patterns over the inner shelf are tidally modulated as a result of both tidal water depth changes and strong tidal flows (∼1 m/s), with the effect from currents being larger. These tidally induced changes in wave refraction result in corresponding variability of the alongshore radiation stress and pressure gradients within the surfzone. Our results indicate that tidal contributions to the surfzone force balance can be significant and important in determining the direction and magnitude of alongshore flow.

California↗
USGS2025

Rupture process of the Mw7.0 December 5, 2024 Offshore Cape Mendocino earthquake

The Mw7.0 December 5, 2024 Offshore Cape Mendocino earthquake ruptured a km long portion of the east-west trending Mendocino fault zone (MFZ). In order to clarify the rupture process, we assemble three-component seismograms from regional seismic stations, horizontal coseismic displacement vectors derived from Global Navigation Satellite System (GNSS) time series, and a Sentinel-1 ascending interferogram. These data are interpreted with a model of slip distributed on two vertical fault planes representative of the eastern MFZ and spanning the ~70 km length of the aftershock zone. Assuming right-lateral strike slip, we find that the rupture initiates in the oceanic mantle at 20-30 km depth and proceeds unilaterally updip and toward the east. Early aftershocks locate adjacent to the peak slip areas, tracking the coseismic rupture propagation from oceanic mantle to shallower depth and implying a significant role of static stress transfer in driving aftershocks in an ocean plate environment.

California↗
USGS2024

Advancing water security in Africa with new high-resolution discharge data

VegDischarge v1 is a comprehensive river discharge across Africa (2000–2021), produced by coupling the agro-hydrologic VegET model and the mizuRoute routing framework. Using remote sensing data and hydrological modeling, the 1-km runoff field simulated by VegET, and routed with mizuRoute, covers over 64,000 river segments in Africa. The VegET model simulates runoff based on vegetation and soil moisture dynamics, while mizuRoute processes this runoff through a detailed river network. Performance metrics show strong model reliability, with R² ranging from 0.5 to 0.9, NSE between 0.6 and 0.9, and KGE from 0.5 to 0.8. The total annual average discharge for Africa is quantified at 3238.1 km³ . year-1, with contributions to various oceanic basins: 989.9 km³ . year-1 to the North Atlantic, primarily from West African rivers like the Senegal, Gambia, Volta, and Niger; 1313.7 km³ . year-1 to the South Atlantic, largely from the Congo River; 212.5 km³ . year-1 to the Mediterranean Sea, predominantly from the Nile River; and 722.0 km³ . year-1 to the Indian Ocean, with substantial inputs from rivers such as the Zambezi. This VegDischarge v1 is valuable for policymakers, stakeholders, and researchers to better understand water availability, its temporal and spatial variations, that impact water-related infrastructure planning, sustainable resource allocation, and the development of climate resilience mitigation strategies.

Scientific Data↗
USGS2018

Spatial variability in ocean-mediated growth potential is linked to Chinook salmon survival

Early ocean survival of Chinook salmon, Oncorhynchus tshawytscha, varies greatly inter-annually and may be the period during which later spawning abundance and fishery recruitment are set. Therefore, identifying environmental drivers related to early survival may inform better models for management and sustainability of salmon in a variable environment. With this in mind, our main objectives were to (a) identify regions of high temporal variability in growth potential over a 23-year time series, (b) determine whether the spatial distribution of growth potential was correlated with observed oceanographic conditions, and (c) determine whether these spatial patterns in growth potential could be used to estimate juvenile salmon survival. We applied this method to the fall run of the Central Valley Chinook salmon population, focusing on the spring and summer period after emigration into central California coastal waters. For the period from 1988 to 2010, juvenile salmon growth potential on the central California continental shelf was described by three spatial patterns. These three patterns were most correlated with upwelling, detrended sea level anomalies, and the strength of onshore/offshore currents, respectively. Using the annual strength of these three patterns, as well as the overall growth potential throughout central California coastal waters, in a generalized linear model we explained 82% of the variation in juvenile salmon survival estimates. We attributed the relationship between growth potential and survival to variability in environmental conditions experienced by juvenile salmon during their first year at sea, as well as potential shifts in predation pressure following out-migration into coastal waters.

Fisheries Oceanography↗
USGS2021

Preliminary assessment of the wave generating potential from landslides at Barry Arm, Prince William Sound, Alaska

We simulated the concurrent rapid motion of landslides on an unstable slope at Barry Arm, Alaska. Movement of landslides into the adjacent fjord displaced fjord water and generated a tsunami, which propagated out of Barry Arm. Rather than assuming an initial sea surface height, velocity, and location for the tsunami, we generated the tsunami directly using a model capable of simulating the dynamics of both water and landslide material. The fjord below most of the landslide source area was occupied by the Barry Glacier until about 2012; therefore, our direct simulation of tsunami generation by landslide motion required new topographic and bathymetric data, which was collected in 2020. The topographic data also constrained landslide geometries and volumes. We considered four scenarios based on two landslide volumes and two landslide mobilities—a more mobile, contractive landslide and a less mobile, noncontractive landslide. The larger of the two volumes is 689 × 10 6 cubic meters (m 3 )—larger than the volume estimate in a previous study—and reflects the largest plausible volume given current observational data. The considered scenario that generated the largest wave heights resulted in forecast wave heights of over 200 meters (m) in the northern part of Barry Arm, adjacent to the landslide source area and runup on the opposite fjord wall in excess of 500 m. Simulated wave heights in excess of 5 m in southern Barry Arm and in Harriman Fjord occurred within 10–15 minutes (min) of landslide motion. The simulated tsunami reached Whittier, Alaska, approximately 20 min after initial rapid landslide motion, with peak heights of just over 2 m in Passage Fjord, 500 m offshore Whittier, occurring 26 min after initial rapid motion. Time of peak wave heights was consistent with previous modeling. Although results are preliminary and can be refined with additional observations and analyses, they provide a refined assessment of the upper bound of the hazard presented by the Barry Arm landslides. The results herein support the National Oceanic and Atmospheric Administration’s National Tsunami Warning Center mission to detect, forecast, and warn for tsunamis in Alaska.

Alaska↗
USGS2013

Marine benthic habitat mapping of the West Arm, Glacier Bay National Park and Preserve, Alaska

Seafloor geology and potential benthic habitats were mapped in West Arm, Glacier Bay National Park and Preserve, Alaska, using multibeam sonar, groundtruthed observations, and geological interpretations. The West Arm of Glacier Bay is a recently deglaciated fjord system under the influence of glacial and paraglacial marine processes. High glacially derived sediment and meltwater fluxes, slope instabilities, and variable bathymetry result in a highly dynamic estuarine environment and benthic ecosystem. We characterize the fjord seafloor and potential benthic habitats using the recently developed Coastal and Marine Ecological Classification Standard (CMECS) by the National Oceanic and Atmospheric Administration (NOAA) and NatureServe. Due to the high flux of glacially sourced fines, mud is the dominant substrate within the West Arm. Water-column characteristics are addressed using a combination of CTD and circulation model results. We also present sediment accumulation data derived from differential bathymetry. These data show the West Arm is divided into two contrasting environments: a dynamic upper fjord and a relatively static lower fjord. The results of these analyses serve as a test of the CMECS classification scheme and as a baseline for ongoing and future mapping efforts and correlations between seafloor substrate, benthic habitats, and glacimarine processes.

Alaska↗
USGS2007

Hydrogeology and Simulated Ground-Water Flow in the Salt Pond Region of Southern Rhode Island

The Salt Pond region of southern Rhode Island extends from Westerly to Narragansett Bay and forms the natural boundary between the Atlantic Ocean and the shallow, highly permeable freshwater aquifer of the South Coastal Basin. Large inputs of fresh ground water coupled with the low flushing rates to the open ocean make the salt ponds particularly susceptible to eutrophication and bacterial contamination. Ground-water discharge to the salt ponds is an important though poorly quantified source of contaminants, such as dissolved nutrients. A ground-water-flow model was developed and used to delineate the watersheds to the salt ponds, including the areas that contribute ground water directly to the ponds and the areas that contribute ground water to streams that flow into ponds. The model also was used to calculate ground-water fluxes to these coastal areas for long-term average conditions. As part of the modeling analysis, adjustments were made to model input parameters to assess potential uncertainties in model-calculated watershed delineations and in ground-water discharge to the salt ponds. The results of the simulations indicate that flow to the salt ponds is affected primarily by the ease with which water is transmitted through a glacial moraine deposit near the regional ground-water divide, and by the specified recharge rate used in the model simulations. The distribution of the total freshwater flow between direct ground-water discharge and ground-water-derived surface-water (streamflow) discharge to the salt ponds is affected primarily by simulated stream characteristics, including the streambed-aquifer connection and the stream stage. The simulated position of the ground-water divide and, therefore, the model-calculated watershed delineations for the salt ponds, were affected only by changes in the transmissivity of the glacial moraine. Selected changes in other simulated hydraulic parameters had substantial effects on total freshwater discharge and the distribution of direct ground-water discharge and ground-water-derived surface-water (streamflow) discharge to the salt ponds, but still provided a reasonable match to the hydrologic data available for model calibration. To reduce the uncertainty in predictions of watershed areas and ground-water discharge to the salt ponds, additional hydrogeologic data would be required to constrain the model input parameters that have the greatest effect on the simulation results.

Scientific Investigations Report↗
USGS2021

Complex vulnerabilities of the water and aquatic carbon cycles to permafrost thaw

The spatial distribution and depth of permafrost are changing in response to warming and landscape disturbance across northern Arctic and boreal regions. This alters the infiltration, flow, surface and subsurface distribution, and hydrologic connectivity of inland waters. Such changes in the water cycle consequently alter the source, transport, and biogeochemical cycling of aquatic carbon (C), its role in the production and emission of greenhouse gases, and C delivery to inland waters and the Arctic Ocean. Responses to permafrost thaw across heterogeneous boreal landscapes will be neither spatially uniform nor synchronous, thus giving rise to expressions of low to medium confidence in predicting hydrologic and aquatic C response despite very high confidence in projections of widespread near-surface permafrost disappearance as described in the 2019 Intergovernmental Panel on Climate Change Special Report on the Ocean and Cryosphere in a Changing Climate: Polar Regions. Here, we describe the state of the science regarding mechanisms and factors that influence aquatic C and hydrologic responses to permafrost thaw. Through synthesis of recent topical field and modeling studies and evaluation of influential landscape characteristics, we present a framework for assessing vulnerabilities of northern permafrost landscapes to specific modes of thaw affecting local to regional hydrology and aquatic C biogeochemistry and transport. Lastly, we discuss scaling challenges relevant to model prediction of these impacts in heterogeneous permafrost landscapes.

Arctic↗
USGS1984

Origin and geochemistry of Cretaceous deep-sea black shales and multicolored claystones, with emphasis on Deep Sea Drilling Project Site 530, southern Angola Basin

Deep-water sedimentary sequences of mid-Cretaceous age, rich in organic carbon, have been recovered at many DSDP sites in the Atlantic Ocean. Most of these sequences have a marked cyclicity in amount of organic carbon resulting in interbedded multicolored shale, marlstone, and (or) limestone that have cycle periods of 20,000 to 100,000 years and average 40,000 to 50,000 years. These cycles may be related to some climatic control on influx of terrigenous organic matter and sediment, rates of upwelling and sea-surface production of organic matter, and preservation of organic matter related to deeper-water dissolved oxygen concentration. These variations in supply of organic matter had pronounced effects on the potential of the sediment for subsequent diagenetic changes and geochemical partitioning in adjacent beds. Many trace elements are enriched in organic-carbon-rich lithologies relative to interbedded organic-carbon-poor lithologies. Elements that are most commonly enriched are Cr, Ni, V, Cu, Zn, and Mo. The association of high traceelement concentrations with organic matter may be the result of concentration of these elements by organisms or by chemical sorption and precipitation processes under anoxic conditions. Detailed trace-element profiles from organiccarbon-rich strata at Site 530 suggest that there may be differential mobility of trace elements, with diffusion of some elements over distances of at least tens of meters. The sequence of trace-element mobility, from highest to lowest, is approximately Ba, Mn, Pb, Ni, Co, Cr, Cu, Zn, V, Cd, and Mo. Slowly deposited, oxidized clays directly overlying some black shale sequences are enriched in some metals, particularly Fe, Mn, Zn, and Cu, relative to normal pelagic clays, and this enrichment may be the result of upward migration of metals in pore waters during compaction or diffusion from the underlying black shale. Most depositional models that have been used to explain the accumulation of the organic-carbon-rich strata imply that reducing conditions in the sediments (and therefore the increased degree of preservation of organic matter) were the result of anoxic or near-anoxic conditions in oceanic bottom waters, or in a midwater oxygen-minimum zone. Evidence from several DSDP sites in the Atlantic, however, indicate that some of these middle Cretaceous "black shale" beds may be the result of variations in rate of supply of organic matter that produced anoxia or near-anoxia within midwater oxygen-minimum zones and possibly, under extreme conditions, throughout much of the bottomwater mass. Although bottom-water anoxia may have occurred during periods of organic-carbon-rich strata, it was not necessarily the only cause for accumulation of these strata. The main reason for the accumulation of organic-carbonrich strata was an increase in the relative amount of organic debris being deposited. Some of this organic debris was derived from continental-margin areas of increased production, accumulation, and preservation of organic matter from marine, terrestrial, or mixed sources and transported to slope and basinal sites by turbidity currents.

Initial Reports of the D.S.D.P.↗
USGS2013

Application of a hydrodynamic and sediment transport model for guidance of response efforts related to the Deepwater Horizon oil spill in the Northern Gulf of Mexico along the coast of Alabama and Florida

U.S. Geological Survey (USGS) scientists have provided a model-based assessment of transport and deposition of residual Deepwater Horizon oil along the shoreline within the northern Gulf of Mexico in the form of mixtures of sand and weathered oil, known as surface residual balls (SRBs). The results of this USGS research, in combination with results from other components of the overall study, will inform operational decisionmaking. The results will provide guidance for response activities and data collection needs during future oil spills. In May 2012 the U.S. Coast Guard, acting as the Deepwater Horizon Federal on-scene coordinator, chartered an operational science advisory team to provide a science-based review of data collected and to conduct additional directed studies and sampling. The goal was to characterize typical shoreline profiles and morphology in the northern Gulf of Mexico to identify likely sources of residual oil and to evaluate mechanisms whereby reoiling phenomena may be occurring (for example, burial and exhumation and alongshore transport). A steering committee cochaired by British Petroleum Corporation (BP) and the National Oceanic and Atmospheric Administration (NOAA) is overseeing the project and includes State on-scene coordinators from four States (Alabama, Florida, Louisiana, and Mississippi), trustees of the U.S. Department of the Interior (DOI), and representatives from the U.S. Coast Guard. This report presents the results of hydrodynamic and sediment transport models and developed techniques for analyzing potential SRB movement and burial and exhumation along the coastline of Alabama and Florida. Results from these modeling efforts are being used to explain the complexity of reoiling in the nearshore environment and to broaden consideration of the different scenarios and difficulties that are being faced in identifying and removing residual oil. For instance, modeling results suggest that larger SRBs are not, under the most commonly observed low-energy wave conditions, likely to move very far alongshore. This finding suggests that SRBs from one source location may not (outside of storm conditions) be redistributed to other up or down coast locations. This information can guide operational response decisions. In addition, because SRBs are less mobile compared with sand, they are likely to become buried and unburied under normal sand transport processes thereby lengthening the time SRBs may take to move onshore. The rate of onshore movement was not specifically addressed by this study, yet the results resolve the cross-shore domain and cross-shore variations in alongshore transport that are relevant to achieving the primary objectives. Furthermore, during infrequent events (for example, winter storms and severe meteorological events such as Hurricane Isaac of August 2012), energy is shown to be sufficient to move a greater range of SRB sizes and potentially expose and break up submerged oil mats. When SRBs do move alongshore, the models indicate that there are regions that are more conducive to accumulation of SRB material than others. Accumulation can occur where there are reversals and decelerations in alongshore currents and where forces created by shear stress drops below critical thresholds to maintain or initiate SRB movement. In addition, flow and SRB mobility patterns around inlets indicate patterns in hydrodynamic forces that influence redistribution of SRBs and the surface oil that mixed with sediment to form oil mats in the first place.

Alabama;Florida↗
USGS2017

Lithospheric density structure beneath the Tarim basin and surroundings, northwestern China, from the joint inversion of gravity and topography

Intraplate strain generally focuses in discrete zones, but despite the profound impact of this partitioning on global tectonics, geodynamics, and seismic hazard, the processes by which deformation becomes localized are not well understood. Such heterogeneous intraplate strain is exemplified in central Asia, where the Indo-Eurasian collision has caused widespread deformation while the Tarim block has experienced minimal Cenozoic shortening. The apparent stability of Tarim may arise either because strain is dominantly accommodated by pre-existing faults in the continental suture zones that bound it—essentially discretizing Eurasia into microplates—or because the lithospheric-scale strength (i.e., viscosity) of the Tarim block is greater than its surroundings. Here, we jointly analyze seismic velocity, gravity, topography, and temperature to develop a 3-D density model of the crust and upper mantle in this region. The Tarim crust is characterized by high density, v s , v p , and v p /v s , consistent with a dominantly mafic composition and with the presence of an oceanic plateau beneath Tarim. Low-density but high-velocity mantle lithosphere beneath southern (southwestern) Tarim underlies a suite of Permian plume-related mafic intrusions and A-type granites sourced in previously depleted mantle lithosphere; we posit that this region was further depleted, dehydrated, and strengthened by Permian plume magmatism. The actively deforming western and southern margins of Tarim—the Tien Shan, Kunlun Shan, and Altyn Tagh fault—are underlain by buoyant upper mantle with low velocity; we hypothesize that this material has been hydrated by mantle-derived fluids that have preferentially migrated along Paleozoic continental sutures. Such hydrous material should be weak, and herein strain focuses there because of lithospheric-scale variations in rheology rather than the pre-existence of faults in the brittle crust. Thus this world-class example of strain partitioning arises not simply from the pre-existence of brittle faults but from the thermo-chemical and therefore rheological variations inherited from prior tectonism.

Earth and Planetary Science Letters↗
USGS2018

Stress rotation across the Cascadia megathrust requires a weak subduction plate boundary at seismogenic depths

The Mendocino Triple Junction region is the most seismically active part of the Cascadia Subduction Zone. The northward moving Pacific plate collides with the subducting Gorda plate causing intense internal deformation within it. Here we show that the stress field rotates rapidly with depth across the thrust interface from a strike-slip regime within the subducting plate, reflecting the Pacific plate collision, to a thrust regime in the overriding plate. We utilize a dense focal mechanism dataset, including observations from the Cascadia Initiative ocean bottom seismograph experiment, to constrain the stress orientations. To quantify the implications of this rotation for the strength of the plate boundary, we designed an inversion that solves for the absolute stress tensors in a three-layer model subject to assumptions about the strength of the subducting mantle. Our results indicate that the shear stress on the plate boundary fault is likely no more than about ∼50 MPa at ∼20 km depth. Regardless of the assumed mantle strength, we infer a relatively weak megathrust fault with an effective friction coefficient of ∼0 to 0.2 at seismogenic depths. Such a low value for the effective friction coefficient requires a combination of high fluid pressures and/or fault-zone minerals with low inherent friction in the region where a great earthquake is expected in Cascadia.

California↗
USGS2008

Constraints on the Pleistocene chronology of sediments from the Lomonosov Ridge

Despite its importance in the global climate system, age-calibrated marine geologic records reflecting the evolultion of glacial cycles through the Pleistocene are largely absent from the central Arctic Ocean. This is especially true for sediments older than 200 ka. Three sites cored during the Integrated Ocean Drilling Program's Expedition 302, the Arctic Coring Expedition (ACEX), provide a 27 m continuous sedimentary section from the Lomonosov Ridge in the central Arctic Ocean. Two key biostratigraphic datums and constraints from the magnetic inclination data are used to anchor the chronology of these sediments back to the base of the Cobb Mountain subchron (1215 ka). Beyond 1215 ka, two best fitting geomagnetic models are used to investigate the nature of cyclostratigraphic change. Within this chronology we show that bulk and mineral magnetic properties of the sediments vary on predicted Milankovitch frequencies. These cyclic variations record "glacial" and "interglacial" modes of sediment deposition on the Lomonosov Ridge as evident in studies of ice-rafted debris and stable isotopic and faunal assemblages for the last two glacial cycles and were used to tune the age model. Potential errors, which largely arise from uncertainties in the nature of downhole paleomagnetic variability, and the choice of a tuning target are handled by defining an error envelope that is based on the best fitting cyclostratigraphic and geomagnetic solutions. Copyright 2008 by the American Geophysical Union.

Paleoceanography↗