USGS Science⌕ Search

SEARCH · USGS Science

Results for “Ocean Modelling”

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

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

Stabilization of volcanic flanks by dike intrusion: An example from Kilauea

Dike propagation and dilation increases the compression of adjacent rocks. On volcanoes, especially oceanic shields, dikes are accordingly thought to be structurally destabilizing. As compression is incremented, volcanic flanks are driven outward or downslope and thus increase their susceptibility to destructive earthquakes and giant landslides. We show, however, that the 2-m-thick dike emplaced along the east rift zone of Kilauea in 1983 actually stabilized that volcano's flank. Specifically, production of flank earthquakes dropped more than twofold after 1983 as maximum downslope motion slowed to 6 cm ?? year-1 from approximately 40 cm ?? year-1 during 1980-1982. As much as 65 cm of deflationary subsidence above Kilauea's summit and upper rift zones accompanied the dike intrusion. According to recent estimates, this deflation corresponds to a reduction in magma-reservoir pressure of approximately 4 MPa, probably about as much as the driving pressure of the 1983 dike. The volume of the dike, approximately 0.10-0.15 km3, is orders of magnitude less than the estimated 200- to 250-km3 volume of Kilauea's reservoir of magma and nearby hot, mushy rock. Thus, deflation of that reservoir reduces the compressional load on the flank over a much larger area than intrusion of the dike adds to it, particularly at the dominant depth of seismicity, 8-9 km. A Coulomb block model for flank motion during intervals between major earthquakes requires the low-angle fault beneath Kilauea's flank to exhibit slip weakening, conducive to earthquake instability. Accordingly, the triggering mechanism of destructive earthquakes, several of which have struck Hawaii during the past 150 years, need not require stresses accumulated by dike intrusions.

Bulletin of Volcanology↗

Metal induced chitinozoan malformations link the Mid-Ludfordian OAE to the venting of metalliferous hydrothermal basin brines into the Silurian Ocean

The mechanisms driving Silurian Ocean Anoxic Events (OAEs) and their impact on marine biota remain poorly understood. Previous studies linked chitinozoan (fossil zooplankton) malformations to heavy metal toxicity under expanding reducing oceanic conditions during the onset of some of these events. Importantly, these malformations represent in vivo reactions to changing paleoenvironmental stress and may serve as a sensitive indicator of marine metal loading and the geochemical processes driving Silurian extinction events. In the current study this hypothesis is tested for the prominent Mid- Ludfordian (upper Silurian) OAE. The onset of the event is examined using the thick and well preserved stratigraphic records in outcrop at Bodudd 1 (Gotland, Sweden). A bed-by-bed collection of 101 samples were analyzed for their whole rock geochemistry with chitinozoan malformations quantified in 66 samples. A total of 587.241 specimens were counted with 133 confirmed malformations and 87 potentially malformed, occurring in two intervals that correlate with Mn enrichments. Notably, malformations appear in the lowermost sample in the section, below the carbon isotope excursion, indicating that biotic stress preceded perturbations to the global carbon cycle. In addition to Mn, anomalies in Zn and Pb are present, alongside minor enrichments in Hg, Co, and Sb. The stratigraphic position of metal anomalies corresponds with an extraordinarily rapid rise in marine Sr isotope compositions. The pulsed nature of trace metal enrichments, the distinct elemental suite, and diagnostic spikes in radiogenic 87 Sr/ 86 Sr support the emerging model that episodic venting of metalliferous hydrothermal brines was the primary trigger for this and other Silurian OAEs.

Gotland↗

Seasat synthetic aperture radar ( SAR) response to lowland vegetation types in eastern Maryland and Virginia

Examination of SEASAT SAR images of eastern Maryland and Virginia reveals botanical distinctions between vegetated lowland areas and adjacent upland areas. Radar returns from the lowland areas can be either brighter or darker than returns from the upland forests. Scattering models and scatterometer measurements predict an increase of 6 dB in backscatter from vegetation over standing water. This agrees with the SCWigital number (DN) increase observed in the digital SEASAT data. The brightest areas in the Chickahominy, Virginia, drainage, containing P. virginica about 0.4 m high, contrast with the brightest areas in the Blackwater, Maryland, marshes, which contain mature loblolly pine in standing water. The darkest vegetated area in the Chickahominy drainage contains a forest of Nyssa aquatica (water tupelo) about 18 m high, while the darkest vegetated area in the Blackwater marshes contains the marsh plant Spartina alterniflora, 0.3 m high. The density, morphology, and relative geometry of the lowland vegetation with respect to standing water can all affect the strength of the return L band signal.

Journal of Geophysical Research - Oceans↗

Application of Wind Fetch and Wave Models for Habitat Rehabilitation and Enhancement Projects

