SEARCH · Search USGS

Results for “Ocean Modeling”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

SourceResearch connectionFrom yearThrough year

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

At least 541 records · Page 30Linked to original sources

USGS2010

Two-dimensional time dependent hurricane overwash and erosion modeling at Santa Rosa Island

A 2DH numerical, model which is capable of computing nearshore circulation and morphodynamics, including dune erosion, breaching and overwash, is used to simulate overwash caused by Hurricane Ivan (2004) on a barrier island. The model is forced using parametric wave and surge time series based on field data and large-scale numerical model results. The model predicted beach face and dune erosion reasonably well as well as the development of washover fans. Furthermore, the model demonstrated considerable quantitative skill (upwards of 66% of variance explained, maximum bias - 0.21 m) in hindcasting the post-storm shape and elevation of the subaerial barrier island when a sheet flow sediment transport limiter was applied. The prediction skill ranged between 0.66 and 0.77 in a series of sensitivity tests in which several hydraulic forcing parameters were varied. The sensitivity studies showed that the variations in the incident wave height and wave period affected the entire simulated island morphology while variations in the surge level gradient between the ocean and back barrier bay affected the amount of deposition on the back barrier and in the back barrier bay. The model sensitivity to the sheet flow sediment transport limiter, which served as a proxy for unknown factors controlling the resistance to erosion, was significantly greater than the sensitivity to the hydraulic forcing parameters. If no limiter was applied the simulated morphological response of the barrier island was an order of magnitude greater than the measured morphological response.

Coastal Engineering↗

USGS2007

Two lithospheric profiles across southern California derived from gravity and seismic data

We present two detailed 2-D density transects for the crust and uppermost mantle across southern California using a linear gravity inversion technique. This technique parameterizes the crust and upper mantle as a set of blocks that are based on published geologic and seismic models. Each block can have a range of densities that are constrained where possible by borehole measurements, seismic velocities, and petrologic data. To further constrain the models, it is assumed that the lithosphere is close to isostatic equilibrium at both ends of the profiles, in the deep ocean and east of the Mojave Desert. We calculate the lithostatic pressure variations field for the whole cross section to rule out the geophysically insignificant solutions. In the linear equation, ρ = a + bV ( V , seismic P-wave velocity; ρ , density), which approximates the mantle density–velocity ( ρ – V ) relationship, different coefficients for b were evaluated. Lower coefficients ( b < 0.2) correspond to an almost purely thermally perturbed mantle, while higher coefficients ( b > 0.3) imply that other effects, such as composition and/or metamorphic changes, play an important role in the mantle. Density models were constructed with the coefficient b ranging from 0 to 0.6. The results indicate that a high b value in the mantle ρ – V relationship is associated with less dense crust in the Mojave block and more dense crust in the Catalina schist block. In the less dense Mojave block, the average density of the whole crust is ∼2.75 g/cm 3 , while that of the lower crust is ∼2.72 g/cm 3 . These densities imply a high silica content in the crust, and a minor fraction of basic rock in the lower crust, or perhaps the absence of a basaltic layer altogether. By comparison, the average density of a typical continental stable platform is ∼2.85 g/cm 3 . Models with higher b coefficients (0.5–0.6) are characterized by a large isostatic imbalance. On the other hand, lower b values (0–0.2) require a consolidated whole crust density in the Mojave Desert of ∼2.78 g/cm 3 , and a lower crust density of ∼2.89 g/cm 3 with mostly basaltic composition. This contradicts the observed, lower V p / V s -ratio in the Mojave Desert associated with mostly felsic and low-density crust. Models with lower b coefficients (0.1–0.2) are characterized by an absence of local Airy compensation beneath the San Gabriel Mountains at the LARSE-1 profile. These, and other non-gravity arguments, suggest optimal solutions to the mantle ρ – V relation of b ∼ 0.2–0.4. This, in turn, means that both thermal and petrological effects occur inside the downwelling of the uppermost mantle high velocity body located beneath the Transverse Ranges. During the development of this mantle downwelling, the basaltic layer of the Mojave block was likely eroded and pulled down into the high velocity body. Those basaltic fragments may have been transformed into eclogites, and this metamorphic change implies a higher b -coefficient density–velocity relationship than would be expected for a purely thermal process.

California↗

USGS2009

Far field tsunami simulations of the 1755 Lisbon earthquake: Implications for tsunami hazard to the U.S. East Coast and the Caribbean

