Trace element diffusion and incorporation in quartz during heating experiments

Heating of quartz crystals in order to study melt and high-temperature fluid inclusions is a common practice to constrain major physical and chemical parameters of magmatic and hydrothermal processes. Diffusion and modification of trace element content in quartz and its hosted melt inclusions have been investigated through step-heating experiments of both matrix-free quartz crystals and quartz crystals associated with sulfides and other minerals using a Linkam TS1500 stage. Magmatic and hydrothermal quartz were successively analyzed after each heating step for Cu, Al, and Ti using electron probe micro-analyzer. After the last heating step, quartz crystals and their hosted melt inclusions were analyzed by laser ablation inductively coupled plasma mass spectrometry and compared to unheated samples. Heated samples reveal modification of Cu, Li, Na, and B contents in quartz and modification of Cu, Li, Ag, and K concentrations in melt inclusions. Our results show that different mechanisms of Cu, Li, and Na incorporation occur in magmatic and hydrothermal quartz. Heated magmatic quartz records only small, up to a few ppm, enrichment in Cu and Na, mostly substituting for Li. By contrast, heated hydrothermal quartz shows enrichment up to several hundreds of ppm in Cu, Li, and Na, which substitute for originally present H. This study reveals that the composition of both quartz and its hosted melt inclusions may be significantly modified upon heating experiments, leading to erroneous quantification of elemental concentrations. In addition, each quartz crystal also becomes significantly enriched in Cu in the sub-surface layer during heating. We propose that sub-surface Cu enrichment is a direct indication of Cu diffusion in quartz externally sourced from both the surrounding sulfides as well as the copper pins belonging to the heating device. Our study shows that the chemical compositions of both heated quartz and its hosted inclusions must be interpreted with great caution to avoid misleading geological interpretations.     

Spatial and seasonal variability of polygonal tundra water balance: Lena River Delta, northern Siberia (Russia)

The summer water balance of a typical Siberian polygonal tundra catchment is investigated in order to identify the spatial and temporal dynamics of its main hydrological processes. The results show that, besides precipitation and evapotranspiration, lateral flow considerably influences the site-specific hydrological conditions. The prominent microtopography of the polygonal tundra strongly controls lateral flow and storage behaviour of the investigated catchment. Intact rims of low-centred polygons build hydrological barriers, which release storage water later in summer than polygons with degraded rims and troughs above degraded ice wedges. The barrier function of rims is strongly controlled by soil thaw, which opens new subsurface flow paths and increases subsurface hydrological connectivity. Therefore, soil thaw dynamics determine the magnitude and timing of subsurface outflow and the redistribution of storage within the catchment. Hydraulic conductivities in the elevated polygonal rims sharply decrease with the transition from organic to mineral layers. This interface causes a rapid shallow subsurface drainage of rainwater towards the depressed polygon centres and troughs. The re-release of storage water from the centres through deeper and less conductive layers helps maintain a high water table in the surface drainage network of troughs throughout the summer.     

Characterizing Major Controls on Spatial and Seasonal Variations in Chemical Composition of Surface and Pore Brine of Maharlu Lake,Southern Iran

Maharlu Lake with Na–Cl water type is the terminal point of a closed basin in southern Iran. A total of 10 water samples from two rivers discharging to the lake and 78 water samples of surface and pore brine of Maharlu Lake have been collected from different depths (surface, 20, 50 and 100 cm) of four sampling stations along the lake during a period of lake water-level fluctuation (November 2014–July 2015). To investigate chemical interaction between lake surface water and shallow pore water and to understand the major factors governing chemical composition of Maharlu brine, concentrations of major and minor (boron, bromide and lithium) solutes, pH and total dissolved solids have been measured in collected water samples. Saturation indices of evaporite minerals in collected water samples have been also calculated. The chemical behavior of dissolved solutes and evaporative evolution of the lake brine during a hydrological period have been simulated using PHREEQC. The results of our investigations indicated that chemical composition of lake surface water and pore brine of Maharlu Lake are mainly connected with lake water-level fluctuations and distance from input rivers (and depth), respectively. Hydrochemical investigations and statistical analysis showed that the brines chemistry of Maharlu is mainly controlled by three processes: brine evaporative evolution, dissolution–precipitation and diagenetic evolution of secondary carbonates.     

Evaluation of the third- and fourth-generation GOCE Earth gravity field models with Australian terrestrial gravity data in spherical harmonics

In March 2013, the fourth generation of European Space Agency’s (ESA) global gravity field models, DIR4 (Bruinsma et al. in Proceedings of the ESA living planet symposium, 28 June–2 July, Bergen, ESA, Publication SP-686, 2010b) and TIM4 (Migliaccio et al. in Proceedings of the ESA living planet symposium, 28 June–2 July, Bergen, ESA, Publication SP-686, 2010), generated from the Gravity field and steady-state Ocean Circulation Explorer (GOCE) gravity observation satellite was released. We evaluate the models using an independent ground truth data set of gravity anomalies over Australia. Combined with Gravity Recovery and Climate Experiment (GRACE) satellite gravity, a new gravity model is obtained that is used to perform comparisons with GOCE models in spherical harmonics. Over Australia, the new gravity model proves to have significantly higher accuracy in the degrees below 120 as compared to EGM2008 and seems to be at least comparable to the accuracy of this model between degree 150 and degree 260. Comparisons in terms of residual quasi-geoid heights, gravity disturbances, and radial gravity gradients evaluated on the ellipsoid and at approximate GOCE mean satellite altitude ( $h=250$  km) show both fourth generation models to improve significantly w.r.t. their predecessors. Relatively, we find a root-mean-square improvement of 39 % for the DIR4 and 23 % for TIM4 over the respective third release models at a spatial scale of 100 km (degree 200). In terms of absolute errors, TIM4 is found to perform slightly better in the bands from degree 120 up to degree 160 and DIR4 is found to perform slightly better than TIM4 from degree 170 up to degree 250. Our analyses cannot confirm the DIR4 formal error of 1 cm geoid height (0.35 mGal in terms of gravity) at degree 200. The formal errors of TIM4, with 3.2 cm geoid height (0.9 mGal in terms of gravity) at degree 200, seem to be realistic. Due to combination with GRACE and SLR data, the DIR models, at satellite altitude, clearly show lower RMS values compared to TIM models in the long wavelength part of the spectrum (below degree and order 120). Our study shows different spectral sensitivity of different functionals at ground level and at GOCE satellite altitude and establishes the link among these findings and the Meissl scheme (Rummel and van Gelderen in Manusrcipta Geodaetica 20:379–385, 1995).