Search USGS⌕ Search

SEARCH · Search USGS

Results for “Journal of Geophysical Research - Planets”

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

Fault interaction and stress triggering of twentieth century earthquakes in Mongolia

A cluster of exceptionally large earthquakes in the interior of Asia occurred from 1905 to 1967: the 1905 M7.9 Tsetserleg and M8.4 Bolnai earthquakes, the 1931 M8.0 Fu Yun earthquake, the 1957 M8.1 Gobi-Altai earthquake, and the 1967 M7.1 Mogod earthquake (sequence). Each of the larger (M ??? 8) earthquakes involved strike-slip faulting averaging more than 5 m and rupture lengths of several hundred kilometers. Available geologic data indicate that recurrence intervals on the major source faults are several thousands of years and distances of about 400 km separate the respective rupture areas. We propose that the occurrences of these and many smaller earthquakes are related and controlled to a large extent by stress changes generated by the compounded static deformation of the preceding earthquakes and subsequent viscoelastic relaxation of the lower crust and upper mantle beneath Mongolia. We employ a spherically layered viscoelastic model constrained by the 1994-2002 GPS velocity field in western Mongolia [Vergnolle et al., 2003]. Using the succession of twentieth century earthquakes as sources of deformation, we then analyze the time-dependent change in Coulomb failure stress (????f). At remote interaction distances, static ????f values are small. However, modeled postseismic stress changes typically accumulate to several tenths of a bar over time intervals of decades. Almost all significant twentieth century regional earthquakes (M ??? 6) with well-constrained fault geometry lie in positive ????f lobes of magnitude about +0.5 bar. Our results suggest that significant stress transfer is possible among continental faults separated by hundreds of kilometers and on timescales of decades. Copyright 2003 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2003

The 2000 Nemuro-Hanto-Oki earthquake, off eastern Hokkaido, Japan, and the high intraslab seismic activity in the southwestern Kuril Trench

The 2000 Nemuro-Hanto-Oki earthquake (Mw6.8) occurred in the southwestern part of the Kuril Trench. The hypocenter was located close to the aftershock region of the great 1994 Kuril earthquake (Mw8.3), named "the 1994 Hokkaido-Toho-Oki earthquake" by the Japan Meteorological Agency, for which the fault plane is still in debate. Analysis of the 2000 event provides a clue to resolve the fault plane issue for the 1994 event. The hypocenters of the 2000 main shock and aftershocks are determined using arrival times from a combination of nearby inland and submarine seismic networks with an improved azimuthal coverage. They clearly show that the 2000 event was an intraslab event occurring on a shallow-dipping fault plane between 55 and 65 km in depth. The well-focused aftershock distribution of the 2000 event, the relative location of the 1994 event with respect to the 2000 event, and the similarity between their focal mechanisms strongly suggest that the faulting of the great 1994 earthquake also occurred on a shallow-dipping fault plane in the subducting slab. The recent hypocenter distribution around the 1994 aftershock region also supports this result. Large intraslab earthquakes occuring to the southeast of Hokkaido may occur due to a strong coupling on the plate boundary, which generates relatively large stress field within the subducting Pacific plate.

Journal of Geophysical Research B: Solid Earth↗
USGS2003

Seismic velocity anisotropy and heterogeneity beneath the Mantle Electromagnetic and Tomography Experiment (MELT) region of the East Pacific Rise from analysis of P and S body waves

