Improving the assessment of ICESat water altimetry accuracy accounting for autocorrelation

Given that water resources are scarce and are strained by competing demands, it has become crucial to develop and improve techniques to observe the temporal and spatial variations in the inland water volume. Due to the lack of data and the heterogeneity of water level stations, remote sensing, and especially altimetry from space, appear as complementary techniques for water level monitoring. In addition to spatial resolution and sampling rates in space or time, one of the most relevant criteria for satellite altimetry on inland water is the accuracy of the elevation data. Here, the accuracy of ICESat LIDAR altimetry product is assessed over the Great Lakes in North America. The accuracy assessment method used in this paper emphasizes on autocorrelation in high temporal frequency ICESat measurements. It also considers uncertainties resulting from both in situ lake level reference data. A probabilistic upscaling process was developed. This process is based on several successive ICESat shots averaged in a spatial transect accounting for autocorrelation between successive shots. The method also applies pre-processing of the ICESat data with saturation correction of ICESat waveforms, spatial filtering to avoid measurement disturbance from the land–water transition effects on waveform saturation and data selection to avoid trends in water elevations across space. Initially this paper analyzes 237 collected ICESat transects, consistent with the available hydrometric ground stations for four of the Great Lakes. By adapting a geostatistical framework, a high frequency autocorrelation between successive shot elevation values was observed and then modeled for 45% of the 237 transects. The modeled autocorrelation was therefore used to estimate water elevations at the transect scale and the resulting uncertainty for the 117 transects without trend. This uncertainty was 8 times greater than the usual computed uncertainty, when no temporal correlation is taken into account. This temporal correlation, corresponding to approximately 11 consecutive ICESat shots, could be linked to low transmitted ICESat GLAS energy and to poor weather conditions. Assuming Gaussian uncertainties for both reference data and ICESat data upscaled at the transect scale, we derived GLAS deviations statistics by averaging the results at station and lake scales. An overall bias of −4.6 cm (underestimation) and an overall standard deviation of 11.6 cm were computed for all lakes. Results demonstrated the relevance of taking autocorrelation into account in satellite data uncertainty assesment.     

Estimating River Conductance from Prior Information to Improve Surface‐Subsurface Model Calibration

Most groundwater models simulate stream‐aquifer interactions with a head‐dependent flux boundary condition based on a river conductance (CRIV). CRIV is usually calibrated with other parameters by history matching. However, the inverse problem of groundwater models is often ill‐posed and individual model parameters are likely to be poorly constrained. Ill‐posedness can be addressed by Tikhonov regularization with prior knowledge on parameter values. The difficulty with a lumped parameter like CRIV, which cannot be measured in the field, is to find suitable initial and regularization values. Several formulations have been proposed for the estimation of CRIV from physical parameters. However, these methods are either too simple to provide a reliable estimate of CRIV, or too complex to be easily implemented by groundwater modelers. This paper addresses the issue with a flexible and operational tool based on a 2D numerical model in a local vertical cross section, where the river conductance is computed from selected geometric and hydrodynamic parameters. Contrary to other approaches, the grid size of the regional model and the anisotropy of the aquifer hydraulic conductivity are also taken into account. A global sensitivity analysis indicates the strong sensitivity of CRIV to these parameters. This enhancement for the prior estimation of CRIV is a step forward for the calibration and uncertainty analysis of surface‐subsurface models. It is especially useful for modeling objectives that require CRIV to be well known such as conjunctive surface water‐groundwater use.     

Nickel speciation in Sebertia acuminata,a plant growing on a lateritic soil of New Caledonia

Nickel speciation in a nickel hyperaccumulating plant (Sebertia acuminata) and its associated soil of southern New Caledonia was studied using various analytical methods. The soil is formed of iron oxides (goethite, hematite), which contain almost all the nickel. The available nickel is probably linked to the organic matter in the litter. Sebertia acuminata, acts as a nickel pump, and concentrates the metal in its leaves. It partitions nickel and silica; nickel is concentrated in the cells (probably in the vacuoles) as organometallic complexes, whereas silica forms the framework of the cells, and the phytolithes. A thorough study of these plants seems essential in order to define the soil–plant relations, and to propose appropriate ways for ecological restoration. To cite this article: N. Perrier et al., C. R. Geoscience 336 (2004).     

Spectral and geological study of the sulfate-rich region of West Candor Chasma, Mars

