Application of 3D empirical travel times to routine event location

The unmodeled effects of lateral heterogeneity are a primary cause of event mislocation in routine methods. Approaches that account for the effect of lateral heterogeneity on travel times are commonly used in specialized studies, but some methods have recently been developed that are suitable for routine location. A variety of changes to the routine procedures have been suggested by other authors. These proposed changes are shown to be compatible to with the use of 3D travel time predictions. 3D empirical travel times are discussed in a routine context. This one step approach, based on direct use of a database of the most reliable previous observations, is shown to be suitable as an additional final step in routine procedures to give a supplementary high accuracy location. The empirical travel time approach is tested on a reference dataset of 155 ground truth events and shown to reduce epicentral mislocation by 45%. The approach is found to be particularly effective at reducing large mislocations. Improvement is demonstrated for every region represented in the dataset but particularly large improvements are shown for events in Europe, USA and Japan.     

Using structure‐from‐motion to create glacier DEMs and orthoimagery from historical terrestrial and oblique aerial imagery

Increased resolution and availability of remote sensing products, and advancements in small‐scale aerial drone systems, allows observations of glacial changes at unprecedented levels of detail. Software developments, such as structure‐from‐motion (SfM), now allow users an easy and efficient method to generate three‐dimensional (3D) models and orthoimages from aerial or terrestrial datasets. While these advancements show promise for current and future glacier monitoring, many regions still suffer a lack of observations from earlier time periods. We report on the use of SfM to extract spatial information from various historic imagery sources. We focus on three geographic regions, the European Alps, high Arctic Norway and the Nepal Himalayas. We used terrestrial field photographs from 1896, high oblique aerial photographs from 1936 and aerial handheld photographs from 1978 to generate digital elevation models (DEMs) and orthophotos of the Rhone glacier, Brøggerhalvøya and the lower Khumbu glacier, respectively. Our analysis shows that applying SfM to historic imagery can generate high quality models using only ground control points. Limited camera/orientation information was largely reproduced using self‐calibrated model data. Using these data, we calculated mean ground sampling distances across each site which demonstrates the high potential resolution of resulting models. Vertical errors for our models are ±5.4 m, ±5.2 m and ±3.3 m. Differencing shows similar patterns of thinning at lower Rhone (European Alps) and Brøggerhalvøya (Norway) glaciers, which have mean thinning rates of 0.31 m a^?1 (1896–2010) to 0.86 m a^?1 (1936–2010) respectively. On these clean ice glaciers thinning is highest in the terminus region and decreasing up‐glacier. In contrast to these glaciers, uneven topography, exposed ice‐cliffs and debris cover on the Khumbu glacier create a highly variable spatial distribution of thinning. The mean thinning rate for the Khumbu study area was found to be 0.54 ± 0.9 m a^?1 (1978–2015). Copyright © 2017 John Wiley & Sons, Ltd.     

A Case Study of Soil-Gas Sampling in Silt and Clay-Rich (Low-Permeability) Soils

Soil-gas sampling and analysis is a common tool used in vapor intrusion assessments; however, sample collection becomes more difficult in fine-grained, low-permeability soils because of limitations on the flow rate that can be sustained during purging and sampling. This affects the time required to extract sufficient volume to satisfy purging and sampling requirements. The soil-gas probe tubing or pipe and sandpack around the probe screen should generally be purged prior to sampling. After purging, additional soil gas must be extracted for chemical analysis, which may include field screening, laboratory analysis, occasional duplicate samples, or analysis for more than one analytical method (e.g., volatile organic compounds and semivolatile organic compounds). At present, most regulatory guidance documents do not distinguish between soil-gas sampling methods that are appropriate for high- or low-permeability soils. This paper discusses permeability influences on soil-gas sample collection and reports data from a case study involving soil-gas sampling from silt and clay-rich soils with moderate to extremely low gas permeability to identify a sampling approach that yields reproducible samples with data quality appropriate for vapor intrusion investigations for a wide range of gas-permeability conditions.     

Parameterisation of incoming longwave radiation over glacier surfaces in the semiarid Andes of Chile

A good understanding of radiation fluxes is important for calculating energy, and hence, mass exchange at glacier surfaces. This study evaluates incoming longwave radiation measured at two nearby glacier stations in the high Andes of the Norte Chico region of Chile. These data are the first published records of atmospheric longwave radiation measurements in this region. Nine previously published optimised parameterisations for clear sky emissivity all produced results with a root mean square error (RMSE) ~20 W?m^?2 and bias within ±5 W m^?2, which is inline with findings from other regions. Six optimised parameterisations for incoming longwave in all sky conditions were trialled for application to this site, five of which performed comparably well with RMSE on daytime data <18 W?m^?2 and bias within ±6 W?m^?2 when applied to the optimisation site and RMSE <20 W?m^?2 and bias within ±10 W m^?2 when applied to the validation site. The parameterisation proposed by Mölg et al. (J Glaciol 55:292-302, 2009) was selected for use in this region. Incorporating the proposed elevation modification into the equation reduced the bias in the modelled incoming longwave radiation for the validation site. It was found that applying the parameterisation optimised in the original work at Kilimanjaro produced good results at both the primary and validation site in this study, suggesting that this formulation may be robust for different high mountain regions.     

