USGS Science⌕ Search

SEARCH · USGS Science

Results for “Ocean Modeling”

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

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

At least 343 records · Page 19Linked to original sources
USGS1980

The age curves of sulfur and oxygen isotopes in marine sulfate and their mutual interpretation

Three hundred new samples of marine evaporite sulfate, of world-wide distribution, were analyzed for δ 34 S, and 60 of these also for δ 18 O in the sulfate ion. Detailed δ 34 S age curves for Tertiary—Cretaceous, Permian—Pennsylvanian, Devonian, Cambrian and Proterozoic times document large variations in δ 34 S. A summary curve forδ 18 O also shows definite variations, some at different times than δ 34 S, and always smaller. The measured δ 34 S and δ 18 O correspond to variations in these isotopes in sulfate of the world ocean surface. The variations of δ 18 O are controlled by input and output fluxes of sulfur in the ocean, three of which are the same that control δ 34 S: deposition and erosion of sulfate, and deposition of sulfide. Erosion of sulfide differs in its effect on the S and O systems. δ 18 O in the sulfate does not seem to be measurably affected by equilibration with either seawater or with subsurface waters after crystallization. In principle, the simultaneous application of both δ 34 S and δ 18 O age curves should help reduce the number of assumptions in calculations of the cycles of sulfur and oxygen through geological time, and a new model involving symmetrical fluxes is introduced here to take advantage of the oxygen data. However, all previously published models as well as this one lead to anomalies, such as unreasonable calcium or oxygen depletions in the ocean—atmosphere system. In addition, most models are incapable of reproducing the sharp rises of the δ 34 S curve in the late Proterozoic, the Devonian and the Triassic which would be the result of unreasonably fast net sulfide deposition. This fast depletion could result from an ocean that has not always been mixed (as previously assumed in all model calculations).

Chemical Geology↗
USGS2011

Bioenergetics in ecosystems

A bioenergetics model for a fish can be defined as a quantitative description of the fish’s energy budget. Bioenergetics modeling can be applied to a fish population in a lake, river, or ocean to estimate the annual consumption of food by the fish population; such applications have proved to be useful in managing fisheries. In addition, bioenergetics models have been used to better understand fish growth and consumption in ecosystems, to determine the importance of the role of fish in cycling nutrients within ecosystems, and to identify the important factors regulating contaminant accumulation in fish from lakes, rivers, and oceans.

Book chapter↗
USGS2010

Cenozoic plate motions and the volcano-tectonic evolution of western Oregon and Washington

A refined northeast Pacific plate-motion model provides a framework for analysis of the Tertiary volcanic and tectonic history of western Oregon and Washington. We examine three possible models for the origin of the allochthonous Paleocene and Eocene oceanic basalt basement of the Coast Range: (1) accretion to the continent of hot spot generated linear seamount chains; (2) accretion of thick oceanic crust and seamounts generated during Farallon-Kula spreading reorganizations between 61 and 48 Ma; and (3) eruption of basalt during oblique rifting of the continental margin as it overrode an active Yellowstone hot spot on the Kula-Farallon ridge. The plate model suggests that microplate rotation and accretion of hot spot generated linear aseismic ridges cannot be easily reconciled with rapid northeast motion of the KuIa and Farallon plates and the well-established paleomagnetic rotations. Following emplacement of the Coast Range basement, changes in the character of forearc, back arc and Cascade arc volcanism correlate with a marked decrease in the rate of Farallon-North America convergence between 43 and and 28 Ma. This slowdown may be responsible for (1) westward stepping of the volcanic arc front from the Challis axis to a Cascade axis at about 42 Ma; (2) a subsequent episode of increased ash flow tuff volcanism and extension in the Cascade arc between 37 and 18 Ma that correlates with the “ignimbrite flare-up” in the Basin and Range; and (3) a period of extensional basaltic and alkalic volcanism and intrusion in the Coast Range between 44 and 28 Ma. Reduction of the convergence rate and westward stepping of the flexure in the subducted slab may have reduced the horizontal compressive stress on the continent, allowing increased injection of magma into the crust, development of large, shallow magma chambers, and the outbreak of extensional volcanism over a large area behind the Farallon-North America subduction zone.

Oregon, Washington↗
USGS2018

Drivers of variability in public‐supply water use across the contiguous United States

