Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Ocean Modeling”

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 1,567 records · Page 87Linked to original sources

A numerical study of wave-driven mean flows and setup dynamics at a coral reef-lagoon system

Two-dimensional mean wave-driven flow and setup dynamics were investigated at a reef-lagoon system at Ningaloo Reef, Western Australia, using the numerical wave-flow model, SWASH. Phase-resolved numerical simulations of the wave and flow fields, validated with highly detailed field observations (including >10 sensors through the energetic surf zone), were used to quantify the main mechanisms that govern the mean momentum balances and resulting mean current and setup patterns, with particular attention to the role of nonlinear wave shapes. Momentum balances from the phase-resolved model indicated that onshore flows near the reef crest were primarily driven by the wave force (dominated by radiation stress gradients) due to intense breaking, whereas the flow over the reef flat and inside the lagoon and channels was primarily driven by a pressure gradient. Wave setup inside the lagoon was primarily controlled by the wave force and bottom stress. The bottom stress reduced the setup on the reef flat and inside the lagoon. Excluding the bottom stress contribution in the setup balance resulted in an over prediction of the wave-setup inside the lagoon by up to 200–370%. The bottom stress was found to be caused by the combined presence of onshore directed wave-driven currents and (nonlinear) waves. Exclusion of the bottom stress contribution from nonlinear wave shapes led to an over prediction of the setup inside the lagoon by approximately 20–40%. The inclusion of the nonlinear wave shape contribution to the bottom stress term was found to be particularly relevant in reef regions that experience a net onshore mass flux over the reef crest.

Journal of Geophysical Research-Oceans↗

The Hardscrabble Creek complex: A newly discovered, mostly buried, Mesoproterozoic mafic-ultramafic pluton in the Wet Mountains, Colorado, USA

The origin of prolific ca. 1.4 Ga ferroan magmatism between the southwestern USA and eastern Canada is enigmatic and has been explained by various models, including extensional, mantle plume, and convergent plate-margin models. Rare mafic plutons are associated with the ferroan plutons, which may help constrain their mantle source and tectonic setting. In the southwestern USA, only two such mafic plutons are known to exist. We present the first evidence for a third, mostly buried, potentially layered, mafic-ultramafic Mesoproterozoic pluton, informally referred to as the Hardscrabble Creek complex, in the central Wet Mountains of Colorado, USA. Recent geophysical data show an elliptical magnetic and gravity high spatially coincident with local gabbroic outcrops. New field and petrographic analyses of these exposed rocks reveal that they consist of ultramafic to mafic cumulates, including orthopyroxenite, olivine norite, norite, and anorthosite. High-precision U-Pb dating of zircon from orthopyroxenite and norite yield weighted mean 206 Pb/ 238 U dates of 1352.36 ± 1.60 Ma and 1352.37 ± 1.71 Ma, respectively. These dates indicate that the complex formed over a narrow timeframe, after the adjacent 1362 ± 7 Ma ferroan San Isabel Granite, and during the waning stages of the regional ca. 1.4 Ga ferroan magmatism. Whole-rock geochemistry and Nd-Sr-Pb isotope compositions of samples from the Hardscrabble Creek complex are similar to those of the San Isabel Granite, suggesting that they were derived from the same or a similar mantle source. The mineral chemistry of the samples is comparable to Proterozoic massif-type anorthosites and related mafic intrusions, indicating that the Hardscrabble Creek complex and San Isabel Granite together represent a rare anorthosite-mangerite-charnockite-granite (AMCG) suite in the southwestern USA. The Hardscrabble Creek complex is unique because it formed ~80 m.y. after the other few mafic plutons in the southwestern USA, and it contains an ultramafic section that is absent from these plutons and rare to the AMCG suite in general. A combination of arc-like whole-rock geochemistry, chondrite uniform reservoir-like Nd-Sr-Pb isotopes, and ocean island basalt (OIB)-like zircon trace element chemistry suggests that the complex was derived from a partial melt of OIB-like mantle and interacted with metasomatically enriched lithospheric mantle. The enriched lithospheric mantle signature, combined with the long ~160 m.y. duration of magmatism in the region, is consistent with a period of protracted convergent tectonism.

Colorado↗

Monitoring the movements of juvenile Pacific Lamprey (Entosphenus tridentatus) in the Yakima River, Washington, using acoustic telemetry, 2019–20

Anthropogenic barriers to main-stem and tributary passage are one of the primary threats associated with declining populations of Pacific Lamprey ( Entosphenus tridentatus ) in the Columbia River Basin. Juvenile lamprey are of special interest because their downstream migration to the ocean may be affected by barriers such as dams or water diversions. Telemetry studies that describe the movement and passage of juvenile lamprey have not been possible until the recent development of a micro-transmitter specifically for use in juvenile lamprey and eels. Through a collaborative research approach, we used these prototype transmitters and acoustic monitoring arrays installed for a juvenile salmon ( Oncorhynchus spp.) migration study to evaluate juvenile lamprey movements in the Yakima River (river kilometer 179 to the river mouth) in 2019 and 2020. We tagged and released 152 juvenile lamprey from April 30 to June 5, 2019, and on June 9, 2020. Lamprey were released 6.9 kilometers (km) upstream from Wapato Dam, 1.2 km upstream from Prosser Dam, and into the canal and tailrace at Prosser Dam. Most tagged lamprey did not initiate downstream movements within the 18 days of tag life, as evidenced by our detections of lamprey in the highest numbers at the first monitoring site downstream from their release site, with limited or no detections at sites farther downstream. There was no evidence of missed detections (lamprey detected at a downstream site without corresponding detections upstream). Overall detections of tagged lamprey were low: 27.0 percent in 2019 and 48.0 percent in 2020. River flows were less than the 10-year average during the monitoring period and water temperatures were variable. Lamprey arrived at detections sites predominantly during periods of darkness (85.3–96.6 percent) following daytime releases. Travel rates through the study area ranged from 0.2 to 45.3 kilometers per day, and lamprey generally remained at each detection station for less than about 20 minutes. Groups of lamprey released together generally had similar travel rates with a small number of fish that moved more quickly or slowly than the remainder of the group. In addition to monitoring the migration and behavior of juvenile lamprey, we also assessed some assumptions of survival models (determining downstream drift of purposely killed fish and empirically measuring transmitter operating life) to benefit future evaluations focused on migration survival of juvenile lamprey.