The great Lisbon earthquake of November 1st, 1755 with an estimated moment magnitude of 8.5-9.0 was the most destructive earthquake in European history. The associated tsunami run-up was reported to have reached 5-15??m along the Portuguese and Moroccan coasts and the run-up was significant at the Azores and Madeira Island. Run-up reports from a trans-oceanic tsunami were documented in the Caribbean, Brazil and Newfoundland (Canada). No reports were documented along the U.S. East Coast. Many attempts have been made to characterize the 1755 Lisbon earthquake source using geophysical surveys and modeling the near-field earthquake intensity and tsunami effects. Studying far field effects, as presented in this paper, is advantageous in establishing constraints on source location and strike orientation because trans-oceanic tsunamis are less influenced by near source bathymetry and are unaffected by triggered submarine landslides at the source. Source location, fault orientation and bathymetry are the main elements governing transatlantic tsunami propagation to sites along the U.S. East Coast, much more than distance from the source and continental shelf width. Results of our far and near-field tsunami simulations based on relative amplitude comparison limit the earthquake source area to a region located south of the Gorringe Bank in the center of the Horseshoe Plain. This is in contrast with previously suggested sources such as Marqu??s de Pombal Fault, and Gulf of C??diz Fault, which are farther east of the Horseshoe Plain. The earthquake was likely to be a thrust event on a fault striking ~ 345?? and dipping to the ENE as opposed to the suggested earthquake source of the Gorringe Bank Fault, which trends NE-SW. Gorringe Bank, the Madeira-Tore Rise (MTR), and the Azores appear to have acted as topographic scatterers for tsunami energy, shielding most of the U.S. East Coast from the 1755 Lisbon tsunami. Additional simulations to assess tsunami hazard to the U.S. East Coast from possible future earthquakes along the Azores-Iberia plate boundary indicate that sources west of the MTR and in the Gulf of Cadiz may affect the southeastern coast of the U.S. The Azores-Iberia plate boundary west of the MTR is characterized by strike-slip faults, not thrusts, but the Gulf of Cadiz may have thrust faults. Southern Florida seems to be at risk from sources located east of MTR and South of the Gorringe Bank, but it is mostly shielded by the Bahamas. Higher resolution near-shore bathymetry along the U.S. East Coast and the Caribbean as well as a detailed study of potential tsunami sources in the central west part of the Horseshoe Plain are necessary to verify our simulation results. ?? 2008 Elsevier B.V.

Marine Geology↗

USGS2008

Complex faulting associated with the 22 December 2003 M w 6.5 San Simeon California, earthquake, aftershocks and postseismic surface deformation

We use data from two seismic networks and satellite interferometric synthetic aperture radar (InSAR) imagery to characterize the 22 December 2003 M w 6.5 San Simeon earthquake sequence. Absolute locations for the mainshock and nearly 10,000 aftershocks were determined using a new three-dimensional (3D) seismic velocity model; relative locations were obtained using double difference. The mainshock location found using the 3D velocity model is 35.704° N, 121.096° W at a depth of 9.7±0.7 km. The aftershocks concentrate at the northwest and southeast parts of the aftershock zone, between the mapped traces of the Oceanic and Nacimiento fault zones. The northwest end of the mainshock rupture, as defined by the aftershocks, projects from the mainshock hypocenter to the surface a few kilometers west of the mapped trace of the Oceanic fault, near the Santa Lucia Range front and the >5 mm postseismic InSAR imagery contour. The Oceanic fault in this area, as mapped by Hall (1991) , is therefore probably a second-order synthetic thrust or reverse fault that splays upward from the main seismogenic fault at depth. The southeast end of the rupture projects closer to the mapped Oceanic fault trace, suggesting much of the slip was along this fault, or at a minimum is accommodating much of the postseismic deformation. InSAR imagery shows ∼72 mm of postseismic uplift in the vicinity of maximum coseismic slip in the central section of the rupture, and ∼48 and ∼45 mm at the northwest and southeast end of the aftershock zone, respectively. From these observations, we model a ∼30-km-long northwest-trending northeast-dipping mainshock rupture surface—called the mainthrust—which is likely the Oceanic fault at depth, a ∼10-km-long southwest-dipping backthrust parallel to the mainthrust near the hypocenter, several smaller southwest-dipping structures in the southeast, and perhaps additional northeast-dipping or subvertical structures southeast of the mainshock plane. Discontinuous backthrust features opposite the mainthrust in the southeast part of the aftershock zone may offset the relic Nacimiento fault zone at depth. The InSAR data image surface deformation associated with both aseismic slip and aftershock production on the mainthrust and the backthrusts at the northwest and southeast ends of the aftershock zone. The well-defined mainthrust at the latitude of the epicenter and antithetic backthrust illuminated by the aftershock zone indicate uplift of the Santa Lucia Range as a popup block; aftershocks in the southeast part of the zone also indicate a popup block, but it is less well defined. The absence of backthrust features in the central part of the zone suggests range-front uplift by fault-propagation folding, or backthrusts in the central part were not activated during the mainshock.

