SEARCH · USGS Science

Results for “Geophysical Journal”

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

SourceResearch connectionFrom yearThrough year

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

At least 541 records · Page 30Linked to original sources

USGS2015

Context of ancient aqueous environments on Mars from in situ geologic mapping at Endeavour Crater

Using the Mars Exploration Rover Opportunity , we have compiled one of the first field geologic maps on Mars while traversing the Noachian terrain along the rim of the 22 km diameter Endeavour Crater (Latitude −2°16′33″, Longitude −5°10′51″). In situ mapping of the petrographic, elemental, structural, and stratigraphic characteristics of outcrops and rocks distinguishes four mappable bedrock lithologic units. Three of these rock units predate the surrounding Burns formation sulfate-rich sandstones and one, the Matijevic Formation, represents conditions on early Mars predating the formation of Endeavour Crater. The stratigraphy assembled from these observations includes several geologic unconformities. The differences in lithologic units across these unconformities record changes in the character and intensity of the Martian aqueous environment over geologic time. Water circulated through fractures in the oldest rocks over periods long enough that texturally and elementally significant alteration occurred in fracture walls. These oldest pre-Endeavour rocks and their network of mineralized and altered fractures were preserved by burial beneath impact ejecta and were subsequently exhumed and exposed. The alteration along joints in the oldest rocks and the mineralized veins and concentrations of trace metals in overlying lithologic units is direct evidence that copious volumes of mineralized and/or hydrothermal fluids circulated through the early Martian crust. The wide range in intensity of structural and chemical modification from outcrop to outcrop along the crater rim shows that the ejecta of large (>8 km in diameter) impact craters is complex. These results imply that geologic complexity is to be anticipated in other areas of Mars where cratering has been a fundamental process in the local and regional geology and mineralogy.

Journal of Geophysical Research E: Planets↗

USGS2015

Evolution of Mars’ Northern Polar Seasonal CO 2 deposits: variations in surface brightness and bulk density

Small scale variations of seasonal ice are explored at different geomorphic units on the Northern Polar Seasonal Cap (NPSC). We use seasonal rock shadow measurements, combined with visible and thermal observations, to calculate density over time. The coupling of volume density and albedo allows us to determine the microphysical state of the seasonal CO 2 ice. We find two distinct endmembers across the NPSC: 1) Snow deposits may anneal to form an overlying slab layer that fractures. These low density deposits maintain relatively constant densities over springtime. 2) Porous slab deposits likely anneal rapidly in early spring and fracture in late spring. These high density deposits dramatically increase in density over time. The endmembers appear to be correlated with latitude.

Journal of Geophysical Research E: Planets↗

USGS2015

Atypical pit craters on Mars: new insights from THEMIS, CTX and HiRISE observations

More than 100 pit craters in the Tharsis region of Mars exhibit morphologies, diameters and thermal behaviors that diverge from the much larger bowl-shaped pit craters that occur in most regions across Mars. These Atypical Pit Craters (APCs) generally have sharp and distinct rims, vertical or overhanging walls that extend down to their floors, surface diameters of ~50-350 m, and high depth-to-diameter (d/D) ratios that are usually greater than 0.3 (which is an upper-range value for impacts and bowl-shaped pit craters), and can exceed values of 1.8. Observations by the Mars Odyssey THermal Emission Imaging System (THEMIS) show that APC floor temperatures are warmer at night, and fluctuate with much lower diurnal amplitudes than nearby surfaces or adjacent bowl-shaped pit craters. Kīlauea volcano, Hawai'i, hosts pit craters that formed through subsurface collapse into active volcanic dikes, resulting in pits that can appear morphologically analogous to either APCs or bowl-shaped pit craters. Partially-drained dikes are sometimes exposed within the lower walls and floors of these terrestrial APC analogs and can form extensive cave systems with unique microclimates. Similar caves in martian pit craters are of great interest for astrobiology. This study uses new observations by the Mars Reconnaissance Orbiter (MRO) High Resolution Imaging Science Experiment (HiRISE) and Context Camera (CTX) to refine previous work where seven APCs were described from lower-resolution THEMIS visible-wavelength (VIS) observations. Here, we identify locations of 115 APCs, map their distribution across the Tharsis region, characterize their internal morphologies with high-resolution observations, and discuss possible formation mechanisms.

Journal of Geophysical Research↗

USGS2015

Examining spectral variations in localized lunar dark mantle deposits

