An automated strategy for calculation of phase diagram sections and retrieval of rock properties as a function of physical conditions

We formulate an algorithm for the calculation of stable phase relations of a system with constrained bulk composition as a function of its environmental variables. The basis of this algorithm is the approximate representation of the free energy composition surfaces of solution phases by inscribed polyhedra. This representation leads to discretization of high variance phase fields into a continuous mesh of smaller polygonal fields within which the composition and physical properties of the phases are uniquely determined. The resulting phase diagram sections are useful for understanding the phase relations of complex metamorphic systems and for applications in which it is necessary to establish the variations in rock properties such as density, seismic velocities and volatile‐content through a metamorphic cycle. The algorithm has been implemented within a computer program that is general with respect to both the choice of variables and the number of components and phases possible in a system, and is independent of the structure of the equations of state used to describe the phases of the system.     

The OSIRIS‐REx target asteroid (101955) Bennu: Constraints on its physical,geological, and dynamical nature from astronomical observations

We review the results of an extensive campaign to determine the physical, geological, and dynamical properties of asteroid (101955) Bennu. This investigation provides information on the orbit, shape, mass, rotation state, radar response, photometric, spectroscopic, thermal, regolith, and environmental properties of Bennu. We combine these data with cosmochemical and dynamical models to develop a hypothetical timeline for Bennu's formation and evolution. We infer that Bennu is an ancient object that has witnessed over 4.5 Gyr of solar system history. Its chemistry and mineralogy were established within the first 10 Myr of the solar system. It likely originated as a discrete asteroid in the inner Main Belt approximately 0.7–2 Gyr ago as a fragment from the catastrophic disruption of a large (approximately 100‐km), carbonaceous asteroid. It was delivered to near‐Earth space via a combination of Yarkovsky‐induced drift and interaction with giant‐planet resonances. During its journey, YORP processes and planetary close encounters modified Bennu's spin state, potentially reshaping and resurfacing the asteroid. We also review work on Bennu's future dynamical evolution and constrain its ultimate fate. It is one of the most Potentially Hazardous Asteroids with an approximately 1‐in‐2700 chance of impacting the Earth in the late 22nd century. It will most likely end its dynamical life by falling into the Sun. The highest probability for a planetary impact is with Venus, followed by the Earth. There is a chance that Bennu will be ejected from the inner solar system after a close encounter with Jupiter. OSIRIS‐REx will return samples from the surface of this intriguing asteroid in September 2023.     

Grain‐scale pressure variations and chemical equilibrium in high‐grade metamorphic rocks

In the classical view of metamorphic microstructures, fast viscous relaxation (and so constant pressure) is assumed, with diffusion being the limiting factor in equilibration. This contribution is focused on the only other possible scenario – fast diffusion and slow viscous relaxation – and brings an alternative interpretation of microstructures typical of high‐grade metamorphic rocks. In contrast to the pressure vessel mechanical model applied to pressure variation associated with coesite inclusions in various host minerals, a multi‐anvil mechanical model is proposed in which strong single crystals and weak grain boundaries can maintain pressure variation at geological time‐scales in a polycrystalline material. In such a mechanical context, exsolution lamellae in feldspar are used to show that feldspar can sustain large differential stresses (>10 kbar) at geological time‐scales. Furthermore, it is argued that the existence of grain‐scale pressure gradients combined with diffusional equilibrium may explain chemical zoning preserved in reaction rims. Assuming zero net flux across the microstructure, an equilibrium thermodynamic method is introduced for inferring pressure variation corresponding to the chemical zoning. This new barometric method is applied to plagioclase rims around kyanite in felsic granulite (Bohemian Massif, Czech Republic), yielding a grain‐scale pressure variation of 8 kbar. In this approach, kinetic factors are not invoked to account for mineral composition zoning preserved in rocks metamorphosed at high grade.     

DHI evaluation by combining rock physics simulation and statistical techniques for fluid identification of Cambrian-to-Cretaceous clastic reservoirs in Pakistan

The use of seismic direct hydrocarbon indicators is very common in exploration and reservoir development to minimise exploration risk and to optimise the location of production wells. DHIs can be enhanced using AVO methods to calculate seismic attributes that approximate relative elastic properties. In this study, we analyse the sensitivity to pore fluid changes of a range of elastic properties by combining rock physics studies and statistical techniques and determine which provide the best basis for DHIs. Gassmann fluid substitution is applied to the well log data and various elastic properties are evaluated by measuring the degree of separation that they achieve between gas sands and wet sands. The method has been applied successfully to well log data from proven reservoirs in three different siliciclastic environments of Cambrian, Jurassic, and Cretaceous ages. We have quantified the sensitivity of various elastic properties such as acoustic and extended elastic (EEI) impedances, elastic moduli (K _sat and K _sat–μ), lambda–mu–rho method (λρ and μρ), P-to-S-wave velocity ratio (V _P/V _S), and Poisson’s ratio (σ) at fully gas/water saturation scenarios. The results are strongly dependent on the local geological settings and our modeling demonstrates that for Cambrian and Cretaceous reservoirs, K _sat–μ, EEI, V _P/V _S, and σ are more sensitive to pore fluids (gas/water). For the Jurassic reservoir, the sensitivity of all elastic and seismic properties to pore fluid reduces due to high overburden pressure and the resultant low porosity. Fluid indicators are evaluated using two metrics: a fluid indicator coefficient based on a Gaussian model and an overlap coefficient which makes no assumptions about a distribution model. This study will provide a potential way to identify gas sand zones in future exploration.     

Global extent and distribution of artificial, residential waterways in estuaries

Artificial residential waterways are now widespread in the world’s estuaries. We used the global mapping tool, Google Earth, to determine that there are nearly 4000 linear km of artificial waterways globally, covering an area of 270 km². Residential waterways constructed as open, flow-through canal estates are at their greatest extent in North America (77% of global linear extent), where systems are typically longer and narrower, with more openings and dead-ends than systems elsewhere. The remaining canal estates are spread across all other continents except Antarctica: Asia (7%), Europe (7%), Oceania (7%), South America (0.9%), and Africa (0.6%). A relatively recent design change from open canals to artificial estuarine lakes with tidal barriers has occurred on all continents except Africa, most extensively in Australia (14 km² area, 57 independent systems). The extremely large expansion in artificial residential waterways aimed at increasing opportunities for waterfront living by humans has also modified and expanded estuarine habitat available to aquatic biota. Research can best underpin planning and management of these types of waterways by focussing on their value as habitat and their provision of other goods and services.