Lake Floor Morphology and Sediment Architecture of Lake Torneträsk,Northern Sweden

Here we present datasets from a hydroacoustic survey in July 2011 at Lake Torneträsk, northern Sweden. Our hydroacoustic data exhibit lake floor morphologies formed by glacial erosion and accumulation processes, insights into lacustrine sediment accumulation since the beginning of deglaciation, and information on seismic activity along the Pärvie Fault. Features of glacial scouring with a high‐energy relief, steep slopes, and relative reliefs of more than 50 m are observed in the large W‐basin. The remainder of the lacustrine subsurface appears to host a broad variety of well preserved formations from glacial accumulation related to the last retreat of the Fennoscandian ice sheet. Deposition of glaciolacustrine and lacustrine sediments is focused in areas situated in proximity to major inlets. Sediment accumulation in distal areas of the lake seldom exceeds 2 m or is not observable. We assume that lack of sediment deposition in the lake is a result of different factors, including low rates of erosion in the catchment, a previously high lake level leading to deposition of sediments in higher elevated paleodeltas, tributaries carrying low suspension loads as a result of sedimentation in upstream lakes, and an overall low productivity in the lake. A clear off‐shore trace of the Pärvie Fault could not be detected from our hydroacoustic data. However, an absence of sediment disturbance in close proximity to the presumed fault trace implies minimal seismic activity since deposition of the glaciolacustrine and lacustrine sediments.     

Sea‐level and ocean‐current control on carbonate‐platform growth,Maldives, Indian Ocean

Multichannel high‐resolution seismic and multibeam data were acquired from the Maldives‐isolated carbonate platform in the Indian Ocean for a detailed characterization of the Neogene bank architecture of this edifice. The goal of the research is to decipher the controlling factors of platform evolution, with a special emphasis on sea‐level changes and changes of the oceanic currents. The stacking pattern of Lower to Middle Miocene depositional sequences, with an evolution of a ramp geometry to a flat‐topped platform, reflects variations of accommodation, which here are proposed to be primarily governed by fluctuations of relative sea level. Easterly currents during this stage of bank growth controlled an asymmetric east‐directed progradation of the bank edge. During the late middle Miocene, this system was replaced by a twofold configuration of bank development. Bank growth continued synchronously with partial bank demise and associated sediment‐drift deposition. This turnover is attributed to the onset and/or intensification of the Indian monsoon and related upwelling and occurrence of currents, locally changing environmental conditions and impinging upon the carbonate system. Mega spill over lobes, shaped by reversing currents, formed as large‐scale prograding complexes, which have previously been interpreted as deposits formed during a forced regression. On a regional scale, a complex carbonate‐platform growth can occur, with a coexistence of bank‐margin progradation and aggradation, as well as partial drowning. It is further shown that a downward shift of clinoforms and offlapping geometries in carbonate platforms are not necessarily indicative for a sea‐level driven forced regression. Findings are expected to be applicable to other examples of Cenozoic platforms in the Indo‐Pacific region.     

Are numerical weather model outputs helpful to reduce tropospheric delay signals in InSAR data?

Interferometric synthetic aperture radar phase data include not only signals due to crustal movements, but also those associated with microwave propagation delay through the atmosphere. In particular, the effect of water vapor can generate apparent signals in the order of a few centimeters or more, and prevent us from detecting such geophysical signals as those due to secular crustal deformation. To examine if and to what extent numerical weather model (NWM) outputs are helpful to reduce the tropospheric delay signals at spatial scales of 5–50 km wavelengths, we compared three approaches of tropospheric signal reduction, using 54 interferograms in central Hokkaido, Japan. The first approach is the conventional topography-correlated delay correction that is based on the regional digital elevation model (DEM). The second approach is based on the Japan Meteorological Agency’s operational meso-scale analysis model (MSM) data, where we compute tropospheric delays and subtract them from the interferogram. However, the MSM data are available at predefined epochs and their spatial resolution is about 10 km; therefore, we need to interpolate both temporally and spatially to match with interferograms. Expecting to obtain a more physically plausible reduction of the tropospheric effects, we ran a 1-km mesh high-resolution numerical weather model WRF (Weather Research and Forecasting model) by ourselves, using the MSM data as the initial and boundary conditions. The third approach is similar to the second approach, except that we make use of the WRF-based tropospheric data. Results show that if the topography-correlated phases are significant, both the conventional DEM-based approach and the MSM-based approach reveal comparable performances. However, when the topography-correlated phases are insignificant, none of the approaches can efficiently reduce the tropospheric phases. Although it could reduce the tropospheric signals in a local area, in none of the case studies did the WRF model produce the “best” performance. Whereas the global atmospheric model outputs are shown to be effective in reducing long-wavelength tropospheric signals, we consider that further improvements are needed for the initial and boundary condition data for high-resolution NWM, so that the NWM-based approach will become more reliable even in the case of a non-stratified troposphere.     

Optimized formulas for the gravitational field of a tesseroid

