Historic emissions from deforestation and forest degradation in Mato Grosso,Brazil: 1) source data uncertainties

Background  

Historic carbon emissions are an important foundation for proposed efforts to Reduce Emissions from Deforestation and forest Degradation and enhance forest carbon stocks through conservation and sustainable forest management (REDD+). The level of uncertainty in historic carbon emissions estimates is also critical for REDD+, since high uncertainties could limit climate benefits from credited mitigation actions. Here, we analyzed source data uncertainties based on the range of available deforestation, forest degradation, and forest carbon stock estimates for the Brazilian state of Mato Grosso during 1990-2008.     

A new multiscale scheme for computing statistical moments in single phase flow in heterogeneous porous media

In this work we develop a new multiscale procedure to compute numerically the statistical moments of the stochastic variables which govern single phase flow in heterogeneous porous media. The technique explores the properties of the log-normally distributed hydraulic conductivity, characterized by power-law or exponential covariances, which shows invariance in its statistical structure upon a simultaneous change of the scale of observation and strength of heterogeneity. We construct a family of equivalent stochastic hydrodynamic variables satisfying the same flow equations at different scales and strengths of heterogeneity or correlation lengths. Within the new procedure the governing equations are solved in a scaled geology and the numerical results are mapped onto the original medium at coarser scales by a straightforward rescaling. The new procedure is implemented numerically within the Monte Carlo algorithm and also in conjunction with the discretization of the low-order effective equations derived from perturbation analysis. Numerical results obtained by the finite element method show the accuracy of the new procedure to approximated the two first moments of the pressure and velocity along with its potential in reducing drastically the computational cost involved in the numerical modeling of both power-law and exponential covariance functions.     

An enhanced ensemble Kalman filter scheme incorporating model error in sequential coupling between flow and geomechanics

In this work, we construct a new methodology for enhancing the predictive accuracy of sequential methods for coupling flow and geomechanics while preserving low computational cost. The new computational approach is developed within the framework of the fixed-stress split algorithm procedure in conjunction with data assimilation based on the ensemble Kalman filter (EnKF). In this context, we identify the high-fidelity model with the two-way formulation where additional source term appears in the flow equation containing the time derivative of total mean stress. The iterative scheme is then interlaced with data assimilation steps, which also incorporate the modeling error inherent to the EnKF framework. Such a procedure gives rise to an “enhanced one-way formulation,” exhibiting substantial improvement in accuracy compared with the classical one-way method. The governing equations are discretized by mixed finite elements, and numerical simulation of a 2D slab problem between injection and production wells illustrate the tremendous achievement of the method proposed herein.     

A new methodology for computing ionic profiles and disjoining pressure in swelling porous media

A new two-scale computational model is proposed to construct the constitutive law of the swelling pressure which appears in the modified form of the macroscopic effective stress principle for expansive clays saturated by an aqueous electrolyte solution containing multivalent ionic species. The microscopic non-local nanoscale model is constructed based on a coupled Poisson-Fredholm integral equation arising from the thermodynamics of inhomogeneous fluids in nanopores (Density Functional Theory), which governs the local electric double layer potential profile coupled with the ion-particle correlation function in an electrolytic solution of finite size ions. The local problem is discretized by invoking the eigenvalue expansion of the convolution kernel in conjunction with the Galerkin method for the Gauss-Poisson equation. The discretization of the Fredholm equation is accomplished by a collocation scheme employing eigenfunction basis. Numerical simulations of the local ionic profiles in rectangular cell geometries are obtained showing considerable discrepancies with those computed with Poisson-Boltzmann based models for point charges, particularly for divalent ions in calcium montmorillonite. The constitutive law for the disjoining pressure is reconstructed numerically by invoking the contact theorem within a post-processing approach. The resultant computational model is capable of capturing ranges of particle attraction characterized by negative values of the disjoining pressure overlooked by the classical electric double layer theory. Such results provide further insight in the role the swelling pressure plays in the modified macroscopic effective stress principle for expansive porous media.     

C-, Sr-isotope and Hg chemostratigraphy of Neoproterozoic cap carbonates of the Sergipano Belt,Northeastern Brazil

Two cap carbonates overlying glaciogenic diamictites crop out extensively in the eastern Vaza Barris Domain of the Sergipano Belt, northeastern Brazil. They are represented by carbonates of the Jacoca Formation, resting on top of diamictite of the Ribeiropolis Formation, and by the Olhos D’Agua Formation (carbonates, organic-rich towards the top), which overlies diamictite of the Palestina Formation. These two sequences were deformed and metamorphosed at sub-greenschist facies-conditions during the Brasiliano cycle (650–600 Ma).     

A new computational strategy for solving two‐phase flow in strongly heterogeneous poroelastic media of evolving scales

