USGS Science⌕ Search

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.

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

At least 1,225 records · Page 68Linked to original sources

Complex mean circulation over the inner shelf south of Martha's Vineyard revealed by observations and a high-resolution model

Inner-shelf circulation is governed by the interaction between tides, baroclinic forcing, winds, waves, and frictional losses; the mean circulation ultimately governs exchange between the coast and ocean. In some cases, oscillatory tidal currents interact with bathymetric features to generate a tidally rectified flow. Recent observational and modeling efforts in an overlapping domain centered on the Martha's Vineyard Coastal Observatory (MVCO) provided an opportunity to investigate the spatial and temporal complexity of circulation on the inner shelf. ADCP and surface radar observations revealed a mean circulation pattern that was highly variable in the alongshore and cross-shore directions. Nested modeling incrementally improved representation of the mean circulation as grid resolution increased and indicated tidal rectification as the generation mechanism of a counter-clockwise gyre near the MVCO. The loss of model skill with decreasing resolution is attributed to insufficient representation of the bathymetric gradients (Δh/h), which is important for representing nonlinear interactions between currents and bathymetry. The modeled momentum balance was characterized by large spatial variability of the pressure gradient and horizontal advection terms over short distances, suggesting that observed inner-shelf momentum balances may be confounded. Given the available observational and modeling data, this work defines the spatially variable mean circulation and its formation mechanism—tidal rectification—and illustrates the importance of model resolution for resolving circulation and constituent exchange near the coast. The results of this study have implications for future observational and modeling studies near the MVCO and other inner-shelf locations with alongshore bathymetric variability.

Massachusetts↗

Wave-current interaction in Willapa Bay

This paper describes the importance of wave-current interaction in an inlet-estuary system. The three-dimensional, fully coupled, Coupled Ocean-Atmosphere-Wave-Sediment Transport (COAWST) modeling system was applied in Willapa Bay (Washington State) from 22 to 29 October 1998 that included a large storm event. To represent the interaction between waves and currents, the vortex-force method was used. Model results were compared with water elevations, currents, and wave measurements obtained by the U.S. Army Corp of Engineers. In general, a good agreement between field data and computed results was achieved, although some discrepancies were also observed in regard to wave peak directions in the most upstream station. Several numerical experiments that considered different forcing terms were run in order to identify the effects of each wind, tide, and wave-current interaction process. Comparison of the horizontal momentum balances results identified that wave-breaking-induced acceleration is one of the leading terms in the inlet area. The enhancement of the apparent bed roughness caused by waves also affected the values and distribution of the bottom shear stress. The pressure gradient showed significant changes with respect to the pure tidal case. During storm conditions the momentum balance in the inlet shares the characteristics of tidal-dominated and wave-dominated surf zone environments. The changes in the momentum balance caused by waves were manifested both in water level and current variations. The most relevant effect on hydrodynamics was a wave-induced setup in the inner part of the estuary.

Washington↗

Utilizing multichannel electrical resistivity methods to examine the dynamics of the fresh water–seawater interface in two Hawaiian groundwater systems

Multichannel electrical resistivity (ER) measurements were conducted at two contrasting coastal sites in Hawaii to obtain new information on the spatial scales and dynamics of the fresh water–seawater interface and rates of coastal groundwater exchange. At Kiholo Bay (located on the dry, Kona side of the Big Island) and at a site in Maunalua Bay (Oahu), there is an evidence for abundant submarine groundwater discharge (SGD). However, the hydrologic and geologic controls on coastal groundwater discharge are likely to be different at these two sites. While at Kiholo Bay SGD is predominantly through lava tubes, at the Maunalua Bay site exchange occurs mostly through nearshore submarine springs. In order to calculate SGD fluxes, it is important to understand the spatial and temporal scales of coastal groundwater exchange. From ER time series data, subsurface salinity distributions were calculated using site-specific formation factors. A salinity mass balance box model was then used to calculate rates of point source (i.e., spatially discreet) and total fresh water discharge. From these data, mean SGD rates were calculated for Kiholo Bay (∼9,200 m 3 /d) and for the Maunalua Bay site (∼5,900 m 3 /d). While such results are on the same order of magnitude to geochemical tracer-derived SGD rates, the ER SGD rates provide enhanced details of coastal groundwater exchange that can enable a more cohesive whole watershed perspective.

