Identifying large‐scale erosion and deposition processes from airborne gamma radiometrics and digital elevation models in a weathered landscape

Patterns of erosion and deposition are difficult to identify and measure at catchment and regional scale but it may be possible to infer their distribution from remote sensing using easily measured surrogate variables. Airborne geophysical surveys provide data on gamma ray emissions from surface and near‐surface material and allow estimation of K, Th and U content. Gamma ray signatures are largely determined by lithology but also change with weathering and with erosion and deposition, and may be used as a partial surrogate for those processes. Comparison of gamma ray signatures with topographic characteristics closely related to sediment transport capacity, downstream sorting of sediments, and the extent of erosion and deposition shows strong linkages. Studies in four small catchments in southeastern Australia indicate that K consistently increases as slopes become steeper while Th and U may either increase or decrease. This suggests the presence of fresh rock rather than weathered material and implies removal of material by erosion, although some patterns may result from systematic changes in lithology across catchments. Analysis by lithology confirms the increase in K with slope in granites, metamorphosed sediments and basalt, and also shows a tendency for U to decrease with slope in the granites and basalt. Gamma ray emissions vary only slightly with catchment area (a surrogate measure of water discharge) suggesting that water erosion is limited or that discharge is not closely related to area. Gamma radiometric profiles down hillslopes, averaged across the full range of gradients, show that, in most cases, radioelements initially decrease, probably reflecting increased weathering, but then increase close to valley floors because of accumulation of fine sediments. Analysis by lithology confirms the increase in radioelements close to valley floors in granites and metamorphosed sediments but the trend is less clear in basalt. Gamma ray profiles down floodplains are variable and indicate the amount of deposition and accumulation of weathered material in valley openings and exposure of bedrock in valley constrictions. Simple erosion and deposition models, based on the conservation of mass equation, and applied to the four catchments, show that all radioelements increase as the potential for deposition increases. They reach a low point in zones of no net erosion or deposition and subsequently increase as erosion becomes more intense and weathered material is lost from slopes. Analysis by lithology largely confirms this pattern. The usefulness of airborne geophysical survey data is limited by flight line spacing with most data being flown at a 200–400 m spacing in Australia. However, general trends in erosion and deposition can still be distinguished and there is capacity for calibrating long‐term erosion and deposition models once better approaches to interpretation of gamma ray data have been developed. Copyright © 2000 John Wiley & Sons, Ltd.     

Quantifying rainfall controls on catchment‐scale landslide erosion in Taiwan

Landslide erosion is a dominant hillslope process and the main source of stream sediment in tropical, tectonically active mountain belts. In this study, we quantified landslide erosion triggered by 24 rainfall events from 2001 to 2009 in three mountainous watersheds in Taiwan and investigated relationships between landslide erosion and rainfall variables. The results show positive power‐law relations between landslide erosion and rainfall intensity and cumulative rainfall, with scaling exponents ranging from 2·94 to 5·03. Additionally, landslide erosion caused by Typhoon Morakot is of comparable magnitude to landslide erosion caused by the Chi‐Chi Earthquake (M_W = 7·6) or 22–24 years of basin‐averaged erosion. Comparison of the three watersheds indicates that deeper landslides that mobilize soil and bedrock are triggered by long‐duration rainfall, whereas shallow landslides are triggered by short‐duration rainfall. These results suggest that rainfall intensity and watershed characteristics are important controls on rainfall‐triggered landslide erosion and that severe typhoons, like high‐magnitude earthquakes, can generate high rates of landslide erosion in Taiwan. Copyright © 2012 John Wiley & Sons, Ltd.     

Evidence for elevated coastal vulnerability following large‐scale historical oyster bed harvesting