This study explores the relationship between municipal water use and an array of climate, economic, behavioral, and policy variables across the contiguous U.S. The relationship is explored using Bayesian‐hierarchical regression models for over 2,500 counties, 18 covariates, and three higher‐level grouping variables. Additionally, a second analysis is included for 83 cities where water price and water conservation policy information is available. A hierarchical model using the nine climate regions (product of National Oceanic and Atmospheric Administration) as the higher‐level groups results in the best out‐of‐sample performance, as estimated by the Widely Available Information Criterion, compared to counties grouped by urban continuum classification or primary economic activity. The regression coefficients indicate that the controls on water use are not uniform across the nation: e.g., counties in the Northeast and Northwest climate regions are more sensitive to social variables, whereas counties in the Southwest and East North Central climate regions are more sensitive to environmental variables. For the national city‐level model, it appears that arid cities with a high cost of living and relatively low water bills sell more water per customer, but as with the county‐level model, the effect of each variable depends heavily on where a city is located.

Water Resources Research↗
USGS1986

Crustal structure of Oaxaca, Mexico, from seismic refraction measurements

Seismic refraction and gravity data have been analyzed to obtain a model of the compressional-wave velocity structure of the ocean-to-continent transition in the State of Oaxaca in southwestern Mexico. Crustal thickness on the continent at the latitude 18°N is 45 ± 4 km, based on reflected phases from the Moho discontinuity. The crust has been modeled with three layers, with velocities of 4.3 to 4.6, 5.0 to 5.7, and 6.85 to 7.0 km/sec, each with positive velocity gradient. The crust thins to 10 km at the coast near Pinotepa Nacional, where Precambrian metamorphic rocks are exposed 45 km from the mid-America trench. Offshore, the oceanic structure consists of an 8-km-thick crust with a normal crustal velocity structure (Spudich and Orcutt, 1980). The apparent dip of the subducting plate beneath western Mexico is 10°. On the oceanic side, strong reflections suggest a minimum depth of 35 km for the lithosphere-asthenosphere boundary. The asthenosphere has a seismic velocity of 7.6 km/sec, and a thin lid in which the velocity is 8.6 km/sec.

Oaxaca↗
USGS2019

Application of the Stream Salmonid Simulator (S3) to Klamath River fall Chinook salmon (Oncorhynchus tshawytscha), California—Parameterization and calibration

Executive Summary In this report, we describe application of the Stream Salmonid Simulator (S3) to Chinook salmon ( Oncorhynchus tshawytscha ) in the Klamath River between Keno Dam in southern Oregon and the ocean in northern California. S3 is a deterministic life-stage-structured population model that tracks daily growth, movement, and survival of juvenile salmon. It can track different source populations or species, such as major tributary populations that enter a river like the Klamath River. A key theme of the model is that river flow affects habitat availability and capacity, which in turn drives density-dependent population dynamics. To explicitly link population dynamics to habitat quality and quantity, the river environment is constructed as a one-dimensional series of linked habitat units, each of which has an associated daily time series of discharge, water temperature, and useable habitat area or carrying capacity. In turn, the physical characteristics of each habitat unit and the number of fish occupying each unit affect survival and growth within each habitat unit and movement of fish among habitat units. The physical template of the Klamath River was formed by classifying the river into 2,635 mesohabitat units composed of runs, riffles, and pools. This template enabled modeling of the unimpounded Klamath River between the Keno Dam (the uppermost of four dams) and Iron Gate Dam (the lowermost dam) to address dam-removal scenarios. However, in this report, our focus was on parameterizing and calibrating the model under existing conditions, which included 1,706 discrete habitat units over the 312-kilometer (km) section of river between Iron Gate Dam and the ocean. For each habitat unit, we developed a time series of daily flow, water temperature, and amount of available habitat (weighted usable habitat area [WUA]) for spawners, fry, and parr. WUA time series were constructed using habitat suitability criteria for Chinook salmon applied to eight two-dimensional (2-D) hydrodynamic models that represented the geomorphic variability in habitat across the Klamath River. Results from the 2-D models were then extrapolated to unmodeled habitat units by scaling WUA curves for changes in habitat unit length and width. These variables were then used to drive population dynamics such as egg development and survival and juvenile movement, growth, and survival.

California, Oregon↗
USGS2007

