SEARCH · Geology Reports

Results for “Soil Research”

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 1,315 records · Page 73Linked to original sources

USGS1998

The discharge of nitrate-contaminated groundwater from developed shoreline to marsh-fringed estuary

As residential development, on-site wastewater disposal, and groundwater contamination increase in the coastal zone, assessment of nutrient removal by soil and sedimentary processes becomes increasingly important. Nitrogen removal efficiency depends largely on the specific flow paths taken by groundwater as it discharges into nitrogen-limited estuarine waters. Shoreline salinity surveys, hydraulic studies, and thermal infrared imagery indicated that groundwater discharge into the Nauset Marsh estuary (Eastham, Massachusetts) occurred in high-velocity seeps immediately seaward of the upland-fringing salt marsh. Discharge was highly variable spatially and occurred through permeable, sandy sediments during low tide. Seepage chamber monitoring showed that dissolved inorganic nitrogen (principally nitrate) traversed nearly conservatively from the aquifer through shallow estuarine sediments to coastal waters at flux rates of 1–3 mmol m −2 h −1 . A significant relationship between pore water NO 3 -N concentrations and NO 3 -N flux rates may provide a rapid method of estimating nitrogen loading from groundwater to the water column.

Massachusets↗

USGS2013

Influence of management and precipitation on carbon fluxes in greatplains grasslands

Suitable management and sufficient precipitation on grasslands can provide carbon sinks. The net carbon accumulation of a site from the atmosphere, modeled as the Net Ecosystem Productivity (NEP), is a useful means to gauge carbon balance. Previous research has developed methods to integrate flux tower data with satellite biophysical datasets to estimate NEP across large regions. A related method uses the Ecosystem Performance Anomaly (EPA) as a satellite-derived indicator of disturbance intensity (e.g., livestock stocking rate, fire, and insect damage). To better understand the interactions among management, climate, and carbon dynamics, we evaluated the relationship between EPA and NEP data at the 250 m scale for grasslands in the Central Great Plains, USA (ranging from semi-arid to mesic). We also used weekly estimates of NEP to evaluate the phenology of carbon dynamics, classified by EPA (i.e., by level of disturbance impact). Results show that the cumulative carbon balance over these grasslands from 2000 to 2008 was a weak net sink of 13.7 g C m −2 yr −1 . Overall, NEP increased with precipitation (R 2 = 0.39, P < 0.05) from west to east. Disturbance influenced NEP phenology; however, climate and biophysical conditions were usually more important. The NEP response to disturbance varies by ecoregion, and more generally by grassland type, where the shortgrass prairie NEP is most sensitive to disturbance, the mixed-grass prairie displays a moderate response, and tallgrass prairie is the least impacted by disturbance (as measured by EPA). Sustainable management practices in the tallgrass and mixed-grass prairie may potentially induce a period of average net carbon sink until a new equilibrium soil organic carbon is achieved. In the shortgrass prairie, management should be considered sustainable if carbon stocks are simply maintained. The consideration of site carbon balance adds to the already difficult task of managing grasslands appropriately to site conditions. Results clarify the seasonal and interannual dynamics of NEP, specifically the influence of disturbance and moisture availability.

Ecological Indicators↗

USGS2006

Overview of the Microscopic Imager Investigation during Spirit's first 450 sols in Gusev crater

The Microscopic Imager (MI) on the Mars Exploration Rover Spirit has returned images of Mars with higher resolution than any previous camera system, allowing detailed petrographic and sedimentological studies of the rocks and soils at the Gusev landing site. Designed to simulate a geologist's hand lens, the MI is mounted on Spirit's instrument arm and can resolve objects 0.1 mm in size or larger. This paper provides an overview of MI operations, data calibration, processing, and analysis of MI data returned during the first 450 sols (Mars days) of the Spirit landed mission. The primary goal of this paper is to facilitate further analyses of MI data by summarizing the methods used to acquire and process the data, the radiometric and geometric accuracy of MI data products, and the availability of archival products. In addition, scientific results of the MI investigation are summarized. MI observations show that poorly sorted soils are common in Gusev crater, although aeolian bedforms have well‐sorted coarse sand grains on their surfaces. Abraded surfaces of plains rocks show igneous textures, light‐toned veins or fracture‐filling minerals, and discrete coatings. The rocks in the Columbia Hills have a wide variety of granular textures, consistent with volcaniclastic or impact origins. Case hardening and submillimeter veins observed in the rocks as well as soil crusts and cemented clods imply episodic subsurface aqueous fluid movement, which has altered multiple geologic units in the Columbia Hills. The MI also monitored Spirit's solar panels and the magnets on the rover's deck.

