USGS Science⌕ Search

SEARCH · USGS Science

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,729 records · Page 96Linked to original sources
USGS2003

A delayed action oscillator shared by biennial, interannual, and decadal signals in the Pacific Basin

Biennial, interannual, and decadal signals in the Pacific basin are observed to share patterns and evolution in covarying sea surface temperature (SST), 18??C isotherm depth (Z18), zonal surface wind (ZSW), and wind stress curl (WSC) anomalies from 1955 to 1999. Each signal has warm SST anomalies propagating slowly eastward along the equator, generating westerly ZSW anomalies in their wake. These westerly ZSW anomalies produce cyclonic WSC anomalies off the equator which pump baroclinic Rossby waves in the western/central tropical North Pacific Ocean. These Rossby waves propagate westward, taking ???6, ???12, and ???36 months to reach the western boundary near ???7??N, ???12??N, and ???18??N on biennial, interannual, and decadal period scales, respectively. There, they reflect as equatorial coupled waves, propagating slowly eastward in covarying SST, Z18, and ZSW anomalies, taking ???6, ???12, and ???24 months to reach the central/eastern equatorial ocean. These equatorial coupled waves produce a delayed-negative feedback to the warm SST anomalies there. The decrease in Rossby wave phase speed with latitude, the increase in meridional scale of equatorial SST anomalies with period scale, and the associated increase in latitude of Rossby wave forcing are consistent with the delayed action oscillator (DAO) model used to explain El Nin??o. However, this is not true of the western-boundary reflection of Rossby waves into slow equatorial coupled waves. This requires modification of the extant DAO model. We construct a modified DAO model, demonstrating how the various mechanisms and the size and sources of their delays yield the resulting frequency of each signal.

Journal of Geophysical Research C: Oceans↗
USGS2003

Late Pleistocene to Holocene slip rates for the Gurvan Bulag thrust fault (Gobi-Altay, Mongolia) estimated with 10Be dates

We surveyed morphotectonic markers along the central part of the Gurvan Bulag thrust, a fault that ruptured with the Bogd fault during the Gobi-Altay earthquake (1957, M 8.3), to document climatic and tectonic processes along the fault for the late Pleistocene- Holocene period. The markers were dated using 10Be produced in situ. Two major periods of alluviation ended at 131 ?? 20 and 16 ?? 4.8 ka. These appear to be contemporaneous with global climatic changes at the terminations of marine isotope stages (MIS) 6 and 2. The vertical slip rates, determined from offset measurements and surfaces ages, are 0.14 ?? 0.03 mm/yr over the late Pleistocene-Holocene and between 0.44 ?? 0.11 and 1.05 ?? 0.25 mm/yr since the end of the late Pleistocene. The higher of these slip rates for the last ???16 kyr is consistent with paleoseismic investigations along the fault [Prentice et al., 2002], and suggests that, at the end of late Pleistocene, the fault evolved from quiescence to having recurrence intervals of 4.0 ?? 1.2 kyr for surface ruptures with ???4 m vertical offset (similar to that of 1957). The inferred recurrence interval is comparable to that of the Bogd fault (3.7 ?? 1.3 kyr) suggesting that the two faults may have ruptured together also earlier during the last ???16 kyr.

Journal of Geophysical Research B: Solid Earth↗
USGS2003

Apparent break in earthquake scaling due to path and site effects on deep borehole recordings

We reexamine the scaling of stress drop and apparent stress, rigidity times the ratio between seismically radiated energy to seismic moment, with earthquake size for a set of microearthquakes recorded in a deep borehole in Long Valley, California. In the first set of calculations, we assume a constant Q and solve for the corner frequency and seismic moment. In the second set of calculations, we model the spectral ratio of nearby events to determine the same quantities. We find that the spectral ratio technique, which can account for path and site effects or nonconstant Q, yields higher stress drops, particularly for the smaller events in the data set. The measurements determined from spectral ratios indicate no departure from constant stress drop scaling down to the smallest events in our data set (Mw 0.8). Our results indicate that propagation effects can contaminate measurements of source parameters even in the relatively clean recording environment of a deep borehole, just as they do at the Earth's surface. The scaling of source properties of microearthquakes made from deep borehole recordings may need to be reevaluated.

Journal of Geophysical Research B: Solid Earth↗
USGS2003

