J. van Hunen

Small-scale convection at the edge of the Colorado Plateau: Implications for topography, magmatism, and evolution of Proterozoic lithosphere

The Colorado Plateau of the southwestern United States is characterized by a bowl-shaped high elevation, late Neogene–Quaternary magmatism at its edge, large gradients in seismic wave velocity across its margins, and relatively low lithospheric seismic wave velocities. We explain these observations by edge-driven convection following rehydration of Colorado Plateau lithosphere. A rapidly emplaced Cenozoic step in lithosphere thickness between the Colorado Plateau and adjacent extended Rio Grande rift and Basin and Range province causes small-scale convection in the asthenosphere. A lithospheric drip below the plateau is removing lithosphere material from the edge that is heated and metasomatized, resulting in magmatism. Edgedriven convection also drives margin uplift, giving the plateau its characteristic bowl shape. The edge-driven convection model shows good consistency with features resolved by seismic tomography.

Small-scale convection during continental rifting: Evidence from the Rio Grande rift

Recent seismic imaging across the Rio Grande rift, western United States, revealed unexpected structures in the underlying mantle. Low seismic wave velocity anomalies below the Rio Grande rift have been interpreted as being partially of melt origin, and high-velocity structures below the western Great Plains have been proposed to be the result of small-scale convection, i.e., cold downwelling lithospheric material with probably a compositional contribution. We perform a dynamic test of these interpretations using a passive rift model for isochemical convection. The models self-consistently produce a rift localized at approximately the right distance from the border to the nearby thicker Great Plains lithosphere. With realistic upper mantle rheologies, small-scale convection forms, aided by the lithospheric step. The resulting thermal anomalies produce seismic low-velocity anomalies below the rift of amplitudes similar to those imaged seismically, requiring the presence of only small amounts of melt. The lateral extent of the observed low velocities below the Rio Grande rift is as in the models, where it is controlled by the spacing between downwelling limbs of the small-scale convection. The fast velocity structure below the western Great Plains can be produced by cold downwelling lithosphere. The thermal rifting models can predict the amplitudes and size of the main seismic anomalies; compositional heterogeneity may contribute to some of the smaller features observed.

Sublithospheric small-scale convection—A mechanism for collision zone magmatism

We studied the effect of increased water content on the dynamics of the lithosphere-asthenosphere boundary in a postsubduction setting. Results from numerical mantle convection models show that the resultant decrease in mantle viscosity and the peridotite solidus produce small-scale convection at the lithosphere-asthenosphere boundary and magmatism that follows the spatially and temporally scattered style and volumes typical for collision magmatism, such as the late Cenozoic volcanism of the Turkish-Iranian Plateau. An inherent feature in small-scale convection is its chaotic nature that can lead to temporally isolated volcanic centers tens of millions of years after initial continental collision, without evident tectonic cause. We also conclude that water input into the upper mantle during and after subduction under the circum-Mediterranean area and the Tibetan Plateau can account for the observed magmatism in these areas. Only fractions (200–600 ppm) of the water input need to be retained after subduction to induce small-scale convection and magmatism on the scale of those observed from the Turkish-Iranian Plateau.

Plate rotation during continental collision and its relationship with the exhumation of UHP metamorphic terranes: Application to the Norwegian Caledonides

Lateral variation and asynchronous onset of collision during the convergence of continents can significantly affect the burial and exhumation of subducted continental crust. Here we use 3-D numerical models for continental collision to discuss how deep burial and exhumation of high and ultrahigh pressure metamorphic (HP/UHP) rocks are enhanced by diachronous collision and the resulting rotation of the colliding plates. Rotation during collision locally favors eduction, the inversion of the subduction, and may explain the discontinuous distribution of ultra-high pressure (UHP) terranes along collision zones. For example, the terminal (Scandian) collision of Baltica and Laurentia, which formed the Scandinavian Caledonides, resulted in the exhumation of only one large HP/UHP terrane, the Western Gneiss Complex (WGC), near the southern end of the collision zone. Rotation of the subducting Baltica plate during collision may provide an explanation for this distribution. We explore this hypothesis by comparing orthogonal and diachronous collision models and conclude that a diachronous collision can transport continental material up to 60 km deeper, and heat material up to 300°C hotter, than an orthogonal collision. Our diachronous collision model predicts that subducted continental margin material returns to the surface only in the region where collision initiated. The diachronous collision model is consistent with petrological and geochonological observations from the WGC and makes predictions for the general evolution of the Scandinavian Caledonides. We propose the collision between Laurentia and Baltica started at the southern end of the collisional zone, and propagated northward. This asymmetric geometry resulted in the counter clockwise rotation of Baltica with respect to Laurentia, consistent with paleomagnetic data from other studies. Our model may have applications to other orogens with regional UHP terranes, such as the Dabie Shan and Papua New Guinea cases, where block rotation during exhumation has also been recorded.

Application of material balance methods to CO2 storage capacity estimation within selected depleted gas reservoirs

Depleted gas reservoirs are potential sites for CO2 storage; therefore, it is important to evaluate their storage capacity. Historically, there have been difficulties in identifying the reservoir drive mechanism of gas reservoirs using traditional P/z plots, having direct impacts for the estimation of the original gas in place (OGIP) and dependent parameters for both theoretical and effective CO2 storage capacity estimation. Cole plots have previously provided an alternative method of characterization, being derived from the gas material balance equation. We use production data to evaluate the reservoir drive mechanism in four depleted gas reservoirs (Hewett Lower Bunter, Hewett Upper Bunter, and North and South Morecambe) on the UK Continental Shelf. Cole plots suggest that the North Morecambe and Hewett Upper Bunter reservoirs experience moderate water drive. Accounting for cumulative water influx into these reservoirs, the OGIP decreases by up to 20% compared with estimates from P/z plots. The revised OGIP values increase recovery factors within these reservoirs; hence, geometrically based theoretical storage capacity estimates for the North Morecambe and Hewett Upper Bunter reservoirs increase by 4 and 30%, respectively. Material balance approaches yield more conservative estimates. Effective storage capacity estimates are between 64 and 86% of theoretical estimates within the depletion drive reservoirs, and are 53 – 79% within the water drive reservoirs. Supplementary material: A more detailed description of the aquifer modelling is available at https://doi.org/10.6084/m9.figshare.c.3803770.v1