California↗

USGS2015

Robust global ocean cooling trend for the pre-industrial Common Era

The oceans mediate the response of global climate to natural and anthropogenic forcings. Yet for the past 2,000 years — a key interval for understanding the present and future climate response to these forcings — global sea surface temperature changes and the underlying driving mechanisms are poorly constrained. Here we present a global synthesis of sea surface temperatures for the Common Era ( CE ) derived from 57 individual marine reconstructions that meet strict quality control criteria. We observe a cooling trend from 1 to 1800 CE that is robust against explicit tests for potential biases in the reconstructions. Between 801 and 1800 CE , the surface cooling trend is qualitatively consistent with an independent synthesis of terrestrial temperature reconstructions, and with a sea surface temperature composite derived from an ensemble of climate model simulations using best estimates of past external radiative forcings. Climate simulations using single and cumulative forcings suggest that the ocean surface cooling trend from 801 to 1800 CE is not primarily a response to orbital forcing but arises from a high frequency of explosive volcanism. Our results show that repeated clusters of volcanic eruptions can induce a net negative radiative forcing that results in a centennial and global scale cooling trend via a decline in mixed-layer oceanic heat content.

Nature Geoscience↗

USGS2023

Volcanic aquifers of Hawaiʻi—Contributions to assessing groundwater availability on Kauaʻi, Oʻahu, and Maui

The volcanic aquifers of the Hawaiian Islands supply water to 1.46 million residents, diverse industries, and a large component of the U.S. military in the Pacific. Groundwater also supplies fresh water that supports ecosystems in streams and near the coast. Hawaii’s aquifers are remarkably productive given their small size, but the capacity of the islands to store fresh groundwater is limited because each island is surrounded by seawater, and salt water underlies much of the fresh groundwater. The amount of fresh groundwater available for human use from Hawai‘i’s volcanic aquifers is constrained by the consequences of groundwater withdrawal. Restrictions placed on these consequences can translate to limitations on groundwater availability. Changes in recharge resulting from changes in land cover or climate can alter the effect of withdrawals. This study uses numerical models of the volcanic aquifers of the islands of Kaua‘i, O‘ahu, and Maui to quantify the consequences of historical and plausible future withdrawals and changes in recharge. The study compares the results of model simulations of multiple scenarios of historical and projected future withdrawal and recharge. Results of the simulations using the groundwater models of the islands of Kaua‘i, O‘ahu, and Maui have implications for other islands in Hawai‘i. Since the first modern water well was drilled in Hawai‘i in 1879, total groundwater withdrawals on Kaua‘i, O‘ahu, and Maui have risen to nearly 400 million gallons per day. Model simulations indicate that these withdrawals have caused reductions in groundwater discharge to streams and springs, reductions in groundwater discharge to the ocean, changes in subsurface flow between sectors within an island, lowering of groundwater levels, and rise of the interface between fresh water and salt water in the aquifers. Future increases in withdrawals will increase the severity of the consequences. Changes in recharge can alter the effect of withdrawals—increases in recharge can offset the consequences of withdrawals, whereas decreases in recharge can exacerbate the effects of withdrawals. This study quantifies the consequences of withdrawals for past and plausible future circumstances. The models can be used to test other circumstances. Limits placed on the consequences of withdrawals—such as restrictions to protect stream or coastal ecosystems that rely on groundwater discharge and limitations on water-level decline and rise of the freshwater-saltwater interface to protect the productivity of existing wells—can translate to limits on groundwater availability from Hawai‘i’s volcanic aquifers. Setting acceptable limits to the consequences of groundwater withdrawal is also a critical part of assessing groundwater availability. Once these limits are set, numerical models can be used to quantify the amount of water that can be withdrawn within those limits and thereby inform management decisions that seek to balance the need to limit the consequences of groundwater withdrawals with the need to develop water for human use.

