Integration of GIS and Spatio-temporal Simulation Models: Interoperable Components for Different Simulation Strategies

New technologies for distributed object computing (DOC) challenge new approaches to integrate Geographical Information Systems (GIS) and simulation tools for spatio-temporal modeling. Several standardization initiatives work on specifying interfaces and defining service architectures for appropriate components. Currently, the most promising attempts are represented by the OpenGIS and the High Level Architecture. This paper gives a brief survey of the current state of the emerging standards, highlights the advantages and shortcomings of these attempts, and defines a possible framework for GIS-integrated simulation of spatio-temporal processes. Two abstraction layers provide interfaces for different GIS kernels and allow the integration of different simulation models. A prototype implementation of the specified interfaces shows the integration of analysis tools and two different simulation strategies: a mesoscale boundary layer model based on the Navier-Stokes equation and a cellular automata to simulate air pollutant dispersion.     

Uranium-lead isotope systematics of Mars inferred from the basaltic shergottite QUE 94201

Uranium-lead ratios (commonly represented as ²³⁸U/²⁰⁴Pb = μ) calculated for the sources of martian basalts preserve a record of petrogenetic processes that were active during early planetary differentiation and formation of martian geochemical reservoirs. To better define the range of μ values represented by the source regions of martian basalts, we completed U-Pb elemental and isotopic analyses on whole rock, mineral and leachate fractions from the martian meteorite Queen Alexandra Range 94201 (QUE 94201). The whole rock and silicate mineral fractions have unradiogenic Pb isotopic compositions that define a narrow range (²⁰⁶Pb/²⁰⁴Pb = 11.16-11.61). In contrast, the Pb isotopic compositions of weak HCl leachates are more variable and radiogenic. The intersection of the QUE 94201 data array with terrestrial Pb in ²⁰⁶Pb/²⁰⁴Pb-²⁰⁷Pb/²⁰⁴Pb-²⁰⁸Pb/²⁰⁴Pb compositional space is consistent with varying amounts of terrestrial contamination in these fractions. We calculate that only 1-7% contamination is present in the purified silicate mineral and whole rock fractions, whereas the HCl leachates contain up to 86% terrestrial Pb. This terrestrial Pb contamination generated a ²⁰⁶Pb-²⁰⁷Pb array in the QUE fractions that appears to represent an ancient age, which contrasts with a much younger crystallization age of 327 ± 10 Ma derived from Rb-Sr and Sm-Nd isochrons (Borg L. E., Nyquist L. E., Taylor L. A., Wiesmann H. and Shih C. -Y. (1997) Constraints on Martian differentiation processes from Rb-Sr and Sm-Nd isotopic analyses of the basaltic shergottite QUE 94201. Geochim. Cosmochim. Acta61, 4915-4931). Despite the contamination, and accepting 327 ± 10 Ma as the crystallization age, we use the U-Pb data to determine the initial ²⁰⁶Pb/²⁰⁴Pb of QUE 94201 to be 11.086 ± 0.008 and to calculate the μ value of its mantle source to be 1.82 ± 0.01. The μ value calculated for the QUE 94201 source is the lowest determined for any martian basalt source, and, when compared to the highest values determined for martian basalt sources, indicates that μ values in martian source reservoirs vary by at least a factor of two. Additionally, the range of source μ values indicates that the μ value of bulk silicate Mars is approximately three. The amount of variation in the μ values of the mantle sources (μ ∼ 2-4) is greater than can be explained by igneous processes involving silicate phases alone. We suggest the possibility that a small amount of sulfide crystallization may generate greater extents of U-Pb fractionation during formation of the mantle sources of martian basalts.     

Using rainfall-runoff modeling to interpret lake level data

Using water balance computations, the behavior of different kinds of lakes is discussed. Simple analytical expressions relating water level to hydrological conditions and lake bathymetry are given. The importance of knowing the river basin area when analyzing lake levels is stressed. A conceptual rainfall-runoff model including lake routing is used to simulate runoff and lake levels and to compute quasi-steady state conditions and long-term transient situations. It is suggested that models can be used to construct curves relating lake levels to precipitation and lake evaporation. By comparing with paleo-lake levels, the annual precipitation related to these levels can be found, provided information is available about the seasonal distribution of the precipitation.     

Effects of CO2 on benthic biota: An in situ benthic chamber experiment in Storfjorden (Norway)

Carbon capture and storage (CCS) methods, either sub-seabed or in ocean depths, introduces risk of CO₂ leakage and subsequent interaction with the ecosystem. It is therefore important to obtain information on possible effects of CO₂. In situ CO₂ exposure experiments were carried out twice for 10 days during 2005 using a Benthic Chamber system at 400 m depth in Storfjorden, Norway. pCO₂ in the water above the sediment in the chambers was controlled at approximately 500, 5000 and 20,000 μatm, respectively. This article describes the experiment and the results from measured the biological responses within the chamber sediments. The results show effects of elevated CO₂ concentrations on biological processes such as increased nanobenthos density. Methane production and sulphate reduction was enhanced in the approximately 5000 μatm chamber.     

