Species richness,diversity and density of understory vegetation along disturbance gradients in the Himalayan conifer forest

We investigated whether species richness, diversity and density of understory herbaceous plants differed along logging(gap) and grazing(primarily by cattle) disturbance gradients, and sought to identify drivers of richness, diversity and density of understory vegetation of logged sites. A factorial experiment was conducted in the mixed conifer forest of Gidakom in Western Bhutan. Levels of the logging treatment included small(0.15 – 0.24 ha), medium(0.25 – 0.35 ha) and large(0.36 – 1.31 ha) gaps. The grazing treatment included grazed(primarily by cattle) and ungrazed(where herbivores were excluded by a fence) plots nested within each gap. Data were collected from 12 gaps(4 replicates at each level of logging) using the point intercept method. Shannon Weaver Diversity and Margalef's indices were used to estimate species diversity and describe species richness, respectively. Soil samples were analyzed for pH and nutrients. The interaction effect of logging and grazing was significant(p≤0.001) only on species diversity. Relative to ungrazed areas, species diversity was significantly higher(0.01≤p≤0.05) in medium grazed gaps. Under grazed conditions, soil P was negatively correlated with gap size and species diversity. While species diversity was positivelycorrelated(0.01≤p≤0.05) with soil N in grazed plots species richness was positively correlated(0.001≤p≤0.01) with soil N in ungrazed plots. Relative density of Yushania microphylla and Carex nubigena were higher under ungrazed conditions. Our study suggests that the combined effect of cattle grazing and logging results in higher species diversity of understory vegetation in medium and grazed gaps in mixed conifer forests of Bhutan,whereas increase or decrease in relative density of major species is determined primarily by the independent effects of grazing and logging. From management perspective, forest managers must refrain from creating large gaps to avoid loss of nutrients(mainly P and N), which may eventually affect tree regeneration. Managers intending to maintain understory vegetation diversity must consider the combined effects of grazing and logging, ensuring low to moderate grazing pressure.     

Non-parametric Bayesian networks for parameter estimation in reservoir simulation: a graphical take on the ensemble Kalman filter (part I)

Reservoir simulation models are used both in the development of new fields and in developed fields where production forecasts are needed for investment decisions. When simulating a reservoir, one must account for the physical and chemical processes taking place in the subsurface. Rock and fluid properties are crucial when describing the flow in porous media. In this paper, the authors are concerned with estimating the permeability field of a reservoir. The problem of estimating model parameters such as permeability is often referred to as a history-matching problem in reservoir engineering. Currently, one of the most widely used methodologies which address the history-matching problem is the ensemble Kalman filter (EnKF). EnKF is a Monte Carlo implementation of the Bayesian update problem. Nevertheless, the EnKF methodology has certain limitations that encourage the search for an alternative method.For this reason, a new approach based on graphical models is proposed and studied. In particular, the graphical model chosen for this purpose is a dynamic non-parametric Bayesian network (NPBN). This is the first attempt to approach a history-matching problem in reservoir simulation using a NPBN-based method. A two-phase, two-dimensional flow model was implemented for a synthetic reservoir simulation exercise, and initial results are shown. The methods’ performances are evaluated and compared. This paper features a completely novel approach to history matching and constitutes only the first part (part I) of a more detailed investigation. For these reasons (novelty and incompleteness), many questions are left open and a number of recommendations are formulated, to be investigated in part II of the same paper.     

Chemical equilibrium and mass balance relationships associated with the Long Valley hydrothermal system, California, U.S.A.

Recent drilling and sampling of hydrothermal fluids from Long Valley permit an accurate characterization of chemical concentrations and equilibrium conditions in the hydrothermal reservoir. Hydrothermal fluids are thermodynamically saturated with secondary quartz, calcite, and pyrite but are in disequilibrium with respect to aqueous sulfide-sulfate speciation. Hydrothermal fluids are enriched in ¹⁸O by approximately 1‰ relative to recharge waters. ¹⁸O and Cl concentrations in well cuttings and core from high-temperature zones of the reservoir are extensively depleted relative to fresh rhyolitic tuff compositions. Approximately 80% of the Li and 50% of the B are retained in the altered reservoir rock. Cl mass balance and open-system ¹⁸O fractionation models produce similar water-rock ratios of between 1.0 and 2.5 kg kg⁻¹. These water-rock ratios coupled with estimates of reservoir porosity and density produce a minimum fluid residence time of 1.3 ka. The low fluid Cl concentrations in Long Valley correlate with corresponding low rock concentrations. Mass balance calculations indicate that leaching of these reservoir rocks accounts for Cl losses during hydrothermal activity over the last 40 ka.     

