Search USGSSearch

SEARCH · Search USGS

Results for “Journal of Geophysical Research: Atmospheres”

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 181 records · Page 10Linked to original sources

Key components and contrasts in the nitrogen budget across a US-Canadian transboundary watershed

Watershed nitrogen (N) budgets provide insights into drivers and solutions for groundwater and surface water N contamination. We constructed a comprehensive N budget for the transboundary Nooksack River Watershed (British Columbia, Canada, and Washington, USA) using locally derived data, national statistics, and standard parameters. Feed imports for dairy (mainly in the United States) and poultry (mainly in Canada) accounted for 30% and 29% of the total N input to the watershed, respectively. Synthetic fertilizer was the next largest source contributing 21% of inputs. Food imports for humans and pets together accounted for 9% of total inputs, lower than atmospheric deposition (10%). N imported by returning salmon representing marine‐derived nutrients accounted for <0.06% of total N input. Quantified N export was 80% of total N input, driven by ammonia emission (32% of exports). Animal product export was the second largest output of N (31%) as milk and cattle in the United States and poultry products in Canada. Riverine export of N was estimated at 28% of total N export. The commonly used crop nitrogen use efficiency (NUE) metric alone did not provide sufficient information on farming activities but in combination with other criteria such as farm‐gate NUE may better represent management efficiency. Agriculture was the primary driver of N inputs to the environment as a result of its regional importance; the N budget information can inform management to minimize N losses. The N budget provides key information for stakeholders across sectors and borders to create environmentally and economically viable and effective solutions.

Washington, British Columbia

Responses of marginal and intrinsic water-use efficiency to changing aridity using FLUXNET observations

According to classic stomatal optimization theory, plant stomata are regulated to maximize carbon assimilation for a given water loss. A key component of stomatal optimization models is marginal water-use efficiency (mWUE), the ratio of the change of transpiration to the change in carbon assimilation. Although the mWUE is often assumed to be constant, variability of mWUE under changing hydrologic conditions has been reported. However, there has yet to be a consensus on the patterns of mWUE variabilities and their relations with atmospheric aridity. We investigate the dynamics of mWUE in response to vapor pressure deficit (VPD) and aridity index using carbon and water fluxes from 115 eddy covariance towers available from the global database FLUXNET. We demonstrate a non-linear mWUE-VPD relationship at a sub-daily scale in general; mWUE varies substantially at both low and high VPD levels. However, mWUE remains relatively constant within the mid-range of VPD. Despite the highly non-linear relationship between mWUE and VPD, the relationship can be informed by the strong linear relationship between ecosystem-level inherent water-use efficiency (IWUE) and mWUE using the slope, m *. We further identify site-specific m * and its variability with changing site-level aridity across six vegetation types. We suggest accurately representing the relationship between IWUE and VPD using Michaelis–Menten or quadratic functions to ensure precise estimation of mWUE variability for individual sites.

Journal of Geophysical Research Biogeosciences

Analysis of time series of glacier speed: Columbia Glacier, Alaska

During the summers of 1984 and 1985, laser measurements were made of the distance from a reference location to markers on the surface of the lower reach of Columbia Glacier, Alaska. The data have numerous gaps, mostly because of inclement weather. The laser measurements were corrected for variations in atmospheric temperature and pressure and for the trajectory of the marker. The marker speed was calculated by using a combination of cubic splines and digital filters to arrive at a data set with 1-hour time intervals. The speed varies from 7 to 15 m/d and has three noteworthy components: (1) a low-frequency perturbation in speed with a time scale of days related to increased precipitation, (2) semidiurnal and diurnal variations related to sea tides, and (3) diurnal variations related to glacier surface melt. The low frequency and tidal period variations are separated by low-pass filtering and subtracting the filtered data from the original data. The variations from melt runoff and tides have the same frequency range and are examined using harmonic analysis and multiple regression with the various forcing functions. The high-frequency portion of the ice speed signal is dominated by tidal influences and to a lesser extent by meltwater influences. The low-frequency portion of the signal is dominated by the effects of precipitation.

Journal of Geophysical Research Solid Earth

Wind-driven particle mobility on Mars: Insights from Mars Exploration Rover observations at "El Dorado" and surroundings at Gusev Crater

