Time-lapse pressure tomography for characterizing CO2 plume evolution in a deep saline aquifer

A time-lapse pressure tomography inversion approach is applied to characterize the CO2 plume development in a virtual deep saline aquifer. Deep CO2 injection leads to flow properties of the mixed-phase, which vary depending on the CO2 saturation. Analogous to the crossed ray paths of a seismic tomographic experiment, pressure tomography creates streamline patterns by injecting brine prior to CO2 injection or by injecting small amounts of CO2 into the two-phase (brine and CO2) system at different depths. In a first step, the introduced pressure responses at observation locations are utilized for a computationally rapid and efficient eikonal equation based inversion to reconstruct the heterogeneity of the subsurface with diffusivity (D) tomograms. Information about the plume shape can be derived by comparing D-tomograms of the aquifer at different times. In a second step, the aquifer is subdivided into two zones of constant values of hydraulic conductivity (K) and specific storage (Ss) through a clustering approach. For the CO2 plume, mixed-phase K and Ss values are estimated by minimizing the difference between calculated and “true” pressure responses using a single-phase flow simulator to reduce the computing complexity. Finally, the estimated flow property is converted to gas saturation by a single-phase proxy, which represents an integrated value of the plume. This novel approach is tested first with a doublet well configuration, and it reveals a great potential of pressure tomography based concepts for characterizing and monitoring deep aquifers, as well as the evolution of a CO2 plume. Still, field-testing will be required for better assessing the applicability of this approach.

Retrieving the paleoclimatic signal from the deeper part of the EPICA Dome C ice core

An important share of paleoclimatic information is buried within the lowermost layers of deep ice cores. Because improving our records further back in time is one of the main challenges in the near future, it is essential to judge how deep these records remain unaltered, since the proximity of the bedrock is likely to interfere both with the recorded temporal sequence and the ice properties. In this paper, we present a multiparametric study (δD-δ18Oice, δ18Oatm, total air content, CO2, CH4, N2O, dust, high-resolution chemistry, ice texture) of the bottom 60 m of the EPICA (European Project for Ice Coring in Antarctica) Dome C ice core from central Antarctica. These bottom layers were subdivided into two distinct facies: the lower 12 m showing visible solid inclusions (basal dispersed ice facies) and the upper 48 m, which we will refer to as the "basal clean ice facies". Some of the data are consistent with a pristine paleoclimatic signal, others show clear anomalies. It is demonstrated that neither large-scale bottom refreezing of subglacial water, nor mixing (be it internal or with a local basal end term from a previous/initial ice sheet configuration) can explain the observed bottom-ice properties. We focus on the high-resolution chemical profiles and on the available remote sensing data on the subglacial topography of the site to propose a mechanism by which relative stretching of the bottom-ice sheet layers is made possible, due to the progressively confining effect of subglacial valley sides. This stress field change, combined with bottom-ice temperature close to the pressure melting point, induces accelerated migration recrystallization, which results in spatial chemical sorting of the impurities, depending on their state (dissolved vs. solid) and if they are involved or not in salt formation. This chemical sorting effect is responsible for the progressive build-up of the visible solid aggregates that therefore mainly originate "from within", and not from incorporation processes of debris from the ice sheet's substrate. We further discuss how the proposed mechanism is compatible with the other ice properties described. We conclude that the paleoclimatic signal is only marginally affected in terms of global ice properties at the bottom of EPICA Dome C, but that the timescale was considerably distorted by mechanical stretching of MIS20 due to the increasing influence of the subglacial topography, a process that might have started well above the bottom ice. A clear paleoclimatic signal can therefore not be inferred from the deeper part of the EPICA Dome C ice core. Our work suggests that the existence of a flat monotonic ice–bedrock interface, extending for several times the ice thickness, would be a crucial factor in choosing a future "oldest ice" drilling location in Antarctica.

The 1987 debris flows in Switzerland: documentation and analysis