Circulation and physical processes within the San Gabriel River Estuary during summer 2005

The Southern California Coastal Water Research Project (SCCWRP) is developing a hydrodynamic model of the SGR estuary, which is part of the comprehensive water-quality model of the SGR estuary and watershed investigated by SCCWRP and other local agencies. The hydrodynamic model will help understanding of 1) the exchange processes between the estuary and coastal ocean; 2) the circulation patterns in the estuary; 3) upstream natural runoff and the cooling discharge from PGS. Like all models, the SGR hydrodynamic model is only useful after it is fully calibrated and validated. In May 2005, SCCWRP requested the assistance of the U.S. geological Survey (USGS) Coastal and Marine Geology team (CMG) in collecting data on the hydrodynamic conditions in the estuary during the summer dry season. The summer was chosen for field data collection as this was assumed to be the season with the greatest potential for chronic degraded water quality due to low river flow and high thermal stratification within the estuary (due to both higher average air temperature and PGS output). Water quality can be degraded in winter as well, when higher river discharge events bring large volumes of water from the Los Angeles basin into the estuary. The objectives of this project were to 1) collect hydrodynamic data along the SGR estuary; 2) study exchange processes within the estuary through analysis of the hydrodynamic data; and 3) provide field data for model calibration and validation. As the data only exist for the summer season, the results herein only apply to summer conditions.

California↗
USGS2013

Modeling volcano growth on the Island of Hawaii: Deep-water perspectives

Recent ocean-bottom geophysical surveys, dredging, and dives, which complement surface data and scientific drilling at the Island of Hawaii, document that evolutionary stages during volcano growth are more diverse than previously described. Based on combining available composition, isotopic age, and geologically constrained volume data for each of the component volcanoes, this overview provides the first integrated models for overall growth of any Hawaiian island. In contrast to prior morphologic models for volcano evolution (preshield, shield, postshield), growth increasingly can be tracked by age and volume (magma supply), defining waxing alkalic, sustained tholeiitic, and waning alkalic stages. Data and estimates for individual volcanoes are used to model changing magma supply during successive compositional stages, to place limits on volcano life spans, and to interpret composite assembly of the island. Volcano volumes vary by an order of magnitude; peak magma supply also varies sizably among edifices but is challenging to quantify because of uncertainty about volcano life spans. Three alternative models are compared: (1) near-constant volcano propagation, (2) near-equal volcano durations, (3) high peak-tholeiite magma supply. These models define inconsistencies with prior geodynamic models, indicate that composite growth at Hawaii peaked ca. 800–400 ka, and demonstrate a lower current rate. Recent age determinations for Kilauea and Kohala define a volcano propagation rate of 8.6 cm/yr that yields plausible inception ages for other volcanoes of the Kea trend. In contrast, a similar propagation rate for the less-constrained Loa trend would require inception of Loihi Seamount in the future and ages that become implausibly large for the older volcanoes. An alternative rate of 10.6 cm/yr for Loa-trend volcanoes is reasonably consistent with ages and volcano spacing, but younger Loa volcanoes are offset from the Kea trend in age-distance plots. Variable magma flux at the Island of Hawaii, and longer-term growth of the Hawaiian chain as discrete islands rather than a continuous ridge, may record pulsed magma flow in the hotspot/plume source.

Hawai'i↗
USGS2022

Systematic mapping of the ocean-continent transform plate boundary of the Queen Charlotte fault system, southeastern Alaska and western British Columbia—A preliminary bathymetric terrain model

In 2015, U.S. Geological Survey scientists in collaboration with scientists from other institutions began a study of the Queen Charlotte fault—the first systematic study of the fault in more than three decades. The primary goal of the study was to gain a better understanding of the earthquake, tsunami, and underwater-landslide hazards throughout southeastern Alaska, as well as gather data to develop geologic models that can be applied to similar plate boundaries around the globe, such as the San Andreas fault system in southern California, the Alpine fault in New Zealand, and the North Anatolian fault in Turkey. A bathymetric terrain model was compiled from six different multibeam surveys of the previously unmapped Queen Charlotte fault offshore of southeastern Alaska and Haida Gwaii archipelago.

Alaska, British Columbia↗
USGS1989

Origin of the oceanic basalt basement of the Solomon Islands arc and its relationship to the Ontong Java Plateau-insights from Cenozoic plate motion models