The ripple field known as 'El Dorado' was a unique stop on Spirit's traverse where dust-raising, active mafic sand ripples and larger inactive coarse-grained ripples interact, illuminating several long-standing issues of Martian dust mobility, sand mobility, and the origin of transverse aeolian ridges. Strong regional wind events endured by Spirit caused perceptible migration of ripple crests in deposits SSE of El Dorado, erasure of tracks in sandy areas, and changes to dust mantling the site. Localized thermal vortices swept across El Dorado, leaving paths of reduced dust but without perceptibly damaging nearly cohesionless sandy ripple crests. From orbit, winds responsible for frequently raising clay-sized dust into the atmosphere do not seem to significantly affect dunes composed of (more easily entrained) sand-sized particles, a long-standing paradox. This disparity between dust mobilization and sand mobilization on Mars is due largely to two factors: (1) dust occurs on the surface as fragile, low-density, sand-sized aggregates that are easily entrained and disrupted, compared with clay-sized air fall particles; and (2) induration of regolith is pervasive. Light-toned bed forms investigated at Gusev are coarse-grained ripples, an interpretation we propose for many of the smallest linear, light-toned bed forms of uncertain origin seen in high-resolution orbital images across Mars. On Earth, wind can organize bimodal or poorly sorted loose sediment into coarse-grained ripples. Coarse-grained ripples could be relatively common on Mars because development of durable, well-sorted sediments analogous to terrestrial aeolian quartz sand deposits is restricted by the lack of free quartz and limited hydraulic sediment processing. Copyright 2008 by the American Geophysical Union.

Journal of Geophysical Research E: Planets

Aeolian dunes as ground truth for atmospheric modeling on Mars

Martian aeolian dunes preserve a record of atmosphere/surface interaction on a variety of scales, serving as ground truth for both Global Climate Models (GCMs) and mesoscale climate models, such as the Mars Regional Atmospheric Modeling System (MRAMS). We hypothesize that the location of dune fields, expressed globally by geographic distribution and locally by dune centroid azimuth (DCA), may record the long-term integration of atmospheric activity across a broad area, preserving GCM-scale atmospheric trends. In contrast, individual dune morphology, as expressed in slipface orientation (SF), may be more sensitive to localized variations in circulation, preserving topographically controlled mesoscale trends. We test this hypothesis by comparing the geographic distribution, DCA, and SF of dunes with output from the Ames Mars GCM and, at a local study site, with output from MRAMS. When compared to the GCM: 1) dunes generally lie adjacent to areas with strongest winds, 2) DCA agrees fairly well with GCM modeled wind directions in smooth-floored craters, and 3) SF does not agree well with GCM modeled wind directions. When compared to MRAMS modeled winds at our study site: 1) DCA generally coincides with the part of the crater where modeled mean winds are weak, and 2) SFs are consistent with some weak, topographically influenced modeled winds. We conclude that: 1) geographic distribution may be valuable as ground truth for GCMs, 2) DCA may be useful as ground truth for both GCM and mesoscale models, and 3) SF may be useful as ground truth for mesoscale models. Copyright 2009 by the American Geophysical Union.

Journal of Geophysical Research E: Planets

Surface methane concentrations along the mid-Atlantic bight driven by aerobic subsurface production rather than seafloor gas seeps

Relatively minor amounts of methane, a potent greenhouse gas, are currently emitted from the oceans to the atmosphere, but such methane emissions have been hypothesized to increase as oceans warm. Here, we investigate the source, distribution, and fate of methane released from the upper continental slope of the U.S. Mid-Atlantic Bight, where hundreds of gas seeps have been discovered between the shelf-break and ~1600 m water depth. Using physical, chemical, and isotopic analyses, we identify two main sources of methane in the water column: seafloor gas seeps and in situ aerobic methanogenesis which primarily occurs at 100 – 200 m depth in the water column. Stable isotopic analyses reveal that water samples collected at all depths were significantly impacted by aerobic methane oxidation, the dominant methane sink in this region, with more than 50% of the methane being oxidized, on average. Due to methane oxidation in the deeper water column, below 200 m depth, surface concentrations of methane are influenced more by methane sources found near the surface (0 – 10 m depth) and in the subsurface (10 - 200 m depth), rather than seafloor emissions at greater depths.

Atlantic margin

Barrier island breach evolution: Alongshore transport and bay-ocean pressure gradient interactions

