Late Ordovician climbing‐dune cross‐stratification: a signature of outburst floods in proglacial outwash environments?

Climbing dune‐scale cross‐statification is described from Late Ordovician paraglacial successions of the Murzuq Basin (SW Libya). This depositional facies is comprised of medium‐grained to coarse‐grained sandstones that typically involve 0·3 to 1 m high, 3 to 5 m in wavelength, asymmetrical laminations. Most often stoss‐depositional structures have been generated, with preservation of the topographies of formative bedforms. Climbing‐dune cross‐stratification related to the migration of lower‐flow regime dune trains is thus identified. Related architecture and facies sequences are described from two case studies: (i) erosion‐based sandstone sheets; and (ii) a deeply incised channel. The former characterized the distal outwash plain and the fluvial/subaqueous transition of related deltaic wedges, while the latter formed in an ice‐proximal segment of the outwash plain. In erosion‐based sand sheets, climbing‐dune cross‐stratification results from unconfined mouth‐bar deposition related to expanding, sediment‐laden flows entering a water body. Within incised channels, climbing‐dune cross‐stratification formed over eddy‐related side bars reflecting deposition under recirculating flow conditions generated at channel bends. Associated facies sequences record glacier outburst floods that occurred during early stages of deglaciation and were temporally and spatially linked with subglacial drainage events involving tunnel valleys. The primary control on the formation of climbing‐dune cross‐stratification is a combination between high‐magnitude flows and sediment supply limitations, which lead to the generation of sediment‐charged stream flows characterized by a significant, relatively coarse‐grained, sand‐sized suspension‐load concentration, with a virtual absence of very coarse to gravelly bedload. The high rate of coarse‐grained sand fallout in sediment‐laden flows following flow expansion throughout mouth bars or in eddy‐related side bars resulted in high rates of transfer of sands from suspension to the bed, net deposition on bedform stoss‐sides and generation of widespread climbing‐dune cross‐stratification. The later structure has no equivalent in the glacial record, either in the ancient or in the Quaternary literature, but analogues are recognized in some flood‐dominated depositional systems of foreland basins.     

Clay mineralogy as a tool for integrated sequence stratigraphic and palaeogeographic reconstructions: Late Oligocene–Early Aquitanian western internal South Iberian Margin,Spain

The different hydrodynamic behaviour of detrital clay minerals in the marine depositional environment allows assessment of relative sea‐level variations in the sedimentary record. Comparison of smectite and kaolinite:illite (S+K:I) changes with the global sea‐level curves and with the third‐order cycles of the eustatic curve for European basins allows assessment of the influence of global eustasy and local tectonics on sequence stratigraphy. In the South Iberian Margin, sedimentation took place both on open‐marine platforms and in deeper water areas. On this margin during the Late Oligocene to Early Aquitanian, the variations in sedimentation were caused not only by global eustasy but also by compressive tectonics. Correlations were made between the S+K:I cycles and the third‐order cycles for European basins, enabling the definition of four third‐order sedimentary sequences (here called C1, C2, C3, and A1) and two lower‐order sequences within C3 and A1 (here called C3a, C3b, A1a, and A1b) related to tectonic movements. High S+K:I values were observed during episodes of maximum flooding in each sequence and lower‐order sequence and in each succession, enabling changes in palaeocoastal morphology to be considered. Copyright © 2012 John Wiley & Sons, Ltd.     

Spatially distributed tracer‐aided modelling to explore water and isotope transport,storage and mixing in a pristine,humid tropical catchment

Rapidly transforming headwater catchments in the humid tropics provide important resources for drinking water, irrigation, hydropower, and ecosystem connectivity. However, such resources for downstream use remain unstudied. To improve understanding of the behaviour and influence of pristine rainforests on water and tracer fluxes, we adapted the relatively parsimonious, spatially distributed tracer‐aided rainfall–runoff (STARR) model using event‐based stable isotope data for the 3.2‐km² San Lorencito catchment in Costa Rica. STARR was used to simulate rainforest interception of water and stable isotopes, which showed a significant isotopic enrichment in throughfall compared with gross rainfall. Acceptable concurrent simulations of discharge (Kling–Gupta efficiency [KGE] ~0.8) and stable isotopes in stream water (KGE ~0.6) at high spatial (10 m) and temporal (hourly) resolution indicated a rapidly responding system. Around 90% of average annual streamflow (2,099 mm) was composed of quick, near‐surface runoff components, whereas only ~10% originated from groundwater in deeper layers. Simulated actual evapotranspiration (ET) from interception and soil storage were low (~420 mm/year) due to high relative humidity (average 96%) and cloud cover limiting radiation inputs. Modelling suggested a highly variable groundwater storage (~10 to 500 mm) in this steep, fractured volcanic catchment that sustains dry season baseflows. This groundwater is concentrated in riparian areas as an alluvial–colluvial aquifer connected to the stream. This was supported by rainfall–runoff isotope simulations, showing a “flashy” stream response to rainfall with only a moderate damping effect and a constant isotope signature from deeper groundwater (~400‐mm additional mixing volume) during baseflow. The work serves as a first attempt to apply a spatially distributed tracer‐aided model to a tropical rainforest environment exploring the hydrological functioning of a steep, fractured‐volcanic catchment. We also highlight limitations and propose a roadmap for future data collection and spatially distributed tracer‐aided model development in tropical headwater catchments.     

