Improved regional groundwater flow modeling using drainage features: a case study of the central northern karst aquifer system of Puerto Rico (USA)

In northern Puerto Rico (USA), subsurface conduit networks with unknown characteristics, and surface features such as springs, rivers, lagoons and wetlands, drain the coastal karst aquifers. In this study, drain lines connecting sinkholes and springs are used to improve the developed regional model by simulating the drainage effects of conduit networks. Implemented in an equivalent porous media (EPM) approach, the model with drains is able to roughly reproduce the spring discharge hydrographs in response to rainfall. Hydraulic conductivities are found to be scale dependent and significantly increase with higher test radius, indicating scale dependency of the EPM approach. Similar to other karst regions in the world, hydraulic gradients are steeper where the transmissivity is lower approaching the coastline. This study enhances current understanding of the complex flow patterns in karst aquifers and suggests that using a drainage feature improves modeling results where available data on conduit characteristics are minimal.     

Benthic macrofauna and productivity regimes in the Barents Sea — Ecological implications in a changing Arctic

Benthic faunal assemblages were analysed from 47 stations in the central and southern parts of the Barents Sea, together with sedimentary and water column parameters, daily ice records and modelled integrated primary productivity. Sampling spanned areas influenced by Atlantic Water (AW) to those lying under Arctic Water (ArW), and included stations with mixed water masses. Ice cover suppressed water column productivity in the northern areas. Three main faunal groups were identified, based on similarity of numerical faunal composition. The northern and southern faunal groups were separated by the northernmost penetration of AW in the bottom water and the third group, the Hopen group, was influenced by modified bank water. Faunal abundances were significantly higher within the southern faunal group relative to the northern group, but the numbers of taxa present were similar. The particularly rich fauna of the Hopen group reflected sediment heterogeneity and tight pelagic–benthic coupling. These results suggest that a retreat and thinning of the ice cover in the Barents Sea likely will result in the northern parts of the Barents Sea becoming more Atlantic in character, with a higher productivity at the sea floor.     

Detection of DNA damage in mussels and sea urchins exposed to crude oil using comet assay

The single-cell microgel electrophoresis assay or the comet assay was used to evaluate DNA damage of dispersed crude oil on sea urchins (Strongylocentrotus droebachiensis) and mussels (Mytilus edulis L.). Sea urchins were exposed to 0.06 and 0.25 mg/L dispersed crude oil in a continuous flow system, while the mussels were exposed to 0.015, 0.06 and 0.25 mg/L dispersed crude oil. Sea urchin coelomocytes and mussel haemocytes were sampled after 4 and 5 weeks exposure, respectively. In the sea urchin coelomocytes, there was a significant concentration-related increase in the percentage of DNA in comet tail. In mussel haemocytes, there was a significantly higher percentage of DNA in comet tail for all treatments compared to the control. The responses were concentration-related up to 0.06 mg/L oil. The two highest exposure concentrations of mussels were not significantly different from each other. These results indicate that the comet assay can be used for biomonitoring of DNA damage in marine invertebrates following oil contamination.     

Experimental study of the carbonation of partially serpentinized and weathered peridotites

Carbonation of partially serpentinized and weathered peridotites was studied experimentally under hydrothermal conditions (T: 200 °C, P_CO2: 130-180 bars). Experiments were performed in a closed system using whole rock drill core samples (height: 1 cm, diameter: 1 cm) as starting material. The initial samples were composed mainly of meshwork serpentine, relicts of primary olivine and an olivine weathering product (deweylite assemblage). Two types of solutions, each with a total salt content corresponding to that of average seawater (35 g/L dissolved salts), were used: (1) a Na-Ca-Cl solution (12.5 g/L CaCl₂ + 22.5 g/L NaCl) and (2) a NaCl solution (35 g/L NaCl). After 15-25 days of experimental treatment, the samples were partly covered with carbonates. In addition, noticeable carbonation reactions had occurred below the sample surfaces within zones with thicknesses up to 250 μm. In the Na-Ca-Cl solution, both the olivine relicts and the deweylite assemblage were partly replaced by calcite along the surrounding serpentine veins. However, the extent of calcitization was found to be considerably larger for the deweylite assemblage than for the olivine. Bulk fluid analyses show an increase in the Mg and Si concentrations with reaction time. In the NaCl solution, the deweylite assemblage was partly dissolved resulting in large voids within the reaction zone. In contrast, the olivine was replaced by magnesite. Under the conditions of our experiments, the meshwork serpentine was not reactive, but aided fluid infiltration into the rock samples. The experimentally produced microtextures closely resemble those found in natural examples. Our study elucidates the mechanisms by which carbonates form in ultramafic rocks under relatively high P_CO2-T conditions and particularly in the presence of Ca-bearing aqueous solutions. The existence of a serpentine meshtexture and the presence of weathering products formed from primary Mg-silicates may have significant beneficial effects on in situ CO₂ mineral sequestration in ultramafic rocks.     

