SEARCH · USGS Science

Results for “Contribution”

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,711 records · Page 95Linked to original sources

USGS2021

An eddy-resolving numerical model to study turbulent flow, sediment and bed evolution using detached eddy simulation in a lateral separation zone at the field-scale

Turbulence-resolving simulations elucidate key elements of fluid dynamics and sediment transport in fluvial environments. This research presents a feasible strategy for applying state-of-the-art computational fluid mechanics to the study of sediment transport and morphodynamic processes in lateral separation zones, which are common features in canyon rivers where massive lateral flow separation causes large-scale turbulence that controls sediment erosion and deposition. An eddy-resolving model was developed and tested at the field-scale, coupling a viscous flow and sediment transport solver using Detached Eddy Simulation techniques. A morphodynamic model was applied to the viscous flow/sediment solver to calculate erosion and deposition. A simulation of turbulence was performed at the grid resolution for a straight channel to determine the relative contributions of modeled and resolved diffusivity. The time-dependent, energetically important, correlative, non-stationary signals of the simulated quantities were captured at the lateral separation zone. Strong periodic signals featured by high amplitude were found at the separation zone, while low frequency pulsations were observed at the reattachment zone of the lateral separation zone. Interactions between the eddies and the loose bed boundaries resulted in erosion of sediment at the main channel followed by deposition at the primary eddy and eddy bars. tions elucidate key elements of fluid dynamics and sediment transport in fluvial environments. This research presents a feasible strategy for applying state-of-the-art computational fluid mechanics to the study of sediment transport and morphodynamic processes in lateral separation zones, which are common features in canyon rivers where massive lateral flow separation causes large-scale turbulence that controls sediment erosion and deposition. An eddy-resolving model was developed and tested at the field scale, coupling viscous flow and sediment transport solver using Detached Eddy Simulation (DES) techniques. A morphodynamic model was applied to the viscous flow/sediment solver to calculate erosion and deposition. A simulation of turbulence was performed at the grid resolution for a straight channel to determine the relative contributions of modeled and resolved diffusivity. The time-dependent, energetically important, correlative non-stationary signals of the simulated quantities were captured at the lateral separation zone. Strong periodic signals featured by high amplitude were found at the separation zone, while low frequency pulsations were observed at the reattachment zone of the lateral separation zone. Interactions between the eddies and the loose bed boundaries resulted in massive erosion of sediment at the main channel followed by deposition at the primary eddy and eddy bars.

Arizona↗

USGS2021

Spatiotemporal methane emission from global reservoirs

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

Journal of Geophysical Research: Biogeosciences↗

USGS2023

Inland water greenhouse gas budgets for RECCAP2: 2. Regionalization and homogenization of estimates

Inland waters are important sources of the greenhouse gasses (GHGs) carbon dioxide (CO 2 ), methane (CH 4 ) and nitrous oxide (N 2 O) to the atmosphere. In the framework of the 2 nd phase of the REgional Carbon Cycle Assessment and Processes (RECCAP-2) initiative, we synthesize existing estimates of GHG emissions from streams, rivers, lakes and reservoirs, and homogenize them with regard to underlying global maps of water surface area distribution and the effects of seasonal ice cover. We then produce regionalized estimates of GHG emissions over 10 extensive land regions. According to our synthesis, inland water GHG emissions have a global warming potential of an equivalent emission of 13.5 (9.9-20.1) and 8.3 (5.7-12.7) Pg CO 2 -eq. yr -1 at a 20 and 100 year horizon (GWP 20 and GWP 100 ), respectively. Contributions of CO 2 dominate GWP 100 , with rivers being the largest emitter. For GWP 20 , lakes and rivers are equally important emitters, and the warming potential of CH 4 is more important than that of CO 2 . Contributions from N 2 O are about two orders of magnitude lower. Normalized to the area of RECCAP-2 regions, S-America and SE-Asia show the highest emission rates, dominated by riverine CO 2 emissions.

Global Biogeochemical Cycles↗

USGS2022

On the use of high-resolution and deep-learning seismic catalogs for short-term earthquake forecasts: Potential benefits and current limitations

