Geology Reports⌕ Search

SEARCH · Geology Reports

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 1,027 records · Page 57Linked to original sources
USGS2022

Insights on multistage rock avalanche behavior from runout modeling constrained by seismic inversions

Inversion of low-frequency regional seismic records to solve for a time series of bulk forces exerted on the earth by a landslide (a force-time function) is increasingly being used to infer volumes and dynamics of large, highly energetic landslides, such as rock avalanches and flowslides, and to provide calibration information on event dynamics and volumes for numerical landslide runout models. Much of the work to date using landslide runout modeling constrained by seismic data has focused on using single-phase models with frictional or velocity-weakening rheologies. Awareness of multistage landslide initiations is increasing, with discrete failures separated in time contributing to the final impact of an event. Our work utilizes a method for incorporating seismic data as a calibration constraint for landslide runout models, considering variable rheologies and different initiation conditions. This study presents a systematic examination of multiple rheologies and initiation conditions, and shows how these factors affect the force-time function derived from the landslide runout model. Our work confirms that, while rheology and fragmenting or initially coherent initiations affect the force-time function, multiple collapses separated by tens of seconds have the greatest impact on the shape and amplitude. We apply this method to the analysis of three real rock avalanches to better constrain plausible initiation conditions and rheology parameters using both seismic and field data. This study provides insights on how assumptions about the initiation dynamics of the source zone and the runout model definition can aid in the interpretation of seismic inversions for multistage rock avalanches.

Journal of Geophysical Research: Solid Earth↗
USGS2021

Cohesive sediment modeling in a shallow estuary: Model and environmental implications of sediment parameter variation

Numerical models of sediment transport in estuarine systems rely on parameter values that are often poorly constrained and can vary on timescales relevant to model processes. The selection of parameter values can affect the accuracy of model predictions, while environmental variation of these parameters can impact the temporal and spatial ranges of sediment fluxes, erosion, and deposition in the real world. We implemented a numerical model of San Pablo Bay, an embayment within San Francisco Bay, California, for November–December 2014, and compared model outputs to observations of water level, velocity, wave parameters, salinity, and suspended sediment concentration (SSC) in the shallow regions. Idealized model runs show that wind timing relative to the phase of the tides is the strongest control on sediment fluxes and bed erosion. We varied sediment erodibility in the outflow of the Petaluma River; while this causes erosion and deposition to vary strongly through the shallows system, total export from the shallows does not change. Model runs with realistic winds show that wind likely resuspends faster settling particles or allows for more particle flocculation; particle settling velocity controls system-wide sediment accumulation. At the margins of the system, the magnitude of SSC is closely tied to wind direction when winds occur during flood tide, but sediment deposition is less connected: Both bed evolution and SSC need to be considered in the prediction of marsh fate. Spatial patterns of light attenuation due to SSC is strongly tied to assumed settling velocity.

California↗
USGS2021

Distinguishing between regression model fits to global mean sea level reconstructions

Global mean sea level (GMSL) has been rising since the last century, posing a serious challenge for the coastal areas. A variety of regression models have been utilized for determining GMSL rise over the past one hundred years, resulting in a large spread of sea level rise rates and multidecadal variations. In this study, we develop a new nonparametric noise model that is data-dependent and considers overfitting due to regression. The noise model is used to determine whether one regression model has significantly better skill than others over the period 1900–2010. The choices of background noise and GMSL reconstruction influence whether two sea level models can be statistically distinguished. With our new nonparametric noise spectra, the differences of model skills in explaining sea level variance are significant only in 34% of model comparisons. However, stepwise trends with three inflection points are significantly more skillful than the linear, quadratic, or exponential trend for most GMSL reconstructions, suggesting the importance of multidecadal variability of sea level rise in the twentieth century. Nevertheless, stepwise trend models cannot be distinguished from models with a long-term harmonic oscillation, indicating that the shape of multidecadal variability is not conclusive. The multidecadal variability is also significant in the steric and barystatic sea level contributions and is related to both natural and anthropogenic forcings. GMSL predictions based on regression fits in the twentieth century underestimate the sea level rise rate over the period 2011–2020 because the sea level acceleration in the recent decade (2011–2020) is not well represented.

Journal of Geophysical Research- Oceans↗
USGS2021

Active Mars: A dynamic world