Hawai'i↗

Use of vertical temperature gradients for prediction of tidal flat sediment characteristics

Sediment characteristics largely govern tidal flat morphologic evolution; however, conventional methods of investigating spatial variability in lithology on tidal flats are difficult to employ in these highly dynamic regions. In response, a series of laboratory experiments was designed to investigate the use of temperature diffusion toward sediment characterization. A vertical thermistor array was used to quantify temperature gradients in simulated tidal flat sediments of varying compositions. Thermal conductivity estimates derived from these arrays were similar to measurements from a standard heated needle probe, which substantiates the thermistor methodology. While the thermal diffusivities of dry homogeneous sediments were similar, diffusivities for saturated homogeneous sediments ranged approximately one order of magnitude. The thermal diffusivity of saturated sand was five times the thermal diffusivity of saturated kaolin and more than eight times the thermal diffusivity of saturated bentonite. This suggests that vertical temperature gradients can be used for distinguishing homogeneous saturated sands from homogeneous saturated clays and perhaps even between homogeneous saturated clay types. However, experiments with more realistic tidal flat mixtures were less discriminating. Relationships between thermal diffusivity and percent fines for saturated mixtures varied depending upon clay composition, indicating that clay hydration and/or water content controls thermal gradients. Furthermore, existing models for the bulk conductivity of sediment mixtures were improved only through the use of calibrated estimates of homogeneous end-member conductivity and water content values. Our findings suggest that remotely sensed observations of water content and thermal diffusivity could only be used to qualitatively estimate tidal flat sediment characteristics.

Journal of Geophysical Research C: Oceans↗

The Moon Mineralogy Mapper (M3) imaging spectrometer for lunar science: Instrument description, calibration, on‐orbit measurements, science data calibration and on‐orbit validation

The NASA Discovery Moon Mineralogy Mapper imaging spectrometer was selected to pursue a wide range of science objectives requiring measurement of composition at fine spatial scales over the full lunar surface. To pursue these objectives, a broad spectral range imaging spectrometer with high uniformity and high signal-to-noise ratio capable of measuring compositionally diagnostic spectral absorption features from a wide variety of known and possible lunar materials was required. For this purpose the Moon Mineralogy Mapper imaging spectrometer was designed and developed that measures the spectral range from 430 to 3000 nm with 10 nm spectral sampling through a 24 degree field of view with 0.7 milliradian spatial sampling. The instrument has a signal-to-noise ratio of greater than 400 for the specified equatorial reference radiance and greater than 100 for the polar reference radiance. The spectral cross-track uniformity is >90% and spectral instantaneous field-of-view uniformity is >90%. The Moon Mineralogy Mapper was launched on Chandrayaan-1 on the 22nd of October. On the 18th of November 2008 the Moon Mineralogy Mapper was turned on and collected a first light data set within 24 h. During this early checkout period and throughout the mission the spacecraft thermal environment and orbital parameters varied more than expected and placed operational and data quality constraints on the measurements. On the 29th of August 2009, spacecraft communication was lost. Over the course of the flight mission 1542 downlinked data sets were acquired that provide coverage of more than 95% of the lunar surface. An end-to-end science data calibration system was developed and all measurements have been passed through this system and delivered to the Planetary Data System (PDS.NASA.GOV). An extensive effort has been undertaken by the science team to validate the Moon Mineralogy Mapper science measurements in the context of the mission objectives. A focused spectral, radiometric, spatial, and uniformity validation effort has been pursued with selected data sets including an Earth-view data set. With this effort an initial validation of the on-orbit performance of the imaging spectrometer has been achieved, including validation of the cross-track spectral uniformity and spectral instantaneous field of view uniformity. The Moon Mineralogy Mapper is the first imaging spectrometer to measure a data set of this kind at the Moon. These calibrated science measurements are being used to address the full set of science goals and objectives for this mission.

