One Million Years tephra record at IODP Sites U1436 and U1437: Insights into explosive volcanism from the Japan and Izu arcs

The 1 Myr tephra records of IODP (International Ocean Discovery Program) Holes U1436A and U1437B in the Izu‐Bonin fore‐ and reararc were investigated in order to assess provenance and eruptive volumes, respectively. In total, 304 tephra samples were examined and 260 primary tephra layers were identified. Tephra provenance was determined by means of major and trace element compositions of glass shards and distinguished between Japan and Izu‐Bonin arc origin of the tephra layers. A total of 33 marine tephra compositions were correlated to the Japan arc and 227 to the Izu arc. Twenty marine tephra layers were correlated between the two drilling sites. Additionally, we defined eleven correlations of marine tephra deposits to major widespread Japanese eruptions; from the 1.05 Ma Shishimuta‐Pink Tephra to the 30 ka Aira‐Tn Tephra, both from Kyushu Island. These eruptions provide independent time markers within the sediment record and six correlations were used to date tephra layers from Japan in Hole U1436A to establish an alternative age model for this hole. Furthermore, the minimum distal tephra volumes of all detected events were calculated, which enabled the comparison of the tephra volumes that derived from the Japan and the Izu‐Bonin arcs. For some of the major Japanese eruptions these are the first volume estimations that also include distal deposits. All of the Japanese tephras derived from events with eruption magnitude M_v ≥ 5.6 and three of the investigated eruptions reach magnitudes M_v ≥ 7. Volcanic events of the Izu‐Bonin arc have mostly eruption magnitudes M_v ≤ 5.     

Oxidized and reduced mineral assemblages in greenstone belt rocks of the St. Ives gold camp, Western Australia: vectors to high-grade ore bodies in Archaean gold deposits?

Hydrothermal sulfide–oxide–gold mineral assemblages in gold deposits in the Archaean St. Ives gold camp in Western Australia indicate extremely variable redox conditions during hydrothermal alteration and gold mineralization in space and time. Reduced alteration assemblages (pyrrhotite–pyrite) occur in deposits in the southwest of the camp (e.g., Argo, Junction deposits) and moderately to strongly oxidized assemblages (magnetite–pyrite, hematite–pyrite) occur in deposits in the Central Corridor in the northeast (e.g., North Orchin, Revenge deposits). Reduced mineral assemblages flank the Central Corridor of oxidized deposits and, locally, cut across it along E–W trending faults. Oxidized mineral assemblages in the Central Corridor are focused on gravity lows which are interpreted to reflect abundant felsic porphyritic intrusions at about 1,000 m below present surface. Hydrothermal magnetite predates and is synchronous with early phases of gold-associated albite–carbonate–pyrite–biotite–chlorite hydrothermal alteration. Later-stage, gold-associated pyrite is in equilibrium with hematite. The spatial distribution and temporal sequence of iron sulfides and oxides with gold indicate the presence of at least two spatially restricted but broadly synchronous hydrothermal fluids with contrasting redox states. Sulfur isotope constraints support the argument that the different mineral assemblages reflect differences in redox conditions. The δ ³⁴S values for pyrite for the St. Ives gold camp range between −8.4‰ and +5.1‰ with the negative values occurring in oxidized magnetite-rich domains and slightly negative or positive values occurring in reduced, pyrrhotitic domains. Preliminary spatial and paragenetic analysis of the distribution of iron sulfides and oxides in the St. Ives camp suggests that gold grades are highest where the redox state of the hydrothermal alteration assemblages switches from relatively reduced pyrrhotite–pyrite to relatively oxidized magnetite–pyrite and hematite–pyrite both in space and time. Gold deposition is inferred to have occurred where fluids of contrasting redox state mixed.     

A priori fully populated covariance matrices in least-squares adjustment—case study: GPS relative positioning

