SEARCH · USGS Science

Results for “One Earth”

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,423 records · Page 79Linked to original sources

USGS2021

The making of the NEAM Tsunami Hazard Model 2018 (NEAMTHM18)

The NEAM Tsunami Hazard Model 2018 (NEAMTHM18) is a probabilistic hazard model for tsunamis generated by earthquakes. It covers the coastlines of the North-eastern Atlantic, the Mediterranean, and connected seas (NEAM). NEAMTHM18 was designed as a three-phase project. The first two phases were dedicated to the model development and hazard calculations, following a formalized decision-making process based on a multiple-expert protocol. The third phase was dedicated to documentation and dissemination. The hazard assessment workflow was structured in Steps and Levels. There are four Steps: Step-1) probabilistic earthquake model; Step-2) tsunami generation and modeling in deep water; Step-3) shoaling and inundation; Step-4) hazard aggregation and uncertainty quantification. Each Step includes a different number of Levels. Level-0 always describes the input data; the other Levels describe the intermediate results needed to proceed from one Step to another. Alternative datasets and models were considered in the implementation. The epistemic hazard uncertainty was quantified through an ensemble modeling technique accounting for alternative models’ weights and yielding a distribution of hazard curves represented by the mean and various percentiles. Hazard curves were calculated at 2,343 Points of Interest (POI) distributed at an average spacing of ∼20 km. Precalculated probability maps for five maximum inundation heights (MIH) and hazard intensity maps for five average return periods (ARP) were produced from hazard curves. In the entire NEAM Region, MIHs of several meters are rare but not impossible. Considering a 2% probability of exceedance in 50 years (ARP≈2,475 years), the POIs with MIH >5 m are fewer than 1% and are all in the Mediterranean on Libya, Egypt, Cyprus, and Greece coasts. In the North-East Atlantic, POIs with MIH >3 m are on the coasts of Mauritania and Gulf of Cadiz. Overall, 30% of the POIs have MIH >1 m. NEAMTHM18 results and documentation are available through the TSUMAPS-NEAM project website (http://www.tsumaps-neam.eu/), featuring an interactive web mapper. Although the NEAMTHM18 cannot substitute in-depth analyses at local scales, it represents the first action to start local and more detailed hazard and risk assessments and contributes to designing evacuation maps for tsunami early warning.

Frontiers in Earth Science↗

USGS1991

Numerical simulations of hydrothermal circulation resulting from basalt intrusions in a buried spreading center

A two-dimensional, one by two-kilometer section through the seafloor was simulated with a numerical model to investigate coupled fluid and heat flow resulting from basalt intrusions in a buried spreading center. Boundary and initial conditions and physical properties of both sediments and basalt were constrained by field surveys and drilling in the Guaymas Basin, central Gulf of California. Parametric variations in these studies included sediment and basalt permeability, anisotropy in sediment permeability, and the size of heat sources. Faults were introduced through new intrusions both before and after cooling. Background heat input caused fluid convection at velocities ≤ 3 cm a −1 through shallow sediments. Eighty to ninety percent of the heat introduced at the base of the simulations exited through the upper, horizontal surface, even when the vertical boundaries were made permeable to fluid flow. The simulated injection of a 25–50 m thick basalt intrusion at a depth of 250 m resulted in about 10 yr of pore-fluid expulsion through the sea-floor in all cases, leaving the sediments above the intrusions strongly underpressured. A longer period of fluid recharge followed, sometimes accompanied by reductions in total seafloor heat output of 10% in comparison to pre-intrusion values. Additional discharge-recharge events were dispersed chaotically through the duration of the cooling period. These cycles in heat and fluid flow resulted from the response of the simulated system to a thermodynamic shock, the sudden emplacement of a large heat source, and not from mechanical displacement of sediments and pore fluids, which was not simulated. Water/rock mass ratios calculated from numerical simulations are in good agreement with geochemical estimates from materials recovered from the Guaymas Basin, assuming a bulk basalt permeability value of at least 10 −17 m 2 /(10 −2 mD). The addition of faults through intrusions and sediments in these simulations did not facilitate continuous, rapid venting. Increased heat input at the base of the faults resulted in temporarily greater fluid discharge, but the flow could not be sustained because the modeled system could not recharge cold fluid quickly enough to remove sufficient heat through the vents.

Earth and Planetary Science Letters↗

USGS2012

Heat flow and energetics of the San Andreas fault zone

Approximately 100 heat flow measurements in the San Andreas fault zone indicate (1) there is no evidence for local factional heating of the main fault trace at any latitude over a 1000-km length from Cape Mendocino to San Bernardino, (2) average heat flow is high (∼2 HFU, ∼80 mW m −2 ) throughout the 550-km segment of the Coast Ranges that encloses the San Andreas fault zone in central California; this broad anomaly falls off rapidly toward the Great Valley to the east, and over a 200-km distance toward the Mendocino Triple Junction to the northwest. As others have pointed out, a local conductive heat flow anomaly would be detectable unless the frictional resistance allocated to heat production on the main trace were ≲100 bars. Frictional work allocated to surface energy of new fractures is probably unimportant, and hydrologic convection is not likely to invalidate the conduction assumption, since the heat discharge by thermal springs near the fault is negligible. Explanations for the low dynamic friction fall into two intergradational classes: those in which the fault is weak all of the time and those in which it is weak only during earthquakes (possibly just large ones). The first class includes faults containing anomalously weak gouge materials and faults containing materials with normal frictional properties under near-lithostatic steady state fluid pressures. In the second class, weakening is caused by the event (for example, a thermally induced increase in fluid pressure, dehydration of clay minerals, or acoustic fluidization). In this class, unlike the first, the average strength and ambient tectonic shear stress may be large, ∼1 kbar, but the stress allocated to elastic radiation (the apparent stress) must be of similar magnitude, an apparent contradiction with seismic estimates. Unless seismic radiation is underestimated for large earthquakes, it is difficult to justify average tectonic stresses on the main trace of the San Andreas fault in excess of ∼200 bars. The development of the broad Coast Range heat flow anomaly southward from Cape Mendocino suggests that heat flow increases by a factor of 2 within 4 m.y. after the passage of the Mendocino Triple Junction. This passage leaves the San Andreas transform fault zone in its wake; the depth of the anomalous sources cannot be much greater than the depth of the seismogenic layer. Some of the anomalous heat may be supplied by conduction from the warmer mantle that must occur south of the Mendocino transform (where there is no subducting slab), and some might be supplied by shear heating in the fault zone. With no contribution from shear heating, extreme mantle upwelling would be required, and asthenosphere conditions should exist today at depths of only ∼20 km in the northernmost Coast Ranges. If there is an appreciable contribution from shear heating, the heat flow constraint implies that the seismogenic layer is partially decoupled at its base and that the basal traction is in the sense that resists right lateral motion on the fault(s). As a result of these basal tractions, the average shearing stress in the seismogenic layer would increase with distance from the main fault, and the seismogenic layer would offer substantial resistance to plate motion even though resistance on the main fault might be negligible. These speculative models have testable consequences.

Journal of Geophysical Research Solid Earth↗

USGS2020

Changing suspended sediment in United States rivers and streams: Linking sediment trends to changes in land use/cover, hydrology and climate

Sediment is one of the leading pollutants in rivers and streams across the United States (US) and the world. Between 1992 and 2012, concentrations of annual mean suspended sediment decreased at over half of the 137 stream sites assessed across the contiguous US. Increases occurred at less than 25 % of the sites, and the direction of change was uncertain at the remaining 25 %. Sediment trends were characterized using the Weighted Regressions on Time, Discharge, and Season (WRTDS) model, and decreases in sediment ranged from −95 % to −8.5 % of the 1992 concentration. To explore potential drivers of these changes, the sediment trends were (1) parsed into two broad contributors of change, changes in land management versus changes in the streamflow regime, and (2) grouped by land use of the watershed and correlated to concurrent changes in land use or land cover (land use/cover), hydrology and climate variables and static/long-term watershed characteristics. At 83 % of the sites, changes in land management (captured by changes in the concentration–streamflow relationship over time; C–Q relationship) contributed more to the change in the sediment trend than changes in the streamflow regime alone (i.e., any systematic change in the magnitude, frequency or timing of flows). However, at >50 % of the sites, changes in the streamflow regime contributed at least a 5 % change in sediment, and at 11 sites changes in the streamflow regime contributed over half the change in sediment, indicating that at many sites changes in streamflow were not the main driver of changes in sediment but were often an important supporting factor. Correlations between sediment trends and concurrent changes in land use/cover, hydrology and climate were often stronger at sites draining watersheds with more homogenous, human-related land uses (i.e., agricultural and urban lands) compared to mixed-use or undeveloped lands. At many sites, decreases in sediment occurred despite small-to-moderate increases in the amount of urban or agricultural land in the watershed, suggesting conservation efforts and best-management practices (BMPs) used to reduce sediment runoff to streams may be successful, up to a point, as lands are converted to urban and agricultural uses.

Hydrology and Earth System Sciences↗

USGS2009

Optical satellite data volcano monitoring: a multi-sensor rapid response system

In this chapter, the use of satellite remote sensing to monitor active geological processes is described. Specifically, threats posed by volcanic eruptions are briefly outlined, and essential monitoring requirements are discussed. As an application example, a collaborative, multi-agency operational volcano monitoring system in the north Pacific is highlighted with a focus on the 2007 eruption of Kliuchevskoi volcano, Russia. The data from this system have been used since 2004 to detect the onset of volcanic activity, support the emergency response to large eruptions, and assess the volcanic products produced following the eruption. The overall utility of such integrative assessments is also summarized. The work described in this chapter was originally funded through two National Aeronautics and Space Administration (NASA) Earth System Science research grants that focused on the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) instrument. A skilled team of volcanologists, geologists, satellite tasking experts, satellite ground system experts, system engineers and software developers collaborated to accomplish the objectives. The first project, Automation of the ASTER Emergency Data Acquisition Protocol for Scientific Analysis, Disaster Monitoring, and Preparedness , established the original collaborative research and monitoring program between the University of Pittsburgh (UP), the Alaska Volcano Observatory (AVO), the NASA Land Processes Distributed Active Archive Center (LP DAAC) at the U.S. Geological Survey (USGS) Earth Resources Observation and Science (EROS) Center, and affiliates on the ASTER Science Team at the Jet Propulsion Laboratory (JPL) as well as associates at the Earth Remote Sensing Data Analysis Center (ERSDAC) in Japan. This grant, completed in 2008, also allowed for detailed volcanic analyses and data validation during three separate summer field campaigns to Kamchatka Russia. The second project, Expansion and synergistic use of the ASTER Urgent Request Protocol (URP) for natural disaster monitoring and scientific analysis , has expanded the project to other volcanoes around the world and is in progress through 2011. The focus on ASTER data is due to the suitability of the sensor for natural disaster monitoring and the availability of data. The instrument has several unique facets that make it especially attractive for volcanic observations (Ramsey and Dehn, 2004). Specifically, ASTER routinely collects data at night, it has the ability to generate digital elevation models using stereo imaging, it can collect data in various gain states to minimize data saturation, it has a cross-track pointing capability for faster targeting, and it collects data up to ±85° latitude for better global coverage. As with any optical imaging-based remote sensing, the viewing conditions can negatively impact the data quality. This impact varies across the optical and thermal infrared wavelengths as well as being a function of the specific atmospheric window within a given wavelength region. Water vapor and cloud formation can obscure surface data in the visible and near infrared (VNIR)/shortwave infrared (SWIR) region due mainly to non-selective scattering of the incident photons. In the longer wavelengths of the thermal infrared (TIR), scattering is less of an issue, but heavy cloud cover can still obscure the ground due to atmospheric absorption. Thin clouds can be optically-transparent in the VNIR and TIR regions, but can cause errors in the extracted surface reflectance or derived surface temperatures. In regions prone to heavy cloud cover, optical remote sensing can be improved through increased temporal resolution. As more images are acquired in a given time period the chances of a clear image improve dramatically. The Advanced Very High Resolution Radiometer (AVHRR) routine monitoring, which commonly collects 4-6 images per day of any north Pacific volcano, takes advantage of this fact. The rapid response program described in this chapter also improves the temporal resolution of the ASTER instrument. ASTER has been acquiring images of volcanic eruptions since soon after its launch in December 1999. An early example included the observations of the large pyroclastic flow deposit emplaced at Bezymianny volcano in Kamchatka, Russia. The first images in March 2000, just weeks after the eruption, revealed the extent, composition, and cooling history of this large deposit and of the active lava dome (Ramsey and Dehn, 2004). The initial results from these early datasets spurred interest in using ASTER data for expanded volcano monitoring in the north Pacific. It also gave rise to the multi-year NASA-funded programs of rapid response scheduling and imaging throughout the Aleutian, Kamchatka and Kurile arcs. Since the formal establishment of the programs, the data have provided detailed descriptions of the eruptions of Augustine, Bezymianny, Kliuchevskoi and Sheveluch volcanoes over the past nine years (Wessels et al., in press; Carter et al., 2007, 2008; Ramsey et al., 2008; Rose and Ramsey, 2009). The initial research focus of this rapid response program was specifically on automating the ASTER sensor’s ability for targeted observational scheduling using the expedited data system. This urgent request protocol is one of the unique characteristics of ASTER. It provides a limited number of emergency observations, typically at a much-improved temporal resolution and quicker turnaround with data processing in the United States rather than in Japan. This can speed the reception of the processed data by several days to a week. The ongoing multi-agency research and operational collaboration has been highly successful. AVO serves as the primary source for status information on volcanic activity, working closely with the National Weather Service (NWS), Federal Aviation Administration (FAA), military and other state and federal emergency services. Collaboration with the Russian Institute of Volcanology and Seismology (IVS)/Kamchatka Volcanic Eruption Response Team (KVERT) is also maintained. Once a volcano is identified as having increased thermal output, ASTER is automatically tasked and the volcano is targeted at the next available opportunity. After the data are acquired, scientists at all the agencies have access to the images, with the primary science analysis carried out at the University of Pittsburgh and AVO. Results are disseminated to the responsible monitoring agencies and the global community through e-mail mailing lists.

