Towards an integrated arid zone water management using simulation-based optimisation

For ensuring both optimal sustainable water resources management and long-term planning in a changing arid environment, we propose an integrated Assessment-, Prognoses-, Planning- and Management tool (APPM). The new APPM integrates the complex interactions of the strongly nonlinear meteorological, hydrological and agricultural phenomena, considering the socio-economic aspects. It aims at achieving best possible solutions for water allocation, groundwater storage and withdrawals including saline water management together with a substantial increase of the water use efficiency employing novel optimisation strategies for irrigation control and scheduling. To obtain a robust and fast operation of the water management system, it unites process modeling with artificial intelligence tools and evolutionary optimisation techniques for managing both water quality and quantity. We demonstrate some key components of our methodology by an exemplary application to the south Al-Batinah region in the Sultanate of Oman which is affected by saltwater intrusion into a coastal aquifer due to excessive groundwater withdrawal for irrigated agriculture. We show the effectiveness and functionality of a new simulation-based water management system for the optimisation and evaluation of different irrigation practices, crop pattern and resulting abstraction scenarios. The results of several optimisation runs indicate that due to contradicting objectives, such as profit-oriented agriculture versus aquifer sustainability only a multi-objective optimisation can provide sustainable solutions for the management of the water resources in respect of the environment as well as the socio-economic development.     

Chaos formation by sublimation of volatile-rich substrate: Evidence from Galaxias Chaos, Mars

Galaxias Chaos deviates significantly from other chaotic regions due to the lack of associated outflow channels, lack of big elevation differences between the chaos and the surrounding terrain and due to gradual trough formation.A sequence of troughs in different stages is observed, and examples of closed troughs within blocks suggest that the trough formation is governed by a local stress field rather than a regional stress field. Moreover, geomorphic evidence suggests that Galaxias Chaos is capped by Elysium lavas, which superpose an unstable subsurface layer that causes chaotic tilting of blocks and trough formation.Based on regional mapping we suggest a formation model, where Vastitas Borealis Formation embedded between Elysium lavas is the unstable subsurface material, because gradual volatile loss causes shrinkage and differential substrate movement. This process undermines the lava cap, depressions form and gradually troughs develop producing a jigsaw puzzle of blocks due to trough coalescence.Observations of chaos west of Elysium Rise indicate that this process might have been widespread along the contact between Vastitas Borealis Formation and Elysium lavas. However, the chaotic regions have probably been superposed by Elysium/Utopia flows to the NW of Elysium Rise, and partly submerged with younger lavas to the west.     

Autocompaction of shallow silty salt marsh clay

Empirically derived algorithms are suggested by means of which it is possible to evaluate relations between physical properties in silty salt marsh clay. A relation between loss on ignition, sand content and bulk dry density in the topmost 5 cm was found to explain 84% of the variation of the latter. The bulk dry density in the topmost ≈ 0.5 m was found to vary as a logarithmic function of depth below the surface controlled by two constants directly dependent on the bulk dry density in the topmost 5 cm. Integration of this function illustrates that mass depth (MSD_z) down core can be directly related to the bulk dry density of the surface layer. For sediments with a surface bulk dry density between 400 kg m^? 3 and 800 kg m^? 3 the porosity was found to vary (+/?10%) with MSD_z as another logarithmic function. The results allow for an evaluation of the use of marker horizons in the topmost layers. Right after such a layer has been spread on the surface, sedimentation measured in length per time (accretion), will be significantly larger than for subsequent periods. Using marker horizons, it is therefore not trivial to know the level of the marker, as deeper located markers will indicate successively smaller accretion rates with the same sediment input.     

Geochemical controls on sediment reactivity and buffering processes in a heterogeneous aquifer

The injection and recovery of oxic water into deep anoxic aquifers may help to alleviate short- and long-term imbalance between freshwater supply and demand. The extent and structure of physical and geochemical heterogeneity of the aquifer will impact the water quality evolution during injection, storage and recovery. Water–sediment interactions within the most permeable parts of the aquifer, where the bulk of the injectant will penetrate, may dominate, however, water quality may also be impacted by interactions within the finer-grained, less permeable but potentially highly reactive media. In this study, the heterogeneity of the reductive capacity of an aquifer selected for water reuse projects was characterised, the amount, type and reactivity of the sedimentary reductants present determined, and the relationship between reductive capacity and sedimentary lithologies quantified. The average potential reductive capacities (PRC_TOT), based on total organic C and pyrite concentrations of the sediment, were quantified for sands (382 μmol O₂ g⁻¹), clays (1522 μmol O₂ g⁻¹), and silts (1957 μmol O₂ g⁻¹). Twenty-seven samples, spanning the three different lithologies, were then incubated for 50 days and the measured reductive capacities (MRC) determined for the sands (29.2 μmol O₂ g⁻¹), silts (136 μmol O₂ g⁻¹), and clays (143 μmol O₂ g⁻¹). On average, the MRC were 10% of the PRC_TOT. The main consumers of O₂ were pyrite (20–100%), sedimentary organic matter (SOM; 3–56%), siderite (3–28%) and Fe(II)-aluminosilicates (8–55%). The incubation data plus hydrogeochemical modelling, indicated that pH-buffering was controlled firstly by dissolution of trace level carbonates, followed by dissolution of feldspars. Zinc, Co, Ni, Cd and Pb were readily mobilized during incubation.     