Journal of Geophysical Research E: Planets↗

USGS2022

Predicting flood damage probability across the conterminous United States

Floods are the leading cause of natural disaster damages in the United States, with billions of dollars incurred every year in the form of government payouts, property damages, and agricultural losses. The Federal Emergency Management Agency oversees the delineation of floodplains to mitigate damages, but disparities exist between locations designated as high risk and where flood damages occur due to land use and climate changes and incomplete floodplain mapping. We harnessed publicly available geospatial datasets and random forest algorithms to analyze the spatial distribution and underlying drivers of flood damage probability caused by excessive rainfall and overflowing water bodies across the conterminous United States. From this, we produced the first spatially complete map of flood damage probability for the nation, along with spatially explicit standard errors for four selected cities. We trained models using the locations of historical reported flood damage events ( n = 71,434) and a suite of geospatial predictors (e.g., flood severity, climate, socio-economic exposure, topographic variables, soil properties, and hydrologic characteristics). We developed independent models for each hydrologic unit code level 2 watershed and generated a flood damage probability for each 100-m pixel. Our model classified damage or no damage with an average area under the curve accuracy of 0.75; however, model performance varied by environmental conditions, with certain land cover classes (e.g., forest) resulting in higher error rates than others (e.g., wetlands). Our results identified flood damage probability hotspots across multiple spatial and regional scales, with high probabilities common in both inland and coastal regions. The highest flood damage probabilities tended to be in areas of low elevation, in close proximity to streams, with extreme precipitation, and with high urban road density. Given rapid environmental changes, our study demonstrates an efficient approach for updating flood damage probability estimates across the nation.

Environmental Research Letters↗

USGS2024

Investigating the atmospheric conditions associated with impactful shallow landslides in California (USA)

Shallow landslides are often triggered during rainfall events, which can increase subsurface soil water pressure and destabilize hillslopes. The likelihood of regional shallow landslide initiation is often assessed through a comparison of rainfall intensity and duration to pre-established thresholds. While informative for landslide warning, this exclusive focus on rainfall exceeding thresholds does not consider the meteorological conditions producing the rainfall. Here, we ask the question, are there common meteorological characteristics that lead to landslide-triggering precipitation? We develop a catalog of 18 post-1995 widespread, impactful shallow landslide events occurring within 13 storms across California, USA, where initiation time could be constrained to a ≤6-h window. We examine storm characteristics during the landslide initiation window using atmospheric reanalysis products, radar observations, and quantitative precipitation estimates. We find that, while there are some common atmospheric characteristics across landslide events, they can occur under a range of atmospheric conditions. For example, all Northern California landslide events assessed are associated with moderate to strong atmospheric rivers (ARs), while Southern California landslides feature non-AR to strong AR conditions. The storm events evaluated herein share many characteristics of hydrologically important storms in California that did not necessarily result in landslides; thus, atmospheric characteristics alone may not be sufficient to determine whether landslides will occur. However, documenting the characteristics of landslide-triggering storms defines the conditions under which landslides tend to occur, provides analog events that can be useful in forecast applications, helps define future research directions relating to atmospheric conditions and landslides, and supports interdisciplinary research efforts.

California↗

USGS2002

Controls on old and new water contributions to stream flow at some nested catchments in Vermont, USA

