Search USGS⌕ Search

SEARCH · Search USGS

Results for “Journal of Geophysical Research - Oceans”

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 613 records · Page 34Linked to original sources

Wind-driven particle mobility on Mars: Insights from Mars Exploration Rover observations at "El Dorado" and surroundings at Gusev Crater

The ripple field known as 'El Dorado' was a unique stop on Spirit's traverse where dust-raising, active mafic sand ripples and larger inactive coarse-grained ripples interact, illuminating several long-standing issues of Martian dust mobility, sand mobility, and the origin of transverse aeolian ridges. Strong regional wind events endured by Spirit caused perceptible migration of ripple crests in deposits SSE of El Dorado, erasure of tracks in sandy areas, and changes to dust mantling the site. Localized thermal vortices swept across El Dorado, leaving paths of reduced dust but without perceptibly damaging nearly cohesionless sandy ripple crests. From orbit, winds responsible for frequently raising clay-sized dust into the atmosphere do not seem to significantly affect dunes composed of (more easily entrained) sand-sized particles, a long-standing paradox. This disparity between dust mobilization and sand mobilization on Mars is due largely to two factors: (1) dust occurs on the surface as fragile, low-density, sand-sized aggregates that are easily entrained and disrupted, compared with clay-sized air fall particles; and (2) induration of regolith is pervasive. Light-toned bed forms investigated at Gusev are coarse-grained ripples, an interpretation we propose for many of the smallest linear, light-toned bed forms of uncertain origin seen in high-resolution orbital images across Mars. On Earth, wind can organize bimodal or poorly sorted loose sediment into coarse-grained ripples. Coarse-grained ripples could be relatively common on Mars because development of durable, well-sorted sediments analogous to terrestrial aeolian quartz sand deposits is restricted by the lack of free quartz and limited hydraulic sediment processing. Copyright 2008 by the American Geophysical Union.

Journal of Geophysical Research E: Planets↗

Light-toned salty soils and co-existing Si-rich species discovered by the Mars Exploration Rover Spirit in Columbia Hills

Light-toned soils were exposed, through serendipitous excavations by Spirit Rover wheels, at eight locations in the Columbia Hills. Their occurrences were grouped into four types on the basis of geomorphic settings. At three major exposures, the light-toned soils are hydrous and sulfate-rich. The spatial distributions of distinct types of salty soils vary substantially: with centimeter-scaled heterogeneities at Paso Robles, Dead Sea, Shredded, and Champagne-Penny, a well-mixed nature for light-toned soils occurring near and at the summit of Husband Hill, and relatively homogeneous distributions in the two layers at the Tyrone site. Aeolian, fumarolic, and hydrothermal fluid processes are suggested to be responsible for the deposition, transportation, and accumulation of these light-toned soils. In addition, a change in Pancam spectra of Tyrone yellowish soils was observed after being exposed to current Martian surface conditions for 175 sols. This change is interpreted to be caused by the dehydration of ferric sulfates on the basis of laboratory simulations and suggests a relative humidity gradient beneath the surface. Si-rich nodules and soils were observed near the major exposures of S-rich soils. They possess a characteristic feature in Pancam visible near-infrared (Vis-NIR) spectra that may be diagnostic of hydrated species, and this spectral feature can be used to search for additional Si-rich species. The exposures of hydrated salty soils within various geomorphic settings imply the potential existence of hydrous minerals in similar settings over a much wider area. Hydrous sulfates represent one of the candidates that may contribute the high level of water equivalent hydrogen in equatorial regions detected by the Neutron Spectrometer on Mars Odyssey.

Journal of Geophysical Research↗

Degradation of Victoria crater, Mars

