Geology Reports⌕ Search

SEARCH · Geology Reports

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 775 records · Page 43Linked to original sources
USGS1969

Heat flow in the Arctic

Defines heat flow as the flux at the earth's solid surface of heat conducted from the interior; the heat-flow-unit (hfu) is on the order of 1-millionth calorie through each sq cm of the surface/sec, which is enough to melt a 4-mm layer of ice over the earth's surface/yr. Earth heat originates from radioactive decay of U, Th and K in the crust and mantle. Although land heat-flow measurements in the Arctic are too few for regional interpretation, those from Cape Thompson, Barrow and Cape Simpson, Northern Alaska are discussed and figured to show what they contribute to understanding of permafrost, climatic change and shoreline movements. Measuring thermal conductivity and gradient is much simpler in ocean basins than on land. Locations of such measurements are mapped, the results for the Alaskan quadrant in more detail. The sharp change in heat flow at the edge of the Alpha Cordillera, shown in a geothermal model, suggests that this feature is a huge accumulation of basalt, rather than mantle material or remnant of a foundering continent as previously postulated. Future Arctic heat flow studies are discussed.

Arctic↗
USGS2021

Watersheds and drainage networks

This topic is an overview of basic concepts about how the distribution of water on the Earth, with specific regard to watersheds, stream and river networks, and waterbodies are represented by geographic data. The flowing and non-flowing bodies of water on the earth’s surface vary in extent largely due to seasonal and annual changes in climate and precipitation. Consequently, modeling the detailed representation of surface water using geographic information is important. The area of land that collects surface runoff and other flowing water and drains to a common outlet location defines a watershed. Terrain and surface features can be naturally divided into watersheds of various sizes. Drainage networks are important data structures for modeling the distribution and movement of surface water over the terrain. Numerous tools and methods exist to extract drainage networks and watersheds from digital elevation models (DEMs). The cartographic representations of surface water are referred to as hydrographic features and consist of a snapshot at a specific time. Hydrographic features can be assigned general feature types, such as lake, pond, river, and ocean. Hydrographic features can be stored, maintained, and distributed for use through vector geospatial databases, such as the National Hydrography Dataset (NHD) for the United States.

Book chapter↗
USGS2017

Deglacial sea level history of the East Siberian Sea and Chukchi Sea margins

Deglacial (12.8–10.7 ka) sea level history on the East Siberian continental shelf and upper continental slope was reconstructed using new geophysical records and sediment cores taken during Leg 2 of the 2014 SWERUS-C3 expedition. The focus of this study is two cores from Herald Canyon, piston core SWERUS-L2-4-PC1 (4-PC1) and multicore SWERUS-L2-4-MC1 (4-MC1), and a gravity core from an East Siberian Sea transect, SWERUS-L2-20-GC1 (20-GC1). Cores 4-PC1 and 20-GC were taken at 120 and 115 m of modern water depth, respectively, only a few meters above the global last glacial maximum (LGM; ∼ 24 kiloannum or ka) minimum sea level of ∼ 125–130 meters below sea level (m b.s.l.). Using calibrated radiocarbon ages mainly on molluscs for chronology and the ecology of benthic foraminifera and ostracode species to estimate paleodepths, the data reveal a dominance of river-proximal species during the early part of the Younger Dryas event (YD, Greenland Stadial GS-1) followed by a rise in river-intermediate species in the late Younger Dryas or the early Holocene (Preboreal) period. A rapid relative sea level rise beginning at roughly 11.4 to 10.8 ka ( ∼ 400 cm of core depth) is indicated by a sharp faunal change and unconformity or condensed zone of sedimentation. Regional sea level at this time was about 108 m b.s.l. at the 4-PC1 site and 102 m b.s.l. at 20-GC1. Regional sea level near the end of the YD was up to 42–47 m lower than predicted by geophysical models corrected for glacio-isostatic adjustment. This discrepancy could be explained by delayed isostatic adjustment caused by a greater volume and/or geographical extent of glacial-age land ice and/or ice shelves in the western Arctic Ocean and adjacent Siberian land areas.

Chukchi Sea, East Siberian Sea↗
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↗
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↗