‘Modelling the Arctic Boundary Layer: An Evaluation of Six Arcmip Regional-Scale Models using Data from the Sheba Project’

A primary climate change signal in the central Arctic is the melting of sea ice. This is dependent on the interplay between the atmosphere and the sea ice, which is critically dependent on the exchange of momentum, heat and moisture at the surface. In assessing the realism of climate change scenarios it is vital to know the quality by which these exchanges are modelled in climate simulations. Six state-of-the-art regional-climate models are run for one year in the western Arctic, on a common domain that encompasses the Surface Heat Budget of the Arctic Ocean (SHEBA) experiment ice-drift track. Surface variables, surface fluxes and the vertical structure of the lower troposphere are evaluated using data from the SHEBA experiment. All the models are driven by the same lateral boundary conditions, sea-ice fraction and sea and sea-ice surface temperatures. Surface pressure, near-surface air temperature, specific humidity and wind speed agree well with observations, with a falling degree of accuracy in that order. Wind speeds have systematic biases in some models, by as much as a few metres per second. The surface radiation fluxes are also surprisingly accurate, given the complexity of the problem. The turbulent momentum flux is acceptable, on average, in most models, but the turbulent heat fluxes are, however, mostly unreliable. Their correlation with observed fluxes is, in principle, insignificant, and they accumulate over a year to values an order of magnitude larger than observed. Typical instantaneous errors are easily of the same order of magnitude as the observed net atmospheric heat flux. In the light of the sensitivity of the atmosphere–ice interaction to errors in these fluxes, the ice-melt in climate change scenarios must be viewed with considerable caution.     

The response of flow duration curves to afforestation

The hydrologic effect of replacing pasture or other short crops with trees is reasonably well understood on a mean annual basis. The impact on flow regime, as described by the annual flow duration curve (FDC) is less certain. A method to assess the impact of plantation establishment on FDCs was developed. The starting point for the analyses was the assumption that rainfall and vegetation age are the principal drivers of evapotranspiration. A key objective was to remove the variability in the rainfall signal, leaving changes in streamflow solely attributable to the evapotranspiration of the plantation. A method was developed to (1) fit a model to the observed annual time series of FDC percentiles; i.e. 10th percentile for each year of record with annual rainfall and plantation age as parameters, (2) replace the annual rainfall variation with the long term mean to obtain climate adjusted FDCs, and (3) quantify changes in FDC percentiles as plantations age. Data from 10 catchments from Australia, South Africa and New Zealand were used. The model was able to represent flow variation for the majority of percentiles at eight of the 10 catchments, particularly for the 10–50th percentiles. The adjusted FDCs revealed variable patterns in flow reductions with two types of responses (groups) being identified. Group 1 catchments show a substantial increase in the number of zero flow days, with low flows being more affected than high flows. Group 2 catchments show a more uniform reduction in flows across all percentiles. The differences may be partly explained by storage characteristics. The modelled flow reductions were in accord with published results of paired catchment experiments. An additional analysis was performed to characterise the impact of afforestation on the number of zero flow days (N_zero) for the catchments in group 1. This model performed particularly well, and when adjusted for climate, indicated a significant increase in N_zero. The zero flow day method could be used to determine change in the occurrence of any given flow in response to afforestation. The methods used in this study proved satisfactory in removing the rainfall variability, and have added useful insight into the hydrologic impacts of plantation establishment. This approach provides a methodology for understanding catchment response to afforestation, where paired catchment data is not available.     

Ultrahigh pressure metamorphic rocks from the Chinese Continental Scientific Drilling Project: I. Petrology and geochemistry of the main hole (0–2,050 m)

