Search USGSSearch

SEARCH · Search USGS

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 361 records · Page 20Linked to original sources

Recent volcanism in the Siqueiros transform fault: Picritic basalts and implications for MORB magma genesis

Small constructional volcanic landforms and very fresh-looking lava flows are present along one of the inferred active strike-slip faults that connect two small spreading centers (A and B) in the western portion of the Siqueiros transform domain. The most primitive lavas (picritic and olivine-phyric basalts), exclusively recovered from the young-looking flows within the A-B strike-slip fault, contain millimeter-sized olivine phenocrysts (up to 20 modal%) that have a limited compositional range (Fo91.5-Fo89.5) and complexly zoned Cr-Al spinels. High-MgO (9.5-10.6 wt%) glasses sampled from the young lava flows contain 1-7% olivine phenocrysts (Fo90.5-Fo89) that could have formed by equilibrium crystallization from basaltic melts with Mg# values between 71 and 74. These high MgO (and high Al2O3) glasses may be near-primary melts from incompatible-element depleted oceanic mantle and little modified by crustal mixing and/or fractionation processes. Phase chemistry and major element systematics indicate that the picritic basalts are not primary liquids and formed by the accumulation of olivine and minor spinel from high-MgO melts (10% < MgO < 14%). Compared to typical N-MORB from the East Pacific Rise, the Siqueiros lavas are more primitive and depleted in incompatible elements. Phase equilibria calculations and comparisons with experimental data and trace element modeling support this hypothesis. They indicate such primary mid-ocean ridge basalt magmas formed by 10-18% accumulative decompression melting in the spinel peridotite field (but small amounts of melting in the garnet peridotite field are not precluded). The compositional variations of the primitive magmas may result from the accumulation of different small batch melt fractions from a polybaric melting column.

Earth and Planetary Science Letters

The use of sea ice habitat by female polar bears in the Beaufort Sea

Polar bears ( Ursus maritimus ) depend on ice-covered seas to satisfy life history requirements. Modern threats to polar bears include oil spills in the marine environment and changes in ice composition resulting from climate change. Managers need practical models that explain the distribution of bears in order to assess the impacts of these threats. We used stepwise procedures to create resource selection models of habitat use for radio-collared female polar bears in the Beaufort Sea. Sea ice characteristics and ocean depths at known polar bear locations were compared to the same features at randomly selected locations. Models generated for each of four seasons confirmed complexities of habitat use by polar bears and their response to numerous factors. Bears preferred shallow water areas where ice concentrations were > 80 % and different ice types intersected. Variation among seasons was reflected mainly in differential selection of ice stages, floe sizes, and their interactions. Water depth, total ice concentration and distance to the nearest interface between different ice types were significant terms in models for most seasons. Variation in ice stage and form also appeared in three models, and several interaction effects were identified. Habitat selection by polar bears is likely related to prey abundance and availability. Use of habitats in shallow water possibly reflects higher productivity in those areas. Habitat use in close proximity to ice edges is probably related to greater access of prey in those habitats.

Conference Paper

Delineation of marsh types and marsh-type change in coastal Louisiana for 2007 and 2013

The Bureau of Ocean Energy Management researchers often require detailed information regarding emergent marsh vegetation types (such as fresh, intermediate, brackish, and saline) for modeling habitat capacities and mitigation. In response, the U.S. Geological Survey in cooperation with the Bureau of Ocean Energy Management produced a detailed change classification of emergent marsh vegetation types in coastal Louisiana from 2007 and 2013. This study incorporates two existing vegetation surveys and independent variables such as Landsat Thematic Mapper multispectral satellite imagery, high-resolution airborne imagery from 2007 and 2013, bare-earth digital elevation models based on airborne light detection and ranging, alternative contemporary land-cover classifications, and other spatially explicit variables. An image classification based on image objects was created from 2007 and 2013 National Agriculture Imagery Program color-infrared aerial photography. The final products consisted of two 10-meter raster datasets. Each image object from the 2007 and 2013 spatial datasets was assigned a vegetation classification by using a simple majority filter. In addition to those spatial datasets, we also conducted a change analysis between the datasets to produce a 10-meter change raster product. This analysis identified how much change has taken place and where change has occurred. The spatial data products show dynamic areas where marsh loss is occurring or where marsh type is changing. This information can be used to assist and advance conservation efforts for priority natural resources.

