Evaluating bedrock outcrop mapping algorithms across diverse landscapes

<p>            Mapping bedrock outcrops is useful across disciplines, but is challenging in environments where ground surface visibility is obscured. The presence of soil or bedrock affects sediment production and transport, local ecology, and runoff generation. The distribution of bedrock outcrops in an area reflects the interplay between regolith production and sediment removal. Outcrop classification methods from Terrestrial-lidar produce millimeter or centimeter resolution DEMs that are highly successful because lidar penetrates through vegetation to the ground surface. However, data availability at such high resolution is limited, and the associated computational complexity required for identifying outcrop, or other surface features, is often impractical for landscape-scale analysis. Aerial lidar datasets at ~1-m resolution (e.g., moderate resolution) are more widely available and less computationally expensive than higher resolution datasets. With increasing accessibility of moderate resolution surface data, there is a need to develop outcrop classification methods and understand the efficacy of these methods across diverse environments. Our objectives are to present a simplified technique that builds on existing methods, and to examine the success of current outcrop identification methods in a variety of landscapes.</p><p>            At moderate resolution, the two most cited metrics to differentiate bedrock from soil-mantled surfaces are based on gradient (e.g., DiBiase et al., 2012) or on surface roughness (e.g., Milodowski et al., 2015). We developed a method that simplifies and combines both metrics, and that improves overall accuracy. We applied all three methods to six landscapes in the USA. For each site, we delineated ground truth from high-resolution orthoimagery for 7-10 test patches with visible ground surface, that evenly spanned 0-100% exposed outcrop. Overall accuracy, true positive rate, and false positive rate for each patch were calculated by comparing the ground truth grids to each lidar-derived outcrop grids on a cell-by-cell basis. Metric success was evaluated for each landscape by assessing the mean and distribution of performance measures across patches. Our combined metric had the highest overall accuracy in an arid, horst and graben landscape (Canyonlands National Park, Utah). It also performed well in a vegetated, high sediment load, active volcano (Mount Rainier, Washington), a canyon carved by channel incision (Boulder Canyon, Colorado), and a chaparral mixed bedrock canyon environment (Mission Trails, San Diego, California). All three methods systematically failed for portions of the landscape in glacially carved canyons (Southern Wind River Range, Wyoming) and on terraced sea cliffs (Santa Cruz County, California). These environments have significant outcrop that is both smooth and low gradient, and therefore cannot be identified using a slope or roughness-based algorithm.</p><p>            Our work highlights the importance of tailoring DEM-based bedrock mapping algorithms to its geomorphic context, and of the need for ground truth. Such data provides the basis for developing more robust methods for error evaluation. In addition, new methods are needed to identify bedrock outcrop from surface DEMs in smooth and low gradient, yet rocky landscapes.</p>

Slope Stability of a Scree Slope Based on Integrated Characterisation and Monitoring

Three years of geotechnical seasonal field monitoring including soil temperature, suction and volumetric water content plus geophysical measurements, lead to a preliminary ground model and assessment of slope stability for a steep scree slope in the Meretschibach catchment, near Agarn village in the Swiss Alps. Building on data reported in a previous paper, which focused on preliminary ground characterisation and seasonal field monitoring, this current research aims to understand whether a surficial failure in the scree slope, triggered by rainfall and depending on bedrock conditions, would represent a relevant natural hazard for Agarn village. A final year of field data is included as well as site-specific sensor calibration, a Ground Penetrating Radar (GPR) profile, and laboratory triaxial testing to provide strength parameters. A bedrock map is presented, based on GPR, with a realistic ground model of the entire scree slope. Furthermore, a preliminary numerical analysis, performed using SEEP-SLOPE/W, shows the influence of a bedrock outcrop observed in the field, for a specific soil thickness, strength parameters and rain intensity. The stability of a gravelly slope decreases with groundwater flow over a step in the bedrock, and the location of the failure will tend to move uphill of a bedrock outcrop at a shallow depth as groundwater flow increases.

Estimating unconsolidated sediment cover thickness by using the horizontal distance to a bedrock outcrop as secondary information

Unconsolidated sediment cover thickness (D) above bedrock was estimated by using a publicly available well database from Norway, GRANADA. General challenges associated with such databases typically involve clustering and bias. However, if information about the horizontal distance to the nearest bedrock outcrop (L) is included, does the spatial estimation of D improve? This idea was tested by comparing two cross-validation results: ordinary kriging (OK) where L was disregarded; and co-kriging (CK) where cross-covariance between D and L was included. The analysis showed only minor differences between OK and CK with respect to differences between estimation and true values. However, the CK results gave in general less estimation variance compared to the OK results. All observations were declustered and transformed to standard normal probability density functions before estimation and back-transformed for the cross-validation analysis. The semivariogram analysis gave correlation lengths for D and L of approx. 10 and 6 km. These correlations reduce the estimation variance in the cross-validation analysis because more than 50 % of the data material had two or more observations within a radius of 5 km. The small-scale variance of D, however, was about 50 % of the total variance, which gave an accuracy of less than 60 % for most of the cross-validation cases. Despite the noisy character of the observations, the analysis demonstrated that L can be used as secondary information to reduce the estimation variance of D.