Physical processes controlling repeated openings and closures of a barrier island breach between a bay and the open ocean are studied using aerial photographs and atmospheric and hydrodynamic observations. The breach site is located on Pea Island along the Outer Banks, separating Pamlico Sound from the Atlantic Ocean. Wind direction was a major control on the pressure gradients between the bay and the ocean to drive flows that initiate or maintain the breach opening. Alongshore sediment flux was found to be a major contributor to breach closure. During the analysis period from 2011 to 2016, three hurricanes had major impacts on the breach. First, Hurricane Irene opened the breach with wind-driven flow from bay to ocean in August 2011. Hurricane Sandy in October 2012 quadrupled the channel width from pressure gradient flows due to water levels that were first higher on the ocean side and then higher on the bay side. The breach closed sometime in Spring 2013, most likely due to an event associated with strong alongshore sediment flux but minimal ocean-bay pressure gradients. Then, in July 2014, Hurricane Arthur briefly opened the breach again from the bay side, in a similar fashion to Irene. In summary, opening and closure of breaches are shown to follow a dynamic and episodic balance between along-channel pressure gradient driven flows and alongshore sediment fluxes.

Pea Island

A comparison of Martian north and south polar cold spots and the long‐term effects of the 2001 global dust storm

In the 1970s, Viking and Mariner observed areas in the polar regions of Mars with winter brightness temperatures below the expected kinetic temperatures for CO 2 ice sublimation. These areas have since been termed “cold spots” and have been identified as surface deposits of CO 2 atmospheric condensates and, occasionally, active CO 2 storms. Three Mars years of data from the Mars Global Surveyor Thermal Emission Spectrometer were used to observe autumn and winter cold spot activity. In this study, cold spots that occur near and on the southern perennial cap were compared to those found near or on the northern perennial cap. On the southern perennial cap, cold spots associated with topographic features (induced by orographic lifting) were less common than cold spots independent of topography, similar to the north. However, the cold spots in the south lasted longer than those observed in the north. There is also evidence that cold spot formation in the south was affected by the global dust storm of 2001, even though the dust storm occurred during the southern spring and summer seasons. Prior to the dust storm, the amount of overall cold spot activity closer to the perennial cap increased and the average CO 2 grain size for most of the cold spots increased as well. Following the dust storm, the majority of cold spots in the south increased in size and duration but they did not form north of 62°S latitude, whereas, in other years, cold spots formed as far north as 48°S.

Journal of Geophysical Research E: Planets

Mars: Abundant recurring slope lineae (RSL) following the planet-encircling dust event (PEDE) of 2018

Recurring slope lineae (RSL) are dark linear markings on Mars that regrow annually and likely originate from the flow of either liquid water or granular material. Following the great dust storm (or planet-encircling dust event, PEDE) of Mars year (MY) 34, Mars Reconnaissance Orbiter/High Resolution Imaging Science Experiment has seen many more candidate RSL than in typical Mars years. They have been imaged at more than 285 unique locations from August 2018 (when the atmosphere was clearing as the PEDE decayed) to August 2019, about half (157) of which are locations where RSL have not been documented previously. In MY34, 150 active RSL sites were identified in the southern middle latitudes (SML, -60° to -30°), whereas an average of 36 active sites were observed in each previous year (MY28–33). Post-PEDE RSL are also present during southern summer over a wider range of latitude, slope aspect, and L s (areocentric longitude of the sun) than in prior years. These RSL sites usually show evidence for recent dust deposition: obscuration of relatively dark areas, an overall brighter and redder surface than in prior years, and dust devil tracks, which indicate dust lifting by several mechanisms. We speculate that dust-lifting processes may initiate and sustain RSL activity. The RSL may form from flows of dust (perhaps clumped) and/or sand that is destabilized by dust movement or directly mobilized by dust devils. If this is the case, then the otherwise puzzling recurrence and year-to-year variability of RSL activity can be at least partly explained. The dust replenishment varies from year to year, which could explain interannual variations in RSL activity.

Journal of Geophysical Research: Planets

Variation of energy and carbon fluxes from a restored temperate freshwater wetland and implications for carbon market verification protocols