Evaluating the salinity distribution of a shallow coastal aquifer by vertical multielectrode profiling (Denmark)

A monitoring system, including five groups of piezometers and five vertical multielectrode profiling probes (VMEP), has been installed in an aquifer beneath a coastal dune in Denmark. In order to assess the salinity distribution within the aquifer, geoelectrical data were gathered in March, June and September 2008, by measuring a dipole-dipole and gradient array using multielectrode profiling. Interpretation of the processed resistivity data was performed by regularized inversion using a one-dimensional, horizontally layered model of formation resistivity. The standard deviation on estimated layer log-resistivity was 0.01–0.03. By estimating two parameters of a power function, observed fluid conductivities derived from samples of porewater were related to corresponding estimated formation resistivities. The conductivity profiles correlate with a winter situation in March with high sea level, active recharge and significant wave activity, causing increased hydraulic heads, a thicker freshwater lens and salt water overlying freshwater close to the sea. In June, the thickness of the freshwater lens is reduced due to less recharge and prevailing offshore winds, imposing density-stable conditions and a sharper transition between fresh and brackish water. During the autumn, aquifer recharge is enhanced and hydraulic heads increase, resulting in a thicker freshwater lens.     

An “ESA-Affordable” Laue-lens

With ESA’s INTEGRAL mission gamma-ray astronomy has advanced to the point where major scientific advances must be expected from detailed studies of the many new point sources. The interest in developing focusing telescopes operating in the soft gamma-ray regime up to 1 MeV is therefore mounting rapidly. Telescopes based on Laue diffraction of gamma-rays from crystals appear as one promising route, although the practical difficulties of realizing a large scale Laue lens are certainly not small. In this paper I have attempted to develop an optimized lens design considering the size and mass constraints of a specific medium size launch vehicle. The introduction of the lens mass as a primary design driver has some surprising effects for the choice of material for the crystals and new tradeoff considerations are introduced.     

Tsunami generation and propagation from the Mjølnir asteroid impact

Abstract— In the late Jurassic period, about 142 million years ago, an asteroid hit the shallow paleo‐Barents Sea, north of present‐day Norway. The geological structure resulting from the impact is today known as the Mjølnir crater. The present work attempts to model the generation and the propagation of the tsunami from the Mjølnir impact. A multi‐material hydrocode SOVA is used to model the impact and the early stages of tsunami generation, while models based on shallow‐water theories are used to study the subsequent wave propagation in the paleo‐Barents Sea. We apply several wave models of varying computational complexity. This includes both three‐dimensional and radially symmetric weakly dispersive and nonlinear Boussinesq equations, as well as equations based on nonlinear ray theory. These tsunami models require a reconstruction of the bathymetry of the paleo‐Barents Sea. The Mjølnir tsunami is characteristic of large bolides impacting in shallow sea; in this case the asteroid was about 1.6 km in diameter and the water depth was around 400 m. Contrary to earthquake‐ and slide‐generated tsunamis, this tsunami featured crucial dispersive and nonlinear effects: a few minutes after the impact, the ocean surface was formed into an undular bore, which developed further into a train of solitary waves. Our simulations indicate wave amplitudes above 200 m, and during shoaling the waves break far from the coastlines in rather deep water. The tsunami induced strong bottom currents, in the range of 30–90 km/h, which presumably caused a strong reworking of bottom sediments with dramatic consequences for the marine environment.     

The geography of global urban greenhouse gas emissions: an exploratory analysis

The purpose of this paper is to describe global urban greenhouse gas emissions by region and sector, examine the distribution of emissions through the urban-to-rural gradient, and identify covariates of emission levels for our baseline year, 2000. We use multiple existing spatial databases to identify urban extent, greenhouse gas emissions (CO₂, N₂O, CH₄ and SF₆) and covariates of emissions in a “top-down” analysis. The results indicate that urban activities are significant sources of total greenhouse gas emissions (36.8 and 48.6 % of total). The urban energy sector accounts for between 41.5 and 66.3 % of total energy emissions. Significant differences exist in the urban share of greenhouse gas emissions between developed and developing countries as well as among source sectors for geographic regions. The 50 largest urban emitting areas account for 38.8 % of all urban greenhouse gas emissions. We find that greenhouse gas emissions are significantly associated with population size, density, growth rates, and per capita income. Finally, comparison of our results to “bottom-up” estimates suggest that this research’s data and techniques are best used at the regional and global scales.