Hawaii↗

USGS1979

Geohydrology and digital-simulation model of the Farrington aquifer in the northern coastal plain of New Jersey

A two-dimensional digital-computer flow model was developed to simulate the Farrington aquifer in the northern part of the Coastal Plain of New Jersey. The area of detailed study includes approximately 500 square miles in Middlesex and Monmouth Couties where the aquifer provides a large part of the municipal and industrial water supply. The area modeled is much larger, extending seaward as well as northeastward into Long Island. The aquifer consists chiefly of the Farrington Sand Member of the Raritan Formation and is composed of sand and some gravel. It thickens from a featheredge in outcrop to more than 170 feet, 11 miles to the southeast. The confining unit between the Farrington and the overlying Old Bridge Sand Member of the Magothy Formation consists primarily of the Woodbridge Clay Member of the Raritan Formation and has a maximum thickness of 244 feet. The model simulates both water-table and artesian conditions. The confining unit overlying the Farrington aquifer is simulated as having a variable thickness and vertical hydraulic conductivity. The effect of a declining water level in the overlying Old Bridge aquifer on the Farrington aquifer is also simulated by the model. Values used to describe the hydraulic properties of the Farrington aquifer are: a hydraulic conductivity of 105 feet per day, a storage coefficient of 1.6 x 10 -4 for artesian conditions, and a specific yield of 0.25 for water-table conditions. Values for the overlying confining unit are: a vertical hydraulic conductivity ranging from 4.2 x 10 -7 to 1.0 x 10 -10 feet per second and a specific storage of 4 x 10 -5 feet -1 . Aquifer simulation for the 15-year period, 1959-73, was used to calibrage the model. The model was calibrated by comparing the observed potentiometric surface of November 1973 with the simulated potentiometric surface. In addition, hydrographs for selected wells were compared with model results. Ground-water withdrawals for 1959 and 1973 were 12.1 and 28.5 milion gallons per day, respectively. Potentiometric surfaces for 1985 and 2000 were computed based on a linear projection of ground-water withdrawals (39.5 and 56.9 million gallons per day in 1985 and 2000, respectively) of the period 1959 through 1973. These surfaces are deeper than that of November, 1973, and the cone of depression is wider. The potentiometric head projected by the model in the vicinity of Sayreville will be more than 150 feet below mean sea level by 2000; the head in this area was 70 feet below sea level in 1973. The model calculated ground-water budgets for steady-state and transient conditions for the entire modeled area and for several rectangular subareas. Ground-water flow into the modeled Farrington aquifer under steady-state conditions before ground-water development was 16 cubic feet per second for the entire area. Recharge in the outcrop area and vertical leakage from the Old Bridge was 8 cubic feet per second each. Approximately 75 percent of the discharge occurred as seepage into surface-water bodies in and near the outcrop and as lateral flow southwestward into Burlington County near the outcrop area. The remaining 25 percent occurred southeast of the outcrop as vertical leakage into the overlying Old Bridge aquifer and as lateral flow to the south into Ocean and Burlington Counties. A transient water budget for 1973 was calculated for a subarea consisting mainly of Middlesex County. The model indicates that 48 percent (14.3 cubic feet per second) of the total inflow to the subareas was through its boundaries. Other sources of water include direct recharge within the subarea (5.4 cubic feet per second), vertical leakage (mainly from the Old Bridge) within the subarea (2.6 cubic feet per second), and water released from storage (3.4 cubic feet per second). Discharge from the subarea consisted mainly of withdrawals (26.5 cubic feet per second). It also included vertical leakage to the Old Bridge and discharge to surface-water bodies simulated by constant-head nodes (3.2 cubic feet per second).

New Jersey↗

USGS1999

Trace-element geochemistry of metabasaltic rocks from the Yukon-Tanana Upland and implications for the origin of tectonic assemblages in east-central Alaska