Various tasks in geodesy, geophysics, and related geosciences require precise information on the impact of mass distributions on gravity field-related quantities, such as the gravitational potential and its partial derivatives. Using forward modeling based on Newton’s integral, mass distributions are generally decomposed into regular elementary bodies. In classical approaches, prisms or point mass approximations are mostly utilized. Considering the effect of the sphericity of the Earth, alternative mass modeling methods based on tesseroid bodies (spherical prisms) should be taken into account, particularly in regional and global applications. Expressions for the gravitational field of a point mass are relatively simple when formulated in Cartesian coordinates. In the case of integrating over a tesseroid volume bounded by geocentric spherical coordinates, it will be shown that it is also beneficial to represent the integral kernel in terms of Cartesian coordinates. This considerably simplifies the determination of the tesseroid’s potential derivatives in comparison with previously published methodologies that make use of integral kernels expressed in spherical coordinates. Based on this idea, optimized formulas for the gravitational potential of a homogeneous tesseroid and its derivatives up to second-order are elaborated in this paper. These new formulas do not suffer from the polar singularity of the spherical coordinate system and can, therefore, be evaluated for any position on the globe. Since integrals over tesseroid volumes cannot be solved analytically, the numerical evaluation is achieved by means of expanding the integral kernel in a Taylor series with fourth-order error in the spatial coordinates of the integration point. As the structure of the Cartesian integral kernel is substantially simplified, Taylor coefficients can be represented in a compact and computationally attractive form. Thus, the use of the optimized tesseroid formulas particularly benefits from a significant decrease in computation time by about 45 % compared to previously used algorithms. In order to show the computational efficiency and to validate the mathematical derivations, the new tesseroid formulas are applied to two realistic numerical experiments and are compared to previously published tesseroid methods and the conventional prism approach.     

Lettau’s Contribution to the Obukhov Length Scale: A Scientific Historical Study

Boundary-Layer Meteorology - It has been repeatedly assumed that Heinz Lettau found the Obukhov length in 1949 independently of Obukhov in 1946. However, it was not the characteristic length scale,...     

The nature and genesis of marginal Cu–PGE–Au sulphide mineralisation in Paleogene Macrodykes of the Kangerlussuaq region, East Greenland

The Kangerlussuaq region of East Greenland hosts a variety of early Tertiary extrusive and intrusive igneous rocks related to continental break up and the passage of the ancestral Iceland plume. These intrusive bodies include a number of gabbroic macrodykes, two of which—the Miki Fjord Macrodyke, and the newly discovered Togeda Macrodyke—contain Cu–PGE–Au sulphide mineralisation along their margins. Sulphides occur as disseminated interstitial blebs and rounded globules of chalcopyrite and pyrrhotite with some Fe–Ti oxides and platinum-group minerals, comprising largely Pd bismuthides and tellurides. The globules are interpreted to have formed from fractionation of trapped droplets of an immiscible Cu- and Pd-rich sulphide melt and show geopetal indicators. Sulphur isotopes imply a local crustal source of S in these from pyritic sediments of the Kangerlussuaq Basin. Thus, generation of these sulphide occurrences was controlled by local country rock type. Low Ni/Cu and Pt/Pd ratios, also present in the Platinova reefs in the Skaergaard Intrusion, indicate that early fractionation of olivine may have depleted the magma of Ni and suggest the likely presence of a large magma chamber at depth. Xenoliths of Ni-rich olivine cumulates in the Miki Fjord Macrodyke may have been sourced from such a body. The location of thus far unidentified conduit or feeder zones to the macrodykes beneath the present day surface may represent potential targets for more massive sulphide orebodies.     

Fluid inclusion study of the Wunugetu Cu–Mo deposit,Inner Mongolia,China

Hydrothermal alteration and mineralization at the Wunugetu porphyry Cu–Mo deposit, China, include four stages, i.e., the early stage characterized by quartz, K-feldspar and minor mineralization, followed by a molybdenum mineralization stage associated with potassic alteration, copper mineralization associated with sericitization, and the last Pb–Zn mineralization stage associated with carbonation. Hydrothermal quartz contains three types of fluid inclusions, namely aqueous (W-type), daughter mineral-bearing (S-type) and CO₂-rich (C-type) inclusion, with the latter two types absent in the late stage. Fluid inclusions in the early stage display homogenization temperatures above 510°C, with salinities up to 75.8 wt.% NaCl equivalent. The presence of S-type inclusions containing anhydrite and hematite daughter minerals and C-type inclusions indicates an oxidizing, CO₂-bearing environment. Fluid inclusions in the Mo- and Cu-mineralization stages yield homogenization temperatures of 342–508°C and 241–336°C, and salinities of 8.6–49.4 and 6.3–35.7 wt.% NaCl equivalent, respectively. The presence of chalcopyrite instead of hematite and anhydrite daughter minerals in S-type inclusions indicates a decreasing of oxygen fugacity. In the late stage, fluid inclusions yield homogenization temperatures of 115–234°C and salinities lower than 12.4 wt.% NaCl equivalent. It is concluded that the early stage fluids were CO₂ bearing, magmatic in origin, and characterized by high temperature, high salinity, and high oxygen fugacity. Phase separation occurred during the Mo- and Cu-mineralization stages, resulting in CO₂ release, oxygen fugacity decrease and rapid precipitation of sulfides. The late-stage fluids were meteoric in origin and characterized by low temperature, low salinity, and CO₂ poor.     