Temperate freshwater wetlands are among the most productive terrestrial ecosystems, stimulating interest in using restored wetlands as biological carbon sequestration projects for greenhouse gas reduction programs. In this study, we used the eddy covariance technique to measure surface energy carbon fluxes from a constructed, impounded freshwater wetland during two annual periods that were 8&thinsp;years apart: 2002&ndash;2003 and 2010&ndash;2011. During 2010&ndash;2011, we measured methane (CH 4 ) fluxes to quantify the annual atmospheric carbon mass balance and its concomitant influence on global warming potential (GWP). Peak growing season fluxes of latent heat and carbon dioxide (CO 2 ) were greater in 2002&ndash;2003 compared to 2010&ndash;2011. In 2002, the daily net ecosystem exchange reached as low as &minus;10.6&thinsp;g&thinsp;C&thinsp;m &minus;2 &thinsp;d &minus;1 , which was greater than 3 times the magnitude observed in 2010 (&minus;2.9&thinsp;g&thinsp;C&thinsp;m &minus;2 &thinsp;d &minus;1 ). CH 4 fluxes during 2010&ndash;2011 were positive throughout the year and followed a strong seasonal pattern, ranging from 38.1&thinsp;mg&thinsp;C&thinsp;m &minus;2 &thinsp;d &minus;1 in the winter to 375.9&thinsp;mg&thinsp;C&thinsp;m &minus;2 &thinsp;d &minus;1 during the summer. The results of this study suggest that the wetland had reduced gross ecosystem productivity in 2010&ndash;2011, likely due to the increase in dead plant biomass (standing litter) that inhibited the generation of new vegetation growth. In 2010&ndash;2011, there was a net positive GWP (675.3&thinsp;g&thinsp;C&thinsp;m &minus;2 &thinsp;yr &minus;1 ), and when these values are evaluated as a sustained flux, the wetland will not reach radiative balance even after 500&thinsp;years.

California

The new GFDL global atmosphere and land model AM2-LM2: Evaluation with prescribed SST simulations

The configuration and performance of a new global atmosphere and land model for climate research developed at the Geophysical Fluid Dynamics Laboratory (GFDL) are presented. The atmosphere model, known as AM2, includes a new gridpoint dynamical core, a prognostic cloud scheme, and a multispecies aerosol climatology, as well as components from previous models used at GFDL. The land model, known as LM2, includes soil sensible and latent heat storage, groundwater storage, and stomatal resistance. The performance of the coupled model AM2-LM2 is evaluated with a series of prescribed sea surface temperature (SST) simulations. Particular focus is given to the model's climatology and the characteristics of interannual variability related to El Nin??o-Southern Oscillation (ENSO). One AM2-LM2 integration was perfor med according to the prescriptions of the second Atmospheric Model Intercomparison Project (AMIP II) and data were submitted to the Program for Climate Model Diagnosis and Intercomparison (PCMDI). Particular strengths of AM2-LM2, as judged by comparison to other models participating in AMIP II, include its circulation and distributions of precipitation. Prominent problems of AM2-LM2 include a cold bias to surface and tropospheric temperatures, weak tropical cyclone activity, and weak tropical intraseasonal activity associated with the Madden-Julian oscillation. An ensemble of 10 AM2-LM 2 integrations with observed SSTs for the second half of the twentieth century permits a statistically reliable assessment of the model's response to ENSO. In general, AM2-LM2 produces a realistic simulation of the anomalies in tropical precipitation and extratropical circulation that are associated with ENSO. ?? 2004 American Meteorological Society.

Journal of Climate

Hydrodynamics and sediment mobility processes over a degraded senile coral reef

Coral reefs can influence hydrodynamics and morphodynamics by dissipating and refracting incident wave energy, modifying circulation patterns, and altering sediment transport pathways. In this study, the sediment and hydrodynamic response of a senile (dead) barrier reef (Crocker Reef, located in the upper portion of the Florida Reef Tract) to storms and quiescent conditions was evaluated using field observations and the Coupled Ocean-Atmosphere-Wave-Sediment Transport model. Waves, circulation, and resultant sediment mobility were modeled across different reef zones. Sediment mobility during quiescent periods and the passage of far-field storms are driven by nonbreaking waves and, to a lesser degree, regional circulation. Spatial variability in these processes produces the present-day distribution of sediment grain size at Crocker Reef, wherein finer-grain material along a shallow central ridge is frequently mobilized (43% to 62% of the time), winnowed away, and deposited along the lower-energy flanks and in the fore reef where sand mobility occurs less frequently (32% to 43% and 1% to 22% of the time, respectively). Analysis of wave conditions for the period of 2006–2014 supports that wave heights rarely exceed the threshold for breaking (0.1% and 0.3% at the reef crest and at the reef flat, respectively), predominantly during the passage of tropical storms. There is a shift to a wave-breaking regime during near-field storms, creating the potential for mobilization of larger material and enhanced reef degradation. Sediment mobility can be enhanced due to wave skewness or the generation of free infragravity waves during periods of depth-induced wave breaking.

Florida

