The Lilstock outcrop in the southern Bristol Channel provides exceptional exposures of several limestone beds displaying stratabound fracture networks, providing the opportunity to create a very large, complete, and ground-truthed fracture model. Here we present the result of automated fracture extraction of high-resolution photogrammetric images (0.9 cm/pixel) of the full outcrop, obtained using an unmanned aerial vehicle, to obtain a spatially extensive, full-resolution map of the complete fracture network with nearly 350,000 ground-truthed fractures. We developed graph-based functions to resolve some common issues that arise in automatic fracture tracing such as incomplete traces, incorrect topology, artificial fragmentation, and linking of fracture segments to generate geologically significant trace interpretations. The fracture networks corresponding to different regions within the outcrop are compared using several network metrics and the results indicate both inter- and intra-network (layer to layer) structural variabilities. The dataset is a valuable benchmark in the study of large-scale natural fracture networks and its extension to stochastic network generation in geomodelling. The dataset also highlights the intrinsic spatial variation in natural fracture networks that can occur even in weakly-deformed rocks over relatively short length scales of tens of metres.

The 100 000 m 2 wave-cut pavement in the Bristol Channel near Lilstock, UK, is a world-class outcrop, perfectly exposing a very large fracture network in several thin limestone layers. We present an analysis based on manual interpretation of fracture generations in selected domains and compare it with automated fracture tracing. Our dataset of high-resolution aerial photographs of the complete outcrop was acquired by an unmanned aerial vehicle, using a survey altitude optimized to resolve all fractures. We map fractures and identify fracture generations based on abutting and overprinting criteria, and we present the fracture networks of five selected representative domains. Each domain is also mapped automatically using ridge detection based on the complex shearlet transform method. The automatic fracture detection technique provides results close to the manually traced fracture networks in shorter time but with a bias towards closely spaced Y over X nodes. The assignment of fractures into generations cannot yet be done automatically, because the fracture traces extracted by the automatic method are segmented at the nodes, unlike the manual interpretation in which fractures are traced as a path from fracture tip to fracture tip and consist of several connected segments. This segmentation makes an interpretation of relative age impossible, because the identification of correct abutting relationships requires the investigation of the complete fracture trace by following a clearly defined set of rules. Generations 1 and 2 are long fractures that traverse all domains. Generation 3 is only present in the southwestern domains. Generation 4 follows an ENE-WSW striking trend, is suborthogonal to generations 1 and 2, and abuts on them and generation 3, if present. Generations 5 is the youngest fracture set with a range of orientations, creating polygonal patterns by abutting at all other fracture generations. Our mapping results show that the northeastern domains only contain four fracture generations; thus, the five generations of the outcrop identified in the southwestern domains are either not all present in each of the five domains or vary locally in their geometry, preventing the interpreter from linking the fractures to their respective generation over several spatially separate mapping domains. Fracture intensities differ between domains where the lowest is in the NE with 7.3 m -1 and the highest is in the SW with 10 m -1, coinciding with different fracture orientations and distributions of abutting relationships. Each domain has slightly different fracture network characteristics, and greater connectivity occurs where the development of later shorter fractures is not affected by the stress shadowing of pre-existing longer fractures.