In this contribution, using the example of the Mátern covariance matrices, we study systematically the effect of apriori fully populated variance covariance matrices (VCM) in the Gauss–Markov model, by varying both the smoothness and the correlation length of the covariance function. Based on simulations where we consider a GPS relative positioning scenario with double differences, the true VCM is exactly known. Thus, an accurate study of parameters deviations with respect to the correlation structure is possible. By means of the mean-square error difference of the estimates obtained with the correct and the assumed VCM, the loss of efficiency when the correlation structure is missspecified is considered. The bias of the variance of unit weight is moreover analysed. By acting independently on the correlation length, the smoothness, the batch length, the noise level, or the design matrix, simulations allow to draw conclusions on the influence of these different factors on the least-squares results. Thanks to an adapted version of the Kermarrec–Schön model, fully populated VCM for GPS phase observations are computed where different correlation factors are resumed in a global covariance model with an elevation dependent weighting. Based on the data of the EPN network, two studies for different baseline lengths validate the conclusions of the simulations on the influence of the Mátern covariance parameters. A precise insight into the impact of apriori correlation structures when the VCM is entirely unknown highlights that both the correlation length and the smoothness defined in the Mátern model are important to get a lower loss of efficiency as well as a better estimation of the variance of unit weight. Consecutively, correlations, if present, should not be neglected for accurate test statistics. Therefore, a proposal is made to determine a mean value of the correlation structure based on a rough estimation of the Mátern parameters via maximum likelihood estimation for some chosen time series of observations. Variations around these mean values show to have little impact on the least-squares results. At the estimates level, the effect of varying the parameters of the fully populated VCM around these approximated values was confirmed to be nearly negligible (i.e. a mm level for strong correlations and a submm level otherwise).     

Deriving incline values for street networks from voluntarily collected GPS traces

ABSTRACT

When producing optimal routes through an environment, considering the incline of surfaces can be of great benefit in a number of use cases. For instance, steep segments need to be avoided for energy-efficient routes and for routes that are suitable for mobility-restricted people. Such incline information may be derived from digital elevation models (DEMs). However, the corresponding data capturing methods (e.g. airborne LiDAR, photogrammetry, and terrestrial surveying) are expensive. Current low-cost and open-licensed DEM (e.g. Shuttle Radar Topography Mission [SRTM] and Advanced Spaceborne Thermal Emission and Reflection Radiometer [ASTER]) generally do not have sufficient horizontal resolution or vertical accuracy, and lack a global coverage. Therefore, we have investigated an alternative low-cost approach which derives street incline values from GPS traces that have been voluntarily collected by the OpenStreetMap contributors. Despite the poor absolute accuracy of this data, the relative accuracy of traces seems to be sufficient enough to compute incline values with reasonable accuracy. A validation shows that the accuracy of incline values calculated from GPS traces slightly outperforms incline values derived from SRTM-1 DEM, though results depend on how many traces per street segment are used for computation.     

Atmospheric Mercury Concentrations from Several Observatory Sites in the Northern Hemisphere

In an effort to investigate both large-scale (spatial) and short/long-term (temporal) distribution characteristics of atmospheric mercury, we have combined and analyzed the Hg concentration data sets collected continuously by four different scientific groups for the areas and periods covering (1) America (three sites near the Canadian Great Lakes (CGL): 1997–2000), (2) Asia (Seoul, Korea (SEL): 1997–2002), (3) Arctic (Alert, Canada (ALT): 1995–2001), and (4) Europe (Mace Head, Ireland (MH): 1996–2002). The mean concentrations of Hg data from those widely dispersed monitoring stations were computed to be (1) 1.58 ± 0.23, 1.69 ± 0.32, and 1.93 ± 0.44 (three sites in CGL), (2) 5.06 ± 2.46 ng m^–3 (SEL), (3) 1.55 ± 0.41 (ALT), and (4) 1.76 ± 0.31 (MH). Intersite relationships were investigated among all different stations using the data groups divided into different temporal intervals. The analysis of diurnal variation patterns of Hg indicated differences in regional source/sink characteristics, with increasing amplitudes of variability toward areas under the strong influence of anthropogenic sources. When the analysis was made over different seasons, the patterns contrasted greatly between the Arctic and the other areas. It was found that the relative enhancement of Hg concentrations was dominant during winter/spring in most areas due to direct or indirect influences of anthropogenic emissions. However, the pattern for the Arctic area was distinguished pronouncedly from others with the spring minimum and summer maximum both of which reflect the potent effects of mercury depletion phenomenon (MDP). By contrast, no long-term trend, either being an increase or decrease, was evident from any of the stations during each respective study period. Although our initial attempts to examine the distribution characteristics of Hg analyzed by different scientific groups were successful, we feel that these efforts should be continued further to extend the compatibility of the global database of Hg.     