Enhanced earthquake catalogs provide detailed images of evolving seismic sequences. Currently, these data sets take some time to be released but will soon become available in real time. Here, we explore whether and how enhanced seismic catalogs feeding into established short-term earthquake forecasting protocols may result in higher predictive skill. We consider three enhanced catalogs for the 2016–2017 Central Italy sequence, featuring a bulk completeness lower by at least two magnitude units compared to the real-time catalog and an improved hypocentral resolution. We use them to inform a set of physical Coulomb Rate-and-State (CRS) and statistical Epidemic-Type Aftershock Sequence (ETAS) models to forecast the space-time occurrence of M3+ events during the first 6 months of the sequence. We track model performance using standard likelihood-based metrics and compare their skill against the best-performing CRS and ETAS models among those developed with the real-time catalog. We find that while the incorporation of the triggering contributions from new small magnitude detections of the enhanced catalogs is beneficial for both types of forecasts, these models do not significantly outperform their respective near real-time benchmarks. To explore the reasons behind this result, we perform targeted sensitivity tests that show how (a) the typical spatial discretizations of forecast experiments ( 2 km) hamper the ability of models to capture highly localized secondary triggering patterns and (b) differences in earthquake parameters (i.e., magnitude and hypocenters) reported in different catalogs can affect forecast evaluation. These findings will contribute toward improving forecast model design and evaluation strategies for next-generation seismic catalogs.

Journal of Geophysical Research--Solid Earth↗

USGS2023

Characterizing performance of freshwater wetland methane models across time scales at FLUXNET-CH4 sites using wavelet analyses

Process-based land surface models are important tools for estimating global wetland methane (CH 4 ) emissions and projecting their behavior across space and time. So far there are no performance assessments of model responses to drivers at multiple time scales. In this study, we apply wavelet analysis to identify the dominant time scales contributing to model uncertainty in the frequency domain. We evaluate seven wetland models at 23 eddy covariance tower sites. Our study first characterizes site-level patterns of freshwater wetland CH 4 fluxes (FCH 4 ) at different time scales. A Monte Carlo approach was developed to incorporate flux observation error to avoid misidentification of the time scales that dominate model error. Our results suggest that (a) significant model-observation disagreements are mainly at multi-day time scales (<15 days); (b) most of the models can capture the CH 4 variability at monthly and seasonal time scales (>32 days) for the boreal and Arctic tundra wetland sites but have significant bias in variability at seasonal time scales for temperate and tropical/subtropical sites; (c) model errors exhibit increasing power spectrum as time scale increases, indicating that biases at time scales <5 days could contribute to persistent systematic biases on longer time scales; and (d) differences in error pattern are related to model structure (e.g., proxy of CH 4 production). Our evaluation suggests the need to accurately replicate FCH 4 variability, especially at short time scales, in future wetland CH 4 model developments.

JGR Biogeosciences↗

USGS2023

Multi-year tracing of spatial and temporal dynamics of post-fire aeolian sediment transport using rare earth elements provide insights into grassland management

Aeolian sediment transport occurs as a function of, and with feedback to ecosystem changes and disturbances. Many desert grasslands are undergoing rapid changes in vegetation, including the encroachment of woody plants, which alters fire regimes and in turn can change the spatial and temporal patterns of aeolian sediment transport. We investigated aeolian sediment transport and spatial distribution of sediment in the surface soil for 7 years following a prescribed fire using a multiple rare earth element (REE) tracer-based approach in a shrub-encroached desert grassland in the northern Chihuahuan desert. Results indicate that even though the aeolian horizontal sediment mass flux increased approximately three-fold in the first windy season in the burned areas compared to control areas, there were no significant differences after three windy seasons. The soil surface of bare microsites was the major contributor of aeolian sediments in unburned areas (87%), while the shrub microsites contributed the least (<2%) during the observation period. However, after the prescribed fire, the contribution of aeolian sediments from shrub microsites increased considerably (∼40%), indicating post-fire microsite-scale sediment redistribution. The findings of this study, which is the first to use multiple REE tracers for multi-year analysis of the spatial and temporal dynamics of aeolian sediment transport, illustrate how disturbance by prescribed fire can influence aeolian processes and alters dryland soil geomorphology in which distinct soils develop over time at very fine spatial scales of individual plants.

New Mexico↗

USGS2024

Quantifying compound and nonlinear effects of hurricane-induced flooding using a dynamically coupled hydrological-ocean model

