The MARID (mica-amphibole-rutile-ilmenite-diopside) suite of xenoliths in kimberlite

Within the ‘glimmerite’ nodules occurring within kimberlite pipes we recognize the MARID suite consisting of varying proportions of mica, amphibole, rutile, ilmenite and diopside. Banding of some specimens is interpreted as cumulate layering. All specimens were deformed either before incorporation into the host kimberlite or during intrusion. Compared with minerals in peridotite xenoliths, the MARID ones are lower in Al₂O₃ and Cr₂O₃, but richer in total iron. The MARID micas, amphiboles, diopsides, ilmenites and probably rutiles contain substantial Fe₂O₃ indicative of oxidizing conditions. The amphibole is potassic richterite. Micas of the megacryst suite in kimberlite have less total iron and Fe₂O₃ than micas of the MARID suite. We suggest that the rocks of the MARID suite crystallized under oxidizing conditions from a magma, chemically similar to kimberlite, within the higher parts of the upper mantle: the presence of amphibole restricts the depth to less than ~ 100 km. A xenolith containing olivine and orthopyroxene as well as minerals similar to but not the same compositionally as MARID-types is interpreted as a metasomite, possibly representing wall-rock of a magma body from which MARID-suite rocks crystallized.     

A numerical approach to study the properties of solutions of the diffusive wave approximation of the shallow water equations

In this paper, we study the properties of approximate solutions to a doubly nonlinear and degenerate diffusion equation, known in the literature as the diffusive wave approximation of the shallow water equations (DSW), using a numerical approach based on the Galerkin finite element method. This equation arises in shallow water flow models when special assumptions are used to simplify the shallow water equations and contains as particular cases the porous medium equation and the p-Laplacian. Diverse numerical schemes have been implemented to approximately solve the DSW equation and have been successfully applied as suitable models to simulate overland flow and water flow in vegetated areas such as wetlands; yet, no formal mathematical analysis has been carried out in order to study the properties of approximate solutions. In this study, we propose a numerical approach as a means to understand some properties of solutions to the DSW equation and, thus, to provide conditions for which the use of the DSW equation may be inappropriate from both the physical and the mathematical points of view, within the context of shallow water modeling. For analysis purposes, we propose a numerical method based on the Galerkin method and we obtain a priori error estimates between the approximate solutions and weak solutions to the DSW equation under physically consistent assumptions. We also present some numerical experiments that provide relevant information about the accuracy of the proposed numerical method to solve the DSW equation and the applicability of the DSW equation as a model to simulate observed quantities in an experimental setting.     

National-scale Assessment of Current and Future Flood Risk in England and Wales

In recent years, through the availability of remotely sensed data and other national datasets, it has become possible to conduct national-scale flood risk assessment in England and Wales. The results of this type of risk analysis can be used to inform policy-making and prioritisation of resources for flood management. It can form the starting point for more detailed strategic and local-scale flood risk assessments. The national-scale risk assessment methodology outlined in this paper makes use of information on the location, standard of protection and condition of flood defences in England and Wales, together with datasets of floodplain extent, topography, occupancy and asset values. The flood risk assessment was applied to all of England and Wales in 2002 at which point the expected annual damage from flooding was estimated to be approximately £1 billion. This figure is comparable with records of recent flood damage. The methodology has subsequently been applied to examine the effects of climate and socio-economic change 50 and 80 years in the future. The analysis predicts increasing flood risk unless current flood management policies, practices and investment levels are changed – up to 20-fold increase in real terms economic risk by the 2080s in the scenario with highest economic growth. The increase is attributable primarily to a combination of climate change (in particular sea level rise and increasing precipitation in parts of the UK) and increasing economic vulnerability.     

Therm odynamics of Diagenetic Fluid and Fluid／Mineral Reactions in the Eogene Xingouzui Formation,Oil Field T,Jianghan Basin

This study focuses on the thermodynamics of diagenetic fluid from the Eogene Xingouzui Formation which represents the most important reservoir in Field Oil T in the Jianghan Basin. The measured homogenization temperatures (110–139 °C) of fluid inclusions in diagenetic minerals fell within the range of 67 –155 °C at the middle diagenetic stage. The pressure of diagenetic fluid is estimated at 10.2 –56 MPa. The activity of ions in the fluid shows a tendency of Ca²⁺ > Mg²⁺ > Na⁺ > K⁺ > Fe³⁺ > Fe²⁺ for cations, and HCO ₃ ⁻ > SO ₄ ²⁻ > F⁻ > Cl⁻ > CO ₃ ²⁻ for anions. For the gaseous facies, there is a tendency of CO₂> CO> H₂S> CH₄> H₂. According to the thermodynamic calculations, the pH and Eh of the fluid are 5.86–6.47 and −0.73–−0.64V, respectively. As a result of the interaction between such a diagenetic fluid and minerals in the sediments, feldspars were dissolved or alterated by other minerals. The clay mineral kaolinite was instable and hence was replaced by illite and chloritoid. This project was jointly funded by the National Natural Science Foundation of China (49133080) and the Open Laboratory of Ore Deposit Geochemistry, Institute of Geochemistry, Chinese Academy of Sciences.     

Oxidation of sulfur dioxide by ozone in highly dispersed water droplets

