Inexact methods for sequential fully implicit (SFI) reservoir simulation

Computational Geosciences - The sequential fully implicit (SFI) scheme was introduced (Jenny et al. J. Comput. Phys. 217(2), 627–641 2006) for solving coupled flow and transport problems....     

The Fe–C–O–H–N system at 6.3–7.8 GPa and 1200–1400 °C: implications for deep carbon and nitrogen cycles

Interactions in a Fe–C–O–H–N system that controls the mobility of siderophile nitrogen and carbon in the Fe⁰-saturated upper mantle are investigated in experiments at 6.3–7.8 GPa and 1200–1400 °C. The results show that the γ-Fe and metal melt phases equilibrated with the fluid in a system unsaturated with carbon and nitrogen are stable at 1300 °C. The interactions of Fe₃C with an N-rich fluid in a graphite-saturated system produce the ε-Fe₃N phase (space group P6₃/mmc or P6₃22) at subsolidus conditions of 1200–1300 °C, while N-rich melts form at 1400 °C. At IW- and MMO-buffered hydrogen fugacity (fH₂), fluids vary from NH₃- to H₂O-rich compositions (NH₃/N₂?>?1 in all cases) with relatively high contents of alkanes. The fluid derived from N-poor samples contains less H₂O and more carbon which mainly reside in oxygenated hydrocarbons, i.e., alcohols and esters at MMO-buffered fH₂ and carboxylic acids at unbuffered fH₂ conditions. In unbuffered conditions, N₂ is the principal nitrogen host (NH₃/N₂?≤?0.1) in the fluid equilibrated with the metal phase. Relatively C- and N-rich fluids in equilibrium with the metal phase (γ-Fe, melt, or Fe₃N) are stable at the upper mantle pressures and temperatures. According to our estimates, the metal/fluid partition coefficient of nitrogen is higher than that of carbon. Thus, nitrogen has a greater affinity for iron than carbon. The general inference is that reduced fluids can successfully transport volatiles from the metal-saturated mantle to metal-free shallow mantle domains. However, nitrogen has a higher affinity for iron and selectively accumulates in the metal phase, while highly mobile carbon resides in the fluid phase. This may be a controlling mechanism of the deep carbon and nitrogen cycles.     

Changes in the Probability of Heavy Precipitation: Important Indicators of Climatic Change

A simple statistical model of daily precipitation based on the gamma distribution is applied to summer (JJA in Northern Hemisphere, DJF in Southern Hemisphere) data from eight countries: Canada, the United States, Mexico, the former Soviet Union, China, Australia, Norway, and Poland. These constitute more than 40% of the global land mass, and more than 80% of the extratropical land area. It is shown that the shape parameter of this distribution remains relatively stable, while the scale parameter is most variable spatially and temporally. This implies that the changes in mean monthly precipitation totals tend to have the most influence on the heavy precipitation rates in these countries. Observations show that in each country under consideration (except China), mean summer precipitation has increased by at least 5% in the past century. In the USA, Norway, and Australia the frequency of summer precipitation events has also increased, but there is little evidence of such increases in any of the countries considered during the past fifty years. A scenario is considered, whereby mean summer precipitation increases by 5% with no change in the number of days with precipitation or the shape parameter. When applied in the statistical model, the probability of daily precipitation exceeding 25.4 mm (1 inch) in northern countries (Canada, Norway, Russia, and Poland) or 50.8 mm (2 inches) in mid-latitude countries (the USA, Mexico, China, and Australia) increases by about 20% (nearly four times the increase in mean). The contribution of heavy rains (above these thresholds) to the total 5% increase of precipitation is disproportionally high (up to 50%), while heavy rain usually constitutes a significantly smaller fraction of the precipitation events and totals in extratropical regions (but up to 40% in the tropics, e.g., in southern Mexico). Scenarios with moderate changes in the number of days with precipitation coupled with changes in the scale parameter were also investigated and found to produce smaller increases in heavy rainfall but still support the above conclusions. These scenarios give changes in heavy rainfall which are comparable to those observed and are consistent with the greenhouse-gas-induced increases in heavy precipitation simulated by some climate models for the next century. In regions with adequate data coverage such as the eastern two-thirds of contiguous United States, Norway, eastern Australia, and the European part of the former USSR, the statistical model helps to explain the disproportionate high changes in heavy precipitation which have been observed.     

Far-zone contributions to the gravity field quantities by means of Molodensky’s truncation coefficients

To reduce the numerical complexity of inverse solutions to large systems of discretised integral equations in gravimetric geoid/quasigeoid modelling, the surface domain of Green’s integrals is subdivided into the near-zone and far-zone integration sub-domains. The inversion is performed for the near zone using regional detailed gravity data. The farzone contributions to the gravity field quantities are estimated from an available global geopotential model using techniques for a spherical harmonic analysis of the gravity field. For computing the far-zone contributions by means of Green’s integrals, truncation coefficients are applied. Different forms of truncation coefficients have been derived depending on a type of integrals in solving various geodetic boundary-value problems. In this study, we utilise Molodensky’s truncation coefficients to Green’s integrals for computing the far-zone contributions to the disturbing potential, the gravity disturbance, and the gravity anomaly. We also demonstrate that Molodensky’s truncation coefficients can be uniformly applied to all types of Green’s integrals used in solving the boundaryvalue problems. The numerical example of the far-zone contributions to the gravity field quantities is given over the area of study which comprises the Canadian Rocky Mountains. The coefficients of a global geopotential model and a detailed digital terrain model are used as input data.     

On a physical mechanism for extra mixing in globular cluster red giants