Louisiana

On the global variations of terrestrial heat-flow

Over 3 500 measurements of surface heat-flux have been catalogued and analyzed to study the large-scale variations of terrestrial heat-flow. It was found that heat-flow values are correlated with major geologic provinces: higher averages and scattered values in active tectonic regions, and lower averages and more uniform values in stable areas. Analyzing the data in the light of new global tectonics shows that the variations of heat-flow are consistent with the hypotheses of sea-floor spreading and plate tectonics. The observed heat-flow across the mid-oceanic ridges can be accounted for by a simple model of a spreading sea floor.

Physics of the Earth and Planetary Interiors

Calcareous plankton biostratigraphic fidelity and species richness during the last 10 m.y. of the Cretaceous at Blake Plateau, subtropical North Atlantic

Species distributions of well-preserved and diverse assemblages of planktonic foraminifera and calcareous nannofossils spanning the last 10 m.y. of the Cretaceous (middle Campanian through Maastrichtian) are analyzed from samples taken across a 1400 m depth transect at Blake Nose in the western subtropical North Atlantic (Ocean Drilling Program Sites 1049, 1050 and 1052). Age models constructed by integrating foraminiferal, calcareous nannofossil, and magnetic polarity datum events provide a reliable framework for temporal correlation of the sites. This framework enables comparisons of species richness and abundance among sites and evaluation of the reliability of first and last appearance datums for regional and global correlation. Among the standard primary zonal marker datums, six of nine planktonic foraminifer and six of seven calcareous nannofossil events are considered reliable for constraining the age-depth models. Secondary datum ages calculated for 17 planktonic foraminiferal events suggest correlation offsets among the three sites of <0.1 m.y. for four species (including Pseudoguembelina praehariaensis Tur and Huber n. sp.), 0.1–0.5 m.y. for six events, and 0.6–1.0 m.y. for eight species. Secondary calcareous nannofossil datum ages calculated for six species show less reliability, with offsets of <0.1 m.y. for one species, 0.7 m.y. for one species, and 1.0–1.5 m.y. for the remaining four species. The distinctly identifiable new species Trinitella suturis Tur and Huber has lowest occurrences that are diachronous by as much as 1.73 m.y. among the Blake Nose sites. Occurrence rarity is the most likely explanation for the age offsets of this and other diachronous species. Planktonic foraminiferal assemblages show no significant differences in species composition and relative abundance among the three sites, suggesting the sites were all located below a single oligotrophic surface water mass. Species richness counts pooled in 200 kyr bins for the Blake Nose sites reveals high species origination rates from the late Campanian through early Maastrichtian and highest species richness (51–63 species) during the Maastrichtian. The only significant extinction pulse during the mid-Campanian through Maastrichtian occurs between 66.2 and 66.4 Ma with loss of five species representing ∼9% of the assemblage. These extinctions occur at the same time as a globally recognized warming event correlated with a pulse of eruptions at the Deccan Traps in India. Calcareous nannofossil assemblages show no significant change in relative abundance among the three sites. Two significant extinction events are documented: one from 75.20 to 75.40 Ma with a loss of six species representing ∼6% of the assemblage and one from 66.2 to 66.4 Ma with a loss of nine species representing ∼9% of the assemblage. The former event is associated with a hiatus at the base of Zone CC24, and the latter corresponds to Deccan Trap warming. Hiatuses are identified at all three Blake Nose sites near the base of the Maastrichtian (∼71.5 Ma; lowermost Pseudoguembelina palpebra /CC24 Zone) and only at the deeper Sites 1049 and 1050 in the mid-Maastrichtian (∼67.2 Ma; Pseudoguembelina hariaensis /CC26b Zone) and the mid-Campanian (∼75.9 Ma; base of Radotruncana calcarata /CC22 Zone). Slumping from across the shelf and slope could have caused the early Maastrichtian hiatus while changes in the pattern and strength of deep-water circulation may have been responsible for the mid-Maastrichtian and mid-Campanian hiatuses.