The oxidation of S(IV) by ozone in aqueous solution has been studied in highly dispersed droplets in a laboratory could chamber at 20°C. Urban air was used. Ozone concentrations ranged from 15 to 120 ppb, orders of magnitude lower than most previous works. Rate constants were derived from differential product determination, rather than from rates of decay of bulk reactants. Comparison of the results obtained here with those from previous work indicates that the reaction rate is indeed first order with respect to ozone. The second-order reaction rate constants found in this work lay between, or close to, other recent results, and could be represented by:k=1.23×10⁴ [H⁺]^–0.51 M^–1 s^–1.     

Heated mine room and pillar-secondary creep response

Heated salt mine room and pillar simulations have been performed to provide information regarding parameters affecting room closure rates to designers of radioactive waste isolation facilities. A coupled secondary creep and heat transfer formulation with large strain capabilities was used to assess the effects of variations in creep law parameters, thermal properties, imposed boundary conditions, temporal integration and meshing resolution on room closure rates. Results indicate that the greatest effect results from variations in parameters appearing in the creep constitutive equation.     

Climate change and river flooding: part 1 classifying the sensitivity of British catchments

Effective national and regional policy guidance on climate change adaptation relies on robust scientific evidence. This two-part series of papers develops and implements a novel scenario-neutral framework enabling an assessment of the vulnerability of flood flows in British catchments to climatic change, to underpin the development of guidance for the flood management community. In this first part, the sensitivity of the 20-year return period flood peak (RP20) to changes in precipitation (P), temperature (T) and potential evapotranspiration (PE) is systematically assessed for 154 catchments. A sensitivity domain of 4,200 scenarios is applied combining 525 and 8 sets of P and T/PE mean monthly changes, respectively, with seasonality incorporated using a single-phase harmonic function. Using the change factor method, the percentage change in RP20 associated with each scenario of the sensitivity domain is calculated, giving flood response surfaces for each catchment. Using a clustering procedure on the response surfaces, the 154 catchments are divided into nine groups: flood sensitivity types. These sensitivity types show that some catchments are (very) sensitive to changes in P but others buffer the response, while the location of catchments of the same type does not show any strong geographical pattern. These results reflect the range of hydrological processes found in Britain, and demonstrate the potential importance of catchment properties (physical and climatic) in the propagation of change in climate to change in floods, and so in characterising the sensitivity types (covered in the companion paper).     

The regularizing Levenberg–Marquardt scheme for history matching of petroleum reservoirs

In this paper, we study on a history matching approach that consists of finding stable approximations to the problem of minimizing the weighted least-squares functional that penalizes the misfit between the reservoir model predictions G(u) and noisy observations y ^η . In other words, we are interested in computing an approximation to the minimizer of $\frac {1}{2}\vert \vert \Gamma ^{-1/2}(y^{\eta }-G(u))\vert \vert _{Y}^{2} $ where Γ is the measurements error covariance, Y is the observation space, and X is a set of admissible parameters. This is an ill-posed nonlinear inverse problem that we address by means of the regularizing Levenberg–Marquardt scheme developed by Hanke (Inverse Probl. 13:79–95, 1997; J. Integr. Equ. Appl. 22(2):259–283, 2010). Under certain conditions on G, the theory of Hanke (Inverse Probl. 13:79–95, 1997; J. Integr. Equ. Appl. 22(2):259–283, 2010) ensures the convergence of the scheme to stable approximations to the inverse problem. We propose an implementation of the regularizing Levenberg–Marquardt scheme that enforces prior knowledge on the geologic properties. In particular, the prior mean $\overline {u}$ is incorporated in the initial guess of the algorithm, and the prior error covariance C is enforced through the definition of the parameter space X. Our main goal is to numerically show that the proposed implementation of the regularizing Levenberg–Marquardt scheme of Hanke is a robust method capable of providing accurate estimates of the geologic properties for small noise measurements. In addition, we provide numerical evidence of the convergence and regularizing results predicted by the theory of Hanke (Inverse Probl. 13:79–95, 1997; J. Integr. Equ. Appl. 22(2):259–283, 2010) for a prototypical oil–water reservoir model. The performance for recovering the true permeability with the regularizing Levenberg–Marquardt scheme is compared to the typical approach of computing the maximum a posteriori (MAP) estimator. In particular, we compare the proposed application of the regularizing Levenberg–Marquardt (LM) scheme against the standard LM approach of Li et al. (SPE J. 8(4):328–340, 2003) and Reynolds et al. (2008) for computing the MAP. Our numerical experiments suggest that the history matching approach based on iterative regularization is robust and could potentially be used to improve further on various methodologies already proposed as effective tools for history matching in petroleum reservoirs.     

A local discontinuous Galerkin scheme for Darcy flow with internal jumps

We present a new version of the local discontinuous Galerkin method which is capable of dealing with jump conditions along a submanifold Γ_LG (i.e., Henry’s Law) in instationary Darcy flow. Our analysis accounts for a spatially and temporally varying, non-linear permeability tensor in all estimates which is also allowed to have a jump at Γ_LG and gives a convergence order result for the primary and the flux unknowns. In addition to this, different approximation spaces for the primary and the flux unknowns are investigated. The results imply that the most efficient choice is to choose the degree of the approximation space for the flux unknowns one less than that of the primary unknown. The only stabilization in the proposed scheme is represented by a penalty term in the primary unknown.