Journal of Geophysical Research E: Planets↗

Stability of Mg-sulfates at-10C and the rates of dehydration/rehydration processes under conditions relevant to Mars

We report the results of low temperature (−10°C) experiments on the stability fields and phase transition pathways of five hydrous Mg-sulfates. A low temperature form of MgSO 4 ·7H 2 O (LT-7w) was found to have a wide stability field that extends to low relative humidity (∼13% RH at −10°C). Using information on the timing of phase transitions, we extracted information on the reaction rates of five important dehydration and rehydration processes. We found that the temperature dependencies of rate constants for dehydration processes differ from those of rehydration, which reflect differences in reaction mechanisms. By extrapolating these rate constants versus T correlations into the T range relevant to Mars, we can evaluate the possibility of occurrence of specific processes and the presence of common Mg-sulfate species present on Mars in different periods and locations. We anticipate in a moderate obliquity period, starkeyite and LH-MgSO 4 ·H 2 O should be two common Mg-sulfates at the surface, another polymorph MH-MgSO 4 ·H 2 O can exist at the locations where hydrothermal processes may have occurred. In polar regions or within the subsurface of other regions, meridianiite (coexisting with water ice, near 100% RH) and LT-7w (over a large RH range) are the stable phases. During a high obliquity period, meridianiite and LT-7w should exhibit widespread occurrence. The correlations of reaction rates versus temperature found in this study imply that dehydration and rehydration of hydrous Mg-sulfates would always be slower than the sublimation and crystallization of water ice, which would be supported by mission observations from Odyssey and by Mars Exploration Rovers.

Journal of Geophysical Research E: Planets↗

Characterization of previously unidentified lunar pyroclastic deposits using Lunar Reconnaissance Orbiter Camera (LROC) data

We used a Lunar Reconnaissance Orbiter Camera (LROC) global monochrome Wide-angle Camera (WAC) mosaic to conduct a survey of the Moon to search for previously unidentified pyroclastic deposits. Promising locations were examined in detail using LROC multispectral WAC mosaics, high-resolution LROC Narrow Angle Camera (NAC) images, and Clementine multispectral (ultraviolet-visible or UVVIS) data. Out of 47 potential deposits chosen for closer examination, 12 were selected as probable newly identified pyroclastic deposits. Potential pyroclastic deposits were generally found in settings similar to previously identified deposits, including areas within or near mare deposits adjacent to highlands, within floor-fractured craters, and along fissures in mare deposits. However, a significant new finding is the discovery of localized pyroclastic deposits within floor-fractured craters Anderson E and F on the lunar farside, isolated from other known similar deposits. Our search confirms that most major regional and localized low-albedo pyroclastic deposits have been identified on the Moon down to ~100 m/pix resolution, and that additional newly identified deposits are likely to be either isolated small deposits or additional portions of discontinuous, patchy deposits.

Journal of Geophysical Research E: Planets↗

Advances in the simulation and automated measurement of well-sorted granular material: 1. Simulation

1. In this, the first of a pair of papers which address the simulation and automated measurement of well-sorted natural granular material, a method is presented for simulation of two-phase (solid, void) assemblages of discrete non-cohesive particles. The purpose is to have a flexible, yet computationally and theoretically simple, suite of tools with well constrained and well known statistical properties, in order to simulate realistic granular material as a discrete element model with realistic size and shape distributions, for a variety of purposes. The stochastic modeling framework is based on three-dimensional tessellations with variable degrees of order in particle-packing arrangement. Examples of sediments with a variety of particle size distributions and spatial variability in grain size are presented. The relationship between particle shape and porosity conforms to published data. The immediate application is testing new algorithms for automated measurements of particle properties (mean and standard deviation of particle sizes, and apparent porosity) from images of natural sediment, as detailed in the second of this pair of papers. The model could also prove useful for simulating specific depositional structures found in natural sediments, the result of physical alterations to packing and grain fabric, using discrete particle flow models. While the principal focus here is on naturally occurring sediment and sedimentary rock, the methods presented might also be useful for simulations of similar granular or cellular material encountered in engineering, industrial and life sciences.