We use teleseismic P and S delay times and shear wave splitting measurements to constrain isotropic and anisotropic heterogeneity in the mantle beneath the southern East Pacific Rise (SEPR). The data comprise 462 P and S delay times and 18 shear wave splitting observations recorded during the Mantle Electromagnetic and Tomography (MELT) Experiment. We estimate the mantle melt content (F) and temperature (T) variation from the isotropic velocity variation. Our results indicate that the maximum variation in F beneath our array is between zero and ???1.2%, and maximum variation in T is between zero and ???100 K. We favor an explanation having partial contributions from both T and F. We approximate the seismic anisotropy of the upper mantle with hexagonal symmetry, consistent with the assumption of two dimensionality of mantle flow. Our new tomographic technique uses a nonlinear inversion of P and slow S polarization delay times to simultaneously solve for coupled VP and VS heterogeneity throughout the model and for the magnitude of anisotropy within discrete domains. The domain dimensions and the dip of the anisotropy are fixed for each inversion but are varied in a grid search, obtaining the misfit of the models to the body wave delay data and to split times of vertically propagating S waves. The data misfit and the isotropic heterogeneity are sensitive to domain dimensions and dip of anisotropy. In a region centered beneath the SEPR the best average dip of the hexagonal symmetry axis is horizontal or dipping shallowly (<30??) west. Given the resolution of our data, a subaxial region characterized by vertically aligned symmetry axes may exist but is limited to be <80 km deep. We infer that the mantle flow beneath the SEPR is consistent with shallow asthenospheric return flow from the direction of the South Pacific superswell.

Journal of Geophysical Research B: Solid Earth↗
USGS2003

The strength and rheology of methane clathrate hydrate

Methane clathrate hydrate (structure I) is found to be very strong, based on laboratory triaxial deformation experiments we have carried out on samples of synthetic, high-purity, polycrystalline material. Samples were deformed in compressional creep tests (i.e., constant applied stress, ??), at conditions of confining pressure P = 50 and 100 MPa, strain rate 4.5 ?? 10-8 ??? ?? ??? 4.3 ?? 10-4 s-1, temperature 260 ??? T ??? 287 K, and internal methane pressure 10 ??? PCH4 ??? 15 MPa. At steady state, typically reached in a few percent strain, methane hydrate exhibited strength that was far higher than expected on the basis of published work. In terms of the standard high-temperature creep law, ?? = A??ne-(E*+PV*)/RT the rheology is described by the constants A = 108.55 MPa-n s-1, n = 2.2, E* = 90,000 J mol-1, and V* = 19 cm3 mol-1. For comparison at temperatures just below the ice point, methane hydrate at a given strain rate is over 20 times stronger than ice, and the contrast increases at lower temperatures. The possible occurrence of syntectonic dissociation of methane hydrate to methane plus free water in these experiments suggests that the high strength measured here may be only a lower bound. On Earth, high strength in hydrate-bearing formations implies higher energy release upon decomposition and subsequent failure. In the outer solar system, if Titan has a 100-km-thick near-surface layer of high-strength, low-thermal conductivity methane hydrate as has been suggested, its interior is likely to be considerably warmer than previously expected.

Journal of Geophysical Research B: Solid Earth↗
USGS2003

Slip distribution of the 1952 Tokachi-Oki earthquake (M 8.1) along the Kuril Trench deduced from tsunami waveform inversion

We inverted 13 tsunami waveforms recorded in Japan to estimate the slip distribution of the 1952 Tokachi-Oki earthquake (M 8.1), which occurred southeast off Hokkaido along the southern Kuril subduction zone. The previously estimated source area determined from tsunami travel times [Hatori, 1973] did not coincide with the observed aftershock distribution. Our results show that a large amount of slip occurred in the aftershock area east of Hatori's tsunami source area, suggesting that a portion of the interplate thrust near the trench was ruptured by the main shock. We also found more than 5 m of slip along the deeper part of the seismogenic interface, just below the central part of Hatori's tsunami source area. This region, which also has the largest stress drop during the main shock, had few aftershocks. Large tsunami heights on the eastern Hokkaido coast are better explained by the heterogeneous slip model than previous uniform-slip fault models. The total seismic moment is estimated to be 1.87 ?? 1021 N m, giving a moment magnitude of Mw = 8.1. The revised tsunami source area is estimated to be 25.2 ?? 103 km2, ???3 times larger than the previous tsunami source area. Out of four large earthquakes with M ??? 7 that subsequently occurred in and around the rupture area of the 1952 event, three were at the edges of regions with relatively small amount of slip. We also found that a subducted seamount near the edge of the rupture area possibly impeded slip along the plate interface.