We present major- and trace- element geochemical data for 27 amphibolites and six greenstones from three structural packages in the Yukon-Tanana Upland of east-central Alaska: the Lake George assemblage (LG) of Devono-Mississippian augen gneiss, quartz-mica schist, quartzite, and amphibolite; the Taylor Mountain assemblage (TM) of mafic schist and gneiss, marble, quartzite, and metachert; and the Seventymile terrane of greenstone, serpentinized peridotite, and Mississippian to Late Triassic metasedimentary rocks. Most LG amphibolites have relatively high Nb, TiO2, Zr, and light rare earth element contents, indicative of an alkalic to tholeiitic, within-plate basalt origin. The within-plate affinities of the LG amphibolites suggest that their basaltic parent magmas developed in an extensional setting and support a correlation of these metamorphosed continental-margin rocks with less metamorphosed counterparts across the Tintina fault in the Selwyn Basin of the Canadian Cordillera. TM amphibolites have a tholeiitic or calc-alkalic composition, low normalized abundances of Nb and Ta relative to Th and La, and Ti/V values of <20, all indicative of a volcanic-arc origin. Limited results from Seventymile greenstones indicate a tholeiitic or calc-alkalic composition and intermediate to high Ti/V values (27-48), consistent with either a within-plate or an ocean-floor basalt origin. Y-La-Nb proportions in both TM and Seventymile metabasalts indicate the proximity of the arc and marginal basin to continental crust. The arc geochemistry of TM amphibolites is consistent with a model in which the TM assemblage includes arc rocks generated above a west-dipping subduction zone outboard of the North American continental margin in mid-Paleozoic through Triassic time. The ocean-floor or within-plate basalt geochemistry of the Seventymile greenstones supports the correlation of the Seventymile terrane with the Slide Mountain terrane in Canada and the hypothesis that these oceanic rocks originated in a basin between the continental margin and an arc to the west.

Alaska↗

USGS2017

Skagit River coho salmon life history model—Users’ guide

Natural resource management is conducted in the context of multiple anthropogenic stressors and is further challenged owing to changing climate. Experiments to determine the effects of climate change on complex ecological systems are nearly impossible. However, using a simulation model to synthesize current understanding of key ecological processes through the life cycle of a fish population can provide a platform for exploring potential effects of and management responses to changing conditions. Potential climate-change scenarios can be imposed, responses can be observed, and the effectiveness of potential actions can be evaluated. This approach is limited owing to future conditions likely deviating in range and timing from conditions used to create the model so that the model is expected to become obsolete. In the meantime, however, the modeling process explicitly states assumptions, clarifies information gaps, and provides a means to better understand which relationships are robust and which are vulnerable to changing climate by observing whether and why model output diverges from actual observations through time. The purpose of the model described herein is to provide such a decision-support tool regarding coho ( Oncorhynchus kisutch ) salmon for the Sauk-Suiattle Indian Tribe of Washington State. The Skagit coho salmon model is implemented in a system dynamics format and has three primary stocks—(1) predicted smolts, (2) realized smolts, and (3) escapement. “Predicted smolts” are the number of smolts expected based on the number of spawners in any year and the Ricker production curve. Pink salmon ( Oncorhynchus gorbuscha ) return to the Skagit River in odd years, and when they overlap with juvenile rearing coho salmon, coho smolt production is substantially higher than in non-pink years. Therefore, the model uses alternative Ricker equations to predict smolts depending on whether their juvenile year was a pink or non-pink year. The stock “realized smolts” is calculated based on the expected effect of streamflow conditions to alter the productivity predicted by the Ricker curve. Adverse conditions include scouring flow events that occur when redds are present; high-flow events during winter on juveniles, which can cause fish displacement and adverse water turbidity; and extremely low flows in summer. The stock “escapement” represents the fish remaining after accounting for ocean mortality and harvest. Ocean mortality has been linked with indices of ocean conditions, which are related to ocean biological productivity. Ocean survival also may have a density-dependent component such that lower survival is associated with higher numbers of smolts. The model allows the user to change certain model parameters and inputs, and choose among alternative predictors for certain modeled relations.

Washington↗

USGS2022

Quantifying the effects of tides, river flow, and barriers on movements of Chinook Salmon smolts at junctions in the Sacramento–San Joaquin River Delta using multistate models

