古翼手龙类的德国翼手龙新标本

The pterodactyloid genus Germanodactylus is relatively poorly known. Two species have been described, although it has been recently suggested (Wang et al., 2008) that they might pertain to different genera.     

Cross-covariances and localization for EnKF in multiphase flow data assimilation

The ensemble Kalman filter has been successfully applied for data assimilation in very large models, including those in reservoir simulation and weather. Two problems become critical in a standard implementation of the ensemble Kalman filter, however, when the ensemble size is small. The first is that the ensemble approximation to cross-covariances of model and state variables to data can indicate the presence of correlations that are not real. These spurious correlations give rise to model or state variable updates in regions that should not be updated. The second problem is that the number of degrees of freedom in the ensemble is only as large as the size of the ensemble, so the assimilation of large amounts of precise, independent data is impossible. Localization of the Kalman gain is almost universal in the weather community, but applications of localization for the ensemble Kalman filter in porous media flow have been somewhat rare. It has been shown, however, that localization of updates to regions of non-zero sensitivity or regions of non-zero cross-covariance improves the performance of the EnKF when the ensemble size is small. Localization is necessary for assimilation of large amounts of independent data. The problem is to define appropriate localization functions for different types of data and different types of variables. We show that the knowledge of sensitivity alone is not sufficient for determination of the region of localization. The region depends also on the prior covariance for model variables and on the past history of data assimilation. Although the goal is to choose localization functions that are large enough to include the true region of non-zero cross-covariance, for EnKF applications, the choice of localization function needs to balance the harm done by spurious covariance resulting from small ensembles and the harm done by excluding real correlations. In this paper, we focus on the distance-based localization and provide insights for choosing suitable localization functions for data assimilation in multiphase flow problems. In practice, we conclude that it is reasonable to choose localization functions based on well patterns, that localization function should be larger than regions of non-zero sensitivity and should extend beyond a single well pattern.     

Do S-type granites commonly sample infracrustal sources? New results from an integrated O, U–Pb and Hf isotope study of zircon

In contrast to I-type granites, which commonly comprise infracrustal and supracrustal sources, S-type granites typically incorporate predominantly supracrustal sources. The initial aim of this study was to identify the sources of three Scottish Caledonian (~460 Ma) S-type granites (Kemnay, Cove and Nigg Bay) by conducting oxygen, U–Pb and Hf isotope analyses in zircon in order to characterise one potential end-member magma involved in the genesis of the voluminous late Caledonian (~430–400 Ma) I-type granites. Field, whole-rock geochemical and isotopic data are consistent with the generation of the S-type granites by melting their Dalradian Supergroup country rocks. While Hf isotope compositions of magmatic zircon, U–Pb data of inherited zircons, and high mean zircon δ¹⁸O values of 9.0 ± 2.7‰ (2SD) and 9.8 ± 2.0‰ for the Kemnay and Cove granites support this model, the Nigg Bay Granite contains zircons with much lower δ¹⁸O values (6.8 ± 2.1‰), similar to those found in Scottish I-type granites. This suggests that the Nigg Bay Granite contains low-δ¹⁸O material representing either altered supracrustal material, or more likely, an infracrustal source component with mantle-like δ¹⁸O. Mixing trends in plots of δ¹⁸O vs. εHf for S-type granite zircons indicate involvement of at least two sources in all three granites. This pilot study of Scottish Caledonian S-type granites demonstrates that, while field and whole-rock geochemical data are consistent with local melting of only supracrustal sources, the oxygen isotopic record stored in zircon reveals a much more complex petrogenetic evolution involving two or more magma sources.     

Weathering, dust, and biocycling effects on soil silicon isotope ratios

Silicon isotope ratios (δ³⁰Si) of bulk mineral materials in soil integrate effects from both silicon sources and processing. Here we report δ³⁰Si values from a climate gradient of Hawaiian soils developed on 170 ka basalt and relate them to patterns of soil chemistry and mineralogy. The results demonstrate informative relationships between the mass fraction of soil Si depletion and δ³⁰Si. In upper (<1 m deep) soil horizons along the climate gradient, Si depletion correlates with decreases of residual δ³⁰Si values in low rainfall soils and increases in high rainfall soils. Strong positive correlation between soil δ³⁰Si and dust-derived quartz and mica content show that both trends are largely controlled by the abundance of these weathering-resistant minerals. The data also lend support to the idea that fractionation of Si isotopes in secondary phases is controlled by partitioning of silicon between dissolved and precipitated products during the initial weathering of primary basalt. Secondary mineral δ³⁰Si values from lower (>1 m deep) soil horizons generally correlate with the isotope fractionation predicted by a study of dissolved Si in basalt-watershed rivers and driven by preferential ²⁸Si removal from the dissolved phase during precipitation. In contrast, after correcting for the influence of dust, secondary mineral Si depletion and δ³⁰Si values in shallow (<1 m deep) soil horizons showed evidence of biocycling induced Si redistribution and substantially lower δ³⁰Si values than predicted. Low δ³⁰Si values in shallow soil horizons compared to predictions can be attributed to repeated fractionation as secondary minerals undergo additional cycles of dissolution and precipitation. Primary mineral weathering, secondary mineral weathering, dust accumulation, and biocycling are major processes in terrestrial Si cycling and these results demonstrate that each can be traced by δ³⁰Si values interpreted in conjunction with mineralogy and measures of Si depletion.     

