Mapping block-and-ash flow hazards based on Titan 2D simulations: a case study from Mt. Taranaki,NZ

Numerical models for simulation of mass flows are typically focussed upon accurately predicting the paths, travel times and inundation from a single flow or collapse event. When considering catchment-based hazards from a volcano, this is complicated by often being faced with several possible scenarios. Over the last 800 years at Mt. Taranaki/Egmont, a number of dome growth and collapse events have resulted in the genesis and emplacement of block-and-ash flows (BAFs). Each BAF was directed northwestward by a breach in the crater rim. The latest dome collapse events in the AD 1880s and AD 1755 inundated the northwestern flank and had run-out lengths 10 km from source. Future activity of this type could have a devastating effect on the Taranaki region’s communities, infrastructure and economy. Hazard planning has involved constructing volcanic hazard maps based upon the areas inundated by past volcanic flows, with little consideration of present-day topography. Here, a numerical geophysical mass flow modelling approach is used to forecast the hazards of future comparable BAF events on NW Mt. Taranaki. The Titan2D programme encompasses a “shallow water”, continuum solution-based, granular flow model. Flow mechanical properties needed for this approach include estimates of internal and basal friction as well as the physical dimensions of the initial collapse. Before this model can be applied to Taranaki BAFs, the input parameters must be calibrated by simulating a range of past collapse events. By using AD 1860 and AD 1755 scenarios, initial collapse volumes can be well constrained and internal and basal friction angles can be evaluated through an iterative approach from previous run-out lengths. A range of possible input parameters was, therefore, determined to produce a suite of potentially inundated areas under present-day terrain. A suite of 10 forecasts from a uniformly distributed range were combined to create a map of relative probabilities of inundation by future BAF events. These results were combined in a GIS package to produce hazard zones related to user-specified hazard thresholds. Using these input parameter constraints, future hazard forecasts for this scale and type of event can also take into account changing summit and topographic configurations following future eruptive or collapse events.     

Using simulated hydrologic response to revisit the 1973 Lerida Court landslide

Hydrologically driven mass wasting in the form of landslides on steep slopes is a worldwide occurrence. High-profile events in, for example, Brazil, Chile, the Philippines, Puerto Rico, and Venezuela during the last three decades all clearly illustrate, based upon significant losses of life and property, that hydrologically driven slope instability in developed (urban) areas can be a major geologic/environmental hazard. The focus of this study is the 1973 hydrologically driven Lerida Court landslide in Portola Valley, CA, USA. Physics-based hydrologic-response simulation, with the comprehensive Integrated Hydrology Model, was employed to forensically estimate the spatiotemporal pore pressure distributions for the Lerida Court site. Slope stability, driven by the simulated pore pressure dynamics, was estimated for the Lerida Court site with the infinite slope/Factor of Safety approach. The pore pressure dynamics for the Lerida Court site were reasonably captured by the hydrologic-response simulation. The estimated time of slope failure for the Lerida Court site compares well with field observations. A recommendation is made that hydrologically driven slope stability estimates including variably saturated subsurface flow be standard protocol for development sites in steep urban settings.     

Application of drone-captured thermal imagery in aiding in the recovery of meteorites within a snow-covered strewn field

The rapid recovery of meteorites mitigates the exposure of astromaterials to the terrestrial environment and subsequent contamination. Modern fireball observatories have enabled the more accurate triangulation of fireball trajectories, which has aided in the location of strewn fields, in the case of meteorite-producing events. Despite this advancement, most meteorite searches still use manual searching to locate any meteorite falls, which is often labor-intensive and has a slow coverage rate (km² day⁻¹). Recent work has begun exploring the application of drone technology to the recovery of meteorites; however, most of this work has focused on falls in arid environments. Our study examines the utilization of drones with thermal imaging technology to aid in the recovery of meteorites that have fallen on a snow-covered field. We created a simulated strewn field that included meteorite specimens as well as Earth rocks with similar properties (“meteowrongs”). Thermal imagery was utilized to determine whether the thermal contrast between meteorites and snow could aid in the identification of meteorites. We found that the thermal contrast was significant enough that meteorites were readily identifiable within thermal images; however, it was not significant enough to distinguish between the meteorites and the meteowrongs. The utilization of thermal imagery in conjunction with visible imagery has the potential to aid in the rapid recovery of meteorites in snow-covered landscapes.     

Truncated dunes as evidence of the 2004 tsunami in North Sumatra and environmental recovery post‐tsunami

The 2004 tsunami transformed the coast of Indonesia. This research investigates a sand dune area in Lampuuk, Sumatra, that was scoured by tsunami flow. We assessed geomorphology one‐year post‐event and examine the timescale of vegetation recovery. Ground Penetrating Radar (GPR) evidence shows an eroded succession of thin dipping units, overlain by aeolian layers 0 to 50 cm thick. Incipient dunes were absent, indicating limited dune rebuilding at one‐year post‐tsunami, possibly resulting from channelised airflow and the absence of vegetation. Recolonisation by vegetation was initially limited but progressed rapidly between 2005 and 2011, highlighting the temporal non‐linearity of recovery processes.     

Teaching physical geography at university with cartoons and comic strips: Motivation,construction and usage

