Simulating Oxygen Intrusion into Highly Heterogeneous Fractured Media Using High Performance Computing

Fractured crystalline rocks are under consideration by several countries as host formations for high-level nuclear waste repositories. The redox evolution in these host rock formations is an important issue for the stability and safety of these disposal sites. If, for example, during a glaciation/deglaciation event, oxygen-rich glacial meltwater penetrates to the depth of the planned repository, some of the engineered barriers would be adversely affected. Moreover, oxidizing conditions would increase the solubility and mobility of many radionuclides. Reactive transport simulations, which are typically used to assess the redox buffering capacity of these host rock formations, are computationally demanding, and thus, calculations for the evaluation of oxygen penetration are usually carried out over simplified geometries and the heterogeneity of the site, both physical (e.g., variability in the groundwater flow field and the kinematic porosity) and mineralogical (e.g., variability in the abundance of Fe(II)-bearing minerals), is usually represented in a simplified fashion. Here, it is shown how a recently developed numerical framework, combined with high performance computing technologies, allows the full geometrical, physical and mineralogical complexity of the site under study to be efficiently included in these types of reactive transport calculations. A synthetically generated realistic three-dimensional fractured medium is used to assess oxygen penetration patterns and their dependence on both the hydrogeological conditions and the availability of Fe(II)-bearing minerals. The results of the calculations point out the significant influence of both physical and mineralogical heterogeneity on the oxygen penetration patterns, thus highlighting the importance of a model parameterisation consistent with the site complexity.     

Identification of the Fugloyarbanki tephra in the NGRIP ice core: a key tie‐point for marine and ice‐core sequences during the last glacial period

A visible tephra horizon in the NGRIP ice core has been identified by geochemical analysis as the Fugloyarbanki Tephra, a widespread marker horizon in marine cores from the Faroe Islands area and the northern North Atlantic. An age of 26 740 ± 390 yr b2k (1σ uncertainty) is derived for this tephra according to the new Greenland Ice Core Chronology (GICC05) based on multi‐parameter counting of annual layers. Detection of this tephra for the first time within the NGRIP ice core provides a key tie‐point between marine and ice‐core records during the transition between MIS 3 and 2. Identification of this volcanic event within the Greenland records demonstrates the future potential of using tephrochronology to precisely correlate palaeoarchives in widely separated localities that span the last glacial period, as well as providing a potential method for examining the extent of the radiocarbon marine reservoir effect at this time. Copyright © 2008 John Wiley & Sons, Ltd.     

Reconstructing the Younger Dryas ice dammed lake in the Baltic Basin: Bathymetry, area and volume

A digital 3D-reconstruction of the Baltic Ice Lake's (BIL) configuration during the termination of the Younger Dryas cold phase (ca. 11 700 cal. yr BP) was compiled using a combined bathymetric–topographic Digital Terrain Model (DTM), Scandinavian ice sheet limits, Baltic Sea Holocene bottom sediment thickness information, and a paleoshoreline database maintained at the Lund University. The bathymetric–topographic DTM, assembled from publicly available data sets, has a resolution of 500 × 500 m on Lambert Azimuthal Equal Area projection allowing area and volume calculations of the BIL to be made with an unprecedented accuracy. When the damming Scandinavian ice sheet margin eventually retreated north of Mount Billingen, the high point in terrain of Southern central Sweden bordering to lower terrain further to the north, the BIL was catastrophically drained resulting in a 25 m drop of the lake level. With our digital reconstruction, we estimate that approximately 7800 km³ of water drained during this event and that the ice dammed lake area was reduced by ca. 18%. Building on previous results suggesting drainage over 1 to 2 years, our lake volume calculations imply that the freshwater flux to the contemporaneous sea in the west was between about 0.12 and 0.25 Sv. The BIL reconstruction provides new detailed information on the paleogeography in the area of southern Scandinavia, both before and after the drainage event, with implications for interpretations of geological records concerning the post-glacial environmental development.     

The Third GABLS Intercomparison Case for Evaluation Studies of Boundary-Layer Models. Part B: Results and Process Understanding