The localized lunar dark mantle deposits (DMDs) in Alphonsus, J. Herschel, and Oppenheimer craters were analyzed using visible-near-infrared spectroscopy data from the Moon Mineralogy Mapper. Spectra of these localized DMDs were analyzed for compositional and mineralogical variations within the deposits and were compared with nearby mare basalt units. Spectra of the three localized DMDs exhibited mafic absorption features indicating iron-rich compositions, although the DMDs were spectrally distinct from nearby mare basalts. All of the DMDs contained spectral signatures of glassy materials, suggesting the presence of volcanic glass in varying concentrations across the individual deposits. In addition, the albedo and spectral signatures were variable within the Alphonsus and Oppenheimer crater DMDs, suggesting variable deposit thickness and/or variations in the amount of mixing with the local substrate. Two previously unidentified localized DMDs were discovered to the northeast of Oppenheimer crater. The identification of high concentrations of volcanic glass in multiple localized DMDs in different locations suggests that the distribution of volcanic glass across the lunar surface is much more widespread than has been previously documented. The presence of volcanic glass implies an explosive, vulcanian eruption style for localized DMDs, as this allows volcanic glass to rapidly quench, inhibiting crystallization, compared to the larger hawaiian-style eruptions typical of regional DMD emplacement where black beads indicate a higher degree of crystallization. Improved understanding of the local and global distributions of volcanic glass in lunar DMDs will further constrain lunar degassing and compositional evolution throughout lunar volcanic history.

Journal of Geophysical Research E: Planets↗

USGS2015

Limited role for thermal erosion by turbulent lava in proximal Athabasca Valles, Mars

The Athabasca Valles flood lava is among the most recent (<50 Ma) and best preserved effusive lava flows on Mars and was probably emplaced turbulently. The Williams et al. (2005) model of thermal erosion by lava has been applied to what we term “proximal Athabasca,” the 75 km long upstream portion of Athabasca Valles. For emplacement volumes of 5000 and 7500 km 3 and average flow thicknesses of 20 and 30 m, the duration of the eruption varies between ~11 and ~37 days. The erosion of the lava flow substrate is investigated for three eruption temperatures (1270°C, 1260°C, and 1250°C), and volatile contents equivalent to 0–65 vol % bubbles. The largest erosion depths of ~3.8–7.5 m are at the lava source, for 20 m thick and bubble-free flows that erupted at their liquidus temperature (1270°C). A substrate containing 25 vol % ice leads to maximum erosion. A lava temperature 20°C below liquidus reduces erosion depths by a factor of ~2.2. If flow viscosity increases with increasing bubble content in the lava, the presence of 30–50 vol % bubbles leads to erosion depths lower than those relative to bubble-free lava by a factor of ~2.4. The presence of 25 vol % ice in the substrate increases erosion depths by a factor of 1.3. Nevertheless, modeled erosion depths, consistent with the emplacement volume and flow duration constraints, are far less than the depth of the channel (~35–100 m). We conclude that thermal erosion does not appear to have had a major role in excavating Athabasca Valles.

Journal of Geophysical Research E: Planets↗

USGS2014

Residual shear strength variability as a primary control on movement of landslides reactivated by earthquake-induced ground motion: Implications for coastal Oregon, U.S.

Most large seismogenic landslides are reactivations of preexisting landslides with basal shear zones in the residual strength condition. Residual shear strength often varies during rapid displacement, but the response of residual shear zones to seismic loading is largely unknown. We used a ring shear apparatus to perform simulated seismic loading tests, constant displacement rate tests, and tests during which shear stress was gradually varied on specimens from two landslides to improve understanding of coseismic landslide reactivation and to identify shear strength models valid for slow gravitational failure through rapid coseismic failure. The landslides we studied represent many along the Oregon, U.S., coast. Seismic loading tests resulted in (1) catastrophic failure involving unbounded displacement when stresses represented those for the existing landslides and (2) limited to unbounded displacement when stresses represented those for hypothetical dormant landslides, suggesting that coseismic landslide reactivation may be significant during future great earthquakes occurring near the Oregon Coast. Constant displacement rate tests indicated that shear strength decreased exponentially during the first few decimeters of displacement but increased logarithmically with increasing displacement rate when sheared at 0.001 cm s −1 or greater. Dynamic shear resistance estimated from shear strength models correlated well with stresses observed during seismic loading tests, indicating that displacement rate and amount primarily controlled failure characteristics. We developed a stress-based approach to estimate coseismic landslide displacement that utilizes the variable shear strength model. The approach produced results that compared favorably to observations made during seismic loading tests, indicating its utility for application to landslides.