Factors controlling the partitioning of old and new water contributions to stream flow were investigated for three events in four catchments (three of which were nested) at Sleepers River Research Watershed in Danville, Vermont. In the 1993 snowmelt period, two-component isotopic hydrograph separations showed that new water (meltwater) inputs to the stream ranged widely from 41 to 74%, and increased with catchment size (41 to 11 125 ha) (with one exception) and with open land cover (0-73%). Peak dissolved organic carbon concentrations and relative alkalinity dilution in stream water ranked in the same order among catchments as the new water fractions, suggesting that new water followed shallow flow paths. During the 1994 snowmelt, despite similar timing and magnitude of melt inputs, the new-water contribution to stream flow ranged only from 30 to 36% in the four catchments. We conclude that the uncommonly high and variable new water fractions in streamwater during the 1993 melt were caused by direct runoff of meltwater over frozen ground, which was prevalent in open land areas during the 1993 winter. In a high-intensity summer rainstorm in 1993, new water fractions were smaller relative to the 1993 snowmelt, ranging from 28 to 46%, but they ranked in the identical catchment order. Reconciliation of the contrasting patterns of new-old water partitioning in the three events appears to require an explanation that invokes multiple processes and effects, including: 1 topographically controlled increase in surface-saturated area with increasing catchment size; 2 direct runoff over frozen ground; 3 low infiltration in agriculturally compacted soils; 4 differences in soil transmissivity, which may be more relevant under dry antecedent conditions. These data highlight some of the difficulties faced by catchment hydrologists in formulating a theory of runoff generation at varying basin scales. Copyright ?? 2002 John Wiley and Sons, Ltd.

Hydrological Processes↗

USGS2015

Heat flux from magmatic hydrothermal systems related to availability of fluid recharge

Magmatic hydrothermal systems are of increasing interest as a renewable energy source. Surface heat flux indicates system resource potential, and can be inferred from soil CO 2 flux measurements and fumarole gas chemistry. Here we compile and reanalyze results from previous CO 2 flux surveys worldwide to compare heat flux from a variety of magma-hydrothermal areas. We infer that availability of water to recharge magmatic hydrothermal systems is correlated with heat flux. Recharge availability is in turn governed by permeability, structure, lithology, rainfall, topography, and perhaps unsurprisingly, proximity to a large supply of water such as the ocean. The relationship between recharge and heat flux interpreted by this study is consistent with recent numerical modeling that relates hydrothermal system heat output to rainfall catchment area. This result highlights the importance of recharge as a consideration when evaluating hydrothermal systems for electricity generation, and the utility of CO 2 flux as a resource evaluation tool.

Journal of Volcanology and Geothermal Research↗

USGS2012

Heat flow in vapor dominated areas of the Yellowstone Plateau volcanic field: implications for the thermal budget of the Yellowstone Caldera

Characterizing the vigor of magmatic activity in Yellowstone requires knowledge of the mechanisms and rates of heat transport between magma and the ground surface. We present results from a heat flow study in two vapor dominated, acid-sulfate thermal areas in the Yellowstone Caldera, the 0.11 km 2 Obsidian Pool Thermal Area (OPTA) and the 0.25 km 2 Solfatara Plateau Thermal Area (SPTA). Conductive heat flux through a low permeability layer capping large vapor reservoirs is calculated from soil temperature measurements at >600 locations and from laboratory measurements of soil properties. The conductive heat output is 3.6 ± 0.4 MW and 7.5 ± 0.4 MW from the OPTA and the SPTA, respectively. The advective heat output from soils is 1.3 ± 0.3 MW and 1.2 ± 0.3 MW from the OPTA and the SPTA, respectively and the heat output from thermal pools in the OPTA is 6.8 ± 1.4 MW. These estimates result in a total heat output of 11.8 ± 1.4 MW and 8.8 ± 0.4 MW from OPTA and SPTA, respectively. Focused zones of high heat flux in both thermal areas are roughly aligned with regional faults suggesting that faults in both areas serve as conduits for the rising acid vapor. Extrapolation of the average heat flux from the OPTA (103 ± 2 W·m −2 ) and SPTA (35 ± 3 W·m −2 ) to the ~35 km 2 of vapor dominated areas in Yellowstone yields 3.6 and 1.2 GW, respectively, which is less than the total heat output transported by steam from the Yellowstone Caldera as estimated by the chloride inventory method (4.0 to 8.0 GW).

Journal of Geophysical Research B: Solid Earth↗

USGS2018

Geospatial data for developing nutrient SPARROW models for the Midcontinental region of Canada and the United States

