Computation of a combined spherical-elastic and viscous-half-space earth model for ice sheet simulation

This report starts by describing the continuum model used by Lingle & Clark (1985) to approximate the deformation of the earth under changing ice sheet and ocean loads. That source considers a single ice stream, but we apply their underlying model to continent-scale ice sheet simulation. Their model combines Farrell's (1972) elastic spherical earth with a viscous half-space overlain by an elastic plate lithosphere. The latter half-space model is derivable from calculations by Cathles (1975). For the elastic spherical earth we use Farrell's tabulated Green's function, as do Lingle & Clark. For the half-space model, however, we propose and implement a significantly faster numerical strategy, a spectral collocation method (Trefethen 2000) based directly on the Fast Fourier Transform. To verify this method we compare to an integral formula for a disc load. To compare earth models we build an accumulation history from a growing similarity solution from (Bueler, et al.~2005) and and simulate the coupled (ice flow)-(earth deformation) system. In the case of simple isostasy the exact solution to this system is known. We demonstrate that the magnitudes of numerical errors made in approximating the ice-earth system are significantly smaller than pairwise differences between several earth models, namely, simple isostasy, the current standard model used in ice sheet simulation (Greve 2001, Hagdorn 2003, Zweck & Huybrechts 2005), and the Lingle & Clark model. Therefore further efforts to validate different earth models used in ice sheet simulations are, not surprisingly, worthwhile.Comment: 36 pages, 16 figures, 3 Matlab program

Clumped isotope evidence for episodic, rapid flow of fluids in a mineralized fault system in the Peak District, UK

We have used clumped isotope thermometry to study a fault-hosted hydrothermal calcite vein associated with the Mississippi Valley Type (MVT) mineralization on the Derbyshire Platform in the southern Pennines, UK. This is the first published dataset obtained using a new mass spectrometer, MIRA, optimized for clumped isotope analysis and an associated clumped isotope-temperature calibration. We analysed multiple generations of vein growth at high spatial resolution in two transects across the vein. The vein grew episodically at temperatures between 40 and 100°C. We interpret each episode of growth as being associated with an increasing flux of formation waters from deep sedimentary basins next to the mineralized platform and an accompanying increase in the precipitation temperatures. Heat is conserved in the fluid as it ascends along the fault surface and, thus, flow must have been fast and restricted to short-lived pulses. The flux could have been driven by high pore pressures associated with rapid sedimentation, hydrocarbon generation and diagenesis in the basinal facies of the Visean Bowland-Hodder group. Natural hydraulic fracturing of shale units and failure of capillary seals, possibly triggered by uplift, allowed the release of fluids into aquifers within the sediment pile. The transmission of high pore fluid pressures from the shales to the fault zone, aided by the compressibility of the gas phase in two-phase pore fluids, may have resulted in fault rupture, accompanied by enhanced fracture permeability and rapid fluid flow. Vein growth ceased as pore pressures dissipated. Such behaviour is closely related to a seismic valve type model for mineralization

EIB Working Paper 2020/06 - Digital technologies and firm performance

As the productivity of the European economy shows signs of slowing down, many hopes are pinned on digital technologies to reverse this trend. This study uses data from the EIBIS 2019 survey to examine whether the adoption of different digital technologies (such as advanced robotics, 3D printing, or Internet of Things) by firms in the EU have different impacts on productivity. It also examines whether these different technologies have different implications for employment growth, and whether there are complementarities between technologies when it comes to firm performance

Quantification of CO2 generation in sedimentary basins through carbonate/clays reactions with uncertain thermodynamic parameters

We develop a methodological framework and mathematical formulation which yields estimates of the uncertainty associated with the amounts of CO2generated by Carbonate-Clays Reactions (CCR) in large-scale subsurface systems to assist characterization of the main features of this geochemical process. Our approach couples a one-dimensional compaction model, providing the dynamics of the evolution of porosity, temperature and pressure along the vertical direction, with a chemical model able to quantify the partial pressure of CO2resulting from minerals and pore water interaction. The modeling framework we propose allows (i) estimating the depth at which the source of gases is located and (ii) quantifying the amount of CO2generated, based on the mineralogy of the sediments involved in the basin formation process. A distinctive objective of the study is the quantification of the way the uncertainty affecting chemical equilibrium constants propagates to model outputs, i.e., the flux of CO2. These parameters are considered as key sources of uncertainty in our modeling approach because temperature and pressure distributions associated with deep burial depths typically fall outside the range of validity of commonly employed geochemical databases and typically used geochemical software. We also analyze the impact of the relative abundancy of primary phases in the sediments on the activation of CCR processes. As a test bed, we consider a computational study where pressure and temperature conditions are representative of those observed in real sedimentary formation. Our results are conducive to the probabilistic assessment of (i) the characteristic pressure and temperature at which CCR leads to generation of CO2in sedimentary systems, (ii) the order of magnitude of the CO2generation rate that can be associated with CCR processes