The ∼750 m diameter and ∼75 m deep Victoria crater in Meridiani Planum, Mars, is a degraded primary impact structure retaining a ∼5 m raised rim consisting of 1–2 m of uplifted rocks overlain by ∼3 m of ejecta at the rim crest. The rim is 120–220 m wide and is surrounded by a dark annulus reaching an average of 590 m beyond the raised rim. Comparison between observed morphology and that expected for pristine craters 500–750 m across indicates that the original, pristine crater was close to 600 m in diameter. Hence, the crater has been erosionally widened by ∼150 m and infilled by ∼50 m of sediments. Eolian processes are responsible for most crater modification, but lesser mass wasting or gully activity contributions cannot be ruled out. Erosion by prevailing winds is most significant along the exposed rim and upper walls and accounts for ∼50 m widening across a WNW–ESE diameter. The volume of material eroded from the crater walls and rim is ∼20% less than the volume of sediments partially filling the crater, indicating eolian infilling from sources outside the crater over time. The annulus formed when ∼1 m deflation of the ejecta created a lag of more resistant hematite spherules that trapped <10–20 cm of darker, regional basaltic sands. Greater relief along the rim enabled meters of erosion. Comparison between Victoria and regional craters leads to definition of a crater degradation sequence dominated by eolian erosion and infilling over time.

Journal of Geophysical Research E: Planets↗

Relative age of interior layered deposits in southwest Candor Chasma based on high-resolution structural mapping

High‐resolution topography generated from stereo HiRISE (High‐Resolution Imaging Science Experiment) imagery reveals the meter‐scale structure of interior layered deposits (ILD) in southwest Candor Chasma. This study seeks to determine the age of the local ILD relative to any normal faults that can be attributed to chasma formation. The study area is located near the contact of these ILD and the wall rock and is in an area where chasma‐forming normal faults have been proposed. We find that while normal faults are found in the study area, these faults are not sufficiently large nor appropriately located or oriented to accommodate the roughly northeast‐southwest extension that is required for normal faults that can be attributed to chasma formation. Additionally, bedding exposed in the local ILD generally dips toward the center of Candor Chasma, consistent with sediment deposition in a preexisting basin. Further, pit craters of Tithonia Catena, presumed to predate or be contemporaneous with the formation of west Candor Chasma, do not cut into the ILD within the study area. These independent lines of evidence support a postchasma age for the ILD exposed within the study area. Chasma‐related normal faults may exist within these ILD at depth but are not exposed at the surface. Approximately 2 km of conformable stratigraphy is exposed in the study area, and therefore at least several kilometers of the local ILD were deposited subsequent to any chasma‐related normal faulting that may have occurred in this part of Candor Chasma.

Journal of Geophysical Research E: Planets↗

Spirit Mars Rover Mission to the Columbia Hills, Gusev Crater: Mission overview and selected results from the Cumberland Ridge to Home Plate

This paper summarizes the Spirit rover operations in the Columbia Hills of Gusev Crater from sols 513 to 1476 and provides an overview of selected findings that focus on synergistic use of the Athena Payload and comparisons to orbital data. Results include discovery of outcrops (Voltaire) on Husband Hill that are interpreted to be altered impact melt deposits that incorporated local materials during emplacement. Evidence for extensive volcanic activity and aqueous alteration in the Inner Basin is also detailed, including discovery and characterization of accretionary lapilli and formation of sulfate, silica, and hematite‐rich deposits. Use of Spirit's data to understand the range of spectral signatures observed over the Columbia Hills by the Mars Reconnaissance Orbiter's Compact Reconnaissance Imaging Spectrometer (CRISM) hyperspectral imager (0.4–4 μ m) is summarized. We show that CRISM spectra are controlled by the proportion of ferric‐rich dust to ferrous‐bearing igneous minerals exposed in ripples and other wind‐blown deposits. The evidence for aqueous alteration derived from Spirit's data is associated with outcrops that are too small to be detected from orbital observations or with materials exposed from the shallow subsurface during rover activities. Although orbital observations show many other locations on Mars with evidence for minerals formed or altered in an aqueous environment, Spirit's data imply that the older crust of Mars has been altered even more extensively than evident from orbital data. This result greatly increases the potential that the surface or shallow subsurface was once a habitable regime.

Journal of Geophysical Research E: Planets↗

Thorium abundances on the Aristarchus plateau: Insights into the composition of the Aristarchus pyroclastic glass deposits