Living coastal barriers, such as coral reefs, tidal marshes, mangroves and shellfish beds are widely recognized for their potential role in mitigating flood risk. Limited data exists, however, for assessing the effectiveness of these natural defenses as forms of flood mitigation. In particular, very few mature shellfish beds exist today for modern study due to their destruction in the past few centuries. As an alternative method of study, we present here sedimentary reconstructions of storm overwash from coastal ponds internal to New York Harbor. We use these reconstructions to show that the initial degradation of oyster beds following European settlement of the area coincides with a significant increase in wave‐derived overwash deposition at all three of our field sites. Numerical simulations of two flood events of record in the harbor (Hurricane Sandy and a severe winter storm in 1992) were run without and with oyster beds of varying heights (1 m above the seafloor‐to‐intertidal). Simulations show that the removal of these oyster beds increases wave energy directly off‐shore of our field sites by between 30% and 200%. Sedimentary reconstructions and wave modeling experiments therefore both support oyster beds serving as a significant form of coastal protection prior to European disturbance. Copyright © 2016 John Wiley & Sons, Ltd.     

Removing sun glint from optical remote sensing images of shallow rivers

Sun glint is the specular reflection of light from the water surface, which often causes unusually bright pixel values that can dominate fluvial remote sensing imagery and obscure the water‐leaving radiance signal of interest for mapping bathymetry, bottom type, or water column optical characteristics. Although sun glint is ubiquitous in fluvial remote sensing imagery, river‐specific methods for removing sun glint are not yet available. We show that existing sun glint‐removal methods developed for multispectral images of marine shallow water environments over‐correct shallow portions of fluvial remote sensing imagery resulting in regions of unreliable data along channel margins. We build on existing marine glint‐removal methods to develop a river‐specific technique that removes sun glint from shallow areas of the channel without over‐correction by accounting for non‐negligible water‐leaving near‐infrared radiance. This new sun glint‐removal method can improve the accuracy of spectrally‐based depth retrieval in cases where sun glint dominates the at‐sensor radiance. For an example image of the gravel‐bed Snake River, Wyoming, USA, observed‐versus‐predicted R² values for depth retrieval improved from 0.66 to 0.76 following sun glint removal. The methodology presented here is straightforward to implement and could be incorporated into image processing workflows for multispectral images that include a near‐infrared band. Copyright © 2016 John Wiley & Sons, Ltd.     

SCALES: a large‐scale assessment model of soil erosion hazard in Basse‐Normandie (northern‐western France)

The cartography of erosion risk is mainly based on the development of models, which evaluate in a qualitative and quantitative manner the physical reproduction of the erosion processes (CORINE, EHU, INRA). These models are mainly semi‐quantitative but can be physically based and spatially distributed (the Pan‐European Soil Erosion Risk Assessment, PESERA). They are characterized by their simplicity and their applicability potential at large temporal and spatial scales. In developing our model SCALES (Spatialisation d'éChelle fine de l'ALéa Erosion des Sols/large‐scale assessment and mapping model of soil erosion hazard), we had in mind several objectives: (1) to map soil erosion at a regional scale with the guarantee of a large accuracy on the local level, (2) to envisage an applicability of the model in European oceanic areas, (3) to focus the erosion hazard estimation on the level of source areas (on‐site erosion), which are the agricultural parcels, (4) to take into account the weight of the temporality of agricultural practices (land‐use concept). Because of these objectives, the nature of variables, which characterize the erosion factors and because of its structure, SCALES differs from other models. Tested in Basse‐Normandie (Calvados 5500 km²) SCALES reveals a strong predisposition of the study area to the soil erosion which should require to be expressed in a wet year. Apart from an internal validation, we tried an intermediate one by comparing our results with those from INRA and PESERA. It appeared that these models under estimate medium erosion levels and differ in the spatial localization of areas with the highest erosion risks. SCALES underlines here the limitations in the use of pedo‐transfer functions and the interpolation of input data with a low resolution. One must not forget however that these models are mainly focused on an interregional comparative approach. Therefore the comparison of SCALES data with those of the INRA and PESERA models cannot result on a convincing validation of our model. For the moment the validation is based on the opinion of local experts, who agree with the qualitative indications delivered by our cartography. An external validation of SCALES is foreseen, which will be based on a thorough inventory of erosion signals in areas with different hazard levels. Copyright © 2010 John Wiley & Sons, Ltd.     