The main hole (MH) of the Chinese Continental Scientific Drilling Project (CCSD) in southern Sulu has penetrated into an ultrahigh-pressure (UHP) metamorphic rock slice which consists of orthogneiss, paragneiss, eclogite, ultramafic rock and minor schist. Recovered eclogites have a UHP metamorphic mineral assemblage of garnet + omphacite + rutile ± phengite ± kyanite ± coesite ± epidote. Ultramafic rocks contain garnet + olivine + clinopyroxene + orthopyroxene ± Ti-clinohumite ± phlogopite. Gneisses and schists contain an amphibolite-facies paragenesis, but their zircons have coesite, garnet, omphacite (or jadeite) and phengite inclusions, indicating that eclogites and gneisses have been subjected to in situ UHP metamorphism. Using available geothermobarometers, P–T estimates of 3.1–4.4 GPa and 678–816°C for eclogites were obtained. If surface outcrops and neighboring shallow drill holes are considered together, we suggest that a huge supracrustal rock slab (> 50 km long × 100 km wide × 5 km deep) was subducted to a depth > 100 km and then exhumed to the surface. The depth interval (0–2,050 m) of the CCSD-MH can be divided into six lithological units. Unit 1 consists of alternating layers of quartz-rich and rutile-rich eclogites, with thin interlayers of gneiss and schist. Eclogites of unit 1 are characterized by Nb, Ta, Sr and Ti depletions, low Mg number and general LREE enrichment. Unit 2 comprises rutile- and ilmenite-rich eclogite and minor “normal” eclogite and is characterized by high TiO₂, total Fe, V, Co and Sr, and very low SiO₂, alkali, Zr, Ba, Nb, Ta and total REE contents, and LREE-depleted REE patterns with slightly positive Eu anomalies. Unit 3 contains ultramafic rock and minor MgO-rich eclogite. Protoliths of UHP rocks from units 1, 2 and 3 represent a layered mafic to ultramafic intrusion at crustal depth. Units 4 and 6 consist of interlayered eclogite and paragneiss; the eclogites are characterized by Th, U, Nb, Ta and Ti depletion and K enrichment and LREE-enriched REE patterns. Paragneisses show Nb, Ta, Sr and Ti depletions and LREE-enriched REE patterns occasionally with slightly negative Eu anomalies, indicating that their protoliths represent metamorphic supracrustal series. Unit 5 consists mainly of orthogneisses, showing distinct Nb, Ta, Sr and Ti depletions, and LREE-enriched REE patterns with pronounced negative Eu anomalies, suggesting granitic protoliths. In conclusion it is proposed that the southern Sulu UHP belt consists of a series of meta-supracrustal rocks, a layered mafic–ultramafic complex and granites.     

Acapulcoite-lodranite meteorites: Ultramafic asteroidal partial melt residues

Acapulcoites (most ancient Hf-W ages are 4,563.1?±?0.8 Ma), lodranites (most ancient Hf-W ages are 4,562.6?±?0.9 Ma) and rocks transitional between them are ancient residues of different degrees of partial melting of a chondritic source lithology (e.g., as indicated by the occurrence of relict chondrules in 9 acapulcoites), although the precise chondrite type is unknown. Acapulcoites are relatively fine- grained (～150–230?μm) rocks with equigranular, achondritic textures and consist of olivine, orthopyroxene, Ca-rich clinopyroxene, plagioclase, metallic Fe,Ni, troilite, chromite and phosphates. Lodranites are coarser grained (540–700?μm), with similar equigranular, recrystallized textures, mineral compositions and contents, although some are significantly depleted in eutectic Fe,Ni-FeS and plagioclase- clinopyroxene partial melts. The acapulcoite-lodranite clan is most readily distinguished from other groups of primitive achondrites (e.g., winoanites/IAB irons) by oxygen isotopic compositions, although more than 50% of meteorites classified as acapulcoites currently lack supporting oxygen isotopic data. The heat source for melting of acapulcoites-lodranites was internal to the parent body, most likely ²⁶Al, although some authors suggest it was shock melting. Acapulcoites experienced lower temperatures of ～980–1170?°C and lower degrees of partial melting (～1–4?vol.%) and lodranites higher temperatures of ～1150–1200?°C and higher degrees (～5?≥?10?vol.%) of partial melting. Hand-specimen and thin section observations indicate movement of Fe,Ni-FeS, basaltic, and phosphate melts in veins over micrometer to centimeter distances. Mineralogical, chemical and isotopic properties, Cosmic Ray Exposure (CRE) ages which cluster around 4–6 Ma and the occurrence of some meteorites consisting of both acapulcoite and lodranite material, indicate that these meteorites come from one parent body and were most likely ejected in one impact event. Whereas the precise parent asteroid of these meteorites is unknown, there is general agreement that it was an S-type object. There is nearly total agreement that the acapulcoite-lodranite parent body was <～100?km in radius and, based on the precise Pb–Pb age for Acapulco of 4555.9?±?0.6 Ma, combined with the Hf/W and U/Pb records and cooling rates deduced from mineralogical and other investigations, that the parent body was fragmented during its cooling which the U/Pb system dates at precisely 4556?±?1 Ma. Hf-W chronometry suggests that the parent body of the acapulcoites-lodranites and, in fact, the parent bodies of all “primitive achondrites” accreted slightly later than those of the differentiated achondrites and, thus, had lower contents of ²⁶Al, the heat producing radionuclide largely responsible for heating of both primitive and differentiated achondrites. Thus, the acapulcoite-lodranite parent body never experienced the high degrees of melting responsible for the formation of the differentiated meteorites, but arrested its melting history at relatively low degrees of ～15?vol.%.     