Book chapter↗

USGS2012

Wandering terranes in southern Alaska: The Aleutia Microplate and implications for the Bering Sea

Paleomagnetic and geological data suggest that much of southern Alaska is a collage of tectonostratigraphic terranes which originated in Mesozoic time at paleolatitudes far south of their present position. The time of ‘docking’ of the terranes against cratonic Alaska is critical to defining their amalgamated size and extent during their northward motion as well as their role in the evolution of the Bering Sea. One of the largest of the tectonostratigraphic terranes, the Peninsular terrane of south central and southwestern Alaska, extends offshore along the outer Bering Sea continental margin (Beringia). Paleomagnetic data suggest that this terrane has moved northward through all of Cenozoic time, but geologic data imply that the terrane had accreted to Alaska by the end of the Mesozoic. In early Cenozoic time the eastern part of the Aleutian arc appears to have been superimposed on the Peninsular terrane, and postulated northward Cenozoic motion of the terrane would therefore have required northward motion of the arc. Two accretion models, based on docking times for terranes in Alaska, are proposed, and they illustrate that large areas of the abyssal Bering Sea, the Alaska Peninsula, the Aleutian arc, and the Beringian continental margin may be part of a superterrane or microplate called Aleutia (microplate as defined by Beck et al. (1980), i.e., a microplate is a displaced segment of lithosphere that has crustal roots, whereas a superterrane is an amalgamation of terranes which may or may not be rootless). Model A implies that the Aleutian arc developed in situ on the southern edge of Aleutia after the microplate had docked. In model B, the final docking time of the Peninsular terrane is late Cenozoic, which implies that the Aleutia microplate encompasses a mammoth area that includes parts of southern Alaska, the Alaska Peninsula, the southern Beringian margin, the abyssal Bering Sea (Kula plate), and the Aleutian arc. If model A is correct, the docking time of the Peninsular terrane is late Mesozoic or earliest Tertiary. The Aleutia microplate in this model is made up solely of the abyssal Bering Sea (Kula plate), which presumably docked at the same time or slightly after the Peninsular terrane accreted against Alaska. If model B is correct, that is, if the Aleutia collided with nuclear Alaska during the Cenozoic, then a late Cenozoic suture zone, the vestige of a large open sea that must have closed between Aleutia and Alaska, must exist in south central and southwest Alaska. Either evidence for Cenozoic closure and suturing has been obliterated in Alaska or the inferences of Cenozoic terrane motion derived from paleomagnetic data are suspect.

