Locating some types of random errors in Digital Terrain Models

The increasing use of Geographical Information System applications has generated a strong interest in the assessment of data quality. As an example of quantitative raster data, we analysed errors in Digital Terrain Models (DTM). Errors might be classified as systematic (strongly dependent on the production methodology) and random. The present work attempts to locate some types of randomly distributed, weakly spatially correlated errors by applying a new methodology based on Principal Components Analysis. The Principal Components approach presented is very different from the typical scheme used in image processing. A prototype implementation has been conducted using MATLAB, and the overall procedure has been numerically tested using a Monte Carlo approach. A DTM of Stockholm, with integer-valued heights varying from 0 to 59 m has been used as a testbed.The model was contaminated by adding randomly located errors, distributed uniformly within 4 m and 4m. The procedure has been applied using both spike shaped (isolated errors) and pyramid-like errors. The preliminary results show that for the former, roughly half of the errors have been located with a Type I error probability of 4.6 per cent on average, checking up to 1 per cent of the dataset. The associated Type II error of the larger errors (of exactly 4m or 4 m) drops from an initial value of 1.21 per cent down to 0.63 per cent. By checking another 1 per cent of the dataset, such error drops to 0.34 per cent implying that about 71 per cent of the 4m errors have been located; Type I error was below 11.27 per cent. The results for pyramid-like errors are slightly worse, with a Type I error of 25.80 per cent on average for the first 1 per cent effort, and a Type II error drop from an initial value of 0.81 per cent down to 0.65 per cent. The procedure can be applied both for error detection during the DTM generation and by end users. It might also be used for other types of quantitative raster data.     

P-wave velocity anisotropy in crystalline rocks

Summary. Compressional wave velocities and anisotropy coefficients determined at high hydrostatic pressures are compiled from the data published for the main types of crystalline rocks. The crack-free elastic anisotropy of igneous crustal rocks is generally very low, between 1 and 3 per cent on average. The anisotropy of metamorphic rocks is higher (up to 22 per cent), but very variable. The average anisotropy coefficients in schists and amphibolites are about 10 per cent, in gneisses between 3 and 7 per cent, and in granulites less than 3 per cent. The average anisotropy of olivine ultramafites is between 7 and 12 per cent, whereas in pyroxenites and eclogites it is usually less than 4 per cent. A comparison of ranges of average velocities and average anisotropies for the individual rock groups suggests that, whereas in the crust the lateral velocity variations are mainly due to compositional changes, in the olivine of the uppermost mantle the velocity variations due to anisotropic structures could be of the same magnitude as the variations due to inhomogeneities.     

Estimating rupture scenario likelihood based on dynamic rupture simulations: the example of the segmented Middle Durance fault, southeastern France

The Middle Durance fault system, southeastern France, is a slow active fault that produced moderate-size historical seismic events and shows evidence of at least one   M _w ≳ 6.5  event in the last 29 000 yr. Based on dynamic rupture simulation, we propose earthquake scenarios that are constrained by knowledge of both the tectonic stress field and of the 3-D geometry of the Durance fault system. We simulate dynamic rupture interaction among several fault segmentations of different strikes, dips and rakes, using a 3-D boundary integral equation method. 50 combinations of reasonable stress field orientations, stress field amplitudes and hypocentre locations are tested. The probability of different rupture evolutions is then computed. Each segment ruptures mainly as a single event (44 per cent of the 50 simulations test in this paper). However, the probability that an event triggers simultaneously along three segments is high (26 per cent), leading to a potential rupture length of 45 km. Finally, 2 per cent of the simulations occur along four adjacent segments, producing the greatest total rupture length of 55 km. The simulation results show that the southernmost segment is most easily ruptured (40 per cent), because of its favourable orientation with respect to the tectonic stress and of its favourable location for interaction with the other segments. South-bound unilateral propagation is slightly preferable (41 per cent), compared to north-bound unilateral and bilateral propagation modes. Although, these rupture scenarios cannot be directly translated into probabilities of occurrence, they do provide a better insight as to which rupture scenarios are more likely, an important element to better estimate near-field strong ground motion and seismic hazard.     

