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.

SourceResearch connectionFrom yearThrough year

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

At least 1,369 records · Page 76Linked to original sources

USGS1998

Structural and kinematic evolution of the Yukon-Tanana upland tectonites, east-central Alaska: A record of late Paleozoic to Mesozoic crustal assembly

The Yukon-Tanana terrane, the largest tectonostratigraphic terrane in the northern North American Cordillera, is polygenetic and not a single terrane. Lineated and foliated (L-S) tectonites, which characterize the Yukon-Tanana terrane, record multiple deformations and formed at different times. We document the polyphase history recorded by L-S tectonites within the Yukon-Tanana upland, east-central Alaska. These upland tectonites compose a heterogeneous assemblage of deformed igneous and metamorphic rocks that form the Alaskan part of what has been called the Yukon-Tanana composite terrane. We build on previous kinematic data and establish the three-dimensional architecture of the upland tectonites through kinematic and structural analysis of more than 250 oriented samples, including quartz c-axis fabric analysis of 39 samples. Through this study we distinguish allochthonous tectonites from parautochthonous tectonites within the Yukon-Tanana upland. The upland tectonites define a regionally coherent stacking order: from bottom to top, they are lower plate North American parautochthonous attenuated continental margin; continentally derived marginal-basin strata; and upper plate ocean-basin and island-arc rocks, including some continental basement rocks. We delineate three major deformation events in time, space, and structural level across the upland from the United States-Canada border to Fairbanks, Alaska: (1) pre-Early Jurassic (>212 Ma) northeast-directed, apparent margin-normal contraction that affected oceanic rocks; (2) late Early to early Middle Jurassic (>188-185 Ma) northwest-directed, apparent margin-parallel contraction and imbrication that resulted in juxtaposition of the allochthonous tectonites with parautochthonous continental rocks; and (3) Early Cretaceous (135-110 Ma) southeast-directed crustal extension that resulted in exposure of the structurally deepest, parautochthonous continental rocks. The oldest event represents deformation within a west-dipping (present coordinates) Permian-Triassic subduction zone. The second event records Early to Middle Jurassic collision of the arc and subduction complex with North American crust, and the third event reflects mid-Cretaceous southeast-directed crustal extension. Events one and two can be recognized and correlated through southern Yukon, even though this region was affected by mid-Cretaceous dextral shear along steep northwest-striking faults. Our data support a model of crustal assembly originally proposed by D. Tempelman-Kluit in which previously deformed allochthonous rocks were thrust over parautochthonous rocks of the attenuated North American margin in Middle Jurassic time. Approximately 50 m.y. after tectonic accretion, east-central Alaska was dissected by crustal extension, exposing overthrust parautochthonous strata.

Alaska↗

USGS2016

Sedimentologic characteristics of recent washover deposits from Assateague Island, Maryland

The U.S. Geological Survey has a long history of responding to and documenting the impacts of storms along the Nation’s coasts and incorporating these data into storm impact and coastal change vulnerability assessments. Although physical changes caused by tropical and extratropical storms to the sandy beaches and dunes fronting barrier islands are generally well documented, the interaction between sandy shoreline erosion and overwash with the back-barrier wetland and estuarine environments is poorly constrained. The goal of the Barrier Island and Estuarine Wetland Physical Change Assessment project is to integrate a wetland-change assessment with existing coastal-change assessments for the adjacent sandy dunes and beaches, initially focusing on Assateague Island along the Maryland and Virginia coastline. Assateague Island was impacted by waves and storm surge associated with the passage of Hurricane Sandy in October 2012, causing erosion and overwash along the ocean-facing sandy shoreline as well as erosion and overwash deposition in the back-barrier and estuarine bay environments. This report describes sediment data collected using sand augers in active overwash zones on Assateague Island in Maryland. Samples were collected by the U.S. Geological Survey (USGS) during two surveys in March/April and October 2014 (USGS Field Activity Numbers [FAN] 2014-301-FA and 2014-322-FA, respectively). The physical characteristics (for example, sediment texture or bedding structure) of and spatial differences among these deposits will provide information about overwash processes and sediment transport from the sandy barrier-island reaches to the back-barrier environments. Metrics derived from these data, such as mean grain size or deposit thicknesses, can be used to ground-truth remote sensing and geophysical data and can also be incorporated into sediment transport models. Data products, including sample location tables, descriptive core logs, core photographs and x-radiographs, the results of sediment grain-size analyses, and Geographic Information System (GIS) data files with accompanying formal Federal Geographic Data Committee (FGDC) metadata can be downloaded from the Data Downloads page.

