Tight dual models of pore spaces

The pore throats in a porous medium control permeability, drainage, and straining through their pore scale geometry and through the way they are connected via pore bodies on the macroscale. Likewise, imbibition is controlled through the geometry of the pore bodies (pore scale) and through the way the pore bodies are connected via pore throats on the macroscale. In an effort to account for both scales at the same time we recently introduced an image-based model for pore spaces that consists of two parts related by duality: (1) a decomposition of a polyhedral pore space into polyhedral pore bodies separated by polygonal pore throats and (2) a polygonal pore network that is homotopy equivalent to the pore space. In this paper we stick to the dual concept while amending the definition of the pore throats and, as a consequence, the other elements of the dual model. Formerly, the pore throats consisted of single two-dimensional Delaunay cells, while they now usually consist of more than one two-dimensional Delaunay cell and extend all the way into the narrowing ends of the pore channel cross sections. This is the first reason for naming the amended dual model “tight”. The second reason is that the formation of the pore throats is now guided by an objective function that always attains its global optimum (tight optimization). At the end of the paper we report on simulations of drainage performed on tight dual models derived from simulated sphere packings and 3D gray-level images. The C-code for the generation of the tight dual model and the simulation of drainage is publicly available at https://jshare.johnshopkins.edu/mhilper1/public_html/tdm.html.     

The gravitational potential and its derivatives of a right rectangular prism with depth-dependent density following an n-th degree polynomial

The direct gravity problem and its solution belong to the basis of the gravimetry. The solutions of this problem are well known for wide class of the source bodies with the constant density contrast. The non-uniform density approximation leads to the relatively complicated mathematical formalism. The analytical solutions for this type of sources are rare and currently these bodies are very useful in the gravimetrical modeling. The solution for the vertical component of the gravitational attraction vector for the 3D right rectangular prism is known in the geophysical literature for the density variations described by the 3-rd degree polynomial. We generalized this solution for an n-th degree, not only for the vertical component, but for the horizontal components, the second-order derivatives and the potential as well. The 2D modifications of all given formulae are presented, too. The presented general solutions, which involve a hypergeometric functions, can be used as they are, or as an auxiliary tool to derive desired solution for the given degree of the density polynomial as a sum of the elementary functions. The pros-and-cons of these approaches (the complexity of the programming codes, runtimes) are discussed, too.     

Experiences with the use of mass-density maps in residual gravity forward modelling

In many modern local and regional gravity field modelling concepts, the short-wavelength gravitational signal modeled by the residual terrain modelling (RTM) technique is used to augment global geopotential models, or to smooth observed gravity prior to data gridding. In practice, the evaluation of RTM effects mostly relies on a constant density assumption, because of the difficulty and complexity of obtaining information on the actual distribution of density of topographic masses. Where the actual density of topographic masses deviates from the adopted value, errors are present in the RTM mass-model, and hence, in the forward-modelled residual gravity field. In this paper we attempt to overcome this problem by combining the RTM technique with a high-resolution mass-density model. We compute RTM gravity quantities over New Zealand, with different combinations of elevation models and mass-density assumptions using gravity and GPS/levelling measurements, precise terrain and bathymetry models, a high-resolution mass-density model and constant density assumptions as main input databases. Based on gravity observations and the RTM technique, optimum densities are detected for North Island of ~2500 kg m^?3, South Island of ~2600 kg m^?3, and the whole New Zealand of ~2590 kg m^?3. Comparison among the three sets of residual gravity disturbances computed from different mass-density assumptions show that, together with a global potential model, the high-resolution New Zealand density model explains ~89.5% of gravitational signals, a constant density assumption of 2670 kg m^?3 explains ~90.2%, while a regionally optimum mass-density explains ~90.3%. Detailed comparison shows that the New Zealand density model works best over areas with small residual heights. Over areas with larger residual heights, subsurface density variations appear to affect the residual gravity disturbance. This effect is found to reach about 30 mGal over Southern Alpine Fault. In order to improve the RTM modelling with mass-density maps, a higher-quality mass-density model that provides radially varying mass-density data would be desirable.     

Can sand dunes be used to study historic storm events?