Alaska↗

USGS2022

Barkley Canyon gas hydrates: A synthesis based on two decades of seafloor observation and remote sensing

Barkley Canyon is one of the few known sites worldwide with the occurrence of thermogenic gas seepage and formation of structure-II and structure-H gas hydrate mounds on the seafloor. This site is the location of continuous seafloor monitoring as part of the Ocean Networks Canada (ONC) cabled observatory off the west coast off Vancouver Island, British Columbia, Canada. We combine repeat remotely operated vehicle (ROV) seafloor video observations, mapping with an autonomous underwater vehicle (AUV), ship-, ROV-, and AUV-based identification of gas flares, as well as seismic and Chirp data to investigate the distribution of fluid migration pathways. Geologically, the site with the prominent gas hydrate mounds and associated fluid seepage is covering an area of ∼0.15 km 2 and is situated on a remnant of a rotated fault block that had slipped off the steep flanks of the north-east facing canyon wall. The gas hydrate mounds, nearly constant in dimension over the entire observation period, are associated with gas and oil seepage and surrounded by debris of chemosynthetic communities and authigenic carbonate. The formation of gas hydrate at and near the seafloor requires additional accommodation space created by forming blisters at the seafloor that displace the regular sediments. An additional zone located centrally on the rotated fault block with more diffuse seepage (∼0.02 km 2 in extent) has been identified with no visible mounds, but with bacterial mats, small carbonate concretions, and clam beds. Gas venting is seen acoustically in the water column up to a depth of ∼300 m. However, acoustic water-column imaging during coring and ROV dives showed rising gas bubbles to much shallower depth, even <50 m, likely a result of degassing of rising oil droplets, which themselves cannot be seen acoustically. Combining all observations, the location of the gas hydrate mounds is controlled by a combination of fault-focused fluid migration from a deeper reservoir and fluid seepage along more permeable strata within the rotated slope block. Fluids must be provided continuously to allow the sustained presence of the gas hydrate mounds at the seafloor.

Barkley Canyon↗

USGS2012

A multi-source satellite data approach for modelling Lake Turkana water level: Calibration and validation using satellite altimetry data