Maryland↗

USGS2001

Wavefield properties of a shallow long-period event and tremor at Kilauea Volcano, Hawaii

The wavefields of tremor and a long-period (LP) event associated with the ongoing eruptive activity at Kilauea Volcano, Hawaii, are investigated using a combination of dense small-aperture (300 m) and sparse large-aperture (5 km) arrays deployed in the vicinity of the summit caldera. Measurements of azimuth and slowness for tremor recorded on the small-aperture array indicate a bimodal nature of the observed wavefield. At frequencies below 2 Hz, the wavefield is dominated by body waves impinging the array with steep incidence. These arrivals are attributed to the oceanic microseismic noise. In the 2-6 Hz band, the wavefield is dominated by waves propagating from sources located at shallow depths (<1 km) beneath the eastern edge of the Halemaumau pit crater. The hypocenter of the LP event, determined from frequency-slowness analyses combined with phase picks, appears to be located close to the source of tremor but at a shallower depth (<0.1 km). The wavefields of tremor and LP event are characterized by a complex composition of body and surface waves, whose propagation and polarization properties are strongly affected by topographic and structural features in the summit caldera region. Analyses of the directional properties of the wavefield in the 2-6 Hz band point to the directions of main scattering sources, which are consistent with pronounced velocity contrasts imaged in a high-resolution three-dimensional velocity model of the caldera region. The frequency and Q of the dominant peak observed in the spectra of the LP event may be explained as the dominant oscillation mode of a crack with scale length 20-100 m and aperture of a few centimeters filled with bubbly water. The mechanism driving the shallow tremor appears to be consistent with a sustained excitation originating in the oscillations of a bubbly cloud resulting from vesiculation and degassing in the magma. ?? 2001 Elsevier Science B.V. All rights reserved.

Journal of Volcanology and Geothermal Research↗

USGS2017

Population genetic structure and gene flow of Adélie penguins (Pygoscelis adeliae) breeding throughout the western Antarctic Peninsula

Adélie penguins ( Pygoscelis adeliae ) are responding to ocean–climate variability throughout the marine ecosystem of the western Antarctic Peninsula (WAP) where some breeding colonies have declined by 80%. Nuclear and mitochondrial DNA (mtDNA) markers were used to understand historical population genetic structure and gene flow given relatively recent and continuing reductions in sea ice habitats and changes in numbers of breeding adults at colonies throughout the WAP. Genetic diversity, spatial genetic structure, genetic signatures of fluctuations in population demography and gene flow were assessed in four regional Adélie penguin colonies. The analyses indicated little genetic structure overall based on bi-parentally inherited microsatellite markers ( F ST =-0.006–0.004). No significant variance was observed in overall haplotype frequency (mtDNA Φ ST =0.017; P =0.112). Some comparisons with Charcot Island were significant, suggestive of female-biased philopatry. Estimates of gene flow based on a two-population coalescent model were asymmetrical from the species’ regional core to its northern range. Breeding Adélie penguins of the WAP are a panmictic population and hold adequate genetic diversity and dispersal capacity to be resilient to environmental change.

Antarctic Science↗

USGS2021

Distinguishing between regression model fits to global mean sea level reconstructions