Mars exhibits diverse surface changes at all latitudes and all seasons. Active processes include impact cratering, aeolian sand and dust transport, a variety of slope processes, changes in polar ices, and diverse effects of seasonal CO 2 frost. The extent of surface change has been surprising and indicates that the present climate is capable of reshaping the surface. Activity has important implications for the Amazonian history of Mars: understanding processes is a necessary step before we can understand their implications and variations over time.

Journal of Geophysical Research: Planets↗
USGS2022

Orbital and in-situ investigation of periodic bedrock ridges in Glen Torridon, Gale Crater, Mars

Wind has been the dominant agent of landscape modification on Mars for the past ~3 billion years. Among the diversity of features formed by aeolian abrasion on the surface of Mars are periodic bedrock ridges (PBRs), a relatively recently recognized class of erosional bedforms on Mars for which Earth analogues are rare. Gale crater, the field site for NASA’s Mars Science Laboratory Curiosity rover since it landed there in 2012, contains a diverse and extensive record of aeolian deposition and erosion. This study focuses on a series of periodic, linear bedrock ridges that occur within the Fe/Mg-smectite clay-bearing Glen Torridon region of Aeolis Mons (informally Mount Sharp). During Curiosity’s exploration of the Glen Torridon region between sols ~2300-3080, the rover drove through this field of ridges, providing the first opportunity for the in situ observation of these enigmatic erosional features. This study characterizes the Glen Torridon ridges using orbiter and rover data to determine their morphology, spatial distribution, compositional and material properties, and association with other aeolian features in the area. Based on these observations, the Glen Torridon ridges are interpreted to be consistent with an origin as wind-eroded periodic bedrock ridges carved during the most recent exhumation of Mount Sharp into the present-day mound. Although there is evidence for multidirectional winds in the Glen Torridon region based on the orientation of modern ripples, megaripples, TARs and other bedrock indicators, the consistent orientation of the Glen Torridon ridges, coupled with morphologic asymmetries within the ridges, support formation and elongation of the Glen Torridon PBRs forms parallel to a net regional northerly wind direction in and around Gale crater.

Journal of Geophysical Research E: Planets↗
USGS2022

Mars Science Laboratory CheMin data from the Glen Torridon region and the significance of lake-groundwater interactions in interpreting mineralogy and sedimentary history

The Glen Torridon (GT) region is positioned in terrains with strong clay mineral signatures, as inferred from orbital spectroscopy. The GT campaign confirmed orbital distinctions with in situ measurements by the Mars Science Laboratory rover, Curiosity , and the CheMin X-ray diffraction instrument with of some of the highest clay mineral abundances to date. Additionally, GT is unique because of distinct phase identifications for the first time by CheMin, including: (i) Fe-carbonates, and (ii) a novel peak in the XRD patterns of some GT samples, with an interplanar spacing of at 9.2 Å. Fe-carbonates have been previously suggested from other instruments onboard, but this is the first definitive reporting by CheMin of multiple samples with Fe-carbonate. This new phase has never been observed in Gale crater with the CheMin instrument and may be a new mineral for Mars, but discrete identification still remains enigmatic because no single phase on Earth is able to account the mineralogical, geochemical, and sedimentological constraints in the GT region. Here, we modeled XRD profiles and propose an interstratified clay mineral, specifically greenalite-minnesotaite (G-M), as a reasonable candidate. The coexistence of Fe-carbonate and Fe-rich clay minerals in the GT samples, supports a conceptual model of a lacustrine groundwater mixing environment in the ancient Gale crater lake. Groundwater interaction with percolating lake waters in the sediments is common in terrestrial lacustrine settings, and the diffusion of two distinct water bodies within the subsurface can create a geochemical gradient and unique mineral front in the sediments. Ultimately, the proximity to this mixing zone controlled the secondary minerals preserved in the Jura, Knockfarril Hill, and Glasgow members of the Mt. Sharp group exposed in GT.

Journal of Geophysical Research E: Planets↗
USGS2022

Bedrock geochemistry and alteration history of the clay-bearing Glen Torridon region of Gale crater, Mars

