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.     

Gravimetric and magnetic fabric study of the Sintra Igneous complex: laccolith-plug emplacement in the Western Iberian passive margin

International Journal of Earth Sciences - The geometry and emplacement of the ~ 96 km2, Late Cretaceous Sintra Igneous complex (SIC, ca. 80 Ma) into the West Iberian passive...     

The Juchatengo complex: an upper-level ophiolite assemblage of late Paleozoic age in Oaxaca,southern Mexico

International Journal of Earth Sciences - The Juchatengo complex (JC) suite is located between the Proterozoic Oaxacan complex to the north and the Xolapa complex to the south, and was amalgamated...     

Rainfall thresholds for the triggering of landslides considering previous soil moisture conditions (Asturias,NW Spain)

Given its geological and climatic conditions and its rugged orography, Asturias is one of the most landslide prone areas in the North of Spain. Most of the landslides occur during intense rainfall episodes. Thus, precipitation is considered the main triggering factor in the study area, reaching average annual values of 960 mm. Two main precipitation patterns are frequent: (i) long-lasting periods of moderate rainfall during autumn and winter and (ii) heavy short rainfall episodes during spring and early summer. In the present work, soil moisture conditions in the locations of 84 landslides are analysed during two rainfall episodes, which represent the most common precipitation patterns: October–November 2008 and June 2010. Empirical data allowed the definition of available water capacity percentages of 99–100% as critical soil moisture conditions for the landslide triggering. Intensity-duration rainfall thresholds were calculated for each episode, considering the periods with sustained high soil moisture levels before the occurrence of each analysed landslide event. For this purpose, data from daily water balance models and weather stations were used. An inverse relationship between the duration of the precipitation and its intensity, consistent with published intensity-duration thresholds, was observed, showing relevant seasonal differences.     

Apatite-hosted melt inclusions from the Panzhihua gabbroic-layered intrusion associated with a giant Fe–Ti oxide deposit in SW China: insights for magma unmixing within a crystal mush

The ～?2-km-thick Panzhihua gabbroic-layered intrusion in SW China is unusual because it hosts a giant Fe–Ti oxide deposit in its lower zone. The deposit consists of laterally extensive net-textured and massive Fe–Ti oxide ore layers, the thickest of which is ～?60 m. To examine the magmatic processes that resulted in the Fe enrichment of parental high-Ti basaltic magma and the formation of thick, Fe–Ti oxide ore layers, we carried out a detailed study of melt inclusions in apatite from a ～?500-m-thick profile of apatite-bearing leucogabbro in the middle zone of the intrusion. The apatite-hosted melt inclusions are light to dark brown in color and appear as polygonal, rounded, oval and negative crystal shapes, which range from ～?5 to ～?50 µm in width and from ～?5 to ～?100 µm in length. They have highly variable compositions and show a large and continuous range of SiO₂ and FeO_t with contrasting end-members; one end-member being Fe-rich and Si-poor (40.2 wt% FeO_t and 17.7 wt% SiO₂) and the other being Si-rich and Fe-poor (74.0 wt% SiO₂ and 1.20 wt% FeO_t). This range in composition may be attributed to entrapment of the melt inclusions over a range of temperature and may reflect the presence of µm-scale and immiscible Fe-rich and Si-rich components in different proportions. Simulating results for the motion of Si-rich droplets within a crystal mush indicate that Si-rich droplets would be separated from Fe-rich melt and migrate upward due to density differences in the interstitial liquid when the magma unmixed. Migration of the Si-rich, immiscible liquid component from the interstitial liquid caused the remaining Fe-rich melt in the lower part to react with plagioclase primocrysts (An_59–60), as evidenced by fine-grained lamellar intergrowth of An-rich plagioclase (An_79–84)?+?clinopyroxene in the oxide gabbro of the lower zone. Therefore, magma unmixing within a crystal mush, combined with gravitationally driven loss of the Si-rich component, resulted in the formation of Fe-rich, melagabbro and major Fe–Ti oxide ores in the lower part and Si-rich, leucogabbro in the upper part of the intrusion.     

Carbon costs and benefits of France’s biomass energy production targets

Background

Concern about climate change has motivated France to reduce its reliance on fossil fuel by setting targets for increased biomass-based renewable energy production. This study quantifies the carbon costs and benefits for the French forestry sector in meeting these targets. A forest growth and harvest simulator was developed for French forests using recent forest inventory data, and the wood-use chain was reconstructed from national wood product statistics. We then projected wood production, bioenergy production, and carbon balance for three realistic intensification scenarios and a business-as-usual scenario. These intensification scenarios targeted either overstocked, harvest-delayed or currently actively managed stands.

Results

All three intensification strategies produced 11.6–12.4 million tonnes of oil equivalent per year of wood-based energy by 2026, which corresponds to the target assigned to French wood-energy to meet the EU 2020 renewable energy target. Sustaining this level past 2026 will be challenging, let alone further increasing it. Although energy production targets can be reached, the management intensification required will degrade the near-term carbon balance of the forestry sector, compared to continuing present-day management. Even for the best-performing intensification strategy, i.e., reducing the harvest diameter of actively managed stands, the carbon benefits would only become apparent after 2040. The carbon balance of a strategy putting abandoned forests back into production would only break even by 2055; the carbon balance from increasing thinning in managed but untended stands would not break even within the studied time periods, i.e. 2015–2045 and 2046–2100. Owing to the temporal dynamics in the components of the carbon balance, i.e., the biomass stock in the forest, the carbon stock in wood products, and substitution benefits, the merit order of the examined strategies varies over time.

Conclusions

No single solution was found to improve the carbon balance of the forestry sector by 2040 in a way that also met energy targets. We therefore searched for the intensification scenario that produces energy at the lowest carbon cost. Reducing rotation time of actively managed stands is slightly more efficient than targeting harvest-delayed stands, but in both cases, each unit of energy produced has a carbon cost that only turns into a benefit between 2060 and 2080.

     

Revisit the calibration errors on experimental slant total electron content (TEC) determined with GPS

The calibration errors on experimental slant total electron content (TEC) determined with global positioning system (GPS) observations is revisited. Instead of the analysis of the calibration errors on the carrier phase leveled to code ionospheric observable, we focus on the accuracy analysis of the undifferenced ambiguity-fixed carrier phase ionospheric observable determined from a global distribution of permanent receivers. The results achieved are: (1) using data from an entire month within the last solar cycle maximum, the undifferenced ambiguity-fixed carrier phase ionospheric observable is found to be over one order of magnitude more accurate than the carrier phase leveled to code ionospheric observable and the raw code ionospheric observable. The observation error of the undifferenced ambiguity-fixed carrier phase ionospheric observable ranges from 0.05 to 0.11 total electron content unit (TECU) while that of the carrier phase leveled to code and the raw code ionospheric observable is from 0.65 to 1.65 and 3.14 to 7.48 TECU, respectively. (2) The time-varying receiver differential code bias (DCB), which presents clear day boundary discontinuity and intra-day variability pattern, contributes the most part of the observation error. This contribution is assessed by the short-term stability of the between-receiver DCB, which ranges from 0.06 to 0.17 TECU in a single day. (3) The remaining part of the observation errors presents a sidereal time cycle pattern, indicating the effects of the multipath. Further, the magnitude of the remaining part implies that the code multipath effects are much reduced. (4) The intra-day variation of the between-receiver DCB of the collocated stations suggests that estimating DCBs as a daily constant can have a mis-modeling error of at least several tenths of 1 TECU.