Global mean sea level (GMSL) has been rising since the last century, posing a serious challenge for the coastal areas. A variety of regression models have been utilized for determining GMSL rise over the past one hundred years, resulting in a large spread of sea level rise rates and multidecadal variations. In this study, we develop a new nonparametric noise model that is data-dependent and considers overfitting due to regression. The noise model is used to determine whether one regression model has significantly better skill than others over the period 1900–2010. The choices of background noise and GMSL reconstruction influence whether two sea level models can be statistically distinguished. With our new nonparametric noise spectra, the differences of model skills in explaining sea level variance are significant only in 34% of model comparisons. However, stepwise trends with three inflection points are significantly more skillful than the linear, quadratic, or exponential trend for most GMSL reconstructions, suggesting the importance of multidecadal variability of sea level rise in the twentieth century. Nevertheless, stepwise trend models cannot be distinguished from models with a long-term harmonic oscillation, indicating that the shape of multidecadal variability is not conclusive. The multidecadal variability is also significant in the steric and barystatic sea level contributions and is related to both natural and anthropogenic forcings. GMSL predictions based on regression fits in the twentieth century underestimate the sea level rise rate over the period 2011–2020 because the sea level acceleration in the recent decade (2011–2020) is not well represented.

Journal of Geophysical Research- Oceans↗

USGS2004

Across-shelf sediment transport since the Last Glacial Maximum, southern California margin

Correlation of continental shelf-slope stratigraphy in Santa Monica Bay (southern California) with Ocean Drilling Program records for nearby slope-basin sites has illuminated the timing and scale of terrigenous sediment dispersal on margin since the Last Glacial Maximum (LGM). Marine flooding surfaces preserved in a transgressive sequence on the Santa Monica Shelf provide a key link between base-level elevation and sediment transport across shelf. Sediment-accumulation rates at slope-basin sites were maximal ca. 15-10 ka, well after the LGM, decreased during the 12-9 ka transition from fluvial-estuarine to fully marine conditions on the shelf, and decelerated throughout the Holocene to 30%-75% of their values at the LGM. The deceleration is interpreted to manifest a landward shift in the margin depocenter with the onset of transgressive sedimentation beginning when sea level surmounted the shelf edge ca. 13 ka, as predicted by sequence-stratigraphic models. However, the records make clear that factors other than base level modulated slope-basin accumulation rates during the deglaciation. ?? 2004 Geological Society of America.

Geology↗

USGS2007

Highstand fans in the California borderland: The overlooked deep-water depositional systems

Contrary to widely used sequence-stratigraphic models, lowstand fans are only part of the turbidite depositional record; our analysis reveals that a comparable volume of coarse-grained sediment has been deposited in California borderland deep-water basins regardless of sea level. Sedimentation rates and periods of active sediment transport have been determined for deep-water canyon-channel systems contributing to the southeastern Gulf of Santa Catalina and San Diego Trough since 40 ka using an extensive grid of high-resolution and deep-penetration seismic-reflection data. A regional seismic-reflection horizon (40 ka) has been correlated across the study area using radiocarbon age dates from the Mohole borehole and U.S. Geological Survey piston cores. This study focused on the submarine fans fed by the Oceanside, Carlsbad, and La Jolla Canyons, all of which head within the length of the Ocean-side littoral cell. The Oceanside Canyon–channel system was active from 45 to 13 ka, and the Carlsbad system was active from 50 (or earlier) to 10 ka. The La Jolla system was active over two periods, from 50 (or earlier) to 40 ka, and from 13 ka to the present. One or more of these canyon-channel systems have been active regardless of sea level. During sea-level fluctuation, shelf width between the canyon head and the littoral zone is the primary control on canyon-channel system activity. Highstand fan deposition occurs when a majority of the sediment within the Oceanside littoral cell is intercepted by one of the canyon heads, currently La Jolla Canyon. Since 40 ka, the sedimentation rate on the La Jolla highstand fan has been >2 times the combined rates on the Oceanside and Carlsbad lowstand fans.

California↗

USGS2003