For the first time we propose a real physical mechanism for 'extra mixing' in red giants that can quantitatively interpret all the known star-to-star abundance variations in globular clusters. This is Zahn's mechanism. It considers extra mixing in a radiative zone of a rotating star as a result of the joint operation of meridional circulation and turbulent diffusion. It is shown that the only free parameter, the angular velocity at the base of the convective envelope, can be so adjusted as to fit the observed abundance correlations without leading to a conflict with available data on rotation velocities of blue horizontal branch stars in the same cluster. There are two critical assumptions in our model, that the top of the radiative zone is not in synchronous rotation with the stellar surface but rotates significantly faster and that the criterion for shear instability takes a particular form. These will eventually be tested by three-dimensional hydrodynamical simulations.     

Changes in the radial growth of trees in relation to biogeomorphic processes in an old-growth forest on flysch,Czechia

Tree radial growth is influenced by individual tree abilities, climate, competition, disturbance regimes, as well as biogeomorphic processes – including biomechanical interactions between trees and soil. Trees are actively involved in hillslope dynamics, both responding to and affecting many (bio)geomorphic processes. Using dendrochronology, we studied feedbacks associated with tree–soil–landscape formation, specifically relationships between hillslope processes, biomechanical effects of trees in soils, tree microhabitat conditions and their morphological adaptations, in the flysch zone of the Carpathians. We visually evaluated stem shape, microhabitat conditions and the biomechanical effects of 1663 trees. Cores were taken in four growing directions from 224 individuals of European beech (Fagus sylvatica L.). In a set of 193 cross-dated beeches, average tree-ring widths and tree eccentricities in all directions were calculated and analysed in relation to the biogeomorphic impacts of trees. Some significant drivers of tree radial growth and sources of stem eccentricity were detected. The radial growth of trees on which deadwood was leaning was markedly limited. In contrast, trees with exposed roots expressed the highest growth rates. This clearly suggests that root exposure may not be an effect of ‘exogenous’ soil creep, but may rather result from individually intensifying tree growth due to fine-scale disturbance dynamics. The response of biomechanical tree–soil interactions in tree radial growth weakened with increasing stem diameter, reflecting the stabilizing role of larger trees. The significance of calendar year on radial growth suggests seasonality in the dynamic component of soil creep. Tree eccentricity was observed mainly in the downslope direction, which suggests a relatively complex effect of biomechanics on stem tilting. © 2020 John Wiley & Sons, Ltd.     

The day-to-day effects of strong geomagnetic disturbances on the North Atlantic Oscillation in the winter periods of years 1951–2003

The day-to-day effects of the strong geomagnetic disturbances on geopotential heights (GPH) in the winter lower atmosphere were described in many papers in the beginning of 1970s. These works focused on the North-East Pacific, while the North Atlantic was until now omitted. Our aim is therefore to investigate the possible effect of strong geomagnetic disturbances on the lower atmosphere GPH changes over the winter North Atlantic on the day-to-day time scale, represented by the daily index of the North Atlantic Oscillation (NAO). The investigated intervals are winter periods (December-March) of 1951–2003. The daily NAO average values in 3-day intervals before and after the disturbance onsets are compared. The graphs of NAO differences are complemented by the maps of GPH differences. The NAO response to geomagnetic disturbance, as registered on the day-to-day time scale, also shows a change in its behaviour around the year 1970. This response reaches its highest values in the years 1951–1969, usually 2–5 days following the onset of geomagnetic disturbances. Intensity of the response depends on the disturbance intensity (the largest differences were associated with extremely strong disturbances).

     

Seismic activity and deep conductivity structure of the Eastern Carpathians

We present results of a study of the seismicity and the geoelectric structure of the Eastern Carpathians. After the evaluation of the seismicity, new methods of processing and analyzing seismic data are developed, which allow constructing an averaged horizontal-layered velocity model of the crust in the Carpathian region of Ukraine, tracing the seismic active faults and localizing the seismic events both in horizontal and in vertical direction with a better precision. For the study of the conductivity structure beneath the Eastern Carpathians, the collected magnetovariation and magnetotelluric data are used. The depth of electrical conductivity anomalies are estimated and the resulting quasi-3D model of the conductivity structure beneath the Carpathians is compared with the seismicity in the depth range of 10 ± 2.5 km. The comparison suggests possible geological mechanisms: the seismic events occur mainly in resistive solid rock domains which surround aseismic high-conductivity zones, consisting of at least partially melted material. Aqueous fluids, or a joint effect of several mechanisms, may also play an active role in this distribution.     

Spherical Harmonic Analysis of Gravitational Curvatures and Its Implications for Future Satellite Missions

In this study we assume that a gravitational curvature tensor, i.e. a tensor of third-order directional derivatives of the Earth’s gravitational potential, is observable at satellite altitudes. Such a tensor is composed of ten different components, i.e. gravitational curvatures, which may be combined into vertical–vertical–vertical, vertical–vertical–horizontal, vertical–horizontal–horizontal and horizontal–horizontal-horizontal gravitational curvatures. Firstly, we study spectral properties of the gravitational curvatures. Secondly, we derive new quadrature formulas for the spherical harmonic analysis of the four gravitational curvatures and provide their corresponding analytical error models. Thirdly, requirements for an instrument that would eventually observe gravitational curvatures by differential accelerometry are investigated. The results reveal that measuring third-order directional derivatives of the gravitational potential imposes very high requirements on the accuracy of deployed accelerometers which are beyond the limits of currently available sensors. For example, for orbital parameters and performance similar to those of the GOCE mission, observing third-order directional derivatives requires accelerometers with the noise level of \({\sim}10^{-17}\,\hbox {m}\,\hbox {s}^{-2}\) Hz\(^{-1/2}\).