Numerical analysis of wave-induced fluid flow effects on seismic data: Application to monitoring of CO2 storage at the Sleipner Field

In this work we analyze how patchy distributions of CO2 and brine within sand reservoirs may lead to significant attenuation and velocity dispersion effects, which in turn may have a profound impact on surface seismic data. The ultimate goal of this paper is to contribute to the understanding of these processes within the framework of the seismic monitoring of CO2 sequestration, a key strategy to mitigate global warming. We first carry out a Monte Carlo analysis to study the statistical behavior of attenuation and velocity dispersion of compressional waves traveling through rocks with properties similar to those at the Utsira Sand, Sleipner field, containing quasi-fractal patchy distributions of CO2 and brine. These results show that the mean patch size and CO2 saturation play key roles in the observed wave-induced fluid flow effects. The latter can be remarkably important when CO2 concentrations are low and mean patch sizes are relatively large. To analyze these effects on the corresponding surface seismic data, we perform numerical simulations of wave propagation considering reservoir models and CO2 accumulation patterns similar to the CO2 injection site in the Sleipner field. These numerical experiments suggest that wave-induced fluid flow effects may produce changes in the reservoir's seismic response, modifying significantly the main seismic attributes usually employed in the characterization of these environments. Consequently, the determination of the nature of the fluid distributions as well as the proper modeling of the seismic data constitute important aspects that should not be ignored in the seismic monitoring of CO2 sequestration problems.

Reduction of symmetry in grandite solid solution

The results of a crystal structure refinement of an anisotropic grandite garnet specimen with composition Gro36-4 And63-6 are given. The structure obtained has orthorrombic symmetry (space group Fddd) and is compared with similar results obtained by other authors. In all cases the reduction of symmetry is due to the ordering of Fe3+ and Al in octahedral sites. Non cubic structures of grandites are discussed in connection with optical, morphological an grou-th features of these minerals.

Numerical analysis of the impact of charcoal production on soil hydrological behavior, runoff response and erosion susceptibility

The impact of charcoal production on soil hydraulic properties, runoff response and erosion susceptibility were studied in both field and simulation experiments. Core and composite samples, from 12 randomly selected sites within the catchment of Kotokosu were taken from the 0-10 cm layer of a charcoal site soil (CSS) and adjacent field soils (AFS). These samples were used to determine saturated hydraulic conductivity (Ksat), bulk density, total porosity, soil texture and color. Infiltration, surface albedo and soil surface temperature were also measured in both CSS and AFS. Measured properties were used as entries in a rainfall runoff simulation experiment on a smooth (5 % slope) plot of 25 x 25 m grids with 10 cm resolutions. Typical rainfall intensities of the study watershed (high, moderate and low) were applied to five different combinations of Ks distributions that could be expected in this landscape. The results showed significantly (p < 0.01) higher flow characteristics of the soil under charcoal kilns (increase of 88 %). Infiltration was enhanced and runoff volume reduced significantly. The results showed runoff reduction of about 37 and 18 %, and runoff coefficient ranging from 0.47-0.75 and 0.04-0.39 or simulation based on high (200 mm h-1) and moderate (100 mm h-1) rainfall events over the CSS and AFS areas, respectively. Other potential impacts of charcoal production on watershed hydrology were described. The results presented, together with watershed measurements, when available, are expected to enhance understanding of the hydrological responses of ecosystems to indiscriminate charcoal production and related activities in this region.

Bootstrapping pairs in Distance-Based Regression

La regressió basada en distàncies és un mètode de predicció que consisteix en dos passos: a partir de les distàncies entre observacions obtenim les variables latents, les quals passen a ser els regressors en un model lineal de mínims quadrats ordinaris. Les distàncies les calculem a partir dels predictors originals fent us d'una funció de dissimilaritats adequada. Donat que, en general, els regressors estan relacionats de manera no lineal amb la resposta, la seva selecció amb el test F usual no és possible. En aquest treball proposem una solució a aquest problema de selecció de predictors definint tests estadístics generalitzats i adaptant un mètode de bootstrap no paramètric per a l'estimació dels p-valors. Incluim un exemple numèric amb dades de l'assegurança d'automòbils.