Cretaceous Research

Heat-producing elements in the lunar mantle: Insights from ion microprobe analyses of lunar pyroclastic glasses

We provide new estimates for the abundance of heat-producing elements in the lunar mantle by using SIMS techniques to measure the concentrations of thorium and samarium in lunar pyroclastic glasses. Lunar pyroclastic glasses are utilized in this study because they represent quenched products of near-primary melts from the lunar mantle and as such, they provide compositional information about the mantle itself. Thorium and samarium were measured because: (1) Th is not significantly fractionated from Sm during partial melting of the pyroclastic glass source regions, which are dominated by olivine and pyroxene. Therefore, the Th/Sm ratios that we measure in the pyroclastic glasses reflect the Th/Sm ratio of the pyroclastic glass source regions. (2) Strong correlations between Th, U, and K on the Moon allow us to use measured Th concentrations to estimate the concentrations of U and K in the pyroclastic glasses. (3) Th, Sm, U, and K are radioactive elements and as such, their concentrations can be used to investigate heat production in the lunar mantle. The results from this study show that the lunar mantle is heterogeneous with respect to heat-producing elements and that there is evidence for mixing of a KREEP component into the source regions of some of the pyroclastic glasses. Because the source regions for many of the glasses are deep (⩾400 km), we propose that a KREEP component was transported to the deep lunar mantle. KREEP enriched sources produce 138% more heat than sources that do not contain KREEP and therefore, could have provided a source of heat for extended periods of nearside basaltic magmatism. Data from this study, in conjunction with models for the fractional crystallization of a lunar magma ocean, are used to show that the average lunar mantle contains 0.15 ppm Th, 0.54 ppm Sm, 0.039 ppm U, and 212 ppm K. This is a greater enrichment in radiogenic elements than some earlier estimates, suggesting a more prolonged impact of radiogenic heat on nearside basaltic volcanism.

Geochimica et Cosmochimica Acta

Nonlinear refraction and reflection travel time tomography

We develop a rapid nonlinear travel time tomography method that simultaneously inverts refraction and reflection travel times on a regular velocity grid. For travel time and ray path calculations, we apply a wave front method employing graph theory. The first-arrival refraction travel times are calculated on the basis of cell velocities, and the later refraction and reflection travel times are computed using both cell velocities and given interfaces. We solve a regularized nonlinear inverse problem. A Laplacian operator is applied to regularize the model parameters (cell slownesses and reflector geometry) so that the inverse problem is valid for a continuum. The travel times are also regularized such that we invert travel time curves rather than travel time points. A conjugate gradient method is applied to minimize the nonlinear objective function. After obtaining a solution, we perform nonlinear Monte Carlo inversions for uncertainty analysis and compute the posterior model covariance. In numerical experiments, we demonstrate that combining the first arrival refraction travel times with later reflection travel times can better reconstruct the velocity field as well as the reflector geometry. This combination is particularly important for modeling crustal structures where large velocity variations occur in the upper crust. We apply this approach to model the crustal structure of the California Borderland using ocean bottom seismometer and land data collected during the Los Angeles Region Seismic Experiment along two marine survey lines. Details of our image include a high-velocity zone under the Catalina Ridge, but a smooth gradient zone between. Catalina Ridge and San Clemente Ridge. The Moho depth is about 22 km with lateral variations. Copyright 1998 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth

The role of Northeast Pacific meltwater events in deglacial climate change