We recently developed a dynamically coupled hydrological-ocean modeling system that provides seamless coverage across the land-ocean continuum during hurricane-induced compound flooding. This study introduced a local inertial equation and a diagonal flow algorithm to the overland routing of the coupled system’s hydrology model (WRF-Hydro). Using Hurricane Florence (2018) as a test case, the performance of the coupled model was significantly improved, evidenced by its enhanced capability of capturing backwater and increased water level simulation accuracy and stability. With four model experiments, we present a framework to detangle, define, and quantify compound and nonlinear effects. The results revealed that the flood peaks in the lower Cape Fear River Basin and the coastal waters were contributed by inland flooding and storm surge, respectively. These two processes had comparable contributions to the flooding in the Cape Fear River Estuary. The compound effect was identified when the flood levels resulting from the combination of land and ocean processes surpassed those caused by an individual process alone. The compound effect during Hurricane Florence exhibited limited impact on flood peaks, primarily due to the time lag between the peaks of the storm surge and the inland flooding. In the period between the two peaks, the compound effect was salient and significantly impacted the magnitude and variation of the flood level. The nonlinear effect, defined as the difference between the compound flood level and the superposition of storm surge and inland flooding water levels, reduced flood levels in the river channels while increasing flood levels on the floodplain.

Water Resources Research↗

USGS2025

Light absorbing particles deposited to snow cover across the Upper Colorado River Basin, Colorado Rocky Mountains, 2013-16: Interannual variations from multiple natural and anthropogenic sources

Atmospheric particulate matter (PM) as light-absorbing particles (LAPs) deposited to snow cover can result in early onset and rapid snow melting, challenging management of downstream water resources. We identified LAPs in 38 snow samples (water years 2013–2016) from the mountainous Upper Colorado River basin by comparing among laboratory-measured spectral reflectance, chemical, physical, and magnetic properties. Dust sample reflectance, averaged over the wavelength range of 0.35–2.50 μm, varied by a factor of 1.9 (range, 0.2300–0.4444) and was suppressed mainly by three components: (a) carbonaceous matter measured as total organic carbon (1.6–22.5 wt. %) including inferred black carbon, natural organic matter, and carbon-based synthetic, black road-tire-wear particles, (b) dark rock and mineral particles, indicated by amounts of magnetite (0.11–0.37 wt. %) as their proxy, and (c) ferric oxide minerals identified by reflectance spectroscopy and magnetic properties. Fundamental compositional differences were associated with different iron oxide groups defined by dominant hematite, goethite, or magnetite. These differences in iron oxide mineralogy are attributed to temporally varying source-area contributions implying strong interannual changes in regional source behavior, dust-storm frequency, and (or) transport tracks. Observations of dust-storm activity in the western U.S. and particle-size averages for all samples (median, 25 μm) indicated that regional dust from deserts dominated mineral-dust masses. Fugitive contaminants, nevertheless, contributed important amounts of LAPs from many types of anthropogenic sources.

Colorado↗

USGS2025

Organic matter availability and the production of methane-derived dissolved organic carbon at methane seeps

Methane seeps located along continental margins and slopes export sedimentary methane and dissolved organic carbon (DOC) into the ocean. The flux of these reduced carbon molecules from the seafloor into the ocean impacts ocean chemistry and supports deep-sea life. While significant effort has been made to understand how the anaerobic oxidation of methane (AOM) regulates the release of methane from the seafloor, little is known about the production of DOC in association with AOM or its flux and fate in the ocean. We hypothesize a mechanism for methane incorporation into DOC at seeps and investigate the relationship between sediment total organic carbon (TOC) availability and the incorporation of methane-derived carbon into DOC at four methane seep regions along the Cascadia margin, with a range of microbial and thermogenic methane sources. At sites with <2.0 wt.% TOC (Hydrate Ridge and Bullseye Vent), up to 60%–80% of carbon in DOC is methane-carbon, much more than sites with >2.0 wt.% TOC (Astoria Canyon and Barkley Canyon). We attribute the greater methane contribution at the more TOC-limited sites to a greater role of AOM in the carbon cycle, whereas at the organic matter-rich sites, microbial competition for sulfate as an electron acceptor for organic matter decomposition limits AOM and hence the transfer of carbon from methane to DOC. We estimate that the global diffusive flux of methane-derived DOC from the seafloor is 0.07–10.1 Tg C/yr, contributing to the stock of DOC present in the deep ocean and/or fueling the deep-sea microbial loop.

