Computing derivative information of sequentially coupled subsurface models

A generic framework for the computation of derivative information required for gradient-based optimization using sequentially coupled subsurface simulation models is presented. The proposed approach allows for the computation of any derivative information with no modification of the mathematical framework. It only requires the forward model Jacobians and the objective function to be appropriately defined. The flexibility of the framework is demonstrated by its application in different reservoir management studies. The performance of the gradient computation strategy is demonstrated in a synthetic water-flooding model, where the forward model is constructed based on a sequentially coupled flow-transport system. The methodology is illustrated for a synthetic model, with different types of applications of data assimilation and life-cycle optimization. Results are compared with the classical fully coupled (FIM) forward simulation. Based on the presented numerical examples, it is demonstrated how, without any modifications of the basic framework, the solution of gradient-based optimization models can be obtained for any given set of coupled equations. The sequential derivative computation methods deliver similar results compared to FIM methods, while being computationally more efficient.     

Precursory slope distress prior to the 2010 Mount Meager landslide,British Columbia

In 2010, the south flank of Mount Meager failed catastrophically, generating the largest (53 ± 3.8 × 10⁶ m³) landslide in Canadian history. We document the slow deformation of the edifice prior to failure using archival historic aerial photographs spanning the period 1948–2006. All photos were processed using Structure from Motion (SfM) photogrammetry. We used the SfM products to produce pre-and post-failure geomorphic maps that document changes in the volcanic edifice and Capricorn Glacier at its base. The photographic dataset shows that the Capricorn Glacier re-advanced from a retracted position in the 1980s then rapidly retreated in the lead-up to the 2010 failure. The dataset also documents 60 years of progressive development of faults, toe bulging, and precursory failures in 1998 and 2009. The 2010 collapse was conditioned by glacial retreat and triggered by hot summer weather that caused ice and snow to melt. Meltwater increased pore water pressures in colluvium and fractured rocks at the base of the slope, causing those materials to mobilize, which in turn triggered several secondary failures structurally controlled by lithology and faults. The landslide retrogressed from the base of the slope to near the peak of Mount Meager involving basement rock and the overlying volcanic sequence. Elsewhere on the flanks of Mount Meager, large fractures have developed in recently deglaciated areas, conditioning these slopes for future collapse. Potential failures in these areas have larger volumes than the 2010 landslide. Anticipated atmospheric warming over the next several decades will cause further loss of snow and glacier ice, likely producing additional slope instability. Satellite- and ground-based monitoring of these slopes can provide advanced warning of future landslides to help reduce risk in populated regions downstream.     

Hydrogeochemical Features of Thermal Waters of South Trungbo (Central Vietnam)

The results of studying the features of the hydrogeological structure and chemical and isotope composition of thermal waters from the central part of Vietnam that are characterized by intense manifestations of intrusive magmatism are presented. It is established that low–and high–thermal waters with temperature varying within 30–85°C are developed in the area under study. The value of total mineralization of the hydrotherms ranges from 0.05 to 10.05 g/dm³. It is assumed that the circulation of thermal waters that are different in temperature and chemical composition occurs at two levels. The regular change of the hydrotherm composition in the direction from mineralized chloride sodium, including with increased Ca content, to fresh sodium bicarbonate is revealed. The ratio of δ¹⁸O–δ²H isotopes indicates that the water component is based on meteoric water. In the coastal areas, there is an isotope shift towards the ocean waters, which is also confirmed by the hydrogeochemical data. The key factors for forming the chemical composition of the thermal waters in South Trungbo are their genetic type, the interaction processes in the “water–rock–gas–organic substance” system, and their equilibrium–nonequilibrium state.     

Ecotoxicological evaluation of dredged sediments from Rodrigo de Freitas Lagoon (Rio de Janeiro State,Brazil) using bioassays with earthworms and collembolans