Journal of Geophysical Research F: Earth Surface↗

Advances in the simulation and automated measurement of well-sorted granular material: 2. Direct measures of particle properties

1. In this, the second of a pair of papers on the structure of well-sorted natural granular material (sediment), new methods are described for automated measurements from images of sediment, of: 1) particle-size standard deviation (arithmetic sorting) with and without apparent void fraction; and 2) mean particle size in material with void fraction. A variety of simulations of granular material are used for testing purposes, in addition to images of natural sediment. Simulations are also used to establish that the effects on automated particle sizing of grains visible through the interstices of the grains at the very surface of a granular material continue to a depth of approximately 4 grain diameters and that this is independent of mean particle size. Ensemble root-mean squared error between observed and estimated arithmetic sorting coefficients for 262 images of natural silts, sands and gravels (drawn from 8 populations) is 31%, which reduces to 27% if adjusted for bias (slope correction between observed and estimated values). These methods allow non-intrusive and fully automated measurements of surfaces of unconsolidated granular material. With no tunable parameters or empirically derived coefficients, they should be broadly universal in appropriate applications. However, empirical corrections may need to be applied for the most accurate results. Finally, analytical formulas are derived for the one-step pore-particle transition probability matrix, estimated from the image's autocorrelogram, from which void fraction of a section of granular material can be estimated directly. This model gives excellent predictions of bulk void fraction yet imperfect predictions of pore-particle transitions.

Journal of Geophysical Research F: Earth Surface↗

Scaling the Teflon Peaks: Rock type and the generation of extreme relief in the glaciated western Alaska Range

Parts of the Alaska Range (Alaska, USA) stand in prominent exception to the “glacial buzzsaw hypothesis,” which postulates that terrain raised above the ELA is rapidly denuded by glaciers. In this paper, we discuss the role of a strong contrast in rock type in the development of this exceptional terrain. Much of the range is developed on pervasively fractured flysch, with local relief of 1000–1500 m, and mean summit elevations that are similar to modern snow line elevations. In contrast, Cretaceous and Tertiary plutons of relatively intact granite support the range's tallest mountains (including Mt. McKinley, or Denali, at 6194 m), with 2500–5000 m of local relief. The high granitic peaks protrude well above modern snow lines and support many large glaciers. We focus on the plutons of the Denali massif and the Kichatna Mountains, to the west. We use field observations, satellite photos, and digital elevation data to demonstrate how exhumation of these plutons affects glacier longitudinal profiles, the glacial drainage network, and the effectiveness of periglacial processes. In strong granite, steep, smooth valley walls are maintained by detachment of rock slabs along sheeting joints. These steep walls act as low-friction surfaces (“Teflon”), efficiently shedding snow. Simple scaling calculations show that this avalanching may greatly enhance the health of the modern glaciers. We conclude that, in places such as Denali, unusual combinations of rapid tectonic uplift and great rock strength have created the highest relief in North America by enhancing glacial erosion in the valleys while preserving the peaks.

Alaska↗

Analysis of low-frequency seismic signals generated during a multiple-iceberg calving event at Jakobshavn Isbræ, Greenland

