The anisotropy of magnetic susceptibility of lava flows: an experimental approach

Measurements of the anisotropy of magnetic susceptibility (AMS) of natural lavas have shown that AMS varies with depth within a lava flow. We have investigated the reasons for such variation by studying the effects of temperature and strain rate on the AMS of recent lava in the laboratory. Samples of lava from Kilauea were melted and subjected to a range of strain rate and cooling histories. The results show that the degree of anisotropy is a function of both the thermal and shearing history of a sample. High degrees of anisotropy were found only in samples that were deformed at temperatures close to those encountered during eruption and then rapidly quenched. Lavas subjected to similar shear stresses at high temperatures had low degrees of anisotropy if allowed to cool down slowly without further deformation. Additionally, lava subjected to complex shearing yield a lower degree of anisotropy even when high strain rates were imposed on it. These results lead to the conclusion that only the last phase of deformation is detectable using AMS and that high strain rates will not result in high degrees of anisotropy if either deformation ends while lava is still fluid or if the orientation of the maximum shear stress varies with time. The relation between the orientation of the principal susceptibilities and that of shear is less sensitive to variation on shear with time. Consequently, flow directions can be inferred confidently with this type of measurements.     

Thick lava flows of Karisimbi Volcano, Rwanda: insights from SIR-C interferometric topography

 We use a digital elevation model (DEM) derived from interferometrically processed SIR-C radar data to estimate the thickness of massive trachyte lava flows on the east flank of Karisimbi Volcano, Rwanda. The flows are as long as 12 km and average 40–60 m (up to >140 m) in thickness. By calculating and subtracting a reference surface from the DEM, we derived a map of flow thickness, which we used to calculate the volume (up to 1 km³ for an individual flow, and 1.8 km³ for all the identified flows) and yield strength of several flows (23–124 kPa). Using the DEM we estimated apparent viscosity based on the spacing of large folds (1.2×10¹² to 5.5×10¹² Pa s for surface viscosity, and 7.5×10¹⁰ to 5.2×10¹¹ Pa s for interior viscosity, for a strain interval of 24 h). We use shaded-relief images of the DEM to map basic flow structures such as channels, shear zones, and surface folds, as well as flow boundaries. The flow thickness map also proves invaluable in mapping flows where flow boundaries are indistinct and poorly expressed in the radar backscatter and shaded-relief images. Received: 6 September 1997 / Accepted: 15 May 1998     

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.     

Multi-instrument remote and in situ observations of the Erebus Volcano (Antarctica) lava lake in 2005: A comparison with the Pele lava lake on the jovian moon Io

The stable, persistent, active lava lake at Erebus volcano (Ross Island, Antarctica) provides an excellent thermal target for analysis of spacecraft observations, and for testing new technology. In the austral summer of 2005 visible and infrared observations of the Erebus lava lake were obtained with sensors on three space vehicles Terra (ASTER, MODIS), Aqua (MODIS) and EO-1 (Hyperion, ALI). Contemporaneous ground-based observations were obtained with hand-held infrared cameras. This allowed a quantitative comparison of the thermal data obtained from different instruments, and of the analytical techniques used to analyze the data, both with and without the constraints imposed by ground-truth. From the thermal camera data, in December 2005 the main Erebus lava lake (Ray Lake) had an area of ≈ 820 m². Surface colour temperatures ranged from 575 K to 1090 K, with a broad peak in the distribution from 730 K to 850 K. Total heat loss was estimated at 23.5 MW. The flux density was ≈ 29 kW m^− 2. Mass flux was estimated at 64 to 93 kg s^− 1. The best correlation between thermal emission and emitting area was obtained with ASTER, which has the best combination of spatial resolution and wavelength coverage, especially in the thermal infrared. The high surface temperature of the lava lake means that Hyperion data are for the most part saturated. Uncertainties, introduced by the need to remove incident sunlight cause the thermal emission from the Hyperion data to be a factor of about two greater than that measured by hand-held thermal camera. MODIS also over-estimated thermal output from the lava lake by the same factor of two because it was detecting reflected sunlight from the rest of the pixel area. The measurement of the detailed temperature distribution on the surface of an active terrestrial lava lake will allow testing of thermal emission models used to interpret remote-sensing data of volcanism on Io, where no such ground-truth exists. Although the Erebus lava lake is four orders of magnitude smaller than the lava lake at Pele on Io, the shape of the integrated thermal emission spectra are similar. Thermal emission from this style of effusive volcanism appears to be invariant. Excess thermal emission in most Pele spectra (compared to Erebus) at short wavelengths (< 3 μm) is most likely due to disruption of the surface on the lava lake by escaping volatiles.     