Cenozoic global plate motion models based on a hotspot reference frame may provide a useful framework for analyzing the tectonic evolution of the Solomon Islands convergent margin. A postulated late Miocene collision of the Ontong Java Plateau (OJP) with a NE-facing arc is consistent with the predicted path of the OJP across the Pacific Basin and its Miocene arrival at the trench. Late-stage igneous activity (65-30 Ma) predicted for the OJP as it rode over the Samoan hotspot occurred in correlative stratigraphic sections on Malaita, the supposed accreted flake of OJP in the Solomon Islands arc. Convergence similar to the present velocities between Australia and the Pacific plates was characteristic of the last 43 million years. Prior to 43 Ma Pacific-Australia plate motions were divergent, seemingly at odds with geologic evidence for early Tertiary convergence, particularly in Papua New Guinea. A postulated South Pacific plate may have existed between Australia and the Pacific plate and would have allowed implied northward subduction along the northeastern Australia plate boundary that lasted into the early Eocene. Subsequent reorganization of plate motions in the middle Eocene correlates with middle Eocene marginal basin formation along ridges oblique to the main plate boundary. Cessation of spreading on the Pacific-South Pacific Ridge and its subsequent subduction beneath Asia followed the change in Pacific plate motion at 43 Ma. A trapped remnant of the extinct, NW-trending ridge may still lie beneath the western Philippine Sea. The terminal deformation, metamorphism and ophiolite obduction in the Eocene orogen of the southwest Pacific also correlates with the major change in Pacific plate motion at 43 Ma and the subsequent compression of the dying Eocene arc against outlying continental and oceanic crustal blocks of the Australian plate. The Solomon Islands oceanic basement may represent juxtaposition of oceanic plateaus of the Australian plate beneath overthrust, dismembered ophiolite derived from adjacent marginal basin crust.

Tectonophysics↗
USGS2018

Global and Arctic climate sensitivity enhanced by changes in North Pacific heat flux

Arctic amplification is a consequence of surface albedo, cloud, and temperature feedbacks, as well as poleward oceanic and atmospheric heat transport. However, the relative impact of changes in sea surface temperature (SST) patterns and ocean heat flux sourced from different regions on Arctic temperatures are not well constrained. We modify ocean-to-atmosphere heat fluxes in the North Pacific and North Atlantic in a climate model to determine the sensitivity of Arctic temperatures to zonal heterogeneities in northern hemisphere SST patterns. Both positive and negative ocean heat flux perturbations from the North Pacific result in greater global and Arctic surface air temperature anomalies than equivalent magnitude perturbations from the North Atlantic; a response we primarily attribute to greater moisture flux from the subpolar extratropics to Arctic. Enhanced poleward latent heat and moisture transport drive sea-ice retreat and low-cloud formation in the Arctic, amplifying Arctic surface warming through the ice-albedo feedback and infrared warming effect of low clouds. Our results imply that global climate sensitivity may be dependent on patterns of ocean heat flux in the northern hemisphere.

Nature Communications↗
USGS2024

Constraints from GPS measurements on plate coupling within the Makran subduction zone and tsunami scenarios in the western Indian Ocean

Plate-coupling estimates and previous seismicity indicate that portions of the Makran megathrust of southern Pakistan and Iran are partially coupled and have the potential to produce future magnitude 7+ earthquakes. However, the GPS observations needed to constrain coupling models are sparse and lead to an incomplete understanding of regional earthquake and tsunami hazard. In this study, we assess GPS velocities for plate coupling of the Makran subduction zone with specific attention to model resolution and the accretionary prism rheology. We use finite element model-derived Green's functions to invert for the interseismic slip deficit under both elastic and viscoelastic Earth assumptions. We use the model resolution matrix to characterize plate-coupling scenarios that are consistent with the limited spatial resolution afforded by GPS observations. We then forward model the corresponding tsunami responses at major coastal cities within the western Indian Ocean basin. Our plate-coupling results show potential segmentation of the megathrust with varying coupling from west to east, but do not rule out a scenario where the entire length of the megathrust could rupture in a single earthquake. The full subduction zone rupture scenarios suggest that the Makran may be able to produce earthquakes up to M w 9.2. The corresponding tsunami model from the largest earthquake event ( M w 9.2) estimates maximum wave heights reaching 2–5 m at major port cities in the northern Arabian Sea region. Cities on the west coast of India are less affected (1–2 m). Coastlines bounding eastern Africa, and the Strait of Hormuz, are the least affected (<1 m).