Simultaneous multielement extraction with the Mehlich-1 solution for Southern Brazilian soils determined by ICP-OES and the effects on the nutrients recommendations to crops

The amounts of macro (P, K, Ca and Mg) and micronutrients (Cu and Zn) extracted with the Mehlich-1 (M1) solution, by the 1.0 mol L-1 KCl (KCl) and with the 0.1 mol L-1 HCl (HCl) for representative soil types of the Rio Grande do Sul state (Brazil) were compared with those extracted with the Mehlich-1 solution determined with the inductively coupled plasma optical emission spectroscopy (ICP). The amounts of nutrients extracted by the different methods showed high correlation coefficients. On average, the Mehlich-1 solution extracted similar amounts of P, determined with colorimetric and ICP methods, and, K determined with emission and ICP. The amounts of Ca and Mg extracted with the Mehlich-1 solution, determined by ICP, were similar to those extracted with the KCl solution determined by the atomic absorption spectrophotometry. The amounts of Cu and Zn extracted with the Mehlich-1 solution, determined by the ICP, were higher than those extracted with the 0.1 mol L-1 HCl determined by the atomic absorption spectrophotometry. The results indicate that the Mehlich-1 solution and ICP can be used for simultaneous multielement extraction and determination for Southern Brazilian soils. However, a conversion factor for values interpretation is needed. The use of the conversion factor to determine the K availability index in soils is adequate and does not affect the K recommendations for crops in southern Brazilian soils.

Nitrate and ammonium in soil solution in tobacco management systems

Tobacco farmers of southern Brazil use high levels of fertilizers, without considering soil and environmental attributes, posing great risk to water resources degradation. The objective of this study was to monitor nitrate and ammonium concentrations in the soil solution of an Entisol in and below the root zone of tobacco under conventional tillage (CT), minimum tillage (MT) and no-tillage (NT). The study was conducted in the small-watershed Arroio Lino, in Agudo, State of Rio Grande do Sul, Brazil. A base fertilization of 850 kg ha-1 of 10-18-24 and topdressing of 400 kg ha-1 of 14-0-14 NPK fertilizer were applied. The soil solution was sampled during the crop cycle with a tension lysimeter equipped with a porous ceramic cup. Ammonium and nitrate concentrations were analyzed by the distillation and titration method. Nitrate concentrations, ranging from 8 to 226 mg L-1, were highest after initial fertilization and decreased during the crop cycle. The average nitrate (N-NO3-) concentration in the root zone was 75 in NT, 95 in MT, and 49 mg L-1 in CT. Below the root zone, the average nitrate concentration was 58 under NT, 108 under MT and 36 mg L-1 under CT. The nitrate and ammonium concentrations did not differ significantly in the management systems. However, the nitrate concentrations measured represent a contamination risk to groundwater of the watershed. The ammonium concentration (N-NH4+) decreased over time in all management systems, possibly as a result of the nitrification process and root uptake of part of the ammonium by the growing plants.

Ruin Theory with Excess of Loss Reinsurance and Reinstatements

The present paper studies the probability of ruin of an insurer, if excess of loss reinsurance with reinstatements is applied. In the setting of the classical Cramer-Lundberg risk model, piecewise deterministic Markov processes are used to describe the free surplus process in this more general situation. It is shown that the finite-time ruin probability is both the solution of a partial integro-differential equation and the fixed point of a contractive integral operator. We exploit the latter representation to develop and implement a recursive algorithm for numerical approximation of the ruin probability that involves high-dimensional integration. Furthermore we study the behavior of the finite-time ruin probability under various levels of initial surplus and security loadings and compare the efficiency of the numerical algorithm with the computational alternative of stochastic simulation of the risk process. (C) 2011 Elsevier Inc. All rights reserved.

Attributes of irrigated rice and soil solution as affected by salinity levels of the water layer