A great number of debris flows occurred during the flood catastrophes of the summer of 1987 in the Swiss Alps. Aerial photography, field investigations and eyewitness accounts documented and analysed the events. As an example of the reconstructed major events, the large debris flow in the Varuna valley involved an estimated peak discharge between 400 and 800 m3/s and an event magnitude of 200,000 m3. Several single pulses were observed; the duration of each of them appeared to be not more than a few minutes. Apart from incision into weak bedrock, the maximum erosion depth seemed to depend on the channel gradient. Based on approximately 600 events, typical starting zones and rainfall conditions are discussed with regard to the triggering conditions. Existing and new empirical formulae are proposed to estimate the most important flow parameters. These values are compared to debris flow data from Canada and Japan.

Reconstruction of in-situ porosity and porewater compositions of low-permeability crystalline rocks: Magnitude of artefacts induced by drilling and sample recovery

Geological site characterisation programmes typically rely on drill cores for direct information on subsurface rocks. However, porosity, transport properties and porewater composition measured on drill cores can deviate from in-situ values due to two main artefacts caused by drilling and sample recovery: (1) mechanical disruption that increases porosity and (2) contamination of the porewater by drilling fluid. We investigated the effect and magnitude of these perturbations on large drill core samples (12–20 cm long, 5 cmdiameter) of high-grade, granitic gneisses obtained from 350 to 600 m depth in a borehole on Olkiluoto Island (SW Finland). The drilling fluid was traced with sodium–iodide. By combining out-diffusion experiments, gravimetry, UV-microscopy and iodide mass balance calculations, we successfully quantified the magnitudes of the artefacts: 2–6% increase in porosity relative to the bulk connected porosity and 0.9 to 8.9 vol.% contamination by drilling fluid. The spatial distribution of the drilling-induced perturbations was revealed by numerical simulations of 2D diffusion matched to the experimental data. This showed that the rims of the samples have a mechanically disrupted zone 0.04 to 0.22 cm wide, characterised by faster transport properties compared to the undisturbed centre (1.8 to 7.7 times higher pore diffusion coefficient). Chemical contamination was shown to affect an even wider zone in all samples, ranging from 0.15 to 0.60 cm, inwhich iodide enrichmentwas up to 180 mg/kgwater, compared to 0.5 mg/kgwater in the uncontaminated centre. For all samples in the present case study, it turned out that the magnitude of the artefacts caused by drilling and sample recovery is so small that no correction is required for their effects. Therefore, the standard laboratory measurements of porosity, transport properties and porewater composition can be taken as valid in-situ estimates. However, it is clear that the magnitudes strongly depend on site- and drilling-specific factors and therefore our results cannot be transferred simply to other locations. We recommend the approach presented in this study as a route to obtain reliable values in future drilling campaigns aimed at characterising in-situ bedrock properties.

Provenance analysis using Raman spectroscopy of carbonaceous material: A case study in the Southern Alps of New Zealand

Detrital provenance analyses in orogenic settings, in which sediments are collected at the outlet of a catchment, have become an important tool to estimate how erosion varies in space and time. Here we present how Raman Spectroscopy on Carbonaceous Material (RSCM) can be used for provenance analysis. RSCM provides an estimate of the peak temperature (RSCM-T) experienced during metamorphism. We show that we can infer modern erosion patterns in a catchment by combining new measurements on detrital sands with previously acquired bedrock data. We focus on the Whataroa catchment in the Southern Alps of New Zealand and exploit the metamorphic gradient that runs parallel to the main drainage direction. To account for potential sampling biases, we also quantify abrasion properties using flume experiments and measure the total organic carbon content in the bedrock that produced the collected sands. Finally, we integrate these parameters into a mass-conservative model. Our results first demonstrate that RSCM-T can be used for detrital studies. Second, we find that spatial variations in tracer concentration and erosion have a first-order control on the RSCM-T distributions, even though our flume experiments reveal that weak lithologies produce substantially more fine particles than do more durable lithologies. This result implies that sand specimens are good proxies for mapping spatial variations in erosion when the bedrock concentration of the target mineral is quantified. The modeling suggests that highest present-day erosion rates (in Whataroa catchment) are not situated at the range front but around 10 km into the mountain belt.

Bedload composition, transport and modification in rivers of Westland, New Zealand, with implications for the distribution of alluvial pounamu (jade)