Washington↗

Lidar-revised geologic map of the Uncas 7.5' quadrangle, Clallam and Jefferson Counties, Washington

In 2000 and 2001, the Puget Sound Lidar Consortium obtained 1 pulse/m 2 lidar data for about 65 percent of the Uncas 7.5' quadrangle. For a brief description of LIDAR (LIght Detection And Ranging) and this data acquisition program, see Haugerud and others (2003). This map combines geologic interpretation (mostly by Haugerud and Tabor) of the 6-ft (2-m) lidar-derived digital elevation model (DEM) with the geology depicted on the Preliminary Geologic Map of the Uncas 7.5' Quadrangle, Clallam and Jefferson Counties, Washington, by Peter J. Haeussler and others (1999). The Uncas quadrangle in the northeastern Olympic Peninsula covers the transition from the accreted terranes of the Olympic Mountains on the west to the Tertiary and Quaternary basin fills of the Puget Lowland to the east. Elevations in the map area range from sea level at Port Discovery to 4,116 ft (1,255 m) on the flank of the Olympic Mountains to the southwest. Previous geologic mapping within and marginal to the Uncas quadrangle includes reports by Cady and others (1972), Brown and others (1960), Tabor and Cady (1978a), Yount and Gower (1991), and Yount and others (1993). Paleontologic and stratigraphic investigations by University of Washington graduate students (Allison, 1959; Thoms, 1959; Sherman, 1960; Hamlin, 1962; Spencer, 1984) also encompass parts of the Uncas quadrangle. Haeussler and Wells mapped in February 1998, following preliminary mapping by Yount and Gower in 1976 and 1979. The description of surficial map units follows Yount and others (1993) and Booth and Waldron (2004). Bedrock map units are modified from Yount and Gower (1991) and Spencer (1984). We used the geologic time scale of Gradstein and others (2005). The Uncas quadrangle lies in the forearc of the Cascadia subduction zone, about 6.25 mi (10 km) east of the Cascadia accretionary complex exposed in the core of the Olympic Mountains (Tabor and Cady, 1978b). Underthrusting of the accretionary complex beneath the forearc uplifted and tilted eastward the Coast Range basalt basement and overlying marginal basin strata, which comprise most of the rocks of the Uncas quadrangle. The Eocene submarine and subaerial tholeiitic basalt of the Crescent Formation on the Olympic Peninsula is thought to be the exposed mafic basement of the Coast Range, which was considered by Snavely and others (1968) to be an oceanic terrane accreted to the margin in Eocene time. In this interpretation, the Coast Range basalt terrane may have originated as an oceanic plateau or by oblique marginal rifting, but its subsequent emplacement history was complex (Wells and others, 1984). Babcock and others (1992) and Haeussler and others (2003) favor the interpretation that the basalts were the product of an oceanic spreading center interacting with the continental margin. Regardless of their origin, onlapping strata in southern Oregon indicate that the Coast Range basalts were attached to North America by 50 Ma; but on southern Vancouver Island, where the terrane-bounding Leech River Fault is exposed, Brandon and Vance (1992) concluded that suturing to North America occurred in the broad interval between 42 and 24 Ma. After emplacement of the Coast Range basalt terrane, the Cascadia accretionary wedge developed by frontal accretion and underplating (Tabor and Cady, 1978b; Clowes and others, 1987). Domal uplift of the part of the accretionary complex beneath the Olympic Mountains occurred after ~18 Ma (Brandon and others, 1998). Continental and alpine glaciation during Quaternary time reshaped the uplifted rocks of the Olympic Mountains.

Washington↗

Preface—Evaluating the response of critical zone processes to human impacts with sediment source fingerprinting