This article presents the cartoons and comic strips that the author draws and has used in his teaching of physical geography at the undergraduate level since 2011. In the context of an image‐based culture, this article discusses the pedagogic goals that cartoons and comic strips fulfil: enhancing learning and creativity, associating pleasure with learning, pushing students to think ‘outside the box’ and relating the students' learning experience to a media framework popular with students. Cartoons and comic strips also answer particular necessities related to the teaching of physical geography. Using characters placed in hypothetical situations, they explain the process of doing geography and being a physical geographer.     

The effects of frequency-dependent quasar variability on the celestial reference frame

We examine the relationship between source position stability and astrophysical properties of radio-loud quasars making up the International Celestial Reference Frame (ICRF2). Understanding this relationship is important for improving quasar selection and analysis strategies, and therefore reference frame stability. We construct flux density time series, known as light curves, for 95 of the most frequently observed ICRF2 quasars at both the 2.3 and 8.4 GHz geodetic very long baseline interferometry (VLBI) observing bands. Because the appearance of new quasar components corresponds to an increase in quasar flux density, these light curves alert us about potential changes in source structure before they appear in VLBI images. We test how source position stability depends on three astrophysical parameters: (1) flux density variability at X band; (2) time lag between flares in S and X bands; (3) spectral index root-mean-square (rms), defined as the variability in the ratio between S and X band flux densities. We find that the time lag between S and X band light curves provides a good indicator of position stability: sources with time lags $<$ 0.06 years are significantly more stable ( $>$ 20 % improvement in weighted rms) than sources with larger time lags. A similar improvement is obtained by observing sources with low $(<$ 0.12) spectral index variability. On the other hand, there is no strong dependence of source position stability on flux density variability in a single frequency band. These findings can be understood by interpreting the time lag between S and X band light curves as a measure of the size of the source structure. Monitoring of source flux density at multiple frequencies therefore appears to provide a useful probe of quasar structure on scales important to geodesy. The observed astrometric position of the brightest quasar component (the core) is known to depend on observing frequency. We show how multi-frequency flux density monitoring may allow the dependence on frequency of the relative core positions along the jet to be elucidated. Knowledge of the position–frequency relation has important implications for current and future geodetic VLBI programs, as well as the alignment between the radio and optical celestial reference frames.     

An analysis and comparison of observed Pleistocene South Carolina (USA) shoreline elevations with predicted elevations derived from Marine Oxygen Isotope Stages

Geological maps of South Carolina, covering > 6800 km², confirm the existence of eight preserved Pleistocene shorelines above current sea level: Marietta (+ 42.6 m), Wicomico (+ 27.4 m), Penholoway (+ 21.3 m), Ladson (+ 17.4 m), Ten Mile Hill (+ 10.7 m), Pamlico (+ 6.7 m), Princess Anne (+ 5.2 m), and Silver Bluff (+ 3 m). Current geochronologic data suggest that these eight shorelines correlate with Marine Oxygen Isotope Stages (MIS) as follows: Marietta—older than MIS 77; Wicomico—MIS 55–45; Penholoway—MIS 19 or 17; Ladson—MIS 11; Ten Mile Hill—MIS 7; Pamlico—MIS 5; Princess Anne—MIS 5; and Silver Bluff—MIS 5 or 3. Except for the MIS 5e Pamlico, and possibly the MIS 11 Ladson, the South Carolina elevations are higher than predicted by isotope proxy-based reconstructions. The < 4 m of total relief from the Pamlico to the Silver Bluff shoreline in South Carolina, for which other reconstructions suggest an expected relief of ~ 80 m, illustrates the lack of match. Our results suggest that processes affecting either post-depositional changes in shoreline elevations or the creation of proxy sea-level estimates must be considered before using paleo sea-level position on continental margins.     

Observations on the crystal structures of lueshite

Laboratory powder XRD patterns of the perovskite-group mineral lueshite from the type locality (Lueshe, Kivu, DRC) and pure NaNbO₃ demonstrate that lueshite does not adopt the same space group (Pbma; #57) as the synthetic compound. The crystal structures of lueshite (2 samples) from Lueshe, Mont Saint-Hilaire (Quebec, Canada) and Sallanlatvi (Kola, Russia) have been determined by single-crystal CCD X-ray diffraction. These room temperature X-ray data for all single-crystal samples can be satisfactorily refined in the orthorhombic space group Pbnm (#62). Cell dimensions, atomic coordinates of the atoms, bond lengths and octahedron tilt angles are given for four crystals. Conventional neutron diffraction patterns for Lueshe lueshite recorded over the temperature range 11–1,000 K confirm that lueshite does not adopt space group Pbma within these temperatures. Neutron diffraction indicates no phase changes on cooling from room temperature to 11 K. None of these neutron diffraction data give satisfactorily refinements but suggest that this is the space group Pbnm. Time-of-flight neutron diffraction patterns for Lueshe lueshite recorded from room temperature to 700 °C demonstrate phase transitions above 550 °C from Cmcm through P4/mbm to \(Pm\overline{3} m\) above 650 °C. Cell dimensions and atomic coordinates of the atoms are given for the three high-temperature phases. The room temperature to 400 °C structures cannot be satisfactorily resolved, and it is suggested that the lueshite at room temperature consists of domains of pinned metastable phases with orthorhombic and/or monoclinic structures. However, the sequence of high-temperature phase transitions observed is similar to those determined for synthetic NaTaO₃, suggesting that the equilibrated room temperature structure of lueshite is orthorhombic Pbnm.