Flood Fatalities in Australia, 1788-1996

Information on fatalities from flooding in Australia has been compiled from newspapers, historical accounts, and government and scientific reports. Records covering 1788 to 1996 indicate that at least 2213 persons have been killed in floods in Australia. The overall decadal death rate has decreased from 23.98 per 100 000 population in the 1800s to 0.04 per 100 000 in the 1990s. The greatest number of fatalities has occurred in New South Wales but, taking population into account, Queensland has been the most vulnerable state historically as, more recently, has been the Northern Territory. Most fatalities have occurred during February, and amongst males (80.6 per cent). The overall male:female death and death-rate ratios are both 4:1. Whereas the male:female death-rate ratio has fluctuated over time from 10:1 to 1:1, it does not approach equality. Most fatalities (38.5 per cent) have occurred through attempts to cross creeks, bridges or roads in times of flood. Another 31.5 per cent of victims were in houses-the majority awaiting rescue, or simply unaware of the flood. Of the total known deaths, work-related fatalities have accounted for 12.4 per cent, and recreational pursuits for 5.7 per cent.     

Phase-array decomposition of a seismic section

A seismic re fraction/wide-angle reflection profile is analysed for the presence of correlated events ('phases'). The correlation problem is formulated in terms of temporally, spatially and frequency-local complex covariances. For robustness, the method concentrates on phase rather than amplitude information. This allows a computationally efficient algorithm that can make allowance for signal correlation length and can model curved wavefronts. A statistical test based on residual phase misfit across the analysed subarray is used to assess the probability that a detected event represents a real correlated signal.
With our chosen analysis parameters and confidence level (over 99.9 per cent). 1222 events were detected in the data. Using simple techniques based on 1-D earth models, detected events are associated with a small number of particular wave types. In this way, we have succeeded in classifying almost 95 per cent of the detected events. Those that remain describe those components of the data that are inconsistent with our simple ray paths in the 1-D assumption and with our prescribed tolerance. These include reverberations, near-surface guided waves and reflected waves from strongly laterally inhomogeneous structures. According to our modelling, about 25 per cent of the detected events are consistent with simple P -wave reflected energy, and these are to a very large extent (over 85 per cent) distinct from all the other wave-type models we have used. A direct mapping of the detected events into the offset-depth domain reveals dear internal and external consistencies among the detections for the various wave types. Estimated earth structure is consistent with models from previous analyses based on much larger data sets.
We have thus succeeded in extracting correlated events from the data and decomposing these, approximately but meaningfully, into distinct classes (ray paths)     

Time of day and telephone interviews in Rural‐Urban Australia

Private household telephone ownership in Australian capital cities ranges between 91.1 per cent (Canberra) and 81.9 per cent (Brisbane). Consequently, telephone ownership need not be viewed as a major factor distinguishing households. As the characteristics of households without telephones are also known, it is possible to supplement telephone interviews with face to face interviews and thereby ensure that such a sampling frame represents the population. Such techniques, coupled with the ease and relative inexpense of telephone interviews, may promote the use of the medium. Using a letter to forewarn potential respondents and a three call‐back strategy (twice during working hours, once at night), a response rate of 77.4 per cent was achieved; 8.7 per cent of respondents could not be contacted and the balance refused. Reasons for refusals related almost exclusively to the subject matter of the interview — bushfire hazard. The first day time calls yielded 166 interviews out of 378 calls and the second 67 out of 197 calls. No particular time of day or day of week resulted in significantly more interviews; however, 55.4 per cent of respondents were female. Night calls yielded the highest response rates, with 86 interviews from 127 calls, and 63.9 per cent of respondents were male.     