Identification of unrecognized tundra fire events on the north slope of Alaska

Characteristics of the natural fire regime are poorly resolved in the Arctic, even though fire may play an important role cycling carbon stored in tundra vegetation and soils to the atmosphere. In the course of studying vegetation and permafrost-terrain characteristics along a chronosequence of tundra burn sites from AD 1977, 1993, and 2007 on the North Slope of Alaska, we discovered two large, previously unrecognized tundra fires. The Meade River fire burned an estimated 500 km 2 and the Ketik River fire burned an estimated 1200 km 2 . Based on radiocarbon dating of charred twigs, analysis of historic aerial photography, and regional climate proxy data, these fires likely occurred between AD 1880 and 1920. Together, these events double the estimated burn area on the North Slope of Alaska over the last ~100 to 130 years. Assessment of vegetation succession along the century-scale chronosequence of tundra fire disturbances demonstrates for the first time on the North Slope of Alaska that tundra fires can facilitate the invasion of tundra by shrubs. Degradation of ice-rich permafrost was also evident at the fire sites and likely aided in the presumed changes of the tundra vegetation postfire. Other previously unrecognized tundra fire events likely exist in Alaska and other Arctic regions and identification of these sites is important for better understanding disturbance regimes and carbon cycling in Arctic tundra.

Alaska

Physical response of a back-barrier estuary to a post-tropical cyclone

This paper presents a modeling investigation of the hydrodynamic and sediment transport response of Chincoteague Bay (VA/MD, USA) to Hurricane Sandy using the Coupled Ocean-Atmosphere-Wave-Sediment-Transport (COAWST) modeling system. Several simulation scenarios with different combinations of remote and local forces were conducted to identify the dominant physical processes. While 80% of the water level increase in the bay was due to coastal sea level at the peak of the storm, a rich spatial and temporal variability in water surface slope was induced by local winds and waves. Local wind increased vertical mixing, horizontal exchanges, and flushing through the inlets. Remote waves (swell) enhanced southward flow through wave setup gradients between the inlets, and increased locally generated wave heights. Locally generated waves had a negligible effect on water level but reduced the residual flow up to 70% due to enhanced apparent roughness and breaking-induced forces. Locally generated waves dominated bed shear stress and sediment resuspension in the bay. Sediment transport patterns mirrored the interior coastline shape and generated deposition on inundated areas. The bay served as a source of fine sediment to the inner shelf, and the ocean-facing barrier island accumulated sand from landward-directed overwash. Despite the intensity of the storm forcing, the bathymetric changes in the bay were on the order of centimeters. This work demonstrates the spectrum of responses to storm forcing, and highlights the importance of local and remote processes on back-barrier estuarine function.

Chincoteague Bay

Links between atmosphere, ocean, and cryosphere from two decades of microseism observations on the Antarctic Peninsula

The lack of landmasses, climatological low pressure, and strong circumpolar westerly winds between the latitudes of 50°S to 65°S produce exceptional storm‐driven wave conditions in the Southern Ocean. This combination makes the Antarctic Peninsula one of Earth's most notable regions of high‐amplitude wave activity and thus, ocean‐swell‐driven microseism noise in both the primary (direct wave‐coastal region interactions) and secondary (direct ocean floor forcing due to interacting wave trains) period bands. Microseism observations are examined across 23 years (1993–2015) from Palmer Station (PMSA), on the west coast of the Antarctic Peninsula, and from East Falkland Island (EFI). These records provide a spatially integrative measure of both Southern Ocean wave amplitudes and the interactions between ocean waves and the solid Earth in the presence of sea ice, which can reduce wave coupling with the continental shelf. We utilize a spatiotemporal correlation‐based approach to illuminate how the distribution of sea ice influences seasonal microseism power. We characterize primary and secondary microseism power due to variations in sea ice and find that primary microseism energy is both more sensitive to sea ice and more capable of propagating across ocean basins than secondary microseism energy. During positive phases of the Southern Annular Mode, sea ice is reduced in the Bellingshausen Sea and overall storm activity in the Drake Passage increases, thus strongly increasing microseism power levels.

Antarctic Penninsula

Representing the effects of alpine grassland vegetation cover on the simulation of soil thermal dynamics by ecosystem models applied to the Qinghai-Tibetan Plateau