1) Background: Critical Zone Processes in the Anthropocene The Earth’s Critical Zone encompasses a suite of interconnected processes in the near-surface lithosphere, pedosphere, biosphere, atmosphere, and hydrosphere (Brantley et al., 2007; Lin, 2010) (Fig. 1). Processes and interactions both within and between these various Critical Zone components supports life-sustaining ecosystem services and resources that establish the foundation for humanity (NRC, 2001). This includes the formation production of fertile soils, flourishing vegetation, productive rivers, lakes and oceans, and our life-sustaining atmosphere (Gaillardet, 2014; Guo and Lin, 2016). Rapid population growth, land use intensification, and global environmental change are disturbing many of these fundamental Critical Zone processes. More than half of the Earth’s terrestrial surface is now impacted by anthropogenic activities (e.g., clearing, grazing, plowing, mining, and logging) (Hooke et al., 2012; Richter and Mobley, 2009). These changes are so widespread and pervasive that the great acceleration of socioeconomic development that occurred around 1950 (Fig. 2) has been recommended to delineate the dawn of the Anthropocene (Waters et al., 2016). Although the utility of adopting and delineating the Anthropocene as the current epoch is subject to debate (Crutzen, 2002; Ruddiman et al., 2015; Smith and Zeder, 2013), the concept effectively highlights both the nature and the extent of our global impact on Earth’s Critical Zone. Soil forming processes and ecosystem services provided by the pedosphere are central to the Critical Zone (Banwart et al., 2011; Lin, 2010). Many of these processes have been disturbed by the agricultural intensification that coincided with the great acceleration resulting in unsustainable land use practices now outpacing soil formation processes (Brantley et al., 2007). As agricultural landscapes now cover an area equivalent to what was scoured during the last glacial maximum (Amundson et al., 2007), the broad-scale intensification of anthropogenic activities has resulted in significant on- and off-site impacts. On-site, soil loss has resulted in decreases in soil fertility and agricultural yields (Ladha et al., 2009) threatening the ability to feed the world’s growing population (Brantley et al., 2007). Off-site, the excess delivery of particulate matter downstream is degrading riverine, lacustrine, and estuarine ecosystems (Bilotta and Brazier, 2008; Clark, 1985; Owens et al., 2005). The challenge, as noted by Brantley et al., (2007), is that despite our society having over 10,000 years of experience working with soils, our conceptual and quantitative models remain inadequate at predicting Critical Zone dynamics under current conditions. Notwithstanding growing pressure for improved environmental management, we still have a limited capacity to predict changes in the Critical Zone in response to anthropogenic activities owing to the multiple spatial and temporal scales at which these complex processes and feedbacks are manifest. As river basin systems are impacted by many of these processes, a deep understanding of soil-sediment continuum dynamics may provide a valuable framework for evaluating the disturbance response of Critical Zone processes. Understanding these processes may also provide land and resource managers with the information necessary to manage both the on-site and off-site effects of accelerated soil erosion.

Journal of Soils and Sediments↗

Mechanisms of wave‐driven water level variability on reef‐fringed coastlines

Wave‐driven water level variability (and runup at the shoreline) is a significant cause of coastal flooding induced by storms. Wave runup is challenging to predict, particularly along tropical coral reef‐fringed coastlines due to the steep bathymetric profiles and large bottom roughness generated by reef organisms, which can violate assumptions in conventional models applied to open sandy coastlines. To investigate the mechanisms of wave‐driven water level variability on a reef‐fringed coastline, we performed a set of laboratory flume experiments on an along‐shore uniform bathymetric profile with and without bottom roughness. Wave setup and waves at frequencies lower than the incident sea‐swell forcing (infragravity waves) were found to be the dominant components of runup. These infragravity waves were positively correlated with offshore wave groups, signifying they were generated in the surf zone by the oscillation of the breakpoint. On the reef flat and at the shoreline, the low‐frequency waves formed a standing wave pattern with energy concentrated at the natural frequencies of the reef flat, indicating resonant amplification. Roughness elements used in the flume to mimic large reef bottom roughness reduced low frequency motions on the reef flat and reduced wave run up by 30% on average, compared to the runs over a smooth bed. These results provide insight into sea‐swell and infragravity wave transformation and wave setup dynamics on steep‐sloped coastlines, and the effect that future losses of reef bottom roughness may have on coastal flooding along reef‐fringed coasts.

Journal of Geophysical Research C: Oceans↗

Status Assessment of Laysan and Black-Footed Albatrosses, North Pacific Ocean, 1923-2005

Over the past century, Laysan (Phoebastria immutabilis) and black-footed (Phoebastria nigripes) albatrosses have been subjected to high rates of mortality and disturbance at the breeding colonies and at sea. Populations were greatly reduced and many colonies were extirpated around the turn of the 20th century as a result of feather hunting. Populations were recovering when military occupation of several breeding islands during World War II led to new population declines at these islands and additional colony extirpations. At sea, thousands of Laysan and black-footed albatrosses were killed each year in high-seas driftnet fisheries, especially from 1978 until the fisheries were banned in 1992. Through the 1990s, there was a growing awareness of the large numbers of albatrosses that were being killed in longline fisheries. During the 1990s, other anthropogenic factors, such as predation by non-native mammals and exposure to contaminants, also were documented to reduce productivity or increase mortality. In response to the growing concerns over the impacts of these threats on albatross populations, the U.S. Fish and Wildlife Service contracted with the U.S. Geological Survey to conduct an assessment of Laysan and black-footed albatross populations. This assessment includes a review of the taxonomy, legal status, geographic distribution, natural history, habitat requirements, threats, and monitoring and management activities for these two species. The second part of the assessment is an analysis of population status and trends from 1923 to 2005. Laysan and black-footed albatrosses forage throughout the North Pacific Ocean and nest on tropical and sub-tropical oceanic islands from Mexico to Japan. As of 2005, 21 islands support breeding colonies of one or both species. The core breeding range is the Hawaiian Islands, where greater than 99 percent of the World's Laysan albatrosses and greater than 95 percent of the black-footed albatrosses nest on the small islands and atolls of the Northwestern Hawaiian Islands. These islands are all protected as part of the Papahanaumokuakea Marine National Monument. Albatrosses are long-lived seabirds with deferred maturity, low fecundity, and high rates of adult survival. Their life history characteristics make populations especially vulnerable, to small increases in adult mortality. The primary threats to Laysan and black-footed albatrosses include interactions with commercial fisheries, predation by introduced mammals, reduced reproductive output due to contaminants, nesting habitat loss and degradation due to human development and invasive plant species, and potential loss and degradation of habitat due to climate change and sea-level rise. Incidental mortality (bycatch) in commercial fisheries is the greatest anthropogenic source of mortality (post-fledging) for both species. We found that longline fishing effort prior to the 1980s was greater than previously estimated and a very significant source of mortality. Regulations to minimize and monitor albatross mortality have been enacted in most U.S. and Canadian longline fisheries, but monitoring of bycatch rates and regulations to minimize seabird mortality are extremely limited in the much larger multinational longline fleets. Management to address threats at the breeding colonies is ongoing and includes eradication or control of non-native species, habitat management, and abatement programs to reduce impacts of contaminants. Effective long-term conservation and management of the Laysan and black-footed albatrosses require management and monitoring at the breeding colonies and at sea and continued assessment of population status and trends. We evaluated the status and trends of Laysan and black-footed albatross populations using linear regression, population viability analysis (PVA), and age-structured matrix models. Analyses were predominantly based on nest-count data gathered at French Frigate Shoals, Laysan Island, and Midw

