Phase equilibrium modelling of the amphibolite to granulite facies transition in metabasic rocks (Ivrea Zone,NW Italy)

The development of thermodynamic models for tonalitic melt and the updated clinopyroxene and amphibole models now allow the use of phase equilibrium modelling to estimate P–T conditions and melt production for anatectic mafic and intermediate rock types at high‐T conditions. The Permian mid‐lower crustal section of the Ivrea Zone preserves a metamorphic field gradient from mid amphibolite facies to granulite facies, and thus records the onset of partial melting in metabasic rocks. Interlayered metabasic and metapelitic rocks allows the direct comparison of P–T estimates and partial melting between both rock types with the same metamorphic evolution. Pseudosections for metabasic compositions calculated in the Na₂O–CaO–K₂O–FeO–MgO–Al₂O₃–SiO₂–H₂O–TiO₂–O (NCKFMASHTO) system are presented and compared with those of metapelitic rocks calculated with consistent end‐member data and a–x models. The results presented in this study show that P–T conditions obtained by phase equilibria modelling of both metabasic and metapelitic rocks give consistent results within uncertainties, allowing integration of results obtained for both rock types. In combination, the calculations for both metabasic and metapelitic rocks allows an updated and more precisely constrained metamorphic field gradient for Val Strona di Omegna to be defined. The new field gradient has a slightly lower dP/dT which is in better agreement with the onset of crustal thinning of the Adriatic margin during the Permian inferred in recent studies.     

New constraints on granulite facies metamorphism and melt production in the Lewisian Complex,northwest Scotland

In this study, we investigate the metamorphic history of the Assynt and Gruinard blocks of the Archean Lewisian Complex, northwest Scotland, which are considered by some to represent discrete crustal terranes. For samples of mafic and intermediate rocks, phase diagrams were constructed in the Na₂O–CaO–K₂O–FeO–MgO–Al₂O₃–SiO₂–H₂O–TiO₂–O₂ (NCKFMASHTO) system using whole‐rock compositions. Our results indicate that all samples equilibrated at similar peak metamorphic conditions of ~8–10 kbar and ~900–1,000°C, consistent with field evidence for in situ partial melting and the classic interpretation of the central region of the Lewisian Complex as representing a single crustal block. Melt‐reintegration modelling was employed in order to estimate probable protolith compositions. Phase equilibria calculated for these modelled undepleted precursors match well with those determined for a subsolidus amphibolite from Gairloch in the southern region of the Lewisian Complex. Both subsolidus lithologies exhibit similar phase relations and potential melt fertility, with both expected to produce orthopyroxene‐bearing hornblende granulites, with or without garnet, at the conditions inferred for the Badcallian metamorphic peak. For fully hydrated protoliths, prograde melting is predicted to first occur at ~620°C and ~9.5 kbar, with up to 45% partial melt predicted to form at peak conditions in a closed‐system environment. Partial melts calculated for both compositions between 610 and 1,050°C are mostly trondhjemitic. Although the melt‐reintegrated granulite is predicted to produce more potassic (granitic) melts at ~700–900°C, the modelled melts are consistent with the measured compositions of felsic sheets from the central region Lewisian Complex.     

Recent advances in offshore geotechnics for deep water oil and gas developments

The paper presents an overview of recent developments in geotechnical analysis and design associated with oil and gas developments in deep water. Typically the seabed in deep water comprises soft, lightly overconsolidated, fine grained sediments, which must support a variety of infrastructure placed on the seabed or anchored to it. A particular challenge is often the mobility of the infrastructure either during installation or during operation, and the consequent disturbance and healing of the seabed soil, leading to changes in seabed topography and strength. Novel aspects of geotechnical engineering for offshore facilities in these conditions are reviewed, including: new equipment and techniques to characterise the seabed; yield function approaches to evaluate the capacity of shallow skirted foundations; novel anchoring systems for moored floating facilities; pipeline and steel catenary riser interaction with the seabed; and submarine slides and their impact on infrastructure. Example results from sophisticated physical and numerical modelling are presented.     

Numerical studies of methane production from Class 1 gas hydrate accumulations enhanced with carbon dioxide injection