Lake Turkana is one of the largest desert lakes in the world and is characterized by high degrees of inter- and intra-annual fluctuations. The hydrology and water balance of this lake have not been well understood due to its remote location and unavailability of reliable ground truth datasets. Managing surface water resources is a great challenge in areas where in-situ data are either limited or unavailable. In this study, multi-source satellite-driven data such as satellite-based rainfall estimates, modelled runoff, evapotranspiration, and a digital elevation dataset were used to model Lake Turkana water levels from 1998 to 2009. Due to the unavailability of reliable lake level data, an approach is presented to calibrate and validate the water balance model of Lake Turkana using a composite lake level product of TOPEX/Poseidon, Jason-1, and ENVISAT satellite altimetry data. Model validation results showed that the satellite-driven water balance model can satisfactorily capture the patterns and seasonal variations of the Lake Turkana water level fluctuations with a Pearson's correlation coefficient of 0.90 and a Nash-Sutcliffe Coefficient of Efficiency (NSCE) of 0.80 during the validation period (2004–2009). Model error estimates were within 10% of the natural variability of the lake. Our analysis indicated that fluctuations in Lake Turkana water levels are mainly driven by lake inflows and over-the-lake evaporation. Over-the-lake rainfall contributes only up to 30% of lake evaporative demand. During the modelling time period, Lake Turkana showed seasonal variations of 1–2 m. The lake level fluctuated in the range up to 4 m between the years 1998–2009. This study demonstrated the usefulness of satellite altimetry data to calibrate and validate the satellite-driven hydrological model for Lake Turkana without using any in-situ data. Furthermore, for Lake Turkana, we identified and outlined opportunities and challenges of using a calibrated satellite-driven water balance model for (i) quantitative assessment of the impact of basin developmental activities on lake levels and for (ii) forecasting lake level changes and their impact on fisheries. From this study, we suggest that globally available satellite altimetry data provide a unique opportunity for calibration and validation of hydrologic models in ungauged basins.

Lake Turkana↗

USGS2012

The Steens Mountain (Oregon) geomagnetic polarity transition, 2. Field intensity variations and discussion of reversal models

We carried out an extensive paleointensity study of the 15.5±0.3 m.y. Miocene reversed‐to‐normal polarity transition recorded in lava flows from Steens Mountain (south central Oregon). One hundred eighty‐five samples from the collection whose paleodirectional study is reported by Mankinen et al. (this issue) were chosen for paleointensity investigations because of their low viscosity index, high Curie point and reversibility, or near reversibility, of the strong field magnetization curve versus temperature. Application of the Thellier stepwise double heating method was very successful, yielding 157 usable paleointensity estimates corresponding to 73 distinct lava flows. After grouping successive lava flows that did not differ significantly in direction and intensity, we obtained 51 distinguishable, complete field vectors of which 10 are reversed, 28 are transitional, and 13 are normal. The record is complex, quite unlike that predicted by simple flooding or standing nondipole field models. It begins with an estimated several thousand years of reversed polarity with an average intensity of 31.5±8.5 μT, about one third lower than the expected Miocene intensity. This difference is interpreted as a long‐term reduction of the dipole moment prior to the reversal. When site directions and intensities are considered, truly transitional directions and intensities appear almost at the same time at the beginning of the transition, and they disappear simultaneously at the end of the reversal. Large deviations in declination occur during this approximately 4500±1000 year transition period that are compatible with roughly similar average magnitudes of zonal and nonzonal field components at the site. The transitional intensity is generally low, with an average of 10.9±4.9 μT for directions more than 45° away from the dipole field and a minimum of about 5 μT. The root‐mean‐square of the three field components X , Y , and Z are of the same order of magnitude for the transitional field and the historical nondipole field at the site latitude. However, a field intensity increase to pretransitional values occurs when the field temporarily reaches normal directions, which suggests that dipolar structure could have been briefly regenerated during the transition in an aborted attempt to reestablish a stationary field. Changes in the field vector are progressive but jerky, with at least two, and possibly three, large swings at astonishingly high rates. Each of those transitional geomagnetic impulses occurs when the field intensity is low (less than 10 μT) and is followed by an interval of directional stasis during which the magnitude of the field increases greatly. For the best documented geomagnetic impulse the rapid directional change corresponds to a vectorial intensity change of 6700±2700 nT yr −1 , which is about 15–50 times larger than the maximum rate of change of the nondipole field observed during the last centuries. The occurrence of geomagnetic impulses seems to support reversal models assuming an increase in the level of turbulence within the liquid core during transitions. The record closes with an estimated several thousand years of normal polarity with an average intensity of 46.7±20.1 μT, agreeing with the expected Miocene value. However, the occurrence of rather large and apparently rapid intensity fluctuations accompanied by little change in direction suggests that the newly reestablished dipole was still somewhat unstable.

Journal of Geophysical Research B: Solid Earth↗

USGS2022

Assembling a safe and effective toolbox for integrated flea control and plague mitigation: Fipronil experiments with prairie dogs

Background Plague, a widely distributed zoonotic disease of mammalian hosts and flea vectors, poses a significant risk to ecosystems throughout much of Earth. Conservation biologists use insecticides for flea control and plague mitigation. Here, we evaluate the use of an insecticide grain bait, laced with 0.005% fipronil (FIP) by weight, with black-tailed prairie dogs (BTPDs, Cynomys ludovicianus ). We consider safety measures, flea control, BTPD body condition, BTPD survival, efficacy of plague mitigation, and the speed of FIP grain application vs. infusing BTPD burrows with insecticide dusts. We also explore conservation implications for endangered black-footed ferrets ( Mustela nigripes ), which are specialized predators of Cynomys . Principal findings During 5- and 10-day laboratory trials in Colorado, USA, 2016–2017, FIP grain had no detectable acute toxic effect on 20 BTPDs that readily consumed the grain. During field experiments in South Dakota, USA, 2016–2020, FIP grain suppressed fleas on BTPDs for at least 12 months and up to 24 months in many cases; short-term flea control on a few sites was poor for unknown reasons. In an area of South Dakota where plague circulation appeared low or absent, FIP grain had no detectable effect, positive or negative, on BTPD survival. Experimental results suggest FIP grain may have improved BTPD body condition (mass:foot) and reproduction (juveniles:adults). During a 2019 plague epizootic in Colorado, BTPDs on 238 ha habitat were protected by FIP grain, whereas BTPDs were nearly eliminated on non-treated habitat. Applications of FIP grain were 2–4 times faster than dusting BTPD burrows. Significance Deltamethrin dust is the most commonly used insecticide for plague mitigation on Cynomys colonies. Fleas on BTPD colonies exhibit the ability to evolve resistance to deltamethrin after repeated annual treatments. Thus, more tools are needed. Accumulating data show orally-delivered FIP is safe and usually effective for flea control with BTPDs, though potential acute toxic effects cannot be ruled out. With continued study and refinement, FIP might be used in rotation with, or even replace deltamethrin, and serve an important role in Cynomys and black-footed ferret conservation. More broadly, our stepwise approach to research on FIP may function as a template or guide for evaluations of insecticides in the context of wildlife conservation.

South Dakota↗

USGS2010

Geologic history of Mars

