Global land mapping of satellite-observed CO2 total columns using spatio-temporal geostatistics

ABSTRACT

This study presents an approach for generating a global land mapping dataset of the satellite measurements of CO₂ total column (XCO₂) using spatio-temporal geostatistics, which makes full use of the joint spatial and temporal dependencies between observations. The mapping approach considers the latitude-zonal seasonal cycles and spatio-temporal correlation structure of XCO₂, and obtains global land maps of XCO₂, with a spatial grid resolution of 1° latitude by 1° longitude and temporal resolution of 3 days. We evaluate the accuracy and uncertainty of the mapping dataset in the following three ways: (1) in cross-validation, the mapping approach results in a high correlation coefficient of 0.94 between the predictions and observations, (2) in comparison with ground truth provided by the Total Carbon Column Observing Network (TCCON), the predicted XCO₂ time series and those from TCCON sites are in good agreement, with an overall bias of 0.01?ppm and a standard deviation of the difference of 1.22?ppm and (3) in comparison with model simulations, the spatio-temporal variability of XCO₂ between the mapping dataset and simulations from the CT2013 and GEOS-Chem are generally consistent. The generated mapping XCO₂ data in this study provides a new global geospatial dataset in global understanding of greenhouse gases dynamics and global warming.     

Isotopically diverse rhyolites coeval with the Columbia River Flood Basalts: evidence for mantle plume interaction with the continental crust

The Columbia River Flood Basalts (CRB) of the northwestern USA are coeval with eruptions of several thousand km³ of rhyolite. A broad survey of major phenocryst oxygen isotopes and of O and Hf isotopes in zircons from these rhyolites reveals significant diversity in inferred δ¹⁸O_melt values, ranging from +1.9 to +10.5‰ (SMOW), and in zircon Hf isotope compositions, which range from ε_Hf = ?39 to +9. This newly identified isotopic diversity shows that the syn‐CRB rhyolites were derived from high‐percentage melting of the crust. Low‐δ¹⁸O rhyolites, which fingerprint the melting of hydrothermally altered crust, are concentrated at the edge of the North American craton. This suggests that the conditions of crustal heating, faulting, and hydrothermal alteration required for the production of these rhyolites were concentrated there by the contrasts in crustal thickness and rheology associated with the boundary between the North American craton and younger accreted terranes.     

The use of ablation‐dominated medial moraines as samplers for 10Be‐derived erosion rates of glacier valley walls,Kichatna Mountains,AK

We use cosmogenic ¹⁰Be concentrations in amalgamated rock samples from active, ice‐cored medial moraines to constrain glacial valley sidewall backwearing rates in the Kichatna Mountains, Alaska Range, Alaska. This dramatic landscape is carved into a small ～65 Ma granitic pluton about 100 km west of Denali, where kilometer‐tall rock walls and ‘cathedral’ spires tower over a radial array of over a dozen valley glaciers. These supraglacial landforms erode primarily by rockfall, but erosion rates are difficult to determine. We use cosmogenic ¹⁰Be to measure rockwall backwearing rates on timescales of 10³–10⁴ years, with a straightforward sampling strategy that exploits ablation‐dominated medial moraines. A medial moraine and its associated englacial debris serve as a conveyor system, bringing supraglacial rockfall debris from accumulation‐zone valley walls to the moraine crest in the ablation zone. We discuss quantitatively several factors that complicate interpretation of cosmogenic concentrations in this material, including the complex scaling of production rates in very steep terrain, the stochastic nature of the rockfall erosion process, the unmixed nature of the moraine sediment, and additional cosmogenic accumulation during transport of the sediment. We sampled medial moraines on each of three glaciers of different sizes and topographic aspects. All three moraines are sourced in areas with identical rock and similar sidewall relief of ～1 km. Each sample was amalgamated from 25 to 35 clasts collected over a 1‐km longitudinal transect of each moraine. Two of the glaciers yield similar ¹⁰Be concentrations (～1·6–2·2 × 10⁴ at/g) and minimum sidewall slope‐normal erosion rates (～0·5–0·7 mm/yr). The lowest ¹⁰Be concentrations (8 × 10³ at/g) and the highest erosion rates (1·3 mm/yr) come from the largest glacier in the range with the lowest late‐summer snowline. These rates are reasonable in an alpine glacial setting, and are much faster than long‐term exhumation rates of the western Alaska Range as determined by thermochronometric studies. Copyright © 2010 John Wiley & Sons, Ltd.     

Preliminary estimates of regolith generation and mobility in the Susquehanna Shale Hills Critical Zone Observatory,Pennsylvania, using meteoric 10Be

This study seeks to quantify the rate and timing of regolith generation in the Critical Zone at the Susquehanna Shale Hills Critical Zone Observatory (SSHO). Meteoric ¹⁰Be depth profiles were determined using measurements from 30 hillslope soil and bedrock core samples in an effort to constrain ¹⁰Be inventories. The SSHO is located in the temperate climate zone of central Pennsylvania and comprises a first-order watershed developed entirely on a Fe-rich, organic-poor, Silurian-aged shale. Two major perturbations to the landscape have occurred at SSHO in the geologically recent past, including significant and sustained periglacial activity until after the retreat of the Laurentide ice sheet (～21 ka) and deforestation during early colonial land-use. Bulk soil samples (n = 16) were collected at three locations along a planar hillslope on the southern ridge of the catchment, representing the ridge top, mid-slope and valley floor. Rock chip samples (n = 14) were also collected from a 24 m deep core drilled into the northern ridge top. All meteoric ¹⁰Be concentration profiles show a declining trend with depth, with most of the ¹⁰Be retained in the uppermost decimeters of the soil. Meteoric ¹⁰Be inventories are higher at the mid-slope and valley floor sample sites, at 3.71 ± 0.02 × 10¹⁰ at/cm² and 3.69 ± 0.02 × 10¹⁰ at/cm², than at the ridge top site (1.90 ± 0.01 × 10¹⁰ at/cm²). The ¹⁰Be inventory at the convex ridge top site implies a minimum residence time of ～10.6 ka, or if erosion is steady, an erosion rate of 19.4 ± 0.2 m/My.     