The toxicity and metal bioavailability were studied in dredged sediments from Rodrigo de Freitas Lagoon (Rio de Janeiro, RJ, Brazil) using acute and avoidance tests with Eisenia andrei, and reproduction tests with Folsomia candida. The sediment was mixed with an artificial soil, and two natural soils (ferralsol and chernosol—representative Brazilian tropical soils) to obtain the following doses: 1, 3, 6, 12, 18, 24 and 30%. Total metal concentrations were determined in the sediment to support the interpretation of ecotoxicological data. Metal concentrations in the mixtures were in agreement with the threshold limits established by Brazilian law. However, significant avoidance responses were found on doses ≥?3% and were the most sensitive endpoint. Earthworm mortality found in artificial soil mixtures (LC50?=?3.9) suggests higher toxicity levels than those obtained in ferralsol (LC50?=?7.6%) and chernosol (11.0%) treatments. Earthworm mortality, avoidance responses and collembolan reproduction levels found in ferralsol mixtures (LC50?=?9.2; avoidance EC50?=?2.3%; reproduction EC50?=?2.8%) were higher compared to chernosol treatments (LC50?=?11.0%; avoidance EC50?=?4.3%; reproduction EC50?=?4.9%). The reduction of toxicity levels in chernosol mixtures is probably due to the abundance of expansive clay minerals in chernosols with capacity of adsorbing metals and other xenobiotic substances from soil pore water, decreasing metal bioavailability. Finally, threshold limits defined by Brazilian legislation for soil quality and land disposal of dredged sediments are not sufficient to prevent noxious effects on soil fauna and should be complemented with a preliminary ecotoxicological evaluation.     

Petrogenesis of low-δ<Superscript>18</Superscript>O quartz porphyry dykes,Koegel Fontein complex,South Africa

This paper investigates the origin of low-δ¹⁸O quartz porphyry dykes associated with the 144–133 Ma Koegel Fontein Igneous Complex, which was intruded during the initial phase of breakup of Africa and South America. The 25-km diameter Rietpoort Granite is the largest and youngest phase of activity, and is roofed by a 10-km diameter pendant of gneiss. Quartz porphyry (QP) dykes, up to 15 m in width, strike NW–SE across the complex. The QP dykes that intruded outside the granite have similar quartz phenocryst δ¹⁸O values (average 8.0‰, ± 0.7, n?=?33) to the granite (average 8.3?±?1.0, n?=?7). The QP dykes that intruded the roof pendant have quartz phenocrysts with more variable δ¹⁸O values (average 1.6‰, ± 2.1, n?=?55). In some cases quartz phenocrysts have δ¹⁸O values as low as ? 2.5‰. The variation in δ¹⁸O value within the quartz crystal population of individual dykes is small relative to the overall range, and core and rim material from individual quartz phenocrysts in three samples are identical within error. There is no evidence that quartz phenocryst δ¹⁸O values have been affected by fluid–rock interaction. Based on a ?_quartz?magma value of 0.6‰, magma δ¹⁸O values must have been as low as ? 3.1‰. Samples collected along the length of the two main QP dykes that traverse the roof pendant have quartz phenocryst δ¹⁸O values that range from +?1.1 to +?4.6‰, and ? 2.3 to +?5.6‰, respectively. These δ¹⁸O values correlate negatively (r = ? 0.96) with initial ⁸⁷Sr/⁸⁶Sr, which can be explained by the event that lowered δ¹⁸O values of the source being older than the dykes. We suggest that the QP dykes were fed by magma produced by partial melting of gneiss, which had been variably altered at high temperature by ¹⁸O-depleted meteoric water during global glaciation at ~?550 Ma. The early melts had variable δ¹⁸O value but as melt pockets interconnected during melting, the δ¹⁸O values approached that of average gneiss. Variable quartz phenocryst δ¹⁸O values in the same dyke can be explained by vertical emplacement, at variable rates of ascent along the dyke. The lateral variation in quartz, and hence magma δ¹⁸O value at a particular point along a single dyke would depend on the rate of ascent of magma at that point along the dyke, and the ‘age’ of the particular magma batch.     

Crystal nucleation and growth produced by continuous decompression of Pinatubo magma