Mars accumulated and differentiated into crust, mantle and core within a few tens of millions of years of Solar System formation. Formation of Hellas, which has been adopted as the base of the Noachian period, is estimated to have occurred around 4.1 to 3.8 Gyr ago, depending on whether or not the planet experienced a late cataclysm. Little is known of the pre-Noachian period except that it was characterized by a magnetic field, subject to numerous large basin-forming impacts, probably including one that formed the global dichotomy. The Noachian period, which ended around 3.7 Gyr ago, was characterized by high rates of cratering, erosion, and valley formation. Most of Tharsis formed and surface conditions were at least episodically such as to cause widespread production of hydrous weathering products such as phyllosilicates. Extensive sulfate deposits accumulated late in the era. Average erosion rates, though high compared with later epochs, fell short of the lowest average terrestrial rates. The record suggests that warm, wet conditions necessary for fluvial activity were met only occasionally, such as might occur if caused by large impacts or volcanic eruptions. At the end of the Noachian, rates of impact, valley formation, weathering, and erosion all dropped precipitously but volcanism continued at a relatively high average rate throughout the Hesperian, resulting in the resurfacing of at least 30% of the planet. Large water floods formed episodically, possibly leaving behind large bodies of water. The canyons formed. The observations suggest the change at the end of the Noachian suppressed most aqueous activity at the surface other than large floods, and resulted in growth of a thick cryosphere. However, presence of discrete sulfate rich deposits, sulfate concentrations in soils, and occasional presence of Hesperian valley networks indicates that water activity did not decline to zero. After the end of the Hesperian around 3 Gyr ago the pace of geologic activity slowed further. The average rate of volcanism during the Amazonian was approximately a factor of ten lower than in the Hesperian and activity was confined largely to Tharsis and Elysium. The main era of water flooding was over, although small floods occurred episodically until geologically recent times. Canyon development was largely restricted to formation of large landslides. Erosion and weathering rates remained extremely low. The most distinctive characteristic of the Amazonian is formation of features that have been attributed to the presence, accumulation, and movement of ice. Included are the polar layered deposits, glacial deposits on volcanoes, ice-rich veneers at high latitudes, and a variety of landforms in the 30–55° latitude belts, including lobate debris aprons, lineated valley fill and concentric crater fill. Most of the gullies on steep slopes also formed late in this era. The rate of formation of the ice-related features and the gullies probably varied as changes in obliquity affected the ice stability relations.

Earth and Planetary Science Letters↗

USGS2002

Digital data grids for the magnetic anomaly map of North America

The digital magnetic anomaly database and map for the North American continent is the result of a joint effort by the Geological Survey of Canada (GSC), U. S. Geological Survey (USGS), and Consejo de Recursos Minerales of Mexico (CRM). This integrated, readily accessible, modern digital database of magnetic anomaly data is a powerful tool for further evaluation of the structure, geologic processes, and tectonic evolution of the continent and may also be used to help resolve societal and scientific issues that span national boundaries. The North American magnetic anomaly map derived from the digital database provides a comprehensive magnetic view of continental-scale trends not available in individual data sets, helps link widely separated areas of outcrop, and unifies disparate geologic studies. This open-file report presents three unique, gridded data sets used to make the magnetic anomaly map of North America. Subsets of these three grids that span only the United States were also created, giving a total of six grids. Details on the data processing and compilation procedures used to produce the grids are described in the booklet that accompanies the North American magnetic anomaly map. All three grids have 1-km spacing and are projected to the DNAG projection (spherical transverse mercator, central meridian of 100 o W, base latitude of 0o, scale factor of 0.926 and Earth radius of 6,371,204 m.) More details are given in the metadata files that accompany the gridded data files. These grids are presented in Geosoft binary grid format, with two files describing each of the six grids (suffixes .grd and .gi). This format can be easily converted to numerous other formats using the free conversion software offered by this company at http://www.geosoft.com/. The first grids (NAmag_origmrg.grd and USmag_origmrg.grd) show the magnetic field at 305 m. above terrain. For the second grids (NAmag_hp500.grd and USmag_hp500.grd) we removed long-wavelength anomalies (500 km and greater) from the first grid. This grid was used for the published map. Although the North American merged grid represents a significant upgrade to older compilations, the existing patchwork of surveys is inherently unable to accurately represent anomalies with long (greater than roughly 150 km) wavelengths, particularly in the US and Canada (U.S. Magnetic-Anomaly Data Set Task Group, 1994). The lack of information about long wavelength anomalies is primarily related to datum shifts between merged surveys, caused by data acquisition at widely different times and by differences in merging procedures. Therefore, we removed anomalies with wavelengths greater than 500 km from the merged grid to reduce the effects caused by the spurious long wavelengths but still maintain the continuity of anomalies. The correction was accomplished by transforming the merged grid to the frequency domain, filtering the transformed data with a long-wavelength cutoff at 500 km, and subtracting the long-wavelength data grid from the merged grid. In addition to the 500-km high pass filter, an equivalent source method, based on long-wavelength characterization using satellite data (CHAMP satellite anomalies, Maus and others, 2002), was also used to correct for spurious shifts in the original magnetic anomaly grid (Ravat and others, 2002). These results are presented in the third grids (NAmag_CM.grd and USmag_CM.grd), in which the wavelengths longer than 500 km have been replaced by downward-continued satellite data. The steps used to create the third long-wavelength-corrected grid are: 0. The North American 1-km merged grid was decimated to 5 km. 1. This 5-km grid was converted to a 0.05 degree grid and was low-pass filtered using a Gaussian filter with a 500-km cutoff, then decimated to 1 degree. 2. A joint inversion of this 1-degree low-pass aeromagnetic grid and satellite data, with the aeromagnetic data weighted very low, was used to produce a stabilized downward continuation of the satellite data. 3. The inverted data were interpolated to 0.05 degrees and again low-pass filtered using the same Gaussian 500-km filter to remove short-wavelength artifacts. 4. The low-pass grid from step 1 was subtracted from the original 0.05-degree aeromagnetic grid to create a 500-km high-pass aeromagnetic grid. This grid was added to the low-pass inverted grid from step 3 to get a corrected 0.05-degree aeromagnetic grid. 5. The corrected 0.05-degree aeromagnetic grid was projected to the DNAG projection and regridded to 5 km. This was subtracted from the decimated 5-km aeromagnetic grid to generate a 5-km correction grid. A matched filter was used to remove short-wavelength artifacts resulting from the projection and regridding process. 6. The resulting 5-km correction grid was regridded to the original 1-km grid and subtracted from the original 1-km aeromagnetic grid to generate the final 1-km corrected aeromagnetic grid. The six grids described in this report are available for download. Two metadata files, one for the North American grids and one for the United States grids, are also included with the gridded data.