Sulfates have been discovered by the OMEGA spectrometer in different locations of the planet Mars. They are strongly correlated to light toned layered deposits in the equatorial regions. West Candor Chasma is the canyon with the thickest stack of layers and one with the largest area covered by sulfates. A detailed study coupling mineralogy derived from OMEGA spectral data and geology derived from HRSC imager and other datasets leads to some straightforward issues. The monohydrated sulfate kieserite is found mainly over heavily eroded scarps of light toned material. It likely corresponds to a mineral present in the initial rock formed either during formation and diagenesis of sediments, or during hydrothermal alteration at depth, because it is typically found on outcrops that are eroded and steep. Polyhydrated sulfates, that match any Ca-, Na-, Fe-, or Mg-sulfates with more than one water molecule, are preferentially present on less eroded and darker outcrops than outcrops of kieserite. These variations can be the result of a diversity in the composition and/or of the rehydration of kieserite on surfaces with longer exposure. The latter possibility of rehydration in the current, or recent, atmosphere suggests the low surface temperatures preserve sulfates from desiccation, and, also can rehydrate part of them. Strong signatures of iron oxides are present on sulfate-rich scarps and at the base of layered deposits scarps. They are correlated with TES gray hematite signature and might correspond to iron oxides present in the rock as sand-size grains, or possibly larger concretions, that are eroded and transported down by gravity at the base of the scarp. Pyroxenes are present mainly on sand dunes in the low lying terrains. Pyroxene is strongly depleted or absent in the layered deposits. When mixed with kieserite, local observations favor a spatial mixing with dunes over layered deposits. Sulfates such as those detected in the studied area require the presence of liquid water to form by precipitation, either in an intermittent lacustrine environment or by hydrothermal fluid circulation. Both possibilities require the presence of sulfur-rich groundwater to explain fluid circulation. The elevation of the uppermost sulfate signatures suggests the presence of aquifers up to 2.5 km above datum, only 1 km below the plateau surface.     

Mass‐dependent fractionation of nickel isotopes in meteoritic metal

Abstract— We measured nickel isotopes via multicollector inductively coupled plasma mass spectrometry (MC‐ICPMS) in the bulk metal from 36 meteorites, including chondrites, pallasites, and irons (magmatic and non‐magmatic). The Ni isotopes in these meteorites are mass fractionated; the fractionation spans an overall range of ～0.4‰ amu^?1. The ranges of Ni isotopic compositions (relative to the SRM 986 Ni isotopic standard) in metal from iron meteorites (～0.0 to ～0.3‰ amu^?1) and chondrites (～0.0 to ～0.2‰ amu^?1) are similar, whereas the range in pallasite metal (～–0.1 to 0.0‰ amu^?1) appears distinct. The fractionation of Ni isotopes within a suite of fourteen IIIAB irons (～0.0 to ～0.3‰ amu^?1) spans the entire range measured in all magmatic irons. However, the degree of Ni isotopic fractionation in these samples does not correlate with their Ni content, suggesting that core crystallization did not fractionate Ni isotopes in a systematic way. We also measured the Ni and Fe isotopes in adjacent kamacite and taenite from the Toluca IAB iron meteorite. Nickel isotopes show clearly resolvable fractionation between these two phases; kamacite is heavier relative to taenite by ～0.4‰ amu^?1. In contrast, the Fe isotopes do not show a resolvable fractionation between kamacite and taenite. The observed isotopic compositions of kamacite and taenite can be understood in terms of kinetic fractionation due to diffusion of Ni during cooling of the Fe‐Ni alloy and the development of the Widmanstätten pattern.     

Sedimentary evidence of deglacial megafloods in the northern Gulf of Mexico (Pigmy Basin)