Thorium (Th) data from the Lunar Prospector gamma ray spectrometer (LP‐GRS) are used to constrain the composition of lunar pyroclastic glass deposits on top of the Aristarchus plateau. Our goal is to use forward modeling of LP‐GRS Th data to measure the Th abundances on the plateau and then to determine if the elevated Th abundances on the plateau are associated with the pyroclastic deposits or with thorium‐rich ejecta from Aristarchus crater. We use a variety of remote sensing data to show that there is a large, homogenous portion of the pyroclastics on the plateau that has seen little or no contamination from the Th‐rich ejecta of Aristarchus crater. Our results show that the uncontaminated pyroclastic glasses on Aristarchus plateau have an average Th content of 6.7 ppm and ∼7 wt % TiO 2 . These Th and Ti values are consistent with Th‐rich, intermediate‐Ti yellow glasses from the lunar sample suite. On the basis of this information, we use petrologic equations and interelement correlations for the Moon to estimate the composition of the source region from which the Aristarchus glasses were derived. We find that the source region for the Aristarchus glasses contained high abundances of heat‐producing elements, which most likely served as a thermal driver for the prolonged volcanic activity in this region of the Moon.

Journal of Geophysical Research E: Planets↗

Observations of periglacial landforms in Utopia Planitia with the High Resolution Imaging Science Experiment (HiRISE)

The region of western Utopia Planitia (80-105 degreesE, 40-55 degrees N) displays several types of landforms similar to Earth periglacial features, including scallop-shaped depressions and networks of polygonal terrains. The scalloped depressions have been proposed to originate from thermokarstic processes such as sublimation and/or melting of near-surface ground ice. Using HiRISE imagery, we characterize these depressions and several associated, distinct polygon networks in unprecedented morphologic and topographic detail and investigate support for an ice-based degradation process. The scalloped depressions and interior polygons and ridges are found to evolve together, mainly influenced by sublimation, local proximity of ground ice to the surface, and obliquity variations.

Journal of Geophysical Research↗

Constraints on deep moonquake focal mechanisms through analyses of tidal stress

[1] A relationship between deep moonquake occurrence and tidal forcing is suggested by the monthly periodicities observed in the occurrence times of events recorded by the Apollo Passive Seismic Experiment. In addition, the typically large S wave to P wave arrival amplitude ratios observed on deep moonquake seismograms are indicative of shear failure. Tidal stress, induced in the lunar interior by the gravitational influence of the Earth, may influence moonquake activity. We investigate the relationship between tidal stress and deep moonquake occurrence by searching for a linear combination of the normal and shear components of tidal stress that best approximates a constant value when evaluated at the times of moonquakes from 39 different moonquake clusters. We perform a grid search at each cluster location, computing the stresses resolved onto a suite of possible failure planes, to obtain the best fitting fault orientation at each location. We find that while linear combinations of stresses (and in some cases stress rates) can fit moonquake occurrence at many clusters quite well; for other clusters, the fit is not strongly dependent on plane orientation. This suggests that deep moonquakes may occur in response to factors other than, or in addition to, tidal stress. Several of our inferences support the hypothesis that deep moonquakes might be related to transformational faulting, in which shear failure is induced by mineral phase changes at depth. The occurrence of this process would have important implications for the lunar interior. Copyright 2009 by the American Geophysical Union.

Journal of Geophysical Research E: Planets↗

Moon meteoritic seismic hum: Steady state prediction

We use three different statistical models describing the frequency of meteoroid impacts on Earth to estimate the seismic background noise due to impacts on the lunar surface. Because of diffraction, seismic events on the Moon are typically characterized by long codas, lasting 1 h or more. We find that the small but frequent impacts generate seismic signals whose codas overlap in time, resulting in a permanent seismic noise that we term the "lunar hum" by analogy with the Earth's continuous seismic background seismic hum. We find that the Apollo era impact detection rates and amplitudes are well explained by a model that parameterizes (1) the net seismic impulse due to the impactor and resulting ejecta and (2) the effects of diffraction and attenuation. The formulation permits the calculation of a composite waveform at any point on the Moon due to simulated impacts at any epicentral distance. The root-mean-square amplitude of this waveform yields a background noise level that is about 100 times lower than the resolution of the Apollo long-period seismometers. At 2 s periods, this noise level is more than 1000 times lower than the low noise model prediction for Earth's microseismic noise. Sufficiently sensitive seismometers will allow the future detection of several impacts per day at body wave frequencies. Copyright 2009 by the American Geophysical Union.

Journal of Geophysical Research E: Planets↗

Characterization of phyllosilicates observed in the central Mawrth Vallis region, Mars, their potential formational processes, and implications for past climate