SHRIMP U–Pb ages of detrital zircons from the Early Cretaceous Nakdong Formation,South East Korea: Timing of initiation of the Gyeongsang Basin and its provenance

To constrain the depositional age of the lowermost Nakdong Formation in the Early Cretaceous Gyeongsang Basin, SHRIMP U–Pb age determination was carried out on zircon separates. The U–Pb compositions of detrital zircons from the Nakdong Formation yield a wide range of ages from the Archean to the Cretaceous but show a marked contrast in age distribution according to the geographical locations within the basin. The provenance of the southern Nakdong Formation is dominantly the surrounding Yeongnam Massif, which is composed of Paleoproterozoic metamorphic rocks and Triassic to Jurassic plutonic rocks, whereas the central to northern Nakdong Formation records significant sediment derivation from the Okcheon Metamorphic Belt, which is distributed to the northwest, in addition to the contribution from the Yeongnam Massif. It is suggested that the maximum depositional age of the Nakdong Formation is ca 127 Ma, based on its youngest detrital zircon age population. The onset of its deposition at 127 Ma coincided with the tectonic inversion in East Asia from a compressional to an extensional geodynamic setting, probably due to the contemporaneous change in the drift direction of the Izanagi Plate and its subsequent oblique subduction.     

Application of Hilbert‐like transforms for enhanced processing of full tensor magnetic gradient data

Commonly, geomagnetic prospection is performed via scalar magnetometers that measure values of the total magnetic intensity. Recent developments of superconducting quantum interference devices have led to their integration in full tensor magnetic gradiometry systems consisting of planar‐type first‐order gradiometers and magnetometers fabricated in thin‐film technology. With these systems measuring directly the magnetic gradient tensor and field vector, a significantly higher magnetic and spatial resolution of the magnetic maps is yield than those produced via conventional magnetometers. In order to preserve the high data quality in this work, we develop a workflow containing all the necessary steps for generating the gradient tensor and field vector quantities from the raw measurement data up to their integration into high­resolution, low­noise, and artefactless two‐dimensional maps of the magnetic field vector. The gradient tensor components are processed by superposition of the balanced gradiometer signals and rotation into an Earth‐centred Earth‐fixed coordinate frame. As the magnetometers have sensitivity lower than that of gradiometers and the total magnetic intensity is not directly recorded, we employ Hilbert‐like transforms, e.g., integration of the gradient tensor components or the conversion of the total magnetic intensity derived by calibrated magnetometer readings to obtain these values. This can lead to a better interpretation of the measured magnetic anomalies of the Earth's magnetic field that is possible from scalar total magnetic intensity measurements. Our conclusions are drawn from the application of these algorithms on a survey acquired in South Africa containing full tensor magnetic gradiometry data.     

An analytical solution for describing the transient temperature distribution in an aquifer thermal energy storage system

A mathematical model is developed for predicting the temperature distribution in an aquifer thermal energy storage (ATES) system, which consists of a confined aquifer bounded from above and below by the rocks of different geological properties. The main transfer processes of heat include the conduction and advection in the aquifer and the conduction in the rocks. The semi‐analytical solution in dimensionless form for the model is developed by Laplace transforms and its corresponding time‐domain solution is evaluated by the modified Crump method. Field geothermal property data are used to simulate the temperature distribution in an ATES system. The results show that the heat transfer in the aquifer is fast and has a vast effect on the vicinity of the wellbore. However, the aquifer temperature decreases with increasing radial and vertical distances. The temperature in the aquifer may be overestimated when ignoring the effect of thermal conductivity. The temperature distribution in an ATES system depends on the vertical thermal conduction in the rocks and the horizontal advection and thermal conduction in the aquifer. The present solution is useful in designing and simulating the heat injection facility in the ATES systems. Copyright © 2010 John Wiley & Sons, Ltd.     