Geophysical Journal International↗
USGS2019

A seascape-scale habitat model to support management of fishing impacts on benthic ecosystems

Minimizing fishing impacts on seafloor ecosystems is a growing focus of ocean management; however, few quantitative tools exist to guide seascape-scale habitat management. To meet these needs, we developed a model to assess benthic ecosystem impacts from fishing gear contact. The habitat impacts model is cast in discrete time and can accommodate overlapping fisheries as well as incorporate gear-specific contact dynamics. We implemented the model in the North Pacific using fishing data from 2003 to 2017, estimating that habitat in 3.1% of the 1.2 million km 2 study area was disturbed at the end of the simulation period. A marked decline in habitat disturbance was evident since 2010, attributable to a single regulatory gear change that lifted trawl gear components off the seafloor. Running scenarios without these gear modifications showed these policies might have contributed to a 24% reduction in habitat disturbance since their implementation. Ultimately, model outputs provide direct estimates of the spatial and temporal trends of habitat effects from fishing — a key component of regulatory policies for many of the world’s fisheries.

Canadian Journal of Fisheries and Aquatic Sciences↗
USGS2005

River nutrient loads and catchment size

We have used a total of 496 sample sites to calibrate a simple regression model for calculating dissolved inorganic nutrient fluxes via runoff to the ocean. The regression uses the logarithms of runoff and human population as the independent variables and estimates the logarithms of dissolved inorganic nitrogen and phosphorus loading with R 2 values near 0.8. This predictive capability is about the same as has been derived for total nutrient loading with process-based models requiring more detailed information on independent variables. We conclude that population and runoff are robust proxies for the more detailed application, landscape modification, and in-stream processing estimated by more process-based models. The regression model has then been applied to a demonstration data set of 1353 river catchments draining to the sea from the North American continent south of the Canadian border. The geographic extents of these basins were extracted from a 1-km digital elevation model for North America, and both runoff and population were estimated for each basin. Most of the basins (72% of the total) are smaller than 103 km2, and both runoff and population density are higher and more variable among small basins than among larger ones.While total load to the ocean can probably be adequately estimated from large systems only, analysis of the geographic distribution of nutrient loading requires consideration of the small basins, which can exhibit significant hydrologic and demographic heterogeneity between systems over their range even within the same geographic region. High-resolution regional and local analysis is necessary for environmental assessment and management. ?? Springer 2005.

Biogeochemistry↗
USGS2019

3D geologic framework for use with the U.S. Geological Survey National Crustal Model, Phase 1—Western United States

A 3D geologic framework is presented here as part of the U.S. Geological Survey National Crustal Model for the western United States, which will be used to improve seismic hazard assessment. The framework is based on 1:250,000 to 1:1,000,000-scale state geologic maps and depths of multiple subsurface unit boundaries. The geology at or near the Earth’s surface is based on published maps with modifications to remove discontinuities across state borders. Extrapolation of rock type and age in the subsurface is achieved by iterative stripping of units of a given age, nearest neighbor interpolation of the remaining units, and constraints on basement geology. The subsurface depth of the interfaces between units is determined by a range of models with varying quantity and quality of constraints. Bedrock depth is derived primarily from a proxy model with added geophysical constraint in some areas. The depths to the base of Cenozoic and Phanerozoic sedimentary and extrusive volcanic rocks are constrained by geophysical methods in many areas. Elsewhere, a simple method is used to estimate their subsurface depth based on the distance to the edge of the geologic units. The remaining continental units are evenly distributed above, below, and between depending on age. The oceanic crust is treated as a simple four-layer model with the added complexity of subduction beneath the North American plate along the Cascadia subduction zone. Refinements to this technique may be accomplished in future versions of the model with more specific information including the location of faults to produce discontinuities in geologic structure and additional information obtained from boreholes and geophysical studies. Further improvements to the geologic framework may be made by incorporating information from more local studies, for example, hydrogeologic studies.

Open-File Report↗
USGS1993

Modeling the tides of Massachusetts and Cape Cod Bays