The influence of alongshore and cross‐shore wave energy flux on large‐ and small‐scale coastal erosion patterns

A strong low‐pressure system traveled along the Japanese main island Honshu in October 2006. High waves and storm surge attacked the Kashima Coast resulting in huge erosion over the area. Airborne laser data measured in October 2005 and November 2006 were analyzed to estimate cross‐sectional changes within the subaerial zone. The results of the alongshore distribution of the changes of cross‐sectional area indicate that the amount of erosion of the 38 km‐long northern and 15 km‐long southern parts decreased toward the south in each part and that the amount of erosion was smaller in protected areas with artificial headlands than in unprotected areas. The local alongshore variation of the erosion and accretion patterns showed wavy fluctuations of several hundreds of meters. The total amounts of the estimated eroded volume of the subaerial zone over the northern and southern parts were 620 000 m³ and 600 000 m³, respectively. The Simulating Waves Nearshore (SWAN) wave model was applied to estimate wave conditions along the coast during the storm. The computational results were verified, and then the alongshore distribution of wave energies, expressed as the alongshore and cross‐shore components of the wave energy flux, was compared with the alongshore distribution of cross‐sectional change. The results show that the distribution of energy flux explains the distribution of erosion well: The alongshore variability in the cross‐shore energy flux is responsible for the large‐scale variability in erosion, and shorter‐scale variability is due to gradients in the alongshore energy fluxes, especially for the areas without coastal works. Copyright © 2011 John Wiley & Sons, Ltd.     

Computation of bridge seismic fragility by large‐scale simulation for probabilistic resilience analysis

Seismic resilience of structures and infrastructure systems is a fast developing concept in the field of disaster management, promoting communities that are resistant and quickly recoverable in case of an extreme event. In this contest, probabilistic seismic demand and fragility analyses are two key elements of the seismic resilience assessment in the majority of the proposed methodologies. Several techniques are available to calculate fragility curves for different types of structures. In particular, to assess the seismic performance of the regional transportation infrastructure, methods for the fragility curve estimation for entire classes of bridges are required. These methods usually rely on a set of assumptions, partially because of the limited information. Other assumptions were introduced at the time when computational resources were inadequate for a purely numerical approach and closed‐form solutions were a convenient alternative. For instance, some of these popular assumptions are aimed at simplifying the model of the engineering demand. In this paper, a simulation‐based methodology is proposed, to take advantage of the computational resources widely available today and avoid such assumptions on the demand. The resulting increase in accuracy is estimated on a typical class of bridges (multi‐span simply supported). Most importantly, the quantitative impact of the assumptions is assessed in the context of a life‐cycle loss estimation analysis and resilience analysis. The results show that some assumptions preserve an acceptable level of accuracy, but others introduce a considerable error in the fragility curves and, in turn, in the expected resilience and life‐cycle losses of the structure. Copyright © 2015 John Wiley & Sons, Ltd.     

Basin‐scale stream‐flow forecasting using the information of large‐scale atmospheric circulation phenomena