High-temperature decompression experiments demonstrate that crystal textures preserve a record of the style and rate of magmatic ascent. To reinforce this link, we performed a suite of isothermal decompression experiments using starting material from the climactic 1991 Pinatubo eruption. We decompressed experiments from 220 MPa to final, quench pressures of 75 or 30 MPa using continuous decompression rates of 100, 30, 10, 3, 1, and 0.3 MPa h^?1. Amphibole, clinopyroxene, and plagioclase crystallized during the experiments, with plagioclase microlites dominating the assemblage. Total microlite number densities range from 10^7.6±0.4 up to 10^8.2±0.2 cm^?3, with plagioclase accounting for up to 65% of the total number. Plagioclase microlite area increased systematically from 19?±?8 to 937?±?487 µm² with increasing experiment duration. Our textures provide time-integrated records of crystal kinetics. Average nucleation and areal growth rates of plagioclase are highest in the fastest decompressions (~?10^7.5 cm^?3 h^?1 and 10.1?±?4.1 µm² h^?1, respectively) and more than an order of magnitude lower in the slowest experiments (~?10^5.5 cm^?3 h^?1 and 0.8?±?0.2 µm² h^?1, respectively). Both nucleation and growth rates are highest at high degrees of disequilibrium. We find that peak supersaturation-dependent instantaneous rates are generally more than an order of magnitude faster than average rates. We use those instantaneous nucleation and growth rates to introduce an iterative model to evaluate the effects of different decompression rates, decompression paths (continuous, single-step or multistep), and the presence of phenocrysts on final crystallinity and microlite size distribution.     

Zoning and exsolution in alkali feldspars from Laacher See volcano (Western Germany): constraints on temperature history prior to eruption

Compositional zoning and exsolution patterns of alkali feldspars in carbonatite-bearing cognate syenites from the 6.3 km³ (D.R.E) phonolitic Laacher See Tephra (LST) deposit in western Germany (12.9 ka) are reported. These rocks represent the cooler outer portion and crystal-rich products of a cooling magma reservoir at upper crustal levels. Major and trace-element difference between cores and rims in sanidine crystals represent two generations of crystal growth separated by unmixing of a carbonate melt. Trace-element differences measured by LA–ICP–MS are in accordance with silicate–carbonate unmixing. Across the core–rim boundary, we extracted gray-scale profiles from multiple accumulations of back-scattered electron images. Gray scales directly represent K/Na ratios owing to low concentrations of Ba and Sr (<?30 ppm). Diffusion gradients are modeled to solve for temperature using known pre-eruptive U–Th zircon ages (0–20 ky) of each sample (Schmitt et al., J Petrol 51:1053–1085, 2010). Estimated temperatures range from 630 °C to 670 °C. For the exsolution boundaries, a diffusive homogenization model is constrained by the solvus temperature of ~ 712^_725 °C and gives short time scales of only 15–50 days. Based on our results, we present a model for the temperature–time history of these rocks. The model also constrains the thermal variation across the cooling crystal-rich carapace of the magma reservoir over 20 ka and suggests a thermal reactivation of cumulates, the cooling carapace, and probably the entire system only a few years prior to the explosive eruption of the remaining molten core of the phonolitic magma reservoir.     

Estimation of commercial diamond grades based on microdiamonds: a case study of the Koidu diamond mine,Sierra Leone

This paper documents the application of a microdiamond-based approach to the estimation of diamond grade in the Pipe 1 kimberlite at the Koidu mine in Sierra Leone. A geological model of Pipe 1 was constructed to represent the distribution and volume of the dominant kimberlite units within the pipe. Bulk samples, along with representative microdiamond samples, were collected from these units at surface and were used to define the ratio between microdiamond stone frequency (+212 μm stones per kilogram) and recoverable macrodiamond grade (+1.2 mm carats per tonne; 1 carat = 0.2 g). These ratios were applied to a comprehensive, spatially representative microdiamond sample dataset and were combined with a spatial model of country-rock xenolith dilution within the pipe to estimate +1.2 mm recoverable grades. The resource estimate was reconciled with subsequent production results in the elevation range 160 to 100 m above sea level. Production results for each of the six 10 m benches covering this elevation range were compared to the estimated average grades for these zones in the pipe. For the five cases where most of the kimberlite mass on a given bench is represented in the production data, the results show a maximum discrepancy of 6% between predicted and reported production grade with no indication of any consistent bias. This indicates that, when supported by a sound geological model and suitable microdiamond and macrodiamond data, the microdiamond-based estimation approach can provide reliable constraints on macrodiamond grade, even in the case of geologically complex bodies such as Koidu Pipe 1.