Class 1 gas hydrate accumulations are characterized by a permeable hydrate-bearing interval overlying a permeable interval with mobile gas, sandwiched between two impermeable intervals. Depressurization-induced dissociation is currently the favored technology for producing gas from Class 1 gas hydrate accumulations. The depressurization production technology requires heat transfer from the surrounding environment to sustain dissociation as the temperature drops toward the hydrate equilibrium point and leaves the reservoir void of gas hydrate. Production of gas hydrate accumulations by exchanging carbon dioxide with methane in the clathrate structure has been demonstrated in laboratory experiments and proposed as a field-scale technology. The carbon dioxide exchange technology has the potential for yielding higher production rates and mechanically stabilizing the reservoir by maintaining hydrate saturations. We used numerical simulation to investigate the advantages and disadvantages of using carbon dioxide injection to enhance the production of methane from Class 1 gas hydrate accumulations. Numerical simulations in this study were primarily concerned with the mechanisms and approaches of carbon dioxide injection to investigate whether methane production could be enhanced through this approach. To avoid excessive simulation execution times, a five-spot well pattern with a 500-m well spacing was approximated using a two-dimensional domain having well boundaries on the vertical sides and impermeable boundaries on the horizontal sides. Impermeable over- and under burden were included to account for heat transfer into the production interval. Simulation results indicate that low injection pressures can be used to reduce secondary hydrate formation and that direct contact of injected carbon dioxide with the methane hydrate present in the formation is limited due to bypass through the higher permeability gas zone.     

Multiseason evaluation of the MM5, COAMPS and WRF over southeast United States

Three models, MM5, COAMPS, and WRF, have been applied for the warm season in 2003 and the cool season in 2003?C2004 to evaluate their performances. All models run over the same domain area covering the north Gulf Mexico and southeastern United States (US) region with the same spatial resolution of 27?km. It was found that the temporal variations of the mean error distribution and strength at 24 and 36?h were rather weak for surface temperature, sea level pressure, and surface wind speed for all models. A warm bias in surface temperature forecasts dominated over land during the warm season, whereas a cool bias existed during the cool season. The MM5 and WRF produced negative biases of sea level pressure during the warm season and positive biases during the cool season while the COAMPS yielded a similar distribution of sea level pressure biases during both seasons. During both seasons, similar surface wind speed biases produced by each model included a high wind speed forecast over most areas by MM5 while the COAMPS and WRF yielded weak surface winds over the western Plains and stronger surface winds over the eastern Plains. Root-mean-squared errors revealed that the forecast of surface temperature, sea level pressure, and surface wind speed were degraded with the increase of forecast time. For rainfall evaluation, it was found that the MM5 underpredicted seasonal precipitation while the COAMPS and WRF overpredicted. The bias scores revealed that the MM5 yielded an underprediction of the coverage of precipitation areas, especially for heavier rainfall events. The MM5 presented the lower threat score at lighter rainfall events compared to the COAMPS and WRF. For moderate and heavier thresholds, all models lacked forecast accuracy. The WRF accuracy in predicting precipitation was heavily dependent upon the performance of the selected cumulus parameterization scheme. Use of the Grell?CDevenyi and Bette?CMiller?CJanjic schemes helps suppress precipitation overprediction.     

Well tie for broadband seismic data

The seismic industry is increasingly acquiring broadband data in order to reap the benefits of extra low‐ and high‐frequency contents. At the low end, as the sharp low‐cut decay gets closer to zero frequency, it becomes harder for a well tie to estimate the low‐frequency response correctly. The fundamental difficulty is that well logs are too short to allow accurate estimation of the long‐period content of the data. Three distinctive techniques, namely parametric constant phase, frequency‐domain least squares with multi‐tapering, and Bayesian time domain with broadband priors, are introduced in this paper to provide a robust solution to the wavelet estimation problem for broadband seismic data. Each of these techniques has a different mathematical foundation that would enable one to explore a wide range of solutions that could be used on a case‐by‐case basis depending on the problem at hand. A case study from the North West Shelf Australia is used to analyse the performance of the proposed techniques. Cross‐validation is proposed as a robust quality control measure for evaluating well‐tie applications. It is observed that when the seismic data are carefully processed, then the constant phase approach would likely offer a good solution. The frequency‐domain method does not assume a constant phase. This flexibility makes it prone to over‐fitting when the phase is approximately constant. Broadband priors for the time‐domain least‐squares method are found to perform well in defining low‐frequency side lobes to the wavelet.     

High‐grade metamorphism and partial melting in Archean composite grey gneiss complexes