Through the International Watersheds Initiative of the International Joint Commission (IJC), the SpatiallyReferenced Regressions on Watershed attributes (SPARROW) model developed by the U.S. Geological Survey (USGS) is being applied to the Great Lakes, Rainy River – Lake of the Woods and Red-Assiniboine basins. The objective of this binational application of the SPARROW model is to better understand and quantify the sources of phosphorus (P) and nitrogen (N) that contribute to regional water-quality issues like algal blooms and eutrophication in Lake Erie and other parts of the Great Lakes, as well as Lake of the Woods. Led by the IJC, a team of researchers from the National Research Council of Canada – Ocean, Coastal and River Engineering Research Centre, USGS, and IJC are extending the SPARROW modelling work previously completed for the Red-Assiniboine basin and the U.S. portions of the Great Lakes, Ohio, Upper Mississippi, and Souris-Red-Rainy river basins to cover all of the Great Lakes, Rainy River – Lake of the Woods and RedAssiniboine basins. The current effort is termed the Midcontinent SPARROW modelling study. This report describes the data used to develop the Midcontinent SPARROW models, specifically the sources of original data, assembling the data, and the processing and harmonization required between the U.S. and Canada data needed to produce these models. Details provided include the: development of a digital stream network and related catchments – most significantly in the Canadian regions of the Great Lakes and Rainy River – Lake of the Woods basins where these data were not available to create a seamless binational network across the model domain; calculation of variables to aid in the determination of in-stream and in-reservoir decay of P and N; quantification of Canadian diversions within the Midcontinental region (i.e., Lake St. Joseph, Long Lake and the Ogoki Reservoir); development of binational input nutrient sources considered for model development (i.e., land cover, inorganic farm fertilizer, manure, atmospheric deposition, point-source pollution from wastewatertreatment plants and contribution from non-modelled watersheds); and development of delivery variables considered to be most predominant (i.e., temperature, precipitation and ensuing runoff, soil permeability and clay content, slope of the catchments, and tile drainage). The majority of the geospatial data collection and processing was required for Canadian datasets because many of the U.S. datasets were already assembled for previous SPARROW model applications in the U.S. The task of harmonizing data between the U.S. and Canada was important to ensure consistency of the datasets used in the models. The harmonized digital stream network, delineated catchments and input data for each catchment (i.e., source and delivery variables), created for the Midcontinent SPARROW models, are available for download at url: https://doi.org/10.4224/300.0001.

Technical Report↗

USGS2014

Influence of stocking, site quality, stand age, low-severity canopy disturbance, and forest composition on sub-boreal aspen mixedwood carbon stocks

Low-severity canopy disturbance presumably influences forest carbon dynamics during the course of stand development, yet the topic has received relatively little attention. This is surprising because of the frequent occurrence of such events and the potential for both the severity and frequency of disturbances to increase as a result of climate change. We investigated the impacts of low-severity canopy disturbance and average insect defoliation on forest carbon stocks and rates of carbon sequestration in mature aspen mixedwood forests of varying stand age (ranging from 61 to 85 years), overstory composition, stocking level, and site quality. Stocking level and site quality positively affected the average annual aboveground tree carbon increment (C AAI ), while stocking level, site quality, and stand age positively affected tree carbon stocks (C TREE ) and total ecosystem carbon stocks (C TOTAL ). Cumulative canopy disturbance (DIST) was reconstructed using dendroecological methods over a 29-year period. DIST was negatively and significantly related to soil carbon (C SOIL ), and it was negatively, albeit marginally, related to C TOTAL . Minima in the annual aboveground carbon increment of trees (C AI ) occurred at sites during defoliation of aspen ( Populus tremuloides Michx.) by forest tent caterpillar ( Malacosoma disstria Hubner), and minima were more extreme at sites dominated by trembling aspen than sites mixed with conifers. At sites defoliated by forest tent caterpillar in the early 2000s, increased sequestration by the softwood component ( Abies balsamea (L.) Mill. and Picea glauca (Moench) Voss) compensated for overall decreases in C AI by 17% on average. These results underscore the importance of accounting for low-severity canopy disturbance events when developing regional forest carbon models and argue for the restoration and maintenance of historically important conifer species within aspen mixedwoods to enhance stand-level resilience to disturbance agents and maintain site-level carbon stocks.

Minnesota↗

USGS1998

Sedimentation of prairie wetlands