Mode of occurrence of Cu and Zn in till over three mineralized areas in Finland

The variation of Cu and Zn contents with grain size was studied in the size fractions: 2000-500 μm, 500-64 μm, 64-2 μm and minus 2 μm. The fixation to different minerals was investigated in the 500-64 μm and minus 64 μm fractions with microscope, microprobe, electron microscope and X-ray diffraction after heavy-liquid separation.The distribution of Cu and Zn between the four size fractions is very similar in all target areas. The contents tend to be highest in the minus 2 μm fraction, and to decrease with increasing grain size except for an increase in the coarsest material. This feature is typical in mineralized areas with till having anomalously high Cu and Zn contents.The anomalous Cu and Zn in the till of Talvivaara are mainly fixed to chalcopyrite and sphalerite, respectively, which indicates that the main part of the material originates from fresh bedrock and that postglacial alterations are slight. In the till of Outokumpu and Maaselkä sulphides are rare. We assume that in Outokumpu the reason for this is the postglacial weathering of sulphides and that in Maaselkä a great deal of till derives from highly weathered, sulphide-poor, preglacial regolith. In both target areas Cu in till seems to fix to limonite and also to clay minerals, particularly in Maaselkä. Zinc seems mainly to be hosted by limonite in Outokumpu.     

Derivation Of Refractive Index And Temperature Gradients From Optical Scintillometry To Correct Atmospherically Induced Errors For Highly Precise Geodetic Measurements

Refraction effects of optical beams are generally caused by an inhomogeneous propagation medium and are a major source of systematic errors in the precise optical determination of angles and distances in the atmospheric surface layer. In this contribution a method for deriving vertical temperature and refractive index gradients from optical scintillation is presented. Knowledge of these gradients is required for the compensation of atmospherically induced errors for highly precise terrestrial geodetic measurements, like direct transfer and levelling. The advantage of the present optical method is, that temperature and refractive index gradients can be derived as line-averaged values over the propagation path, which is not possible by meteorological point measurements. Field observations have been carried out with a displaced-beam scintillometer over flat terrain and under different atmospheric conditions in order to verify this method. The experiments show, that this method allows to derive accurate correction values for precise terrestrial geodetic measurements.     

Enforcing smoothness and assessing uncertainty in non-linear one-dimensional prestack seismic inversion

Estimation of elastic properties of rock formations from surface seismic amplitude measurements remains a subject of interest for the exploration and development of hydrocarbon reservoirs. This paper develops a global inversion technique to estimate and appraise 1D distributions of compressional‐wave velocity, shear‐wave velocity and bulk density, from normal‐moveout‐corrected PP prestack surface seismic amplitude measurements. Specific objectives are: (a) to evaluate the efficiency of the minimization algorithm (b) to appraise the impact of various data misfit functions, and (c) to assess the effect of the degree and type of smoothness criterion enforced by the inversion. Numerical experiments show that very fast simulated annealing is the most efficient minimization technique among alternative approaches considered for global inversion. It is also found that an adequate choice of data misfit function is necessary for a reliable and efficient match of noisy and sparse seismic amplitude measurements. Several procedures are considered to enforce smoothness of the estimated 1D distributions of elastic parameters, including predefined quadratic measures of length, flatness and roughness. Based on the general analysis of global inversion techniques, we introduce a new stochastic inversion algorithm that initializes the search for the minimum with constrained random distributions of elastic parameters and enforces predefined autocorrelation functions (semivariograms). This strategy readily lends itself to the assessment of model uncertainty. The new global inversion algorithm is successfully tested on noisy synthetic amplitude data. Moreover, we present a feasibility analysis of the resolution and uncertainty of prestack seismic amplitude data to infer 1D distributions of elastic parameters measured with wireline logs in the deepwater Gulf of Mexico. The new global inversion algorithm is computationally more efficient than the alternative global inversion procedures considered here.