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.

SourceResearch connectionFrom yearThrough year

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

At least 1,513 records · Page 84Linked to original sources

USGS1997

Organic carbon balance and net ecosystem metabolism in Chesapeake Bay

The major fluxes of organic carbon associated with physical transport and biological metabolism were compiled, analyzed and compared for the mainstem portion of Chesapeake Bay (USA). In addition, 5 independent methods were used to calculate the annual mean net ecosystem metabolism (NEM = production - respiration) for the integrated Bay. These methods, which employed biogeochemical models, nutrient mass-balances anti summation of individual organic carbon fluxes, yielded remarkably similar estimates, with a mean NEM of +50 g C m-2 yr-1 (?? SE = 751, which is approximately 8% of the estimated annual average gross primary production. These calculations suggest a strong cross-sectional pattern in NEM throughout the Bay, wherein net heterotrophic metabolism prevails in the pelagic zones of the main channel, while net autotrophy occurs in the littoral zones which flank the deeper central area. For computational purposes, the estuary was separated into 3 regions along the land-sea gradient: (1) the oligohaline Upper Bay (11% of total area); (2) the mesohaline Mid Bay (36% of area); and (3) the polyhaline Lower Bay (53% of area). A distinct regional trend in NEM was observed along this salinity gradient, with net here(atrophy (NEM = 87 g C m-2 yr-1) in the Upper Bay, balanced metabolism in the Mid Bay and net autotrophy (NEM = +92 g C m-2 yr-1) in the Lower Bay. As a consequence of overall net autotrophy, the ratio of dissolved inorganic nitrogen (DIN) to total organic nitrogen (TON) changed from DIN:TON = 5.1 for riverine inputs to DIN:TON = 0.04 for water exported to the ocean. A striking feature of this organic C mass-balance was the relative dominance of biologically mediated metabolic fluxes compared to physical transport fluxes. The overall ratio of physical TOC inputs (1) to biotic primary production (P) was 0.08 for the whole estuary, but varied dramatically from 2.3 in the Upper Bay to 0.03 in the Mid and Lower Bay regions. Similarly, ecosystem respiration was some 6-fold higher than the sum of all physical carbon sinks. This general negative correspondence between I:P ratio and NEM, which occurred among Bay regions, was also evident in data available for organic C fluxes in other coastal ecosystems. An inverse relationship between NEM and P, postulated in a previous study, did not apply to Chesapeake Bay, and closer examination of available data revealed the importance of the loading ratio of DIN:TOC as a key control on coastal NEM. It is proposed here that the general global trend of coastal eutrophication will lead to increasing values of NEM in estuaries worldwide. The management implications of this trend are complex, involving both increased potential fisheries harvest and decreased demersal habitat.

Marine Ecology Progress Series↗

USGS1991

A comparison of Sr-Nd-Pb isotopes in young and old continental lithospheric mantle: Patagonia and eastern China

It is commonly accepted that beneath the continental crust lies a keel of lithospheric mantle, which extends 50–200 kilometres downward to a transition zone into the asthenosphere. The chemical and physical properties of this reservoir are best known through studies of the basalts and xenoliths that provide samples of the subcrustal mantle. Although sharing many characteristics with oceanic island basalts, some continental basalts become increasingly distinct isotopically as crustal age increases, strongly supporting a permanent association between crust and mantle. Consequently, the distinctive trace element and isotope composition of the lithospheric mantle is able to give important clues to its origin and evolution. The mantle under newly‐created crust is typified by a radiogenic isotope variability that emphasizes the materials from which the continental lithosphere is assembled. Old lithospheric mantle, on the other hand, exhibits more evolved isotopic patterns that attest to the existence of long‐lived, chemically complex systems. A comparison of the Pb, Sr and Nd isotopes in alkalic to sub‐alkalic basalt derived from Phanerozoic (Patagonia) and Middle Archaean to Early Proterozoic (eastern China) subcrustal mantle is useful for identifying ‘end‐member’ components of the lithosphere. One component, having an isotopic composition close to PREMA, either continues to evolve virtually unchanged after incorporation into the lithosphere or is, itself, a relatively new addition even to old lithosphere. Another component, beginning with the isotopic composition of BSE, undergoes significant reduction in U/Pb and Sm/Nd (but not Rb/Sr) upon incorporation into the lithosphere and, with time, shows an increasingly retarded evolution of 206 Pb/ 204 Pb and negative ε Nd‐values approaching the isotopic composition of EMI. Five models are discussed that relate the isotopic composition of the continental lithospheric mantle to that of other parts of the terrestrial system, which may be involved in its origin and evolution. The potential locations of the contributing components and the mechanisms and timing of their assembly into lithosphere are considered. Current knowledge, however, does not allow us to distinguish unequivocally among the various scenarios for the creation and evolution of this reservoir.

Australian Journal of Earth Sciences↗

USGS2017

Emulation of long-term changes in global climate: application to the late Pliocene and future

Multi-millennial transient simulations of climate changes have a range of important applications, such as for investigating key geologic events and transitions for which high-resolution palaeoenvironmental proxy data are available, or for projecting the long-term impacts of future climate evolution on the performance of geological repositories for the disposal of radioactive wastes. However, due to the high computational requirements of current fully coupled general circulation models (GCMs), long-term simulations can generally only be performed with less complex models and/or at lower spatial resolution. In this study, we present novel longterm “continuous” projections of climate evolution based on the output from GCMs, via the use of a statistical emulator. The emulator is calibrated using ensembles of GCM simulations, which have varying orbital configurations and atmospheric CO2 concentrations and enables a variety of investigations of long-term climate change to be conducted, which would not be possible with other modelling techniques on the same temporal and spatial scales. To illustrate the potential applications, we apply the emulator to the late Pliocene (by modelling surface air temperature – SAT), comparing its results with palaeo-proxy data for a number of global sites, and to the next 200 kyr (thousand years) (by modelling SAT and precipitation). A range of CO2 scenarios are prescribed for each period. During the late Pliocene, we find that emulated SAT varies on an approximately precessional timescale, with evidence of increased obliquity response at times. A comparison of atmospheric CO2 concentration for this period, estimated using the proxy sea surface temperature (SST) data from different sites and emulator results, finds that relatively similar CO2 concentrations are estimated based on sites at lower latitudes, whereas higher-latitude sites show larger discrepancies. In our second illustrative application, spanning the next 200 kyr into the future, we find that SAT oscillations appear to be primarily influenced by obliquity for the first ∼ 120 kyr, whilst eccentricity is relatively low, after which precession plays a more dominant role. Conversely, variations in precipitation over the entire period demonstrate a strong precessional signal. Overall, we find that the emulator provides a useful and powerful tool for rapidly simulating the long-term evolution of climate, both past and future, due to its relatively high spatial resolution and relatively low computational cost. However, there are uncertainties associated with the approach used, including the inability of the emulator to capture deviations from a quasi-stationary response to the forcing, such as transient adjustments of the deep-ocean temperature and circulation, in addition to its limited range of fixed ice sheet configurations and its requirement for prescribed atmospheric CO2 concentrations.

Climate of the Past↗

USGS2020

Plate boundary localization, slip-rates and rupture segmentation of the Queen Charlotte Fault based on submarine tectonic geomorphology

Linking fault behavior over many earthquake cycles to individual earthquake behavior is a primary goal in tectonic geomorphology, particularly across an entire plate boundary. Here, we examine the 1150-km-long, right-lateral Queen Charlotte-Fairweather fault system using comprehensive multibeam bathymetry data acquired along the Queen Charlotte Fault (QCF) offshore southeastern Alaska and western British Columbia. Fine-scale analysis of tectonic geomorphology allowed us to identify and reconstruct 184 strike-slip piercing points over a 630 km stretch of the QCF. Age constraints from glacial recession and offshore sedimentation patterns yield a consistent slip-rate of ∼50–57 mm/yr since ∼17–12 ka, the fastest rate for a continent-ocean strike-slip fault on Earth. These slip-rates equal or exceed estimates of Pacific-North America (PA-NA) relative motion from global plate reconstructions, indicating that PA-NA motion is highly localized. The QCF cuts the seafloor along a narrow and unusually straight trace for its entire length and multiple fault traces are observed only at local step-overs. The geometry and behavior of the QCF over many earthquake cycles is simple and typical of mature faults with relatively homogeneous stress fields. Since the QCF is the primary PA-NA plate boundary, we used the trace of the QCF to define the small circle path for relative plate motion and computed the associated Euler pole. Predicted along-strike obliquity variations based on the new pole agree with observed tectonic geomorphology and suggest that previous global plate reconstructions overestimated the degree of oblique convergence along the QCF. We also find that subtle, long-wavelength (75–150 km) bends and discrete step-overs appear to define the endpoints of M>7 earthquakes, suggesting that obliquity and resultant fault geometry may control rupture segmentation and asperity development. Lastly, the agreement between predicted obliquity and tectonic geomorphology along the entire length of QCF compelled a reevaluation of regional tectonic models. In the north, the eastern Yakatat Terrane appears to be translating northwest with the Pacific plate, and slip transferred from the QCF to the Fairweather Fault results in ∼20 mm/yr of convergence along the southern St. Elias mountains. In the south, we predict a reduced rate of convergence along the QCF west of Haida Gwaii (∼5–6 mm/yr of shortening, on average) relative to previous studies. Our results support a model for transpression and strike-slip partitioning along the edge of a hot and weak Pacific Plate, leading to crustal thickening and growth of the Queen Charlotte Terrace to the west of Haida Gwaii.

Alaska, British Columbia↗

USGS2014

Dissolved-solids sources, loads, yields, and concentrations in streams of the conterminous United States

Recent studies have shown that excessive dissolved-solids concentrations in water can have adverse effects on the environment and on agricultural, domestic, municipal, and industrial water users. Such effects motivated the U.S. Geological Survey’s National Water Quality Assessment Program to develop a SPAtially-Referenced Regression on Watershed Attributes (SPARROW) model that has improved the understanding of sources, loads, yields, and concentrations of dissolved solids in streams of the conterminous United States. Using the SPARROW model, long-term mean annual dissolved-solids loads from 2,560 water-quality monitoring stations were statistically related to several spatial datasets that are surrogates for dissolved-solids sources and land-to-water delivery processes. Specifically, sources in the model included variables representing geologic materials, road deicers, urban lands, cultivated lands, and pasture lands. Transport of dissolved solids from these sources was modulated by land-to-water delivery variables that represent precipitation, streamflow, soil, vegetation, terrain, population, irrigation, and artificial drainage characteristics. Where appropriate, the load estimates, source variables, and transport variables were statistically adjusted to represent conditions for the base year 2000. The nonlinear least-squares estimated SPARROW model was used to predict long-term mean annual conditions for dissolved-solids sources, loads, yields, and concentrations in a digital hydrologic network representing nearly 66,000 stream reaches and their corresponding incremental catchments that drain the Nation. Nationwide, the predominant source of dissolved solids yielded from incremental catchments and delivered to local streams is geologic materials in 89 percent of the catchments, road deicers in 5 percent of the catchments, pasture lands in 3 percent of the catchments, urban lands in 2 percent of the catchments, and cultivated lands in 1 percent of the catchments. Whereas incremental catchments with dissolved solids that originated predominantly from geologic sources or from urban lands are found across much of the Nation, incremental catchments with dissolved solids yields that originated predominantly from road deicers are largely found in the Northeast, and incremental catchments with dissolved solids that originated predominantly from cultivated or pasture lands are largely found in the West. The total amount of dissolved solids delivered to the Nation’s streams is 271.9 million metric tons (Mt) annually, of which 194.2 million Mt (71.4%) come from geologic sources, 37.7 million Mt (13.9%) come from road deicers, 18.2 million Mt (6.7%) come from pasture lands, 13.9 million Mt (5.1%) come from urban lands, and 7.9 million Mt (2.9%) come from cultivated lands. Nationwide, the median incremental-catchment yield delivered to local streams is 26 metric tons per year per square kilometer [(Mt/yr)/km 2 ]. Ten percent of the incremental catchments yield less than 4 (Mt/yr)/km 2 , and 10 percent yield more than 90 (Mt/yr)/km 2 . Incremental-catchment yields greater than 50 (Mt/yr)/km 2 mostly occur along the northern part of the West Coast and in a crescent shaped band south of the Great Lakes. For example, the median incremental-catchment yield is 81 (Mt/yr)/km 2 for the Great Lakes, 78 (Mt/yr)/km 2 for the Ohio, and 74 (Mt/yr)/km 2 for the Upper Mississippi water-resources regions. Incremental-catchment yields less than 10 (Mt/yr)/km 2 mostly occur in a wide band across the arid lowland of the interior West that excludes areas along the coast and the extensive, higher mountain ranges. For example, the median incremental-catchment yield is 3 (Mt/yr)/km 2 for the Lower Colorado, 5 (Mt/yr)/km 2 for the Rio Grande, and 8 (Mt/yr)/km 2 for the Great Basin water-resources regions. Predicted incremental loads were cascaded down through the reach network, with loads accumulating from reach to reach. For most stream reaches, the entire incremental load of dissolved solids delivered to the reach was transported to either the ocean or to one of the large streams flowing along the U.S. international boundary without losses occurring along the way. The exceptions to this include streams in the southwestern part of the country, such as the Colorado River, Rio Grande, and streams of internally drained drainages in the Great Basin, where dissolved-solids loads decreased through streamflow diversion for off-stream use, or by infiltration through the streambed. Long-term mean annual flow-weighted concentrations were derived from the predicted accumulated-load and stream-discharge data. Widespread low concentrations, generally less than 100 milligrams per liter (mg/L), occur in many reaches of the New England, South Atlantic-Gulf, and Pacific Northwest water-resources regions as a result of moderate dissolved-solids yields and high runoff rates. Widespread moderate concentrations, generally between 100 and 500 mg/L, occur in many reaches of the Great Lakes, Ohio, and Upper Mississippi River water-resources regions. Whereas dissolved-solids yields are generally high in these regions, runoff rates are also high, which helps moderate concentrations in these regions. Widespread higher concentrations, generally greater than 500 mg/L, occur across a belt of reaches that extends almost continuously from Canada to Mexico in the Midwest, cutting through the Souris-Red-Rainy, Missouri, Arkansas-White-Red, Texas-Gulf, and Rio Grande water-resources regions. Although dissolved-solids yields are moderate to low in these areas, low runoff rates result in the high concentrations for these areas. In 12.6 percent of the Nation’s stream reaches, predicted concentrations of dissolved solids exceed 500 mg/L, the U.S. Environmental Protection Agency’s secondary, nonenforceable drinking water standard. While this standard provides a metric for evaluating predicted concentrations in the context of drinking-water supplies, it should be noted that it only applies to drinking water actually served to customers by water utilities, and it does not apply to all stream reaches in the Nation nor does it apply during times when water is not being withdrawn for use. Exceedance of 500 mg/L is more pronounced in certain water-resources regions than others. For example, about half of the reaches in the Souris-Red-Rainy region have concentrations predicted to exceed 500 mg/L, and between 25 and 37 percent of the reaches in the Missouri, Arkansas-White-Red, Texas-Gulf, Rio Grande, and Lower Colorado regions are predicted to exceed 500 mg/L. Development of stream-load data for use in the SPARROW model also provided long-term temporal trend information in dissolved-solids concentrations at the monitoring stations for their period of record, which was constrained between 1980 and 2009. For the 2,560 monitoring stations used in this study, long-term trends in flow-adjusted dissolved-solids concentrations increased over time at 23 percent of the stations, decreased at 18 percent of the stations, and did not change over time at 59 percent of the stations. Long-term trends show a strong regional spatial pattern where from the western parts of the Great Plains to the West Coast, concentrations mostly either did not change or decreased over time, and from the eastern parts of the Great Plains to the East Coast, concentrations mostly either did not change or increased over time. Results from the trend analysis and from the SPARROW model indicate that, compared to monitoring stations with no trends or decreasing trends, stations with increasing trends are associated with a smaller percentage of the predicted dissolved-solids load originating from geologic sources, and a larger percentage originating from urban lands and road deicers. Conversely, compared to stations with increasing trends or no trends, stations with decreasing trends have a larger percentage of the predicted dissolved-solids load originating from geologic sources and a smaller percentage originating from urban lands and road deicers. Stations with decreasing trends also have larger percentages of predicted dissolved-solids load originating from cultivated lands and pasture lands, compared to stations with increasing trends or no trends.

Scientific Investigations Report↗

USGS2004

Long-Term Oceanographic Observations in Western Massachusetts Bay Offshore of Boston, Massachusetts: Data Report for 1989-2002

This data report presents long-term oceanographic observations made in western Massachusetts Bay at two locations: (1) 42 deg 22.6' N., 70 deg 47.0' W. (Site A, 33 m water depth) from December 1989 through December 2002 (figure 1), and (2) 42 deg 9.8' N., 70 deg 38.4' W. (Site B, 21 m water depth) from October 1997 through December 2002. Site A is approximately 1 km south of the new ocean outfall that began discharging treated sewage effluent from the Boston metropolitan area into Massachusetts Bay on September 6, 2000. These long-term oceanographic observations have been collected by the U.S. Geological Survey (USGS) in partnership with the Massachusetts Water Resources Authority (MWRA) and with logistical support from the U.S. Coast Guard (USCG - http://www.uscg.mil). This report presents time series data through December 2002, updating a similar report that presented data through December 2000 (Butman and others, 2002). In addition, the Statistics and Mean Flow sections include some new plots and tables and the format of the report has been streamlined by combining yearly figures into single .pdfs. Figure 1 (PDF format) The long-term measurements are planned to continue at least through 2005. The long-term oceanographic observations at Sites A and B are part of a USGS study designed to understand the transport and long-term fate of sediments and associated contaminants in the Massachusetts bays. (See http://woodshole.er.usgs.gov/project-pages/bostonharbor/ and Butman and Bothner, 1997.) The long-term observations document seasonal and inter-annual changes in currents, hydrography, and suspended-matter concentration in western Massachusetts Bay, and the importance of infrequent catastrophic events, such as major storms or hurricanes, in sediment resuspension and transport. They also provide observations for testing numerical models of circulation. This data report presents a description of the field program and instrumentation, an overview of the data through summary plots and statistics, and the data in NetCDF and ASCII format for the period December 1989 through December 2002 for Site A and October 1997 through December 2002 for Site B. The objective of this report is to make the data available in digital form and to provide summary plots and statistics to facilitate browsing of the long-term data set.

Data Series↗

USGS2021

A survey of storm-induced seaward-transport features observed during the 2019 and 2020 hurricane seasons

Hurricanes are known to play a critical role in reshaping coastlines, but often only impacts on the open ocean coast are considered, ignoring seaward-directed forces and responses. The identification of subaerial evidence for storm-induced seaward transport is a critical step towards understanding its impact on coastal resiliency. The visual features, found in the National Oceanic and Atmospheric Administration, National Geodetic Survey Emergency Response Imagery (ERI) collected after recent hurricanes on the U.S. East Atlantic and Gulf of Mexico coasts, include scours and channelized erosion, but also deposition on the shoreface or in the nearshore as deltas and fans of various sizes. We catalog all available ERI and describe recently formed features found on the North Core Banks, North Carolina, after Hurricane Dorian (2019); the Carolina coasts after Hurricane Isaias (2020); the Isles Dernieres, Louisiana, after Hurricane Zeta (2020); and the southwest coast of Louisiana, after Hurricanes Laura and Delta (2020). Hundreds of features were identified over nearly 200 km of coastline with the density of features exceeding 20 per km in some areas. Individual features range in size from 5 m to 500 m in the alongshore, with similar dimensions in the cross-shore direction, including the formation or reactivation of outlets. The extensive occurrence of these storm-induced return-flow and seawardflow morphologic features demonstrates that their role in coastal evolution and resilience may be more prominent than previously thought. Based on these observations we propose clarifying terms for return- and seaward-flow features to distinguish them from more frequently documented landward-flow features and advocate for their inclusion in coastal change hazards classification schemes and coastal evolution morphodynamic models.

southeast United States↗

USGS1985

Origin and evolution of the alkalic ultramafic rocks in the Coyote Peak diatreme, Humboldt County, California

Instrumental-neutron-activation analyses are reported for two uncontaminated rocks, a phlogopite-rich clot, and two contaminated rocks from the Coyote Peak diatreme, northwestern California. These data, combined with Nd, Sr, and Pb isotopic evidence, have been modeled to a multi-stage evolution for the uncontaminated rocks. Fertile mantle material (refractory elements 2.5× chondritic abundances; Rb / Sr = 0.029 by weight) was depleted about 900 m.y. ago by congruent melting and removal of ~4% basaltic liquid; this depleted residue provided the source rock from which the Coyote Peak magma was ultimately derived. About 66 m.y. ago, the depleted mantle residue was incongruently melted in the presence of H 2 O and CO 2 at a total pressure > 26 kb to yield ~0.5% of a Si-poor, Ca-rich melt. This melt then metasomatized depleted garnet-free harzburgite in the upper mantle at about 26 kb to produce a rock similar to phlogopite-bearing wehrlite. About 29 m.y. ago, this rock was subjected to an increase in pressure to >26 kb and incongruently melted to give ~0.5% of a second-stage melt resembling olivine melilitite in composition. Enroute to the surface, about 28% olivine and 2% titanomagnetite were lost from the highly fluid melt. Coarse-grained phlogopite-rich clots in the uncontaminated rocks apparently crystallized from a latestage liquid derived from the uncontaminated melt. Contaminated rocks appear to be the result of partial assimilation of, and dilution by, ~14% Franciscan graywacke country rock. The diatreme was emplaced near a converging plate margin where young hot oceanic mantle and crust of the Juan de Fuca plate was probably subducting obliquely beneath a thin lip of the North American plate. The unusual chemistry of the rocks may be the result of this complex tectonic setting which could also have included local strike-slip and extensional environments within the two plates pierced by the diatreme.

Geochimica et Cosmochimica Acta↗

USGS1990

Geometry of the September 1971 eruptive fissure at Kilauea volcano, Hawaii

A three-dimensional model has been used to estimate the location and dimensions of the eruptive fissure for the 24-29 September 1971 eruption along the southwest rift zone of Kilauea volcano, Hawaii. The model is an inclined rectangular sheet embedded in an elastic half-space with constant displacement on the plane of the sheet. The set of "best" model parameters suggests that the sheet is vertical, extends from a depth of about 2 km to the surface, and has a length of about 14 km. Because this sheet intersects the surface where eruptive vents and extensive ground cracking formed during the eruption, this sheet probably represents the conduit for erupted lava. The amount of displacement perpendicular to the sheet is about 1.9 m, in the middle range of values measured for the amount of opening across the September 1971 eruptive fissure. The thickness of the eruptive fissure associated with the January 1983 east rift zone eruption was determined in an earlier paper to be 3.6 m, about twice the thickness determined here for the September 1971 eruption. Because the lengths (12 km for 1983 and 14 km for 1971) and heights (about 2 km) of the sheet models derived for the January 1983 and September 1971 rift zone eruptions are nearly identical, the greater thickness for the January 1983 eruptive fissure implies that the magma pressure was about a factor of two greater to form the January 1983 eruptive fissure. Because the September 1971 and January 1983 eruptive fissures extent to depths of only a few kilometers, the region of greatest compressive stress produced along the volcano's flank by either of these eruptive fissures would also be within a few kilometers of the surface. Previous work has shown that rift eruptions and intrusions contribute to the buildup of compressive stress along Kilauea's south flank and that this buildup is released by increased seismicity along the south flank. Because south flank earthquakes occur at significantly greater depths, i.e., from 5 to 13 km, than the vertical extent of the 1971 and 1983 eruptiv fissures, the depth of emplacement of these eruptive fissures cannot be the main factor in controlling the hypocentral depths of south flank earthquakes. Two possible explanations for the occurrence of south flank earthquakes in the depth range of 5-13 km are (1) a deeper pressure source, possibly related to deeper magma storage within the rift zone, and (2) a lowstrength region located between 5 and 13 km beneath Kilauea's south flank, possibly at the interface between oceanic sediments and the base of the Hawaiian volcanics. ?? 1990 Springer-Verlag.

Bulletin of Volcanology↗

USGS2018

Executive summary. In Second State of the Carbon Cycle Report (SOCCR2): A Sustained Assessment Report

Central to life on Earth, carbon is essential to the molecular makeup of all living things and plays a key role in regulating global climate. To understand carbon’s role in these processes, researchers measure and evaluate carbon stocks and fluxes. A stock is the quantity of carbon contained in a pool or reservoir in the Earth system (e.g., carbon in forest trees), and a flux is the direction and rate of carbon’s transfer between pools (e.g., the movement of carbon from the atmosphere into forest trees during photosynthesis). This document, the Second State of the Carbon Cycle Report (SOCCR2), examines the patterns of carbon stocks and fluxes—collectively called the “carbon cycle.” Emphasis is given to these patterns in specific sectors (e.g., agriculture and energy) and ecosystems (e.g., forests and coastal waters) and to the response of the carbon cycle to human activity. The purpose of SOCCR2 is to assess the current state of the North American carbon cycle and to present recent advances in understanding the factors that influence it. Concentrating on North America—Canada, the United States, and Mexico—the report describes carbon cycling for air, land, inland waters (streams, rivers, lakes, and reservoirs), and coastal waters (see Figure ES.1, p. 23). The questions framing the publication A U.S. Carbon Cycle Science Plan (Michalak et al., 2011) inspired development of three slightly modified questions that guide SOCCR2’s content and focus on North America in a global context: How have natural processes and human actions affected the global carbon cycle on land, in the atmosphere, in the ocean and other aquatic systems, and at ecosystem interfaces (e.g., coastal, wetland, and urban-rural)? How have socioeconomic trends affected atmospheric levels of the primary carbon-containing gases, carbon dioxide (CO2) and methane (CH4)? How have species, ecosystems, natural resources, and human systems been impacted by increasing greenhouse gas (GHG) concentrations, associated changes in climate, and carbon management decisions and practices? SOCCR2 synthesizes the most recent understanding of carbon cycling in North America, assessing new carbon cycle findings and information, the state of knowledge regarding core methods used to study the carbon cycle, and future research needed to best inform carbon management and policy options. Focusing on scientific developments in the decade since the First State of the Carbon Cycle Report (SOCCR1; CCSP 2007), SOCCR2 summarizes the past, current, and projected state of carbon sources, sinks, and natural processes, as well as contributions by human activities. In addition to CO2 and CH4, the report sometimes discusses nitrous oxide (N2O), a GHG associated with activities and processes that affect fluxes of carbon gases.1 SOCCR2 also describes improvements in analysis tools; developments in decision support; and new insights into ecosystem carbon cycling, human causes of changes in the carbon cycle, and social science perspectives on carbon. Since publication of SOCCR1, coordinated research from agencies in the three North American countries has enabled innovative observational, analytical, and modeling capabilities to further advance understanding of the North American carbon cycle (see Appendix D: Carbon Measurement Approaches and Accounting Frameworks, p. 834). Some of the report’s main conclusions, based on the Key Findings of each chapter, are highlighted in Box ES.1, Main Findings of SOCCR2, p. 24.

Report↗

USGS2020

Challenges in quantifying air-water carbon dioxide flux using estuarine water quality data: Case study for Chesapeake Bay

Estuaries play an uncertain but potentially important role in the global carbon cycle via CO 2 outgassing. The uncertainty mainly stems from the paucity of studies that document the full spatial and temporal variability of estuarine surface water partial pressure of carbon dioxide ( p CO 2 ). Here, we explore the potential of utilizing the abundance of pH data from historical water quality monitoring programs to fill the data void via a case study of the mainstem Chesapeake Bay (eastern United States). We calculate p CO 2 and the air-water CO 2 flux at monthly resolution from 1998 to 2018 from tidal fresh to polyhaline waters, paying special attention to the error estimation. The biggest error is due to the pH measurement error, and errors due to the gas transfer velocity, temporal sampling, the alkalinity mixing model, and the organic alkalinity estimation are 72%, 27%, 15%, and 5%, respectively, of the error due to pH. Seasonal, interannual, and spatial variability in the air-water flux and surface p CO 2 is high, and a correlation analysis with oxygen reveals that this variability is driven largely by biological processes. Averaged over 1998–2018, the mainstem bay is a weak net source of CO 2 to the atmosphere of 1.2 (1.1, 1.4) mol m −2 yr −1 (best estimate and 95% confidence interval). Our findings suggest that the abundance of historical pH measurements in estuaries around the globe should be mined in order to constrain the large spatial and temporal variability of the CO 2 exchange between estuaries and the atmosphere.

Maryland, Virginia↗

USGS2016

Characterizing supraglacial meltwater channel hydraulics on the Greenland Ice Sheet from in situ observations

Supraglacial rivers on the Greenland ice sheet (GrIS) transport large volumes of surface meltwater toward the ocean, yet have received relatively little direct research. This study presents field observations of channel width, depth, velocity, and water surface slope for nine supraglacial channels on the southwestern GrIS collected between 23 July and 20 August, 2012. Field sites are located up to 74 km inland and span 494-1485 m elevation, and contain measured discharges larger than any previous in situ study: from 0.006 to 23.12 m 3 /s in channels 0.20 to 20.62 m wide. All channels were deeply incised with near vertical banks, and hydraulic geometry results indicate that supraglacial channels primarily accommodate greater discharges by increasing velocity. Smaller streams had steeper water surface slopes (0.74-8.83%) than typical in terrestrial settings, yielding correspondingly high velocities (0.40-2.60 m/s) and Froude numbers (0.45-3.11) with supercritical flow observed in 54% of measurements. Derived Manning's n values were larger and more variable than anticipated from channels of uniform substrate, ranging from 0.009 to 0.154 with a mean value of 0.035 +/- 0.027 despite the absence of sediment, debris, or other roughness elements. Ubiquitous micro-depressions in shallow sections of the channel bed may explain some of these roughness values. However, we find that other, unobserved sources of flow resistance likely contributed to these elevated n values: future work should explicitly consider additional sources of flow resistance beyond bed roughness in supraglacial channels. We conclude that hydraulic modelling for these channels must allow for both sub- and supercritical flow, and most importantly must refrain from assuming that all ice-substrate channels exhibit similar hydraulic behavior, especially for Froude numbers and Manning's n. Finally, this study highlights that further theoretical and empirical work on supraglacial channel hydraulics is necessary before broad scale understanding of ice sheet hydrology can be achieved. This article is protected by copyright. All rights reserved.

Earth Surface Processes and Landforms↗

USGS2004

Three-dimensional distribution of gas hydrate beneath southern Hydrate Ridge: Constraints from ODP Leg 204

Large uncertainties about the energy resource potential and role in global climate change of gas hydrates result from uncertainty about how much hydrate is contained in marine sediments. During Leg 204 of the Ocean Drilling Program (ODP) to the accretionary complex of the Cascadia subduction zone, we sampled the gas hydrate stability zone (GHSZ) from the seafloor to its base in contrasting geological settings defined by a 3D seismic survey. By integrating results from different methods, including several new techniques developed for Leg 204, we overcome the problem of spatial under-sampling inherent in robust methods traditionally used for estimating the hydrate content of cores and obtain a high-resolution, quantitative estimate of the total amount and spatial variability of gas hydrate in this structural system. We conclude that high gas hydrate content (30-40% of pore space or 20-26% of total volume) is restricted to the upper tens of meters below the seafloor near the summit of the structure, where vigorous fluid venting occurs. Elsewhere, the average gas hydrate content of the sediments in the gas hydrate stability zone is generally <2% of the pore space, although this estimate may increase by a factor of 2 when patchy zones of locally higher gas hydrate content are included in the calculation. These patchy zones are structurally and stratigraphically controlled, contain up to 20% hydrate in the pore space when averaged over zones ???10 m thick, and may occur in up to ???20% of the region imaged by 3D seismic data. This heterogeneous gas hydrate distribution is an important constraint on models of gas hydrate formation in marine sediments and the response of the sediments to tectonic and environmental change. ?? 2004 Published by Elsevier B.V.

Earth and Planetary Science Letters↗

USGS2004

Twentieth century demographic changes in cirio and cardón in Baja California, Mexico

Aim&ensp; Our purpose was to discern long-term large-scale patterns of survivorship and recruitment of two dominant, charismatic and useful desert plants, cirio [Fouquieriaceae; Fouquieria columnaris (Kell.) Kell. ex Curran] and cardón [Cactaceae; Pachycereus pringlei (S.Watson) Britton & Rose], and to evaluate the effects of physical and human factors. Location&ensp; The study included 77 sites distributed over c. 13,200 km 2 in the northern Vizcaíno Region of the Sonoran Desert, state of Baja California, México. Method&ensp; Censuses used identified individuals ( n = 3780 cirio and 2246 cardón) and were based on repeat photography. Time intervals between censuses ranged from 29 to 96 years, ending in 1996–2002. Earlier repetition (1963) of old photographs at 14 sites allowed analysis of temporal variation. The effect of elapsed time was modelled with Weibull functions for survivorship and linear functions for recruitment. To explain the residual variations, we tested categorical and continuous variables related to substrate, geomorphology, climate and geography. Human impacts were weighed with historical and quantitative analysis. Results&ensp; Elapsed time accounted for 45% of the variation of survivorship in cirio but only 35% in cardón. The fitted Weibull curves were used to estimate longevity: all-size cohorts of 1000 individuals would die out in 388 (cirio) and 723 years (cardón). Recruitment in cirio was linearly related to time ( r 2 = 0.71) and averaged < 0.003 year −1 but was less than mortality at 90% of the sites. For cardón, recruitment averaged > 0.005 year −1 , was linearly related to time ( r 2 = 0.66) and was less than mortality at only 33% of the sites. Vital rates were not strongly correlated between mid- and late-century but were mostly similar to the long-term regional rates. Removing the factor of elapsed time, the residual variations showed that survivorship was greater on sites protected from winds for cirio and was less on very coarse soils for cardón. Recruitment increased with latitude in cirio and was greater on southern exposures for both species. The residual variations were not clearly affected, at a regional scale, by other factors thought to be important, including elevation, distance to the Pacific Ocean, geology, slope gradient, soil stability, older vs. young alluvial soils, and soil Ca : Mg and Na : K ratios. Human impacts have been sporadic and heterogeneous but locally strong; our quantitative indices of accessibility did not show regionally significant effects. Blowdown by hurricanes is a sporadic cause of major mortality for cirio but not necessarily for cardón. Main conclusions&ensp; At our scales, effects of time outweigh those of substrate, and human impact was scant. Large patterns were pronounced: cirio experienced a prolonged, widespread decline in the 1900s, while cardón fluctuated in different directions and degree among local populations. Cirio was found to be inherently much shorter-lived than cardón. We also suggest that recruitment in cirio was low due to long periods with relatively dry winters that affected the entire region, while spatial heterogeneity of cardón recruitment seemed more related to the variation in summer rains.

Journal of Biogeography↗

USGS1985

Accretion tectonics and crustal structure in Alaska

The entire width of the North American Cordillera in Alaska is made up of "suspect terranes". Pre-Late Cretaceous paleogeography is poorly constrained and the ultimate origins of the many fragments which make up the state are unclear. The Prince William and Chugach terranes accreted since Late Cretaceous time and represent the collapse of much of the northeast Pacific Ocean swept into what today is southern Alaska. Greater Wrangellia, a composite terrane now dispersed into fragments scattered from Idaho to southern Alaska, apparently accreted into Alaska in Late Cretaceous time crushing an enormous deep-marine flysch basin on its inboard side. Most of interior eastern Alaska is the Yukon Tanana terrane, a very large entirely fault-bounded metamorphic-plutonic assemblage covering thousands of square kilometers in Canada as well as Alaska. The original stratigraphy and relationship to North America of the Yukon-Tanana terrane are both obscure. A collapsed Mesozoic flysch basin, similar to the one inboard of Wrangellia, lies along the northern margin. Much of Arctic Alaska was apparently a vast expanse of upper Paleozoic to Early Mesozoic deep marine sediments and mafic volcanic and plutonic rocks now scattered widely as large telescoped sheets and Klippen thrust over the Ruby geanticline and the Brooks Range, and probably underlying the Yukon-Koyukuk basin and the Yukon flats. The Brooks Range itself is a stack of north vergent nappes, the telescoping of which began in Early Cretaceous time. Despite compelling evidence for thousands of kilometers of relative displacement between the accreted terranes, and large amounts of telescoping, translation, and rotation since accretion, the resulting new continental crust added to North America in Alaska carries few obvious signatures that allow application of currently popular simple plate tectonic models. Intraplate telescoping and strike-slip translations, delamination at mid-crustal levels, and large-scale lithospheric wedging were important processes in northern Cordilleran tectonic evolution.

Alaska↗

USGS2014

Passage and survival probabilities of juvenile Chinook salmon at Cougar Dam, Oregon, 2012

This report describes studies of juvenile-salmon dam passage and apparent survival at Cougar Dam, Oregon, during two operating conditions in 2012. Cougar Dam is a 158-meter tall rock-fill dam used primarily for flood control, and passes water through a temperature control tower to either a powerhouse penstock or to a regulating outlet (RO). The temperature control tower has moveable weir gates to enable water of different elevations and temperatures to be drawn through the dam to control water temperatures downstream. A series of studies of downstream dam passage of juvenile salmonids were begun after the National Oceanic and Atmospheric Administration determined that Cougar Dam was impacting the viability of anadromous fish stocks. The primary objectives of the studies described in this report were to estimate the route-specific fish passage probabilities at the dam and to estimate the survival probabilities of fish passing through the RO. The first set of dam operating conditions, studied in November, consisted of (1) a mean reservoir elevation of 1,589 feet, (2) water entering the temperature control tower through the weir gates, (3) most water routed through the turbines during the day and through the RO during the night, and (4) mean RO gate openings of 1.2 feet during the day and 3.2 feet during the night. The second set of dam operating conditions, studied in December, consisted of (1) a mean reservoir elevation of 1,507 ft, (2) water entering the temperature control tower through the RO bypass, (3) all water passing through the RO, and (4) mean RO gate openings of 7.3 feet during the day and 7.5 feet during the night. The studies were based on juvenile Chinook salmon (Oncorhynchus tshawytscha) surgically implanted with radio transmitters and passive integrated transponder (PIT) tags. Inferences about general dam passage percentage and timing of volitional migrants were based on surface-acclimated fish released in the reservoir. Dam passage and apparent survival probabilities were estimated using the Route-Specific-Survival Model with data from surface-acclimated fish released near the water surface directly upstream of the temperature control tower (treatment group) and slightly downstream of the dam (control group). In this study, apparent survival is the joint probability of surviving and migrating through the study area during the life of the transmitters. Two rearing groups were used to enable sufficient sample sizes for the studies. The groups differed in feed type, and for the December study only, the rearing location. Fish from each group were divided nearly equally among all combinations of release sites, release times, and surgeons. The sizes, travel times, and survivals of the two rearing groups were similar. There were statistical differences in fish lengths and travel times of the two groups, but they were small and likely were not biologically meaningful. There also was evidence of a difference in single-release estimates of survival between the rearing groups during the December study, but the differences had little effect on the relative survival estimates so the analyses of passage and survival were based on data from the rearing groups pooled. Conditions during the December study were more conducive to passing volitionally migrating fish than conditions during the November study. The passage percentage of the fish released in the reservoir was similar between studies (about 70 percent), but the passage occurred in a median of 1.0 day during the December study and a median of 9.3 days during the November study. More than 93 percent of the dam passage of volitionally migrating fish occurred at night during each study. This finding corroborates results of previous studies at Cougar Dam and suggests that the operating conditions at night are most important to volitionally migrating fish, given the current configuration of the dam. Most fish released near the temperature control tower passed through the RO. A total of 92.2 percent of the treatment group passed through the RO during the November study and the RO was the only route open during the December study. The assumptions of the survival model were either met or adjusted for during each study. There was little evidence that tagger skill or premature failure of radio transmitters had an effect on survival estimates. There were statistically significant differences in travel times between treatment and control groups through several of the river reaches they had in common, but the differences were typically only a few hours, and the two groups likely experienced the same in-river conditions. There was direct evidence of bias due to detection of euthanized fish with live transmitters released as part of the study design. The bias was ameliorated by adjusting the survival estimates for the probability of detecting dead fish with live transmitters, which reduced the estimated survival probabilities by about 0.02. The data and models indicated that the treatment effect was not fully expressed until the study reach terminating with Marshall Island Park on the Willamette River, a distance of 105.8 kilometers downstream of Cougar Dam. This was the first reach in which the 95-percent confidence interval of the estimated reach-specific relative survival overlapped 1.0, indicating similar survival of treatment and control groups. The median travel time of the treatment group from release to Marshall Island Park was 1.64 days during the November study and 1.36 days during the December study. The survival probability of fish that passed into the RO was greater during the December study than during the November study. The relative survival probability of fish passing through the RO was 0.4594 (standard error [SE] 0.0543) during the November study and 0.7389 (SE 0.1160) during the December study. These estimates represent relative survival probabilities from release near Cougar Dam to the Marshall Island site. The estimated survival probability of RO passage was lower than previous studies based on balloon and PIT tags, but higher than a similar study based on radio transmitters. We suggest that, apart from dam operations, the differences in survival primarily are due to the release location. We hypothesize that the balloon- and PIT-tagged fish released through a hose at a point near the RO gate opening experienced more benign conditions than the radio-tagged fish passing the RO volitionally. This hypothesis could be tested with further study. An alternative hypothesis is that some live fish remained within the study area beyond the life of their radio transmitter. The results from these and previous studies indicate that entrainment and survival of juvenile salmonids passing Cougar Dam varies with dam operating conditions. The condition most conducive to dam passage has been the discharge and low pool elevation condition tested during December 2012. That condition included large RO gate openings and was the condition with the highest dam passage survival.

Oregon↗

USGS2023

Camera trap distance sampling survey design, Andersen Airforce Base, Guam

Reliable population estimates of animal density is one of the most elementary needs for the control and management of wildlife, particularly for introduced ungulates on oceanic islands. On Guam, Philippine deer (Rusa marianna) and wild pigs (Sus scrofa; wild boar and descendants of domestic pigs) cause agricultural and ecological damage and are hunted for recreational, nutritional, and cultural uses. Most common population estimation methods are based on capture-recapture and related methods that require marking or uniquely identifying individuals. Capturing, marking, and either recapturing or resighting individuals repeatedly is labor intensive and expensive. In many situations marking or individually distinguishing animals is not feasible, necessitating estimating densities and abundance from unmarked animal populations. Motion-triggered camera traps are a relatively low-cost approach that can be used to generate presence/pseudo-absence and indices of relative abundance on multiple species simultaneously. We used distance sampling with camera traps to estimate deer and pig densities from non-independent observations of unmarked animals while accounting for imperfect detection where some present individuals are not detected. We present methods to (1) process the digital imagery data automatically for species detection and species categorization using a machine learning algorithm, (2) automatically estimate distance to detected species using a separate machine learning algorithm, and (3) estimate densities using distance sampling with camera trap methods. We compare accuracy statistics and results of ungulate densities estimated from automated methods to those estimated from manual assessment. We collected 7,695 videos: 381 videos contained deer and 377 contained pigs. The object detection and identification model performed well with overall accuracy above 80% and F1 scores above 0.9. The hazard-rate key detection function was chosen for deer and pigs based on Akaike’s information criterion accounting for overdispersion. Deer density estimates were 0.53 ± 0.20 deer/ha with higher density in the Plateau area than the Tarague area of Guam. Pig density estimates were 0.53 ± 0.32 pigs/ha, also with higher densities in the Plateau area than the Tarague area. Coefficients of variation ranged from 0.38 to 1.15, and greater numbers of camera traps would be required for pigs than deer to achieve desired coefficients of variation. On average, 101.9 ± 82.3 deer and 131.6 ± 118.8 pigs were detected per day. Microsite heterogeneity affected densities where orientation-specific estimates were less precise than estimates made with the full dataset. We developed a camera trap survey design based on standard camera trapping sampling protocols using motion-activated, digital cameras and determined that distance sampling methods using camera traps produce reliable densities of unmarked deer and pigs on Guam. Our camera trap survey design is based on a regularly sized trapping grid that is generalizable and can be expanded to survey other areas of Guam.

Guam↗

USGS1996

Comparison of the partitioning behaviours of yttrium, rare earth elements, and titanium between hydrogenetic marine ferromanganese crusts and seawater

In order to evaluate details of the partitioning behaviours of Y, rare earth elements (REEs), and Ti between inorganic metal oxide surfaces and seawater, we studied the distribution of these elements in hydrogenetic marine ferromanganese (Fe-Mn) crusts from the Central Pacific Ocean. Nonphosphatized Fe-Mn crusts display shale-normalized rare earths and yttrium (REY SN ) patterns (Y inserted between Dy and Ho) that are depleted in light REEs (LREEs) and which show negative anomalies for Y sn , and positive anomalies for La SN , Eu SN , Gd SN , and in most cases, Ce sn . They show considerably smaller Y/ Ho ratios than seawater or common igneous and clastic rocks, indicating that Y and Ho are fractionated in the marine environment. Compared to P-poor crusts, REY SN patterns of phosphatized Fe-Mn crusts are similar, but yield pronounced positive Y sn anomalies, stronger positive La SN anomalies, and enrichment of the HREEs relative to the MREEs. The data suggest modification of REY during phosphatization and indicate that studies requiring primary REY distributions or isotopic ratios should be restricted to non-phosphatized (layers of) Fe-Mn crusts. Apparent bulk coefficients, K d m , describing trace metal partitioning between nonphosphatized hydrogenetic Fe-Mn crusts and seawater, are similar for Pr to Eu and decrease for Eu to Yb. Exceptionally high values of K D Ce , which are similar to those of Ti, result from oxidative scavenging of Ce and support previous suggestions that Ce (IV) is a hydroxide-dominated element in seawater. Yttrium and Gd show lower K D values than their respective neighbours in the REY series. Results of modelling the exchange equilibrium between REY dissolved in seawater and REY sorbed on hydrous Fe-Mn oxides corroborate previous studies that suggested the surface complexation of REY can be approximated by their first hydroxide binding constant. Negative “anomalies” occur for stabilities of bulk surface complexes of Gd, La, and particularly Y. The differences in inorganic surface complex stability between Y and Ho and between Gd and its REE neighbours are similar to those shown by the stabilities of complexes with aminocarboxylic acids and are significantly larger than those shown by stabilities of complexes with carboxylic acids. Hence, sorption of Y and REEs onto hydrous Fe-Mn oxides may contribute significantly to the positive Y SN and Gd SN anomalies in seawater.

Geochimica et Cosmochimica Acta↗