Many wetlands in the prairie pothole region are embedded within an agricultural landscape where they are subject to varying degrees of siltation. Cultivation of wetland catchment areas has exacerbated soil erosion; wetlands in agricultural fields receive more sediment from upland areas than wetlands in grassland landscapes and hence are subject to premature filling (i.e., they have shorter topographic lives). Associated impacts from increased turbidity, sediment deposition, and increased surface water input likely have impaired natural wetland functions. Although trapping of sediments by wetlands is often cited as a water quality benefit, sediment input from agricultural fields has potential to completely fill wetlands and shorten their effective life-span. Thus, the value placed on wetlands to trap sediments is in conflict with maximizing the effective topographic life of wetlands. Herein, we provide an overview of sedimentation, identify associated impacts on wetlands, and suggest remedial management strategies. We also highlight the need to evaluate the impact of agricultural practices on wetland functions from an interdisciplinary approach to facilitate development of best management practices that benefit both wetland and agricultural interests.

Great Plains Research↗

USGS2015

Plant phenolics and absorption features in vegetation reflectance spectra near 1.66 μm

Past laboratory and field studies have quantified phenolic substances in vegetative matter from reflectance measurements for understanding plant response to herbivores and insect predation. Past remote sensing studies on phenolics have evaluated crop quality and vegetation patterns caused by bedrock geology and associated variations in soil geochemistry. We examined spectra of pure phenolic compounds, common plant biochemical constituents, dry leaves, fresh leaves, and plant canopies for direct evidence of absorption features attributable to plant phenolics. Using spectral feature analysis with continuum removal, we observed that a narrow feature at 1.66 μm is persistent in spectra of manzanita, sumac, red maple, sugar maple, tea, and other species. This feature was consistent with absorption caused by aromatic C- H bonds in the chemical structure of phenolic compounds and non-hydroxylated aromatics. Because of overlapping absorption by water, the feature was weaker in fresh leaf and canopy spectra compared to dry leaf measurements. Simple linear regressions of feature depth and feature area with polyphenol concentration in tea resulted in high correlations and low errors (% phenol by dry weight) at the dry leaf ( r 2 = 0.95, RMSE = 1.0%, n = 56), fresh leaf ( r 2 = 0.79, RMSE = 2.1%, n = 56), and canopy ( r 2 = 0.78, RMSE = 1.0%, n = 13) levels of measurement. Spectra of leaves, needles, and canopies of big sagebrush and evergreens exhibited a weak absorption feature centered near 1.63 μm, short ward of the phenolic compounds, possibly consistent with terpenes. This study demonstrates that subtle variation in vegetation spectra in the shortwave infrared can directly indicate biochemical constituents and be used to quantify them. Phenolics are of lesser abundance compared to the major plant constituents but, nonetheless, have important plant functions and ecological significance. Additional research is needed to advance our understanding of the spectral influences of plant phenolics and terpenes relative to dominant leaf biochemistry (water, chlorophyll, protein/nitrogen, cellulose, and lignin).

International Journal of Applied Earth Observation↗

USGS2002

Analysis of δ15N and δ18O to differentiate NO3− sources in runoff at two watersheds in the Catskill Mountains of New York

To quantify the movement of atmospheric nitrogen deposition through two forested watersheds in the Catskill Mountains of New York, dual‐isotope analysis (δ 15 N and δ 18 O) was used to differentiate NO 3 − derived from precipitation from NO 3 − derived by microbial nitrification and to quantify the contributions of these sources to NO 3 − in drainage waters. Samples of stream water, soil water, precipitation, snowmelt, and O‐horizon soil were collected during the March and April snowmelt period of 1994 and throughout an 18‐month period from August 1995 through February 1997. The mean δ 18 O‐NO 3 − value of precipitation was +50.5‰, whereas the mean values for stream water and soil water were +17.7‰ and +23.6‰, respectively. The mean δ 15 N‐NO 3 − of precipitation was −0.2‰, that of soil water was +1.4‰, and that of stream water was +2.3‰; these values showed greater overlap among the three different waters than did the δ 18 O‐NO 3 − values, indicating that δ 15 N‐NO 3 − was not as useful for source separation. Soil water δ 18 O‐NO 3 − values decreased, and δ 15 N‐NO 3 − values increased, from the O to the B and C horizons, but most of the differences among horizons were not statistically significant. Nitrate derived by nitrification in incubated soil samples had a wide range of δ 15 N‐NO 3 − values, from +1.5‰ to +16.1‰, whereas δ 18 O‐NO 3 − values ranged more narrowly, from +13.2‰ to +16.0‰. Values of δ 18 O‐NO 3 − indicated that NO 3 − in stream water is mainly derived from nitrification. Only during a high‐flow event that exceeded the annual flood was precipitation a major contributor to stream water NO 3 − . Values of δ 18 O‐NO 3 − and δ 15 N‐NO 3 − changed at differing rates as NO 3 − cycled through these watersheds because δ 18 O‐NO 3 − values change sharply through the incorporation of oxygen from ambient water and gas during nitrification, whereas δ 15 N‐NO 3 − values change only incrementally through fractionation during biocycling processes. The results of this study show that most NO 3 − is first cycled through the biota and nitrified before entering the stream.