California, Oregon, Washington↗

USGS2025

Wetlands, groundwater and seasonality influence the spatial distribution of stream chemistry in a low-relief catchment

Evaluating stream water chemistry patterns provides insight into catchment ecosystem and hydrologic processes. Spatially distributed patterns and controls of stream solutes are well-established for high-relief catchments where solute flow paths align with surface topography. However, the controls on solute patterns are poorly constrained for low-relief catchments where hydrogeologic heterogeneities and river corridor features, like wetlands, may influence water and solute transport. Here, we provide a data set of solute patterns from 58 synoptic surveys across 28 sites and over 32 months in a low-relief wetland-rich catchment to determine the major surface and subsurface controls along with wetland influence across the catchment. In this low-relief catchment, the expected wetland storage, processing, and transport of solutes is only apparent in solute patterns of the smallest subcatchments. Meanwhile, downstream seasonal and wetland influence on observed chemistry can be masked by large groundwater contributions to the main stream channel. These findings highlight the importance of incorporating variable groundwater contributions into catchment-scale studies for low-relief catchments, and that understanding the overall influence of wetlands on stream chemistry requires sampling across various spatial and temporal scales. Therefore, in low-relief wetland-rich catchments, given the mosaic of above and below ground controls on stream solutes, modeling efforts may need to include both surface and subsurface hydrological data and processes.

Michigan↗

USGS2012

A model for the generation of two-dimensional surf beat

A finite difference model predicting group-forced long waves in the nearshore is constructed with two interacting parts: an incident wave model providing time-varying radiation stress gradients across the nearshore, and a long-wave model which solves the equations of motion for the forcing imposed by the incident waves. Both shallow water group-bound long waves and long waves generated by a time-varying breakpoint are simulated. Model-generated time series are used to calculate the cross correlation between wave groups and long waves through the surf zone. The cross-correlation signal first observed by Tucker (1950) is well predicted. For the first time, this signal is decomposed into the contributions from the two mechanisms of leaky mode forcing. Results show that the cross-correlation signal can be explained by bound long waves which are amplified, though strongly modified, through the surf zone before reflection from the shoreline. The breakpoint-forced long waves are added to the bound long waves at a phase of pi/2 and are a secondary contribution owing to their relatively small size.

Journal of Geophysical Research C: Oceans↗

USGS2008

Preliminary analysis of the hydrologic and geochemical controls on acid‐neutralizing capacity in two acidic seepage lakes in Florida

In late 1988, parallel studies of Lake Five‐O ( p H 5.14) in the Florida panhandle and Lake Barco ( p H 4.50) in north central Florida were initiated to develop hydrologic and major ion budgets of these lakes as part of an overall effort to improve understanding of the hydrologic, depositional, and biogeochemical factors that control acid‐neutralizing capacity (ANC) in seepage lakes. Preliminary findings from these studies indicate that earlier perceptions of lake hydrology and mechanisms of ANC regulation in Florida seepage lakes may have to be revised. The traditional perspective of seepage lakes in the Florida panhandle views these systems as dominated by precipitation inputs and that ANC regulation is due largely to in‐lake processes. Our results for Lake Five‐O show modest to steep hydraulic gradients almost entirely around the lake. In addition, the horizontal hydraulic conductivity of the surficial aquifer is high (8–74 m day −1 ), indicating that large quantities of groundwater flow into Lake Five‐O. Calculations of net groundwater flow from hydrologic budgets also indicate that groundwater may comprise at least 38 to 46% of the total inflow. For Lake Barco, net flow estimates of the minimum groundwater inflow range from 5 to 14% of total inflow. Enrichment factor and ion flux calculations for Lake Five‐O and Lake Barco indicate that terrestrial as well as in‐lake processes contribute significantly to ANC regulation. The extent that terrestrial processes contribute to ANC generation is directly related to the magnitude of groundwater inflow as well as the degree of ion enrichment or depletion that occurs in the surficial aquifer. Net ANC generation in both study lakes was dominated by anion retention (NO 3 − and SO 4 2− ). Where previous studies concluded that in‐lake reduction was the primary sink for SO 4 2− , our preliminary calculations show that adsorption of SO 4 2− within the watershed is perhaps twice as important as in‐lake reduction as a source of ANC. Net base cation enrichment in both lakes was negligible.