Glen Torridon is a topographic trough located on the slope of Aeolis Mons, Gale crater, Mars. It corresponds to what was previously referred to as the “clay-bearing unit”, due to the relatively strong spectral signatures of clay minerals (mainly ferric smectites) detected from orbit. Starting in January 2019, the Curiosity rover explored Glen Torridon for more than 700 sols (Martian days). The objectives of this campaign included acquiring a detailed understanding of the geologic context in which the clay minerals were formed and determining the intensity of aqueous alteration experienced by the sediments. Here, we present the major-element geochemistry of the bedrock as analyzed by the ChemCam instrument. Our results reveal that the two main types of bedrock exposures identified in the lower part of Glen Torridon are associated with distinct chemical compositions (K-rich and Mg-rich), for which we are able to propose mineralogical interpretations. Moreover, the topmost stratigraphic member exposed in the region displays a stronger diagenetic overprint, especially at two locations close to the unconformable contact with the overlying Stimson formation, where the bedrock composition significantly deviates from the rest of Glen Torridon. Overall, the values of the Chemical Index of Alteration determined with ChemCam are elevated by Martian standards, suggesting the formation of clay minerals through open system weathering. However, there is no indication that the alteration was stronger than in some terrains previously visited by Curiosity , which in turn implies that the enhanced orbital signatures are mostly controlled by non-compositional factors.

Journal of Geophysical Research E: Planets↗
USGS2021

Postwildfire soil‐hydraulic recovery and the persistence of debris flow hazards

Deadly and destructive debris flows often follow wildfire, but understanding of changes in the hazard potential with time since fire is poor. We develop a simulation‐based framework to quantify changes in the hydrologic triggering conditions for debris flows as postwildfire infiltration properties evolve through time. Our approach produces time‐varying rainfall intensity‐duration thresholds for runoff‐ and infiltration‐generated debris flows with physics‐based hydrologic simulations that are parameterized with widely available hydroclimatic, vegetation reflectance, and soil texture data. When we apply our thresholding protocol to a test case in the San Gabriel Mountains (California, USA), the results are consistent with existing regional empirical thresholds and rainstorms that caused runoff‐ and infiltration‐generated debris flows soon after and three years following a wildfire, respectively. We find that the hydrologic triggering mechanisms for the two observed debris flow types are coupled with the effects of fire on the soil saturated hydraulic conductivity. Specifically, the rainfall intensity needed to generate debris flows via runoff increases with time following wildfire while the rainfall duration needed to produce debris flows via subsurface pore‐water pressures decreases. We also find that variations in soil moisture, rainfall climatology, median grain size, and root reinforcement could impact the median annual probability of postwildfire debris flows. We conclude that a simulation‐based method for calculating rainfall thresholds is a tractable approach to improve situational awareness of debris flow hazard in the years following wildfire. Further development of our framework will be important to quantify postwildfire hazard levels in variable climates, vegetation types, and fire regimes.

Journal of Geophysical Research: Earth Surface↗
USGS2021

Multi-model comparison of computed debris flow runout for the 9 January 2018 Montecito, California post-wildfire event

Hazard assessment for post-wildfire debris flows, which are common in the steep terrain of the western United States, has focused on the susceptibility of upstream basins to generate debris flows. However, reducing public exposure to this hazard also requires an assessment of hazards in downstream areas that might be inundated during debris flow runout. Debris flow runout models are widely available, but their application to hazard assessment for post-wildfire debris flows has not been extensively tested. Necessary inputs to these models include the total volume of the mobilized flow, flow properties (either inherent material properties or calibration coefficients), and site topography. Estimates of volume are possible in post-event (“back calculation”) studies, yet before an event, volume is an uncertain quantity. We simulated debris flow runout for the well-constrained 9 January 2018 Montecito event using three models (RAMMS, FLO2D, and D-Claw) to determine the relative importance of volume and flow properties. We broke the impacted area into three domains, and for each model-domain combination, we performed a numerical sampling study in which volume and flow properties varied within a wide, but plausible range. We assessed model performance based on inundation patterns and peak flow depths. We found all models could simulate the event with comparable results. Simulation performance was most sensitive to flow volume and less sensitive to flow properties. Our results emphasize the importance of reducing uncertainty in pre-event estimates of flow volume for hazard assessment.

California↗
USGS2022

Vadose zone thickness limits pore-fluid pressures and acceleration in a large, slow-moving landslide