A comparison of the vegetation response to rainfall in the Sahel and East Africa,using normalized difference vegetation index from NOAA AVHRR

This article presents the results of a study of the relationship between rainfall and Normalized Difference Vegetation Index (NDVI) in East Africa and the Sahel. Monthly data for the years 1982 to 1985 have been analyzed. We have evaluated NDVI-rainfall relationships by vegetation type, using the major formations described by White (1983). In the article, a comparison of the differential response of vegetation growth to rainfall in the two study regions is emphasized. The most important conclusions of our research are as follows:

1. The spatial patterns of annually-integrated NDVI closely reflect mean annual rainfall.
2. There is a good relationship between rainfall variations and NDVI on seasonal and interannual time scales for areas where mean annual rainfall ranges from approximately 200 to 1200 mm.
3. In most cases, NDVI is best correlated with the rainfall total for the concurrent plus two antecedent months; the correlation is better in the Sahel than in East Africa.
4. The ratios of NDVI to rainfall are considerably higher in East Africa than in the Sahel.
5. Mean annually-integrated NDVI is linearly related to mean annual rainfall in the Sahel. In East Africa the relationship is approximately log-linear; above some threshold value of rainfall, NDVI values level off and vary minimally with rainfall.

Two possible explanations of this last conclusion are suggested: above this threshold, rainfall is no longer the limiting factor in vegetation growth and/or NDVI is not a good indicator of vegetation growth. The latter is a likely possibility since NDVI directly reflects photosynthetic activity and becomes a poor indicator of biomass (i.e., growth) as high canopy densities are reached. The NDVI-rainfall relationship for East Africa is markedly similar to the relationship between NDVI and Leaf Area Index demonstrated by Sellers (1985) and Asrar et al. (1984).     

Petrogenesis of the Merensky Reef in the Rustenburg section of the Bushveld Complex

A petrogenetic model for the Merensky Reef in the Rustenburg section of the Bushveld Complex has been developed based on detailed field and petrographic observations and electron microprobe data. The model maintains that the reef formed by reaction of hydrous melt and a partially molten cumulate assemblage. The model is devised to account for several key observations: (1) Although the dominant rock type in the Rusterburg sections is pegmatoidal feldspathic pyroxenite, there is a continous range of reef lithology from pyroxenite to pegmatoidal harzburgite and dunite, and small amounts of olivine are present in nearly all pegmatoids. (2) The pegmatoid is usually bounded above and below by chromitite seams and the basal chromitite separated from underlying norite by a centimeter-thick layer of anorthosite. The thicknesses of the two layers exhibit a well-defined, positive correlation. (3) Inclusion of pyroxenite identical to the hanging wall and of leuconorite identical to the footwall are present in the pegmatoid. The leuconorite inclusions are surrounded by thin anorthosite and chromitite layers in the same sequence as that at the base of the reef. (4) Chromite in seams adjacent to plagioclase-rich rocks is characterized by higher Mg/Mg+Fe and Al/R3 and lower Cr/R3 than that in seams adjacent to pyroxene-rich rocks. Similar variations in mineral compositions are observed across individual chromitite seams where the underlying and overlying rock types differ. The chromite compositional variations cannot be rationalized in terms of either fractional crystallization or reequilibration with surrounding silicates. It is proposed that the present reef was originally a melt-rich horizon in norite immediately overlain by relatively crystallized pyroxenite. Magmatic vapor generated by crystallization of intercumulus melt migrated upward through fractures in the cumulate pile below the protoreef. The melt-rich protoreef became hydrated because fractures were unable to propagate through it and because the melt itself was water-undersaturated. Hydration of the intercumulus melt was accompanied by melting, and the hydration/melting front migrated downward into the footwall and upward into the hanging wall. In the footwall melting resulted first in the dissolution of orthopyroxene and then of plagioclase. With continued hydration chromite was stabilized as melt alumina content increased. The regular variations in chromite compositions reflect the original gradients in melt composition at the hydration front. The stratigraphic sequence downward through the base of the reef or pegmatoid (melt)-chromitite-anorthosite-norite represents the sequence of stable mineral assemblages across the hydration/melting front. The sequence is shown to be consistent with knowledge gained from experiments on melting of hydrous mafic systems at crustal pressures. With cooling the hydrated mixture from partial melting of norite footwall and more mafic hanging wall crystallized in the sequence chromite-olivine-pyroxene-plagioclase, with peritectic loss of some olivine. Calculations of mass balance indicate that a significant proportion of the melt was lost from the melt-rich horizon. Variations in the development of the pegmatoid and associated lithologies and amount of modal olivine in the pegmatoids along the strike of the Merensky Reef resulted because the processes of hydration, melting and melt loss operated to varying extents.     