Ordinary People and East–West Symbolism

This chapter explores the possibility that east-west symbolism is an enduring aspect of Maya culture that was initially developed in farmers' community centers and homes and later appropriated by society's nobility and encoded in hieroglyphic texts and images and enshrined in the monuments of major civic-centers. This chapter explores these ideas through the analysis of the central religious complex at the ancient Maya farming community of Chan in Belize. Chan’s central religious complex follows a construction sequence reminiscent of other E Groups, particularly the Cenote-style E Group, across the Maya area. It begins initially as a linear bedrock outcrop upon which architecture construction begins in the Late Preclassic (350 BCE-CE 0). The distance between the east and west structures of Chan’s central religious complex was always maintained throughout its architectural reconstruction history. If one of the functions of the west structure was viewing sunrise over the tripartite east structure, then maintaining a relative distance between the two structures could have been integral in maintaining lines of sight. Chan’s central religious complex was also a location for the burial of venerated ancestors, and those individuals from the community selected for veneration included men, women, and children.

Exploring the sensitivity on a soil area-slope-grading relationship to changes in process parameters using a pedogenesis model

This paper generalises the physical dependence of the relationship between contributing area, local slope, and the surface soil grading using a pedogenesis model and allows an exploration of soilscape self-organisation. A parametric study was carried out using different parent materials, erosion, and weathering mechanisms. These simulations confirmed the generality of the area-slope-d50 relationship. The relationship is also true for other statistics of soil grading (e.g. d10,d90) and robust for different depths within the profile. For small area-slope regimes (i.e. hillslopes with small areas and/or slopes) only the smallest particles can be mobilised by erosion and the area-slope-d50 relationship appears to reflect the erosion model and its Shield's Stress threshold. For higher area-slope regimes, total mobilization of the entire soil grading occurs and self-organisation reflects the relative entrainment of different size fractions. Occasionally the interaction between the in-profile weathering and surface erosion draws the bedrock to the surface and forms a bedrock outcrop. The study also shows the influence on different depth-dependent in-profile weathering functions in the formation of the equilibrium soil profile and the grading characteristics of the soil within the profile. We outline the potential of this new model and its ability to numerically explore soil and landscape properties.

Exploring the sensitivity on a soil area-slope-grading relationship to changes in process parameters using a pedogenesis model

This paper generalises the physical dependence of the relationship between contributing area, local slope, and the surface soil grading first described by Cohen et al, [2009, 2010] using their mARM1D and mARM3D pedogenesis models. A more general computational model, SSSPAM5D, extending the conceptualisation of mARM3D has been developed to further our exploration of soilscape self-organisation. A parametric study was carried out using different parent materials, erosion, and weathering mechanisms. These simulations confirmed the generality of the area-slope-d50 relationship. The relationship is also true for other statistics of soil grading (e.g. d10, d90) and robust for different depths within the profile. For small area-slope regimes (i.e. hillslopes with small areas and/or slopes) only the smallest particles can be mobilised by erosion and the area-slope-d50 relationship appears to reflect the erosion model and its Shields Stress threshold. For higher area-slope regimes, total mobilization of the entire soil grading occurs and self-organisation reflects the relative entrainment of different size fractions. Occasionally the interaction between the in-profile weathering and surface erosion draws the bedrock to the surface and forms a bedrock outcrop. The study also shows the influence on different depth dependent in-profile weathering functions in the formation of the equilibrium soil profile and the grading characteristics of the soil within the profile.

Late Wisconsin ice-flow history in the Buffalo Head Hills kimberlite field, north-central Alberta

Ice flow of the last glaciation in the Buffalo Head Hills kimberlite field of northern Alberta is reconstructed from landform interpretations and clast orientations for the purpose of aiding kimberlite exploration in the region. The paucity of bedrock outcrop and the absence of preserved striae and other erosional ice-flow indicators on the soft Cretaceous marine sediments inhibit detailed interpretations on glacial flow chronology. Poorly developed bedrock drumlins on the Buffalo Head Hills and erosional ice-flow indicators preserved on the kimberlite outcrops indicate southwestward ice flow during the maximum extent of ice during the last glaciation. During the deglaciation of northern Alberta, later phases of ice flow were controlled by lobes of surging ice, which surged into proglacial lakes. West of the Buffalo Head Hills, the maximum phase of southwest flow was followed by southeastward ice movement of the Peace River ice lobe. Similarly, east of the Buffalo Head Hills, the maximum phase of ice flow was superceded by a south-southwest ice advance of the Wasbasca ice lobe.