Columbia River megafloods occurred repeatedly during the last deglaciation, but the impacts of this fresh water on Pacific hydrography are largely unknown. To reconstruct changes in ocean circulation during this period, we used a numerical model to simulate the flow trajectory of Columbia River megafloods and compiled records of sea surface temperature, paleo-salinity, and deep-water radiocarbon from marine sediment cores in the Northeast Pacific. The North Pacific sea surface cooled and freshened during the early deglacial (19.0-16.5 ka) and Younger Dryas (12.9-11.7 ka) intervals, coincident with the appearance of subsurface water masses depleted in radiocarbon relative to the sea surface. We infer that Pacific meltwater fluxes contributed to net Northern Hemisphere cooling prior to North Atlantic Heinrich Events, and again during the Younger Dryas stadial. Abrupt warming in the Northeast Pacific similarly contributed to hemispheric warming during the Bølling and Holocene transitions. These findings underscore the importance of changes in North Pacific freshwater fluxes and circulation in deglacial climate events.

Science Advances

Two tickets to paradise: multiple dispersal events in the founding of hoary bat populations in Hawai'i

The Hawaiian islands are an extremely isolated oceanic archipelago, and their fauna has long served as models of dispersal in island biogeography. While molecular data have recently been applied to investigate the timing and origin of dispersal events for several animal groups including birds, insects, and snails, these questions have been largely unaddressed in Hawai'i's only native terrestrial mammal, the Hawaiian hoary bat, Lasiurus cinereus semotus . Here, we use molecular data to test the hypotheses that (1) Hawaiian L . c . semotus originated via dispersal from North American populations of L . c . cinereus rather than from South American L . c . villosissimus , and (2) modern Hawaiian populations were founded from a single dispersal event. Contrary to the latter hypothesis, our mitochondrial data support a biogeographic history of multiple, relatively recent dispersals of hoary bats from North America to the Hawaiian islands. Coalescent demographic analyses of multilocus data suggest that modern populations of Hawaiian hoary bats were founded no more than 10 kya. Our finding of multiple evolutionarily significant units in Hawai'i highlights information that should be useful for re-evaluation of the conservation status of hoary bats in Hawai'i.

PLoS ONE

Paleoceanographic perspectives on Arctic Ocean change

The Arctic Ocean is presently experiencing changes in ocean temperature and sea ice extent that are unprecedented in the observational time period (satellite observations: 1979-Present). To provide context for the current changes, scientists turn to paleo records of past climate to document and study natural variability in the Arctic system. Paleoceanographic records that extend limited Arctic instrumental measurements are central to improving our understanding of sea ice dynamics and ocean warming and for enhancing the predictive capability of models. By coupling paleoceanographic records with modern observations, scientists can also contextualize the rate and magnitude of modern change with the deep past.

Arctic Report Card

The Holocene sea-level highstand in the equatorial Pacific: Analysis of the insular paleosea-level database

A review of the literature provides 92 estimates of the middle to late Holocene sea-level highstand on Pacific Islands. These data generally support geophysical model calculations that predict a +1 to 3 m relative sea-level highstand on oceanic islands due to the Earth’s rheological response to the melting of the last continental ice sheets and subsequent redistribution of meltwater. Both predictions and observations indicate sea level was higher than present in the equatorial Pacific between 5000 and 1500 y B.P. A non-linear relationship exists between the age and elevation of the highstand peak, suggesting that different rates of isostatic adjustment may occur in the Pacific, with the highest rates of sea-level fall following the highstand near the equator. It is important to resolve detailed sea-level histories from insular sites to test and refine models of climatic, oceanographic, and geophysical processes including hydroisostasy, equatorial ocean siphoning, and lithospheric flexure that are invoked as mechanisms affecting relative sea-level position. We use a select subset of the available database meeting specific criteria to examine model relationships of paleosea-surface topography. This new evaluated database of paleosea-level positions is also validated for testing and constraining geophysical model predictions of past and present sea-level variations.

Coral Reefs

Dissolved Fe supply to the central Gulf of Alaska is inferred to be derived from Alaskan glacial dust that is not resolved by dust transport models