Stratiform manganese mineralisation near Inverness,Scotland: A Devonian sublacustrine hot-spring deposit?

A stratiform manganiferous horizon, up to 4 m thick but localised in aerial extent, occurs 10 km south-east of Inverness in north-east Scotland. It overlies Precambrian Grampian Group schists and is itself overlain by Middle Old Red Sandstone (Devonian) conglomerate, which elsewhere rests unconformably on the schists. The deposit was mined briefly in the early 1920s and is composed principally of fine-grained braunite with minor bustamite. Lower sections of the deposit contain cavities, infilled by late baryte, calcite and rhodochrosite, which resemble stromatactis. This horizon is underlain by a hematite-quartz layer, thought to represent a laterite formed by pre-MORS weathering of the underlying schists. The mineralogy and geochemistry of the manganese ore (enriched in Ag, As, Ba, Bi, Pb, Sb, Sr and Zn) are characteristic of a hydrothermal origin. Such an origin, when considered in conjunction with the outcrop, textures, mineralogy and stratigraphy of the deposit together with the palaeogeography of the region, can only be explained by deposition from hot-springs in a sublacustrine environment. Nearby baryte mineralisation may also have been deposited from the same geothermal system, which had the Caledonian 407 Ma Moy granitoid as heat source. The managanese ore is enriched in gold pathfinder elements and minor detrital gold has been found in the area. Gold mineralisation may therefore have also been deposited in the upper sections of this Devonian geothermal system.     

Modelling the hydrodynamics of offshore sandbanks

This paper describes the details of a quasi-three-dimensional model (3DBANK), which has been developed to investigate medium and long-term morphological evolution and development of offshore sandbanks. The model is based on a three-dimensional tidal module using the Galerkin-eigenfunction method, but also includes four sub-modules to compute: the instantaneous bedform characteristics from which the temporal and spatial variations of the shear stresses at the sea bed can be derived; the suspended sediment concentration through the water column; the bed-load and suspended sediment fluxes at a point-in-plan; and the resulting morphological changes, respectively. The model also includes the effects of the wind and waves at the sea surface, as well as the wave–current interaction (WCI), and operates with full hydrodynamic and morphodynamic interaction. The components of the model were tested against laboratory and field data, and the complete model was then applied to Middlekerke Bank off the Flemish coast where extensive field measurements were taken during the European Community (EC) funded Circulation and Sediment Transport Around Banks (CSTAB) Project using various advanced instrumentation including STABLE and HF OSCR. Comparisons of the model results with the field measurements and observations show that the model is capable of reproducing the current and wave-induced bedforms, bed roughness, tidal currents and tidal residuals around the sandbank satisfactorily, and can be used to study the long-term sandbank evolution under various offshore conditions. This paper, however, focuses on the hydrodynamic aspects of the model, while the details of the morphological components will be given in a companion paper.     

A Semi-Quantitative, Regional Precipitation Data Set for Studying African Climates of the Nineteenth Century, Part I. Overview of the Data Set

This article describes a historical archive of proxy and actual precipitation data that extends the African climate record back to the early nineteenth century. The `proxy' archive includes verbal, documentary references which contain information related to rainfall conditions, such as references to famine, drought, agriculture or the nature of the rainy season. The precipitation archive includes all observations made in Africa during the nineteenth century. It consists of records for 60 stations in Algeria, 87 stations in South Africa and 304 stations scattered over the rest of Africa. Information is particularly plentiful for the 1880s and 1890s. The two parts have been be combined into a semi-quantitative regional data set indicating annual rainfall conditions in terms of anomaly classes (e.g., normal, dry, wet). This data set extends from the early nineteenth century to 1900 and distinguishes seven anomaly classes, using numbers ranging from –3 to +3 to represent very wet, wet, good rains, normal, dry, drought, and severe drought. The regionalization is based on 90 geographical regions shown via studies of the modern precipitation record to be climatically homogeneous with respect to the interannual variability of rainfall. The regional aggregation allows the voluminous fragmentary information available in historical sources to be used systematically to produce multi-year time series that can be directly integrated into the modern record for each region. The resultant time series can also be subjected to statistical analysis, in order to investigate nineteenth century climate over Africa. Spatial detail is added to the data set by utilizing a unique methodology based on climatic teleconnections established from studies of rainfall variability over Africa. The historical information and station records have been combined into a file containing a regional anomaly value for up to 90 geographic regions and the years 1801–1900. Gaps necessarily remain in the matrix, but as early as the 1820s over 40 regions are represented. By the 1880s generally around 70 regions or more are represented.