Effect of the Dry Ashing Method on Cadmium Isotope Measurements in Soil and Plant Samples

A dry ashing method is commonly used to remove organic material from samples prior to geochemical analysis. In the course of this study, the Cd isotope ratios of a series of soil and plant reference materials and samples were studied to evaluate the effect of the dry ashing method on measurement results of Cd isotope ratios. The samples were pre‐treated using the dry ashing method and high‐pressure bomb for comparison. The results show that the digestion using high‐pressure bombs did not lead to Cd loss, but using the dry ashing method would cause different proportions of Cd loss. The whole range of Cd isotope difference between two methods was from ?0.07‰ to 3.01‰. There was also an obvious difference in measured Cd isotope ratios from the same leaf sample pre‐treated independently by the dry ashing method, indicating that the amount of Cd loss and the effect on Cd isotope measurement during dry ashing is related to the properties of the samples. Therefore, dry ashing may not be appropriate for the removal of organic material in Cd isotope ratio measurement, especially for samples with high organic contents. The δ^114/110Cd values of reference materials NIST SRM 1573a and GSD‐30 are reported for the first time in this study.     

Parallelization of the distinct lattice spring model

The distinct lattice spring model (DLSM) is a newly developed numerical tool for modeling rock dynamics problems, i.e. dynamic failure and wave propagation. In this paper, parallelization of DLSM is presented. With the development of parallel computing technologies in both hardware and software, parallelization of a code is becoming easier than before. There are many available choices now. In this paper, Open Multi‐Processing (OpenMP) with multicore personal computer (PC) and message passing interface (MPI) with cluster are selected as the environments to parallelize DLSM. Performances of these parallel DLSM codes are tested on different computers. It is found that the parallel DLSM code with OpenMP can reach a maximum speed‐up of 4.68× on a quad‐core PC. The parallel DLSM code with MPI can achieve a speed‐up of 40.886× when 256 CPUs are used on a cluster. At the end of this paper, a high‐resolution model with four million particles, which is too big to handle by the serial code, is simulated by using the parallel DLSM code on a cluster. It is concluded that the parallelization of DLSM is successful. Copyright © 2011 John Wiley & Sons, Ltd.     

Evidence of ameliorated Middle Weichselian climate and sub‐arctic environment in the western Baltic region: coring lake sediments at Klintholm,Møn,Denmark

Coring through glaciotectonically stacked Quaternary sediments situated below sea level on the island of Møn, Denmark, recovered a succession of interstadial sediments of Middle Weichselian age. Plant and animal remains including insects found in laminated sand and mud indicate deposition in a lake surrounded by dwarf shrubs, herbs, mosses and rare trees. The insect fauna indicates a mean July temperature of 8–12 °C, suggesting an arctic to sub‐arctic environment, while winter temperatures around ?8 to ?22 °C suggest periglacial conditions with permafrost. Luminescence dating of sediment samples gave ages from 48–29 ka, and radiocarbon dating indicates deposition of plant fragments between 45 and 36 ka BP. The fossil assemblage from Møn shows close resemblance to those from other sites with similar ages found in the vicinity of the western Baltic Basin.     

Process length variation in cysts of the dinoflagellate Protoceratium reticulatum,from surface sediments of the Baltic–Kattegat–Skagerrak estuarine system: a regional salinity proxy

Mertens, K. N., Dale, B., Ellegaard, M., Jansson, I.‐M., Godhe, A., Kremp, A. & Louwye, S. 2010: Process length variation in cysts of the dinoflagellate Protoceratium reticulatum, from surface sediments of the Baltic–Kattegat–Skagerrak estuarine system: a regional salinity proxy. Boreas, 10.1111/j.1502‐3885.2010.00193.x. ISSN 0300‐9483. Results are presented from a regional comparison of average process length variation in cysts of Protoceratium reticulatum and Lingulodinium polyedrum, extracted from surface sediments in the Skagerrak–Kattegat–Baltic estuarine system, with the environmental variables of seawater temperature and salinity. Although too few cysts of Lingulodinium polyedrum were recovered from the sediments to make reliable correlations, cysts of Protoceratium reticulatum were well represented, and average process length was correlated significantly with both salinity and temperature. Owing to dominant summer surface production, and regional covariation between salinity and density, we propose the use of the significant correlation with summer sea surface salinity (SSS_summer) by the equation SSS_summer=3.16 × average process length ?0.84 (R²=0.8). Application of this equation down‐core in Limfjord (northern Denmark) shows its usefulness as a regional palaeosalinity proxy.