In the State of Rio Grande do Sul, Brazil, flooded rice fields using Patos Lagoon as the source of water for irrigation are subject to be damaged by salinity, since this source is bound to the sea on its southern end. The sensitivity of rice is variable during plant development, being higher in the seedling and reproductive periods. However, there is not enough information about the behavior of plants under salt stress during the course of its development, especially in the vegetative stage. This study evaluated the effect of different levels of salinity of irrigation water on the salinity of soil solution over time and on some plant attributes, during the vegetative stage of rice. The study was conducted in a greenhouse, where seeds of the variety IRGA 424 were sown in pots and irrigated with water with electrical conductivity (ECi) levels of: 0.3, 0.75, 1.5, 3.0 and 4.5 dS m-1; from the tillering initiation (V4) until the panicle initiation (PI). The evaluations made were the electrical conductiviy of soil solution (ECe), the dry biomass of plants and stems, tillering, height and the transpiration of plants. The ECe increased with the ECi over time, and was determined by water transpiration flux in pots. The ECe values at the end of the experiment were high and, in most cases, higher than the critical values for flooded rice. The growth attributes of rice were negatively affected from ECi of 2.0 dS m-1 and ECe of 4.0 dS m-1.

Rice grown in nutrient solution with doses of manganese and silicon

Although silicon is not recognized as a nutrient, it may benefit rice plants and may alleviate the Mn toxicity in some plant species. The dry matter yield (root, leaf, sheaths and leaf blade) and plant architecture (angle of leaf insertion and leaf arc) were evaluated in rice plants grown in nutrient solutions with three Mn doses, with and without Si addition. The treatments were arranged in a 2 x 3 factorial [with and without (2 mmol L-1) Si; three Mn doses (0.5; 2.5 and 10 µmol L-1)], in a randomized block design with 4 replications. The experimental unit was a 4 L plastic vase with 4 rice (Metica-1 cultivar) plants. Thirty nine days after keeping the seedlings in the nutrient solution the plant dry matter yield was determined; the angle of leaf insertion in the sheath and the leaf arc were measured; and the Si and Mn concentrations in roots, sheaths and leaves were determined. The analysis of variance (F test at 5 and 1 % levels) and the regression analysis (for testing plant response to Mn with the Si treatments) were performed. The Si added to the nutrient solution increased the dry matter yield of roots, sheaths and leaf blades and also decreased the angle of leaf blade insertion into the sheath and the foliar arc in the rice plant. Additionally, it ameliorated the rice plant architecture which allowed an increase in the dry matter yield. Similarly, the addition of Mn to the solution improved the architecture of the rice plants with gain in dry matter yield. As Si was added to the nutrient solution, the concentration of Mn in leaves decreased and in roots increased thus alleviating the toxic effects of Mn on the plants.

Numerical algorithm for Ginzburg-Landau equations with multiplicative noise: Application to domain growth

We consider stochastic partial differential equations with multiplicative noise. We derive an algorithm for the computer simulation of these equations. The algorithm is applied to study domain growth of a model with a conserved order parameter. The numerical results corroborate previous analytical predictions obtained by linear analysis.

Characterization of gravitational rock slope deformations at different spatial scales based on field, remote sensing and numerical approaches