We investigated seismic signals generated during a large-scale, multiple iceberg calving event that occurred at Jakobshavn Isbræ, Greenland, on 21 August 2009. The event was recorded by a high-rate time-lapse camera and five broadband seismic stations located within a few hundred kilometers of the terminus. During the event two full-glacier-thickness icebergs calved from the grounded (or nearly grounded) terminus and immediately capsized; the second iceberg to calve was two to three times smaller than the first. The individual calving and capsize events were well-correlated with the radiation of low-frequency seismic signals (<0.1 Hz) dominated by Love and Rayleigh waves. In agreement with regional records from previously published ‘glacial earthquakes’, these low-frequency seismic signals had maximum power and/or signal-to-noise ratios in the 0.05–0.1 Hz band. Similarly, full waveform inversions indicate that these signals were also generated by horizontal single forces acting at the glacier terminus. The signals therefore appear to be local manifestations of glacial earthquakes, although the magnitudes of the signals (twice-time integrated force histories) were considerably smaller than previously reported glacial earthquakes. We thus speculate that such earthquakes may be a common, if not pervasive, feature of all full-glacier-thickness calving events from grounded termini. Finally, a key result from our study is that waveform inversions performed on low-frequency, calving-generated seismic signals may have only limited ability to quantitatively estimate mass losses from calving. In particular, the choice of source time function has little impact on the inversion but dramatically changes the earthquake magnitude. Accordingly, in our analysis, it is unclear whether the smaller or larger of the two calving icebergs generated a larger seismic signal.

Journal of Geophysical Research↗

Displacement fields from point cloud data: Application of particle imaging velocimetry to landslide geodesy

Acquiring spatially continuous ground-surface displacement fields from Terrestrial Laser Scanners (TLS) will allow better understanding of the physical processes governing landslide motion at detailed spatial and temporal scales. Problems arise, however, when estimating continuous displacement fields from TLS point-clouds because reflecting points from sequential scans of moving ground are not defined uniquely, thus repeat TLS surveys typically do not track individual reflectors. Here, we implemented the cross-correlation-based Particle Image Velocimetry (PIV) method to derive a surface deformation field using TLS point-cloud data. We estimated associated errors using the shape of the cross-correlation function and tested the method's performance with synthetic displacements applied to a TLS point cloud. We applied the method to the toe of the episodically active Cleveland Corral Landslide in northern California using TLS data acquired in June 2005–January 2007 and January–May 2010. Estimated displacements ranged from decimeters to several meters and they agreed well with independent measurements at better than 9% root mean squared (RMS) error. For each of the time periods, the method provided a smooth, nearly continuous displacement field that coincides with independently mapped boundaries of the slide and permits further kinematic and mechanical inference. For the 2010 data set, for instance, the PIV-derived displacement field identified a diffuse zone of displacement that preceded by over a month the development of a new lateral shear zone. Additionally, the upslope and downslope displacement gradients delineated by the dense PIV field elucidated the non-rigid behavior of the slide.

California↗

Elementary theory of bed-sediment entrainment by debris flows and avalanches

Analyses of mass and momentum exchange between a debris flow or avalanche and an underlying sediment layer aid interpretations and predictions of bed-sediment entrainment rates. A preliminary analysis assesses the behavior of a Coulomb slide block that entrains bed material as it descends a uniform slope. The analysis demonstrates that the block's momentum can grow unstably, even in the presence of limited entrainment efficiency. A more-detailed, depth-integrated continuum analysis of interacting, deformable bodies identifies mechanical controls on entrainment efficiency, and shows that entrainment rates satisfy a jump condition that involves shear-traction and velocity discontinuities at the flow-bed boundary. Explicit predictions of the entrainment rate E result from making reasonable assumptions about flow velocity profiles and boundary shear tractions. For Coulomb-friction tractions, predicted entrainment rates are sensitive to pore fluid pressures that develop in bed sediment as it is overridden. In the simplest scenario the bed sediment liquefies completely, and the entrainment-rate equation reduces to E = 2 μ 1 gh 1 cos θ (1 − λ 1 )/ , where θ is the slope angle, μ 1 is the flow's Coulomb friction coefficient, h 1 is its thickness, λ 1 is its degree of liquefaction, and is its depth-averaged velocity. For values of λ 1 ranging from 0.5 to 0.8, this equation predicts entrainment rates consistent with rates of 0.05 to 0.1 m/s measured in large-scale debris-flow experiments in which wet sediment beds liquefied almost completely. The propensity for bed liquefaction depends on several factors, including sediment porosity, permeability, and thickness, and rates of compression and shear deformation that occur when beds are overridden.