Oregon↗

USGS2014

Characterizing riverbed sediment using high-frequency acoustics 1: spectral properties of scattering

Bed-sediment classification using high-frequency hydro-acoustic instruments is challenging when sediments are spatially heterogeneous, which is often the case in rivers. The use of acoustic backscatter to classify sediments is an attractive alternative to analysis of topography because it is potentially sensitive to grain-scale roughness. Here, a new method is presented which uses high-frequency acoustic backscatter from multibeam sonar to classify heterogeneous riverbed sediments by type (sand, gravel,rock) continuously in space and at small spatial resolution. In this, the first of a pair of papers that examine the scattering signatures from a heterogeneous riverbed, methods are presented to construct spatially explicit maps of spectral properties from geo-referenced point clouds of geometrically and radiometrically corrected echoes. Backscatter power spectra are computed to produce scale and amplitude metrics that collectively characterize the length scales of stochastic measures of riverbed scattering, termed ‘stochastic geometries’. Backscatter aggregated over small spatial scales have spectra that obey a power-law. This apparently self-affine behavior could instead arise from morphological- and grain-scale roughnesses over multiple overlapping scales, or riverbed scattering being transitional between Rayleigh and geometric regimes. Relationships exist between stochastic geometries of backscatter and areas of rough and smooth sediments. However, no one parameter can uniquely characterize a particular substrate, nor definitively separate the relative contributions of roughness and acoustic impedance (hardness). Combinations of spectral quantities do, however, have the potential to delineate riverbed sediment patchiness, in a data-driven approach comparing backscatter with bed-sediment observations (which is the subject of part two of this manuscript).

Journal of Geophysical Research F: Earth Surface↗

USGS2014

Characterizing riverbed sediment using high-frequency acoustics 2: scattering signatures of Colorado River bed sediment in Marble and Grand Canyons

In this, the second of a pair of papers on the statistical signatures of riverbed sediment in high-frequency acoustic backscatter, spatially explicit maps of the stochastic geometries (length- and amplitude-scales) of backscatter are related to patches of riverbed surfaces composed of known sediment types, as determined by geo-referenced underwater video observations. Statistics of backscatter magnitudes alone are found to be poor discriminators between sediment types. However, the variance of the power spectrum, and the intercept and slope from a power-law spectral form (termed the spectral strength and exponent, respectively) successfully discriminate between sediment types. A decision-tree approach was able to classify spatially heterogeneous patches of homogeneous sands, gravels (and sand-gravel mixtures), and cobbles/boulders with 95, 88, and 91% accuracy, respectively. Application to sites outside the calibration, and surveys made at calibration sites at different times, were plausible based on observations from underwater video. Analysis of decision trees built with different training data sets suggested that the spectral exponent was consistently the most important variable in the classification. In the absence of theory concerning how spatially variable sediment surfaces scatter high-frequency sound, the primary advantage of this data-driven approach to classify bed sediment over alternatives is that spectral methods have well understood properties and make no assumptions about the distributional form of the fluctuating component of backscatter over small spatial scales.

Arizona↗

USGS2015

Glaciological and marine geological controls on terminus dynamics of Hubbard Glacier, southeast Alaska

Hubbard Glacier, located in southeast Alaska, is the world's largest non-polar tidewater glacier. It has been steadily advancing since it was first mapped in 1895; occasionally, the advance creates an ice or sediment dam that blocks a tributary fjord (Russell Fiord). The sustained advance raises the probability of long-term closure in the near-future, which will strongly impact the ecosystem of Russell Fiord and the nearby community of Yakutat. Here, we examine a 43-year record of flow speeds and terminus position to understand the large-scale dynamics of Hubbard Glacier. Our long-term record shows that the rate of terminus advance has increased slightly since 1895, with the exception of a slowed advance between approximately 1972 and 1984. The short-lived closure events in 1986 and 2002 were not initiated by perturbations in ice velocity or environmental forcings, but were likely due to fluctuations in sedimentation patterns at the terminus. This study points to the significance of a coupled system where short-term velocity fluctuations and morainal shoal development control tidewater glacier terminus position.

Alaska↗

USGS2015

Radar attenuation and temperature within the Greenland Ice Sheet