Journal of Geophysical Research B: Solid Earth↗
USGS2003

Shallow soil CO2 flow along the San Andreas and Calaveras Faults, California

We evaluate a comprehensive soil CO2 survey along the San Andreas fault (SAF) in Parkfield, and the Calaveras fault (CF) in Hollister, California, in the context of spatial and temporal variability, origin, and transport of CO2 in fractured terrain. CO2 efflux was measured within grids with portable instrumentation and continously with meteorological parameters at a fixed station, in both faulted and unfaulted areas. Spatial and temporal variability of surface CO2 effluxes was observed to be higher at faulted SAF and CF sites, relative to comparable background areas. However, ??13C (-23.3 to - 16.4???) and ??14C (75.5 to 94.4???) values of soil CO2 in both faulted and unfaulted areas are indicative of biogenic CO2, even though CO2 effluxes in faulted areas reached values as high as 428 g m-2 d-1. Profiles of soil CO2 concentration as a function of depth were measured at multiple sites within SAF and CF grids and repeatedly at two locations at the SAF grid. Many of these profiles suggest a surprisingly high component of advective CO2 flow. Spectral and correlation analysis of SAF CO2 efflux and meteorological parameter time series indicates that effects of wind speed variations on atmospheric air flow though fractures modulate surface efflux of biogenic CO2. The resulting areal patterns in CO2 effluxes could be erroneously attributed to a deep gas source in the absence of isotopic data, a problem that must be addressed in fault zone soil gas studies.

Journal of Geophysical Research B: Solid Earth↗
USGS2003

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

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

Journal of Geophysical Research C: Oceans↗
USGS2003

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

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

Journal of Geophysical Research B: Solid Earth↗
USGS2003

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

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

Journal of Geophysical Research B: Solid Earth↗
USGS2003

Tsunamis generated by subaerial mass flows

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

Journal of Geophysical Research B: Solid Earth↗
USGS2003

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

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

Journal of Geophysical Research B: Solid Earth↗
USGS2003

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

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

Journal of Geophysical Research B: Solid Earth↗
USGS2003

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

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

Journal of Geophysical Research B: Solid Earth↗
USGS2003

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

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

Journal of Geophysical Research B: Solid Earth↗
USGS2003

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

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

Journal of Geophysical Research B: Solid Earth↗
USGS2014

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

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

New York↗
USGS1991

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

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

Hawaii↗
USGS1999

Carbon cycling in boreal wetlands: A comparison of three approaches

Three independent methods were used to measure net ecosystem production (NEP) in four wetlands near Thompson, Manitoba, Canada. The first method calculated NEP by subtracting heterotrophic respiration from net primary productivity, using both measurements and estimates derived from the literature. The second method used radiocarbon data from cores to derive long-term NEP averaged over the past several decades. The third method used direct measurement of NEP combined with a model to fill in for days with no data. The three methods, with their independently derived uncertainties, all show the same magnitude and pattern of NEP variation across four different wetland types. However, direct measurement yielded distinctly lower estimates of NEP in the most productive sites. Highest NEP (31–180 gC m −2 yr −1 ) was observed in the two wetlands with the highest proportion of sedge vegetation. A bog collapse scar and a nutrient-rich fen had NEP values not significantly different from zero. The maximum NEP at sites with intermediate nutrient status is due to slower overall decomposition and is likely associated with greater allocation of production below ground by sedges. The three methods for estimating NEP differ in the effort required, the sources of error, and in the timescale over which they apply. Used in combination, they allow estimation of parameters such as below- ground production and the contribution of heterotrophic decomposition to total soil respiration. Using the radiocarbon method, we also derived estimates of the rate of N accumulation in the four wetland types.

Manitoba↗