Journal of Geophysical Research F: Earth Surface↗

Monitoring glacier surface seismicity in time and space using Rayleigh waves

Sliding glaciers and brittle ice failure generate seismic body and surface wave energy characteristic to the source mechanism. Here we analyze continuous seismic recordings from an array of nine short-period passive seismometers located on Bench Glacier, Alaska (USA) (61.033°N, 145.687°W). We focus on the arrival-time and amplitude information of the dominant Rayleigh wave phase. Over a 46-hour period we detect thousands of events using a cross-correlation based event identification method. Travel-time inversion of a subset of events (7% of the total) defines an active crevasse, propagating more than 200 meters in three hours. From the Rayleigh wave amplitudes, we estimate the amount of volumetric opening along the crevasse as well as an average bulk attenuation ( = 42) for the ice in this part of the glacier. With the remaining icequake signals we establish a diurnal periodicity in seismicity, indicating that surface run-off and subglacial water pressure changes likely control the triggering of these surface events. Furthermore, we find that these events are too weak (i.e., too noisy) to locate individually. However, stacking individual events increases the signal-to-noise ratio of the waveforms, implying that these periodic sources are effectively stationary during the recording period.

Journal of Geophysical Research F: Earth Surface↗

Sediment entrainment by debris flows: In situ measurements from the headwaters of a steep catchment

Debris flows can dramatically increase their volume, and hence their destructive potential, by entraining sediment. Yet quantitative constraints on rates and mechanics of sediment entrainment by debris flows are limited. Using an in situ sensor network in the headwaters of a natural catchment we measured flow and bed properties during six erosive debris-flow events. Despite similar flow properties and thicknesses of bed sediment entrained across all events, time-averaged entrainment rates were significantly faster for bed sediment that was saturated prior to flow arrival compared with rates for sediment that was dry. Bed sediment was entrained from the sediment-surface downward in a progressive fashion and occurred during passage of dense granular fronts as well as water-rich, inter-surge flow. En masse failure of bed sediment along the sediment-bedrock interface was never observed. Large-magnitude, high-frequency fluctuations in total normal basal stress were dissipated within the upper 5 cm of bed sediment. Within this near surface layer, concomitant fluctuations in Coulomb frictional resistance are expected, irrespective of the influence of pore fluid pressure or fluctuations in shear stress. If the near-surface sediment was wet as it was overridden by a flow, additional large-magnitude, high-frequency pore pressure fluctuations were measured in the near-surface bed sediment. These pore pressure fluctuations propagated to depth at subsonic rates and in a diffusive manner. The depth to which large excess pore pressures propagated was typically less than 10 cm, but scaled as ( D / f i ) 0.5 , in which D is the hydraulic diffusivity and f i is the frequency of a particular pore pressure fluctuation. Shallow penetration depths of granular-normal-stress fluctuations and excess pore pressures demonstrate that only near-surface bed sediment experiences the full dynamic range of effective-stress fluctuations, and as a result, can be more easily entrained than deeper sediment. These data provide robust tests for mechanical models of entrainment and demonstrate that a debris flow over wet bed sediment will be larger than the same flow over dry bed sediment.

Colorado↗

Probabilistic prediction of barrier-island response to hurricanes