High-precision high field strength element partitioning between garnet, amphibole and alkaline melt from Kakanui, New Zealand

The high field strength elements (HFSE: Zr, Hf, Nb, Ta, and W) are an important group of chemical tracers that are increasingly used to investigate magmatic differentiation processes. Successful modeling of these processes requires the availability of accurate mineral-melt partition coefficients (D). To date, these have largely been determined by ion microprobe or laser ablation-ICP-MS analyses of the run products of high-pressure, high-temperature experiments. Since HFSE are (highly) incompatible, relatively immobile, high-charge, and difficult to ionize, these experiments and their analysis are challenging. Here we explore whether high-precision analyses of natural mineral-melt systems can provide additional constraints on HFSE partitioning.The HFSE concentrations in natural garnet and amphibole and their alkaline host melt from Kakanui, New Zealand are determined with high precision isotope dilution on a multi-collector-ICP-MS. Major and trace element compositions combined with Lu-Hf isotopic systematics and detailed petrographic sample analysis are used to assess mineral-melt equilibrium and to provide context for the HFSE D measurements. The whole-rock nephelinite, ∼1 mm sized amphiboles in the nephelinite, and garnet megacrysts have similar initial Hf isotope ratios with a mean initial ¹⁷⁶Hf/¹⁷⁷Hf_{(34 Ma)} = 0.282900 ± 0.000026 (2σ). In contrast, the amphibole megacrysts are isotopically distinct (¹⁷⁶Hf/¹⁷⁷Hf_{(34 Ma)} = 0.282830 ± 0.000011). Rare earth element D values for garnet megacryst-nephelinite melt and ∼1 mm amphibole-nephelinite melt plotted as a function of ionic radii show classic near-parabolic trends that are in excellent agreement with crystal lattice-strain models. These observations are consistent with equilibrium between the whole-rock nephelinite, the ∼1 mm amphibole grains within the nephelinite and the garnet megacrysts.High-precision isotope dilution results for Zr and Hf in garnet (D_Zr = 0.220 ± 0.007 and D_Hf = 0.216 ± 0.005 [2σ]), and for all HFSE in amphibole are consistent with previous experimental findings. However, our measurements for Nb and Ta in garnet (D_Nb = 0.0007 ± 0.0001 and D_Ta = 0.0011 ± 0.0006 [2σ]) show that conventional methods may overestimate Nb and Ta concentrations, thereby overestimating both Nb and Ta absolute D values for garnet by up to 3 orders of magnitude and underestimating D_Nb/D_Ta by greater than a factor of 100. As a consequence, the role of residual garnet in imposing Nb/Ta fractionation may be less important than previously thought. Moreover, garnet D_Hf/D_W = 17 and D_Nb/D_Zr = 0.003 imply fractionation of Hf from W and Nb from Zr upon garnet crystallization, which may have influenced short-lived ¹⁸²Hf-¹⁸²W and ⁹²Nb-⁹²Zr isotopic systems in Hadean time.     

Amphibole and phlogopite in “hybrid” metasomatic bands monitor trace element transfer at the interface between felsic and ultramafic rocks (Eastern Alps,Italy)