Rock slope instabilities such as rock slides, rock avalanche or deep-seated gravitational slope deformations are widespread in Alpine valleys. These phenomena represent at the same time a main factor that control the mountain belts erosion and also a significant natural hazard that creates important losses to the mountain communities. However, the potential geometrical and dynamic connections linking outcrop and slope-scale instabilities are often unknown. A more detailed definition of the potential links will be essential to improve the comprehension of the destabilization processes and to dispose of a more complete hazard characterization of the rock instabilities at different spatial scales. In order to propose an integrated approach in the study of the rock slope instabilities, three main themes were analysed in this PhD thesis: (1) the inventory and the spatial distribution of rock slope deformations at regional scale and their influence on the landscape evolution, (2) the influence of brittle and ductile tectonic structures on rock slope instabilities development and (3) the characterization of hazard posed by potential rock slope instabilities through the development of conceptual instability models. To prose and integrated approach for the analyses of these topics, several techniques were adopted. In particular, high resolution digital elevation models revealed to be fundamental tools that were employed during the different stages of the rock slope instability assessment. A special attention was spent in the application of digital elevation model for detailed geometrical modelling of past and potential instabilities and for the rock slope monitoring at different spatial scales. Detailed field analyses and numerical models were performed to complete and verify the remote sensing approach. In the first part of this thesis, large slope instabilities in Rhone valley (Switzerland) were mapped in order to dispose of a first overview of tectonic and climatic factors influencing their distribution and their characteristics. Our analyses demonstrate the key influence of neotectonic activity and the glacial conditioning on the spatial distribution of the rock slope deformations. Besides, the volumes of rock instabilities identified along the main Rhone valley, were then used to propose the first estimate of the postglacial denudation and filling of the Rhone valley associated to large gravitational movements. In the second part of the thesis, detailed structural analyses of the Frank slide and the Sierre rock avalanche were performed to characterize the influence of brittle and ductile tectonic structures on the geometry and on the failure mechanism of large instabilities. Our observations indicated that the geometric characteristics and the variation of the rock mass quality associated to ductile tectonic structures, that are often ignored landslide study, represent important factors that can drastically influence the extension and the failure mechanism of rock slope instabilities. In the last part of the thesis, the failure mechanisms and the hazard associated to five potential instabilities were analysed in detail. These case studies clearly highlighted the importance to incorporate different analyses and monitoring techniques to dispose of reliable and hazard scenarios. This information associated to the development of a conceptual instability model represents the primary data for an integrated risk management of rock slope instabilities. - Les mouvements de versant tels que les chutes de blocs, les éboulements ou encore les phénomènes plus lents comme les déformations gravitaires profondes de versant représentent des manifestations courantes en régions montagneuses. Les mouvements de versant sont à la fois un des facteurs principaux contrôlant la destruction progressive des chaines orogéniques mais aussi un danger naturel concret qui peut provoquer des dommages importants. Pourtant, les phénomènes gravitaires sont rarement analysés dans leur globalité et les rapports géométriques et mécaniques qui lient les instabilités à l'échelle du versant aux instabilités locales restent encore mal définis. Une meilleure caractérisation de ces liens pourrait pourtant représenter un apport substantiel dans la compréhension des processus de déstabilisation des versants et améliorer la caractérisation des dangers gravitaires à toutes les échelles spatiales. Dans le but de proposer un approche plus globale à la problématique des mouvements gravitaires, ce travail de thèse propose trois axes de recherche principaux: (1) l'inventaire et l'analyse de la distribution spatiale des grandes instabilités rocheuses à l'échelle régionale, (2) l'analyse des structures tectoniques cassantes et ductiles en relation avec les mécanismes de rupture des grandes instabilités rocheuses et (3) la caractérisation des aléas rocheux par une approche multidisciplinaire visant à développer un modèle conceptuel de l'instabilité et une meilleure appréciation du danger . Pour analyser les différentes problématiques traitées dans cette thèse, différentes techniques ont été utilisées. En particulier, le modèle numérique de terrain s'est révélé être un outil indispensable pour la majorité des analyses effectuées, en partant de l'identification de l'instabilité jusqu'au suivi des mouvements. Les analyses de terrain et des modélisations numériques ont ensuite permis de compléter les informations issues du modèle numérique de terrain. Dans la première partie de cette thèse, les mouvements gravitaires rocheux dans la vallée du Rhône (Suisse) ont été cartographiés pour étudier leur répartition en fonction des variables géologiques et morphologiques régionales. En particulier, les analyses ont mis en évidence l'influence de l'activité néotectonique et des phases glaciaires sur la distribution des zones à forte densité d'instabilités rocheuses. Les volumes des instabilités rocheuses identifiées le long de la vallée principale ont été ensuite utilisés pour estimer le taux de dénudations postglaciaire et le remplissage de la vallée du Rhône lié aux grands mouvements gravitaires. Dans la deuxième partie, l'étude de l'agencement structural des avalanches rocheuses de Sierre (Suisse) et de Frank (Canada) a permis de mieux caractériser l'influence passive des structures tectoniques sur la géométrie des instabilités. En particulier, les structures issues d'une tectonique ductile, souvent ignorées dans l'étude des instabilités gravitaires, ont été identifiées comme des structures très importantes qui contrôlent les mécanismes de rupture des instabilités à différentes échelles. Dans la dernière partie de la thèse, cinq instabilités rocheuses différentes ont été étudiées par une approche multidisciplinaire visant à mieux caractériser l'aléa et à développer un modèle conceptuel trois dimensionnel de ces instabilités. A l'aide de ces analyses on a pu mettre en évidence la nécessité d'incorporer différentes techniques d'analyses et de surveillance pour une gestion plus objective du risque associée aux grandes instabilités rocheuses.