Tsunamis generated by subaerial mass flows

Tsunamis generated in lakes and reservoirs by subaerial mass flows pose distinctive problems for hazards assessment because the domain of interest is commonly the "near field," beyond the zone of complex splashing but close enough to the source that wave propagation effects are not predominant. Scaling analysis of the equations governing water wave propagation shows that near-field wave amplitude and wavelength should depend on certain measures of mass flow dynamics and volume. The scaling analysis motivates a successful collapse (in dimensionless space) of data from two distinct sets of experiments with solid block "wave makers." To first order, wave amplitude/water depth is a simple function of the ratio of dimensionless wave maker travel time to dimensionless wave maker volume per unit width. Wave amplitude data from previous laboratory investigations with both rigid and deformable wave makers follow the same trend in dimensionless parameter space as our own data. The characteristic wavelength/water depth for all our experiments is simply proportional to dimensionless wave maker travel time, which is itself given approximately by a simple function of wave maker length/water depth. Wave maker shape and rigidity do not otherwise influence wave features. Application of the amplitude scaling relation to several historical events yields "predicted" near-field wave amplitudes in reasonable agreement with measurements and observations. Together, the scaling relations for near-field amplitude, wavelength, and submerged travel time provide key inputs necessary for computational wave propagation and hazards assessment.

Journal of Geophysical Research B: Solid Earth↗
USGS2003

Evidence for a Battle Mountain-Eureka crustal fault zone, north-central Nevada, and its relation to Neoproterozoic-Early Paleozoic continental breakup

Combined evidence from gravity, radiogenic isotope, and magnetotelluric (MT) data indicates a crustal fault zone that coincides with the northwest-trending Battle Mountain-Eureka (BME) mineral trend in north-central Nevada, USA. The BME crustal fault zone likely originated during Neoproterozoic-Early Paleozoic rifting of the continent and had a large influence on subsequent tectonic events, such as emplacement of allochthons and episodic deformation, magmatism, and mineralization throughout the Phanerozoic. MT models show the fault zone is about 10 km wide, 130-km long, and extends from 1 to 5 km below the surface to deep crustal levels. Isotope data and gravity models imply the fault zone separates crust of fundamentally different character. Geophysical evidence for such a long-lived structure, likely inherited from continental breakup, defies conventional wisdom that structures this old have been destroyed by Cenozoic extensional processes. Moreover, the coincidence with the alignment of mineral deposits supports the assertion by many economic geologists that these alignments are indicators of buried regional structures.

Journal of Geophysical Research B: Solid Earth↗
USGS2003

Slip rate and earthquake recurrence along the central Septentrional fault, North American-Caribbean plate boundary, Dominican Republic

The Septentrional fault zone (SFZ) is the major North American-Caribbean, strike-slip, plate boundary fault at the longitude of eastern Hispaniola. The SFZ traverses the densely populated Cibao Valley of the Dominican Republic, forming a prominent scarp in alluvium. Our studies at four sites along the central SFZ are aimed at quantifying the late Quaternary behavior of this structure to better understand the seismic hazard it represents for the northeastern Caribbean. Our investigations of excavations at sites near Rio Cenovi show that the most recent ground-rupturing earthquake along this fault in the north central Dominican Republic occurred between A.D. 1040 and A.D. 1230, and involved a minimum of ???4 m of left-lateral slip and 2.3 m of normal dip slip at that site. Our studies of offset stream terraces at two locations, Rio Juan Lopez and Rio Licey, provide late Holocene slip rate estimates of 6-9 mm/yr and a maximum of 11-12 mm/yr, respectively, across the Septentrional fault. Combining these results gives a best estimate of 6-12 mm/yr for the slip rate across the SFZ. Three excavations, two near Tenares and one at the Rio Licey site, yielded evidence for the occurrence of earlier prehistoric earthquakes. Dates of strata associated with the penultimate event suggest that it occurred post-A.D. 30, giving a recurrence interval of 800-1200 years. These studies indicate that the SFZ has likely accumulated elastic strain sufficient to generate a major earthquake during the more than 800 years since it last slipped and should be considered likely to produce a destructive future earthquake.

Journal of Geophysical Research B: Solid Earth↗
USGS2003

Reservoir depletion at The Geysers geothermal area, California, shown by four-dimensional seismic tomography