Spatial and temporal variations in thunderstorm casualties over India

Thunderstorms are one of the most dangerous convective weather events. Despite the recent advances in prediction of convective storms worldwide, thousands of casualties occur annually. In the present study, the authors highlight specifically, thunderstorm‐related casualties reported in India from 1978–2012. Analysis of long‐term data have revealed about 16 308 casualties resulting from 1381 thunderstorm events with an average of 465 casualties occurring annually. The maximum number of casualties were concentrated in north‐eastern and central north‐eastern states. About 80 per cent of total casualties were recorded in West Bengal (23 per cent), Assam (20 per cent), Orissa (14 per cent), Bihar (13 per cent) and Jharkhand (8 per cent) states. The national casualties rate per million population per year and casualties density standardized by area has been found to be 0.50 and 5.07, respectively. Male casualties were found to be more prominent than female and children casualties, probably due to the larger proportion of males performing their work outdoors. The number of thunderstorm events and casualties was observed to be highest during pre‐monsoon season and lowest during winter. It is believed that the findings from this study will help policy makers to draw strategies to cope with the perils of thunderstorms.     

地表物质组成与土地沙漠化

土地沙漠化是在干旱、半干旱地区,由于特定的自然条件和人类不合理地利用自然资源,所引起的环境退化过程。随着这个过程的进一步发展对人类的不利影响越来越大。因此,对土地沙漠化的问题,也越来越得到人们的普遍重视。在内蒙伊克昭盟及宁夏、陕北地区,沙漠化的土地面积,分别占总土地面积的32%、25.2%和45.9%。这说明了在上述地区土地沙漠化是相当严重的。造成土地沙漠化的自然因素,除气候的干旱、多风外,地表物质的组成是其物质基础。而地表物质组成与沉积条件密切相关,因此,不同的沉积环境形成不同的沉积物,其质地性质差异很大。这些不同的沉积物,对土地沙漠化的发生及发展速度、类型和治理措施都有直接关系。因此,研究地表物质组成与土地沙漠化之间的关系是极为重要。     

The additivity of partial thermal remanent magnetization in magnetite

Summary. Experiments were done to test the additivity of partial thermal remanent magnetizations (PTRMs) for prepared samples containing magnetite particles whose sizes range from SD (single domain) to MD (multidomain). The experiments compare the sum of two PTRMs with total-TRM, all produced by the same external field of 0.47 oe. The most significant conclusion of this paper is that, to first order, the additivity of PTRMs is obeyed for the magnetites of this study regardless of particle size. However, small, higher order deviations from additivity occur such that ΣPTRM > TRM by an average of about 1 per cent. Though small, these departures from additivity are significant at the 99 per cent confidence level, and they can be understood in terms of Néel's theory for SD particles. The small departures from additivity are intrinsic to the experimental procedure in which some particles acquire remanence twice, in each of the two PTRM steps. In the limit of small inducing fields additivity should be obeyed exactly for the magnetites of this study and for samples of interest in palaeomagnetism. The deviations from additivity should have no effect on palaeointensity determinations by the Thelliers' version of the Thellier palaeointensity method. For palaeointensity determinations by Coe's version of the Thellier method the effects of deviations from additivity would be very small, less than 4 per cent on the average for a worst-case experimental configuration, and these effects can be minimized by producing PTRMs parallel to the original NRM and by using weak laboratory fields.     

A rationale for beach selection by the public on the coast of Wales, UK

Questionnaires were completed by 2306 beach users at 19 Welsh beaches with respect to beach selection parameters. The modal group of respondents was female aged 30–39. Beach choice was primarily determined by clean litter-free sand and seawater, followed by safety. Refreshment facilities and beach awards were deemed minor considerations by the public when choosing a beach to visit. Approximately 58 per cent of respondents were aware of beach award and rating schemes. Of coastal visitors interviewed for this paper, 67 per cent rated a beach as 'important' or 'very important' to their holiday, with just 2 per cent replying that they were unimportant.     

