Abstract

Theoretical calculations indicate that elastic stresses induced by surface topography may be large enough in some landscapes to fracture rocks, which in turn could influence slope stability, erosion rates, and bedrock hydrologic properties. These calculations typically have involved idealized topographic profiles, with few direct comparisons of predicted topographic stresses and observed fractures at specific field sites. We use a numerical model to calculate the stresses induced by measured topographic profiles and compare the calculated stress field with fractures observed in shallow boreholes. The model uses a boundary element method to calculate the stress distribution beneath an arbitrary topographic profile in the presence of ambient tectonic stress. When applied to a topographic profile across the Susquehanna Shale Hills Critical Zone Observatory in central Pennsylvania, the model predicts where shear fractures would occur based on a Mohr–Coulomb criterion, with considerable differences in profiles of stresses with depth beneath ridgetops and valley floors. We calculate the minimum cohesion required to prevent shear failure, Cmin, as a proxy for the potential for fracturing or reactivation of existing fractures. We compare depth profiles of Cmin with structural analyses of image logs from four boreholes located on the valley floor, and find that fracture abundance declines sharply with depth in the uppermost 15 m of the bedrock, consistent with the modeled profile of Cmin. In contrast, Cmin increases with depth at comparable depths below ridgetops, suggesting that ridgetop fracture abundance patterns may differ if topographic stresses are indeed important. Thus, the present results are consistent with the hypothesis that topography can influence subsurface rock fracture patterns and provide a basis for further observational tests.

Citation

Mirna Slim, J. Taylor Perron, Stephen J. Martel and Kamini Singha (2014): Topographic stress and rock fracture: a two dimensional numerical model for arbitrary topography and preliminary comparison with borehole observations. Earth Surface Processes and Landforms, 40(4):512–529. DOI: 10.1002/esp.3646

This Paper/Book acknowledges NSF CZO grant support.