Soil surface temperature is a critical boundary condition for the simulation of soil temperature by environmental models. It is influenced by atmospheric and soil conditions and by vegetation cover. In sophisticated land surface models, it is simulated iteratively by solving surface energy budget equations. In ecosystem, permafrost, and hydrology models, the consideration of soil surface temperature is generally simple. In this study, we developed a methodology for representing the effects of vegetation cover and atmospheric factors on the estimation of soil surface temperature for alpine grassland ecosystems on the Qinghai-Tibetan Plateau. Our approach integrated measurements from meteorological stations with simulations from a sophisticated land surface model to develop an equation set for estimating soil surface temperature. After implementing this equation set into an ecosystem model and evaluating the performance of the ecosystem model in simulating soil temperature at different depths in the soil profile, we applied the model to simulate interactions among vegetation cover, freeze-thaw cycles, and soil erosion to demonstrate potential applications made possible through the implementation of the methodology developed in this study. Results showed that (1) to properly estimate daily soil surface temperature, algorithms should use air temperature, downward solar radiation, and vegetation cover as independent variables; (2) the equation set developed in this study performed better than soil surface temperature algorithms used in other models; and (3) the ecosystem model performed well in simulating soil temperature throughout the soil profile using the equation set developed in this study. Our application of the model indicates that the representation in ecosystem models of the effects of vegetation cover on the simulation of soil thermal dynamics has the potential to substantially improve our understanding of the vulnerability of alpine grassland ecosystems to changes in climate and grazing regimes.

Qinghai;Tibet

Evolving environmental and geometric controls on Columbia Glacier’s continued retreat

Geometry strongly controls the dynamic behavior of marine‐terminating (tidewater) glaciers, significantly influencing advance and retreat cycles independent of climate. Yet the recent, nearly ubiquitous retreat of tidewater glaciers suggests that changes in atmospheric and oceanic forcing may also drive dynamic change. To isolate the influence of geometry on tidewater glacier dynamics, we analyzed detailed observational time series from 2012 to 2016 for two tidewater glaciers with shared dynamic histories and environmental forcing: Columbia Glacier and its former tributary (Post Glacier) in southcentral Alaska. We find that although terminus retreat has driven decadal‐scale changes in dynamics of the Columbia‐Post system, environmental factors contribute to short‐term (i.e., seasonal) dynamic variability. In particular, analysis of force balance time series indicates that observed variations in speed result from seasonal changes to the subglacial hydrologic system and associated changes in basal drag. Variations in terminus position only drive noticeable speed change when the terminus retreats from regions of relatively high basal drag. In agreement with long‐term analyses of Columbia Glacier, we find that terminus geometry can perturb the timing of seasonal ice flow patterns. Specifically, our data support the idea that retreat of a glacier terminus into deeper water is accompanied by a shift in the primary control on frontal ablation. Although our analysis focuses on two Alaskan glaciers, our data suggest that changes in the relative importance of surface meltwater and buoyancy effects on submarine melting and/or calving may manifest as a shift in terminus change seasonality and offer a mechanism to identify frontal ablation controls.

Alaska

Drivers of methane flux differ between lakes and reservoirs, complicating global upscaling efforts

Methane is an important greenhouse gas with growing atmospheric concentrations. Freshwater lakes and reservoirs contribute substantially to atmospheric methane concentrations, but the magnitude of this contribution is poorly constrained. Uncertainty stems partially from whether the sites currently sampled represent the global population as well as incomplete knowledge of which environmental variables predict methane flux. Thus, determining the main drivers of methane flux across diverse waterbody types will inform more accurate upscaling approaches. Here we use a new database of total, diffusive, and ebullitive areal methane emissions from 313 lakes and reservoirs (ranging in surface area from 6 m 2 to 5,400 km 2 ) to identify the best predictors of methane emission. We found that the best predictors of methane emission differed by waterbody type (lakes vs. reservoirs), and that ecosystem morphometric variables (e.g., surface area and maximum depth) were more important predictors in lakes whereas metrics of autochthonous production (e.g., chlorophyll a ) were more important in reservoirs. We also found that productivity strongly predicted methane ebullition, whereas ecosystem morphometry and waterbody type were more important predictors of diffusive methane flux. Finally, we identify several knowledge gaps that limit upscaling efforts. First, we need more methane emission measurements in small reservoirs, large lakes, and both natural and artificial ponds. Additionally, more accurate upscaling efforts require improved global information about waterbody surface area, waterbody type (lake vs. reservoir), ice phenology, and the distribution of productivity‐related predictor variables such as total phosphorus, DOC, and chlorophyll a .

Journal of Geophysical Research-Biogeosciences