The flow of ice is temperature-dependent, but direct measurements of englacial temperature are sparse. The dielectric attenuation of radio waves through ice is also temperature-dependent, and radar sounding of ice sheets is sensitive to this attenuation. Here we estimate depth-averaged radar-attenuation rates within the Greenland Ice Sheet from airborne radar-sounding data and its associated radiostratigraphy. Using existing empirical relationships between temperature, chemistry, and radar attenuation, we then infer the depth-averaged englacial temperature. The dated radiostratigraphy permits a correction for the confounding effect of spatially varying ice chemistry. Where radar transects intersect boreholes, radar-inferred temperature is consistently higher than that measured directly. We attribute this discrepancy to the poorly recognized frequency dependence of the radar-attenuation rate and correct for this effect empirically, resulting in a robust relationship between radar-inferred and borehole-measured depth-averaged temperature. Radar-inferred englacial temperature is often lower than modern surface temperature and that of a steady state ice-sheet model, particularly in southern Greenland. This pattern suggests that past changes in surface boundary conditions (temperature and accumulation rate) affect the ice sheet's present temperature structure over a much larger area than previously recognized. This radar-inferred temperature structure provides a new constraint for thermomechanical models of the Greenland Ice Sheet.

Journal of Geophysical Research F: Earth Surface↗

USGS2015

Coevolution of bed surface patchiness and channel morphology: 1. Mechanisms of forced patch formation

Riverbeds frequently display a spatial structure where the sediment mixture composing the channel bed has been sorted into discrete patches of similar grain size. Even though patches are a fundamental feature in gravel bed rivers, we have little understanding of how patches form, evolve, and interact. Here we present a two-dimensional morphodynamic model that is used to examine in greater detail the mechanisms responsible for the development of forced bed surface patches and the coevolution of bed morphology and bed surface patchiness. The model computes the depth-averaged channel hydrodynamics, mixed-grain-size sediment transport, and bed evolution by coupling the river morphodynamic model Flow and Sediment Transport with Morphological Evolution of Channels (FaSTMECH) with a transport relation for gravel mixtures and the mixed-grain-size Exner equation using the active layer assumption. To test the model, we use it to simulate a flume experiment in which the bed developed a sequence of alternate bars and temporally and spatially persistent forced patches with a general pattern of coarse bar tops and fine pools. Cross-stream sediment flux causes sediment to be exported off of bars and imported into pools at a rate that balances downstream gradients in the streamwise sediment transport rate, allowing quasi-steady bar-pool topography to persist. The relative importance of lateral gravitational forces on the cross-stream component of sediment transport is a primary control on the amplitude of the bars. Because boundary shear stress declines as flow shoals over the bars, the lateral sediment transport is increasingly size selective and leads to the development of coarse bar tops and fine pools.

Journal of Geophysical Research F: Earth Surface↗

USGS2015

Coevolution of bed surface patchiness and channel morphology: 2. Numerical experiments

In gravel bed rivers, bed topography and the bed surface grain size distribution evolve simultaneously, but it is not clear how feedbacks between topography and grain sorting affect channel morphology. In this, the second of a pair of papers examining interactions between bed topography and bed surface sorting in gravel bed rivers, we use a two-dimensional morphodynamic model to perform numerical experiments designed to explore the coevolution of both free and forced bars and bed surface patches. Model runs were carried out on a computational grid simulating a 200 m long, 2.75 m wide, straight, rectangular channel, with an initially flat bed at a slope of 0.0137. Over five numerical experiments, we varied (a) whether an obstruction was present, (b) whether the sediment was a gravel mixture or a single size, and (c) whether the bed surface grain size feeds back on the hydraulic roughness field. Experiments with channel obstructions developed a train of alternate bars that became stationary and were connected to the obstruction. Freely migrating alternate bars formed in the experiments without channel obstructions. Simulations incorporating roughness feedbacks between the bed surface and flow field produced flatter, broader, and longer bars than simulations using constant roughness or uniform sediment. Our findings suggest that patches are not simply a by-product of bed topography, but they interact with the evolving bed and influence morphologic evolution.

Journal of Geophysical Research F: Earth Surface↗

USGS2014

Differential response of carbon fluxes to climate in three peatland ecosystems that vary in the presence and stability of permafrost