Re-examination of previously published dissolved iron time-series data from Ocean Station Papa in the central Gulf of Alaska (GoA) reveals 33-70% increases in the dissolved iron inventories occurring between September and February of successive years, implying a source of Fe to this region during autumn or early winter. Because I can virtually rule out many possible iron sources at this time of year, I suggest Alaskan glacial dust is the likely iron source. Large plumes of such dust are known to be generated regularly in the autumn by anomalous offshore winds and channelled through mountain gaps, simultaneously from several locations spanning ∼1000 km of the northern Gulf of Alaska coastline. Large dust flux events occur when below-freezing, low-humidity air temperatures persist for many days during the autumn. I suggest that existing state-of-the-art global dust models fail to reproduce this Alaskan dust flux because the model spatial resolution is too coarse to resolve the high winds through the narrow mountain gaps that generate the dust. Future work that could help to confirm this Fe source to the central GoA includes time-series profiles of iron concentrations, and ancillary information from sensor-equipped profiling floats. If this mechanism of Fe supply to the central GoA were confirmed, it would imply this Alaskan dust is transported ≥ 1100 km from the coast, more than twice as far as has been visually documented from satellite observations.

Alaska

Storm surge modeling and applications in coastal areas

This chapter introduces the reader to a wide spectrum of storm surge modeling systems used to assess the impact of tropical cyclones, covering a range of numerical methods, model domains, forcing and boundary conditions, and purposes. New technologies to obtain data such as deployment of temporary sensors and remote sensing practices to support modeling are also presented. Extensive storm surge modeling applications have been made with existing modeling systems and some of them are described in this chapter. The authors recognize the importance of evaluating river-ocean interactions in coastal environments during tropical cyclones. Therefore, the coupling of hydraulic (riverine) and storm surge models is discussed. In addition, results from studies performed in the coast of India are shown which generated maps to help emergency managers and reduce risk due to coastal inundation.

Book chapter

A Systems Thinking approach to post-disaster restoration of maritime transportation systems

A Systems Thinking approach is used to examine elements of a maritime transportation system that are most likely to be impacted by an extreme event. The majority of the literature uses a high-level view that can fail to capture the damage at the sub-system elements. This work uses a system dynamics simulation for a better view and understanding of the Port of San Juan, Puerto Rico, as a whole system and uses Hurricane Georges (1998), as a representative disruptive event. The model focuses on the impacts of natural disasters at the sub-system level with a final goal of determining the sequence needed to restore an ocean-going port to its pre-event state. This work in progress details model development and outlines steps for using real-world information to assist maritime port manager planning and recommendations for best practices to mitigate disaster damage.

Port of San Juan

Simulating post-dam removal effects of hatchery operations and disease on juvenile Chinook salmon (Oncorhynchus tshawytscha) production in the Lower Klamath River, California