We develop a new computational methodology for solving two‐phase flow in highly heterogeneous porous media incorporating geomechanical coupling subject to uncertainty in the poromechanical parameters. Within the framework of a staggered‐in‐time coupling algorithm, the numerical method proposed herein relies on a Petrov–Galerkin postprocessing approach projected on the Raviart–Thomas space to compute the Darcy velocity of the mixture in conjunction with a locally conservative higher order finite volume discretization of the nonlinear transport equation for the saturation and an operator splitting procedure based on the difference in the time‐scales of transport and geomechanics to compute the effects of transient porosity upon saturation. Notable features of the numerical modeling proposed herein are the local conservation properties inherited by the discrete fluxes that are crucial to correctly capture the fingering patterns arising from the interaction between heterogeneity and nonlinear viscous coupling. Water flooding in a poroelastic formation subject to an overburden is simulated with the geology characterized by multiscale self‐similar permeability and Young modulus random fields with power‐law covariance structure. Statistical moments of the poromechanical unknowns are computed within the framework of a high‐resolution Monte Carlo method. Numerical results illustrate the necessity of adopting locally conservative schemes to obtain reliable predictions of secondary recovery and finger growth in strongly heterogeneous deformable reservoirs. Copyright © 2011 John Wiley & Sons, Ltd.     

Nd isotopic composition of Paleoproterozoic volcanic and granitoid rocks of Vila Riozinho: implications for the crustal evolution of the Tapajós gold province, Amazon craton

Nd data from the Paleoproterozoic magmatic rocks of Vila Riozinho and Jamanxim (Tapajós gold province) indicate that original magmas were not produced exclusively by the remelting of Archean sialic crust and point to dominant Paleoproterozoic sources. ε_Nd(T) values preclude derivation from mantle sources for the ca. 2.0 Ga Vila Riozinho volcanics and older São Jorge granite. They may represent a subduction-related magmatic arc with magmas modified by interaction with crust or a post- to late-orogenic remelting of an older Paleoproterozoic juvenile arc with minimal contribution from the Archean crust. The origin of the 1.88 Ga Parauari, Maloquinha, Iriri, and Moraes Almeida igneous associations and the Jamanxim rhyolites has been attributed to large-scale taphrogenesis that marked the breakup of a large Paleoproterozoic continent. Derivation of the original magmas from the remelting of crustal sources older than ca. 1.9 Ga is consistent with geochemical and Nd isotopic data. Archean remnants probably occur between the Paleoproterozoic terrains of the Ventuari-Tapajós province. Archean terrains of the Amazon craton extend from the Xingu to the Itaituba region but have not been identified in the southern Guyana shield. Thus, data reveal that the boundaries between the central Amazon and Ventuari-Tapajós provinces need better definition and more detailed field and geochronological work.     

Geomorphological controls on coastal vegetation at the Virginia Coast Reserve

The establishment and succession of vegetation on migrating, low-profile barrier islands is greatly affected by the physical hydrogeomorphological processes that regulate island topography, saline and fresh groundwater table surfaces. Apart from the physical destruction of plants by overwash processes, fluctuations in water table elevations and variations in groundwater salinity, both spatially and temporally, also appear to have significant impact on the nature and distribution of vegetation on these islands. Species composition, community structure and biodiversity on the Virginia barrier islands are controlled by the same processes that give rise to landforms and maintain their form. These processes include marine water inundations, groundwater salinity variations and changes in depth to the fresh-water table. Land surface elevation, landform morphology and position on the barrier island determine exposure to high tides, storm surges, sand burial, and the extent of the fresh-water reserves.In this article, the underpinnings of a Long-term Ecological Research Program in which 25 geologists, geomorphologists, climatologists, and ecologists have a common research plan is presented and several examples of the product of this research partnership dealing with geomorphological and hydrologic controls on vegetation dynamics are detailed. Among the aspects of ecological dynamics examined in terms of geomorphological processes are vegetation zonation, succession, disturbance, and ecosystem state change.     

Benthic species as mud patrol - modelled effects of bioturbators and biofilms on large-scale estuarine mud and morphology

Sediment-stabilizing and -destabilizing organisms, i.e. microphytobenthos (biofilms) and macrozoobenthos (bioturbators), affect the erodibility of muddy sediments, potentially altering large-scale estuarine morphology. Using a novel eco-morphodynamic model of an idealized estuary, we investigate eco-engineering effects of microphytobenthos and two macrozoobenthic bioturbators. Local mud erodibility is based on species pattern predicted through hydrodynamics, soil mud content, competition and grazing. Mud resuspension and export is enhanced under bioturbation and prevented under biostabilization through respective exposure and protection of the supra- and intertidal. Bioturbation decreases mud thickness and bed elevations, which increases net mud fluxes. Microphytobenthos reduces erosion, leading to a local mud increase of intertidal sediments. In multi-species scenarios, an effective mud-prone bioturbator strongly alters morphology, exceeding that of a more abundant sand-prone moderate species, showing that morphological change depends on species traits as opposed to abundance. Altering their habitat, the effective mud-prone bioturbator facilitates expansion of the sand-prone moderate bioturbator. Grazing and species competition favor species distributions of dominant bioturbators. Consequently, eco-engineering affects habitat conditions while species interactions determine species dominance. Our results show that eco-engineering species determine the mud content of the estuary, which suggests large effects on the morphology of estuaries with aggravating habitat degradation.