Scientific Investigations Report↗

Direct measurements of mean Reynolds stress and ripple roughness in the presence of energetic forcing by surface waves

Direct covariance observations of the mean flow Reynolds stress and sonar images of the seafloor collected on a wave‐exposed inner continental shelf demonstrate that the drag exerted by the seabed on the overlying flow is consistent with boundary layer models for wave‐current interaction, provided that the orientation and anisotropy of the bed roughness are appropriately quantified. Large spatial and temporal variations in drag result from nonequilibrium ripple dynamics, ripple anisotropy, and the orientation of the ripples relative to the current. At a location in coarse sand characterized by large two‐dimensional orbital ripples, the observed drag shows a strong dependence on the relative orientation of the mean current to the ripple crests. At a contrasting location in fine sand, where more isotropic sub‐orbital ripples are observed, the sensitivity of the current to the orientation of the ripples is reduced. Further, at the coarse site under conditions when the currents are parallel to the ripple crests and the wave orbital diameter is smaller than the wavelength of the relic orbital ripples, the flow becomes hydraulically smooth. This transition is not observed at the fine site, where the observed wave orbital diameter is always greater than the wavelength of the observed sub‐orbital ripples. Paradoxically, the dominant along‐shelf flows often experience lower drag at the coarse site than at the fine site, despite the larger ripples, highlighting the complex dynamics controlling drag in wave‐exposed environments with heterogeneous roughness.

Journal of Geophysical Research C: Oceans↗

Evolution of a trench-slope basin within the Cascadia subduction margin: The Neogene Humboldt Basin, California

The Neogene Humboldt (Eel River) Basin is located along the north-eastern margin of the Pacific Ocean within the Cascadia subduction zone. This sedimentary basin originated near the base of the accretionary prism in post-Eocene time. Subduction processes since that time have elevated strata in the south-eastern portion of the basin above sea level. High-resolution chronostratigraphic data from the onshore portion of the Humboldt Basin enable correlation of time-equivalent lithofacies across the palaeomargin, reconstruction of slope-basin evolution, and preliminary delineation of climatic and tectonic influence on lithological variation. Emergent basin fill is divided into five lithofacies which clearly document shoaling of the inner trench slope from deep-water environments in early Miocene time to paralic environments in Pleistocene time. The oldest strata consist of hemipelagic mudstones and minor debris-flow breccias deposited in a deep-water setting during elevated sea level. These strata are overlain by glauconite-rich, fine-grained turbidites which heralded an increasing influx of terrigenous detritus. Water depths shoaled earlier in the eastern basin area as the palaeoshoreline prograded seaward. Turbidite deposition ceased in the eastern basin area at about 2-2 Ma, whereas 22 km to the west, turbidite deposition continued until about 1-8 Ma. Lithofacies at the western study site change abruptly across a middle Pleistocene unconformity from outer shelf to paralic deposits. In the east, a more complete Pleistocene section records transition from outer to inner shelf, beach and fluvial environments. The Humboldt Basin lithofacies sequence is overprinted by eustatic control of sediment source. Comparison of sediment character with palaeoceanographic conditions indicates dominance of hemipelagic facies during periods of elevated sea level in the middle Miocene and early Pliocene when depocentres were isolated from terrigenous sediment. Glauconite-rich facies were mobilized from an upper slope setting following these periods of elevated sea level and redeposited in a deep-marine environment. Pleistocene shoreline lithofacies display glacio-esutatic control of depositional environment by recording several cycles of nearshore to fluvial progressions. General models of accretionary prism behaviour and trench-slope basin evolution are compatible with the overall coarsening-upward lithofacies sequence filling the Humboldt Basin. Early structural barriers precluded deposition of terrigenous material except from locally derived debris flows; subsequent shoaling and burial of deactivated thrust-folds enabled turbidity flows to reach the basin floor. However, late-stage tectonism apparently controlled the onset of coarse-grained deposition in this sequence. Significant sand-rich turbidite deposition began in the middle Pliocene, synchronous with tectonic uplift of the southern basin margin. Conversely, cessation of turbidite deposition in the eastern basin area in latest Pliocene time was synchronous with growth of anticlinal structures which again blocked widespread dispersal of turbidity flows. This middle Pliocene to Holocene period of crustal shortening is synchronous with continued reduction in spreading rate along the southern Juan de Fuca ridge, and probably reflects partial coupling between the subducting lithosphere and the overlying accretionary prism.

California↗

Sora rail studies on the Patuxent River, Maryland