The rate and timing of hydrologically forced landslides is a complex function of precipitation patterns, material properties, topography, and groundwater hydrology. In the simplest form, however, slopes fail when subsurface pore pressure grows large enough to exceed the Mohr-Coulomb failure criterion. The capacity for pore pressure rise in a landslide is determined in part by the thickness of the unsaturated zone above the water table, which itself is set by weathering patterns that should have predictable patterns across different lithologies. To investigate how this structure affects landslide behavior, we exploit a multi-year record of precipitation, pore pressure, and velocity from Oak Ridge earthflow, a slow-moving landslide set in Franciscan mélange, northern California, USA. In conjunction with electrical resistivity tomography and hydraulic conductivity measurements, these data show that Oak Ridge has a thin weathered profile that is comparable in thickness to other mélange landslides in California. We propose that due to the inherently thin vadose zone, mélange landscapes experience an unusually high water table that frequently brings them close to movement; however, the capacity to increase stress is limited by the small amount of dynamic storage available. Instead, excess pore pressure is shed via springs and saturation overland flow once the water table reaches the surface. Linkages between weathering patterns, hydrology, and deformation can explain behavior patterns exhibited by Franciscan mélange earthflows across a large precipitation gradient.

California↗
USGS2022

Multi-decadal simulation of marsh topography evolution under sea level rise and episodic sediment loads

Coastal marsh within Mediterranean climate zones is exposed to episodic watershed runoff and sediment loads that occur during storm events. Simulating future marsh accretion under sea level rise calls for attention to: (a) physical processes acting over the time scale of storm events and (b) biophysical processes acting over time scales longer than storm events. Using the upper Newport Bay in Southern California as a case study, we examine the influence of event-scale processes on simulated change in marsh topography by comparing: (a) a biophysical model that integrates with an annual time step and neglects event-scale processes (BP-Annual), (b) a physical model that resolves event-scale processes but neglects biophysical interactions (P-Event), and (c) a biophysical model that resolves event-scale physical processes and biophysical processes at annual and longer time scales (BP-Event). A calibrated BP-Event model shows that large (>20-year return period) episodic storm events are major drivers of marsh accretion, depositing up to 30 cm of sediment in one event. Greater deposition is predicted near fluvial sources and tidal channels and less on marshes further from fluvial sources and tidal channels. In contrast, the BP-Annual model poorly resolves spatial structure in marsh accretion as a consequence of neglecting event-scale processes. Furthermore, the P-Event model significantly overestimates marsh accretion as a consequence of neglecting marsh surface compaction driven by annual scale biophysical processes. Differences between BP-Event and BP-Annual models translate up to 20 cm per century in marsh surface elevation.

California↗
USGS2021

Storm-scale and seasonal dynamics of carbon export from a nested subarctic watershed underlain by permafrost

Subarctic catchments underlain by permafrost sequester a major stock of frozen organic carbon (C), which may be mobilized as the Arctic warms. Warming can impact C export from thawing soils by altering the depth and timing of runoff related to changing storm and fire regimes and altered soil thaw depths. We investigated C export in a first order headwater stream (West Twin Creek) and its receiving third order river (Beaver Creek) in interior Alaska using discrete sampling of dissolved organic and inorganic C (DOC and DIC) and 15-min collection of specific conductance (SC), fluorescent dissolved organic matter (fDOM) and water discharge (Q). Storm SC-Q relationships displayed negative slopes, indicating solute limitation and limited influence of seasonal soil thaw on storm runoff chemistry. Concurrently, fDOM-Q displayed positive slopes that decreased over the summer, indicating flushing of a limited fDOM pool. Baseflow DIC increased over the season concurrent with soil thaw, with higher DIC at the larger scale indicating greater influence of deeper, mineral-rich flow paths. Storm and seasonal trends were generally similar at both scales. The biggest difference was in fDOM, which displayed higher concentrations and slower depletion in the first order stream. Improved process understanding from this study can be used to better predict carbon export and cycling by stream networks as northern forests and arctic regions continue to warm.

Alaska↗
USGS2021

Event scale relationships of DOC and TDN fluxes in throughfall and stemflow diverge from stream exports in a forested catchment