Knowing the long‐term frequency of high magnitude storm events that cause coastal inundation is critical for present coastal management, especially in the context of rising sea levels and potentially increasing frequency and severity of storm events. Coastal sand dunes may provide a sedimentary archive of past storm events from which long‐term frequencies of large storms can be reconstructed. This study uses novel portable optically stimulated luminescence (POSL) profiles from coastal dunes to reconstruct the sedimentary archive of storm and surge activity for Norfolk, UK. Application of POSL profiling with supporting luminescence ages and particle size analysis to coastal dunes provides not only information of dunefield evolution but also on past coastal storms. In this study, seven storm events, two major, were identified from the dune archive spanning the last 140 years. These appear to correspond to historical reports of major storm surges. Dunes appear to be only recording (at least at the sampling resolution used here) the highest storm levels that were associated with significant flooding. As such the approach seems to hold promise to obtain a better understanding of the frequency of large storms by extending the dune archive records further back to times when documentation of storm surges was sparse. © 2017 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

Turbulent Heating in the Accelerating Region Using a Multishell Model

Solar energetic particles (SEPs) are released into the heliosphere by solar flares and coronal mass ejections (CMEs). They are mostly protons, with smaller amounts of heavy ions from helium to iron, and lesser amounts of species heavier than iron. The spectra of heavy ions have been previously studied mostly by using the fluence of the particles in an event-integrated spectrum in a small number of spectral snapshots. In this article, we analyze the temporal evolution of the heavy-ion spectra using two large SEP events (27 January 2012 and 7 January 2014) from the Solar TErrestrial Relations Observatory (STEREO) era using Advanced Composition Explorer (ACE) Solar Isotope Spectrometer (SIS) and Ultra Low Energy Isotope Spectrometer (ULEIS), Energetic Particles: Acceleration, Composition and Transport (EPACT) onboard Wind, and the STEREO-A (Ahead) and -B (Behind) Low-Energy Telescope (LET) and Suprathermal Ion Telescope (SIT) instruments, taking a large number of snapshots covering the temporal evolution of the event. We find large differences in the spectra of the ions after the main flux enhancement in terms of the grouping of similar species, but also in terms of the location of the instruments. Although it is somewhat less noticeable than in the case of the temporal evolution of protons (Doran and Dalla, Solar Phys. 291, 2071, 2016), we observe a wave-like pattern travelling through the heavy ion spectra from the highest energies to the lowest, creating an “arch” structure that later straightens into a power law after 18 to 24 hours.

     

土体触探静力加载试验的FLAC模拟模型选择

本文将触探静力加载试验作为轴对称问题处理,建立了两种不同FLAC轴对称模型,等荷载模型和等位移模型。研究表明,两种模型的FLAC模拟结果基本一致,一般条件下,应用简化的、运算速度快的等荷载模型可以满足试验要求。     

Dispersant Use and a Bioremediation Strategy as Alternate Means of Reducing Impacts of Large Oil Spills on Mangroves: The Gladstone Field Trials

Over a three-year period (1995–1998), we studied short-term effects of dispersant use and a bioremediation strategy in two consecutive field trials in sub-tropical Australian mangroves. In each case, weathered oil was applied, and a large spill simulated, in mature Rhizophora stylosa trees around 4–9 m tall. In the first trial, we used Gippsland light crude oil with or without dispersant, Corexit 9527. In the second, a bioremediation strategy followed application of Gippsland oil or Bunker C fuel oil. Bioremediation involved forced aeration with supplemental application of nutrients. Dispersant use had an overall positive benefit shown as reduced tree mortality. By contrast, there was no apparent reduction in mortality of trees with bioremediation. However, one year after oiling, leaf densities of surviving trees were greater in bioremediation plots than in controls, and less in oil-only plots. These and other results have been incorporated into spill response management strategies in Australia.     

Fish assemblages as indicators of water quality in the middle Thames estuary, England (1980–1989)