It is well recognized that the time series of hydrologic variables, such as rainfall and streamflow are significantly influenced by various large‐scale atmospheric circulation patterns. The influence of El Niño‐southern oscillation (ENSO) on hydrologic variables, through hydroclimatic teleconnection, is recognized throughout the world. Indian summer monsoon rainfall (ISMR) has been proved to be significantly influenced by ENSO. Recently, it was established that the relationship between ISMR and ENSO is modulated by the influence of atmospheric circulation patterns over the Indian Ocean region. The influences of Indian Ocean dipole (IOD) mode and equatorial Indian Ocean oscillation (EQUINOO) on ISMR have been established in recent research. Thus, for the Indian subcontinent, hydrologic time series are significantly influenced by ENSO along with EQUINOO. Though the influence of these large‐scale atmospheric circulations on large‐scale rainfall patterns was investigated, their influence on basin‐scale stream‐flow is yet to be investigated. In this paper, information of ENSO from the tropical Pacific Ocean and EQUINOO from the tropical Indian Ocean is used in terms of their corresponding indices for stream‐flow forecasting of the Mahanadi River in the state of Orissa, India. To model the complex non‐linear relationship between basin‐scale stream‐flow and such large‐scale atmospheric circulation information, artificial neural network (ANN) methodology has been opted for the present study. Efficient optimization of ANN architecture is obtained by using an evolutionary optimizer based on a genetic algorithm. This study proves that use of such large‐scale atmospheric circulation information potentially improves the performance of monthly basin‐scale stream‐flow prediction which, in turn, helps in better management of water resources. Copyright © 2007 John Wiley & Sons, Ltd.     

Field‐based and spectral indicators for soil erosion mapping in semi‐arid mediterranean environments (Coastal Cordillera of central Chile)

The Coastal Cordillera of central Chile is naturally sensitive to soil erosion due to moderate to steep slopes, intense winter rains when the vegetation cover is scarce, and deeply weathered granitic rocks. In 1965, 60 per cent of its surface was moderately to very severely eroded. Today this process is still largely active, but no data are currently available to evaluate the real extent, distribution and severity of soil degradation on a regional scale. This information is vital to support efficient soil conservation plans. A multi‐scale approach was implemented to produce regional land degradation maps based on remote sensing technologies. Fieldwork has shown that the surface colour or ‘redness’ and the density of coarse fragments are pertinent erosion indicators to describe a typical sequence of soil degradation in the context of mediterranean soil developed on granitic materials and micaschists. Field radiometric experiments concluded that both factors influence the reflectance of natural surfaces and can be modelled using radiometric indices accessible from most satellites operating in the optical domain, i.e. redness index and brightness index. Finally the radiometric indices were successfully applied to SPOT images to produce land degradation maps. Only broad classes of erosion status were discriminated and the detection of the degradation processes was only possible when most of the fertile layer had already been removed. This technology provides decision‐making information required to develop regional soil conservation plans and to prioritize actions between catchment areas, especially in vast inter‐tropical regions where spatialized data are not always readily available. Copyright © 2006 John Wiley & Sons, Ltd.     

Perspectives in using a remotely sensed dryness index in distributed hydrological models at the river‐basin scale

In a previous study a spatially distributed hydrological model, based on the MIKE SHE code, was constructed and validated for the 375 000 km² Senegal River basin in West Africa. The model was constructed using spatial data on topography, soil types and vegetation characteristics together with time‐series of precipitation from 112 stations in the basin. The model was calibrated and validated based on river discharge data from nine stations in the basin for 11 years. Calibration and validation results suggested that the spatial resolution of the input data in parts of the area was not sufficient for a satisfactory evaluation of the modelling performance. The study further examined the spatial patterns in the model input and output, and it was found that particularly the spatial resolution of the precipitation input had a major impact on the model response. In an attempt to improve the model performance, this study examines a remotely sensed dryness index for its relationship to simulated soil moisture and evaporation for six days in the wet season 1990. The index is derived from observations of surface temperature and vegetation index as measured by the NOAA Advanced Very High Resolution Radiometer (AVHRR) sensor. The correlation results between the index and the simulation results are of mixed quality. A sensitivity analysis, conducted on both estimates, reveals significant uncertainties in both. The study suggests that the remotely sensed dryness index with its current use of NOAA AVHRR data does not offer information that leads to a better calibration or validation of the simulation model in a spatial sense. The method potentially may become more suitable with the use of the upcoming high‐resolution temporal Meteosat Second Generation data. Copyright © 2002 John Wiley & Sons, Ltd.     

Flow resistance in gravel‐bed channels with large‐scale roughness