Models based upon coastal engineering equations have been developed to quantify wind fetch length and several physical wave characteristics including significant height, length, peak period, maximum orbital velocity, and shear stress. These models, developed using Environmental Systems Research Institute's ArcGIS 9.2 Geographic Information System platform, were used to quantify differences in proposed island construction designs for three Habitat Rehabilitation and Enhancement Projects (HREPs) in the U.S. Army Corps of Engineers St. Paul District (Capoli Slough and Harpers Slough) and St. Louis District (Swan Lake). Weighted wind fetch was calculated using land cover data supplied by the Long Term Resource Monitoring Program (LTRMP) for each island design scenario for all three HREPs. Figures and graphs were created to depict the results of this analysis. The difference in weighted wind fetch from existing conditions to each potential future island design was calculated for Capoli and Harpers Slough HREPs. A simplistic method for calculating sediment suspension probability was also applied to the HREPs in the St. Paul District. This analysis involved determining the percentage of days that maximum orbital wave velocity calculated over the growing seasons of 2002-2007 exceeded a threshold value taken from the literature where fine unconsolidated sediments may become suspended. This analysis also evaluated the difference in sediment suspension probability from existing conditions to the potential island designs. Bathymetric data used in the analysis were collected from the LTRMP and wind direction and magnitude data were collected from the National Oceanic and Atmospheric Administration, National Climatic Data Center.

Open-File Report↗

Oregon OCS seafloor mapping: Selected lease blocks relevant to renewable energy

In 2014 the U.S. Geological Survey (USGS) and the Bureau of Ocean Energy Management (BOEM) entered into Intra-agency agreement M13PG00037 to map an area of the Oregon Outer Continental Shelf (OCS) off of Coos Bay, Oregon, under consideration for development of a floating wind energy farm. The BOEM requires seafloor mapping and site characterization studies in order to evaluate the impact of seafloor and sub-seafloor conditions on the installation, operation, and structural integrity of proposed renewable energy projects, as well as to assess the potential effects of construction and operations on archaeological resources. The mission of the USGS is to provide geologic, topographic, and hydrologic information that contributes to the wise management of the Nation's natural resources and that promotes the health, safety, and well being of the people. This information consists of maps, databases, and descriptions and analyses of the water, energy, and mineral resources, land surface, underlying geologic structure, and dynamic processes of the earth. For the Oregon OCS study, the USGS acquired multibeam echo sounder and seafloor video data surrounding the proposed development site, which is 95 km2 in area and 15 miles offshore from Coos Bay. The development site had been surveyed by Solmar Hydro Inc. in 2013 under a contract with WindFloat Pacific. The USGS subsequently produced a bathymetry digital elevation model and a backscatter intensity grid that were merged with existing data collected by the contractor. The merged grids were published along with visual observations of benthic geo-habitat from the video data in an associated USGS data release (Cochrane and others, 2015). This report includes the results of analysis of the video data conducted by Oregon State University and the geo-habitat interpretation of the multibeam echo sounder (MBES) data conducted by the USGS. MBES data was published in Cochrane and others (2015). Interpretive data associated with this publication is published in Cochrane (2017). All the data is provided as geographic information system (GIS) files that contain both Esri ArcGIS geotiffs or shapefiles. For those who do not own the full suite of Esri GIS and mapping software, the data can be read using Esri ArcReader, a free viewer that is available at http://www.esri.com/software/arcgis/arcreader/index.html (last accessed August 29, 2016). Web services, which consist of standard implementations of ArcGIS representational state transfer (REST) Service and Open Geospatial Consortium (OGC) GIS web map service (WMS), also are available for all published GIS data. Web services were created using an ArcGIS service definition file, resulting in data layers that are symbolized as shown on the associated report figures. Both the ArcGIS REST Service and OGC WMS Service include all the individual GIS layers. Data layers are bundled together in a map-area web service; however, each layer can be symbolized and accessed individually after the web service is ingested into a desktop application or web map. Web services enable users to download and view data, as well as to easily add data to their own workflows, using any browser-enabled, standalone or mobile device. Though the surficial substrate is dominated by combinations of mud and sand substrate, a diverse assortment of geomorphologic features are related to geologic processes—one anticlinal ridge where bedrock is exposed, a slump and associated scarps, and pockmarks. Pockmarks are seen in the form of fields of small pockmarks, a lineation of large pockmarks with methanogenic carbonates, and areas of large pockmarks that have merged into larger variously shaped depressions. The slump appears to have originated at the pockmark lineation. Video-supervised numerical analysis of the MBES backscatter intensity data and vector ruggedness derived from the MBES bathymetry data was used to produce a substrate model called a seafloor character raster for the study area. The seafloor character raster consists of three substrate classes: soft-flat areas, hard-flat areas, and hard-rugged areas. A Coastal and Marine Ecological Classification Standard (CMECS) geoform and substrate map was also produced using depth, slope, and benthic position index classes to delineate geoform boundaries. Seven geoforms were identified in this process, including ridges, slump scars, slump deposits, basins, and pockmarks. Statistical analysis of the video data for correlations between substrate, depth, and invertebrate assemblages resulted in the identification of seven biomes: three hard-bottom biomes and four softbottom biomes. A similar analysis of vertebrate observations produces a similar set of biomes. The biome between-group dissimilarity was very high or high. Invertebrates alone represent most of the structure of the whole benthic community into different assemblages. A biotope map was generated using the seafloor character raster and the substrate and depth values of the biomes. Hard substrate biotopes were small in size and were located primarily on the ridge and in pockmarks along the pockmark lineation. The soft-bottom bitopes consisted of large contiguous areas delimited by isobaths.

Open-File Report↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