Application of matrix analysis in delineating sinkhole risk areas along highway (I-70 near Frederick, Maryland)

Sinkhole collapse in the area of Maryland Interstate 70 (I-70) and nearby roadways south of Frederick, Maryland, has been posing a threat to the safety of the highway operation as well as other structures. The occurrence of sinkholes is associated with intensive land development. However, the geological conditions that have been developing over the past 200 million years in the Frederick Valley control the locations of the sinkholes. Within an area of approximately 8 km2, 138 sinkholes are recorded and their spatial distribution is irregular, but clustered. The clustering indicates the existence of an interaction between the sinkholes. The point pattern of sinkholes is considered to be a sample of a Gibbsian point process from which the hard-core Strauss Model is developed. The radius of influence is calculated for the recorded sinkholes which are most likely to occur within 30 m of an existing sinkhole. The stochastic analysis of the existing sinkholes is biased toward the areas with intensive land use. This bias is adjusted by considering (1) topography, (2) proximity to topographic depressions, (3) interpreted rock formation, (4) soil type, (5) geophysical anomalies, (6) proximity to geologic structures, and (7) thickness of overburden. Based on the properties of each factor, a scoring system is developed and the average relative risk score for individual 30-m segments of the study area is calculated. The areas designated by higher risk levels would have greater risk of a sinkhole collapse than the areas designated by lower risk levels. This risk assessment approach can be updated as more information becomes available.