Ultramafic rocks in contact to gneisses in the Mt. Hochwart HP mélange (Eastern Italian Alps) preserve a series of metasomatic mineral zones. A phlogopitite with minor tremolite and accessory zircon and apatite forms close to the gneiss (Zone 1). Zone 2 consists of tremolite, phlogopite and anthophyllite followed by Zone 3 with anthophyllite plus minor chlorite and talc. Zone 3 grades into an amphibole–garnet peridotite lens. This reaction zone has been generated by infiltration of hydrous fluids at T of 660–700 °C and P < 1.2 GPa, which occurred during exhumation of coupled continental crust and mantle peridotites.The reaction zone between a trace element-rich (gneiss) and a trace element depleted reservoir (peridotite) allows assessment of local trace element mobility in aqueous fluids. We present the results of in situ LA-ICP-MS trace element analysis of minerals from the reaction zone. Phlogopite is the main host for Large Ion Lithophile Elements (LILE) and contributes significantly to the Li, Ti, Nb, Ta, Pb and Sc budget. Anthophyllite is the main host for Li whereas all other trace elements including Rare Earth Elements (REE) are preferentially incorporated into tremolite. Combined with the abundance of these minerals over the contact zone, the mineral trace element data suggests that the LILE and REE were mobile on a small scale of a few centimetres only. Limited mobility of Ta, which is generally regarded as barely mobile in fluids, is documented in elevated contents of Ta in anthophyllite coupled with low Nb/Ta. The high Li content in minerals throughout the reaction zone suggests that Li was the most mobile element.The studied metasomatic zones mirror geochemical processes occurring in subduction zones at the slab–mantle interface. Phlogopite crystallization at the slab–mantle interface is an efficient mechanism to filter LILE from the aqueous fluid. Thus, such reaction zones, forming at temperatures < 660–700 °C, likely prevents that the typical slab signature with enriched LILE is transported by aqueous fluids over long distances in the mantle wedge. However, if coupled to the downgoing slab, phlogopite- and tremolite-rich rocks from such reaction zones might be able to act as carriers of trace elements and water into deeper parts of the subduction zone.     

Fingerprinting a multistage metamorphic fluid–rock history: Evidence from grain scale Sr,O and C isotopic and trace element variations in high-grade marbles from East Antarctica

Granulite grade marble layers interlayered with metapelitic granulites from Lützow Holm Bay, East Antarctica, provide insight into fluid–rock interactions during burial to and exhumation from lower crustal levels. Sub-millimeter scale strontium, oxygen and carbon isotope variations along with LA-ICPMS trace element geochemistry and mineral chemistry of texturally characterized carbonates and associated minerals helped to reconstruct the multistage metamorphic fluid history.Fluid–rock interaction dating back to prograde metamorphism are still preserved in consistently low oxygen and high strontium isotope compositions (δ¹⁸O = 12‰; ⁸⁷Sr/⁸⁶Sr_(550Ma) = 0.7248) within a massif dolomitic marble layer that escaped significant later metasomatism. In most marbles, total re-crystallization and isotopic resetting occurred in the presence of “externally derived” hyper-saline fluids that circulated along the carbonate layers during the early stages of prograde metamorphism. This leads to a trend of increased radiogenic Sr in marbles towards the value of associated metapelitic rocks that have ⁸⁷Sr/⁸⁶Sr_(550Ma) of 0.764.LA-ICPMS studies on trace elements in carbonate and associated silicate minerals at different textural settings, distinguished using cathodoluminescence microscopy, revealed multiple metasomatic events during retrograde metamorphism. Trace element contents of Ba, Sr, Pb and U gave compelling evidence for metasomatic alteration that postdate the exsolution of carbonate at ~ 600 ºC, which can be correlated with the fluids released from the crystallization of anatectic melts and pegmatites. Subsequently, meteoric fluid infiltration occurred at a shallower level of the crust and caused extreme oxygen isotopic heterogeneity (δ¹⁸O = 14.7 ~ ? 4.9‰) and imprinted high concentration of fluid mobile elements. Taken together our results emphasize the importance of integrating textural and chemical heterogeneities to reveal the multiple episodes of fluid–rock interaction processes in a dynamic continental crust, which has major implications on migration of fluids and material and help in formulating models on the geodynamic evolution of crust.     

Robust methods for the decomposition and interpretation of compound dunes applied to a complex hydromorphological setting

Underwater dunes are a morphological feature that are explored by marine scientists and coastal engineers alike. This study presents new methodologies in order to simplify bedform identification and morphodynamic analyses. Specifically, subaqueous compound dunes are decomposed with a simple yet extensive tracking algorithm, which relies on a repeated evaluation of unfiltered bed elevation profiles according to five predefined length classes. In a second step, morphological trends are assessed in the form of bed migration rates, bed slope asymmetries and net sediment changes, in which all parameters are referred to equidistant sections of the examined fairway stretch. This integrated approach not only avoids the challenges in weighting the varying size and abundance of dunes of different scales but also ensures comparability between dune-specific and areal parameters, which significantly improves the interpretation of the morphological setting as a whole. The developed methods are applied to the Outer Jade fairway, an anthropogenically influenced and regularly maintained waterway in the German Bight, and allow scrutiny of spatio-temporal trends in this region. Based on a unique data set of 100 sequential high-quality echo-sounding surveys, various types of bedforms are identified, comprising large-scale primary as well as superimposing secondary dunes that are assumed to interfere with each other. Temporal trends show a long-term rise of the troughs of major bedforms and constant maximum crest elevations near the official maintenance depth, which matches the observed long-term aggradation of sediments. The spatial distribution of integrated morphodynamic parameters reflects a previously described zone of primary dune convergence and facilitates the precise localization of this geophysical singularity. The presented findings both confirm the robustness of the proposed methodologies and, in return, enhance the understanding of morphological processes in the Outer Jade. © 2020 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