Calcium Carbonate and Phosphate Reference Materials for Monitoring Bulk and Microanalytical Determination of Sr Isotopes

In situ laser ablation analyses rely on the microanalytical homogeneity of reference materials (RMs) and a similar matrix and mass fraction between unknown samples and RMs to obtain reliable results. Suitable carbonate and phosphate RMs for determination of Sr isotope ratios in such materials are limited. Thus, we determined ⁸⁷Sr/⁸⁶Sr ratios of several carbonate (JCt‐1, JCp‐1, MACS‐1, MACS‐3) and phosphate (MAPS‐4, MAPS‐5, NIST SRM 1400, NIST SRM 1486) international RMs using dissolved samples and two different multi‐collector inductively coupled plasma‐mass spectrometers (MC‐ICP‐MS). Our Sr isotope data are in agreement with published data and have an improved measurement precision for some RMs. For MACS‐1, we present the first ⁸⁷Sr/⁸⁶Sr value. We tested the suitability of these materials for microanalytical analyses by LA‐MC‐ICP‐MS, with two different laser ablation systems: a conventional nanosecond laser and a state‐of‐the‐art femtosecond laser. We investigated the RMs micro‐homogeneity and compared the data with our solution data. Both laser ablation systems yielded identical ⁸⁷Sr/⁸⁶Sr ratios within uncertainty to the solution data for RMs with low interferences of REEs. Therefore, these carbonate and phosphate RMs can be used to achieve accurate and precise results for in situ Sr isotope investigations by LA‐MC‐ICP‐MS of similar materials.     

Review paper: What Kelvin might have written about Elasticity

One hundred and fifty five years ago, Kelvin published the first part of a fundamental analysis of the elastic tensor, in which he proposed a coordinate‐free representation through its eigensystem. His thoughts were apparently far ahead of his time, since it took 125 years before the paper elicited a positive reaction (it is now accessible through several modern reviews). Science not only lost track for 125 years of the original paper but also lost the ideas Kelvin might have proposed in the second part, a publication that was never put to paper, presumably in view of the lack of appreciation of the first part. In an attempt to establish what might have been on Kelvin's mind for a second part, one has to ‘forget’ the progress of mathematical physics in the intervening time and base all arguments strictly on the content of the first part and on the state of science in the second half of the 19^thcentury. The theory of elasticity would certainly have developed faster, had Kelvin's paper peen appreciated by his ‘peers’. But a theory based on Kelvin's ideas would be fruitful even today.     

Research note: Seismic attenuation due to wave‐induced fluid flow at microscopic and mesoscopic scales