Cored sediments from the Pigmy Basin, northern Gulf of Mexico, were analyzed in order to better constrain late deglacial and early Holocene paleoenvironmental and sedimentary changes in response to North American climate evolution. Mineralogical and geochemical proxies indicate the succession of two sedimentary regimes: dominantly detrital during the deglaciation (15–12.9 cal ka BP) whereas biogenic contribution relatively increased later on during the Younger Dryas and early Holocene (12.9 and 10 cal ka BP). Geochemical data reveal that the deglacial record mainly reflects variations of terrigenous supply via the Mississippi River rather than modifications of redox conditions in the basin. Specific variations of almost all the parameters measured in this paper are synchronous with the main deglacial meltwater episode (Meltwater Spike) described or modeled in previous marine or continental studies. During this episode, most parameters display “stair-step-like” – pattern variations highlighting three successive steps within the main meltwater flow. Variations in grain-size and clay mineral assemblage recorded in the Pigmy Basin indicate that the erosional regime was very strong on land during the first part of the Meltwater Spike, and then milder, inducing more subtle modifications in the sedimentary regime in this part of the Gulf. Specific geochemical and mineralogical signatures (notably, clay minerals and trace metal geochemistry) pinpoint a dominant origin from NW North America for detrital particles reflecting meltwater outflow from the south-western Laurentide Ice Sheet (LIS) margin during the most intense freshwater discharge. The observed decrease of the sedimentation rate from about 200 to 25 cm/ka at ca 12.9 ka evidenced a drastic decrease of erosional processes during late phase of discharge, consistently with the hypotheses of major reduction of meltwater flow. The major modification at 12.9 cal ka BP is interpreted to result from both modifications of the main Mississippi fluvial regime due to eastward and northward rerouting of meltwater flow at the onset of the Younger Dryas, and the increase of sea-surface temperature linked to insolation. Finally, slight grain-size modifications suggest that some freshwater discharges may have episodically reached the Gulf of Mexico after the Younger Dryas reflecting possible small adjustments of the postglacial hydrological regime.     

Ca isotope cycling in a forested ecosystem

Reports of large Ca isotope fractionations between trees and soils prompted this study of a Boreal forest ecosystem near La Ronge, Saskatchewan, to improve understanding of this phenomenon. The results on five tree species (black spruce, trembling aspen, white spruce, jack pine, balsam poplar) confirm that nutrient Ca uptake by plants favors the light isotopes, thus driving residual Ca in plant available soil pools towards enrichment in the heavy isotopes. Substantial within-tree fraction occurs in tissues formed along the transpiration stream, with low δ⁴⁴Ca values in fine roots (2 mm), intermediate values in stemwood, and high values in foliage. Separation factors between different plant tissues are similar between species, but the initial fractionation step in the tips of the fine roots is species specific, and/or sensitive to the local soil environment. Soil water δ⁴⁴Ca values appear to increase with depth to at least 35 cm below the top of the forest floor, which is close to the deepest level of fine roots. The heavy plant fractionated signature of Ca in the finely rooted upper soils filters downward where it is retained on ion exchange sites, leached into groundwater, and discharged into surface waters.The relationship between Ca uptake by tree fine roots and the pattern of δ⁴⁴Ca enrichment with soil depth was modeled for two Ca pools: the forest floor (litter) and the underlying (upper B) mineral soil. Six study plots were investigated along two hillside toposequences trending upwards from a first order stream. We used allometric equations describing the Ca distribution in boreal tree species to calculate weighted average δ⁴⁴Ca values for the stands in each plot and estimate Ca uptake rates. The δ⁴⁴Ca value of precipitation was measured, and soil weathering signatures deduced, by acid leaching of lower B mineral soils. Steady state equations were used to derive a set of model Ca fluxes and fractionation factors for each plot. The model reproduces the increase in δ⁴⁴Ca with depth found in forest floor and upper B soil waters. Transient model runs show that the forest Ca cycle is sensitive to changes in plant Ca uptake rate, such as would occur during ontogeny or disturbance. Accordingly, secular records of δ⁴⁴Ca in tree ring cellulose have the potential to monitor changes in the forest Ca cycle through time, thus providing a new tool for evaluating natural and anthropogenic impacts on forest health. Another model run shows that by changing the size of the isotope fractionation factor and adjusting for differences in forest productivity, that the range in Ca isotope fractionation in forested ecosystems reported in the literature, thus far, is reproduced. As a quantitative tool, the Ca cycling model produces a reasonable set of relative Ca fluxes for the La Ronge site, consistent with Environment Canada’s measurements for wet deposition in the region and simulated Ca release from soil mineral weathering using the PROFILE model. But the sensitivity of the model is limited by the small range of fractionation observed in this boreal shield setting of ∼1‰, which limits accuracy. If the model were applied to a site with a greater range in δ⁴⁴Ca values among the principal Ca fluxes, it is capable of producing robust and reliable estimations of Ca fluxes that are otherwise difficult to measure in forested ecosystems.