River bedload surveyed at 50 sites in Westland is dominated by Alpine Schist or Torlesse Greywacke from the Alpine Fault hanging wall, with subordinate Pounamu Ultramafics or footwall-derived Western Province rocks. Tumbling experiments found ultramafics to have the lowest attrition rates, compared with greywacke sandstone and granite (which abrade to produce silt to medium-sand), or incompetent schist (which fragments). Arahura has greater total concentrations (103–105 t/km2) and proportions (5–40%) of ultramafic bedload compared with Hokitika and Taramakau catchments (101–104 t/km2, mostly <10%), matching relative areas of mapped Pounamu Ultramafic bedrock, but enriched relative to absolute areal proportions. Western Province rocks downthrown by the Alpine Fault are under-represented in the bedload. Enriched concentrations of ultramafic bedload decrease rapidly with distance downstream from source rock outcrops, changing near prominent ice-limit moraines. Bedload evolution with transport involves both downstream fining and dilution from tributaries, in a sediment supply regime more strongly influenced by tectonics and the imprint of past glaciation. Treasured New Zealand pounamu (jade) is associated with ultramafic rocks. Chances of discovery vary between catchments, are increased near glacial moraines, and are highest near source-rock outcrops in remote mountain headwaters.

Long-bone fractures in persons with spinal cord injury

STUDY DESIGN Retrospective data analysis. OBJECTIVES To document fracture characteristics, management and related complications in individuals with traumatic spinal cord injury (SCI). SETTING Rehabilitation centre for SCI individuals. METHOD Patients' records were reviewed. Patients with traumatic SCI and extremity fractures that had occurred after SCI were included. Patient characteristics, fractured bone, fracture localisation, severity and management (operative/conservative), and fracture-related complications were extracted. RESULTS A total of 156 long-bone fractures in 107 SCI patients (34 women and 73 men) were identified. The majority of patients were paraplegics (77.6%) and classified as American Spinal Injury Association Impairment Scale A (86.0%). Only the lower extremities were affected, whereby the femur (60.9% of all fractures) was fractured more frequently than the lower leg (39.1%). A total of 70 patients (65.4%) had one fracture, whereas 37 patients (34.6%) had two or more fractures. Simple or extraarticular fractures were most common (75.0%). Overall, 130 (83.3%) fractures were managed operatively. Approximately half of the femur fractures (48.2%) were treated with locking compression plates. In the lower leg, fractures were mainly managed with external fixation (48.8%). Conservative fracture management was applied in 16.7% of the cases and consisted of braces or a well-padded soft cast. Fracture-associated complications were present in 13.5% of the cases but did not differ significantly between operative (13.1%) and conservative (15.4%) fracture management. CONCLUSION SCI was associated with simple or extraarticular fractures of the distal femur and the lower leg. Fractures were mainly managed operatively with a low complication rate.

Patterns of landscape form in the upper Rhône basin, Central Swiss Alps, predominantly show lithologic controls despite multiple glaciations and variations in rock uplift rates

The development of topography depends mainly on the interplay between uplift and erosion. These processes are controlled by various factors including climate, glaciers, lithology, seismic activity and short-term variables, such as anthropogenic impact. Many studies in orogens all over the world have shown how these controlling variables may affect the landscape's topography. In particular, it has been hypothesized that lithology exerts a dominant control on erosion rates and landscape morphology. However, clear demonstrations of this influence are rare and difficult to disentangle from the overprint of other signals such as climate or tectonics. In this study we focus on the upper Rhône Basin situated in the Central Swiss Alps in order to explore the relation between topography, possible controlling variables and lithology in particular. The Rhône Basin has been affected by spatially variable uplift, high orographically driven rainfalls and multiple glaciations. Furthermore, lithology and erodibility vary substantially within the basin. Thanks to high-resolution geological, climatic and topographic data, the Rhône Basin is a suitable laboratory to explore these complexities. Elevation, relief, slope and hypsometric data as well as river profile information from digital elevation models are used to characterize the landscape's topography of around 50 tributary basins. Additionally, uplift over different timescales, glacial inheritance, precipitation patterns and erodibility of the underlying bedrock are quantified for each basin. Results show that the chosen topographic and controlling variables vary remarkably between different tributary basins. We investigate the link between observed topographic differences and the possible controlling variables through statistical analyses. Variations of elevation, slope and relief seem to be linked to differences in long-term uplift rate, whereas elevation distributions (hypsometry) and river profile shapes may be related to glacial imprint. This confirms that the landscape of the Rhône Basin has been highly preconditioned by (past) uplift and glaciation. Linear discriminant analyses (LDAs), however, suggest a stronger link between observed topographic variations and differences in erodibility. We therefore conclude that despite evident glacial and tectonic conditioning, a lithologic control is still preserved and measurable in the landscape of the Rhône tributary basins.