The competing models for the origin and internal evolution of granitic pegmatites in the light of melt and fluid inclusion research

In this paper we discuss the main petrogenetic models for granitic pegmatites and how these models have evolved over time. We suggest that the present state of knowledge requires that some aspects of these models to be modified, or absorbed into newer ones. Pegmatite formation and internal evolution have long supposed the need for highly water- and flux-enriched magmas to explain the differences between pegmatites and other intrusives of similar major element composition. Compositions and textural characteristics of fluid and melt inclusions in pegmatite minerals provide strong evidence for such magmas. Furthermore, we show that melt inclusion research has increased the number of potential flux components, which may include H₂O, OH^?, CO₂, HCO ₃ ^? , CO ₃ ^2? , SO ₄ ^2? , PO ₄ ^3? , H₃BO₃, F , and Cl, as well as the elements Li, Na, K, Rb, Cs, and Be, herein described as melt structure modifiers. In this paper we emphasize that the combined effect which these components have on the properties of pegmatite melts is difficult to deduce from experimental studies using only a limited number of these components. The combination and the amount of the different magmatic species, together with differences in the source region, and variations in pressure and temperature cause the great diversity of the pegmatites observed. Some volatile species, such as CO ₃ ^2? and alkalis, have the capacity to increase the solubility of H₂O in silicate melt to an extraordinary degree, to the extent that melt-melt-fluid immiscibility becomes inevitable. It is our view that the formation of pegmatites is connected with the complex interplay of many factors.     

Sixty thousand years of magmatic volatile history before the caldera-forming eruption of Mount Mazama, Crater Lake, Oregon

The well-documented eruptive history of Mount Mazama, Oregon, provides an excellent opportunity to use pre-eruptive volatile concentrations to study the growth of an explosive silicic magmatic system. Melt inclusions (MI) hosted in pyroxene and plagioclase crystals from eight dacitic–rhyodacitic eruptive deposits (71–7.7?ka) were analyzed to determine variations in volatile-element concentrations and changes in magma storage conditions leading up to and including the climactic eruption of Crater Lake caldera. Temperatures (Fe–Ti oxides) increased through the series of dacites, then decreased, and increased again through the rhyodacites (918–968 to ~950 to 845–895?°C). Oxygen fugacity began at nickel–nickel-oxide buffer (NNO) +0.8 (71?ka), dropped slightly to NNO +0.3, and then climbed to its highest value with the climactic eruption (7.7?ka) at NNO +1.1 log units. In parallel with oxidation state, maximum MI sulfur concentrations were high early in the eruptive sequence (~500?ppm), decreased (to ~200?ppm), and then increased again with the climactic eruption (~500?ppm). Maximum MI sulfur correlates with the Sr content (as a proxy for LREE, Ba, Rb, P₂O₅) of recharge magmas, represented by basaltic andesitic to andesitic enclaves and similar-aged lavas. These results suggest that oxidized Sr-rich recharge magmas dominated early and late in the development of the pre-climactic dacite–rhyodacite system. Dissolved H₂O concentrations in MI do not, however, correlate with these changes in dominant recharge magma, instead recording vapor solubility relations in the developing shallow magma storage and conduit region. Dissolved H₂O concentrations form two populations through time: the first at 3–4.6 wt% (with a few extreme values up to 6.1 wt%) and the second at ≤2.4 wt%. CO₂ concentrations measured in a subset of these inclusions reach up to 240?ppm in early-erupted deposits (71?ka) and are below detection in climactic deposits (7.7?ka). Combined H₂O and CO₂ concentrations and solubility models indicate a dominant storage region at 4–7?km (up to 12?km), with drier inclusions that diffusively re-equilibrated and/or were trapped at shallower depths. Boron and Cl (except in the climactic deposit) largely remained in the melt, suggesting vapor–melt partition coefficients and gas fractions were low. Modeled Li, F, and S vapor–melt partition coefficients are higher than those of B and Cl. The decrease in maximum MI CO₂ concentration following the earliest dacitic eruptions is interpreted to result from a broadening of the shallow storage region to greater than the diameter of subjacent feeders, so that greater proportions of reservoir magma were to the side of CO₂-bearing vapor bubbles ascending vertically from the locus of recharge magma injection, thereby escaping recarbonation by streaming vapor bubbles. The Mazama melt inclusions provide a picture of a growing magma storage region, where chemical variations in melt and magma occur due to changes in the nature and supply rate of magma recharge, the timing of degassing, and the possible degree of equilibration with gases from below.