Prediction of barrier-island response to hurricane attack is important for assessing the vulnerability of communities, infrastructure, habitat, and recreational assets to the impacts of storm surge, waves, and erosion. We have demonstrated that a conceptual model intended to make qualitative predictions of the type of beach response to storms (e.g., beach erosion, dune erosion, dune overwash, inundation) can be reformulated in a Bayesian network to make quantitative predictions of the morphologic response. In an application of this approach at Santa Rosa Island, FL, predicted dune-crest elevation changes in response to Hurricane Ivan explained about 20% to 30% of the observed variance. An extended Bayesian network based on the original conceptual model, which included dune elevations, storm surge, and swash, but with the addition of beach and dune widths as input variables, showed improved skill compared to the original model, explaining 70% of dune elevation change variance and about 60% of dune and shoreline position change variance. This probabilistic approach accurately represented prediction uncertainty (measured with the log likelihood ratio), and it outperformed the baseline prediction (i.e., the prior distribution based on the observations). Finally, sensitivity studies demonstrated that degrading the resolution of the Bayesian network or removing data from the calibration process reduced the skill of the predictions by 30% to 40%. The reduction in skill did not change conclusions regarding the relative importance of the input variables, and the extended model's skill always outperformed the original model.

Florida↗

Effects of river regulation on aeolian landscapes, Colorado River, southwestern USA

Connectivity between fluvial and aeolian sedimentary systems plays an important role in the physical and biological environment of dryland regions. This study examines the coupling between fluvial sand deposits and aeolian dune fields in bedrock canyons of the arid to semiarid Colorado River corridor, southwestern USA. By quantifying significant differences between aeolian landscapes with and without modern fluvial sediment sources, this work demonstrates for the first time that the flow- and sediment-limiting effects of dam operations affect sedimentary processes and ecosystems in aeolian landscapes above the fluvial high water line. Dune fields decoupled from fluvial sand supply have more ground cover (biologic crust and vegetation) and less aeolian sand transport than do dune fields that remain coupled to modern fluvial sand supply. The proportion of active aeolian sand area also is substantially lower in a heavily regulated river reach (Marble–Grand Canyon, Arizona) than in a much less regulated reach with otherwise similar environmental conditions (Cataract Canyon, Utah). The interconnections shown here among river flow and sediment, aeolian sand transport, and biologic communities in aeolian dunes demonstrate a newly recognized means by which anthropogenic influence alters dryland environments. Because fluvial–aeolian coupling is common globally, it is likely that similar sediment-transport connectivity and interaction with upland ecosystems are important in other dryland regions to a greater degree than has been recognized previously.

Colorado↗

Characterizing the performance of ecosystem models across time scales: A spectral analysis of the North American Carbon Program site‐level synthesis

[1] Ecosystem models are important tools for diagnosing the carbon cycle and projecting its behavior across space and time. Despite the fact that ecosystems respond to drivers at multiple time scales, most assessments of model performance do not discriminate different time scales. Spectral methods, such as wavelet analyses, present an alternative approach that enables the identification of the dominant time scales contributing to model performance in the frequency domain. In this study we used wavelet analyses to synthesize the performance of 21 ecosystem models at 9 eddy covariance towers as part of the North American Carbon Program's site‐level intercomparison. This study expands upon previous single‐site and single‐model analyses to determine what patterns of model error are consistent across a diverse range of models and sites. To assess the significance of model error at different time scales, a novel Monte Carlo approach was developed to incorporate flux observation error. Failing to account for observation error leads to a misidentification of the time scales that dominate model error. These analyses show that model error (1) is largest at the annual and 20–120 day scales, (2) has a clear peak at the diurnal scale, and (3) shows large variability among models in the 2–20 day scales. Errors at the annual scale were consistent across time, diurnal errors were predominantly during the growing season, and intermediate‐scale errors were largely event driven. Breaking spectra into discrete temporal bands revealed a significant model‐by‐band effect but also a nonsignificant model‐by‐site effect, which together suggest that individual models show consistency in their error patterns. Differences among models were related to model time step, soil hydrology, and the representation of photosynthesis and phenology but not the soil carbon or nitrogen cycles. These factors had the greatest impact on diurnal errors, were less important at annual scales, and had the least impact at intermediate time scales.

Journal of Geophysical Research: Biogeosciences↗