The freshwater marshes of the tidal Patuxent River are well known for their annual fall concentration of migrant soras ( Porzana carolina ) and were formerly the most famous rail hunting grounds in the Chesapeake Bay region. Because of concern over the apparent long-term decline in number of soras and the decline in the quality of the Patuxent marshes, especially the loss of wild rice ( Zizania aquatica ), the Maryland National-Capital Park and Planning Commission (MNCPPC), co-steward of the Jug Bay National Estuarine Research Reserve, sponsored rail-related research beginning in 1987. Past efforts focused on developing efficient trapping techniques, age and sex criteria, and monitoring body mass dynamics. Noted progress was made in developing digital playback systems and trap improvements to enhance sora captures. These improvements increased capture success by over an order of magnitude and resulted in capture of 2,315 soras and 276 Virginia rails ( Rallus limicola ) in the 5 year period, 1993-97. Although these methods demonstrate the efficacy of banding large numbers of soras on migration and possibly winter concentration areas, captures at the Patuxent River site have been 70-90% hatching-year birds and recoveries and recaptures have been virtually nonexistent. With the present effort, this outcome precludes population parameter estimation using traditional capture-recapture or recovery model methodologies. In 1996, studies were initiated to employ radio telemetry methods to investigate length of stay, habitat use, survival, and migration characteristics of fall migrant soras. These studies are ongoing and will be continued through 1998 with a grant from the U.S. Fish and Wildlife Service's Webless Migratory Game Bird Research Program and support from the U.S. Geological Survey's Patuxent Wildlife Research Center. Supplemental funding has also been provided by MNCPPC, FWS Region 5, the Maryland Ornithological Society, Quail Unlimited, and Prince Georges Community College. During 1996-97 we developed a successful radio transmitter attachment technique to secure 1.8g radio transmitters over the synsacrum of migrant soras. We modified Rappole and Tipton's (1991) leg-loop attachment method by addition of a waist loop to prevent soras from slipping transmitters over their short tails. Thin gauge (0.6 mm) elastic thread proved ideal for attachment and allowed for girth expansion associated with fattening during stopover. Sixty instrumented soras have been monitored in two years of study from early September until early November. Only a single mortality was recorded and 41 (68%) were confIrmed and another 13 (total 90%) were believed to have migrated from the study area. Only a single bird slipped a radio transmitter. Most birds demonstrated a sedentary nature in the marsh throughout stopover. Average length of stay was 44 days in 1997 (n=29) with peak departure occurring 20-24 October. Departing migrants were detected using a receiver/data-logger monitoring system placed 4 miles down river from the study site. Thirty-six of 37 (97%) soras departed in a 2-hour window of time, beginning 1 hour after sunset. Departure was synchronized with cold fronts on clear, starlit nights. Twenty-five soras were monitored on migration from 8 km to as far away as 770 km. Findings indicate migration flight speeds of 40 to 50 mph and a direct southward orientation from the study site at least until contact with the Atlantic Ocean west of Cape Lookout, North Carolina. We also attempted to monitor passage of migrant soras radio tagged at Iroquois National Wildlife Refuge in western New York by biologist Soch Lor. None were recorded passing a monitoring unit at Haldeman Island in the middle of the Susquehanna River 10 miles north of Harrisburg, Pennsylvania. This further corroborates our findings that when migrating overland, soras move in a direct southward orientation and are not following land features such as major rivers

Book chapter↗

Mitigating flood risks in urban estuaries: Tidal dynamics, shoreline hardening, nature-based solutions, and floodgates in San Francisco Bay

Hydrodynamic models are valuable tools for understanding the primary factors influencing daily and peak water levels and for guiding discussions on potential adaptation strategies for managing flood risk in coastal areas. This analysis uses the Delft3D San Francisco Bay-Delta Community Model to simulate water levels and incorporates the effects of a number of adaptation measures in the urban San Francisco Bay estuary, California. In particular, we examine the influence of shoreline hardening, nature-based solutions, and subregional floodgates on regional water levels. The result shows that under present conditions, tidal amplification is responsible for generating a wide distribution of extreme water levels across San Francisco Bay. Tidal amplification is found to decrease under sea level rise, thereby producing a relative damping effect on extremes. A comparison of different shoreline scenarios demonstrates that hard frontal shorelines result in higher tidal amplification, whereas restored (soft) shorelines lower amplification. The current shoreline configuration has both hard and soft characteristics and results in an intermediate tidal response. In some areas, wetland restoration reduces extreme water levels by as much as 20 cm, whereas hard-shoreline addition elevates them by as much as 10 cm for 1.5 m of sea level rise. Furthermore, local floodgates can significantly reduce high water levels without major adverse effects elsewhere in San Francisco Bay. These findings point toward the justification for a range of adaptive measures across political boundaries, weighing hard and soft options in addressing the mounting danger of sea level rise.

California↗

Carbon and oxygen isotopic composition of coal and carbon dioxide derived from laboratory coal combustion: A preliminary study