Successful migration of Chinook Salmon ( Oncorhynchus tshawytscha ) smolts seaward in the Sacramento – San Joaquin River Delta (hereafter, Delta) requires navigating a network of numerous branching channels. Within the Delta, several key junctions route smolts either towards more direct paths to the ocean or towards the interior Delta, an area associated with decreased survival. Movements within these junctions that determine route choice can be influenced by numerous behavioral and environmental factors, including the complex interplay between tidal and riverine hydraulics. Here, we apply continuous time multistate Markov models to examine the influence of tidal and riverine hydraulics, behavioral factors, and management actions on smolt movements. These models incorporate more information from acoustic telemetry data compared with previous approaches to modeling smolt movements in the Delta. By decomposing modeled flows into tidal and net flow signals we elucidate how each component influences movements into and out of distributary channels. Increasing net flows generally increased movement rates, while flood tides decreased seaward movement rates. Similarly, ebb tides increased downstream movements as fish go with the flow. We found less support for diel movement behaviors compared to flow metrics. Additionally, we quantify the effects of a large management action, the placement of a physical barrier, which was effective at decreasing entrainment into the interior Delta. Together, these results help inform management of Chinook Salmon and increase our understanding of the major factors driving smolt movements within these key junctions.

California↗

USGS2015

Slab melting beneath the Cascades Arc driven by dehydration of altered oceanic peridotite

Water is returned to Earth’s interior at subduction zones. However, the processes and pathways by which water leaves the subducting plate and causes melting beneath volcanic arcs are complex; the source of the water—subducting sediment, altered oceanic crust, or hydrated mantle in the downgoing plate—is debated; and the role of slab temperature is unclear. Here we analyse the hydrogen-isotope and trace-element signature of melt inclusions in ash samples from the Cascade Arc, where young, hot lithosphere subducts. Comparing these data with published analyses, we find that fluids in the Cascade magmas are sourced from deeper parts of the subducting slab—hydrated mantle peridotite in the slab interior—compared with fluids in magmas from the Marianas Arc, where older, colder lithosphere subducts. We use geodynamic modelling to show that, in the hotter subduction zone, the upper crust of the subducting slab rapidly dehydrates at shallow depths. With continued subduction, fluids released from the deeper plate interior migrate into the dehydrated parts, causing those to melt. These melts in turn migrate into the overlying mantle wedge, where they trigger further melting. Our results provide a physical model to explain melting of the subducted plate and mass transfer from the slab to the mantle beneath arcs where relatively young oceanic lithosphere is subducted.

Cascade Arc↗

USGS2025

Declining marine survival of steelhead trout linked to climate and ecosystem change

Species with complex life cycles, such as anadromous fish that perform spawning migrations between freshwater and the ocean, may be particularly sensitive to global change because freshwater and marine habitats experience distinct shifts in climate and ecosystem dynamics. Abundances of wild steelhead trout ( Oncorhynchus mykiss ) have declined across most of their range over the past 40–50 years. We examined whether declines in steelhead survival can be linked to changing climate conditions and species interactions. A novel hierarchical integrated population model that accounts for the species' complex life history was fitted to data from multiple wild steelhead populations on the Washington coast, U.S.A. The model estimates recruitment residuals and kelt survival rates as time-varying processes, which reflect annual variation in survival before and after first maturation. We found that survival rates of immature steelhead (recruits) and adult steelhead (kelts) have declined over time and that survival trends across populations were strongly associated with climate and ecosystem change, specifically summer sea surface temperature and pink salmon abundance in the North Pacific Ocean, the NPGO index and river flows. Including these drivers in the model reduced unexplained annual variation in shared recruitment and kelt survival anomalies and largely accounted for their negative long-term trends. Our findings provide evidence that rising temperatures and increased interspecific competition at sea have contributed to declines in steelhead survival over the last five decades. Considering projected warming and high pink salmon abundances in the ocean, steelhead will likely continue to experience low marine survival rates.

Washington↗

USGS1996

Crustal and upper mantle velocity structure of the Salton Trough, southeast California