Florida↗

USGS1997

Variation of the 3-μm absorption feature on Mars: observations over eastern Valles Marineris by the mariner 6 infrared spectrometer

A new approach for calibration of the shortest wavelength channel (1.8 to 6.0 μm) of the Mariner 6 infrared spectrometer was derived. This calibration provides a new description of the instrument response function from 1.8 to 3.7 μm and accounts for the thermal contribution to the signal at longer wavelengths. This allows the two segments from 1.8 to 6 μm to be merged into a single spectrum. The broad water of hydration absorption spans these two segments and is examined in these merged spectra using a method of band integration. Unlike previous analyses which rely on ratios at two wavelengths, the integration method can assess the band strength independently from the albedo in the near infrared. Spectra taken over the eastern end of the Valles Marineris are examined for variations of the band-integrated value, and three distinct clusters are found. Within the estimated uncertainty, two clusters (both low and high albedo) have approximately the same integrated band depth. The third cluster (medium albedo) has an integrated band depth about 10% higher. This difference cannot be systematically attributed to either surface or atmospheric parameters and suggests variation in the amount of water either chemically or physically bound in surface materials. Approximately one-half of the high integrated band depth cluster is associated with chaotic terrain at the source of outflow channels, the other half occurs over lower inertia plains adjacent to chasmata. This suggests both surface physical properties and mineralogy as well as water in exchange with the atmosphere contribute to the 3-μm bound water absorption.

Journal of Geophysical Research E: Planets↗

USGS1996

Metastable mantle phase transformations and deep earthquakes in subducting oceanic lithosphere

Earth's deepest earthquakes occur as a population in subducting or previously subducted lithosphere at depths ranging from about 325 to 690 km. This depth interval closely brackets the mantle transition zone, characterized by rapid seismic velocity increases resulting from the transformation of upper mantle minerals to higher-pressure phases. Deep earthquakes thus provide the primary direct evidence for subduction of the lithosphere to these depths and allow us to investigate the deep thermal, thermodynamic, and mechanical ferment inside slabs. Numerical simulations of reaction rates show that the olivine → spinel transformation should be kinetically hindered in old, cold slabs descending into the transition zone. Thus wedge-shaped zones of metastable peridotite probably persist to depths of more than 600 km. Laboratory deformation experiments on some metastable minerals display a shear instability called transformational faulting. This instability involves sudden failure by localized superplasticity in thin shear zones where the metastable host mineral transforms to a denser, finer-grained phase. Hence in cold slabs, such faulting is expected for the polymorphic reactions in which olivine transforms to the spinel structure and clinoenstatite transforms to ilmenite. It is thus natural to hypothesize that deep earthquakes result from transformational faulting in metastable peridotite wedges within cold slabs. This consideration of the mineralogical states of slabs augments the traditional largely thermal view of slab processes and explains some previously enigmatic slab features. It explains why deep seismicity occurs only in the approximate depth range of the mantle transition zone, where minerals in downgoing slabs should transform to spinel and ilmenite structures. The onset of deep shocks at about 325 km is consistent with the onset of metastability near the equilibrium phase boundary in the slab. Even if a slab penetrates into the lower mantle, earthquakes should cease at depths near 700 km, because the seismogenic phase transformations in the slab are completed or can no longer occur. Substantial metastability is expected only in old, cold slabs, consistent with the observed restriction of deep earthquakes to those settings. Earthquakes should be restricted to the cold cores of slabs, as in any model in which the seismicity is temperature controlled, via the distribution of metastability. However, the geometries of recent large deep earthquakes pose a challenge for any such models. Transformational faulting may give insight into why deep shocks lack appreciable aftershocks and why their source characteristics, including focal mechanisms indicating localized shear failure rather than implosive deformation, are so similar to those of shallow earthquakes. Finally, metastable phase changes in slabs would produce an internal source of stress in addition to those due to the weight of the sinking slab. Such internal stresses may explain the occurrence of earthquakes in portions of lithosphere which have foundered to the bottom of the transition zone and/or are detached from subducting slabs. Metastability in downgoing slabs could have considerable geodynamic significance. Metastable wedges would reduce the negative buoyancy of slabs, decrease the driving force for subduction, and influence the state of stress in slabs. Heat released by metastable phase changes would raise temperatures within slabs and facilitate the transformation of spinel to the lower mantle mineral assemblage, causing slabs to equilibrate more rapidly with the ambient mantle and thus contribute to the cessation of deep seismicity. Because wedge formation should occur only for fast subducting slabs, it may act as a “parachute” and contribute to regulating plate speeds. Wedge formation would also have consequences for mantle evolution because the density of a slab stagnated near the bottom of the transition zone would increase as it heats up and the wedge transforms to denser spinel, favoring the subsequent sinking of the slab into the lower mantle.