Glacier mapping in the Cordillera Blanca,Peru, tropical Andes,using Sentinel‐2 and Landsat data

Sentinel‐2 images were used for mapping debris‐free glaciers for the first time in Cordillera Blanca. Landsat‐8 and Sentinel‐2 data were compared for glacier area estimation in 2016, obtaining comparable results. It was observed that normalized difference snow index method for glacier mapping using MSI data is less sensitive to cast shadows and steep terrain compared with Landsat data. Estimated total glacier areas in 1975, 1994, and 2016 were 726 ± 20.3 km², 576.9 ± 15.1 km² and 482.8 ± 7.4 km², respectively. Glacier area in 2016 using Landsat was slightly lower (475.7 ± 16.8 km²) compared to the area estimated using MSI data. Observed glacier shrinkage between 1975 and 2016 was 33.5 per cent, which is lower compared to observed glacier area loss in the eastern cordilleras of Peru. Glacier shrinkage was higher at northern and northeastern slopes (47.9 per cent and 48.1 per cent, respectively) compared to the south‐western slopes (11.1 per cent).     

Effects of Urbanisation on Channel Morphology: The Case of Ekulu River in Enugu, South Eastern Nigeria

This study focuses on the effects of urbanisation on the morphology and hydraulic geometry of the Ekulu River in Enugu, southeastern Nigeria. Measurements of channel properties were taken along a 25 km stretch at 45 sites grouped into three sections: section A, upstream of the town (15); section B, within the town (15); and section C, downstream of the town (15). Spatial interpolation techniques were used to predict channel morphometric properties within and downstream of Enugu on the basis of upstream relationships. The independent variable in the analysis was drainage area, used as a surrogate for discharge. In section A, channel parameters increased systematically downstream with basin area. Within the town, channel parameters exhibited sudden dimensional increases. Based on sectoral relationships, measured channel dimensions showed moderate increases above the expected values, with average enlargement ratio indices of 34 per cent, 91 per cent and 65 per cent in width, channel capacity and depth, respectively, between sections A and B. Between B and C the channel exhibited average reduction ratio indices of 21 per cent and 27 per cent in capacity and depth, respectively, and an average enlargement ratio index of 17 per cent in width relative to the expected value. The increase in channel dimension within the town is explained by: (1) the confluence with a third-order tributary; (2) the poorly structured, loose, incoherent and highly erodible channel bank materials; (3) high volumes of runoff generated from the impermeable urban surfaces; (4) urban hydrological routing due to urban drainage; (5) human traffic across the river banks; and (6) the huge volume of garbage and solid wastes dumped into the channel.     

Regional Revival in the Nord Region of France

The Nord region of northern France, comprising the départements of the Nord and the Pas-de-Calais, epitomises the economic and social problems endured by many areas of the developed world which underwent rapid industrial growth and urbanisation in the nineteenth and early twentieth centuries but have experienced economic contraction and social depression since 1950. In 1968 the Nord region contained 3.8 million inhabitants (7.6 per cent of the French population) on 12,000 km² (2.3 per cent of the national area) at an average density of 308 persons/km², by comparison with the national average of 91/km2 (Roques 1969, p. 7). Eighty-one per cent of the regional population lived in urban areas, by contrast with 67 per cent nationally.     

Tourists’ reviews of weather in five Indian Ocean islands