This paper presents data and modelling results from a crustal and upper mantle wide-angle seismic transect across the Salton Trough region in southeast California. The Salton Trough is a unique part of the Basin and Range province where mid-ocean ridge/transform spreading in the Gulf of California has evolved northward into the continent. In 1992, the U.S. Geological Survey (USGS) conducted the final leg of the Pacific to Arizona Crustal Experiment (PACE). Two perpendicular models of the crust and upper mantle were fit to wide-angle reflection and refraction travel times, seismic amplitudes, and Bouguer gravity anomalies. The first profile crossed the Salton Trough from the southwest to the northeast, and the second was a strike line that paralleled the Salton Sea along its western edge. We found thin crust (∼21–22 km thick) beneath the axis of the Salton Trough (Imperial Valley) and locally thicker crust (∼27 km) beneath the Chocolate Mountains to the northeast. We modelled a slight thinning of the crust further to the northeast beneath the Colorado River (∼24 km) and subsequent thickening beneath the metamorphic core complex belt northeast of the Colorado River. There is a deep, apparently young basin (∼5–6 km unmetamorphosed sediments) beneath the Imperial Valley and a shallower (∼2–3 km) basin beneath the Colorado River. A regional 6.9-km/s layer (between ∼15-km depth and the Moho) underlies the Salton Trough as well as the Chocolate Mountains where it pinches out at the Moho. This lower crustal layer is spatially associated with a low-velocity (7.6–7.7 km/s) upper mantle. We found that our crustal model is locally compatible with the previously suggested notion that the crust of the Salton Trough has formed almost entirely from magmatism in the lower crust and sedimentation in the upper crust. However, we observe an apparently magmatically emplaced lower crust to the northeast, outside of the Salton Trough, and propose that this layer in part predates Salton Trough rifting. It may also in part result from migration of magmatic spreading centers associated with the southern San Andreas fault system. These spreading centers may have existed east of their current locations in the past and may have influenced the lower crust and upper mantle to the east of the current Salton Trough.

California↗

USGS2014

Interactions between waves, sediment, and turbulence on a shallow estuarine mudflat

Measurements were collected on a shallow estuarine mudflat in northern San Francisco Bay to examine the physical processes controlling waves, turbulence, sediment resuspension, and their interactions. Tides alone forced weak to moderate currents of 10–30 cm s -1 in depths of 0–3 m, and maintained a background suspension of 30–50 mg L21 of fine sediment. In the presence of wind waves, bottom orbital velocities spanned 20–30 cm s -1 , suspended-sediment concentrations (SSC) at 15 and 30 cm above the bed (cmab) increased by 1–2 orders of magnitude, and vertical gradients in SSC were strong enough to produce turbulence-limiting stratification, with gradient Richardson numbers exceeding 0.25. Simultaneously, turbulent stresses (decomposed from wave motions) increased by an order of magnitude. The apparent contradiction of energetic turbulence in the presence of strong stratification was reconciled by considering the turbulent kinetic energy (TKE) budget: in general, dissipation and buoyancy flux were balanced by local shear production, and each of these terms increased during wave events. The classic wave-current boundary layer model represented the observations qualitatively, but not quantitatively since the velocity profile could not be approximated as logarithmic. Rather, the mean shear was elevated by the Stokes drift return flow and wind-generated surface stress, which diffused sediment upward and limited stratification. Our findings highlight a pathway for waves to supply energy to both the production and destruction of turbulence, and demonstrate that in such shallow depths, TKE and SSC can be elevated over more of the water column than predicted by traditional models.

California↗

USGS2008

Seasonal variation in sensitivity of larval sea lampreys to the lampricide 3-trifluoromethyl-4-nitrophenol

We evaluated the sensitivity of larval sea lampreys Petromyzon marinus to the lampricide 3-trifluoromethyl-4-nitrophenol (TFM) in a series of toxicity tests in spring and summer. Although noted previously, the seasonal variation in sensitivity to TFM had never been tested as a means of reducing TFM usage in stream treatments. A preliminary study consisted of three spring and four summer static toxicity tests conducted at 12??C. A more comprehensive study consisted of 12 spring and summer paired flow-through toxicity tests conducted both at seasonal water temperatures and at 12??C. The sensitivity of larval sea lampreys to TFM was greater in spring than in summer. The preliminary static toxicity tests indicated that the concentration of TFM needed to kill larval sea lampreys in spring (May and June) was about one-half that required in summer (August); the concentrations lethal to 50% and 99.9% of the test animals (the LC50 and LC99.9 values) were less in spring than in summer. Analysis of variance of the flow-through toxicity data indicated that season significantly affected both the LC50 and LC99.9 values. For all 12 paired flow-through toxicity tests, the spring LC50 and LC99.9 values were less than the corresponding summer values. For 9 of the 12 paired flow-through toxicity tests, the dose-response toxicity lines were parallel and allowed statistical comparison of the LC50 values. The spring LC50 values were significantly lower than the summer values in eight of the nine tests. Verification of a seasonal variation in the sensitivity of larval sea lampreys to TFM will allow inclusion of this factor in the selection model currently used by both the U.S. Fish and Wildlife Service and the Department of Fisheries and Oceans-Canada to schedule lampricide stream treatments. ?? Copyright by the American Fisheries Society 2008.