Sediment mobilization deposits from episodic subsurface fluid flow—A new tool to reveal long-term earthquake records?

Subsurface fluid flow can be affected by earthquakes; increased spring activity, mud vol- cano eruptions, groundwater fluctuations, changes in geyser frequency, and other forms of altered subsurface fluid flow have been documented during, after, or even prior to seismic shaking. Recently discovered giant pockmarks on the bottom of Lake Neuchâtel, Switzerland, are the lake-floor expression of subsurface fluid flow. They discharge groundwater from the Jura Mountains karstic aquifers and experience episodically increased subsurface fluid flow documented by subsurface sediment mobilization deposits at the levees of the pockmarks. In this study, we present the spatio-temporal distribution of event deposits from these phases of sediment expulsion and of multiple time-correlative mass-transport deposits. We report five striking instances of concurrent multiple subsurface sediment deposits and multiple mass- transport deposits since late glacial times, for which we propose past earthquakes as a trigger. Comparison of this new event catalogue with historic earthquakes and other independent paleoseismic records suggests that initiation of sediment expulsion requires a minimum mac- roseismic intensity of VII. Thus, our study presents for the first time sedimentary deposits resulting from increased subsurface fluid flow as a paleoseismic proxy.

Regional surface morphology of comet 67P/Churyumov-Gerasimenko from Rosetta/OSIRIS images

Aims. The OSIRIS camera onboard the Rosetta spacecraft has been acquiring images of the comet 67P/Churyumov-Gerasimenko (67P)'s nucleus at spatial resolutions down to similar to 0.17 m/px ever since Aug. 2014. These images have yielded unprecedented insight into the morphological diversity of the comet's surface. This paper presents an overview of the regional morphology of comet 67P. Methods. We used the images that were acquired at orbits similar to 20-30 km from the center of the comet to distinguish different regions on the surface and introduce the basic regional nomenclature adopted by all papers in this Rosetta special feature that address the comet's morphology and surface processes. We used anaglyphs to detect subtle regional and topographical boundaries and images from close orbit (similar to 10 km from the comet's center) to investigate the fine texture of the surface. Results. Nineteen regions have currently been defined on the nucleus based on morphological and/or structural boundaries, and they can be grouped into distinctive region types. Consolidated, fractured regions are the most common region type. Some of these regions enclose smooth units that appear to settle in gravitational sinks or topographically low areas. Both comet lobes have a significant portion of their surface covered by a dusty coating that appears to be recently placed and shows signs of mobilization by aeolian-like processes. The dusty coatings cover most of the regions on the surface but are notably absent from a couple of irregular large depressions that show sharp contacts with their surroundings and talus-like deposits in their interiors, which suggests that short-term explosive activity may play a significant role in shaping the comet's surface in addition to long-term sublimation loss. Finally, the presence of layered brittle units showing signs of mechanical failure predominantly in one of the comet's lobes can indicate a compositional heterogeneity between the two lobes.

Influence of catchment quality and altitude on the water and sediment composition of 68 small lakes in Central Europe