Comparison of the response of the Landsat 7 Enhanced Thematic Mapper Plus and the Earth Observing-1 Advanced Land Imager over active lava flows

Using near simultaneously acquired Landsat 7 Enhanced Thematic Mapper Plus (ETM+) and Earth Observing-1 Advanced Land Imager (ALI) data we assess the relative radiant responses over active lava flows from the Mt. Etna July/August 2001 flank eruption. By assessing the extent of saturation between the two instruments and using the dual-band method of extracting sub pixel thermal information, we show that the ALI represents an improvement over the ETM+ in the present ability to assess temperatures of hot active lava flows for a number of reasons. (1) The extra spectral channels provided by ALI compliment the current SWIR channels on ETM+ by providing a greater number of paired channel combinations for input into the dual-band method. Thus, dual-band temperature solutions can be determined for a greater range of lava flow types than previously possible using the two paired channel combinations available with the ETM+. (2) The ALI instrument is less susceptible than ETM+ to saturation within the SWIR, especially when using channels 5, 5p and 4p at wavelengths of 1.65, 1.25 and 0.87 μm respectively. (3) The greater radiometric sensitivity of the ALI 12 bit electronics coupled with a significantly higher signal to noise ratio aid in obtaining successful dual-band solutions.     

Post-emplacement tilting of lava flows inferred from magnetic fabric study: the example of Oligocene lavas in the Jeanne d’Arc Peninsula (Kerguelen Islands)

Statistical analysis of the magnetic fabric of samples from several successive lava flows emplaced under similar conditions can allow determination of the mean flow direction when magnetic fabric data from individual flows do not lead to reliable results. A difference between the obtained flow direction and the present dip direction indicates that the flows were tilted after emplacement. For 2 successive series of flows on the Jeanne d’Arc Peninsula presently NNW dipping, this method shows lava emplacement along a SSW–NNE direction. This indicates a gentle tilting acquired during a period of weak deformation in the whole archipelago. Additionally, the magnetic fabric data allow the reconstruction of the different conditions of emplacement of these two series of lava flows and of formation of the local thick conglomerate interbedded between these series.     

Thermal–mechanical modeling of cooling history and fracture development in inflationary basalt lava flows

Thermal–mechanical analyses of isotherms in low-volume basalt flows having a range of aspect ratios agree with inferred isotherm patterns deduced from cooling fracture patterns in field examples on the eastern Snake River Plain, Idaho, and highlight the caveats of analytical models of sheet flow cooling when considering low-volume flows. Our field observations show that low-volume lava flows have low aspect ratios (width divided by thickness), typically < 5. Four fracture types typically develop: column-bounding, column-normal, entablature (all of which are cooling fractures), and inflation fractures. Cooling fractures provide a proxy for isotherms during cooling and produce patterns that are strongly influenced by flow aspect ratio. Inflation fractures are induced by lava pressure-driven inflationary events and introduce a thermal perturbation to the flow interior that is clearly evidenced by fracture patterns around them. Inflation fracture growth occurs incrementally due to blunting of the lower tip within viscoelastic basalt, allowing the inflation fracture to pivot open. The final stage of growth involves propagation beyond the blunted tip towards the stress concentration at the tapered tip of a lava core, resulting in penetration of the core that causes quenching of the lava and the formation of a densely fractured entablature. We present numerical models that include the effects of inflation fractures on lava cooling and which support field-based inferences that inflation fractures depress the isotherms in the vicinity of the fracture, cause a subdivision of the lava core, control the location of the final portion of the lava flow to solidify, and cause significant changes in the local cooling fracture orientations. In addition to perturbing isotherms, inflation fractures cause a lava flow to completely solidify in a shorter amount of time than an identically shaped flow that does not contain an inflation fracture.     

An empirical function to estimate the depths of linear hot sources: Laboratory modeling and field measurements of lava tubes

Estimating depths of buried lava tubes is important for determining the thermal budgets and effusion rates of basaltic volcanic systems. This research used a laboratory experiment scaled to a lava tube system to measure the 3D temperature field surrounding a hot viscous fluid flowing through a buried glass tube while varying conditions such as flow rate and temperature. The depth of the glass tube was changed for different experimental runs. Numerical techniques were applied to model the laboratory experiment. The surface thermal distributions from 166 thermal traverses, constrained to a depth to width ratio of 0.6 to 1.6, were analyzed to empirically derive a depth estimation function using regression techniques. This “Linear Anomaly Surface Transect (LAST)” depth function is a scaleable depth estimation technique which can be solved with thermal imaging data alone. The minimum temperature, maximum temperature and width of a Lorentzian distribution fit to a surface thermal transect, are the only inputs required for the LAST function to estimate the depths of the hot source. The input parameters were then applied to non-laboratory situations including the Kuhio lava tube system in Hawai’i. The LAST function produced depth estimates of ∼ 0.3 m for the Kuhio lava tube in Hawai’i, which did not agree with observations on the ground. This is the result of the complex composition and geometry of an actual lava tube where heat transfer is controlled by more than a simple fluid filling a tube, but also by convection of gasses and fluids in a partially filled passage. Though not effective for lava tubes at this time, the model provides promising results for simple cases applied to engineering and underground coal fires.     

遥感技术在长白山天池火山熔岩流表壳类型结构研究的可行性分析

溢流式熔岩喷发是火山喷发的一种常见形式,容易造成火山灾害.从火口或裂隙呈液态喷发出的溶岩,沿地表向前流动,熔岩温度逐渐降低,黏度随之增大,流动速度逐渐减慢,直至停滞、凝结.炙热熔岩流在流动过程中,会对沿途及周边植被、建筑物和人员等造成破坏和伤亡.Andrew和Rowland(2015)对aa和pahoehoe的力学原理...     

Delineation of riparian vegetation from Landsat multi‐temporal imagery using PCA

A deficiency in crucial digital data, such as vegetation cover, in remote regions is a challenging issue for water management and planning, especially for areas undergoing rapid development, such as mining in the Pilbara, Western Australia. This is particularly relevant to riparian vegetation, which provides important ecological services and, as such, requires regional protection. The objective of this research was to develop an approach to riparian vegetation mapping at a regional scale using remotely sensed data. The proposed method was based on principal component analysis applied to multi‐temporal Normalized Difference Vegetation Index datasets derived from Landsat TM 5 imagery. To delimit the spatial extent of riparian vegetation, a thresholding method was required and various thresholding algorithms were tested. The accuracy of results was estimated for various Normalized Difference Vegetation Index multi‐temporal datasets using available ground‐truth data. The combination of a 14‐dry‐date dataset and Kittler's thresholding method provided the most accurate delineation of riparian vegetation.     

Monitoring turbid freshwater plume characteristics by means of remotely sensed data

Remotely sensed data may provide easy access for monitoring the spatial separation and obtaining the hydrodynamic characteristics of turbid freshwater plumes created by river flow in the marine environment. Traditional methods are time consuming and require great effort to produce sufficient data. In this project, integrated research has been carried out on a river to demonstrate the utility of remote sensing (RS) technology for studying fundamental theoretical properties of turbulent mixtures. The Filyos River mouth, located on the Black Sea coast of Turkey, is the research area. Flow properties, such as the horizontal dispersion coefficient, have been calculated (using Landsat TM sensor images taken on two different dates). The effects of the plume on the morphology of neighbouring beaches are also examined. This study shows the utility of RS technology for generating quantitative data and better defining the hydraulic behaviour of a river with high turbidity. Copyright © 2006 John Wiley & Sons, Ltd.     

Ecohydrological responses to surface flow across borders: Two decades of changes in vegetation greenness and water use in the riparian corridor of the Colorado River delta

Hydrological and bioclimatic processes that lead to drought may stress plants and wildlife, restructure plant community type and architecture, increase monotypic stands and bare soils, facilitate the invasion of non-native plant species and accelerate soil erosion. Our study focuses on the impact of a paucity of Colorado River surface flows from the United States (U.S.) to Mexico. We measured change in riparian plant greenness and water use over the past two decades using remotely sensed measurements of vegetation index (VI), evapotranspiration (ET) and a new annualized phenology assessment metric (PAM) for ET. We measure these long-term (2000–2019) metrics and their short-term (2014–2019) response to an environmental pulse flow in 2014, as prescribed under Minute 319 of the 1944 Water Treaty between the two nations. In subsequent years, small-directed flows were provided to restoration areas under Minute 323. We use 250 m MODIS and 30 m Landsat imagery to evaluate three vegetation indices (NDVI, EVI, EVI2). We select EVI2 to parameterize an optical-based ET algorithm and test the relationship between ET from Landsat and MODIS by regression approaches. Our analyses show significant decreases in VIs and ET for both the 20-year and post-pulse 5-year periods. Over the last 20 years, EVI _Landsat declined 34% (30% by EVI_MODIS) and ET_Landsat-EVI declined 38% (27% by ET_MODIS-EVI), overall ca. 1.61 mm/day or 476 mm/year drop in ET; using PAM ET_Landsat-EVI the drop was from 1130 to 654 mm/year. Over the 5 years since the 2014 pulse flow, EVI_Landsat declined 20% (13% by EVI_MODIS) and ET_Landsat-EVI declined 23% (4% by ET_MODIS-EVI) with a 0.77 mm/day or a 209 mm/year 5-year drop in ET; using PAM ET_Landsat-EVI the drop was from 863 to 654 mm/year. Data and change maps show the pulse flow contributed enough water to slow the rate of loss, but only for the very short-term (1–2 years). These findings are critically important as they suggest further deterioration of biodiversity, wildlife habitat and key ecosystem services due to anthropogenic diversions of water in the U.S. and Mexico and from land clearing, fires and plant-related drought which affect hydrological processes.     

Investigation of meandering in Filyos River by means of satellite sensor data

The most significant morphological property of a river is the meandering process, which is dominated and governed by hydraulic, hydrologic and topographic characteristics of the river and its drainage area. It is possible to obtain reliable data on river morphology in the long term by using remotely sensed data. In this study the Filyos River, located at the Western Black Sea region of Turkey, has been selected as the study area to show the capabilities of remote sensing technology and to define the appropriate techniques for achieving the most reliable information on the river morphology by monitoring the meandering processes. The findings of the study indicate that remotely sensed data can be used successfully in defining some basic characteristics of the meandering process on rivers. Copyright © 2004 John Wiley & Sons, Ltd.     

Vegetation-groundwater dynamics at a former uranium mill site following invasion of a biocontrol agent: A time series analysis of Landsat normalized difference vegetation index data

Because groundwater recharge in dry regions is generally low, arid and semiarid environments have been considered well-suited for long-term isolation of hazardous materials (e.g., radioactive waste). In these dry regions, water lost (transpired) by plants and evaporated from the soil surface, collectively termed evapotranspiration (ET), is usually the primary discharge component in the water balance. Therefore, vegetation can potentially affect groundwater flow and contaminant transport at waste disposal sites. We studied vegetation health and ET dynamics at a Uranium Mill Tailings Radiation Control Act (UMTRCA) disposal site in Shiprock, New Mexico, where a floodplain alluvial aquifer was contaminated by mill effluent. Vegetation on the floodplain was predominantly deep-rooted, non-native tamarisk shrubs (Tamarix sp.). After the introduction of the tamarisk beetle (Diorhabda sp.) as a biocontrol agent, the health of the invasive tamarisk on the Shiprock floodplain declined. We used Landsat normalized difference vegetation index (NDVI) data to measure greenness and a remote sensing algorithm to estimate landscape-scale ET along the floodplain of the UMTRCA site in Shiprock prior to (2000–2009) and after (2010–2018) beetle establishment. Using groundwater level data collected from 2011 to 2014, we also assessed the role of ET in explaining seasonal variations in depth to water of the floodplain. Growing season scaled NDVI decreased 30% (p < .001), while ET decreased 26% from the pre- to post-beetle period and seasonal ET estimates were significantly correlated with groundwater levels from 2011 to 2014 (r² = .71; p = .009). Tamarisk greenness (a proxy for health) was significantly affected by Diorhabda but has partially recovered since 2012. Despite this, increased ET demand in the summer/fall period might reduce contaminant transport to the San Juan River during this period.     

Large‐scale regional delineation of riparian vegetation in the arid and semi‐arid Pilbara region,WA

Multiscene Landsat 5 TM imagery, Principal Component Analysis, and the Normalized Difference Vegetation Index were used to produce the first region‐scale map of riparian vegetation for the Pilbara (230,000 km²), Western Australia. Riparian vegetation is an environmentally important habitat in the arid and desert climate of the Pilbara. These habitats are supported by infrequent flow events and in some locations by groundwater discharge. Our analysis suggests that riparian vegetation covers less than 4% of the Pilbara region, whereas almost 10.5% of this area is composed of groundwater dependent vegetation (GDV). GDV is often associated with open water (river pools), providing refugia for a variety of species. GDV has an extremely high ecological value and are often important Indigenous sites. This paper demonstrates how Landsat data calibrated to Top of Atmosphere reflectance can be used to delineate riparian vegetation across 16 Landsat scenes and two Universal Transverse Mercator spatial zones. The proposed method is able to delineate riparian vegetation and GDV, without the need for Bidirectional Reflectance Distribution Function correction. Results were validated using ground truth data from local and regional scale vegetation surveys.     

Wide‐area estimates of evapotranspiration by red gum (Eucalyptus camaldulensis) and associated vegetation in the Murray–Darling River Basin,Australia

Floodplain red gum forests (Eucalyptus camaldulensis plus associated grasses, reeds and sedges) are sites of high biodiversity in otherwise arid regions of southeastern Australia. They depend on periodic floods from rivers, but dams and diversions have reduced flood frequencies and volumes, leading to deterioration of trees and associated biota. There is a need to determine their water requirements so environmental flows can be administered to maintain or restore the forests. Their water requirements include the frequency and extent of overbank flooding, which recharges the floodplain soils with water, as well as the actual amount of water consumed in evapotranspiration (ET). We estimated the flooding requirements and ET for a 38 134 ha area of red gum forest fed by the Murrumbidgee River in Yanga National Park, New South Wales. ET was estimated by three methods: sap flux sensors placed in individual trees; a remote sensing method based on the Enhanced Vegetation Index from MODIS satellite imagery and a water balance method based on differences between river flows into and out of the forest. The methods gave comparable estimates yet covered different spatial and temporal scales. We estimated flood frequency and volume requirements by comparing Normalized Difference Vegetation Index values from Landsat images with flood history from 1995 to 2014, which included both wet periods and dry periods. ET during wet years is about 50% of potential ET but is much less in dry years because of the trees' ability to control stomatal conductance. Based on our analyses plus other studies, red gum trees at this location require environmental flows of 2000 GL yr^?1 every other year, with peak flows of 20 000 ML d^?1, to produce flooding sufficient to keep them in good condition. However, only about 120–200 GL yr^?1 of river water is consumed in ET, with the remainder flowing out of the forest where it enters the Murray River system. Copyright © 2015 John Wiley & Sons, Ltd.     

Post‐Younger Dryas deglaciation of the Greenland western margin as revealed by spatial analysis of lakes

Lake shapes and their spatial distribution are important geomorphological indicators in previously glaciated areas. Their shapes are influenced by the underlying geological structure and processes of glacial sediment deposition or erosion. Since these processes act on large areas, distribution of lakes can reflect the intensity of glacial erosional/depositional processes and their spatial extent. Landsat imagery was used to extract lake outlines from a selected pilot‐study area on the widest ice‐free coastal margin of the south‐western Greenland north of Kangerlussuaq. Analysis included image classification and spatial analysis of lakes with elevation data using geographic information system (GIS) tools. A morphometric index was applied to extract kettle lakes as indicators of a specific glacial process – ice stagnation. Analysis of their spatial distribution helped in the reconstruction of glacial dynamics in formerly glaciated terrain. Our results show that spatial lake distribution combined with elevation analysis can be used to identify zones of glacial erosion and deposition. The highest concentrations of lakes within the study area occupy the elevation range between 164 and 361 m above sea level (a.s.l.). This zone can be identified as an area where intensive glacial erosion took place in the past. The widespread distribution of modeled kettle lake features within the same elevation range and across the study area suggests that the last deglaciation process was accompanied by abandonment of blocks of stagnant ice. This conclusion is supported by surface exposure ages obtained in the same study area and published elsewhere. Copyright © 2009 John Wiley & Sons, Ltd.     

Monitoring of the sediment dynamics along a sandy shoreline by means of airborne hyperspectral remote sensing and LIDAR: a case study in Belgium

Airborne hyperspectral data and airborne laserscan or LIDAR data were applied to analyse the sediment transport and the beach morphodynamics along the Belgian shoreline. Between 2000 and 2004, four airborne acquisitions were performed with both types of sensor. The hyperspectral data were classified into seven sand type classes following a supervised classification approach, in which feature selection served to reduce the number of bands in the hyperspectral data. The seven classes allowed us to analyse the spatial dynamics of specific sediment volumes. The technique made it possible to distinguish the sand used for berm replenishment works or for beach nourishments from the sand naturally found on the backshore and the foreshore. Subtracting sequential DTMs (digital terrain models) resulted in height difference maps indicating the erosion and accretion zones. The combination of both data types, hyperspectral data and LIDAR data, provides a powerful tool, suited to analyse the dynamics of sandy shorelines. The technique was demonstrated on three sites along the Belgian shoreline: Koksijde, located on the West Coast and characterized by wide accretional beaches, influenced by dry berm replenishment works and the construction of groins; Zeebrugge, on the Middle Coast, where a beach nourishment was executed one year before the acquisitions started and where the dams of the harbour of Zeebrugge are responsible for the formation of a large accretional beach, and Knokke‐Heist, located on the East Coast and characterized by narrow, locally reflective, beaches, heavily influenced by nourishment activities. The methodology applied allowed retrieval of the main sediment transport directions as well as the amount of sediment transported. It proved to be specifically suited to follow up the redistribution and the re‐sorting of the fill in beach nourishment areas. Copyright © 2007 John Wiley & Sons, Ltd.