A general reduction scheme for reactive transport in porous media

We present a method to transform the governing equations of multispecies reactive transport in porous media. The reformulation leads to a smaller problem size by decoupling of equations and by elimination of unknowns, which increases the efficiency of numerical simulations. The reformulation presented here is a generalization of earlier works. In fact, a whole class of transformations is now presented. This class is parametrized by the choice of certain transformation matrices. For specific choices, some known formulations of reactive transport can be retrieved. Hence, the software based on the presented transformation can be used to obtain efficiency comparisons of different solution approaches. For our efficiency tests, we use the MoMaS benchmark problem on reactive transport.     

Fault displacement modelling using 3D vector fields

In history matching and sensitivity analysis, flexibility in the structural modelling is of great importance. The ability to easily generate multiple realizations of the model will have impact both on the updating workflow in history matching and uncertainty studies based on Monte Carlo simulations. The main contribution to fault modelling by the work presented in this paper is a new algorithm for calculating a 3D displacement field applicable to a wide range of faults due to a flexible representation. This gives the possibility to apply this field to change the displacement and thereby moving horizons and fault lines. The fault is modelled by a parametric format where the fault has a reference plane defined by a centre point, dip and strike angles. The fault itself is represented as a surface defined by a function z = f(x,y), where x, y and z are coordinates local to the reference plane, with the z-axis being normal to the plane. The displacement associated with the fault outside the fault surface is described by a 3D vector field. The displacement on the fault surface can be found by identifying the intersection lines between horizons and the fault surface (fault lines), and using kriging techniques to fill in values at all points on the surface. Away from the fault surface the displacement field is defined by a monotonic decreasing function which ensures zero displacement at a specified distance from the fault. An algorithm is developed where the displacement can be increased or decreased according to user-defined parameters. This means that the whole displacement field is changed and points on horizons around the fault can be moved accordingly by applying the modified displacement field on them. The interaction between several faults influencing the same points is taken care of by truncation rules and the ordering of the faults. The method is demonstrated on a realistic synthetic case based on a real reservoir.     

Soft-sediment deformation of Late Pleistocene sediments along the southwestern coast of the Baltic Sea (NE Germany)

A 1,460-m-long profile of a Late Glacial subglacial, glacio-fluvial, glacio-limnic and glacio-deltaic sequence exposed at a cliff section on Usedom Island (SW Baltic Sea coast) is described. The sequence is up to 31 m thick and shows sedimentary structures typical of a glacial setting. Soft-sediment deformation is encountered and is associated with changes in lithology. These deformations include liquefaction, slumping, and faulting. As the most plausible cause, earthquake-induced shaking is discussed. The associated neotectonic activity is seen as a consequence of the postglacial isostatic crustal rebound. As the deglaciation earthquake ratio diminishes with time and as the rebound is phasing out, no large earthquakes are anticipated for northern Germany, although in conclusion the lithosphere of the North German Basin has to be regarded as weakened by repeated ice loading and deloading.     

Quaternary coastline evolution of Lake Ohrid (Macedonia/Albania)

Lake Ohrid (between FYR of Macedonia and Albania), situated in an active tectonic region of the Balkanides, is characterized by N ?? S trending active faults. To reconstruct the Holocene shoreline evolution we investigated the coastline using sediment cores and geophysical methods to image sedimentary and tectonic structures. We revealed areas of differing sedimentation regimes. The plains north and south of the lake are dominated by clastic input related to climate variations and uplift/erosion, whereas the steep western and eastern margins are controlled by recent tectonics. Furthermore, no evidence for a much higher lake-level during the Holocene was found in the plains north and south of the lake, except rare temporary floodings. This is supported by mappings of the limestone cliffs around Lake Ohrid, which yielded no evidence for abrasional platforms or notches as indicators for past highstands.     

Late Quaternary marine-based Kara Sea ice sheets: a review of terrestrial stratigraphic data highlighting their formation

Reconstructions of the Late Quaternary glacial history of the Kara Sea area show repeated build-up of ice-sheet domes over the shallow epicontinental Kara Sea. Inferred ice divides were situated over the central Kara Sea, and the ice sheet repeatedly inundated the surrounding coastal areas of western Siberia. Geological fingerprinting of the Kara Sea ice sheet include end moraine zones, raised beaches, tills, glaciotectonic deformations and coarsening-upward sediment sequences, reflecting isostatic rebound cycles. This paper reviews evidence from several areas along the perimeter of the Kara Sea, suggesting that peripheral sites were critical for the initiation of the large Kara Sea ice sheet. Ice-sheet inception progressed with the formation of local ice caps that later coalesced on the adjacent shelf with globally falling sea levels, eventually merging and growing into a large ice dome.     

Long-term changes of (polar) mesosphere summer echoes