We describe and analyze the results of the third global energy and water cycle experiment atmospheric boundary layer Study intercomparison and evaluation study for single-column models. Each of the nineteen participating models was operated with its own physics package, including land-surface, radiation and turbulent mixing schemes, for a full diurnal cycle selected from the Cabauw observatory archive. By carefully prescribing the temporal evolution of the forcings on the vertical column, the models could be evaluated against observations. We focus on the gross features of the stable boundary layer (SBL), such as the onset of evening momentum decoupling, the 2-m minimum temperature, the evolution of the inertial oscillation and the morning transition. New process diagrams are introduced to interpret the variety of model results and the relative importance of processes in the SBL; the diagrams include the results of a number of sensitivity runs performed with one of the models. The models are characterized in terms of thermal coupling to the soil, longwave radiation and turbulent mixing. It is shown that differences in longwave radiation schemes among the models have only a small effect on the simulations; however, there are significant variations in downward radiation due to different boundary-layer profiles of temperature and humidity. The differences in modelled thermal coupling to the land surface are large and explain most of the variations in 2-m air temperature and longwave incoming radiation among models. Models with strong turbulent mixing overestimate the boundary-layer height, underestimate the wind speed at 200 m, and give a relatively large downward sensible heat flux. The result is that 2-m air temperature is relatively insensitive to turbulent mixing intensity. Evening transition times spread 1.5 h around the observed time of transition, with later transitions for models with coarse resolution. Time of onset in the morning transition spreads 2 h around the observed transition time. With this case, the morning transition appeared to be difficult to study, no relation could be found between the studied processes, and the variation in the time of the morning transition among the models.     

Influence of algal farming on fish assemblages

We examined the influence of algal farming on fish assemblages in two shallow coastal lagoons in Zanzibar, Tanzania. Fish assemblages were visually investigated using a belt transect method and the line-intercept technique was used to examine the substrate composition. 101 species of fish belonging to 31 families were recorded. Algal farming affected the associated fish fauna in terms of abundance, species richness, trophic identity and fish community composition. However, the impact differed between the lagoons. Algal farms in one lagoon hosted a more abundant and diversified fish fauna than controls, whereas farms in the other lagoon exhibited lower fish densities and similar species diversity compared to controls. The discrepancies between lagoons may be an effect of differences in farming intensity and character of the substratum.     

Comparison of DNA damage in the early life stages of cod, Gadus morhua, originating from the Barents Sea and Baltic Sea

DNA adducts in cod embryos and larvae were analysed by ³²P-postlabeling to test the hypothesis that anthropogenic substances, which could form reactive intermediates, are involved in the reproductive failure of cod (Gadus morhua) from the Baltic Sea. A comparison with cod from the Barents Sea was performed. The mean value of DNA adducts in cod embryos/larvae from the Baltic Sea was 2.3 nmol of adducts/mol nucleotides, compared to 0.12 nmol of adducts/mol nucleotides in the embryos/larvae from the Barents Sea.     

Taking varves to bits: Scanning electron microscopy in the study of laminated sediments and varves

Conventional high resolution studies of varved sediments are able to identify clastic and biogenic laminae, but are often unable to resolve the nature of fine-scale lamination contained therein. This intra-annual signal provides us with the highest potential resolution from the sedimentary record and can be resolved using scanning electron microscopy (SEM). Six case studies from lacustrine and marine settings are presented to illustrate the combination of clastic and biogenic fabric types typically found in laminated sediments. Clastic laminae fabrics include those which originate through grain settling and those which are event deposits. The correct identification of event deposits is essential if varves are to be used chronologically. SEM-based biogenic laminae fabric studies have identified seasonal faunal successions where individual laminae may be less than 100 thick and most recently, deep chlorophyll maxima (DCM) summer diatom floras, providing an insight into seasonal scale processes. High resolution lamina fabric studies can provide a basis for generating records of seasonal and inter-annual variability, thus contributing to our understanding of lacustrine and marine processes and palaeoenvironmental interpretation.