Open-File Report↗

USGS2000

Stratigraphic and structural evolution of the Selenga Delta Accommodation Zone, Lake Baikal Rift, Siberia

Seismic reflection profiles from the Lake Baikal Rift reveal extensive details about the sediment thickness, structural geometry and history of extensional deformation and syn-rift sedimentation in this classic continental rift. The Selenga River is the largest single source of terrigenous input into Lake Baikal, and its large delta sits astride the major accommodation zone between the Central and South basins of the lake. Incorporating one of the world's largest lacustrine deltas, this depositional system is a classic example of the influence of rift basin structural segmentation on a major continental drainage. More than 3700 km of deep basin-scale multi-channel seismic reflection (MCS) data were acquired during the 1989 Russian and the 1992 Russian–American field programs. The seismic data image most of the sedimentary section, including pre-rift basement in several localities. The MCS data reveal that the broad bathymetric saddle between these two major half-graben basins is underlain by a complex of severely deformed basement blocks, and is not simply a consequence of long-term deltaic deposition. Maximum sediment thickness is estimated to be more than 9 km in some areas around the Selenga Delta. Detailed stratigraphic analyses of the Selenga area MCS data suggest that modes of deposition have shifted markedly during the history of the delta. The present mode of gravity- and mass-flow sedimentation that dominates the northern and southern parts of the modern delta, as well as the pronounced bathymetric relief in the area, are relatively recent developments in the history of the Lake Baikal Rift. Several episodes of major delta progradation, each extending far across the modern rift, can be documented in the MCS data. The stratigraphic framework defined by these prograding deltaic sequences can be used to constrain the structural as well as depositional evolution of this part of the Baikal Rift. An age model has been established for this stratigraphy, by tying the delta sequences to the site of the Baikal Drilling Project 1993 Drill Hole. Although the drill hole is only 100 m deep, and the base of the cores is only ∼670 ka in age, ages were extrapolated to deeper stratigraphic intervals using the Reflection-Seismic-Radiocarbon method of Cohen et al. (1993). The deep prograding delta sequences now observed in the MCS data probably formed in response to major fluctuations in sediment supply, rather than in response to shifts in lake level. This stratigraphic framework and age model suggest that the deep delta packages developed at intervals of approximately 400 ka and may have formed as a consequence of climate changes affiliated with the northern hemisphere glaciations. The stratigraphic analysis also suggests that the Selenga Basin and Syncline developed as a distinct depocentre only during the past ∼2–3 Ma.

Lake Baikal Rift↗

USGS2016

Bayesian estimation of magma supply, storage, and eruption rates using a multiphysical volcano model: Kīlauea Volcano, 2000–2012

Estimating rates of magma supply to the world's volcanoes remains one of the most fundamental aims of volcanology. Yet, supply rates can be difficult to estimate even at well-monitored volcanoes, in part because observations are noisy and are usually considered independently rather than as part of a holistic system. In this work we demonstrate a technique for probabilistically estimating time-variable rates of magma supply to a volcano through probabilistic constraint on storage and eruption rates. This approach utilizes Bayesian joint inversion of diverse datasets using predictions from a multiphysical volcano model, and independent prior information derived from previous geophysical, geochemical, and geological studies. The solution to the inverse problem takes the form of a probability density function which takes into account uncertainties in observations and prior information, and which we sample using a Markov chain Monte Carlo algorithm. Applying the technique to Kīlauea Volcano, we develop a model which relates magma flow rates with deformation of the volcano's surface, sulfur dioxide emission rates, lava flow field volumes, and composition of the volcano's basaltic magma. This model accounts for effects and processes mostly neglected in previous supply rate estimates at Kīlauea, including magma compressibility, loss of sulfur to the hydrothermal system, and potential magma storage in the volcano's deep rift zones. We jointly invert data and prior information to estimate rates of supply, storage, and eruption during three recent quasi-steady-state periods at the volcano. Results shed new light on the time-variability of magma supply to Kīlauea, which we find to have increased by 35–100% between 2001 and 2006 (from 0.11–0.17 to 0.18–0.28 km 3 /yr), before subsequently decreasing to 0.08–0.12 km 3 /yr by 2012. Changes in supply rate directly impact hazard at the volcano, and were largely responsible for an increase in eruption rate of 60–150% between 2001 and 2006, and subsequent decline by as much as 60% by 2012. We also demonstrate the occurrence of temporal changes in the proportion of Kīlauea's magma supply that is stored versus erupted, with the supply “surge” in 2006 associated with increased accumulation of magma at the summit. Finally, we are able to place some constraints on sulfur concentrations in Kīlauea magma and the scrubbing of sulfur by the volcano's hydrothermal system. Multiphysical, Bayesian constraint on magma flow rates may be used to monitor evolving volcanic hazard not just at Kīlauea but at other volcanoes around the world.

Hawai'i↗

USGS1993

Volcanic eruption of the mid-ocean ridge along the East Pacific Rise crest at 9°45-52'N: direct submersible observations of seafloor phenomena associated with an eruption event in April, 1991