Interpretation of single-well tracer tests using fractional-flow dimensions. Part 1: Theory and mathematical models

Generalized equations using fractional-flow dimensions were derived to estimate the Darcy and seepage velocities obtained from the point-dilution and the single-well injection-withdrawal field tests conducted in fractured-rock aquifers. Seepage velocities can only be estimated from single-well tests if the hydraulic conductivity and the hydraulic gradient are known a priori. However, if a radial-convergent test is also performed between two boreholes, the kinematic porosity can be estimated and be used to estimate the seepage velocity from the single-well test results. To apply the generalized equations, the flow dimension and the extent of the flow region must be known. Therefore, the generalized radial flow (GRF) model of Barker (1988; a generalized radial flow model for hydraulic tests in fractured rock. Water Resour Res 24(10):1796–1804) is used to estimate the flow dimension because of its wide range of applications. A pumping test performed on the boreholes will yield an estimate of the fractional-flow dimension by applying the GRF model. Electronic Publication     

DISCHARGE AND SEDIMENT-LOAD CHARACTERISTICS OF THE HILDA ROCK-GLACIER STREAM,CANADIAN ROCKY MOUNTAINS,ALBERTA

Water discharge and sediment load characteristics of the Hilda rock-glacier stream are described for portions of the 1985 and 1986 ablation seasons. Covering an area of 1.5 km², the Hilda rock glacier yields a relatively consistent flow, generally low in suspended sediments and relatively high in total dissolved solids during the ablation season. The seasonal hydrograph shows decreasing discharge and suspended sediment load and increasing dissolved loads through the ablation season as snowmelt water sources are depleted. Superimposed on this is a diurnal hydrograph showing that daily peaks in discharge lag several hours behind daily temperature maxima. Rain events produce hydrograph anomalies and act to flush sediment, elevating suspended loads, while diluting dissolved loads. These data suggest that rock-glacier hydrological systems are more responsive to short-term and transient meteorological conditions than previously suggested in the literature. Nonetheless, the data indicate that the Hilda rock glacier is more conservative in both water and sediment yield than glacier systems of comparable size and in the same environment.     

Building damage characteristics based on surveyed data and fragility curves of the 2011 Great East Japan tsunami

A large amount of buildings was damaged or destroyed by the 2011 Great East Japan tsunami. Numerous field surveys were conducted in order to collect the tsunami inundation extents and building damage data in the affected areas. Therefore, this event provides us with one of the most complete data set among tsunami events in history. In this study, fragility functions are derived using data provided by the Ministry of Land, Infrastructure and Transportation of Japan, with more than 250,000 structures surveyed. The set of data has details on damage level, structural material, number of stories per building and location (town). This information is crucial to the understanding of the causes of building damage, as differences in structural characteristics and building location can be taken into account in the damage probability analysis. Using least squares regression, different sets of fragility curves are derived to demonstrate the influence of structural material, number of stories and coastal topography on building damage levels. The results show a better resistant performance of reinforced concrete and steel buildings over wood or masonry buildings. Also, buildings taller than two stories were confirmed to be much stronger than the buildings of one or two stories. The damage characteristic due to the coastal topography based on limited number of data in town locations is also shortly discussed here. At the same tsunami inundation depth, buildings along the Sanriku ria coast were much greater damaged than buildings from the plain coast in Sendai. The difference in damage states can be explained by the faster flow velocities in the ria coast at the same inundation depth. These findings are key to support better future building damage assessments, land use management and disaster planning.     

Indicator mineralogy of kimberlite boulders from eskers in the Kirkland Lake and Lake Timiskaming areas, Ontario, Canada