Intensive geothermal exploitation at The Geysers geothermal area, California, induces myriads of small-magnitude earthquakes that are monitored by a dense, permanent, local seismometer network. Using this network, tomographic inversions were performed for the three-dimensional Vp and Vp/Vs structure of the reservoir for April 1991, February 1993, December 1994, October 1996, and August 1998. The extensive low-Vp/Vs anomaly that occupies the reservoir grew in strength from a maximum of 9% to a maximum of 13.4% during the 7-year study period. This is attributed to depletion of pore liquid water in the reservoir and replacement with steam. This decreases Vp by increasing compressibility, and increases Vs because of reduction in pore pressure and the drying of argillaceous minerals, e.g., illite, which increase the shear modulus. These effects serendipitously combine to lower Vp/Vs, resulting in a strong overall effect that provides a convenient tool for monitoring reservoir depletion. Variations in the Vp and Vs fields indicate that water depletion is the dominant process in the central part of the exploited reservoir, and pressure reduction and mineral drying in the northwest and southeast parts of the reservoir. The rate at which the Vp/Vs anomaly grew in strength in the period 1991-1998 suggests most of the original anomaly was caused by exploitation. Continuous monitoring of Vp, Vs, and Vp/Vs is an effective geothermal reservoir depletion monitoring tool and can potentially provide information about depletion in parts of the reservoir that have not been drilled.

Journal of Geophysical Research B: Solid Earth↗
USGS2003

Crustal magnetization and accretion at the Southwest Indian Ridge near the Atlantis II fracture zone, 0-25 Ma

We analyze geophysical data that extend from 0 to 25-Myr-old seafloor on both flanks of the Southwest Indian Ridge (SWIR). Lineated marine magnetic anomalies are consistent and identifiable within the study area, even over seafloor lacking a basaltic upper crust. The full spreading rate of 14 km/Myr has remained nearly constant since at least 20 Ma, but crustal accretion has been highly asymmetric, with half rates of 8.5 and 5.5 km/Myr on the Antarctic and African flanks, respectively. This asymmetry may be unique to a ???400 km wide corridor between large-offset fracture zones of the SWIR. In contrast to the Mid-Atlantic Ridge, crustal magnetization amplitudes correlate directly with seafloor topography along the present-day rift valleys. This pattern appears to be primarily a function of along-axis variations in crustal thickness, rather than magnetic mineralogy. Off-axis, magnetization amplitudes at paleo-segment ends are more positive than at paleo-segment midpoints, suggesting the presence of an induced component of magnetization within the lower crust or serpentinized upper mantle. Alteration of the magnetic source layer at paleo-segment midpoints reduces magnetization amplitudes by 70-80% within 20 Myr of accretion. Magnetic and Ocean Drilling Program (ODP) Hole 735B data suggest that the lower crust cooled quickly enough to lock in a primary thermoremanent magnetization that is in phase with that of the overlying upper crust. Thus magnetic polarity boundaries within the intrusive lower crust may be steeper than envisioned in prior models of ocean crustal magnetization. As the crust ages, the lower crust becomes increasingly important in preserving marine magnetic stripes.

Journal of Geophysical Research B: Solid Earth↗
USGS2003

Three-dimensional crustal structure of Long Valley caldera, California, and evidence for the migration of CO2 under Mammoth Mountain

A temporary network of 69 three-component seismic stations captured a major seismic sequence in Long Valley caldera in 1997. We performed a tomographic inversion for crustal structure beneath a 28 km ?? 16 km area encompassing part of the resurgent dome, the south moat, and Mammoth Mountain. Resolution of crustal structure beneath the center of the study volume was good down to ???3 km below sea level (???5 km below the surface). Relatively high wave speeds are associated with the Bishop Tuff and lower wave speeds characterize debris in the surrounding moat. A low-Vp/Vs anomaly extending from near the surface to ???1 km below sea level beneath Mammoth Mountain may represent a CO2 reservoir that is supplying CO2-rich springs, venting at the surface, and killing trees. We investigated temporal variations in structure beneath Mammoth Mountain by differencing our results with tomographic images obtained using data from 1989/1990. Significant changes in both Vp and Vs were consistent with the migration of CO2 into the upper 2 km or so beneath Mammoth Mountain and its depletion in peripheral volumes that correlate with surface venting areas. Repeat tomography is capable of detecting the migration of gas beneath active silicic volcanoes and may thus provide a useful volcano monitoring tool.