Water Resources Research↗

USGS2008

Veneers, rinds, and fracture fills: Relatively late alteration of sedimentary rocks at Meridiani Planum, Mars

Veneers and thicker rinds that coat outcrop surfaces and partially cemented fracture fills formed perpendicular to bedding document relatively late stage alteration of ancient sedimentary rocks at Meridiani Planum, Mars. The chemistry of submillimeter thick, buff-colored veneers reflects multiple processes at work since the establishment of the current plains surface. Veneer composition is dominated by the mixing of silicate-rich dust and sulfate-rich outcrop surface, but it has also been influenced by mineral precipitation, including NaCl, and possibly by limited physical or chemical weathering of sulfate minerals. Competing processes of chemical alteration (perhaps mediated by thin films of water or water vapor beneath blanketing soils) and sandblasting of exposed outcrop surfaces determine the current distribution of veneers. Dark-toned rinds several millimeters thick reflect more extensive surface alteration but also indicate combined dust admixture, halite precipitation, and possible minor sulfate removal. Cemented fracture fills that are differentially resistant to erosion occur along the margins of linear fracture systems possibly related to impact. These appear to reflect limited groundwater activity along the margins of fractures, cementing mechanically introduced fill derived principally from outcrop rocks. The limited thickness and spatial distribution of these three features suggest that aqueous activity has been rare and transient or has operated at exceedingly low rates during the protracted interval since outcropping Meridiani strata were exposed on the plains surface.

Journal of Geophysical Research E: Planets↗

USGS2010

Estimating the timing and location of shallow rainfall-induced landslides using a model for transient, unsaturated infiltration

Shallow rainfall-induced landslides commonly occur under conditions of transient infiltration into initially unsaturated soils. In an effort to predict the timing and location of such landslides, we developed a model of the infiltration process using a two-layer system that consists of an unsaturated zone above a saturated zone and implemented this model in a geographic information system (GIS) framework. The model links analytical solutions for transient, unsaturated, vertical infiltration above the water table to pressure-diffusion solutions for pressure changes below the water table. The solutions are coupled through a transient water table that rises as water accumulates at the base of the unsaturated zone. This scheme, though limited to simplified soil-water characteristics and moist initial conditions, greatly improves computational efficiency over numerical models in spatially distributed modeling applications. Pore pressures computed by these coupled models are subsequently used in one-dimensional slope-stability computations to estimate the timing and locations of slope failures. Applied over a digital landscape near Seattle, Washington, for an hourly rainfall history known to trigger shallow landslides, the model computes a factor of safety for each grid cell at any time during a rainstorm. The unsaturated layer attenuates and delays the rainfall-induced pore-pressure response of the model at depth, consistent with observations at an instrumented hillside near Edmonds, Washington. This attenuation results in realistic estimates of timing for the onset of slope instability (7 h earlier than observed landslides, on average). By considering the spatial distribution of physical properties, the model predicts the primary source areas of landslides.

Journal of Geophysical Research F: Earth Surface↗

USGS2023

Subsurface porewater flow accelerates talik development under the Alaska Highway, Yukon: A prelude to road collapse and permafrost thaw?

The presence of taliks (perennially unfrozen zones in permafrost areas) adversely affects the thermal stability of infrastructure in cold regions, including roads. The role of heat advection on talik development and feedback on permafrost degradation has not been quantified methodically in this context. We incorporate a surface energy balance model into a coupled groundwater flow and energy transport numerical model (SUTRA-ice). The model, calibrated with long-term observations (1997–2018 on the Alaska Highway), is used to investigate and quantify the role of heat advection on talik initiation and development under a road embankment. Over the 25-year simulation period, the new model is driven by reconstructed meteorological data and has a good agreement with near surface soil temperatures. The model successfully reproduces the increasing depth to the permafrost table (mean absolute error <0.2 m), and talik development. The results demonstrate that heat advection provides an additional energy source that expedites the rate of permafrost thaw and roughly doubles the rate of permafrost table deepening, compared to purely conductive thawing. Talik initially formed and grew over time under the combined effect of water flow, snow insulation, road construction and climate warming. Talik formation creates a new thermal state under the road embankment, resulting in acceleration of underlying permafrost degradation, due to the positive feedback of heat accumulation created by trapped unfrozen water. In a changing climate, mobile water flow will play a more important role in permafrost thaw and talik development under road embankments, and is likely to significantly increase maintenance costs and reduce the long-term stability of the infrastructure.