Fish abundance and environmental data collected over ten years (1980–1989) from the middle Thames estuary, England, were analyzed to detect temporal trends in fish populations and relationship with environmental parameters, and to assess water quality. Fish were collected from the cooling water intake screens of West Thurrock power station, situated 35.5 km below London Bridge, in the mid-estuary. Marine species abundance were highly seasonal, with peaks in December–March for herring (Clupea harengus), sprat (Sprattus sprattus), 3-spined-stickleback (Gasterosteus aculeatus), and poor cod (Trisopterus minutus); July–August for flounder (Platichthys flesus); and September–December for sand goby (Pomatoschistus minutus), whiting (Merlangius merlangus), bass (Dicentrarchus labrax), plaice (Pleuronectes platessa), and dab (Limanda limanda). Bimodal seasonal patterns of peaks or unclear seasonality in abundance characterized marine estuarine-dependent sole (Solea solea), Nilsson's pipefish (Syngnathus rostelattus) (April/May and September/October), and pouting (Trisopterus luscus) (May and November/December); the estuarine smelt (Osmerus eperlanus) (October and January) and the catadromous eel (Anguilla anguilla) (June and October). There was substantial variation in the abundance of common species over the period of ten years, with herring, sand goby, flounder, and plaice showing a stable abundance in 1980–1984, increasing sharply in 1985–1986, and then decreasing successively through the remainder of the decade (1987–1989). The first half of the decade was a period of higher abundance for less tolerant species such as smelt, sprat, and poor cod, while the second half showed higher abundances of species tolerant to harsh environmental conditions such as sand goby, flounder, eel, and plaice. A general pattern of stable fish populations with a slight trend of deterioration was found to emerge over the years, related to the number of species and quantities of common species. Multivariate techniques of principal component and canonical correspondence ordinations were used for assessing relationships between fish populations abundance and environmental variables. The most significant environmental variables correlated with fish species were temperature and dissolved oxygen. High abundances of flounder were associated with high temperature, while high abundance of poor cod, sprat, herring, and 3-spined-stickleback were associated with high dissolved oxygen, flow, ammonical nitrogen, and low temperature. Plaice, whiting, sand goby, bass, and dab were preferentially found in high salinity and suspended solids, while smelt and sole were likely to prefer average values or showed no clear preferences.     

AMG preconditioning for sedimentary basin simulations in Temis calculator

Basin modelling aims at reconstructing the time evolution of a sedimentary basin in order to make quantitative predictions of geological phenomena leading to hydrocarbon accumulations. It accounts for porous medium compaction, heat transfer, hydrocarbon generation and migration. These physical phenomena are modelled by partial differential equations [Schneider, F., Wolf, S., Faille, I., Pot, D. 2000. A 3D basin model for hydrocarbon potential evaluation: application to Congo offshore. Oil and Gas Science Technologie, Rev-IFP 55, 3–13] representing the mass balances of solid, fluids (water and oil) coupled with Darcy's law and a compaction law. These equations are discretized using a cell-centered Finite Volume method in space and an implicit Euler integration in time. At each time step, the resulting nonlinear system is solved using Newton's method ending up at each Newton iteration with the solution of a linear system which represents the most cpu-time-consuming part of the basin simulator.     

Assimilation of preexisting Pleistocene intrusions at Long Valley by periodic magma recharge accelerates rhyolite generation: rethinking the remelting model

Rhyolite flows and tuffs from the Long Valley area of California, which were erupted over a two-million-year time period, exhibit systematic trends in Nd, Hf, and Pb isotopes, trace element composition, erupted volume, and inferred magma residence time that provide evidence for a new model for the production of large volumes of silica-rich magma. Key constraints come from geochronology of zircon crystal populations combined with a refined eruption chronology from Ar–Ar geochronology; together these data give better estimates of magma residence time that can be evaluated in the context of changing magma compositions. Here, we report Hf, Nd, and Sr isotopes, major and trace element compositions, ⁴⁰Ar/³⁹Ar ages, and U–Pb zircon ages that combined with existing data suggest that the chronology and geochemistry of Long Valley rhyolites can be explained by a dynamic interaction of crustal and mantle-derived magma. The large volume Bishop Tuff represents the culmination of a period of increased mantle-derived magma input to the Long Valley volcanic system; the effect of this input continued into earliest postcaldera time. As the postcaldera evolution of the system continued, new and less primitive crustal-derived magmas dominated the system. A mixture of varying amounts of more mafic mantle-derived and felsic crustal-derived magmas with recently crystallized granitic plutonic materials offers the best explanation for the observed chronology, secular shifts in Hf and Nd isotopes, and the apparently low zircon crystallization and saturation temperatures as compared to Fe–Ti oxide eruption temperatures. This scenario in which transient crustal magma bodies remained molten for varying time periods, fed eruptions before solidification, and were then remelted by fresh recharge provides a realistic conceptual framework that can explain the isotopic and geochemical evidence. General relationships between crustal residence times and magma sources are that: (1) precaldera rhyolites had long crustal magma residence times and high crustal affinity, (2) the caldera-related Bishop Tuff and early postcaldera rhyolites have lower crustal affinity and short magma residence times, and (3) later postcaldera rhyolites again have stronger crustal signatures and longer magma residence times.