Refined discrete and empirical horizontal gradients in VLBI analysis

Missing or incorrect consideration of azimuthal asymmetry of troposphere delays is a considerable error source in space geodetic techniques such as Global Navigation Satellite Systems (GNSS) or Very Long Baseline Interferometry (VLBI). So-called horizontal troposphere gradients are generally utilized for modeling such azimuthal variations and are particularly required for observations at low elevation angles. Apart from estimating the gradients within the data analysis, which has become common practice in space geodetic techniques, there is also the possibility to determine the gradients beforehand from different data sources than the actual observations. Using ray-tracing through Numerical Weather Models (NWMs), we determined discrete gradient values referred to as GRAD for VLBI observations, based on the standard gradient model by Chen and Herring (J Geophys Res 102(B9):20489–20502, 1997.  https://doi.org/10.1029/97JB01739) and also for new, higher-order gradient models. These gradients are produced on the same data basis as the Vienna Mapping Functions 3 (VMF3) (Landskron and Böhm in J Geod, 2017.  https://doi.org/10.1007/s00190-017-1066-2), so they can also be regarded as the VMF3 gradients as they are fully consistent with each other. From VLBI analyses of the Vienna VLBI and Satellite Software (VieVS), it becomes evident that baseline length repeatabilities (BLRs) are improved on average by 5% when using a priori gradients GRAD instead of estimating the gradients. The reason for this improvement is that the gradient estimation yields poor results for VLBI sessions with a small number of observations, while the GRAD a priori gradients are unaffected from this. We also developed a new empirical gradient model applicable for any time and location on Earth, which is included in the Global Pressure and Temperature 3 (GPT3) model. Although being able to describe only the systematic component of azimuthal asymmetry and no short-term variations at all, even these empirical a priori gradients slightly reduce (improve) the BLRs with respect to the estimation of gradients. In general, this paper addresses that a priori horizontal gradients are actually more important for VLBI analysis than previously assumed, as particularly the discrete model GRAD as well as the empirical model GPT3 are indeed able to refine and improve the results.     

Constraints on the source of Cu in a submarine magmatic-hydrothermal system,Brothers volcano,Kermadec island arc

Most magmatic-hydrothermal Cu deposits are genetically linked to arc magmas. However, most continental or oceanic arc magmas are barren, and hence new methods have to be developed to distinguish between barren and mineralised arc systems. Source composition, melting conditions, the timing of S saturation and an initial chalcophile element-enrichment represent important parameters that control the potential of a subduction setting to host an economically valuable deposit. Brothers volcano in the Kermadec island arc is one of the best-studied examples of arc-related submarine magmatic-hydrothermal activity. This study, for the first time, compares the chemical and mineralogical composition of the Brothers seafloor massive sulphides and the associated dacitic to rhyolitic lavas that host the hydrothermal system. Incompatible trace element ratios, such as La/Sm and Ce/Pb, indicate that the basaltic melts from L’Esperance volcano may represent a parental analogue to the more evolved Brothers lavas. Copper-rich magmatic sulphides (Cu?>?2 wt%) identified in fresh volcanic glass and phenocryst phases, such as clinopyroxene, plagioclase and Fe–Ti oxide suggest that the surrounding lavas that host the Brothers hydrothermal system represent a potential Cu source for the sulphide ores at the seafloor. Thermodynamic calculations reveal that the Brothers melts reached volatile saturation during their evolution. Melt inclusion data and the occurrence of sulphides along vesicle margins indicate that an exsolving volatile phase extracted Cu from the silicate melt and probably contributed it to the overlying hydrothermal system. Hence, the formation of the Cu-rich seafloor massive sulphides (up to 35.6 wt%) is probably due to the contribution of Cu from a bimodal source including wall rock leaching and magmatic degassing, in a mineralisation style that is hybrid between Cyprus-type volcanic-hosted massive sulphide and subaerial epithermal–porphyry deposits.     

Visibility in airborne volcanic ash: considerations for surface transportation using a laboratory-based method

All modes of surface transportation can be disrupted by visibility degradation caused by airborne volcanic ash. Despite much qualitative evidence of low visibility on roads following historical eruptions worldwide, there have been few detailed studies that have attempted to quantify relationships between visibility conditions and observed impacts on network functionality and safety. In the absence of detailed field observations, such gaps in knowledge can be filled by developing empirical datasets through laboratory investigations. Here, we use historical eruption data to estimate a plausible range of ash-settling rates and ash particle characteristics for Auckland city, New Zealand. We propose and implement a new experimental set-up in controlled laboratory conditions, which incorporates a dual-pass transmissometer and solid aerosol generator, to reproduce these ash-settling rates and calculate visual ranges through the associated airborne volcanic ash. Our findings demonstrate that visibility is most impaired for high ash-settling rates (i.e. > 500 g m^?2 h^?1) and particle size is deemed the most influential ash characteristic for visual range. For the samples tested (all < 320 μm particle diameter), visibility was restricted to ~ 1–2 m when ash settling was replicated for very high rates (i.e. ~ 4000 g m^?2 h^?1) and was especially low when ash particles were fine-grained, more irregular in shape and lighter in colour. Finally, we consider potential implications for disruption to surface transportation in Auckland through comparisons with existing research which investigates the consequences of visual range reduction for other atmospheric hazards such as fog. This includes discussing how our approach might be utilised in emergency and transport management planning. Finally, we summarise strategies available for the mitigation of visibility degradation in environments contaminated with volcanic ash.     

Pennsylvanian glacimarine sedimentation in the Cushamen Formation,western North Patagonian Massif

The metasedimentary sequence of the Cushamen Formation in the western North Patagonian Massif is540 m thick and comprises six sedimentary lithofacies associations related to a glacimarine environment.Four of these lithofacies represent distal glacimarine environments, whereas another one was deposited in proximal glacimarine environments, and the last includes subglacial environments. The organization and configuration of these lithofacies associations represent the advance and retreat of the glacier masses. The maximum glacial advance is correlatable with the G2 glacial interval of the Pennsylvanian Pampa de Tepuel, Las Salinas and Valley Chico, formations of the Extraandean Chubut, and the southern part of Neuquen Cordillera. Contemporaneously, in southern Chile there are marine and glacimarine sediments. The chronostratigraphic relationships between the Silurian to Permian units allow five paleogeographic stages to be distinguished. The middle Silurian-late Devonian igneous rocks represent the first magmatic stage. The second stage, which is transitional to the first, is represented by a marine basin that includes the late Devonian-early Carboniferous Esquel and Rio Pescado formations and the Llanquihue Complex. The third stage(early-late Carboniferous) includes granitoids of the second magmatic event that partially overlapped the first magmatic igneous belt. The fourth stage belongs to the late Carboniferous sedimentation of the Cushamen and equivalent formations. The extended early Permian magmatism was the last Paleozoic event in the studied area.     

3-D crustal-scale gravity model of the San Rafael Block and Payenia volcanic province in Mendoza,Argentina

The San Rafael Block(SRB)is part of one of the main retroarc volcanic provinces in southern Central Andes in Mendoza,Argentina.This block is located in the Andean foothills between the orogenic front and foreland basement uplifts of late Miocene age.In order to analyze the geochronological evolution of the Quaternary volcanism in the region,several geologic and geophysical studies have been conducted.Nevertheless,the crust,where the SRB is located,has not been well characterized yet.Based on gravimetric and magnetic data,together with isostatic and elastic thickness analyses,we modeled the crustal structure of the area.Information obtained has allowed us to understand the crust where the SRB and the Payenia volcanic province are located.Bouguer anomalies indicate that the SRB presents higher densities to the North of Cerro Nevado and Moho calculations suggest depths for this block between 40 and 50 km.Determinations of elastic thickness would indicate that the crust supporting the San Rafael Block presents values of approximately 10 km,being enough to support the block loading.However,in the Payenia region,elastic thickness values are close to zero due to the regional temperature increase.     

CyberConnector: a service-oriented system for automatically tailoring multisource Earth observation data to feed Earth science models

Feeding multisource Earth observation (EO) data into Earth science models (ESM) remains a daunting challenge. This paper presents a service-oriented approach as an alternative solution. It uses geospatial web services to process the EO data and geoprocessing workflow for automation. Different from existing approaches, it takes advantage of virtual data products (VDP) to release modelers from intensive data processing. It can directly connect ESMs to public EO sources via Cyberinfrastructure. A prototype called CyberConnector is implemented. CyberConnector supports intuitive building of VDP, automatic execution of workflows and effortless retrieval of model-ready input files. We used it to stream multiple datasets to several ESMs including finite-volume coastal ocean model (FVCOM) and cloud-resolving model (CRM). The results show that CyberConnector can truly benefit modelers on time saving and effort minimizing.     

Converting scanned images of seismic reflection data into SEG-Y format

Archives across the world contain vast amounts of old or “vintage” seismic reflection data, which are largely inaccessible for geo-scientific research, due to the out-dated media on which they are stored. Despite the age of these data, they often have great potential to be of use in modern day research. It is often the case that seismic reflection data within these archives are only available as a processed stacked section, printed on paper or film. In this study, a method for the conversion (vectorization) of scanned images of stacked reflection seismic data to standard SEG-Y format is presented. The method addresses data displayed with a line denoting the waveform, where areas on one side of the baseline are shaded (i.e. wiggle trace, variable fill). The method provides an improvement on other published methods utilized within currently available academic software. Unlike previous studies, the method used to detect trace baselines and to detect and remove timelines on the seismic image is described in detail. Furthermore, a quantitative analysis of the performance of the method is presented, showing that an average trace-to-trace correlation coefficient of between 0.8 and 0.95 can be achieved for typical plotting styles. Finally, a case study where the method is applied to vectorize over 1700 km of land seismic data from the island of Gotland (Sweden) is presented.     

Standardization of a NaI(Tl) scintillation detector and its use on environments for Geology teaching in a Brazilian university

This paper describes analytical protocols adopted for calibrating a portable gamma rays spectrometer using a 2? × 2? NaI(Tl) scintillation detector for performing exploration and environmental surveys. The steps involved the determination of counting times with good accuracy, calibration in energy and concentration, and data conversion to effective dose (Ed) for monitoring the radioactivity levels in the air, for indoor and outdoor environments. The χ² test was applied to the counting rate data acquired for ten different readings, allowing determine the greater accuracy among them. Some standards yielded values above of the detection limit for calibrating the K, eU, and eTh windows, whereas empirical equations provided satisfactory results. The count rates in these windows were obtained in selected environments related to Geology teaching activities in a Brazilian university. Then, they were converted to Ed, showing the applicability of the adopted procedures for environmental surveys. The radiological monitoring showed two sites exhibiting Ed values above the WHO guideline reference value of 1 mSv/year. The first was in a minerals and rocks museum and the second in a deposit (litoteca) dedicated to the rocks storage. Thus, the members of the public are only occasionally exposed to the natural gamma radiation there. Additionally, the Ed values are much lower than the recommended limit of 20 mSv/year for occupationally exposed individuals.     

Robust system size reduction of discrete fracture networks: a multi-fidelity method that preserves transport characteristics

We propose a multi-fidelity system reduction technique that uses weighted graphs paired with three-dimensional discrete fracture network (DFN) modelling for efficient simulation of subsurface flow and transport in fractured media. DFN models are used to simulate flow and transport in subsurface fractured rock with low-permeability. One method to alleviate the heavy computational overhead associated with these simulations is to reduce the size of the DFN using a graph representation of it to identify the primary flow sub-network and only simulate flow and transport thereon. The first of these methods used unweighted graphs constructed solely on DFN topology and could be used for accurate predictions of first-passage times. However, these techniques perform poorly when predicting later stages of the mass breakthrough. We utilize a weighted-graph representation of the DFN where edge weights are based on hydrological parameters in the DFN that allows us to exploit the kinematic quantities derivable a posteriori from the flow solution obtained on the graph representation of the DFN to perform system reduction and predict the later stages of the breakthrough curve with high fidelity. We also propose and demonstrate the use of an adaptive pruning algorithm with error control that produces a pruned DFN sub-network whose predicted mass breakthrough agrees with the original DFN within a user-specified tolerance. The method allows for the level of accuracy to be a user-controlled parameter.     

Rock Glaciers in the Patagonian Andes: An Inventory for the Monte San Lorenzo (Cerro Cochrane) Massif, 47° S

Although rock glaciers in the Central and Desert Andes of Argentina and Chile have been previously studied in detail, much less attention has been paid to the occurrence of these permafrost forms in Patagonia. Recently, however, the establishment of the Argentinean Glacier Inventory program, which intends to inventory and monitor all ice masses along the Argentinean Andes, has started a large amount of new geocryological research. The project is designed to provide reliable and worldwide comparable results, supported by well established technical procedures and background information. Presented here is the first rock glacier inventory of the Monte San Lorenzo (Cerro Cochrane) region in the southern Patagonian Andes. A total of 130 intact (9.86 km²) and 47 fossil (1.45 km²) landforms were inventoried using two 2.5 m resolution ALOS Panchromatic Remote‐sensing Instruments for Stereo Mapping images. Since the Argentinean federal initiative described above legally protects all rock glaciers in the country as water reserves, and due to the little scientific knowledge concerning rock glaciers in the vast majority of the Patagonian Andes, this inventory provides an important basis for political decision‐making and opens further geocryological research avenues for the Patagonian region in general.