Precision spacecraft navigation using a low-cost GPS receiver

Within the PROBA-2 microsatellite mission, a miniaturized single-frequency GPS receiver based on commercial-off-the-shelf (COTS) technology is employed for onboard navigation and timing. A rapid electronic fuse protects against destructive single-event latch-ups (SEL) and enables a quasi-continuous receiver operation despite the inherent sensitivity to space radiation. While limited to single-frequency C/A-code tracking with a narrow-band frontend, the receiver is able to provide precision navigation services through processing of raw GPS measurements on ground as well as a built-in real-time navigation system. In both cases, ionospheric path delays are eliminated through a combination of L1 pseudorange and carrier phase measurements, which also offers a factor-of-two noise reduction relative to code-only processing. By comparison with satellite laser ranging (SLR) measurements, a 0.3-m (3D rms) accuracy is demonstrated for the PROBA-2 reduced dynamic orbit determinations using post-processed GPS orbit and clock products. Furthermore, the experimental onboard navigation system is shown to provide real-time position information with a 3D rms accuracy of about 1?m, which notably outperforms the specification of the Standard Positioning Service (SPS). In view of their lower hardware complexity, mass budget and power requirements as well as the reduced interference susceptibility, legacy C/A-code receivers can thus provide an attractive alternative to dual-frequency receivers even for demanding navigation applications in low Earth orbit.     

The 2008 DGFI realization of the ITRS: DTRF2008

A new realization of the International Terrestrial System was computed at the ITRS Combination Centre at DGFI as a contribution to ITRF2008. The solution is labelled DTRF2008. In the same way as in the DGFI computation for ITRF2005 it is based on either normal equation systems or estimated parameters derived from VLBI, SLR, GPS and DORIS observations by weekly or session-wise processing. The parameter space of the ITRS realization comprises station positions and velocities and daily resolved Earth Orientation Parameters (EOP), whereby for the first time also nutation parameters are included. The advantage of starting from time series of input data is that the temporal behaviour of geophysical parameters can be investigated to decide whether the parameters can contribute to the datum realization of the ITRF. In the same way, a standardized analysis of station position time series can be performed to detect and remove discontinuities. The advantage of including EOP in the ITRS realization is twofold: (1) the combination of the coordinates of the terrestrial pole—estimated from all contributing techniques—links the technique networks in two components of the orientation, leading to an improvement of consistency of the Terrestrial Reference Frame (TRF) and (2) in their capacity as parameters common to all techniques, the terrestrial pole coordinates enhance the selection of local ties as they provide a measure for the consistency of the combined frame. The computation strategy of DGFI is based on the combination of normal equation systems while at the ITRS Combination Centre at IGN solutions are combined. The two independent ITRS realizations provide the possibility to assess the accuracy of ITRF by comparison of the two frames. The accuracy evaluation was done separately for the datum parameters (origin, orientation and scale) and the network geometry. The accuracy of the datum parameters, assessed from the comparison of DTRF2008 and ITRF2008, is between 2–5?mm and 0.1–0.8?mm/year depending on the technique. The network geometry (station positions and velocities) agrees within 3.2?mm and 1.0?mm/year. A comparison of DTRF2008 and ITRF2005 provides similar results for the datum parameters, but there are larger differences for the network geometry. The internal accuracy of DTRF2008—that means the level of conservation of datum information and network geometry within the combination—was derived from comparisons with the technique-only multi-year solutions. From this an internal accuracy of 0.32?mm for the VLBI up to 3.3?mm for the DORIS part of the network is found. The internal accuracy of velocities ranges from 0.05?mm/year for VLBI to 0.83?mm/year for DORIS. The internal consistency of DTRF2008 for orientation can be derived from the analysis of the terrestrial pole coordinates. It is estimated at 1.5–2.5?mm for the GPS, VLBI and SLR parts of the network. The consistency of these three and the DORIS network part is within 6.5?mm.