Contaminant remediation decision analysis using information gap theory

Decision making under severe lack of information is a ubiquitous situation in nearly every applied field of engineering, policy, and science. A severe lack of information precludes our ability to determine a frequency of occurrence of events or conditions that impact the decision; therefore, decision uncertainties due to a severe lack of information cannot be characterized probabilistically. To circumvent this problem, information gap (info-gap) theory has been developed to explicitly recognize and quantify the implications of a severe lack of information in decision making. This paper presents a decision analysis based on info-gap theory developed for a contaminant remediation scenario. The analysis provides decision support in determining the fraction of contaminant mass to remove from the environment. An info-gap uncertainty model is developed to characterize uncertainty due to a lack of information concerning the contaminant flux. The info-gap uncertainty model groups nested, convex sets of functions defining contaminant flux over time based on their level of deviation from a nominal contaminant flux. The nominal contaminant flux defines a best estimate of contaminant flux over time based on existing, though incomplete, information. A robustness function is derived to quantify the maximum level of deviation from nominal that still ensures compliance for alternative decisions. An opportuneness function is derived to characterize the possibility of meeting a desired contaminant concentration level. The decision analysis evaluates how the robustness and opportuneness change as a function of time since remediation and as a function of the fraction of contaminant mass removed.     

New constraints on the late Cenozoic incision history of the New River, Virginia

The New River crosses three physiogeologic provinces of the ancient, tectonically quiescent Appalachian orogen and is ideally situated to record variability in fluvial erosion rates over the late Cenozoic. Active erosion features on resistant bedrock that floors the river at prominent knickpoints demonstrate that the river is currently incising toward base level. However, thick sequences of alluvial fill and fluvial terraces cut into this fill record an incision history for the river that includes several periods of stalled downcutting and aggradation. We used cosmogenic ¹⁰Be exposure dating, aided by mapping and sedimentological examination of terrace deposits, to constrain the timing of events in this history. ¹⁰Be concentration depth profiles were used to help account for variables such as cosmogenic inheritance and terrace bioturbation. Fill-cut and strath terraces at elevations 10, 20, and 50 m above the modern river yield model cosmogenic exposure ages of 130, 600, and 600–950 ka, respectively, but uncertainties on these ages are not well constrained. These results provide the first direct constraint on the history of alluvial aggradation and incision events recorded by New River terrace deposits. The exposure ages yield a long-term average incision rate of 43 m/my, which is comparable to rates measured elsewhere in the Appalachians. During specific intervals over the last 1 Ma, however, the New River's incision rate reached 100 m/my. Modern erosion rates on bedrock at a prominent knickpoint are between 28 and 87 m/my, in good agreement with rates calculated between terrace abandonment events and significantly faster than 2 m/my rates of surface erosion from ancient terrace remnants. Fluctuations between aggradation and rapid incision operate on timescales of 10⁴− 10⁵ year, similar to those of late Cenozoic climate variations, though uncertainties in model ages preclude direct correlation of these fluctuations to specific climate change events. These second-order fluctuations appear within a longer-term signal of dominant aggradation (until 2 Ma) followed by dominant incision. A similar signal is observed on other Appalachian rivers and may be the result of sediment supply fluctuations driven by the increased frequency of climate changes in the late Cenozoic.     

Deaggregation of multi-hazard damages,losses, risks,and connectivity: an application to the joint seismic-tsunami hazard at Seaside,Oregon

This paper presents a methodology to deaggregate the results of a multi-hazard damage analysis by extending the traditional multi-hazard damage analysis to consider both population characteristics and independent hazards. The methodology is applied to the joint seismic-tsunami hazard at Seaside, Oregon, considering four infrastructure systems: (1) buildings, (2) transportation network, (3) electric power network and (4) water supply network. Damages to all infrastructure systems are evaluated, and the networked infrastructures are used to inform parcel connectivity to critical facilities. US Census data and a probabilistic housing unit allocation method are implemented to assign detailed household demographic characteristics at the parcel level. Six dimensions of deaggregation are introduced: (1) spatial, (2) hazard type, (3) hazard intensity, (4) infrastructure system, (5) infrastructure component, and (6) housing unit characteristics. The damages, economic losses and risks, and connectivity to critical facilities are deaggregated across these six dimensions. The results show that deaggregated economic loss and risk plots can allow community resilience planners the ability to isolate high-risk events, as well as provide insights into the underlying driving forces. Geospatial representation of the results allows for the identification of both vulnerable buildings and areas within a community and is highlighted by the spatial pattern of parcel disconnection from critical facilities. The incorporation of population characteristics provides an understanding of how hazards disproportionately impact population subgroups and can aide in equitable resilience planning.

     

Surface Estimation for Multiple Misaligned Point Sets

Mathematical Geosciences - Two common tasks when processing point cloud data sets are surface estimation and point cloud registration. In this paper, a statistical approach is developed to solve...