Aquatic fluxes of carbon and nutrients link terrestrial and aquatic ecosystems. Within forests, storm events drive both the delivery of carbon and nitrogen to the forest floor and the export of these solutes from the land via streams. To increase understanding of the relationships between hydrologic event character and the relative fluxes of carbon and nitrogen in throughfall, stemflow and streams, we measured dissolved organic carbon (DOC) and total dissolved nitrogen (TDN) concentrations in each flow path for 23 events in a forested watershed in Vermont, USA. DOC and TDN concentrations increased with streamflow, indicating their export was limited by water transport of catchment stores. DOC and TDN concentrations in throughfall and stemflow decreased exponentially with increasing precipitation, suggesting that precipitation removed a portion of available sources from tree surfaces during the events. DOC and TDN fluxes were estimated for 76 events across a 2-year period. For most events, throughfall and stemflow fluxes greatly exceeded stream fluxes, but the imbalance narrowed for larger storms (>30 mm). The largest 10 stream events exported 40% of all stream event DOC whereas those same 10 events contributed 14% of all throughfall export. Approximately 2–5 times more DOC and TDN was exported from trees during rain events than left the catchment via streams annually. The diverging influence of event size on tree versus stream fluxes has important implications for forested ecosystems as hydrological events increase in intensity and frequency due to climate change.

Vermont↗
USGS2021

Aquatic-terrestrial linkages control metabolism and carbon dynamics in a mid-sized, urban stream influenced by snowmelt

Freshwater streams can exchange nutrients and carbon with the surrounding terrestrial environment through various mechanisms including physical erosion, flooding, leaf drop, and snowmelt. These aquatic-terrestrial interactions are crucial in carbon mobilization, transformation, ecosystem productivity, and have important implications for the role of freshwater ecosystems in the global carbon budget. We utilized high-frequency oxygen, temperature, and carbon dioxide (CO 2 ) data to infer watershed connectivity in Boulder Creek, a mid-sized (1160 km 2 ) watershed located in Colorado, USA. Daily modeled gross primary production (GPP), ecosystem respiration (ER), net ecosystem production (NEP), and reaeration coefficients ( K 600 ) were paired with high-frequency, in-situ dissolved CO 2 data to characterize changes in metabolic regime and carbon flux on a stream influenced by seasonal snowmelt. GPP and ER were correlated ( ρ = −0.72, p ≪ 0.001) during the non-snowmelt period and NEP was frequently negative. Mean F CO2 during the non-snowmelt period was approximately 302 (±171) mmol C m −2 d −1 and was primarily supported by watershed CO 2 inputs. During snowmelt, GPP and ER were not significantly correlated ( ρ = −0.22, p = 0.05), and mean NEP was significantly more negative than during non-snowmelt. Watershed connectivity was higher during snowmelt, as evidenced by significantly higher F CO2 (843 ± 338 mmol C m −2 d −1 ) and greater allochthonous CO 2 inputs than during non-snowmelt periods, emphasizing the effects of seasonal differences in aquatic-terrestrial linkages in this stream. We suggest that our understanding of watershed carbon budgets is subject to temporal dynamics which control the degree of connectivity between terrestrial and aquatic ecosystems.

Colorado↗
USGS2021

Spatiotemporal methane emission from global reservoirs

Inland aquatic systems, such as reservoirs, contribute substantially to global methane (CH 4 ) emissions; yet are among the most uncertain components of the total CH 4 budget. Reservoirs have received recent attention as they may generate high CH 4 fluxes. Improved quantification of these CH 4 fluxes, particularly their spatiotemporal distribution, is key to realistically incorporating them in CH 4 modeling and budget studies. Here we report on a new global, gridded (0.25° lat × 0.25° lon) study of reservoir CH 4 emissions, accounting for new knowledge regarding reservoir areal extent and distribution, and spatiotemporal emission patterns influenced by diurnal variability, temperature-dependent seasonality, satellite-derived freeze-thaw dynamics, and eco-climatic zone. The results of this new data set comprise daily CH 4 emissions throughout the full annual cycle and show that reservoirs cover 297 × 10 3 km 2 globally and emit 10.1 Tg CH 4 yr −1 (1σ uncertainty range of 7.2–12.9 Tg CH 4 yr −1 ) from diffusive (1.2 Tg CH 4 yr −1 ) and ebullitive (8.9 Tg CH 4 yr −1 ) emission pathways. This analysis of reservoir CH 4 emission addresses multiple gaps and uncertainties in previous studies and represents an important contribution to studies of the global CH 4 budget. The new data sets and methodologies from this study provide a framework to better understand and model the current and future role of reservoirs in the global CH 4 budget and to guide efforts to mitigate reservoir-related CH 4 emissions.

Journal of Geophysical Research: Biogeosciences↗
USGS2021