North American Journal of Fisheries Management↗

USGS1983

An interpretation of gravity and aeromagnetic surveys of the Greater Jiddah area, Kingdom of Saudi Arabia

A gravity survey of the Jiddah area between lat 21°23' and 21°48' N., long 39°03' and 39°20' E., carried out in early 1977, covers an area of 744 km 2 and includes 215 stations, at an average spacing of 1.9 km. The resulting simple Bouguer gravity anomaly map shows a series of four large (7-15 mgal) gravity anomaly highs along a north-trending gravity gradient located immediately west of outcrops of Precambrian rocks. The regional Bouguer gravity anomaly field of the map area decreases eastward at an approximately constant rate of about 1.1 mgal-km-1. Several anomalies suggest northeast-trending, right-lateral strike-slip faulting throughout the area. Interpretation of the Bouguer gravity anomaly map, combined with geologic, aeromagnetic, and paleomagnetic data, leads to an interpretive model in which the gravity gradient marks the boundary zone between Precambrian continental crust to the east and oceanic sea-floor crust overlain by marine sedimentary rocks to the west. The strike-slip faults in the area are probably the landward extension of offshore Red Sea transform faults, and, because they penetrate and offset continental crust, their study may yield useful information about the mechanism of transform faulting. The gravity anomaly highs are interpreted to result from shallow intrusive bodies, probably gabbroic in composition, that are buried a few hundred meters below the surface and that are approximately 2 km thick.

Report↗

USGS2012

Free tropospheric transport of microorganisms from Asia to North America

Microorganisms are abundant in the troposphere and can be transported vast distances on prevailing winds. This study measures the abundance and diversity of airborne bacteria and fungi sampled at the Mt. Bachelor Observatory (located 2.7 km above sea level in North America) where incoming free tropospheric air routinely arrives from distant sources across the Pacific Ocean, including Asia. Overall deoxyribonucleic acid (DNA) concentrations for microorganisms in the free troposphere, derived from quantitative polymerase chain reaction assays, averaged 4.94 × 10(-5) ng DNA m(-3) for bacteria and 4.77 × 10(-3) ng DNA m(-3) for fungi. Aerosols occasionally corresponded with microbial abundance, most often in the springtime. Viable cells were recovered from 27.4 % of bacterial and 47.6 % of fungal samples (N = 124), with 49 different species identified by ribosomal DNA gene sequencing. The number of microbial isolates rose significantly above baseline values on 22-23 April 2011 and 13-15 May 2011. Both events were analyzed in detail, revealing distinct free tropospheric chemistries (e.g., low water vapor, high aerosols, carbon monoxide, and ozone) useful for ruling out boundary layer contamination. Kinematic back trajectory modeling suggested air from these events probably originated near China or Japan. Even after traveling for 10 days across the Pacific Ocean in the free troposphere, diverse and viable microbial populations, including presumptive plant pathogens Alternaria infectoria and Chaetomium globosum, were detected in Asian air samples. Establishing a connection between the intercontinental transport of microorganisms and specific diseases in North America will require follow-up investigations on both sides of the Pacific Ocean.

Microbial Ecology↗

USGS2010

Modeling downstream fining in sand-bed rivers. II: Application

In this paper the model presented in the companion paper, Wright and Parker (2005) is applied to a generic river reach typical of a large, sand-bed river flowing into the ocean in order to investigate the mechanisms controlling longitudinal profile development and downstream fining. Three mechanisms which drive downstream fining are studied: a delta prograding into standing water, sea-level rise, and tectonic subsidence. Various rates of sea-level rise (typical of the late Holocene) and tectonic subsidence are modeled in order to quantify their effects on the degree of profile concavity and downstream fining. Also, several other physical mechanisms which may affect fining are studied, including the relative importance of the suspended versus bed load, the effect of the loss of sediment overbank, and the influence of the delta bottom slope. Finally, sensitivity analysis is used to show that the grain-size distribution at the interface between the active layer and substrate has a significant effect on downstream fining.

Journal of Hydraulic Research↗