Shale gas production: potential versus actual greenhouse gas emissions

Estimates of greenhouse gas (GHG) emissions from shale gas production and use are controversial. Here we assess the level of GHG emissions from shale gas well hydraulic fracturing operations in the United States during 2010. Data from each of the approximately 4000 horizontal shale gas wells brought online that year are used to show that about 900 Gg CH[subscript 4] of potential fugitive emissions were generated by these operations, or 228 Mg CH[subscript 4] per well—a figure inappropriately used in analyses of the GHG impact of shale gas. In fact, along with simply venting gas produced during the completion of shale gas wells, two additional techniques are widely used to handle these potential emissions: gas flaring and reduced emission 'green' completions. The use of flaring and reduced emission completions reduce the levels of actual fugitive emissions from shale well completion operations to about 216 Gg CH[subscript 4], or 50 Mg CH[subscript 4] per well, a release substantially lower than several widely quoted estimates. Although fugitive emissions from the overall natural gas sector are a proper concern, it is incorrect to suggest that shale gas-related hydraulic fracturing has substantially altered the overall GHG intensity of natural gas production

Submarine gas seepage in a mixed contractional and shear deformation regime: Cases from the Hikurangi oblique-subduction margin

Gas seepage from marine sediments has implications for understanding feedbacks between the global carbon reservoir, seabed ecology and climate change. Although the relationship between hydrates, gas chimneys and seafloor seepage is well established, the nature of fluid sources and plumbing mechanisms controlling fluid escape into the hydrate zone and up to the seafloor remain one of the least understood components of fluid migration systems. In this study we present the analysis of new three-dimensional high-resolution seismic data acquired to investigate fluid migration systems sustaining active seafloor seepage at Omakere Ridge, on the Hikurangi subduction margin, New Zealand. The analysis reveals at high resolution, complex overprinting fault structures (i.e. protothrusts, normal faults from flexural extension, and shallow (<1 km) arrays of oblique shear structures) implicated in fluid migration within the gas hydrate stability zone in an area of 2x7 km. In addition to fluid migration systems sustaining seafloor seepage on both sides of a central thrust fault, the data show seismic evidence for sub-seafloor gas-rich fluid accumulation associated with proto-thrusts and extensional faults. In these latter systems fluid pressure dissipation through time has been favored, hindering the development of gas chimneys. We discuss the elements of the distinct fluid migration systems and the influence that a complex partitioning of stress may have on the evolution of fluid flow systems in active subduction margins

Mantle plume capture, anchoring, and outflow during Galápagos plume-ridge interaction

Compositions of basalts erupted between the main zone of Galápagos plume upwelling and adjacent Galápagos Spreading Center (GSC) provide important constraints on dynamic processes involved in transfer of deep-mantle-sourced material to mid-ocean ridges. We examine recent basalts from central and northeast Galápagos including some that have less radiogenic Sr, Nd, and Pb isotopic compositions than plume-influenced basalts (E-MORB) from the nearby ridge. We show that the location of E-MORB, greatest crustal thickness, and elevated topography on the GSC correlates with a confined zone of low-velocity, high-temperature mantle connecting the plume stem and ridge at depths of ∼100 km. At this site on the ridge, plume-driven upwelling involving deep melting of partially dehydrated, recycled ancient oceanic crust, plus plate-limited shallow melting of anhydrous peridotite, generate E-MORB and larger amounts of melt than elsewhere on the GSC. The first-order control on plume stem to ridge flow is rheological rather than gravitational, and strongly influenced by flow regimes initiated when the plume was on axis (>5 Ma). During subsequent northeast ridge migration material upwelling in the plume stem appears to have remained “anchored” to a contact point on the GSC. This deep, confined NE plume stem-to-ridge flow occurs via a network of melt channels, embedded within the normal spreading and advection of plume material beneath the Nazca plate, and coincides with locations of historic volcanism. Our observations require a more dynamically complex model than proposed by most studies, which rely on radial solid-state outflow of heterogeneous plume material to the ridge.We thank Galápagos National Park authorities and CDRS for permitting fieldwork in Galápagos. D. Villagomez and D. Toomey generously shared their extensive seismic data set for Galápagos, and D. McKenzie kindly provided help with temperature calculations. End-member compositions of Galápagos mantle reservoirs in Figure 4 were estimated from principal component analysis; data related to these calculations are available in the supporting information. We are grateful to Kaj Hoernle and two anonymous reviewers for their constructive comments on an earlier version of this manuscript. The research was funded by the University of Cambridge, Geological Society of London, NERC (RG57434), and NSF (EAR 0838461, EAR 0944229, and EAR-11452711).This is the final published version of the article. It first appeared at http://dx.doi.org/10.1002/2015GC00572