Maps Showing Inundation Depths, Ice-Rafted Erratics, and Sedimentary Facies of Late Pleistocene Missoula Floods in the Willamette Valley, Oregon

Glacial Lake Missoula, impounded by the Purcell Trench lobe of the late Pleistocene Cordilleran Icesheet, repeatedly breached its ice dam, sending floods as large as 2,500 cubic kilometers racing across the Channeled Scabland and down the Columbia River valley to the Pacific Ocean. Peak discharges for some floods exceeded 20 million cubic meters per second. At valley constrictions along the flood route, floodwaters temporarily ponded behind each narrow zone. One such constriction at Kalama Gap-northwest of Portland-backed water 120-150 meters high in the Portland basin, and backflooded 200 km south into Willamette Valley. Dozens of floods backed up into the Willamette Valley, eroding 'scabland' channels, and depositing giant boulder gravel bars in areas of vigorous currents as well as bedded flood sand and silt in backwater areas. Also, large chunks of ice entrained from the breached glacier dam rafted hundreds of 'erratic' rocks, leaving them scattered among the flanking foothills and valley bottom. From several sources and our own mapping, we have compiled information on many of these features and depict them on physiographic maps derived from digital elevation models of the Portland Basin and Willamette Valley. These maps show maximum flood inundation levels, inundation levels associated with stratigraphic evidence of repeated floodings, distribution of flood deposits, and sites of ice-rafted erratics. Accompanying these maps, a database lists locations, elevations, and descriptions of approximately 400 ice-rafted erratics-most compiled from early 20th-century maps and notes of A.M. Piper and I.S. Allison.

Open-File Report↗

USGS2023

Data-limited fishery assessment methods shed light on the exploitation history and population dynamics of Endangered Species Act-listed Yelloweye Rockfish in Puget Sound, Washington

Objective The distinct population segment (DPS) of Yelloweye Rockfish Sebastes ruberrimus inhabiting the Puget Sound/Georgia Basin was listed under the Endangered Species Act (ESA) in 2010, and a formal recovery plan for the DPS was published by National Oceanic and Atmospheric Administration Fisheries in 2017. In this recovery plan, the biological criteria for delisting or downlisting were specified as certain levels of spawning potential ratio (SPR), a commonly used metric of equilibrium stock status for commercially exploited fishes. Although this metric can be estimated from length compositions, the combination of length data with a catch history (which was not previously available for this DPS) improves our understanding of population dynamics over time and allows us to estimate a different measure of stock status, relative (to unfished) spawning stock biomass (SSB), rather than only SPR. Methods To estimate relative SSB and reconstruct the historical dynamics of this DPS, we reconstructed the catch history from fisheries records, collated length data from historical and contemporary hook-and-line surveys, and fitted a data-limited version of a statistical catch-at-age model. Result Despite a high level of uncertainty, we estimated that Yelloweye Rockfish in Puget Sound are above 25% of unfished biomass (a reference point detailed in the recovery criteria) under the assumption of deterministic recruitment, presenting the first direct estimates of Yelloweye Rockfish population status in Puget Sound. Conclusion However, as informed by recent genetic studies, the DPS boundaries of ESA-listed Yelloweye Rockfish extend from South Puget Sound to Queen Charlotte Strait in British Columbia. The Canadian portion of this population is managed separately and is currently estimated to be at 32% of unfished biomass (95% quantiles = 15%–68%). Thus, the disjunction between the biological boundaries of the population and the jurisdictional boundaries between Canada and the United States presents an additional source of uncertainty in assessing recovery that must be addressed to achieve DPS-wide recovery goals.

Washington↗

USGS1989

Salinity variations in submarine hydrothermal systems by layered double-diffusive convection