Strong VHF radar echoes have been observed not only during summer months at polar latitudes (polar mesosphere summer echoes, PMSE) but also at middle latitudes (mesosphere summer echoes, MSE). These echoes are closely connected with small ice particles, thus containing information about mesospheric temperature and water vapour content. But the (P)MSE also depend on the ionisation due to solar wave radiation and precipitating high energetic particles. Observations with VHF radars at Andenes (69.3°N; 16.0°E) since 1994 and at Kühlungsborn (54.6°N; 11.8°E) since 1998 are used for investigations of the solar and geomagnetic control of the (P)MSE as well as of possible long-term changes. The (P)MSE are positively correlated with the solar Lyman α radiation and the geomagnetic activity and have slightly positive trends. Due to the limited measuring period, the significance levels of the detected (P)MSE trends are small. Positive trends in noctilucent clouds (NLC) and polar mesospheric clouds (PMC) are in general agreement with (P)MSE trends.     

Spatio-temporal prediction of daily temperatures using time-series of MODIS LST images

A computational framework to generate daily temperature maps using time-series of publicly available MODIS MOD11A2 product Land Surface Temperature (LST) images (1 km resolution; 8-day composites) is illustrated using temperature measurements from the national network of meteorological stations (159) in Croatia. The input data set contains 57,282 ground measurements of daily temperature for the year 2008. Temperature was modeled as a function of latitude, longitude, distance from the sea, elevation, time, insolation, and the MODIS LST images. The original rasters were first converted to principal components to reduce noise and filter missing pixels in the LST images. The residual were next analyzed for spatio-temporal auto-correlation; sum-metric separable variograms were fitted to account for zonal and geometric space-time anisotropy. The final predictions were generated for time-slices of a 3D space-time cube, constructed in the R environment for statistical computing. The results show that the space-time regression model can explain a significant part of the variation in station-data (84%). MODIS LST 8-day (cloud-free) images are unbiased estimator of the daily temperature, but with relatively low precision (±4.1°C); however their added value is that they systematically improve detection of local changes in land surface temperature due to local meteorological conditions and/or active heat sources (urban areas, land cover classes). The results of 10–fold cross-validation show that use of spatio-temporal regression-kriging and incorporation of time-series of remote sensing images leads to significantly more accurate maps of temperature than if plain spatial techniques were used. The average (global) accuracy of mapping temperature was ±2.4°C. The regression-kriging explained 91% of variability in daily temperatures, compared to 44% for ordinary kriging. Further software advancement—interactive space-time variogram exploration and automated retrieval, resampling and filtering of MODIS images—are anticipated.     

Calibrating adaptive antenna arrays for high-integrity GPS

A major challenge in using GPS guidance for aircraft final approach and landing is to reject interference that can jam reception of the GPS signals. Antenna arrays, which use space–time adaptive processing (STAP), significantly improve the signal to interference plus noise ratio, but at the possible expense of distorting the received signals, leading to timing biases that may degrade navigation performance. Rather than a sophisticated calibration approach to remove biases introduced by STAP, this paper demonstrates that a relatively compact calibration strategy can substantially reduce navigation biases, even under elevated interference conditions. Consequently, this paper develops an antenna bias calibration strategy for two classes of adaptive array algorithm and validates this method using both simulated and experimental data with operational hardware in the loop. A proof-of-concept system and an operational prototype are described, which implement the adaptive antenna algorithms and deterministic corrections. This investigation demonstrates that systems with adaptive antenna arrays can approach the accuracy and integrity requirements for automatic aircraft landing, and in particular for sea-based landing on board aircraft carriers, while simultaneously providing significant attenuation of interference. Evidence suggests that achieving these goals is possible with minimal restrictions on system hardware configuration—specifically, limitations on the permissible level of antenna anisotropy and the use of sufficient analog-to-digital converter resolution.     

A solubility and surface complexation study of a non-stoichiometric hydroxyapatite

The dissolution and surface complexation of a non-stoichiometric hydroxyapatite (Ca_8.4(HPO₄)_1.6(PO₄)_4.4(OH)_0.4), (HAP) was studied in the pH range 3.5-10.5, at 25 °C in 0.1 M Na(Cl). The results from well-equilibrated batch experiments, potentiometric titrations, and zeta-potential measurements were combined with information provided by Attenuated Total Reflectance Fourier Transform Infrared (ATR-FTIR) spectroscopy and X-ray Photoelectron Spectroscopy (XPS). The information from the analyses was used to design an equilibration model that takes into account dissolution, surface potential, solution and surface complexation, as well as possible phase transformations. The results from the XPS measurements clearly show that the surface of the mineral has a different composition than the bulk and that the Ca/P ratio of the surface layer is 1.4 ± 0.1. This ratio was also found in solution in the batches equilibrated at low pH where the dominating reaction is dissolution. In the batches equilibrated at near neutral pH values, however, the Ca/P ratio in solution attains values as high as 25, which is due to re-adsorption of phosphate ions to the HAP surface. The total concentration of protons as well as the total concentration of dissolved calcium and phosphate in solution were used to calculate a model for the dissolution and surface complexation of HAP. The constant capacitance model was applied in designing the following surface complexation model: