Earthen barriers to control lava flows in the 2001 eruption of Mt. Etna

Preceded by four days of intense seismicity and marked ground deformation, a new eruption of Mt. Etna started on 17 July and lasted until 9 August 2001. It produced lava emission and strombolian and phreatomagmatic activity from four different main vents located on a complex fracture system extending from the southeast summit cone for about 4.5 km southwards, from 3000 to 2100 m elevation (a.s.l.). The lava emitted from the lowest vent cut up an important road on the volcano and destroyed other rural roads and a few isolated country houses. Its front descended southwards to about 4 km distance from the villages of Nicolosi and Belpasso. A plan of intervention, including diversion and retaining barriers and possibly lava flow interruption, was prepared but not activated because the flow front stopped as a consequence of a decrease in the effusion rate. Extensive interventions were carried out in order to protect some important tourist facilities of the Sapienza and Mts. Silvestri zones (1900 m elevation) from being destroyed by the lava emitted from vents located at 2700 m and 2550 m elevation. Thirteen earthen barriers (with a maximum length of 370 m, height of 10–12 m, base width of 15 m and volume of 25 000 m³) were built to divert the lava flow away from the facilities towards a path implying considerably less damage. Most of the barriers were oriented diagonally (110–135°) to the direction of the flow. They were made of loose material excavated nearby and worked very nicely, resisting the thrust of the lava without any difficulty. After the interventions carried out on Mt. Etna in 1983 and in 1991–1992, those of 2001 confirm that earthen barriers can be very effective in controlling lava flows.     

Stratigraphy,structure and geochronology of ore leads in the Matsitama Schist belt of Northern Botswana

The Matsitama schist belt in northeastern Botswana comprises an area of metasediments, notably quartzites, limestones, shales and amphibolites that are bounded by granites and gneisses. The belt lies southwest of the Rhodesian cration and north of the Limpopo mobile belt.Stratigraphic, structural and lead isotopic evidence indicates that the Matsitama metasediments are equivalent to the Shashi metasediments in the Limpopo belt. There is strong evidence that the Matsitama and Shashi metasediments stratigraphically underlie volcanic rocks of the Tati belt which have been correlated with Archaean schist belts of about 2700 Ma of Rhodesia. Therefore, the Matsitama and Shashi rocks are at least as old as the schist belts of the Rhodesian craton and may represent a shallow-water facies that occurs only in the Limpopo area.There is no structural evidence that the Matsitama and Shashi metasediments were deposited unconformably on basement rocks, although the presence of gneiss, amphibolite and ironstone pebbles in a Matsitama conglomerate, as well as the presence of orthoquartzites, shows the existence of a basement source region. However, the surrounding granites intrude the Matsitama and Shashi metasediments and all underwent several deformation phases.The structural history of the Matsitama rocks can be described in terms of five phases of deformation. The main cleavage-producing deformation phase, F₂, folded the rocks into a major synform and intensely deformed them. Before this, however, the rocks had been folded and thrust so that part of the succession shows downward-facing F₂ structures and there are possibly repetitions of the stratigraphy due to imbrication. Structures of the F₃ and F₄ phases fold the main cleavage but locally are sufficiently intense to modify the shape of the finite strain ellipsoid. There is a major ductile shear zone of F₄ age, south of which F₄ folds are tight, while to the north, F₄ deformation is negligible. All of these structures can be correlated with deformation phases in the Tati schist belt to the east and in the northern part of the Limpopo mobile belt.Lead isotope evidence suggests that mineralization in the Matsitama metasediments occurred at least 2200 Ma ago, and that leads from Dihudi/Thakadu and Messina, in the centre of the Limpopo belt, underwent a two-stage history of events at 2600–2700 Ma and 2000–2100 Ma ago, agreeing with other geochronological evidence. The leads from Matsitama and Messina are isotopically distinct from leads from the Rhodesian schist belts, which show evidence of transfer to the crust some 3500 Ma ago. The absence of this 3500 Ma-old lead from the Matsitama and Messina environments may indicate different crustal conditions and possibly the absence of the Rhodesian-type early basement.     

The influence of particle composition on thorium scavenging in the NE Atlantic ocean (POMME experiment)

²³⁰Th, ²³²Th and ²³⁴Th were analyzed in sinking particles collected by moored and drifting sediment traps in the NE Atlantic Ocean (POMME experiment) in order to constrain the phase(s) carrying Th isotopes in the water column. It reveals a contrasted behaviour between ²³⁴Th and ²³⁰Th. ²³⁴Th is correlated to the particulate organic carbon suggesting that it is primarily scavenged by organic compounds in the surface waters. ²³⁰Th_xs is correlated with Mn, Ba and the lithogenic fraction that are enriched in small suspended particles and incorporated in the sinking particulate flux throughout the water column. The lack of correlation between ²³⁰Th_xs and CaCO₃ or biogenic silica (bSi) indicates that CaCO₃ and bSi are not responsible for ²³⁰Th scavenging in the deep waters of this oceanic region. ²³⁰Th is generally correlated with the lithogenic content of the trapped material but this correlation disappears in winter during strong atmospheric dust inputs suggesting that lithogenic matter is not directly responsible for ²³⁰Th scavenging in the deep waters or that sufficient time is required to achieve particle–solution equilibration. MnO₂ could be the prevalent ²³⁰Th_xs-bearing phase. The narrow range of K_{d_MnO2}^Th obtained for very contrasted oceanic environments supports a global control of ²³⁰Th_xs scavenging by MnO₂ and raises the possibility that the ²³⁰Th–²³¹Pa fractionation is controlled by the amount of colloidal MnO₂ in seawater.     

Refining the ice flow chronology and subglacial dynamics across the migrating Labrador Divide of the Laurentide Ice Sheet with age constraints on deglaciation

The Laurentide Ice Sheet was characterized by a dynamic polythermal base. However, important data and knowledge gaps have led to contrasting reconstructions in areas such as the Labrador Ice Divide. In this study, detailed fieldwork was conducted at the southeastern edge of a major landform boundary to resolve the relative ice flow chronology and constrain the evolution of the subglacial dynamics, including the migration and collapse of the Labrador Ice Divide. Surficial mapping and analysis of 94 outcrop‐scale ice flow indicators were used to develop a relative ice flow chronology. ¹⁰Be exposure ages were used with optical ages to confine the timing of deglaciation within the study area. Four phases of ice flow were identified. Flow 1 was a northeasterly ice flow preserved under non‐erosive subglacial conditions associated with the development of an ice divide. Flow 2 was a northwest ice flow, which we correlate to the Ungava Bay Ice Stream and led to a westward migration of the ice divide, preserving Flow 2 features and resulting in Flow 3's eastward‐trending indicators. Flow 4 is limited to sparse fine striations within and around the regional uplands. The new optical ages and ¹⁰Be exposure ages add to the regional geochronology dataset, which further constrains the timing of ice margin retreat in the area to around 8.0 ka. Copyright © 2019 The Authors. Journal of Quaternary Science Published by John Wiley & Sons Ltd.     

Observing changes in groundwater resource using hydro-chemical and isotopic parameters: a case study from Bist Doab,Punjab

The continuous and large-scale abstraction of groundwater has created a groundwater depletion problem in several parts of the Punjab state including Bist Doab, the interfluve region of Beas and Satluj rivers. In the present study, a few important parameters, viz. water level, stable isotope, EC, temperature, groundwater age, that can be used to fingerprint the over-exploitation of groundwater have been examined. It has been observed that with the increase in over-exploitation, the yield of shallow aquifer is progressively getting reduced and as a result forcing the farmers to sink their wells to deeper depths. With abstraction of deeper aquifer, the storage of old groundwater at the deeper aquifer is declining and getting replaced by induced accelerated inflow of young water from the recharge zone and the overlying shallow aquifer. The signatures of the modern water have been observed in the data analyzed for isotopic, hydro-chemical facies, electrical conductivity and temperature of water from deeper aquifer. The study has identified the usefulness of these parameters for identifying groundwater over-exploitation in the region. Depleting water resource may stagnate the economic progress of the region. The paper provides suitable water resource management strategies to be adopted to improve the sustainability of water resources and economic growth in the region.     

Rivers turned to rock: Late Quaternary alluvial induration influencing the behaviour and morphology of an anabranching river in the Australian monsoon tropics

Late Quaternary alluvial induration has greatly influenced contemporary channel morphology on the anabranching Gilbert River in the monsoon tropics of the Gulf of Carpentaria. The Gilbert, one of a number of rivers in this region, has contributed to an extensive system of coalescing low-gradient and partly indurated riverine plains. Extensive channel sands were deposited by enhanced flow conditions during marine oxygen isotope (OI) Stage 5. Subsequent flow declined, probably associated with increased aridity, however, enhanced runoff recurred again in OI Stages 4–3 (65–50 ka). Aridity then capped these plains with 4–7 m of mud. A widespread network of sandy distributary channels was incised into this muddy surface from sometime after the Last Glacial Maximum (LGM) to the mid Holocene during a fluvial episode more active than the present but less so than those of OI Stages 5 and 3. This network is still partly active but with channel avulsion and abandonment now occurring largely proximal to the main Gilbert flow path.A tropical climate and reactive catchment lithology have enhanced chemical weathering and lithification of alluvium along the river resulting in the formation of small rapids, waterfalls and inset gorges, features characteristic more of bedrock than alluvial systems. Thermoluminescence (TL) and comparative optically stimulated luminescence (OSL) ages of the sediments are presented along with U/Th ages of pedogenic calcrete and Fe/Mn oxyhydroxide/ oxide accumulations. They show that calcrete precipitated during the Late Quaternary at times similar to those that favoured ferricrete formation, possibly because of an alternating wet–dry climate. Intense chemical alteration of the alluvium leading to induration appears to have prevailed for much of the Late Quaternary but, probably due to exceptional dryness, not during the LGM. The result has been restricted channel migration and a reduced capacity for the channel to adjust and accommodate sudden changes in bedload. Consequent avulsions have caused local stream powers to increase by an order of magnitude, inducing knickpoint erosion, local incision and the sudden influx of additional bedload that has triggered further avulsions. The Gilbert River, while less energetic than its Pleistocene ancestors, is clearly an avulsive system, and emphasizes the importance in some tropical rivers of alluvial induration for reinforcing the banks, generating nickpoints, reworking sediment and thereby developing and maintaining an indurated and anabranching river style.     

Four-dimensional data assimilation method based on SVD: Theoretical aspect

Summary A new method of four-dimensional data assimilation based on Singular Value Decomposition (SVD) is proposed. In it, a set of atmospheric states is obtained by integrating a numerical weather prediction model and simulated observations are taken and calculated from the model variables. Then the SVD technique is used to create the base vectors from this coupled data set. Finally, the analysis is obtained by projecting actual observation data into a space spanned by the base vectors. Using this approach, the four-dimensional data assimilation becomes a simple linear inverse problem the linearization of the nonlinear forward model is avoided, and the developments of the adjoint and background error covariance matrix are no longer needed. Since the SVD technique is used here, the method is simply called 4DSVD.     

Evaluation of the first GOCE static gravity field models using terrestrial gravity,vertical deflections and EGM2008 quasigeoid heights

Recently, four global geopotential models (GGMs) were computed and released based on the first 2 months of data collected by the Gravity field and steady-state Ocean Circulation Explorer (GOCE) dedicated satellite gravity field mission. Given that GOCE is a technologically complex mission and different processing strategies were applied to real space-collected GOCE data for the first time, evaluation of the new models is an important aspect. As a first assessment strategy, we use terrestrial gravity data over Switzerland and Australia and astrogeodetic vertical deflections over Europe and Australia as ground-truth data sets for GOCE model evaluation. We apply a spectral enhancement method (SEM) to the truncated GOCE GGMs to make their spectral content more comparable with the terrestrial data. The SEM utilises the high-degree bands of EGM2008 and residual terrain model data as a data source to widely bridge the spectral gap between the satellite and terrestrial data. Analysis of root mean square (RMS) errors is carried out as a function of (i) the GOCE GGM expansion degree and (ii) the four different GOCE GGMs. The RMS curves are also compared against those from EGM2008 and GRACE-based GGMs. As a second assessment strategy, we compare global grids of GOCE GGM and EGM2008 quasigeoid heights. In connection with EGM2008 error estimates, this allows location of regions where GOCE is likely to deliver improved knowledge on the Earth’s gravity field. Our ground truth data sets, together with the EGM2008 quasigeoid comparisons, signal clear improvements in the spectral band ~160–165 to ~180–185 in terms of spherical harmonic degrees for the GOCE-based GGMs, fairly independently of the individual GOCE model used. The results from both assessments together provide strong evidence that the first 2 months of GOCE observations improve the knowledge of the Earth’s static gravity field at spatial scales between ~125 and ~110 km, particularly over parts of Asia, Africa, South America and Antarctica, in comparison with the pre-GOCE-era.