Reviews of Geophysics↗

USGS1997

The roles of time and displacement in velocity-dependent volumetric strain of fault zones

The relationship between measured friction μ A and volumetric strain during frictional sliding was determined using a rate and state variable dependent friction constitutive equation, a common work balance relating friction and volume change, and two types of experimental faults: initially bare surfaces of Westerly granite and rock surfaces separated by a 1 mm layer of less than 90 μm Westerly granite gouge. The constitutive equation is the sum of a constant term representing the nominal resistance to sliding and two smaller terms: a rate dependent term representing the shear viscosity of the fault surface (direct effect), and a term which represents variations in the area of contact (evolution effect). The work balance relationship requires that μ A differs from the frictional resistance that leads to shear heating by the derivative of fault normal displacement with respect shear displacement, d δ n / d δ s . An implication of this relationship is that the rate dependence of d δ n / d δ s contributes to the rate dependence of μ A . Experiments show changes in sliding velocity lead to changes in both fault strength and volume. Analysis of data with the rate and state equations combined with the work balance relationship preclude the conventional interpretation of the direct effect in the rate and state variable constitutive equations. Consideration of a model bare surface fault consisting of an undeformable indentor sliding on a deformable surface reveals a serious flaw in the work balance relationship if volume change is time-dependent. For the model, at zero slip rate indentation creep under the normal load leads to time-dependent strengthening of the fault surface but, according to the work balance relationship, no work is done because compaction or dilatancy can only be induced by shearing. Additional tests on initially bare surfaces and gouges show that fault normal strain in experiments is time-dependent, consistent with the model. This time-dependent fault normal strain, which is not accounted for in the work balance relationship, explains the inconsistency between the constitutive equations and the work balance. For initially bare surface faults, all rate dependence of volume change is due to time dependence. Similar results are found for gouge. We conclude that μ A reflects the frictional resistance that results in shear heating, and no correction needs to be made for the volume changes. The result that time-dependent volume changes do not contribute to μ A is a general result and extends beyond these experiments, the simple indentor model and particular constitutive equations used to illustrate the principle.

Journal of Geophysical Research B: Solid Earth↗

USGS1997

Comparing nocturnal eddy covariance measurements to estimates of ecosystem respiration made by scaling chamber measurements at six coniferous boreal sites

During the growing season, nighttime ecosystem respiration emits 30–100% of the daytime net photosynthetic uptake of carbon, and therefore measurements of rates and understanding of its control by the environment are important for understanding net ecosystem exchange. Ecosystem respiration can be measured at night by eddy covariance methods, but the data may not be reliable because of low turbulence or other methodological problems. We used relationships between woody tissue, foliage, and soil respiration rates and temperature, with temperature records collected on site to estimate ecosystem respiration rates at six coniferous BOREAS sites at half-hour or 1-hour intervals, and then compared these estimates to nocturnal measurements of CO 2 exchange by eddy covariance. Soil surface respiration was the largest source of CO 2 at all sites (48–71%), and foliar respiration made a large contribution to ecosystem respiration at all sites (25–43%). Woody tissue respiration contributed only 5–15% to ecosystem respiration. We estimated error for the scaled chamber predictions of ecosystem respiration by using the uncertainty associated with each respiration parameter and respiring biomass value. There was substantial uncertainty in estimates of foliar and soil respiration because of the spatial variability of specific respiration rates. In addition, more attention needs to be paid to estimating foliar respiration during the early part of the growing season, when new foliage is growing, and to determining seasonal trends of soil surface respiration. Nocturnal eddy covariance measurements were poorly correlated to scaled chamber estimates of ecosystem respiration ( r 2 =0.06–0.27) and were consistently lower than scaled chamber predictions (by 27% on average for the six sites). The bias in eddy covariance estimates of ecosystem respiration will alter estimates of gross assimilation in the light and of net ecosystem exchange rates over extended periods.