Mawrth Vallis contains one of the largest exposures of phyllosilicates on Mars. Nontronite, montmorillonite, kaolinite, and hydrated silica have been identified throughout the region using data from the Compact Reconnaissance Imaging Spectrometer for Mars (CRISM). In addition, saponite has been identified in one observation within a crater. These individual minerals are identified and distinguished by features at 1.38-1.42, ???1.91, and 2.17-2.41 ??m. There are two main phyllosilicate units in the Mawrth Vallis region. The lowermost unit is nontronite bearing, unconformably overlain by an Al-phyllosilicate unit containing montmorillonite plus hydrated silica, with a thin layer of kaolinite plus hydrated silica at the top of the unit. These two units are draped by a spectrally unremarkable capping unit. Smectites generally form in neutral to alkaline environments, while kaolinite and hydrated silica typically form in slightly acidic conditions; thus, the observed phyllosilicates may reflect a change in aqueous chemistry. Spectra retrieved near the boundary between the nontronite and Al-phyllosilicate units exhibit a strong positive slope from 1 to 2 ??m, likely from a ferrous component within the rock. This ferrous component indicates either rapid deposition in an oxidizing environment or reducing conditions. Formation of each of the phyllosilicate minerals identified requires liquid water, thus indicating a regional wet period in the Noachian when these units formed. The two main phyllosilicate units may be extensive layers of altered volcanic ash. Other potential formational processes include sediment deposition into a marine or lacustrine basin or pedogenesis. Copyright 2009 by the American Geophysical Union.

Journal of Geophysical Research E: Planets↗

Simulating and understanding sand wave variation: A case study of the Golden Gate sand waves

In this paper we present a detailed comparison between measured features of the Golden Gate sand wave field and the results of a nonlinear sand wave model. Because the Golden Gate sand waves exhibit large variation in their characteristics and in their environmental physics, this area gives us the opportunity to study sand wave variation between locations, within one well-measured, large area. The nonlinear model used in this paper is presently the only tool that provides information on the nonlinear evolution of large-amplitude sand waves. The model is used to increase our understanding of the coupling between the variability in environmental conditions and the sand wave characteristics. Results show that the model is able to describe the variation in the Golden Gate sand waves well when both the local oscillating tidal current and the residual current are taken into account. Current and water depth seem to be the most important factors influencing sand wave characteristics. The simulation results give further confidence in the underlying model hypothesis and assumptions. Copyright 2009 by the American Geophysical Union.

Journal of Geophysical Research F: Earth Surface↗

Compositions of modern dust and surface sediments in the Desert Southwest, United States

Modern dusts across southwestern United States deserts are compositionally similar to dust-rich Av soil horizons (depths of 0-0.5 cm and 1-4 cm at 35 sites) for common crustal elements but distinctly different for some trace elements. Chemical compositions and magnetic properties of the soil samples are similar among sites relative to dust sources, geographic areas, and lithologic substrates. Exceptions are Li, U, and W, enriched in Owens Valley, California, and Mg and Sr, enriched in soils formed on calcareous fan gravel in southeast Nevada. The Av horizons are dominated by dust and reflect limited mixing with substrate sediments. Modern dust samples are also similar across the region, except that Owens Valley dusts are higher in Mg, Ba, and Li and dusts both there and at sites to the north on volcanic substrates are higher in Sb and W. Thus, dust and Av horizons consist of contributions from many different sources that are well mixed before deposition. Modern dusts contain significantly greater amounts of As, Cd, Cr, Cu, Ni, Pb, and Sb than do Av horizons, which record dust additions over hundreds to thousands of years. These results suggest that modern dust compositions are influenced by anthropogenic sources and emissions from Owens (dry) Lake after its artificial desiccation in 1926. Both modern dusts and Av horizons are enriched in As, Ba, Cu, Li, Sb, Th, U, and W relative to average crustal composition, which we interpret to indicate that the geologic sources of dust in the southwestern United States are geochemically distinctive.

Journal of Geophysical Research F: Earth Surface↗

Instrumental record of debris flow initiation during natural rainfall: Implications for modeling slope stability