In April, 1991, we witnessed from the submersible Alvin a suite of previously undocumented seafloor phenomena accompanying an in-progress eruption of the mid-ocean ridge on the East Pacific Rise crest at 9°45′N–52′N. The volume of the eruption could not be precisely determined, although comparison of pre- and post-eruption SeaBeam bathymetry indicate that any changes in ridge crest morphology resulting from the eruption were < 10 m high. Effects of the eruption included: (1) increased abundance and redistribution of hydrothermal vents, disappearance of numerous vent communities, and changes in characteristics of vent fauna and mineral deposits within the eruption area since December, 1989; (2) murkiness of bottom waters up to tens of meters above the seafloor due to high densities of suspended mineral and biogenic particulates; (3) destruction of a vent community by lava flows, mass wasting, and possible hydrovolcanic explosion at a site known as ‘Tubeworm Barbecue’ in the axial summit caldera (ASC) at 9°50.6′N; (4) near-critical temperatures of hydrothermal vent fluids, ranging up to 403°C; (5) temporal variations over a 2 week interval in both temperatures and chemical/isotopic compositions of hydrothermal fluids; (6) unusual compositions of end-member vent fluids, with pH values ranging to a record low of 2.5, salinities ranging as low as 0.3 wt% NaCl (one-twelfth that of seawater), and dissolved gases reaching high concentrations (> 65 mmol/l for both CO 2 and H 2 S); (7) venting at temperatures above 380°C of visually detectable white vapor that transformed to plumes of gray smoke a few centimeters above vent orifices; (8) disorganized venting of both high-temperature fluids (black and gray smoke) and large volumes of cooler, diffuse hydrothermal fluids directly from the basaltic seafloor, rather than from hydrothermal mineral constructions; (9) rapid and extensive growth of flocculent white bacterial mats (species unknown) on and under the seafloor in areas experiencing widespread venting of diffuse hydrothermal fluid; and (10) subseafloor downslope migration of magma normal to the ridge axis in a network of small-scale (1–5 m diameter) lava tubes and channels to distances at least 100–200 m outside the ASC. We suggest that, in April, 1991, intrusion of dikes in the eruption area to < 200 m beneath the ASC floor resulted in phase separation of fluids near the tops of the dikes and a large flux of vapor-rich hydrothermal fluids through the overlying rubbly, cavernous lavas. Low salinities and gas-rich compositions of hydrothermal fluids sampled in the eruption area are appropriate for a vapor phase in a seawater system undergoing subcritical liquid-vapor phase separation (boiling) and phase segregation. Hydrothermal fluids streamed directly from fissures and pits that may have been loci of lava drainback and/or hydrovolcanic explosions. These fissures and pits were lined with white mats of a unique fast-growing bacteria that was the only life associated with the brand-new vents. The prolific bacteria, which covered thousands of square meters on the ridge crest and were also abundant in subseafloor voids, may thrive on high levels of gases in the vapor-rich hydrothermal fluids initially escaping the hydrothermal system. White bacterial particulates swept from the seafloor by hydrothermal vents swirled in an unprecedented biogenic ‘blizzard’ up to 50 m above the bottom. The bacterial proliferation of April, 1991 is likely to be a transient bloom that will be checked quickly either by decline of dissolved gas concentrations in the fluids as rapid heat loss brings about cessation of boiling, and/or by grazing as other organisms are re-established in the biologically devastated area.

Earth and Planetary Science Letters↗

USGS2013

Historical methane hydrate project review

In 1995, U.S. Geological Survey made the first systematic assessment of the volume of natural gas stored in the hydrate accumulations of the United States. That study, along with numerous other studies, has shown that the amount of gas stored as methane hydrates in the world greatly exceeds the volume of known conventional gas resources. However, gas hydrates represent both a scientific and technical challenge and much remains to be learned about their characteristics and occurrence in nature. Methane hydrate research in recent years has mostly focused on: (1) documenting the geologic parameters that control the occurrence and stability of gas hydrates in nature, (2) assessing the volume of natural gas stored within various gas hydrate accumulations, (3) analyzing the production response and characteristics of methane hydrates, (4) identifying and predicting natural and induced environmental and climate impacts of natural gas hydrates, and (5) analyzing the effects of methane hydrate on drilling safety. Methane hydrates are naturally occurring crystalline substances composed of water and gas, in which a solid water-‐lattice holds gas molecules in a cage-‐like structure. The gas and water becomes a solid under specific temperature and pressure conditions within the Earth, called the hydrate stability zone. Other factors that control the presence of methane hydrate in nature include the source of the gas included within the hydrates, the physical and chemical controls on the migration of gas with a sedimentary basin containing methane hydrates, the availability of the water also included in the hydrate structure, and the presence of a suitable host sediment or “reservoir”. The geologic controls on the occurrence of gas hydrates have become collectively known as the “methane hydrate petroleum system”, which has become the focus of numerous hydrate research programs. Recognizing the importance of methane hydrate research and the need for a coordinated effort, the U.S. Congress enacted Public Law 106-‐193, the Methane Hydrate Research and Development Act of 2000. This Act called for the Secretary of Energy to begin a methane hydrate research and development program in consultation with other U.S. federal agencies. At the same time a new methane hydrate research program had been launched in Japan by the Ministry of International Trade and Industry to develop plans for a methane hydrate exploratory drilling project in the Nankai Trough. Since this early start we have seen other countries including India, China, Canada, and the Republic of Korea establish large gas hydrate research and development programs. These national led efforts have also included the investment in a long list of important scientific research drilling expeditions and production test studies that have provided a wealth of information on the occurrence of methane hydrate in nature. The most notable expeditions and projects have including the following: -‐Ocean Drilling Program Leg 164 (1995) -‐Japan Nankai Trough Project (1999-‐2000) -‐Ocean Drilling Program Leg 204 (2004) -‐Japan Tokai-‐oki to Kumano-‐nada Project (2004) -‐Gulf of Mexico JIP Leg I (2005) -‐Integrated Ocean Drilling Program Expedition 311 (2005) -‐Malaysia Gumusut-‐Kakap Project (2006) -‐India NGHP Expedition 01 (2006) -‐China GMGS Expedition 01 (2007) -‐Republic of Korea UBGH Expedition 01 (2007) -‐Gulf of Mexico JIP Leg II (2009) -‐Republic of Korea UBGH Expedition 02 (2010) -‐MH-‐21 Nankai Trough Pre-‐Production Expedition (2012-‐2013) -‐Mallik Gas Hydrate Testing Projects (1998/2002/2007-‐2008) -‐Alaska Mount Elbert Stratigraphic Test Well (2007) -‐Alaska Iġnik Sikumi Methane Hydrate Production Test Well (2011-‐2012) Research coring and seismic programs carried out by the Ocean Drilling Program (ODP) and Integrated Ocean Drilling Program (IODP), starting with the ODP Leg 164 drilling of the Blake Ridge in the Atlantic Ocean in 1995, have also contributed greatly to our understanding of the geologic controls on the formation, occurrence, and stability of gas hydrates in marine environments. For the most part methane hydrate research expeditions carried out by the ODP and IODP provided the foundation for our scientific understanding of gas hydrates. The methane hydrate research efforts under ODP-‐IODP have mostly dealt with the assessment of the geologic controls on the occurrence of gas hydrate, with a specific goal to study the role methane hydrates may play in the global carbon cycle. Over the last 10 years, national led methane hydrate research programs, along with industry interest have led to the development and execution of major methane hydrate production field test programs. Two of the most important production field testing programs have been conducted at the Mallik site in the Mackenzie River Delta of Canada and in the Eileen methane hydrate accumulation on the North Slope of Alaska. Most recently we have also seen the completion of the world’s first marine methane hydrate production test in the Nankai Trough in the offshore of Japan. Industry interest in gas hydrates has also included important projects that have dealt with the assessment of geologic hazards associated with the presence of hydrates. The scientific drilling and associated coring, logging, and borehole monitoring technologies developed in the long list of methane hydrate related field studies are one of the most important developments and contributions associated with methane hydrate research and development activities. Methane hydrate drilling has been conducted from advanced scientific drilling platforms like the JOIDES Resolution and the D/V Chikyu, which feature highly advanced integrated core laboratories and borehole logging capabilities. Hydrate research drilling has also included the use of a wide array of industry, geotechnical and multi-‐service ships. All of which have been effectively used to collect invaluable geologic and engineering data on the occurrence of methane hydrates throughout the world. Technologies designed specifically for the collection and analysis of undisturbed methane hydrate samples have included the development of a host of pressure core systems and associated specialty laboratory apparatus. The study and use of both wireline conveyed and logging-‐while-‐drilling technologies have also contributed greatly to our understanding of the in-‐situ nature of hydrate-‐bearing sediments. Recent developments in borehole instrumentation specifically designed to monitor changes associated with hydrates in nature through time or to evaluate the response of hydrate accumulations to production have also contributed greatly to our understanding of the complex nature and evolution of methane hydrate systems. Our understanding of how methane hydrates occur and behave in nature is still growing and evolving – we do not yet know if methane hydrates can be economically produced, nor do we know fully the role of hydrates as an agent of climate change or as a geologic hazard. But it is known for certain that scientific drilling has contributed greatly to our understanding of hydrates in nature and will continue to be a critical source of the information to advance our understanding of methane hydrates.