Sixteen kimberlite boulders were collected from three sites on the Munro and Misema River Eskers in the Kirkland Lake kimberlite field and one site on the Sharp Lake esker in the Lake Timiskaming kimberlite field. The boulders were processed for heavy-mineral concentrates from which grains of Mg-ilmenite, chromite, garnet, clinopyroxene and olivine were picked, counted and analyzed by electron microprobe. Based on relative abundances and composition of these mineral phases, the boulders could be assigned to six mineralogically different groups, five for the Kirkland Lake area and one for the Lake Timiskaming area. Their indicator mineral composition and abundances are compared to existing data for known kimberlites in both the Kirkland Lake and Lake Timiskaming areas. Six boulders from the Munro Esker form a compositionally homogeneous group (I) in which the Mg-ilmenite population is very similar to that of the A1 kimberlite, located 7–12 km N (up-ice), directly adjacent to the Munro esker in the Kirkland Lake kimberlite field. U–Pb perovskite ages of three of the group I boulders overlap with that of the A1 kimberlite. Three other boulders recovered from the same localities in the Munro Esker also show some broad similarities in Mg-ilmenite composition and age to the A1 kimberlite. However, they are sufficiently different in mineral abundances and composition from each other and from the A1 kimberlite to assign them to different groups (II–IV). Their sources could be different phases of the same kimberlite or—more likely—three different, hitherto unknown kimberlites up-ice of the sample localities along the Munro Esker in the Kirkland Lake kimberlite field. A single boulder from the Misema River esker, Kirkland Lake, has mineral compositions that do not match any of the known kimberlites from the Kirkland Lake field. This suggests another unknown kimberlite exists in the area up-ice of the Larder Lake pit along the Misema River esker. Six boulders from the Sharp Lake esker, within the Lake Timiskaming field, form a homogeneous group with distinct mineral compositions unmatched by any of the known kimberlites in the Lake Timiskaming field. U–Pb perovskite age determinations on two of these boulders support this notion. These boulders are likely derived from an unknown kimberlite source up-ice from the Seed kimberlite, 4 km NW of the Sharp Lake pit, since indicator minerals with identical compositions to those of the Sharp Lake boulders have been found in till samples collected down-ice from Seed. Based on abundance and composition of indicator minerals, most importantly Mg-ilmenite, and supported by U–Pb age dating of perovskite, we conclude that the sources of 10 of the 16 boulders must be several hitherto unknown kimberlite bodies in the Kirkland Lake and Lake Timiskaming kimberlite fields.     

Measuring impact crater depth throughout the solar system

One important, almost ubiquitous, tool for understanding the surfaces of solid bodies throughout the solar system is the study of impact craters. While measuring a distribution of crater diameters and locations is an important tool for a wide variety of studies, so too is measuring a crater's “depth.” Depth can inform numerous studies including the strength of a surface and modification rates in the local environment. There is, however, no standard data set, definition, or technique to perform this data-gathering task, and the abundance of different definitions of “depth” and methods for estimating that quantity can lead to misunderstandings in and of the literature. In this review, we describe a wide variety of data sets and methods to analyze those data sets that have been, are currently, or could be used to derive different types of crater depth measurements. We also recommend certain nomenclature in doing so to help standardize practice in the field. We present a review section of all crater depths that have been published on different solar system bodies which shows how the field has evolved through time and how some common assumptions might not be wholly accurate. We conclude with several recommendations for researchers which could help different data sets to be more easily understood and compared.     

A New Appraisal of Sri Lankan BB Zircon as a Reference Material for LA‐ICP‐MS U‐Pb Geochronology and Lu‐Hf Isotope Tracing

A potential zircon reference material (BB zircon) for laser ablation‐inductively coupled plasma‐mass spectrometry (LA‐ICP‐MS) U‐Pb geochronology and Hf isotope geochemistry is described. A batch of twenty zircon megacrysts (0.5–1.5 cm³) from Sri Lanka was studied. Within‐grain rare earth element (REE) compositions are largely homogeneous, albeit with some variation seen between fractured and homogeneous domains. Excluding fractured cathodoluminescence bright domains, the variation in U content for all analysed crystals ranged from 227 to 368 μg g^?1 and the average Th/U ratios were between 0.20 and 0.47. The Hf isotope composition (0.56–0.84 g/100 g Hf) is homogeneous within and between the grains – mean ¹⁷⁶Hf/¹⁷⁷Hf of 0.281674 ± 0.000018 (2s). The calculated alpha dose of 0.59 × 10¹⁸ g^?1 for a number of BB grains falls within the trend of previously studied, untreated zircon samples from Sri Lanka. Aliquots of the same crystal (analysed by ID‐TIMS in four different laboratories) gave consistent U‐Pb ages with excellent measurement reproducibility (0.1–0.4% RSD). Interlaboratory assessment (by LA‐ICP‐MS) from individual crystals returned results that are within uncertainty equivalent to the TIMS ages. Finally, we report on within‐ and between‐grain homogeneity of the oxygen isotope systematic of four BB crystals (13.16‰ VSMOW).