Favourable climate forms a key resource in attracting tourists to tropical coastal destinations. Warm temperatures and clear skies are ideal conditions for tourists to spend long days on the beach and engage in outdoor activities. By contrast, rain, strong winds, high humidity and cloud cover serve as a deterrent. These factors influence tourists’ choice of destination, timing and duration of stay, and enjoyment of their vacation, and are therefore imperative to the economic sustainability of a destination. This study explores tourists’ self-reported sensitivity to and satisfaction of weather and climatic conditions during visits to five tropical Indian Ocean Islands through the analysis of 13 618 TripAdvisor reviews. An average of 12 per cent of reviews mentioned climatic factors, indicating a relatively high sensitivity to climate, with a highest proportion of 15.7 per cent of reviews for Maldives. Rain was mentioned most frequently (21 per cent of climate mentions), followed by sunshine and hot conditions. The majority of the reviews were authored in August and May, and more than 40 per cent were written in a neutral tone, with similar but lower proportions of overtly negative and positive reviews. Thematic analysis of reviews highlight the importance of preparedness amongst tourists, emphasizing the need for accurate information communicated.     

The British Tertiary igneous province: palaeomagnetism of the Arran dykes

Summary. Five hundred and sixteen sites in Arran were sampled. After cleaning by heating and/or alternating fields 87 per cent yielded directions with α₉₅ < 10°, comprising 435 dykes, 12 sills and one granite; polarities could be determined for some of the remainder. Directions are thought to be primary and on average are shallower than corresponds to the mean centred axial dipole, with many shallow reversely-magnetized dykes present.
Seventeen per cent of the dykes are normally magnetized, 5 per cent have intermediate directions. Polarity does not correlate with petrography, thickness or trend, but does depend upon location. In particular, on the NE coast up to 70 per cent are normal, and possibly belong to a different swarm.
Normal polarities were found for both sites in the Northern Granite, for Holy Island and for some of the sills. In conjunction with the known stratigraphy this shows that the rocks were formed during an R-N-R sequence of polarities, or more probably R-N-R-N. Using the radiometric age of 58 Ma for the Northern Granite and the best estimate of the polarity time-scale the duration of activity was probably in the range 0.5 to 3.4 Ma if three polarity periods were involved, and 3.4 to 4.4 Ma if four were involved.     

First-order P-wave ray synthetic seismograms in inhomogeneous weakly anisotropic media

We propose approximate equations for P -wave ray theory Green's function for smooth inhomogeneous weakly anisotropic media. Equations are based on perturbation theory, in which deviations of anisotropy from isotropy are considered to be the first-order quantities. For evaluation of the approximate Green's function, earlier derived first-order ray tracing equations and in this paper derived first-order dynamic ray tracing equations are used.
The first-order ray theory P -wave Green's function for inhomogeneous, weakly anisotropic media of arbitrary symmetry depends, at most, on 15 weak-anisotropy parameters. For anisotropic media of higher-symmetry than monoclinic, all equations involved differ only slightly from the corresponding equations for isotropic media. For vanishing anisotropy, the equations reduce to equations for computation of standard ray theory Green's function for isotropic media. These properties make the proposed approximate Green's function an easy and natural substitute of traditional Green's function for isotropic media.
Numerical tests for configuration and models used in seismic prospecting indicate negligible dependence of accuracy of the approximate Green's function on inhomogeneity of the medium. Accuracy depends more strongly on strength of anisotropy in general and on angular variation of phase velocity due to anisotropy in particular. For example, for anisotropy of about 8 per cent, considered in the examples presented, the relative errors of the geometrical spreading are usually under 1 per cent; for anisotropy of about 20 per cent, however, they may locally reach as much as 20 per cent.     

Effects of Flow Regulation in Flow Regime on the Murrumbidgee River,South Eastern Australia: an assessment using a daily estimation hydrological model

In this paper, modelled hydrological data are used to quantify the effects of regulation on the flow regime of the lower Murrumbidgee River in the period 1970–1998. Although other studies report historical changes in flood frequency and duration, this study uses modelled natural daily flow data rather than pre-regulation records or aggregated modelled monthly data. The comparison of modelled natural and regulated daily flows shows the magnitude of changes to mean and seasonal flows, flood peaks and flow duration. At gauges upstream of major irrigation off-takes, mean flows have been increased by approximately 10 per cent, flood peaks have been reduced by 21–46 per cent, and there has been a seasonal redistribution such that flows in summer and autumn have been increased at the expense of those in winter and spring. At gauges downstream of the major irrigation off-takes, mean flows have been reduced by 8–46 per cent, flood peaks have been reduced by 16–61 per cent, and flows have been decreased in all seasons.     