Journal of Geophysical Research B: Solid Earth↗
USGS2014

An empirical approach to modeling methylmercury concentrations in an Adirondack stream watershed

Inverse empirical models can inform and improve more complex process-based models by quantifying the principal factors that control water quality variation. Here we developed a multiple regression model that explains 81% of the variation in filtered methylmercury (FMeHg) concentrations in Fishing Brook, a fourth-order stream in the Adirondack Mountains, New York, a known “hot spot” of Hg bioaccumulation. This model builds on previous observations that wetland-dominated riparian areas are the principal source of MeHg to this stream and were based on 43 samples collected during a 33 month period in 2007–2009. Explanatory variables include those that represent the effects of water temperature, streamflow, and modeled riparian water table depth on seasonal and annual patterns of FMeHg concentrations. An additional variable represents the effects of an upstream pond on decreasing FMeHg concentrations. Model results suggest that temperature-driven effects on net Hg methylation rates are the principal control on annual FMeHg concentration patterns. Additionally, streamflow dilutes FMeHg concentrations during the cold dormant season. The model further indicates that depth and persistence of the riparian water table as simulated by TOPMODEL are dominant controls on FMeHg concentration patterns during the warm growing season, especially evident when concentrations during the dry summer of 2007 were less than half of those in the wetter summers of 2008 and 2009. This modeling approach may help identify the principal factors that control variation in surface water FMeHg concentrations in other settings, which can guide the appropriate application of process-based models.

New York↗
USGS1991

Fractal hierarchies of magma transport in Hawaii and critical self- organization of tremor

A hierarchical model of magma transport in Hawaii is developed from the seismic records of deep (30-60 km) and intermediate-depth (5-15 km) harmonic tremor between January 1, 1962, and December 31, 1983. A tremor model of magma transport is developed from mass balances of percolation that are proportional to tremor durations. It gives reasonable magma fractions and residence times for a vertical drift velocity of 4 km yr -1 and yields patterns of intermittency that are in accord with singularity analyses of the 22-year time series record. It is suggested that spatiotemporal universality extends from small to large scales in Hawaiian and other magmatic systems. The apparent universal scaling of frequencies may be more than 15 decades in time (1 s to about 60 m.y.) and 10 decades in length (0.1 mm to 10 3 km). -from Authors

Hawaii↗
USGS2022

Combining eddy covariance and chamber methods to better constrain CO2 and CH4 fluxes across a heterogeneous restored tidal wetland

Tidal wetlands play an important role in global carbon cycling by storing carbon in sediment at millennial time scales, transporting dissolved carbon into coastal waters, and contributing significantly to global CH 4 budgets. However, these ecosystems' greenhouse gas monitoring and predictions are challenging due to spatial heterogeneity and tidal flooding. We utilized eddy covariance and chamber measurements to quantify fluxes of CO 2 and CH 4 at a restored tidal saltmarsh across spatial and temporal scales. Eddy covariance data revealed that the site was a strong net sink for CO 2 (−387 g C-CO 2 m −2 yr −1 , SD = 46) and a small net source of CH 4 (0.7 g C-CH 4 m −2 yr −1 , SD = 0.4). After partitioning net ecosystem exchange of CO 2 into gross primary production and ecosystem respiration, we found that high net uptake of CO 2 was due to low respiration emissions rather than high photosynthetic rates. We also found that respiration rates varied between land covers with increased respiration in mudflats compared to vegetated areas. Daytime soil chamber measurements revealed that the greatest CO 2 emission was from higher elevation mudflat soils (0.5 μmol m −2 s −1 , SE = 1.3) and CH 4 emission was greatest from lower elevation Spartina foliosa soils (1.6 nmol m −2 s −1 , SD = 8.2). Overall, these results highlight the importance of the relationships between wetland plant community and elevation, and inundation for CO 2 and CH 4 fluxes. Future research should include the use of high-resolution imagery, automated chambers, and a focus on quantifying carbon exported in tidal waters.

California↗
USGS2019

Basal stress equations for granular debris masses on smooth or discretized slopes