Much of the exposed Archean crust is composed of composite gneiss which includes a large proportion of intermediate to tonalitic material. These gneiss terranes were typically metamorphosed to amphibolite to granulite facies conditions, with evidence for substantial partial melting at higher grade. Recently published activity–composition (a?x) models for partial melting of metabasic to intermediate compositions allows calculation of the stable metamorphic minerals, melt production and melt composition in such rocks for the first time. Calculated P?T pseudosections are presented for six bulk rock compositions taken from the literature, comprising two metabasic compositions, two intermediate/dioritic compositions and two tonalitic compositions. This range of bulk compositions captures much of the diversity of rock types found in Archean banded gneiss terranes, enabling us to present an overview of metamorphism and partial melting in such terranes. If such rocks are fluid saturated at the solidus, they first begin to melt in the upper amphibolite facies. However, at such conditions, very little (< 5%) melt is produced and this melt is granitic in composition for all rocks. The production of greater proportions of melt requires temperatures ～800–850 °C and is associated with the first appearance of orthopyroxene at pressures below 8–9 kbar or with the appearance and growth of garnet at higher pressures. The temperature at which orthopyroxene appears varies little with composition providing a robust estimate of the amphibolite–granulite facies boundary. Across this boundary, melt production is coincident with the breakdown of hornblende and/or biotite. Melts produced at granulite facies range from tonalite–trondhjemite–granodiorite for the metabasic protoliths, granodiorite to granite for the intermediate protoliths and granite for the tonalitic protoliths. Under fluid‐absent conditions the melt fertility of the different protoliths is largely controlled by the relative proportions of hornblende and quartz at high grade, with the intermediate compositions being the most fertile. The least fertile rocks are the most leucocratic tonalites due to their relatively small proportions of hydrous mafic phases such as hornblende or biotite. In the metabasic rocks, melt production becomes limited by the complete consumption of quartz to higher temperatures. The use of phase equilibrium forward‐modelling provides a thermodynamic framework for understanding melt production, melt loss and intracrustal differentiation during the Archean.     

Anisotropic damage of rock joints during cyclic loading: constitutive framework and numerical integration

This work describes a constitutive framework for modeling the behavior of rough joints under cyclic loading. Particular attention is paid to the intrinsic links between dilatancy, surface degradation, and mobilized shear strength. The framework also accounts for the important effect of shear‐induced anisotropy. The resulting approach is fully three‐dimensional and is not restricted to plane‐displacement kinematics. Both the governing formulation and an algorithm for implicit numerical integration are presented. While the proposed methods are general, we also postulate a specific model that is compared with experimental data. It employs relatively few free parameters but shows good agreement with laboratory tests. Copyright © 2013 John Wiley & Sons, Ltd.     

New mineral activity–composition relations for thermodynamic calculations in metapelitic systems

New activity–composition (a–x) relations for minerals commonly occurring in metapelites are presented for use with the internally consistent thermodynamic dataset of Holland & Powell ( 2011 , Journal of Metamorphic Geology, 29 , 333–383). The a–x relations include a broader consideration of Fe₂O₃ in minerals, changes to the formalism of several phases and order–disorder in all ferromagnesian minerals where Fe–Mg mixing occurs on multiple sites. The a–x relations for chlorite, biotite, garnet, chloritoid, staurolite, cordierite, orthopyroxene, muscovite, paragonite and margarite have been substantially reparameterized using the approach outlined in the companion paper in this issue. For the first time, the entire set of a–x relations for the common ferromagnesian minerals in metapelitic rocks is parameterized simultaneously, with attention paid to ensuring that they can be used together to calculate phase diagrams of geologically appropriate topology. The a–x relations developed are for use in the Na₂O–CaO–K₂O–FeO–MgO–Al₂O₃–SiO₂–H₂O–TiO₂–O₂ (NCKFMASHTO) system for both subsolidus and suprasolidus conditions. Petrogenetic grids in KFMASH and KFMASHTO are similar in topology to those produced with earlier end‐member datasets and a–x relations, but with some notable differences. In particular, in subsolidus equilibria, the FeO/(FeO + MgO) of garnet is now greater than in coexisting staurolite, bringing a number of key staurolite‐bearing equilibria into better agreement with inferences from field and petrographic observations. Furthermore, the addition of Fe³⁺ and Ti to a number of silicate phases allows more plausible equilibria to be calculated in relevant systems. Pseudosections calculated with the new a–x relations are also topologically similar to equivalent diagrams using earlier a–x relations, although with many low variance fields shifting in P–T space to somewhat lower pressure conditions.