A previously published mixing length (ML) model for evaluating the Darcy–Weisbach friction factor for a large‐scale roughness condition (depth to sediment height ratio ranging from 1 to 4) is brie?y reviewed and modi?ed (MML). Then the MML model and a modi?ed drag (MD) model are experimentally tested using laboratory measurements carried out for gravel‐bed channels and large‐scale roughness condition. This analysis showed that the MML gives accurate estimates of the Darcy–Weisbach coef?cient and for Froude number values greater than 0·5 the MML model coincides with the ML one. Testing of the MD model shows limited accuracy in estimating ?ow resistance. Finally, the MML and MD models are compared with the performance of a quasi‐theoretical (QT) model deduced applying the P‐theorem of the dimensional analysis and the incomplete self‐similarity condition for the depth/sediment ratio and the Froude number. Using the experimental gravel‐bed data to calibrate the QT model, a constant value of the exponent of the Froude number is determined while two relationships are proposed for estimating the scale factor and the exponent of the depth/sediment ratio. This indirect estimate procedure of the coef?cients (b₀, b₁ and b₂) of the QT model can produce a negligible overestimation or underestimation of the friction factor. Copyright © 2003 John Wiley & Sons, Ltd.     

Damage assessment through structural identification of a three‐story large‐scale precast concrete structure

This paper investigates the damage assessment of a three‐story half‐scale precast concrete building resembling a parking garage through structural identification. The structure was tested under earthquake‐type loading on the NEES large high‐performance outdoor shake table at the University of California San Diego in 2008. The tests provide a unique opportunity to capture the dynamic performance of precast concrete structures built under realistic boundary conditions. The effective modal parameters of the structure at different damage states have been identified from white‐noise and scaled earthquake test data with the assumption that the structure responded in a quasi‐linear manner. Modal identification has been performed using the deterministic‐stochastic subspace identification method based on the measured input–output data. The changes in the identified modal parameters are correlated to the observed damage. In general, the natural frequencies decrease, and the damping ratios increase as the structure is exposed to larger base excitations, indicating loss of stiffness, development/propagation of cracks, and failure in joint connections. The analysis of the modal rotations and curvatures allowed the localization of shear and flexural damages respectively and the checking of the effectiveness of repair actions. Copyright © 2013 John Wiley & Sons, Ltd.     

Cyclic performance of large‐scale corroded reinforced concrete beams

Nine large‐scale beam specimens were constructed. Of which, one was used as the control, whereas the other eight ones were divided into four sets. Each set had two specimens and was subjected to accelerated corrosion using an imposed current for the same time interval. Following the corrosion, a specimen in each set was tested using cyclic loading to examine the seismic performance, whereas the other one was demolished to examine the extent of corrosion. Cyclic loading results indicated that with an increasing corrosion level, the ultimate drift, ductility, plastic rotation capacity, and energy dissipation of the beams initially increased and later decreased. The failure mode switched from flexural failure, largely owing to buckling of the longitudinal reinforcement to flexural‐shear failure, which is mainly caused by fracturing of the transverse reinforcement. Corrosion increased shear deformation and the spread of plasticity of the plastic hinge region. The residual flexural strength, as estimated by an empirical equation based on the maximum pit depth in the longitudinal reinforcement, closely corresponds to experimental values. Furthermore, the residual shear strength estimated based on the minimum reduced cross‐sectional area of transverse reinforcement correlates better with the experimental observations than that based on the weight loss. Copyright © 2011 John Wiley & Sons, Ltd.     

Validation of (not‐historical) large‐event near‐fault ground‐motion simulations for use in civil engineering applications