Loss of ice cover, shifting phenology, and more extreme events in Northern Hemisphere lakes

Long-term lake ice phenological records from around the Northern Hemisphere provide unique sensitive indicators of climatic variations, even prior to the existence of physical meteorological measurement stations. Here, we updated ice phenology records for 60 lakes with time-series ranging from 107–204 years to provide the first re-assessment of Northern Hemispheric ice trends since 2004 by adding 15 additional years of ice phenology records and 40 lakes to our study. We found that, on average, ice-on was 11.0 days later, ice-off was 6.8 days earlier, and ice duration was 17.0 days shorter per century over the entire record for each lake. Trends in ice-on and ice duration were six times faster in the last 25-year period (1992–2016) than previous quarter centuries. More extreme events in recent decades, including late ice-on, early ice-off, shorter periods of ice cover, or no ice cover at all, contribute to the increasing rate of lake ice loss. Reductions in greenhouse gas emissions could limit increases in air temperature and abate losses in lake ice cover that would subsequently limit ecological, cultural, and socioeconomic consequences, such as increased evaporation rates, warmer water temperatures, degraded water quality, and the formation of toxic algal blooms.

Journal of Geophysical Research: Biogeosciences↗
USGS2021

Episodic nutrient addition affects water column nutrient processing rates in river-to-lake transitional zones

Storm-driven nutrient loading from tributaries can fuel eutrophication in nearshore and open water areas of lentic ecosystems. However, nutrient processing in river-to-lake transitional zones can substantially alter the amount and composition of nutrients transported to lakes from upstream surface waters. We measured the removal of nutrients and dissolved organic carbon (DOC) from the water column in the Fox rivermouth (Green Bay, Lake Michigan) to evaluate the response of rivermouth plankton to episodic nutrient enrichment. Light and dark water column incubations (8–12 hr) were conducted on four occasions from April through September to measure changes in dissolved nitrogen (N), phosphorus (P), and DOC concentrations in three locations along the Fox rivermouth. Two incubation experiments were conducted on consecutive days, (a) under ambient nutrient concentrations, and (b) under experimentally enriched N and P concentrations. Spatial and temporal variation was observed in nutrient uptake rates, but light incubations consistently had higher nutrient uptake rates than dark incubations. Nutrient enrichment increased total dissolved P and total dissolved N uptake and DOC release in light incubations, but only increased total dissolved P uptake in dark incubations. Moreover, nutrient uptake ratios (N:P) decreased from ambient to nutrient enriched conditions and indicated preferential P uptake by phytoplankton communities in light conditions. Our study substantiates that rivermouths can process nutrients bound for downstream ecosystems and demonstrates the potential of plankton communities to dynamically increase net uptake rates in response to episodic nutrient enrichment.

Wisconsin↗
USGS2021

Lagged wetland CH4 flux response in a historically wet year

While a stimulating effect of plant primary productivity on soil carbon dioxide (CO 2 ) emissions has been well documented, links between gross primary productivity (GPP) and wetland methane (CH 4 ) emissions are less well investigated. Determination of the influence of primary productivity on wetland CH 4 emissions (FCH 4 ) is complicated by confounding influences of water table level and temperature on CH 4 production, which also vary seasonally. Here, we evaluate the link between preceding GPP and subsequent FCH 4 at two fens in Wisconsin using eddy covariance flux towers, Lost Creek (US-Los) and Allequash Creek (US-ALQ). Both wetlands are mosaics of forested and shrub wetlands, with US-Los being larger in scale and having a more open canopy. Co-located sites with multi-year observations of flux, hydrology, and meteorology provide an opportunity to measure and compare lag effects on FCH 4 without interference due to differing climate. Daily average FCH 4 from US-Los reached a maximum of 47.7 ηmol CH 4 m −2 s −1 during the study period, while US-ALQ was more than double at 117.9 ηmol CH 4 m −2 s −1 . The lagged influence of GPP on temperature-normalized FCH 4 ( T air -FCH 4 ) was weaker and more delayed in a year with anomalously high precipitation than a following drier year at both sites. FCH 4 at US-ALQ was lower coincident with higher stream discharge in the wet year (2019), potentially due to soil gas flushing during high precipitation events and lower water temperatures. Better understanding of the lagged influence of GPP on FCH 4 due to this study has implications for climate modeling and more accurate carbon budgeting.

Journal of Geophysical Research: Biogeosciences↗