Various mechanisms have been proposed to explain the salinity variations in vent fluids of seafloor geothermal systems. New experiments reacting diabase and evolved seawater were carried out to reproduce earlier published observations of Cl depletions attributed to formation of an ephemeral Cl-bearing mineral. The absence of any Cl depletions in the present study suggests that the formation of Cl-bearing minerals is not sufficiently widespread to account for the observed salinity variations in the vent fluids. A re-evaluation of both field and laboratory evidence has led to a new model for subseafloor circulation that accounts for salinity variations as well as other chemical and mineralogic observations. In place of a simple single-pass convection system, we propose that the seafloor systems consist of two vertically nested convection cells in which a brine layer at depth heats and drives an overlying seawater cell. Such layering of salinities, a process known in fluid mechanics as double-diffusive convection, is an expected result when convection is induced in saline fluids. The process provides for stable high-temperature heat transfer upward from the cracking front adjacent to the magma, and for limited chemical exchange of the brine with the overlying seawater to explain salinity variations and high metal contents in the vent fluids. The brine also provides an effective medium to produce the secondary mineral assemblages observed in rocks from the mid-ocean ridges and ophiolites unsuccessfully produced in laboratory studies using seawater. The brine originates from the two-phase separation of seawater during magmatic/tectonic events and accumulates and remains relatively stable in the region immediately above the magma chamber.

Journal of Geology↗

USGS2009

Geologic Map of MTM 35337, 40337, and 45337 Quadrangles, Deuteronilus Mensae Region of Mars

Deuteronilus Mensae, first defined as an albedo feature at lat 35.0 deg N., long 5.0 deg E., by U.S. Geological Survey (USGS) and International Astronomical Union (IAU) nomenclature, is a gradational zone along the dichotomy boundary in the northern mid-latitudes of Mars. The boundary in this location includes the transition from the rugged cratered highlands of Arabia Terra to the northern lowland plains of Acidalia Planitia. Within Deuteronilus Mensae, polygonal mesas are prominent along with features diagnostic of Martian fretted terrain, including lobate debris aprons, lineated valley fill, and concentric crater fill. Lobate debris aprons, as well as the valley and crater fill deposits, are geomorphic indicators of ground ice, and their concentration in Deuteronilus Mensae is of great interest because of their potential association with Martian climate change. The paucity of impact craters on the surfaces of debris aprons and the presence of ice-cemented mantle material imply young (for example, Amazonian) surface ages that are consistent with recent climate change in this region of Mars. North of Deuteronilus Mensae are the northern lowlands, a potential depositional sink that may have had large standing bodies of water or an ocean in the past. The northern lowlands have elevations that are several kilometers below the ancient cratered highlands with significantly younger surface ages. The morphologic and topographic characteristics of the Deuteronilus Mensae region record a diverse geologic history, including significant modification of the ancient highland plateau and resurfacing of low-lying regions. Previous studies of this region have interpreted a complex array of geologic processes, including eolian, fluvial and glacial activity, coastal erosion, marine deposition, mass wasting, tectonic faulting, effusive volcanism, and hydrovolcanism. The origin and age of the Martian crustal dichotomy boundary are fundamental questions that remain unresolved at the present time. Several scenarios for its formation, including single and multiple large impact events, have been proposed and debated in the literature. Endogenic processes whereby crust is thinned by internal mantle convection and tectonic processes have also been proposed. Planetary accretion models and isotopic data from Martian meteorites suggest that the crust formed very early in Martian history. Using populations of quasi-circular depressions extracted from the topography of Mars, other studies suggest that the age difference between the highlands and lowlands could be ~100 m.y.. Furthermore, understanding the origin and age of the dichotomy boundary has been made more complicated due to significant erosion and deposition that have modified the boundary and its adjacent regions. The resulting diversity of terrains and features is likely a combined result of ancient and recent events. Detailed geologic analyses of dichotomy boundary zones are important for understanding the spatial and temporal variations in highland evolution. This information, and comparisons to other highland regions, can help elucidate the scale of potential environmental changes. Previous geomorphic and geologic mapping investigations of the Deuteronilus Mensae region have been completed at local to global scales. The regional geology was first mapped by Lucchitta (1978) at 1:5,000,000 scale using Mariner 9 data. This study concluded that high crater flux early in Martian history formed overlapping craters and basins that were later filled by voluminous lava flows that buried the impacted surface, creating the highlands. After this period of heavy bombardment, fluvial erosion of the highlands formed the canyons and valleys, followed by dissection that created the small mesas and buttes, and later, formation of the steep escarpment marking the present-day northern highland margin. After valley dissection, mass wasting and eolian processes caused lateral retreat of mesas and buttes