Knowledge of basal stresses is essential for analyzing slope stability and modeling the dynamics and erosive potential of debris flows and avalanches. Here we derive and test new algebraic formulas for calculating the shear stress τ and normal stress σ at the base of variable‐thickness granular debris masses in states of static or dynamic equilibrium on slopes. The formulas include a lateral pressure coefficient κ , but use of a fixed value κ = 0.7 yields predictions of σ that on average err by less than 3% and of τ that on average err by less than 13% in matching basal stresses measured in six large‐scale experiments involving wet debris masses with varying geometries and compositions. Much larger prediction errors result from use of infinite‐slope or shallow‐debris approximations. Specialized versions of the new formulas apply if basal topography is discretized and represented by a “staircase” function in a digital elevation model. Use of these formulas to assess static limiting equilibrium conditions shows that the apparent basal Coulomb friction angle ϕ tread of debris that engages friction acting on the horizontal surfaces (or “treads”) of a staircase sloping at an angle θ is generally described by tan ϕ tread = tan ( ϕ − θ )+ κ tan θ , where ϕ is the true basal friction angle of the same debris in contact with a uniformly sloping bed. Differences between the values of ϕ and ϕ tread can greatly influence the results of numerical simulations that use unsmoothed digital elevation model topography to calculate the stability or dynamics of debris masses on slopes.

Journal of Geophysical Research-Earth Surface↗
USGS2021

Temporal seismic velocity variations: Recovery following from the 2019 Mw 7.1 Ridgecrest, California earthquake

We investigated seismic velocity changes ( dv / v ) associated with the 2019 Ridgecrest earthquake sequence with high‐frequency autocorrelations of ambient seismic noise data. Daily autocorrelation functions were computed for the entirety of 2019 and the first quarter of 2020 for broadband stations within the region, including the temporary broadband stations installed during the aftershock deployment. Travel time shifts in the daily autocorrelation functions, relative to the mean autocorrelation waveform, were computed to produce dv / v time series, which are sensitive to the evolving material properties of the shallow crust surrounding the Ridgecrest fault zone (RFZ). A short‐term velocity drop follows the M w 7.1 earthquake at stations in the vicinity of the rupture surface, while those greater than 50 km away showed no such drop. The maximum, absolute changes in seismic velocity are proportional to the logarithm of distance from the fault rupture and to the peak dynamic strain experienced during the earthquake. Near the areas of the highest coseismic slip within the RFZ, seismic velocities recovered over 3 months. However, in the vicinity of the nearby Garlock fault, where triggered slip manifested, and north of the RFZ, seismic velocities recovered within a month. We interpret the seismic velocity changes and their recovery to be largely due to changes in the physical properties of the shallow crust, such as fault zone damage recovery caused by the earthquake rupture process and in response to the large dynamic stresses of passing seismic waves from the mainshock.

California↗
USGS2021

A numerical study of wave-driven mean flows and setup dynamics at a coral reef-lagoon system

Two-dimensional mean wave-driven flow and setup dynamics were investigated at a reef-lagoon system at Ningaloo Reef, Western Australia, using the numerical wave-flow model, SWASH. Phase-resolved numerical simulations of the wave and flow fields, validated with highly detailed field observations (including >10 sensors through the energetic surf zone), were used to quantify the main mechanisms that govern the mean momentum balances and resulting mean current and setup patterns, with particular attention to the role of nonlinear wave shapes. Momentum balances from the phase-resolved model indicated that onshore flows near the reef crest were primarily driven by the wave force (dominated by radiation stress gradients) due to intense breaking, whereas the flow over the reef flat and inside the lagoon and channels was primarily driven by a pressure gradient. Wave setup inside the lagoon was primarily controlled by the wave force and bottom stress. The bottom stress reduced the setup on the reef flat and inside the lagoon. Excluding the bottom stress contribution in the setup balance resulted in an over prediction of the wave-setup inside the lagoon by up to 200–370%. The bottom stress was found to be caused by the combined presence of onshore directed wave-driven currents and (nonlinear) waves. Exclusion of the bottom stress contribution from nonlinear wave shapes led to an over prediction of the setup inside the lagoon by approximately 20–40%. The inclusion of the nonlinear wave shape contribution to the bottom stress term was found to be particularly relevant in reef regions that experience a net onshore mass flux over the reef crest.

Journal of Geophysical Research-Oceans↗
USGS2021