Visualization-based analysis of structural and dynamical properties of simulated hydrous silicate melt

We have explored first-principles molecular dynamics simulation data for hydrous MgSiO₃ liquid (with 10 wt% water) to gain insight into its structural and dynamical behavior as a function of pressure (0–150 GPa) and temperature (2,000–6,000 K). By visualizing/analyzing a number of parameters associated with short- and mid-range orders, we have shown that the melt structure changes substantially on compression. The speciation of the water component at low pressures is dominated by the isolated structures (with over 90% hydrogen participated) consisting of hydroxyls, water molecules, O–H–O bridging and four-atom (O–H–O–H and H–O–H–O) groups, where every oxygen atom may be a part of polyhedron or free (i.e., bound to only magnesium atom). Hydroxyls favor polyhedral sites over magnesium sites whereas molecular water is almost entirely bound to magnesium sites, and also interpolyhedral bridging (Si–O–H–O–Si) dominates other types of bridging. Water content is shown to enhance and suppress, respectively, the proportions of hydroxyls and molecular water. As compression increases, these isolated structures increasingly combine with each other to form extended structures involving a total of five or more O and H atoms and also containing threefold coordination species, which together consume over 80% hydrogen at the highest compression studied. Our results show that water lowers the mean coordination numbers of different types including all cation–anion environments. The hydrous melt tends to be more tetrahedrally coordinated but with the Si–Si network being more disrupted compared to the anhydrous melt. Protons increase the content of non-bridging oxygen and decrease the contents of bridging oxygen as well as oxygen triclusters (present at pressures above 10 GPa). The calculated self-diffusion coefficients of all atomic species are enhanced in the presence of water compared to those of the anhydrous melt. This is consistent with the prediction that water depolymerizes the melt structure at all pressures. Our analysis also suggests that proton diffusion involves two processes—the transfer of H atoms (requiring the rupture and formation of O–H bonds) and the motion of hydroxyls as hydrogen carriers (requiring the rupture and formation of Si–O and/or Mg–O bonds). Both the processes are operative at low compression whereas only the first process is operative at high compression.     

Seismic‐stratigraphic constraints on the age of the Faroe Islands Basalt Group,Faroe–Shetland region,Northeast Atlantic Ocean

Lower Palaeogene extrusive igneous rocks of the Faroe Islands Basalt Group (FIBG) dominate the Faroese continental margin, with flood basalts created at the time of breakup and separation from East Greenland extending eastwards into the Faroe‐Shetland Basin. This volcanic succession was emplaced in connection with the opening of the NE Atlantic; however, consensus on the age and duration of volcanism remains lacking. On the Faroe Islands, the FIBG comprises four main basaltic formations (the pre‐breakup Lopra and Beinisvørð formations, and the syn‐breakup Malinstindur and Enni formations) locally separated by thin intrabasaltic sedimentary and/or volcaniclastic units. Offshore, the distribution of these formations remains ambiguous. We examine the stratigraphic framework of these rocks on the Faroese continental margin combining onshore (published) outcrop information with offshore seismic‐reflection and well data. Our results indicate that on seismic‐reflection profiles, the FIBG can be informally divided into lower and upper seismic‐stratigraphic packages separated by the strongly reflective A‐horizon. The Lower FIBG comprises the Lopra and Beinisvørð formations; the upper FIBG includes the Malinstindur and Enni formations. The strongly reflecting A‐horizon is a consequence of the contrast in properties of the overlying Malinstindur and underlying Beinisvørð formations. Onshore, the A‐horizon is an erosional surface, locally cutting down into the Beinisvørð Formation; offshore, we have correlated the A‐horizon with the Flett unconformity, a highly incised, subaerial unconformity, within the juxtaposed and interbedded sedimentary fill of the Faroe‐Shetland Basin. We refer to this key regional boundary as the A‐horizon/Flett unconformity. The formation of this unconformity represents the transition from the pre‐breakup to the syn‐breakup phase of ocean margin development in the Faroe–Shetland region. We examine the wider implications of this correlation considering existing stratigraphic models for the FIBG, discussing potential sources of uncertainty in the correlation of the lower Palaeogene succession across the Faroe–Shetland region, and implications for the age and duration of the volcanism.     

Frequency locking in a periodically forced thermal two-phase oscillator

Welander (1977) showed that a fluid, heated internally and cooled to freezing from above, can for a certain parameter range and for steady thermal forcing demonstrate oscillatory behaviour with ice-covered and ice-free states alternating. This concept has been extended to encompass the geophysically important case of cyclic cooling. Numerical investigations have revealed that periodically forced either shows frequency locking or quasiperiodic behaviour.