Scientific Investigations Map↗

USGS2024

Oblique contraction along the fastest ocean-continent transform plate boundary focuses rock uplift west of the Fairweather fault, southeast Alaska

Contraction along the Yakutat–North America plate boundary drives 4.6–9.0 mm/year Holocene rock uplift rates along Earth's fastest slipping (≥49 mm/year) ocean–continent transform fault, the Fairweather Fault. Between Icy Point and Lituya Bay, the near-vertical Fairweather fault focuses rock uplift and rapid right-lateral slip by accommodating both vertical and fault-parallel strain during oblique-slip and separate, predominantly strike-slip ruptures. Unusually high uplift rates, indicated by radiocarbon and luminescence dating, result from a 10-km-wide, asymmetric, positive flower structure along a 20°, ∼30-km-long restraining double bend in the Fairweather fault. The principal reverse fault in the flower structure, the offshore, blind Icy Point–Lituya Bay fault, ruptures no more than every 460–1040 years evidenced by uplifted Holocene shorelines. Maximum 3–5 m coseismic uplifts imply 3.1–10 m dip slip per event and earthquake magnitudes of M w 7.0–7.5. The Yakutat block collides obliquely into North America, and our model entails oblique slip on the Fairweather fault with and without corupture on the reverse fault. Oblique slip is evident by vertically offset (>25 m) fluvial and marine terraces and by the primary Fairweather fault strand that strikes >20° to the west of plate-boundary motion.

Alaska↗

USGS1986

Carbon isotope systematics of a mantle "hotspot": A comparison of Loihi Seamount and MORB glasses

The carbon isotope geochemistry of glasses from Loihi Seamount has been compared with that of MORB glasses. Stepped heating shows two carbon components in both sample suites: (1) isotopically light carbon ( avg. δ 13 C = −26.3‰ ) released < 600°C, ascribed to surficial contamination, and (2) isotopically heavy carbon released > 600°C, regarded as indigenous. The high-temperature component in MORB samples varied from 52 to 169 ppm C, average δ 13 C = −6.6‰ , consistent with previous studies (overall MORD average δ 13 C = −6.4 ± 0.9‰ ), and new results for Indian Ocean glasses are similar to Atlantic and Pacific Ocean samples. Carbon release profiles produced by stepped heating may be typical of locality, but there are no significant differences in δ 13 C values between MORB samples from different areas. Lower yields (17–110 ppm C) correlated with depth in the Loihi samples suggest that they are partially degassed. This degassing has not affected δ 13 C values significantly (avg. −5.8‰). Loihi tholeiites have higher δ 13 C (avg. −5.6‰) than the alkali basalts (avg. −7.1‰). Carbon abundances correlate well with He concentration data. Comparison of the δ 13 C values with trace element and He, Sr, Nd, and Pb isotope data from the literature suggests that the Loihi samples with highest δ 13 C have high 3 He/ 4 He and possibly the least depleted 143 Nd/ 144 Nd and 87 Sr/ 86 Sr. The carbon isotope data are consistent with previous models for Loihi involving several mantle sources, lithospheric contamination, and mixing. The slightly higher δ 13 C of Loihi tholeiites suggests that the undegassed “plume” component manifested by high 3 He/ 4 He values might have δ 13 C about 1‰ higher than the MORB average.

Earth and Planetary Science Letters↗

USGS2012

Effects of smectite to illite transformation on the frictional strength and sliding stability of intact marine mudstones