Ground‐motion simulations generated from physics‐based wave propagation models are gaining increasing interest in the engineering community for their potential to inform the performance‐based design and assessment of infrastructure residing in active seismic areas. A key prerequisite before the ground‐motion simulations can be used with confidence for application in engineering domains is their comprehensive and rigorous investigation and validation. This article provides a four‐step methodology and acceptance criteria to assess the reliability of simulated ground motions of not historical events, which includes (1) the selection of a population of real records consistent with the simulated scenarios, (2) the comparison of the distribution of Intensity Measures (IMs) from the simulated records, real records, and Ground‐Motion Prediction Equations (GMPEs), (3) the comparison of the distribution of simple proxies for building response, and (4) the comparison of the distribution of Engineering Demand Parameters (EDPs) for a realistic model of a structure. Specific focus is laid on near‐field ground motions (<10km) from large earthquakes (M_w7), for which the database of real records for potential use in engineering applications is severely limited. The methodology is demonstrated through comparison of (2490) near‐field synthetic records with 5 Hz resolution generated from the Pitarka et al (2019) kinematic rupture model with a population of (38) pulse‐like near‐field real records from multiple events and, when applicable, with NGA‐W2 GMPEs. The proposed procedure provides an effective method for informing and advancing the science needed to generate realistic ground‐motion simulations, and for building confidence in their use in engineering domains.     

Analytical modelling of a large‐scale dynamic testing facility

An analytical model, which aims at reproducing the response of a large‐scale dynamic testing facility, that is a system composed of the specimen/shaking table/reaction‐mass/airbags/dampers/soil is developed. The Lagrangian of the system is derived, under the assumption of large displacements and rotations. A set of four nonlinear differential equations is obtained and solved with numerical methods. Preliminary verifications of the derived model are carried out by reproducing both well‐known results in the literature as well as those of a lumped model employed in the design of an existing dynamic testing facility. The case‐study for validating the nonlinear equations of motion is the shaking table of the EUCENTRE Laboratory. Copyright © 2011 John Wiley & Sons, Ltd.     

The Daily Erosion Project – daily estimates of water runoff,soil detachment,and erosion

Water runoff and sediment transport from agricultural uplands are substantial threats to water quality and sustained crop production. To improve soil and water resources, farmers, conservationists, and policy‐makers must understand how landforms, soil types, farming practices, and rainfall interact with water runoff and soil erosion processes. To that end, the Iowa Daily Erosion Project (IDEP) was designed and implemented in 2003 to inventory these factors across Iowa in the United States. IDEP utilized the Water Erosion Prediction Project (WEPP) soil erosion model along with radar‐derived precipitation data and government‐provided slope, soil, and management information to produce daily estimates of soil erosion and runoff at the township scale (93 km² [36 mi²]). Improved national databases and evolving remote sensing technology now permit the derivation of slope, soil, and field‐level management inputs for WEPP. These remotely sensed parameters, along with more detailed meteorological data, now drive daily WEPP hillslope soil erosion and water runoff estimates at the small watershed scale, approximately 90 km² (35 mi²), across sections of multiple Midwest states. The revisions constitute a substantial improvement as more realistic field conditions are reflected, more detailed weather data are utilized, hill slope sampling density is an order of magnitude greater, and results are aggregated based on surface hydrology enabling further watershed research and analysis. Considering these improvements and the expansion of the project beyond Iowa it was renamed the Daily Erosion Project (DEP). Statistical and comparative evaluations of soil erosion simulations indicate that the sampling density is adequate and the results are defendable. The modeling framework developed is readily adaptable to other regions given suitable inputs. © 2017 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

Regional‐scale mapping of groundwater discharge zones using thermal satellite imagery

Mapping groundwater discharge zones at broad spatial scales remains a challenge, particularly in data sparse regions. We applied a regional scale mapping approach based on thermal remote sensing to map discharge zones in a complex watershed with a broad diversity of geological materials, land cover and topographic variation situated within the Prairie Parkland of northern Alberta, Canada. We acquired winter thermal imagery from the USGS Landsat archive to demonstrate the utility of this data source for applications that can complement both scientific and management programs. We showed that the thermally determined potential discharge areas were corroborated with hydrological (spring locations) and chemical (conservative tracers of groundwater) data. This study demonstrates how thermal remote sensing can form part of a comprehensive mapping framework to investigate groundwater resources over broad spatial scales. Copyright © 2013 John Wiley & Sons, Ltd.     

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.