Effect of NPK fertilization on production and leaf nutrient content of eucalyptus minicuttings in nutrient solution

Adequate nutrient levels in plants vary according to the species or clone, age and management practice. Therefore, adjustments of the nutrient solution are often necessary according to the plant material for multiplication. This study aimed to evaluate the influence of NPK fertilization on production and leaf nutrient contents of eucalyptus cuttings in nutrient solution. The study was conducted from November 2008 to January 2009 in a greenhouse. The experimental design was completely randomized fractional factorial (4 x 4 x 4)½, with a total of 32 treatments with three replications. The treatments consisted of four doses of N (50, 100, 200 and 400 mg L-1) as urea, P (7.5, 15, 30 and 60 mg L-1) in the form of phosphoric acid and K (50, 100, 200 and 400 mg L-1) in the form of potassium chloride in the nutrient solution. Only the effect of N alone was significant for the number and dry weight of minicuttings per ministump, with a linear decreasing effect with increasing N levels. The highest number of cuttings was obtained at a dose of 50, 7.5 and 50 mg L-1 of N, P and K, respectively.

Numerical modeling of the effect of variation of boundary conditions on vadose zone hydraulic properties

An accurate estimation of hydraulic fluxes in the vadose zone is essential for the prediction of water, nutrient and contaminant transport in natural systems. The objective of this study was to simulate the effect of variation of boundary conditions on the estimation of hydraulic properties (i.e. water content, effective unsaturated hydraulic conductivity and hydraulic flux) in a one-dimensional unsaturated flow model domain. Unsaturated one-dimensional vertical water flow was simulated in a pure phase clay loam profile and in clay loam interlayered with silt loam distributed according to the third iteration of the Cantor Bar fractal object Simulations were performed using the numerical model Hydrus 1D. The upper and lower pressure heads were varied around average values of -55 cm for the near-saturation range. This resulted in combinations for the upper and lower constant head boundary conditions, respectively, of -50 and -60 cm, -40 and -70 cm, -30 and -80 cm, -20 and -90 cm, and -10 and -100 cm. For the drier range the average head between the upper and lower boundary conditions was set to -550 cm, resulting in the combinations -500 and -600 cm, -400 and -700 cm, -300 and -800 cm, -200 and -900 cm, and -100 and -1,000 cm, for upper and lower boundary conditions, respectively. There was an increase in water contents, fluxes and hydraulic conductivities with the increase in head difference between boundary conditions. Variation in boundary conditions in the pure phase and interlayered one-dimensional profiles caused significant deviations in fluxes, water contents and hydraulic conductivities compared to the simplest case (a head difference between the upper and lower constant head boundaries of 10 cm in the wetter range and 100 cm in the drier range).

Zenith angle distributions at Super-Kamiokande and SNO and the solution of the solar neutrino problem

We have performed a detailed study of the zenith angle dependence of the regeneration factor and distributions of events at SNO and SK for different solutions of the solar neutrino problem. In particular, we discuss the oscillatory behavior and the synchronization effect in the distribution for the LMA solution, the parametric peak for the LOW solution, etc. A physical interpretation of the effects is given. We suggest a new binning of events which emphasizes the distinctive features of the zenith angle distributions for the different solutions. We also find the correlations between the integrated day-night asymmetry and the rates of events in different zenith angle bins. The study of these correlations strengthens the identification power of the analysis.