Executive Summary The Federal Energy Regulatory Commission has been considering the approval to breach four dams on lower Klamath River in southern Oregon and northern California. Approval of this application would allow for Strikeouts indicate text deletion hereafter. decommissioning and dam removal, beginning as early as 2023. This action would affect Klamath River salmon ( Oncorhynchus ssp.) populations, a critical food source for federally endangered Southern Resident Killer Whales ( Orcinus orca ). In the long run, reintroduction of salmon populations to the upper Klamath River Basin may increase salmon abundance available to Southern Resident Killer Whales, but in the near term, it is uncertain how changes in hatchery management and disease-caused mortality by the myxosporean parasite Ceratonova shasta will influence abundance of salmon populations entering the ocean. To assess this uncertainty, we used the Stream Salmonid Simulator (S3) to simulate population dynamics of juvenile Chinook salmon ( Oncorhynchus tshawytscha ) for nine different population sources that rear and migrate through the Klamath River. S3 is a spatially explicit population model that runs on a daily time-step and simulates daily growth, survival, and movement of juvenile Chinook salmon from the time of spawning through ocean entry. The key features of this model relevant to this report include (1) a C. shasta disease submodel; (2) a temperature-dependent bioenergetics model that calculates daily growth rates; (3) size-dependent movement; (4) density-dependent dynamics that are influenced by the effect of flow on suitable habitat area; and (5) habitat, river flow, and water temperature specific to each scenario. We constructed and ran four scenarios: two scenarios for dams in place (Dams In) and dams removed (Dams Out), and given these dam-removal conditions, a low- and high-spore scenario for C. shasta . Each scenario was run for nine water years representing a range of conditions from dry to wet. Previously published daily river flows and water temperatures for Dams In and Dams Out provided physical inputs for each scenario. Daily spore concentrations were simulated using a three-part mechanistic model that used river discharge, water temperature, and the prevalence of infection (POI) of hatchery-origin Chinook salmon juveniles with C. shasta in the previous year 3 . We constructed two spore scenarios for each Dams In and Dams Out scenario, a “Low Spore” scenario and a “High Spore” scenario resulting in four scenarios for comparison. Spore scenarios were established by setting the prior-year POI of hatchery fish to 0.15 and 0.75 in the estimation of spore concentrations. Hatchery releases under Dams Out differed from those under the current Dams In scenario. Hatchery releases under the Dams Out scenario were modified to emulate changes in hatchery production that would occur under Dams Out conditions. This included moving hatchery production and releases from Iron Gate Dam to a proposed hatchery at Fall Creek, which would be located about 11 kilometers (km) upstream of Iron Gate Dam. It is anticipated that the Fall Creek hatchery would produce fewer fish at smaller and larger sizes at different release timings. For salmon inputs, we used observed historical abundance of main-stem spawners from brood year 2009 and juvenile salmon entering from tributaries in water year 2010, which represented an average return year for the 2005–18 period. Main-stem spawning was allowed to shift upstream from Iron Gate Dam under the Dams Out scenario. We also included hatchery-origin fish as natural spawners that would have otherwise returned to Iron Gate Hatchery in the first 3 years following dam removal. The S3 model simulated considerably higher total abundance for Dams Out relative to the respective Dams In scenarios, and higher abundance for the Low Spore scenario relative to the High Spore scenario. The difference in abundance between the four combinations of the dam-removal and spore scenarios varied among population groups. For main-stem natural production, juvenile abundance at ocean entry was 2–3 times higher for Dams Out scenarios than for Dams In scenarios, and juvenile abundance for High Spore scenarios was lower than that for the Dams Out Low Spores scenario. For hatchery releases, abundance at ocean entry was similar between Dams In and Dams Out scenarios for most water years, despite lower release sizes from Fall Creek Hatchery under Dams Out. For tributary populations, abundance for the High Spore scenarios was consistently lower than for the Low Spore scenarios, but differences between dam-removal scenarios varied among water years, with Dams Out scenarios having similar or higher abundance than Dams In scenarios, and dry water years having the largest difference between Dams In and Dams Out scenarios. We determined that different factors affected the response of each population group. For main-stem natural production, survival from fry emergence to ocean entry was higher under Dams Out scenarios compared to Dams In scenarios because juveniles emerged later and tended to arrive at the ocean sooner and at larger sizes, causing the population to have less time-dependent in-river mortality. Owing to their late release timing, hatchery populations had high disease-caused mortality in Dams In and Dams Out High Spore scenarios. Furthermore, a high proportion of infected fish (those that would be expected to die at some future point) survived to the ocean. Iron Gate Hatchery fish had lower survival rates than releases from Fall Creek Hatchery because the last mid-June release group from the 2010 Iron Gate Hatchery release incurred nearly total mortality in most water years owing to water temperatures exceeding 24 degrees Celsius. Our analysis shows how the S3 model was able to track different populations and provide insights on how the differential response of each population combined to influence the simulated number of juvenile Chinook salmon arriving at the Pacific Ocean where they become available as a food source for Southern Resident Killer Whales.

California

Variational analysis of drifter positions and model outputs for the reconstruction of surface currents in the central Adriatic during fall 2002