Yukon↗

USGS1996

Identifying coseismic subsidence in tidal-wetland stratigraphic sequences at the Cascadia subduction zone of western North America

Tidal-wetland stratigraphy reveals that great plate boundary earthquakes have caused hundreds of kilometers of coast to subside at the Cascadia subduction zone. However, determining earthquake recurrence intervals and mapping the coastal extent of past great earthquake ruptures in this region are complicated by the effects of many sedimentologic, hydrographic, and oceanographic processes that occur on the coasts of tectonically passive as well as active continental margins. Tidal-wetland stratigraphy at many Cascadia estuaries differs little from that at similar sites on passive-margin coasts where stratigraphic sequences form through nonseismic processes unrelated to coseismic land level changes. Methods developed through study of similar stratigraphic sequences in Europe provide a framework for investigating the Cascadia estuarine record. Five kinds of criteria must be evaluated when inferring regional coastal subsidence due to great plate boundary earthquakes: the suddenness and amount of submergence, the lateral extent of submerged tidal-wetland soils, the coincidence of submergence with tsunami deposits, and the degree of synchroneity of submergence events at widely spaced sites. Evaluation of such criteria at the Cascadia subduction zone indicates regional coastal subsidence during at least two great earthquakes. Evidence for a coseismic origin remains equivocal, however, for the many peat-mud contacts in Cascadia stratigraphic sequences that lack (1) contrasts in lithology or fossils indicative of more than half a meter of submergence, (2) well-studied tsunami deposits, or (3) precise ages needed for regional correlation. Paleoecologic studies of fossil assemblages are particularly important in estimating the size of sudden sea level changes recorded by abrupt peat-mud contacts and in helping to distinguish erosional and gradually formed contacts from coseismic contacts. Reconstruction of a history of great earthquakes for the Cascadia subduction zone will require rigorous application of the above criteria and many detailed investigations.

Journal of Geophysical Research B: Solid Earth↗

USGS2005

Reconstructing a mid-Cretaceous landscape from paleosols in western Canada

The Albian Stage of the mid-Cretaceous was a time of equable climate conditions with high sea levels and broad shallow epeiric seas that may have had a moderating affect on continental climates. A Late Albian landscape surface that developed during a regression and subsequent sea-level rise in the Western Canada Foreland Basin is reconstructed on the basis of correlation of paleosols penetrated by cores through the Paddy Member of the Peace River Formation. Reconstruction of this landscape refines chronostratigraphic relationships and will benefit future paleoclimatological studies milizing continental sphaerosiderite proxy records. The paleosols developed in estuarine sandstones and mudstones, and they exhibit evidence of a polygenetic history. Upon initial exposure and pedogenesis, the Paddy Member developed deeply weathered, well-drained cumulative soil profiles. Later stages of pedogenesis were characterized by hydromorphic soil conditions. The stages of soil development interpreted for the Paddy Member correlate with inferred stages of pedogenic development in time-equivalent formations located both basinward and downslope (upper Viking Formation), and landward and upslope (Boulder Creek Formation). On the basis of the genetic similarity among paleosols in these three correlative formations, the paleosols are interpreted as having formed along a single, continuous landscape surface. Results of this study indicate that the catena concept of pedogenesis along sloping landscapes is applicable to ancient successions. Sphaerosiderites in the Paddy Mem ber paleosols are used to provide proxy values for meteoric ??18O values at 52?? N paleolatitude in the Cretaceous Western Interior Basin. The meteoric ??18O values are used to refine existing interpretations about the mid-Cretaceous paleolatitudinal gradient in meteoric ?? 18O values, and the mid-Cretaceous hydrologic cycle. Copyright ?? 2005, SEPM (Society for Sedimentary Geology).

Journal of Sedimentary Research↗