The middle of a hillslope hollow in the Oregon Coast Range failed and mobilized as a debris flow during heavy rainfall in November 1996. Automated pressure transducers recorded high spatial variability of pore water pressure within the area that mobilized as a debris flow, which initiated where local upward flow from bedrock developed into overlying colluvium. Postfailure observations of the bedrock surface exposed in the debris flow scar reveal a strong spatial correspondence between elevated piezometric response and water discharging from bedrock fractures. Measurements of apparent root cohesion on the basal ( C b ) and lateral ( C l ) scarp demonstrate substantial local variability, with areally weighted values of C b = 0.1 and C l = 4.6 kPa. Using measured soil properties and basal root strength, the widely used infinite slope model, employed assuming slope parallel groundwater flow, provides a poor prediction of hydrologic conditions at failure. In contrast, a model including lateral root strength (but neglecting lateral frictional strength) gave a predicted critical value of relative soil saturation that fell within the range defined by the arithmetic and geometric mean values at the time of failure. The 3‐D slope stability model CLARA‐W, used with locally observed pore water pressure, predicted small areas with lower factors of safety within the overall slide mass at sites consistent with field observations of where the failure initiated. This highly variable and localized nature of small areas of high pore pressure that can trigger slope failure means, however, that substantial uncertainty appears inevitable for estimating hydrologic conditions within incipient debris flows under natural conditions.

Oregon↗

Hindcasting of decadal‐timescale estuarine bathymetric change with a tidal‐timescale model

Hindcasting decadal-timescale bathymetric change in estuaries is prone to error due to limited data for initial conditions, boundary forcing, and calibration; computational limitations further hinder efforts. We developed and calibrated a tidal-timescale model to bathymetric change in Suisun Bay, California, over the 1867–1887 period. A general, multiple-timescale calibration ensured robustness over all timescales; two input reduction methods, the morphological hydrograph and the morphological acceleration factor, were applied at the decadal timescale. The model was calibrated to net bathymetric change in the entire basin; average error for bathymetric change over individual depth ranges was 37%. On a model cell-by-cell basis, performance for spatial amplitude correlation was poor over the majority of the domain, though spatial phase correlation was better, with 61% of the domain correctly indicated as erosional or depositional. Poor agreement was likely caused by the specification of initial bed composition, which was unknown during the 1867–1887 period. Cross-sectional bathymetric change between channels and flats, driven primarily by wind wave resuspension, was modeled with higher skill than longitudinal change, which is driven in part by gravitational circulation. The accelerated response of depth may have prevented gravitational circulation from being represented properly. As performance criteria became more stringent in a spatial sense, the error of the model increased. While these methods are useful for estimating basin-scale sedimentation changes, they may not be suitable for predicting specific locations of erosion or deposition. They do, however, provide a foundation for realistic estuarine geomorphic modeling applications.

California↗

Accumulation of atmospheric sulfur in some Costa Rican soils

Sulfur is one of the macronutrient elements whose sources to terrestrial ecosystems should shift from dominance by rock-weathering to atmospheric deposition as soils and underlying substrate undergo progressive weathering and leaching. However, the nature and timing of this transition is not well known. We investigated sources of sulfur to tropical rain forests growing on basalt-derived soils in the Osa Peninsula region of Costa Rica. Sulfur sources were examined using stable isotope ratios (δ 34 S) and compared to chemical indices of soil development. The most weathered soils, and the forests they supported, are dominated by atmospheric sulfur, while a less weathered soil type contains both rock-derived and atmospheric sulfur. Patterns of increasing δ 34 S with increasing soil sulfur concentration across the landscape suggest atmospheric sulfur is accumulating, and little rock-derived sulfur has been retained. Soil sulfur, minus adsorbed sulfate, is correlated with carbon and nitrogen, implying that sulfur accumulation occurs as plants and microbes incorporate sulfur into organic matter. Only the lower depth increments of the more weathered soils contained significant adsorbed sulfate. The evidence suggests a pattern of soil development in which sulfur-bearing minerals in rock, such as sulfides, weather early relative to other minerals, and the released sulfate is leached away. Sulfur added via atmospheric deposition is retained as organic matter accumulates in the soil profile. Adsorbed sulfate accumulates later, driven by changes in soil chemistry and mineralogy. These aspects of sulfur behavior during pedogenesis in this environment may hasten the transition to dominance by atmospheric sources.

Journal of Geophysical Research: Biogeosciences↗

Boreal soil carbon dynamics under a changing climate: A model inversion approach