In this paper we present an application of a variational method for the reconstruction of the velocity field in a coastal flow in the central Adriatic Sea, using in situ data from surface drifters and outputs from the ROMS circulation model. The variational approach, previously developed and tested for mesoscale open ocean flows, has been improved and adapted to account for inhomogeneities on boundary current dynamics over complex bathymetry and coastline and for weak Lagrangian persistence in coastal flows. The velocity reconstruction is performed using nine drifter trajectories over 45 d, and a hierarchy of indirect tests is introduced to evaluate the results as the real ocean state is not known. For internal consistency and impact of the analysis, three diagnostics characterizing the particle prediction and transport, in terms of residence times in various zones and export rates from the boundary current toward the interior, show that the reconstruction is quite effective. A qualitative comparison with sea color data from the MODIS satellite images shows that the reconstruction significantly improves the description of the boundary current with respect to the ROMS model first guess, capturing its main features and its exchanges with the interior when sampled by the drifters. Copyright 2008 by the American Geophysical Union.

Journal of Geophysical Research C: Oceans

Changes in surfzone morphodynamics driven by multi-decadal contraction of a large ebb-tidal delta

The impact of multi-decadal, large-scale deflation (76 million m 3 of sediment loss) and contraction (~ 1 km) of a 150 km 2 ebb-tidal delta on hydrodynamics and sediment transport at adjacent Ocean Beach in San Francisco, CA (USA), is examined using a coupled wave and circulation model. The model is forced with representative wave and tidal conditions using recent (2005) and historic (1956) ebb-tidal delta bathymetry data sets. Comparison of the simulations indicates that along north/south trending Ocean Beach the contraction and deflation of the ebb-tidal delta have resulted in significant differences in the flow and sediment dynamics. Between 1956 and 2005 the transverse bar (the shallow attachment point of the ebb-tidal delta to the shoreline) migrated northward ~ 1 km toward the inlet while a persistent alongshore flow and transport divergence point migrated south by ~ 500 m such that these features now overlap. A reduction in tidal prism and sediment supply over the last century has resulted in a net decrease in offshore tidal current-generated sediment transport at the mouth of San Francisco Bay, and a relative increase in onshore-directed wave-driven transport toward the inlet, accounting for the observed contraction of the ebb-tidal delta. Alongshore migration of the transverse bar and alongshore flow divergence have resulted in an increasing proportion of onshore migrating sediment from the ebb-tidal delta to be transported north along the beach in 2005 versus south in 1956. The northerly migrating sediment is then trapped by Pt. Lobos, a rocky headland at the northern extreme of the beach, consistent with the observed shoreline accretion in this area. Conversely, alongshore migration of the transverse bar and divergence point has decreased the sediment supply to southern Ocean Beach, consistent with the observed erosion of the shoreline in this area. This study illustrates the utility of applying a high-resolution coupled circulation-wave model for understanding coastal response to large-scale bathymetric changes over multi-decadal timescales, common to many coastal systems adjacent to urbanized estuaries and watersheds worldwide.

California

Modeling groundwater-level responses to multiple stresses using transfer-function models and wavelet analysis in a coastal aquifer system

In coastal aquifers, dynamic stresses such as climate forcings, groundwater withdrawals, and ocean tidal fluctuations cause nonlinear responses to groundwater levels. Such responses to the stresses impact groundwater resources and related flooding and infrastructure risks at multiple scales. We used time-series models such as transfer-function models and wavelet analysis to quantify the relative contribution of these stresses to groundwater-level fluctuation in wells from the unconfined and confined aquifers in an Atlantic coastal aquifer. Climate forcings, such as precipitation and temperature, explained most of the groundwater-level variation for wells in the unconfined aquifer, whereas groundwater withdrawals were the dominant driver of groundwater levels for wells in the confined aquifer. The impact of groundwater withdrawals also was detected in several wells in the unconfined aquifer. Although the influence of ocean tides on groundwater levels commonly is observed in coastal aquifers, we found that daily groundwater withdrawals can obscure the semi-diurnal coherence signal of the two series. The magnitude of groundwater-level fluctuation that could be explained solely by tides was minor compared to that explained by climate or withdrawal stresses. Transfer-function modeling showed seasonal withdrawals from wells in confined aquifers had a significant, yet heterogeneous influence on groundwater levels in coastal aquifers, which highlights climate and withdrawals as key compounding stresses in coastal hydrology. This study demonstrates the value of time-series approaches to advance characterization of groundwater systems in areas with limited hydrogeologic parameter information.

Virginia