68 lakes (63 Swiss, 2 French and 3 Italian) located in an altitudinal range between 334 and 2339m spanning a wide range of land-use have been investigated. The aim of the study was to discuss influences of geographic location, vegetation and land-use in the catchment area on the water and sediment chemistry of small lakes. Detailed quantitative description of land-use, vegetation, and climate in the watershed of all lakes was established. Surface and bottom water samples collected from each lake were analyzed for major ions and nutrients. Correlations were interpreted using linear regression analysis. Chemical parameters of water and sediment reflect the characteristics of the catchment areas. All lakes were alkaline since they were situated on calcareous bedrock. Concentrations of nitrogen and phosphorus strongly increase with increasing agricultural land-use. Na and K, however, are positively correlated with the amount of urbanization within the catchment area. These elements as well as dissolved organic carbon (DOC), Mg, Ca, and alkalinity, increase when the catchment is urbanized or used for agriculture. Total nitrogen and organic carbon in the sediments decrease distinctly if large parts of the catchment consist of bare land. No correlations between sediment composition and maximum water depth or altitude of the lakes were found.¶Striking differences in the water compositions of lakes above and below approximately 700 m of altitude were observed. Concentrations of total nitrogen and nitrate, total phosphorus, DOC, Na, K, Mg, Ca, and alkalinity are distinctly higher in most lakes below 700 m than above, and the pH of the bottom waters of these lakes is generally lower. Estimates of total nitrogen concentrations, even in remote areas, indicate that precipitation is responsible for increased background concentrations. At lower altitudes nitrogen concentrations in lakes is explained by the nitrogen loaded rain from urban areas deposited on the catchment, and with high percentages of agricultural land-use in the watershed.

Modification of aviary design reduces incidence of falls, collisions and keel bone damage in laying hens

Non-cage housing systems for laying hens such as aviaries provide greater freedom to perform species-specific behavior and thus are thought to improve welfare of the birds; however, aviaries are associated with a high prevalence of keel bone damage (fractures and deviations), which is a major welfare problem in commercial laying hens. Potential causes of keel bone damage are falls and collisions with internal housing structures that occur as birds move between tiers or perches in the aviary. The aim of this study was to investigate the scope for reducing keel bone damage by reducing falls and collisions through modifications of aviary design. Birds were kept in 20 pens in a laying hen house (225 hens per pen) that were assigned to four different treatments (n = 5 pens per treatment group) including (1) control pens and pens modified by the addition of (2) perches, (3) platforms and (4) ramps. Video recordings at 19, 22, 29, 36 and 43 weeks of age were used to analyze controlled movements and falls (including details on occurrence of collision, cause of fall, height of fall and behavior after fall) during the transitional dusk and subsequent dark phase. Palpation assessments (focusing on fractures and deviations) using 20 focal hens per pen were conducted at 18, 20, 23, 30, 37, 44, 52 and 60 weeks of age. In comparison to the control group, we found 44% more controlled movements in the ramp (P = 0.003) and 47% more controlled movements in the platform treatments (P = 0.014) as well as 45% fewer falls (P = 0.006) and 59% fewer collisions (P < 0.001) in the ramp treatment. There were no significant differences between the control and perch treatments. Also, at 60 weeks of age, 23% fewer fractured keel bones were found in the ramp compared with the control treatment (P = 0.0053). After slaughter at 66 weeks of age, no difference in keel bone damage was found between treatment groups and the prevalence of fractures increased to an average of 86%. As a potential mechanism to explain the differences in locomotion, we suggest that ramps facilitated movement in the vertical plane by providing a continuous path between the tiers and thus supported more natural behavior (i.e. walking and running) of the birds. As a consequence of reducing events that potentially damage keel bones, the installation of ramps may have reduced the prevalence of keel fractures for a major portion of the flock cycle. We conclude that aviary design and installation of specific internal housing structures (i.e. ramps and platforms) have considerable potential to reduce keel bone damage of laying hens in aviary systems.

Isolation and Characterization of a Microorganism from Groundwater that Reduces Arsenate