Spatial and temporal patterns of fires in tropical savannas of Indonesia

This paper examines spatial and temporal patterns of large to medium‐sized fires across the Indonesian Archipelago using remote sensing, particularly the MODIS fire products, over a 14 years period (2000–13 inclusive) in conjunction with climatic and land data in a GIS. Results showed fires (burned area and number of fires) were detected the most in year 2002 and least in 2010. The extent of burning was correlated with the annual Southern Oscillation Index with most burning occurring during times of sustained negative SOI values, which generally means drier conditions across the region. Most of the detected fires occurred in the middle of the dry season and in drier eastern Indonesia. Across Indonesia, approximately 1.5 per cent of available open vegetation area burned, whereas only 0.5 per cent of forest area burned. Most open vegetation burning was detected in the savanna with an area equivalent to 11.1 per cent of its extent being burned over the 14 years. On Sumba Island, where the pronounced dry season extends for at least four months, open vegetation is the most common land cover type, of which 99 per cent burned by area over the study period making it the island with the most burned areas in the Indonesian Archipelago.     

Optical and radar remote sensing data for forest cover mapping in Peninsular Malaysia

This study aims to map forest cover in Peninsular Malaysia using satellite images as deforestation is of concern in the recent decades, and is an important environmental issue for the future too. The Carnegie Landsat Analysis System‐Lite (CLASlite) program was used in this study to detect forest cover in Peninsular Malaysia using Landsat satellite data. The results of the study show that CLASlite algorithm misclassified some oil palm, rubber and urban areas as forest vegetation. A reliable forest cover map was produced by first combining Landsat and ALOS PALSAR images to identify oil palm, rubber and urban areas, and then subsequently removing them. The HH and HV polarization data of ALOS PALSAR (threshold method) could detect oil palm plantations with 85.26 per cent of overall accuracy. For urban area detection, Enhance Build up Index (EBBI) using spectral bands from Landsat provided higher overall accuracy of 94 per cent. These methods produced a forest cover reading of 5 914 421 ha with an overall classification accuracy of 94.5 per cent. The forest cover (including rubber areas) detected in this study is 0.38 per cent higher than the percentage of 2010 forest cover detected by the Forestry Department of Peninsular Malaysia. The technique described in this paper presents an alternative and viable approach for updating forest cover maps in Malaysia.     

Retrieving earthquake signature in grace gravity solutions

The GRACE satellites have been orbiting the Earth since 2002, monitoring the time variable gravity field. Some of the observed fluctuations are due to geodynamic causes, but they are often hidden in the complex signal, composed of hydrology, ocean, atmosphere, and geodynamics, the signal of geodynamic origin being usually the smallest. In addition, dealiasing residuals and noise make the separation of the signal from the different causes more difficult. We proposed a method based on the Empirical Orthogonal Function decomposition to extract the signal of physical origin, under the hypothesis that the physical signal is spatially more consistent than the noise and aliasing incomplete correction. We used synthetic geoid variations associated with earthquakes located at nearly 2000 positions at the Earth surface, based on several examples of large actual subduction events. We show that, with the present day accuracy, we can retrieve the geoid variations associated with more than 98 per cent of the earthquakes of magnitude 9 or above, around 60 per cent for magnitude 8.8, 40 per cent for magnitude 8.6 and 33 per cent for magnitude 8.3. Some events, with the right properties and location, can be detected with magnitude as low as 8. We then applied the method to the GRACE solutions, and retrieved the Hokkaido event (2003) and the Sumatra event (2004), which is in agreement with the retrieval rates mentioned here above.