The concentration of carbon dioxide (CO 2 ) in the atmosphere has dramatically increased from the start of the industrial revolution in the mid-1700s to present levels exceeding 400 ppm. Carbon dioxide derived from fossil fuel combustion is a greenhouse gas and a major contributor to on-going climate change. Carbon and oxygen stable isotope geochemistry is a useful tool to help model and predict the contributions of anthropogenic sources of CO 2 in the global carbon cycle. Surprisingly few studies have addressed the carbon and oxygen isotopic composition of CO 2 derived from coal combustion. The goal of this study is to document the relationships between the carbon and oxygen isotope signatures of coal and signatures of the CO 2 produced from laboratory coal combustion in atmospheric conditions. Six coal samples were selected that represent various geologic ages (Carboniferous to Tertiary) and coal ranks (lignite to bituminous). Duplicate splits of the six coal samples were ignited and partially combusted in the laboratory at atmospheric conditions. The resulting coal-combustion gases were collected and the molecular composition of the collected gases and isotopic analyses of δ 13 C of CO 2 , δ 13 C of CH 4 , and δ 18 O of CO 2 were analysed by a commercial laboratory. Splits (~ 1 g) of the un-combusted dried ground coal samples were analyzed for δ 13 C and δ 18 O by the U.S. Geological Survey Reston Stable Isotope Laboratory. The major findings of this preliminary work indicate that the isotopic signatures of δ 13 C (relative to the Vienna Pee Dee Belemnite scale, VPDB) of CO 2 resulting from coal combustion are similar to the δ 13 C VPDB signature of the bulk coal (− 28.46 to − 23.86 ‰) and are not similar to atmospheric δ 13 C VPDB of CO 2 (~ − 8 ‰, see http://www.esrl.noaa.gov/gmd/outreach/isotopes/c13tellsus.html ). The δ 18 O values of bulk coal are strongly correlated to the coal dry ash yields and appear to have little or no influence on the δ 18 O values of CO 2 resulting from coal combustion in open atmospheric conditions. There is a wide range of δ 13 C values of coal reported in the literature and the δ 13 C values from this study generally follow reported ranges for higher plants over geologic time. The values of δ 18 O (relative to Vienna Standard Mean Ocean Water) of CO 2 derived from atmospheric combustion of coal and other high-carbon fuels (peat and coal) range from + 19.03 to + 27.03‰ and are similar to atmospheric oxygen δ 18 O VSMOW values which average + 23.8‰. Further work is needed on a broader set of samples to better define the relationships between coal composition and combustion-derived gases.

International Journal of Coal Geology↗

Concentrations and Loads of Organic Compounds and Trace Elements in Tributaries to Newark and Raritan Bays, New Jersey

A study was undertaken to determine the concentrations and loads of sediment and chemicals delivered to Newark and Raritan Bays by five major tributaries: the Raritan, Passaic, Rahway, Elizabeth, and Hackensack Rivers. This study was initiated by the State of New Jersey as Study I-C of the New Jersey Toxics Reduction Workplan for the New York-New Jersey Harbor, working under the NY-NJ Harbor Estuary Program (HEP) Contaminant Assessment and Reduction Program (CARP). The CARP is a comprehensive effort to evaluate the levels and sources of toxic contaminants to the tributaries and estuarine areas of the NY-NJ Harbor, including Newark and Raritan Bays. The Raritan and Passaic Rivers are large rivers (mean daily discharges of 1,189 and 1,132 cubic feet per second (ft3/s), respectively), that drain large, mixed rural/urban basins. The Elizabeth and Rahway Rivers are small rivers (mean daily discharges of 25.9 and 49.1 ft3/s, respectively) that drain small, highly urbanized and industrialized basins. The Hackensack River drains a small, mixed rural/urban basin, and its flow is highly controlled by an upstream reservoir (mean daily discharge of 90.4 ft3/s). These rivers flow into urbanized estuaries and ultimately, to the Atlantic Ocean. Each of these tributaries were sampled during two to four storm events, and twice each during low-flow discharge conditions. Samples were collected using automated equipment installed at stations adjacent to U.S. Geological Survey streamflow-gaging stations near the heads-of-tide of these rivers. Large-volume (greater than 50 liters of water and a target of 1 gram of sediment), flow-weighted composite samples were collected for chemical analysis using filtration to collect suspended particulates and exchange resin (XAD-2) to sequester dissolved contaminants. Composite whole-water samples were collected for dissolved polycyclic aromatic hydrocarbons (PAH) and for trace element analysis. Additional discrete grab samples were collected throughout each event for trace-element analysis, and multiple samples were collected for suspended sediment (SS), particulate carbon (POC), and dissolved organic carbon (DOC) analysis. The suspended sediment and exchange resin were analyzed for 114 polychlorinated biphenyls (PCBs, by US EPA method 1668A, modified), seven 2,3,7,8-substituted chlorinated dibenzo-p-dioxins (CDD) and 10 dibenzo-p-difurans (CDF) (by US EPA method 1613), 24 PAHs (by low-resolution isotope dilution/mass-spectral methods), 27 organo-chlorine pesticides (OCPs) (by high resolution isotope dilution/mass-spectral methods), and the trace elements mercury (Hg), methyl-mercury (MeHg), lead (Pb), and cadmium (Cd). Isotope dilution methods using gas chromatography and high-and low-resolution mass spectral (GC/MS) detection were used to accurately identify and quantify organic compounds in the sediment and water phases. Trace elements were measured using inductively coupled plasma-mass spectrometry and cold-vapor atomic fluorescence spectrometry methods. The loads of sediment, carbon, and chemicals were calculated for each storm and low-flow event sampled. Because only a few storm events were sampled, yearly loads of sediment were calculated from rating curves developed using historical SS and POC data. The average annual loads of sediment and carbon were calculated for the period 1975-2000, along with the loads for the selected water years being modeled as part of the New York New Jersey Harbor Estuary Program CARP. Comparison of loads calculated using the rating curve method to loads measured during the sampled storm events indicated that the rating curve method likely underpredicts annual loads. Average annual loads of suspended sediment in the tributaries were estimated to be 395,000 kilograms per year (kg/yr) in the Hackensack River, 417,000 kg/yr in the Elizabeth River, 882,000 kg/yr in the Rahway River, 22,700,000 kg/yr in the Passaic River, and 93,100,000 kg/yr in the Raritan River. Averag

Scientific Investigations Report↗

Across-shelf sediment transport: Interactions between suspended and bed sediment

We use a two-dimensional, time-dependent sediment-transport model to quantify across-shelf transport, deposition, and sorting during wave-driven resuspension events characteristic of those that dominate sediment transport on many continental shelves. Decreases in wave-orbital velocities as water depth increases, and the resulting cross-shelf gradient in bed shear stress favor a net offshore transport of sediment. On wide, flat shelves (slopes ∼0.1%percnt;), these gradients are low, and the depth to which the seabed is reworked depends mainly on bottom shear stress and local sediment availability. On narrow, steep shelves (slopes ∼0.5%percnt;), however, the gradient in bottom stress generates significant cross-shelf suspended sediment flux gradients that create regions of net erosion and deposition. While the magnitude of waves generally determines the water depth to which sediment can be resuspended, erosional and depositional patterns on narrow shelves are sensitive to cross-shelf gradients in wave energy, nonlocal sediment availability, and the direction and magnitude of the cross-shelf current. During energetic waves, cross-shelf divergence of suspended sediment flux can create a coarsened, erosional area on the inner shelf that abuts a region of fine-grained sediment deposition on the mid-to-outer shelf. If currents are strongly shoreward, however, flux divergence leads to erosion over the entire shelf.

Journal of Geophysical Research C: Oceans↗

Impact of Hurricane Irma on coral reef sediment redistribution at Looe Key Reef, Florida, USA

Understanding event-driven sediment transport in coral reef environments is essential to assessing impacts on reef species, habitats, restoration, and mitigation, yet a global knowledge gap remains due to limited quantitative studies. Hurricane Irma made landfall in the Lower Florida Keys with sustained 209 km h −1 winds and waves greater than 8 m on 10 September 2017, directly impacting the Florida Reef Tract (FRT) and providing an opportunity to perform a unique comprehensive, quantitative assessment of its impact on coral reef structure and sediment redistribution. We used lidar and multibeam derived digital elevation models (DEMs) collected before and after the passing of Hurricane Irma over a 15.98 km 2 area along the lower FRT including Looe Key Reef to quantify changes in seafloor elevation, volume, and structure due to storm impacts. Elevation change was calculated at over 4 million point locations across 10 habitat types within this study area for two time periods using data collected (1) approximately 1 year before the passing of Irma and 3 to 6 months following the storm's impact as well as (2) 3 to 6 months after and up to 16.5 months after the storm. Elevation change data were then used to generate triangulated irregular network (TIN) models in ArcMap to calculate changes in seafloor volume during each time period. Our results indicate that Hurricane Irma was primarily a depositional event that increased mean seafloor elevation and volume at this study site by 0.34 m and up to 5.4 Mm 3 , respectively. Sediment was transported primarily west-southwest (WSW) and downslope, modifying geomorphic seafloor features including the migration of sand waves and rubble fields, formation of scour marks in shallow seagrass habitats, and burial of seagrass and coral-dominated habitats. Approximately 16.5 months after Hurricane Irma (during a 13-month period between 2017 and 2019), net erosion was observed across all habitats with mean elevation change of −0.15 m and net volume change up to −2.46 Mm 3 . Rates of elevation change during this post-storm period were 1 to 2 orders of magnitude greater than decadal and multi-decadal rates of change in the same location, and changes showed erosion of approximately 50 % of sediment deposited during the storm event as seafloor sediment distribution began to re-equilibrate to non-storm sea-state conditions. Our results suggest that higher-resolution elevation change data collected over seasonal and annual time periods could enhance characterization and understanding of short-term and long-term rates and processes of seafloor change.

Florida↗

Reexamining the Honolulu Volcanics: Hawai‘i's classic case of rejuvenation volcanism

Rejuvenated volcanism is a worldwide phenomenon occurring on many oceanic islands in all of the major ocean basins. This plume-related volcanism follows the main edifice-building stage after a hiatus of variable duration (e.g. 0.6–2 Myrs in Hawai'i). The Honolulu Volcanics (HV), the classic case of rejuvenated volcanism, involved monogenetic eruptions from at least 48 vent areas. Previous studies inferred these vents were aligned along 3 to 11 rifts oriented orthogonal to the propagation direction of the Hawaiian plume. HV basalts are known for having high MgO contents (greater than 10 wt %) and upper mantle xenoliths. Thus, HV magmas are assumed to be relatively primitive and to have ascended rapidly (less than 1 day) through the crust. However, new analyses of olivine cores in basalts from 24 HV vents are mostly too low in forsterite content (74–86 mol %) to be in equilibrium with mantle melts. Olivine and clinopyroxene in HV basalts commonly show reverse zoning indicating magma mixing prior to eruption. These results are inconsistent with the rapid ascent of HV magmas directly from their mantle source. Many of the HV magmas underwent storage (probably in the lower crust or uppermost mantle), crystal fractionation and magma mixing prior to eruption. New 40 Ar/ 39 Ar dates were determined for 11 HV lavas to evaluate their eruptive history. These ages, 80 to 685 ka, combined with our previous and other 40 Ar/ 39 Ar ages for HV lavas reveal long gaps (greater than 50 kyr) between some eruptions. Our comprehensive, whole-rock major and trace element database (63 XRF analyses, 57 ICPMS analyses) of basalts from 37 vents show remarkable compositional diversity with no obvious spatial pattern or temporal trends. The two most recent eruptive sequences have the greatest diversity (basanite and melilitite compositions). HV basanites show systematic trace element trends that may reflect mixing of multiple source components. The nephelinites and melilitites require a complex source history that may have involved residual accessory minerals during mantle melting and a metasomatic component that was not carbonatitic. The new ages and geochemical data show eruptions along most of the previously proposed rift systems were unrelated (except for the Koko Rift). Therefore, geodynamic models that relate HV volcanism to these rift systems are invalid. Lava volumes for two HV eruptions were estimated at 0.11 and 0.23 km 3 using surface mapping and water well data. Similar size, recent monogenetic eruptions in Auckland, New Zealand, were inferred to have lasted several months. Thus, if another HV eruption were to occur, which is possible given the long hiatus between eruptions, it would be extremely disruptive for the nearly 1 million residents of Honolulu. None of the existing geodynamic models fully explain the age duration, volumes and the locations of Hawai'i's rejuvenated volcanism. Thus, the cause of this secondary volcanism remains enigmatic.

Hawaii↗

Tidal triggering of low frequency earthquakes near Parkfield, California: Implications for fault mechanics within the brittle-ductile transition

Studies of nonvolcanic tremor (NVT) have established the significant impact of small stress perturbations on NVT generation. Here we analyze the influence of the solid earth and ocean tides on a catalog of ∼550,000 low frequency earthquakes (LFEs) distributed along a 150 km section of the San Andreas Fault centered at Parkfield. LFE families are identified in the NVT data on the basis of waveform similarity and are thought to represent small, effectively co-located earthquakes occurring on brittle asperities on an otherwise aseismic fault at depths of 16 to 30 km. We calculate the sensitivity of each of these 88 LFE families to the tidally induced right-lateral shear stress (RLSS), fault-normal stress (FNS), and their time derivatives and use the hypocentral locations of each family to map the spatial variability of this sensitivity. LFE occurrence is most strongly modulated by fluctuations in shear stress, with the majority of families demonstrating a correlation with RLSS at the 99% confidence level or above. Producing the observed LFE rate modulation in response to shear stress perturbations requires low effective stress in the LFE source region. There are substantial lateral and vertical variations in tidal shear stress sensitivity, which we interpret to reflect spatial variation in source region properties, such as friction and pore fluid pressure. Additionally, we find that highly episodic, shallow LFE families are generally less correlated with tidal stresses than their deeper, continuously active counterparts. The majority of families have weaker or insignificant correlation with positive (tensile) FNS. Two groups of families demonstrate a stronger correlation with fault-normal tension to the north and with compression to the south of Parkfield. The families that correlate with fault-normal clamping coincide with a releasing right bend in the surface fault trace and the LFE locations, suggesting that the San Andreas remains localized and contiguous down to near the base of the crust. The deep families that have high sensitivity to both shear and tensile normal stress perturbations may be indicative of an increase in effective fault contact area with depth. Synthesizing our observations with those of other LFE-hosting localities will help to develop a comprehensive understanding of transient fault slip below the “seismogenic zone” by providing constraints on parameters in physical models of slow slip and LFEs.

California↗

Geologic and isostatic gravity map of the Nenana basin area, central Alaska

The Nenana Basin area is a prospective petroleum province in central Alaska, and this geologic and isostatic gravity map is part of a petroleum resource assessment of the area. The geology was compiled from published sources as shown on the index map (map sheet). Map units are organized and presented according to the scheme of lithotectonic terranes proposed by Jones and others (1987) and Silberling and Jones (1984); we recognize, however, that this terrane scheme is controversial and likely to be revised in the future. In some cases, we combined certain terranes because we were unable to match the terrane boundaries given by Jones and others (1987) and Silberling and Jones (1984) with specific faults shown on existing geologic maps. Postaccretion cover deposits represent overlap assemblages that depositionally overlie accreted terranes. Plutonic igneous rocks shown on this map include several plutons that are clearly postaccretionary, based on isotopic ages and (or) field relations. It is possible that some of the plutons predate accretion, but this has not been demonstrated. According to Jones and others (1982), the terranes in the area of our map were assembled during late Mesozoic or earliest Cenozoic time. The gravity contours are derived from data used in earlier compilations that are supplemented by some National Oceanic and Atmospheric Administration data along the Alaska Pipeline level line. The earlier compilations were used for simple Bouguer maps, prepared primarily by non-digital methods, and are superseded by this map. The present map is the result of digital processing that includes the 1967 Geodetic Reference System, the IGSN-71 datum, digital terrain corrections, and conversion to isostatic gravity so that geologic structures on the margin of the Alaska Range are more clearly portrayed. Computation procedures are described in part by Barnes (1972, 1984), Jachens and Roberts (1981), and Barnes and others (1994). The calculations used a crustal density of 2.67 g/cm 3 , a density contrast at the base of the isostatic root of 0.4 g/cm 3 , and a root thickness at sea level of 25 km. The distribution of data within the map area is uneven and locally controls the shape of the computer-generated contours. Altimetry was used for most of the elevation control and its inconsistency is responsible for many of the small contour irregularities. Ninety percent of the measurements are estimated to have an accuracy of about 1.5 mgal or about a quarter of the 5 mgal contour interval. Data collection and analysis were assisted by R.V. Allen, R.C. Jachens, M.A. Fisher, T.R. Bruns, J.G. Blank, J.W. Bader, Z.C. Valin, J.W. Cady, R.L. Morin, and P.V. Woodward. The most promising area for petroleum exploration is a prominent 25 mgal isostatic gravity low north of Nenana (T. 2 S., R. 8 W.). This gravity low probably corresponds to the deepest part of a sedimentary basin filled by Cenozoic strata that includes nonmarine fluvial and lacustrine deposits of the Eocene to Miocene Usibelli Group. Smaller gravity lows are associated with outcrops of these sedimentary rocks north of Suntrana (T. 12 S., R. 6-9 W.) and Sable Pass (T. 16 S., R. 11 W.). A broad low on the north flank of the Alaska Range east of the Wood River (T. 10 S., R. 1 E.) indicates another basin under the Tanana lowland that extends eastward off the map area towards Delta Junction, where its presence was confirmed by both gravity and seismic data. Gravity modelling suggests that the base of the Usibelli Group in the area north of Nenana (T. 2 S., R. 8 W.) is about 3,000 to 3,350 m beneath the ground surface. Organic geochemical studies indicate that mudstones and coals in the Usibelli Group are potential sources of petroleum; calculations based on borehole temperatures suggest that, in the area of the gravity low, these rocks may have been buried deeply enough to generate oil and gas. Two exploratory wells, the Union Nenana No. 1 and the ARCO Totek Hills No. 11, were drilled some distance away from the gravity low in areas where the Usibelli Group is thin. Mudlogs show that both wells were dry holes that bottomed in schist and had gas shows associated with coal beds in the Usibelli Group, but no reported signs of oil.

Alaska↗