“…The rasters are then converted into hill shade layers, with N000 • E, N045 • E, N090 • E, and N135 • E orientation, to avoid misinterpretation of the lineament and fracture network, as explained in [12,13]. The following properties of the fracture network are extracted from GIS views (Table 3, Figures 4 and 5): length, orientation, linear density (P10), areal density (P20), areal intensity (P21), connectivity (C L ) [85], spacing (C V ) [86] and node topology [87,88]. The fracture clusters are extracted from the LiDAR and GIS analyses, with the input parameters for a stochastic distribution modelling, to be implemented in the DFN models (Tables 4 and 5).…”