A time-dependent, three-dimensional numerical modeling study of the tides of Massachusetts and Cape Code Bays, motivated by construction of a new sewage treatment plant and ocean outfall for the city of Boston, has been undertaken by the authors. The numerical model being used is a hybrid version of the Blumberg and Mellor ECOM3D model, modified to include a semi-implicit time-stepping scheme and transport of a non-reactive dissolved constituent. Tides in the bays are dominated by the semi-diurnal frequencies, in particular by the M2 tide, due to the resonance of these frequencies in the Gulf of Maine. The numerical model reproduces, well, measured tidal ellipses in unstratified wintertime conditions. Stratified conditions present more of a problem because tidal-frequency internal wave generation and propagation significantly complicates the structure of the resulting tidal field. Nonetheless, the numerical model reproduces qualitative aspects of the stratified tidal flow that are consistent with observations in the bays.

Conference Paper↗
USGS2017

Inner-shelf ocean dynamics and seafloor morphologic changes during Hurricane Sandy

Hurricane Sandy was one of the most destructive hurricanes in US history, making landfall on the New Jersey coast on Oct 30, 2012. Storm impacts included several barrier island breaches, massive coastal erosion, and flooding. While changes to the subaerial landscape are relatively easily observed, storm-induced changes to the adjacent shoreface and inner continental shelf are more difficult to evaluate. These regions provide a framework for the coastal zone, are important for navigation, aggregate resources, marine ecosystems, and coastal evolution. Here we provide unprecedented perspective regarding regional inner continental shelf sediment dynamics based on both observations and numerical modeling over time scales associated with these types of large storm events. Oceanographic conditions and seafloor morphologic changes are evaluated using both a coupled atmospheric-ocean-wave-sediment numerical modeling system and observation analysis from a series of geologic surveys and oceanographic instrument deployments focused on a region offshore of Fire Island, NY. The geologic investigations conducted in 2011 and 2014 revealed lateral movement of sedimentary structures of distances up to 450 m and in water depths up to 30 m, and vertical changes in sediment thickness greater than 1 m in some locations. The modeling investigations utilize a system with grid refinement designed to simulate oceanographic conditions with progressively increasing resolutions for the entire US East Coast (5-km grid), the New York Bight (700-m grid), and offshore of Fire Island, NY (100-m grid), allowing larger scale dynamics to drive smaller scale coastal changes. Model results in the New York Bight identify maximum storm surge of up to 3 m, surface currents on the order of 2 ms-1 along the New Jersey coast, waves up to 8 m in height, and bottom stresses exceeding 10 Pa. Flow down the Hudson Shelf Valley is shown to result in convergent sediment transport and deposition along its axis. Modeled sediment redistribution along Fire Island showed erosion across the crests of inner shelf sand ridges and sedimentation in adjacent troughs, consistent with the geologic observations.

Continental Shelf Research↗
USGS2017

The interaction of climate change and methane hydrates

Gas hydrate, a frozen, naturally-occurring, and highly-concentrated form of methane, sequesters significant carbon in the global system and is stable only over a range of low-temperature and moderate-pressure conditions. Gas hydrate is widespread in the sediments of marine continental margins and permafrost areas, locations where ocean and atmospheric warming may perturb the hydrate stability field and lead to release of the sequestered methane into the overlying sediments and soils. Methane and methane-derived carbon that escape from sediments and soils and reach the atmosphere could exacerbate greenhouse warming. The synergy between warming climate and gas hydrate dissociation feeds a popular perception that global warming could drive catastrophic methane releases from the contemporary gas hydrate reservoir. Appropriate evaluation of the two sides of the climate-methane hydrate synergy requires assessing direct and indirect observational data related to gas hydrate dissociation phenomena and numerical models that track the interaction of gas hydrates/methane with the ocean and/or atmosphere. Methane hydrate is likely undergoing dissociation now on global upper continental slopes and on continental shelves that ring the Arctic Ocean. Many factors—the depth of the gas hydrates in sediments, strong sediment and water column sinks, and the inability of bubbles emitted at the seafloor to deliver methane to the sea-air interface in most cases—mitigate the impact of gas hydrate dissociation on atmospheric greenhouse gas concentrations though. There is no conclusive proof that hydrate-derived methane is reaching the atmosphere now, but more observational data and improved numerical models will better characterize the climate-hydrate synergy in the future.

Reviews of Geophysics↗