Several fundamental but important factors controlling the feedback of boreal organic carbon (OC) to climate change were examined using a mechanistic model of soil OC dynamics, including the combined effects of temperature and moisture on the decomposition of OC and the factors controlling carbon quality and decomposition with depth. To estimate decomposition rates and evaluate their variations with depth, the model was inverted using a global optimization algorithm. Three sites with different drainage conditions that represent a broad diversity of boreal black spruce ecosystems were modeled. The comparison among the models with different depth patterns of decomposition rates (i.e., constant, linear, and exponential decrease) revealed that the model with constant inherent decomposition rates through the soil profile was able to fit the observed data in the most efficient way. There were also lower turnover times in the wettest site compared to the drier site even after accounting for moisture and temperature differences. Taken together, these results indicate that decomposition (especially for the wetter site) was not accurately represented with standard moisture and temperature controls and that other important protection mechanisms (e.g., limitation of O 2 , redox conditions, and permafrost) rather than low inherent decomposition rates are responsible for the recalcitrance of deep OC. The simulation results also showed that most of the soil CO 2 efflux is generated from subsurface layers of OC because of the large OC stocks and optimal moisture conditions, suggesting that these deeper soil OC stocks are likely to be critically important to the future carbon dynamics.

Journal of Geophysical Research Biogeosciences↗

Linking soil organic matter dynamics and erosion-induced terrestrial carbon sequestration at different landform positions

Recently, the potential for terrestrial carbon (C) sequestration by soil erosion and deposition has received increased interest. Erosion and deposition constitute a sink for atmospheric carbon dioxide relative to a preerosional state or a noneroding scenario, if the posterosion watershed C balance is increased due to (1) partial replacement of eroded C by new photosynthate in the eroded site; and (2) preservation from decomposition of at least some eroded soil organic carbon (SOC) arriving in depositional settings. Little is known, however, about differences in C dynamics at different erosional and depositional landform positions within the same eroding system. We determined the contribution of different landform positions to erosion-induced terrestrial C sequestration by measuring rates of net primary productivity (NPP), replacement of eroded C, and decomposition of organic matter (OM) at four categorically different landform positions within a naturally eroding toposequence in northern California. We found that eroded C is replaced by NPP 15 times over in the summit of the site studied and 5 times over in the slope. Profile-averaged, long-term rate constant for SOM decomposition was 2 to 14 times slower in the depositional settings compared with that in eroding slopes. As a result, the inventory of C in the depositional settings was 2 to 3 times larger than that of the eroding positions. Owing to both C replacement at eroding sites and reduced rates of OM decomposition in depositional sites, soil erosion constitutes a C sink from the atmosphere at our study site.

California↗

Hydrologic support of carbon dioxide flux revealed by whole-lake carbon budgets

Freshwater lakes are an important component of the global carbon cycle through both organic carbon (OC) sequestration and carbon dioxide (CO 2 ) emission. Most lakes have a net annual loss of CO 2 to the atmosphere and substantial current evidence suggests that biologic mineralization of allochthonous OC maintains this flux. Because net CO 2 flux to the atmosphere implies net mineralization of OC within the lake ecosystem, it is also commonly assumed that net annual CO 2 emission indicates negative net ecosystem production (NEP). We explored the relationship between atmospheric CO 2 emission and NEP in two lakes known to have contrasting hydrologie characteristics and net CO 2 emission. We calculated NEP for calendar year 2004 using whole-lake OC and inorganic carbon (IC) budgets, NEP oc and NEP IC , respectively, and compared the resulting values to measured annual CO 2 flux from the lakes. In both lakes, NEP Ic and NEP Ic were positive, indicating net autotrophy. Therefore CO 2 emission from these lakes was apparently not supported by mineralization of allochthonous organic material. In both lakes, hydrologie CO 2 inputs, as well as CO 2 evolved from netcalcite precipitation, could account for the net CO 2 emission. NEP calculated from diel CO 2 measurements was also affected by hydrologie inputs of CO 2 . These results indicate that CO 2 emission and positive NEP may coincide in lakes, especially in carbonate terrain, and that all potential geologic, biogeochemical, and hydrologie sources of CO 2 need to be accounted for when using CO 2 concentrations to infer lake NEP. Copyright 2009 by the American Geophysical Union.

Journal of Geophysical Research G: Biogeosciences↗