At subduction zones, earthquake nucleation and coseismic slip occur only within a limited depth range, known as the “seismogenic zone”. One leading hypothesis for the upper aseismic-seismic transition is that transformation of smectite to illite at ∼100–150°C triggers a change from rate-strengthening frictional behavior that allows only stable sliding, to rate weakening behavior considered a prerequisite for unstable slip. Previous studies on powdered gouges have shown that changes in clay mineralogy alone are unlikely to control this transition, but associated fabric and cementation developed during diagenesis remain possible candidates. We conducted shearing experiments designed specifically to evaluate this hypothesis, by using intact wafers of mudstone from Ocean Drilling Program Site 1174, offshore SW Japan, which have undergone progressive smectite transformation in situ. We sheared specimens along a sawcut in a triaxial configuration, oriented parallel to bedding, at normal stresses of ∼20–150 MPa and a pore pressure of 1 MPa. During shearing, we conducted velocity-stepping tests to measure the friction rate parameter (a-b). Friction coefficient ranges from 0.28–0.40 and values of (a-b) are uniformly positive; both are independent of clay transformation progress. Our work represents the most direct and comprehensive test of the clay transformation hypothesis to date, and suggests that neither illitization, nor accompanying fabric development and cementation, trigger a transition to unstable frictional behavior. We suggest that strain localization, in combination with precipitation of calcite and quartz, is a viable alternative that is consistent with both field observations and recent conceptual models of a heterogeneous seismogenic zone.

Philippine Sea Plate;Eurasian Plate↗

USGS2003

Sequential megafaunal collapse in the North Pacific Ocean: An ongoing legacy of industrial whaling?

Populations of seals, sea lions, and sea otters have sequentially collapsed over large areas of the northern North Pacific Ocean and southern Bering Sea during the last several decades. A bottom-up nutritional limitation mechanism induced by physical oceanographic change or competition with fisheries was long thought to be largely responsible for these declines. The current weight of evidence is more consistent with top-down forcing. Increased predation by killer whales probably drove the sea otter collapse and may have been responsible for the earlier pinniped declines as well. We propose that decimation of the great whales by post-World War II industrial whaling caused the great whales' foremost natural predators, killer whales, to begin feeding more intensively on the smaller marine mammals, thus "fishing-down" this element of the marine food web. The timing of these events, information on the abundance, diet, and foraging behavior of both predators and prey, and feasibility analyses based on demographic and energetic modeling are all consistent with this hypothesis.

Proceedings of the National Academy of Sciences of↗

USGS2009

An evaluation of density-dependent and density-independent influences on population growth rates in Weddell seals

Much of the existing literature that evaluates the roles of density-dependent and density-independent factors on population dynamics has been called into question in recent years because measurement errors were not properly dealt with in analyses. Using state-space models to account for measurement errors, we evaluated a set of competing models for a 22-year time series of mark-resight estimates of abundance for a breeding population of female Weddell seals (Leptonychotes weddellii) studied in Erebus Bay, Antarctica. We tested for evidence of direct density dependence in growth rates and evaluated whether equilibrium population size was related to seasonal sea-ice extent and the Southern Oscillation Index (SOI). We found strong evidence of negative density dependence in annual growth rates for a population whose estimated size ranged from 438 to 623 females during the study. Based on Bayes factors, a density-dependence-only model was favored over models that also included en! vironmental covariates. According to the favored model, the population had a stationary distribution with a mean of 497 females (SD = 60.5), an expected growth rate of 1.10 (95% credible interval 1.08-1.15) when population size was 441 females, and a rate of 0.90 (95% credible interval 0.87-0.93) for a population of 553 females. A model including effects of SOI did receive some support and indicated a positive relationship between SOI and population size. However, effects of SOI were not large, and including the effect did not greatly reduce our estimate of process variation. We speculate that direct density dependence occurred because rates of adult survival, breeding, and temporary emigration were affected by limitations on per capita food resources and space for parturition and pup-rearing. To improve understanding of the relative roles of various demographic components and their associated vital rates to population growth rate, mark-recapture methods can be applied that incorporate both environmental covariates and the seal abundance estimates that were developed here. An improved understanding of why vital rates change with changing population abundance will only come as we develop a better understanding of the processes affecting marine food resources in the Southern Ocean.