Changes in vegetation and soil properties following permafrost degradation and thermokarst development in peatlands may cause changes in net carbon storage. To better understand these dynamics, we established three sites in Alaska that vary in permafrost regime, including a black spruce peat plateau forest with stable permafrost, an internal collapse scar bog formed as a result of thermokarst, and a rich fen without permafrost. Measurements include year-round eddy covariance estimates of carbon dioxide (CO 2 ), water, and energy fluxes, associated environmental variables, and methane (CH 4 ) fluxes at the collapse scar bog. The ecosystems all acted as net sinks of CO 2 in 2011 and 2012, when air temperature and precipitation remained near long-term means. In 2013, under a late snowmelt and late leaf out followed by a hot, dry summer, the permafrost forest and collapse scar bog were sources of CO 2 . In this same year, CO 2 uptake in the fen increased, largely because summer inundation from groundwater inputs suppressed ecosystem respiration. CO 2 exchange in the permafrost forest and collapse scar bog was sensitive to warm air temperatures, with 0.5 g C m −2 lost each day when maximum air temperature was very warm (≥29°C). The bog lost 4981 ± 300 mg CH 4 m −2 between April and September 2013, indicating that this ecosystem acted as a significant source of both CO 2 and CH 4 to the atmosphere in 2013. These results suggest that boreal peatland responses to warming and drying, both of which are expected to occur in a changing climate, will depend on permafrost regime.

Alaska↗

USGS2014

Radiocarbon age-offsets in an arctic lake reveal the long-term response of permafrost carbon to climate change

Continued warming of the Arctic may cause permafrost to thaw and speed the decomposition of large stores of soil organic carbon (OC), thereby accentuating global warming. However, it is unclear if recent warming has raised the current rates of permafrost OC release to anomalous levels or to what extent soil carbon release is sensitive to climate forcing. Here we use a time series of radiocarbon age-offsets ( 14 C) between the bulk lake sediment and plant macrofossils deposited in an arctic lake as an archive for soil and permafrost OC release over the last 14,500 years. The lake traps and archives OC imported from the watershed and allows us to test whether prior warming events stimulated old carbon release and heightened age-offsets. Today, the age-offset (2 ka; thousand of calibrated years before A.D. 1950) and the depositional rate of ancient OC from the watershed into the lake are relatively low and similar to those during the Younger Dryas cold interval (occurring 12.9–11.7 ka). In contrast, age-offsets were higher (3.0–5.0 ka) when summer air temperatures were warmer than present during the Holocene Thermal Maximum (11.7–9.0 ka) and Bølling-Allerød periods (14.5–12.9 ka). During these warm times, permafrost thaw contributed to ancient OC depositional rates that were ~10 times greater than today. Although permafrost OC was vulnerable to climate warming in the past, we suggest surface soil organic horizons and peat are presently limiting summer thaw and carbon release. As a result, the temperature threshold to trigger widespread permafrost OC release is higher than during previous warming events.

Alaska↗

USGS2014

Using dissolved organic matter age and composition to detect permafrost thaw in boreal watersheds of interior Alaska

Recent warming at high latitudes has accelerated permafrost thaw, which can modify soil carbon dynamics and watershed hydrology. The flux and composition of dissolved organic matter (DOM) from soils to rivers are sensitive to permafrost configuration and its impact on subsurface hydrology and groundwater discharge. Here, we evaluate the utility of DOM composition and age as a tool for detecting permafrost thaw in three rivers (Beaver, Birch, and Hess Creeks) within the discontinuous permafrost zone of interior Alaska. We observed strong temporal controls on Δ 14 C content of hydrophobic acid isolates (Δ 14 C-HPOA) across all rivers, with the most enriched values occurring during spring snowmelt (75 ± 8&permil;) and most depleted during winter flow (−21 ± 8&permil;). Radiocarbon ages of winter flow samples ranged from 35 to 445 yr BP, closely tracking estimated median base flow travel times for this region (335 years). During spring snowmelt, young DOM was composed of highly aromatic, high molecular-weight compounds, whereas older DOM of winter flow had lower aromaticity and molecular weight. We observed a significant correlation between Δ 14 C-HPOA and UV absorbance coefficient at 254 nm ( α 254 ) across all study rivers. Using α 254 as an optical indicator for Δ 14 C-HPOA, we also observed a long-term decline in α 254 during maximum annual thaw depth over the last decade at the Hess Creek study site. These findings suggest a shift in watershed hydrology associated with increasing active layer thickness. Further development of DOM optical indicators may serve as a novel and inexpensive tool for detecting permafrost degradation in northern watersheds.