Holocene environmental history in high‐Arctic North Greenland revealed by a combined biomarker and macrofossil approach

In this study, we use a combined biomarker and macrofossil approach to reconstruct the Holocene climate history recorded in Trifna Sø, Skallingen area, eastern North Greenland. Chronological information is derived from comparison of lithological, biogeochemical and macrofossil characteristics with a well‐dated record from nearby Lille Sneha Sø. Following local deglaciation around c. 8 cal. ka BP, the local peak warmth occurred between c. 7.4 and 6.2 cal. ka BP as indicated by maximum macrofossil abundances of warmth‐demanding plants (Salix arctica and Dryas integrifolia) and invertebrates (Daphnia pulex and Chironomidae). Warm conditions were dominated by terrestrial organic matter (OM) sedimentation as implied by the alkane‐based P_aq ratio, but increased aquatic productivity is indicated when temperature was highest around 6.5 cal. ka BP. The n‐C₂₉/n‐C₃₁ alkane ratio shows that vegetation in the catchment was dominated by shrubs after deglaciation, but shifted towards relatively more grassy/herbaceous vegetation during peak warmth. After 5.4 cal. ka BP, the disappearance of warmth‐demanding plant and invertebrate macrofossils indicates cooling in the Skallingen area. This cooling was characterized by a significant shift towards dominance of aquatic OM sedimentation in Trifna Sø as implied by high P_aq ratios. Cooling was also associated with a shift in vegetation type from dwarf‐shrub heaths towards relatively more herbaceous vegetation in the catchment, stronger erosion and more oligotrophic conditions in the lake. Our data show that mean air temperatures inferred using branched glycerol dialkyl glycerol tetraethers (brGDGTs) do not seem to accurately reflect the local climatic history. Irrespective of calibration, methylation of branched tetraethers (MBT) palaeothermometry cannot be reconciled with the macrofossil evidence and seems to be biased by either changing brGDGT sources (in situ vs. soil‐derived) or changing species assemblages and/or an unknown physiological response to changing environmental conditions at high latitude.     

Use of Picard and Newton iteration for solving nonlinear ground water flow equations

This study examines the use of Picard and Newton iteration to solve the nonlinear, saturated ground water flow equation. Here, a simple three-node problem is used to demonstrate the convergence difficulties that can arise when solving the nonlinear, saturated ground water flow equation in both homogeneous and heterogeneous systems with and without nonlinear boundary conditions. For these cases, the characteristic types of convergence patterns are examined. Viewing these convergence patterns as orbits of an attractor in a dynamical system provides further insight. It is shown that the nonlinearity that arises from nonlinear head-dependent boundary conditions can cause more convergence difficulties than the nonlinearity that arises from flow in an unconfined aquifer. Furthermore, the effects of damping on both convergence and convergence rate are investigated. It is shown that no single strategy is effective for all problems and how understanding pitfalls and merits of several methods can be helpful in overcoming convergence difficulties. Results show that Picard iterations can be a simple and effective method for the solution of nonlinear, saturated ground water flow problems.