Extreme precipitation across adjacent burned and unburned watersheds reveals impacts of low severity wildfire on debris-flow processes

In steep landscapes, wildfire-induced changes to soil and vegetation can lead to extreme and hazardous geomorphic responses, including debris flows. The wildfire-induced mechanisms that lead to heightened geomorphic responses, however, depend on many site-specific factors including regional climate, vegetation, soil texture, and soil burn severity. As climate and land use change drive changes in fire regime, there is an increasing need to understand how fire alters geomorphic responses, particularly in areas where fire has been historically infrequent. Here, we examine differences in the initiation, magnitude, and particle-size distribution of debris flows that initiated within the area burned by the 2019 Woodbury Fire in central Arizona, USA, and those that initiated in a nearby unburned area. Despite similar rainfall intensities, unburned watersheds were less likely to produce debris flows. Debris flows in unburned areas initiated from both runoff and shallow landslides, while debris flows only initiated from runoff-related processes in the burned area. The grain-size distribution making up the matrix of debris-flow deposits within the burned area generally had a lower ratio of sand to silt relative to debris flows that initiated in the unburned area, though there were no systematic differences in the coarse fraction of debris-flow sediment between burned and unburned areas. Results help expand our ability to predict postwildfire debris-flow activity in a wider range of settings, specifically the Sonoran Desert ecoregion, and provide general insight into the impact of wildfire on geomorphic processes in steep terrain.

Arizona↗
USGS2022

Synthetic evaluation of infrasonic multipole waveform inversion

Acoustic source inversions estimate the mass flow rate of volcanic explosions or yield of chemical explosions and provide insight into potential source directionality. However, the limitations of applying these methods to complex sources and their ability to resolve a stable solution have not been investigated in detail. We perform synthetic infrasound waveform inversions that use 3-D Green’s functions for a variety of idealized and realistic deployment scenarios using both a flat plane and Yasur volcano, Vanuatu as examples. We investigate the ability of various scenarios to retrieve the input source functions and relative amplitudes for monopole and multipole (monopole and dipole) inversions. Infrasound waveform inversions appear to be a robust method to quantify mass flow rates from simple sources (monopole) using deployments of infrasound sensors placed around a source, but care should be taken when analyzing and interpreting results from more complex acoustic sources (multipole) that have significant directional components. In the examples we consider the solution is stable for monopole inversions with a signal-to-noise ratio greater than five and the dipole component is small. For most scenarios investigated, the vertical dipole component of the multipole explosion source is poorly constrained and can impact the ability to recover the other source term components. Because multipole inversions are ill-posed for many deployments, a low residual does not necessarily mean the proper source vector has been recovered. Synthetic studies can help investigate the limitations and place bounds on information that may be missing using monopole and multipole inversions for potentially directional sources.

Journal of Geophysical Research-Solid Earth↗
USGS2021

Influence of permafrost type and site history on losses of permafrost carbon after thaw

We quantified permafrost peat plateau and post-thaw carbon (C) stocks across a chronosequence in Interior Alaska to evaluate the amount of C lost with thaw. Macrofossil reconstructions revealed three stratigraphic layers of peat: (1) a base layer of fen/marsh peat, (2) peat from a forested peat plateau (with permafrost) and, (3) collapse-scar bog peat (at sites where permafrost thaw has occurred). Radiocarbon dating revealed that peat initiated within the last 2,500 years and that permafrost aggraded during the Little Ice Age (ca. 250 – 575 years ago) and degraded within the last several decades. The timing of permafrost thaw within each feature was not related to thaw bog size. Their rate of expansion may be more influenced by local factors, such as ground ice content and subsurface water inputs. We found C losses due to thaw over the past century were up to 46% of the C available, but the absolute amount of C lost (kg m -2 ) was over 50% lower than losses previously described in other Alaskan peatland chronosequences. We hypothesize that this difference stems from the process by which permafrost aggraded, with sites that formed permafrost epigenetically (significantly later than most peat accumulation) experiencing less absolute C loss with thaw than sites that formed syngenetically (simultaneously with peat accumulation). Epigenetic peat from our site had lower C:N ratios as compared to Alaskan sites that have syngenetic peat. This difference could help predict the magnitude of C loss with thaw across a range or permafrost types and histories.

Alaska↗