Wave‐induced fluid flow at microscopic and mesoscopic scales arguably constitutes the major cause of intrinsic seismic attenuation throughout the exploration seismic and sonic frequency ranges. The quantitative analysis of these phenomena is, however, complicated by the fact that the governing physical processes may be dependent. The reason for this is that the presence of microscopic heterogeneities, such as micro‐cracks or broken grain contacts, causes the stiffness of the so‐called modified dry frame to be complex‐valued and frequency‐dependent, which in turn may affect the viscoelastic behaviour in response to fluid flow at mesoscopic scales. In this work, we propose a simple but effective procedure to estimate the seismic attenuation and velocity dispersion behaviour associated with wave‐induced fluid flow due to both microscopic and mesoscopic heterogeneities and discuss the results obtained for a range of pertinent scenarios.     

Multi-stage Ag–Bi–Co–Ni–U and Cu–Bi vein mineralization at Wittichen,Schwarzwald, SW Germany: geological setting,ore mineralogy,and fluid evolution

The Wittichen Co–Ag–Bi–U mining area (Schwarzwald ore district, SW Germany) hosts several unconformity-related vein-type mineralizations within Variscan leucogranite and Permian to Triassic redbeds. The multistage mineralization formed at the intersection of two fault systems in the last 250 Ma. A Permo-Triassic ore stage I with minor U–Bi–quartz–fluorite mineralization is followed by a Jurassic to Cretaceous ore stage II with the main Ag and Co mineralization consisting of several generations of gangue minerals that host the sub-stages of U–Bi, Bi–Ag, Ni–As–Bi and Co–As–Bi. Important ore minerals are native elements, Co and Ni arsenides, and pitchblende; sulphides are absent. The Miocene ore stage III comprises barite with the Cu–Bi sulfosalts emplectite, wittichenite and aikinite, and the sulphides anilite and djurleite besides native Bi, chalcopyrite, sphalerite, galena and tennantite. The mineral-forming fluid system changed from low salinity (<5 wt.% NaCl) at high temperature (around 300°C) in Permian to highly saline (around 25 wt.% NaCl + CaCl₂) at lower temperatures (50–150°C) in Triassic to Cretaceous times. Thermodynamic calculations and comparison with similar mineralizations worldwide show that the Mesozoic ore-forming fluid was alkaline with redox conditions above the hematite–magnetite buffer. We suggest that the precipitation mechanism for native elements, pitchblende and arsenides is a decrease in pH during fluid mixing processes. REE patterns in fluorite and the occurrence of Bi in all stages suggest a granitic source of some ore-forming elements, whereas, e.g. Ag, Co and Ni probably have been leached from the redbeds. The greater importance of Cu and isotope data indicates that the Miocene ore stage III is more influenced by fluids from the overlying redbeds and limestones than the earlier mineralization stages.     

Origin and extent of fresh groundwater, salty paleowaters and recent saltwater intrusions in Red River flood plain aquifers, Vietnam

A model has been established on the origin and extent of fresh groundwater, salty paleowaters and saltwater from recent seawater intrusions in the Red River flood plain in Vietnam. This was done with geological observations, geophysical borehole logging and transient electromagnetic methods. Salt paleowater is present up to 50–75?km from the coastline, with occurrence controlled by the Holocene transgression. A density-driven leaching of salty porewater has occurred from high-permeability Holocene sediments into underlying Pleistocene deposits, whereas diffusion has dominated in low-permeability layers. In the Pleistocene aquifer, the highest content of dissolved solids is found below two intrinsic valleys with Holocene marine sediments and along the coastline. Recent intrusion of saltwater from the South China Sea is observed in shallow groundwater 35?km inland, probably a result of transport of salty water inland in rivers or leaching of paleowaters from very young near-coast marine sediments. The observed inverted salinity profile, with high saline water overlying fresher groundwater, has been formed due to the global eustatic sea-level changes during the last 8,000–9,000?years. The proposed model may therefore be applicable to other coastal aquifers, with a proper incorporation of the local geological environments.