Report↗

USGS2010

Analysis of nonvolcanic tremor on the San Andreas Fault near Parkfield, CA using U.S. Geological Survey Parkfield Seismic Array

Reports by Nadeau and Dolenc (2005) that tremor had been detected near Cholame Valley spawned an effort to use UPSAR (U. S. Geological Survey Parkfield Seismic Array) to study characteristics of tremor. UPSAR was modified to record three channels of velocity at 40–50 sps continuously in January 2005 and ran for about 1 month, during which time we recorded numerous episodes of tremor. One tremor, on 21 January at 0728, was recorded with particularly high signal levels as well as another episode 3 days later. Both events were very emergent, had a frequency content between 2 and 8 Hz, and had numerous high-amplitude, short-duration arrivals within the tremor signal. Here using the first episode as an example, we discuss an analysis procedure, which yields azimuth and apparent velocity of the tremor at UPSAR. We then provide locations for both tremor episodes. The emphasis here is how the tremor episode evolves. Twelve stations were operating at the time of recording. Slowness of arrivals was determined using cross correlation of pairs of stations; the same method used in analyzing the main shock data from 28 September 2004. A feature of this analysis is that 20 s of the time series were used at a time to calculate correlation; the longer windows resulted in more consistent estimates of slowness, but lower peak correlations. These values of correlation (peaks of about 0.25), however, are similar to that obtained for the S wave of a microearthquake. Observed peaks in slowness were traced back to source locations assumed to lie on the San Andreas fault. Our inferred locations for the two tremor events cluster near the locations of previously observed tremor, south of the Cholame Valley. Tremor source depths are in the 14–24 km range, which is below the seismogenic brittle zone, but above the Moho. Estimates of error do not preclude locations below the Moho, however. The tremor signal is very emergent but contains packets that are several times larger than the background tremor signal and lasts about 5 s. These impulsive wavelets are similar to low-frequency earthquakes signals seen in Japan but appear to be broader band rather than just higher in low-frequency energy. They may be more appropriately called high-energy tremor (HET). HET signals at UPSAR correlate well with the record of this event from station GHIB of the HRSN borehole array at Parkfield and HETs typically have a higher cross-correlation coefficient than the rest of the tremor event. The amplitudes of a large HET are consistent with a magnitude of 0.1 when compared with a M2.3 event that had about the same epicenter. Polarizations of the tremor episode at UPSAR are mostly just north of east. Both linearity and azimuth evolve over time suggesting a change in tremor source location over time and linearity is typically higher at the HETs.

California↗

USGS1998

Chemistry, isotopic composition, and origin of a methane-hydrogen sulfide hydrate at the Cascadia subduction zone

Although the presence of extensive gas hydrate on the Cascadia margin, offshore from the western U.S. and Canada, has been inferred from marine seismic records and pore water chemistry, solid gas hydrate has only been found at one location. At Ocean Drilling Program (ODP) Site 892, offshore from central Oregon, gas hydrate was recovered close to the sediment-water interface at 2–19 m below the seafloor (mbsf) at 670 m water depth. The gas hydrate occurs as elongated platy crystals or crystal aggregates, mostly disseminated irregularly, with higher concentrations occurring in discrete zones, thin layers, and/or veinlets parallel or oblique to the bedding. A 2- to 3-cm thick massive gas hydrate layer, parallel to bedding, was recovered at ∼ 17 mbsf. Gas from a sample of this layer was composed of both CH 4 and H 2 S. This sample is the first mixed-gas hydrate of CH 4 H 2 S documented in ODP; it also contains ethane and minor amounts of CO 2 . Measured temperatures of the recovered core ranged from 2 to −1.8°C and are 6 to 8 degrees lower than in-situ temperatures. These temperature anomalies were caused by the partial dissociation of the CH 4 H 2 S hydrate during recovery without a pressure core sampler.During this dissociation, toxic levels of H 2 S (δ 34 S, +27.4‰) were released. The δ 13 C values of the CH 4 in the gas hydrate, −64.5 to −67.5‰ (PDB) , together with δD values of −197 to −199‰ (SMOW) indicate a primarily microbial source for the CH 4 . The δ 18 O value of the hydrate H 2 O is +2.9‰ (SMOW) , comparable with the experimental fractionation factor for sea-ice. The unusual composition (CH 4 H 2 S) and depth distribution (2–19 mbsf) of this gas hydrate indicate mixing between a methane-rich fluid with a pore fluid enriched in sulfide; at this site the former is advecting along an inclined fault into the active sulfate reduction zone. The facts that the CH 4 H 2 S hydrate is primarily confined to the present day active sulfate reduction zone (2–19 mbsf), and that from here down to the BSR depth (19–68 mbsf) the gas hydrate inferred to exist is a ≥99% CH 4 hydrate, suggest that the mixing of CH 4 and H 2 S is a geologically young process. Because the existence of a mixed CH 4 H 2 S hydrate is indicative of moderate to intense advection of a methane-rich fluid into a near surface active sulfate reduction zone, tectonically active (faulted) margins with organic-rich sediments and moderate to high sedimentation rates are the most likely regions of occurrence. The extension of such a mixed hydrate below the sulfate reduction zone should reflect the time-span of methane advection into the sulfate reduction zone.

Earth and Planetary Science Letters↗