Alaska↗

USGS2014

Uncertainty in the fate of soil organic carbon: A comparison of three conceptually different soil decomposition models

Conventional Q10 soil organic matter decomposition models and more complex microbial models are available for making projections of future soil carbon dynamics. However, it is unclear (1) how well the conceptually different approaches can simulate observed decomposition and (2) to what extent the trajectories of long-term simulations differ when using the different approaches. In this study, we compared three structurally different soil carbon (C) decomposition models (one Q10 and two microbial models of different complexity), each with a one- and two-horizon version. The models were calibrated and validated using 4 years of measurements of heterotrophic soil CO 2 efflux from trenched plots in a Dahurian larch ( Larix gmelinii Rupr.) plantation. All models reproduced the observed heterotrophic component of soil CO 2 efflux, but the trajectories of soil carbon dynamics differed substantially in 100 year simulations with and without warming and increased litterfall input, with microbial models that produced better agreement with observed changes in soil organic C in long-term warming experiments. Our results also suggest that both constant and varying carbon use efficiency are plausible when modeling future decomposition dynamics and that the use of a short-term (e.g., a few years) period of measurement is insufficient to adequately constrain model parameters that represent long-term responses of microbial thermal adaption. These results highlight the need to reframe the representation of decomposition models and to constrain parameters with long-term observations and multiple data streams. We urge caution in interpreting future soil carbon responses derived from existing decomposition models because both conceptual and parameter uncertainties are substantial.

Journal of Geophysical Research G: Biogeosciences↗

USGS2014

Nitrogen cycling processes and microbial community composition in bed sediments in the Yukon River at Pilot Station

Information on the contribution of nitrogen (N)-cycling processes in bed sediments to river nutrient fluxes in large northern latitude river systems is limited. This study examined the relationship between N-cycling processes in bed sediments and N speciation and loading in the Yukon River near its mouth at the Bering Sea. We conducted laboratory bioassays to measure N-cycling processes in sediment samples collected over distinct water cycle seasons. In conjunction, the microbial community composition in the bed sediments using genes involved in N-cycling ( narG , napA , nosZ , and amoA ) and 16S rRNA gene pyrosequences was examined. Temporal variation was observed in net N mineralization, nitrate uptake, and denitrification rate potentials and correlated strongly with sediment carbon (C) and extractable N content and microbial community composition rather than with river water nutrient concentrations. The C content of the bed sediment was notably impacted by the spring flood, ranging from 1.1% in the midst of an ice-jam to 0.1% immediately after ice-out, suggesting a buildup of organic material (OM) prior to scouring of the bed sediments during ice break up. The dominant members of the microbial community that explained differences in N-processing rates belonged to the genera Crenothrix , Flavobacterium , and the family of Comamonadaceae. Our results suggest that biogeochemical processing rates in the bed sediments appear to be more coupled to hydrology, nutrient availability in the sediments, and microbial community composition rather than river nutrient concentrations at Pilot Station.

Alaska↗

USGS2014

Temporal variability of carbon and nutrient burial, sediment accretion, and mass accumulation over the past century in a carbonate platform mangrove forest of the Florida Everglades.

The objective of this research was to measure temporal variability in accretion and mass sedimentation rates (including organic carbon (OC), total nitrogen (TN), and total phosphorous (TP)) from the past century in a mangrove forest on the Shark River in Everglades National Park, USA. The 210 Pb Constant Rate of Supply model was applied to six soil cores to calculate annual rates over the most recent 10, 50, and 100 year time spans. Our results show that rates integrated over longer timeframes are lower than those for shorter, recent periods of observation. Additionally, the substantial spatial variability between cores over the 10 year period is diminished over the 100 year record, raising two important implications. First, a multiple-decade assessment of soil accretion and OC burial provides a more conservative estimate and is likely to be most relevant for forecasting these rates relative to long-term processes of sea level rise and climate change mitigation. Second, a small number of sampling locations are better able to account for spatial variability over the longer periods than for the shorter periods. The site average 100 year OC burial rate, 123 ± 19 (standard deviation) g m -2 yr -1 , is low compared with global mangrove values. High TN and TP burial rates in recent decades may lead to increased soil carbon remineralization, contributing to the low carbon burial rates. Finally, the strong correlation between OC burial and accretion across this site signals the substantial contribution of OC to soil building in addition to the ecosystem service of CO 2 sequestration.

Florida↗