This is an investigation into the microbially mediated processes involved in the transformation of arsenic. With the recent change in the Federal Maximum Contaminant Level for arsenic in drinking water, an increasing amount of resources are being devoted to understanding the mechanisms involved in the movement of arsenic. Arsenic in drinking water typically comes from natural sources, but the triggers that result in increased release of arsenic from parent material are poorly understood. Knowledge of these processes is necessary in order to make sound engineering decisions regarding drinking water management practices. Recent years have brought forth the idea that bacteria play a significant role in arsenic cycling. Groundwater is a major source of potable water in this and many other countries. To date, no reports have been made indicating the presence and activity of arsenate reducing bacteria in groundwater settings, which may increase dissolved arsenic concentrations. This research was designed to address this question and has shown that these bacteria are present in Maine groundwater. Two Maine wells were sampled in order to culture resident bacteria that are capable of dissimilatory arsenate reduction. Samples were collected using anaerobic techniques fiom wells in Northport and Green Lake. These samples were amended with specific compounds to enrich the resident population of arsenate utilizing bacteria. These cultures were monitored over time to establish rates of arsenate reduction. Cultures fiom both sites exhibited arsenate reduction in initial enrichment cultures. Isolates obtained fiom the Green Lake enrichments, however, did not reduce arsenate. This indicates either that a symbiotic relationship was required for the observed arsenate reduction or that fast-growing fermentative organisms that could survive in high arsenate media were picked in the isolation procedure. The Northport cultures exhibited continued arsenate reduction after isolation and successive transfers into fiesh media. The cultured bacteria reduced the majority of 1 a arsenate solutions in less than one week, accompanied by a corresponding oxidation of lactate. The 16s rRNA fiom the isolate was arnplifled and sequenced. The results of the DNA sequence analysis indicate that the rRNA sequence of the bacteria isolated at the Northport site is unique. This means that this strain of bacteria has not been reported before. It is in the same taxonomic subgroup as two previously described arsenate respirers. The implications of this study are significant. The fact that resident bacteria are capable of reducing arsenate has implications for water management practices. Reduction of arsenate to arsenite increases the mobility of the compound, as well as the toxicity. An understanding of the activity of these types of organisms is necessary in order to understand the contribution they are making to arsenic concentrations in drinking water. The next step in this work would be to quantitj the actual loading of dissolved arsenic present in aquifers because of these organisms.

AN ASSESSMENT OF THE RELATION OF CANCER MORTALITY RATES AND SOURCES OF DRINKING WATER IN TEXAS

The geographic distribution of average annual age-adjusted mortality rates (1964-1976) for four types of cancer (all cancer sites combined, gastrointestinal, urinary, and lung cancer) were compared by sources of drinking water for 254 Texas counties and county rural areas and 301 Texas cities. Exposure variables considered were surface versus ground water, public water supplies versus individuals wells, and trihalomethane levels in municipal water supplies. Each general source of "surface" and "ground" water was further divided by aggregating ground water using areas by aquifers and surface water using study areas by river basins. Potential confounding variables taken into account included median education, employment in cancer risk industries, population mobility, ethnicity, and urbanicity. A pattern of higher and lower cancer mortality rates was found for populations using some aquifers and river basins. Further study is required to determine whether the differences in cancer mortality rates that were observed are related to drinking water content or are coincidental with differences in personal characteristics which could not be taken into account in this ecologic study design. ^

Alpine permafrost index map

The objective of this study is the production of an Alpine Permafrost Index Map (APIM) covering the entire European Alps. A unified statistical model that is based on Alpine-wide permafrost observations is used for debris and bedrock surfaces across the entire Alps. The explanatory variables of the model are mean annual air temperatures, potential incoming solar radiation and precipitation. Offset terms were applied to make model predictions for topographic and geomorphic conditions that differ from the terrain features used for model fitting. These offsets are based on literature review and involve some degree of subjective choice during model building. The assessment of the APIM is challenging because limited independent test data are available for comparison and these observations represent point information in a spatially highly variable topography. The APIM provides an index that describes the spatial distribution of permafrost and comes together with an interpretation key that helps to assess map uncertainties and to relate map contents to their actual expression in terrain. The map can be used as a first resource to estimate permafrost conditions at any given location in the European Alps in a variety of contexts such as research and spatial planning. Results show that Switzerland likely is the country with the largest permafrost area in the Alps, followed by Italy, Austria, France and Germany. Slovenia and Liechtenstein may have marginal permafrost areas. In all countries the permafrost area is expected to be larger than the glacier-covered area.