Ecology↗

USGS2017

Overview of the U.S. Nuclear Regulatory Commission collaborative research program to assess tsunami hazard for nuclear power plants on the Atlantic and Gulf Coasts

In response to the 2004 Indian Ocean Tsunami, the United States Nuclear Regulatory Commission (US NRC) initiated a long-term research program to improve understanding of tsunami hazard levels for nuclear facilities in the United States. For this effort, the US NRC organized a collaborative research program with the United States Geological Survey (USGS) and the National Oceanic and Atmospheric Administration (NOAA) with a goal of assessing tsunami hazard on the Atlantic and Gulf Coasts of the United States. Necessarily, the US NRC research program includes both seismic- and landslide-based tsunamigenic sources in both the near and the far fields. The inclusion of tsunamigenic landslides, an important category of sources that impact tsunami hazard levels for the Atlantic and Gulf Coasts is a key difference between this program and most other tsunami hazard assessment programs. The initial phase of this work consisted of collection, interpretation, and analysis of available offshore data, with significant effort focused on characterizing offshore near-field landslides and analyzing their tsunamigenic potential and properties. In the next phase of research, additional field investigations will be conducted in key locations of interest and additional analysis will be undertaken. Simultaneously, the MOST tsunami generation and propagation model used by NOAA will first be enhanced to include landslide-based initiation mechanisms and then will be used to investigate the impact of the tsunamigenic sources identified and characterized by the USGS. The potential for probabilistic tsunami hazard assessment will also be explore in the final phases of the program.

Conference Paper↗

USGS2023

Constraints on near-ridge magmatism using 40Ar/39Ar geochronology of enriched MORB from the 8°20' N seamount chain

Our understanding of the spatial-temporal-compositional relationships between off-axis magmatism and mid-ocean ridge spreading centers is limited. Determining the 40 Ar/ 39 Ar ages of mid-ocean ridge basalt (MORB) lavas erupting near mid-ocean ridges (MOR) has been a challenge due to the characteristically low K 2 O contents in incompatible element-depleted normal MORB (NMORB). High-precision 40 Ar/ 39 Ar geochronology is used here to determine ages of young, basaltic lavas erupted along the 8°20' N seamount chain west of the East Pacific Rise (EPR) axis that have a range of incompatible element enrichments (EMORB) suitable for 40 Ar/ 39 Ar geochronology (e.g., K 2 O contents > 0.3 wt%). 40 Ar/ 39 Ar ages were determined in 29 well-characterized basalts sampled using HOV Alvin and dredging. Detailed geochronology and geochemical analyses provide important constraints on the timing, distribution, and origins of lavas that constructed this extensive volcanic lineament relative to magmatism beneath the adjacent EPR axis. Seamount eruption ages are up to ∼1.6 Ma younger than the underlying lithosphere, supporting a model of prolonged off-axis magmatism for at least 2 Myrs at distances as great as ∼90 km from the ridge axis. Increasing geochemical heterogeneity with eruption distance reflects the diminishing effect of sub-ridge melt focusing. The range of geochemically distinct lavas erupted at given distances from the ridge highlights the dynamic nature of the near-ridge magmatic environment over Myr timescales. Linear ridge-like (EPR-parallel) morphotectonic features erupt the youngest and most incompatible element-enriched lavas of the entire seamount chain, indicating there is a recent change in the influence of mantle heterogeneity and off-axis melt metasomatism on the near-ridge lithospheric mantle. Changes in seamount morphologies are attributed to counter-clockwise rotation and southward migration of the nearby Siqueiros transform over the last few million years.

Earth and Planetary Science Letters↗