Journal of Geophysical Research D: Atmospheres↗

USGS1997

Case study modeling of turbulent and mesoscale fluxes over the BOREAS region

Results from aircraft and surface observations provided evidence for the existence of mesoscale circulations over the Boreal Ecosystem-Atmosphere Study (BOREAS) domain. Using an integrated approach that included the use of analytical modeling, numerical modeling, and data analysis, we have found that there are substantial contributions to the total budgets of heat over the BOREAS domain generated by mesoscale circulations. This effect is largest when the synoptic flow is relatively weak, yet it is present under less favorable conditions, as shown by the case study presented here. While further analysis is warranted to document this effect, the existence of mesoscale flow is not surprising, since it is related to the presence of landscape patches, including lakes, which are of a size on the order of the local Rossby radius and which have spatial differences in maximum sensible heat flux of about 300 W m-2. We have also analyzed the vertical temperature profile simulated in our case study as well as high-resolution soundings and we have found vertical profiles of temperature change above the boundary layer height, which we attribute in part to mesoscale contributions. Our conclusion is that in regions with organized landscapes, such as BOREAS, even with relatively strong synoptic winds, dynamical scaling criteria should be used to assess whether mesoscale effects should be parameterized or explicitly resolved in numerical models of the atmosphere.

Journal of Geophysical Research D: Atmospheres↗

USGS2012

Heat flow and energetics of the San Andreas fault zone

Approximately 100 heat flow measurements in the San Andreas fault zone indicate (1) there is no evidence for local factional heating of the main fault trace at any latitude over a 1000-km length from Cape Mendocino to San Bernardino, (2) average heat flow is high (∼2 HFU, ∼80 mW m −2 ) throughout the 550-km segment of the Coast Ranges that encloses the San Andreas fault zone in central California; this broad anomaly falls off rapidly toward the Great Valley to the east, and over a 200-km distance toward the Mendocino Triple Junction to the northwest. As others have pointed out, a local conductive heat flow anomaly would be detectable unless the frictional resistance allocated to heat production on the main trace were ≲100 bars. Frictional work allocated to surface energy of new fractures is probably unimportant, and hydrologic convection is not likely to invalidate the conduction assumption, since the heat discharge by thermal springs near the fault is negligible. Explanations for the low dynamic friction fall into two intergradational classes: those in which the fault is weak all of the time and those in which it is weak only during earthquakes (possibly just large ones). The first class includes faults containing anomalously weak gouge materials and faults containing materials with normal frictional properties under near-lithostatic steady state fluid pressures. In the second class, weakening is caused by the event (for example, a thermally induced increase in fluid pressure, dehydration of clay minerals, or acoustic fluidization). In this class, unlike the first, the average strength and ambient tectonic shear stress may be large, ∼1 kbar, but the stress allocated to elastic radiation (the apparent stress) must be of similar magnitude, an apparent contradiction with seismic estimates. Unless seismic radiation is underestimated for large earthquakes, it is difficult to justify average tectonic stresses on the main trace of the San Andreas fault in excess of ∼200 bars. The development of the broad Coast Range heat flow anomaly southward from Cape Mendocino suggests that heat flow increases by a factor of 2 within 4 m.y. after the passage of the Mendocino Triple Junction. This passage leaves the San Andreas transform fault zone in its wake; the depth of the anomalous sources cannot be much greater than the depth of the seismogenic layer. Some of the anomalous heat may be supplied by conduction from the warmer mantle that must occur south of the Mendocino transform (where there is no subducting slab), and some might be supplied by shear heating in the fault zone. With no contribution from shear heating, extreme mantle upwelling would be required, and asthenosphere conditions should exist today at depths of only ∼20 km in the northernmost Coast Ranges. If there is an appreciable contribution from shear heating, the heat flow constraint implies that the seismogenic layer is partially decoupled at its base and that the basal traction is in the sense that resists right lateral motion on the fault(s). As a result of these basal tractions, the average shearing stress in the seismogenic layer would increase with distance from the main fault, and the seismogenic layer would offer substantial resistance to plate motion even though resistance on the main fault might be negligible. These speculative models have testable consequences.

Journal of Geophysical Research Solid Earth↗