The use of remote sensing techniques is now common practice in different working environments, including engineering geology. Moreover, in recent years the development of structure from motion (SfM) methods, together with rapid technological improvement, has allowed the widespread use of cost-effective remotely piloted aircraft systems (RPAS) for acquiring detailed and accurate geometrical information even in evolving environments, such as mining contexts. Indeed, the acquisition of remotely sensed data from hazardous areas provides accurate 3-D models and high-resolution orthophotos minimizing the risk for operators. The quality and quantity of the data obtainable from RPAS surveys can then be used for inspection of mining areas, audit of mining design, rock mass characterizations, stability analysis investigations and monitoring activities. Despite the widespread use of RPAS, its potential and limitations still have to be fully understood.
In this paper a case study is shown where a RPAS was used for the engineering geological investigation of a closed marble mine area in Italy; direct ground-based techniques could not be applied for safety reasons. In view of the re-activation of mining operations, high-resolution images taken from different positions and heights were acquired and processed using SfM techniques to obtain an accurate and detailed 3-D model of the area. The geometrical and radiometrical information was subsequently used for a deterministic rock mass characterization, which led to the identification of two large marble blocks that pose a potential significant hazard issue for the future workforce. A preliminary stability analysis, with a focus on investigating the contribution of potential rock bridges, was then performed in order to demonstrate the potential use of RPAS information in engineering geological contexts for geohazard identification, awareness and reduction.
In open-pit or quarry areas, personnel and equipment involved in mining operations can be exposed to different types of slope instability processes. Rock collapses can be due to a series of predisposing and triggering factors, mostly depending on relationships between localized geological conditions and mining activities. According to Zajc et al. (2014), for example, hazardous situations may occur when unfavourable sedimentological characteristics and geological discontinuities (e.g. joints, faults) of rock masses are made even more critical by extraction of the resource or ore material. In addition, Zheng et al. (2015) highlight the crucial role played by morphological features, such as sharp cuts and steep slopes, for potential triggering of rockfalls in mining areas. As widely demonstrated in the literature, the understanding of geometric relationships between geological discontinuities and slope morphology is essential to evaluate the potential occurrence of rock failures, since orientation of joint sets may influence both size and failure mechanisms of rock blocks prone to collapse (e.g. Stead and Wolter, 2015).
Generally, discontinuity characterization is carried out in the field by traditional engineering geological surveys (Priest, 1993). Measurements may be subjected to different source of errors which can result in either under- or overestimation of the discontinuity geometrical properties (Tuckey and Stead, 2016). In order to avoid this deficiency Sturzenegger and Stead (2009) suggest coupling traditional field measurements with remote sensing techniques. Indeed, techniques such as terrestrial laser scanning (TLS) and digital terrestrial photogrammetry (DTP) for rock mass characterization are increasingly being used, especially in open-pit mines where rock slopes subjected to excavation are analysed (e.g. Kovanič and Blišťan, 2014; Salvini et al., 2015; Tuckey and Stead, 2016). TLS and DTP allow accurate representation of rock outcrops by means of 3-D point clouds or interpolated models. A limitation of ground-based remote sensing is related to the survey of complex topography from sub-optimal camera or scanner positions, resulting in occlusion zones (Passalacqua et al., 2015). A solution to this problem is provided through use of remotely piloted aircraft systems (RPAS), which can be used as a platform to acquire light detection and ranging (lidar) or photogrammetric data. Indeed, RPAS photogrammetry for engineering geological investigations has become widespread mainly because it is a cost-efficient, highly flexible and safe technique (Remondino et al., 2011; Siebert and Teizer, 2014; Tannant, 2015). There are several photogrammetric studies using RPAS for the geomorphic feature characterization or mapping of the surface extent in open-pit mines (Lamb, 2000; Chen et al., 2015; Shahbazi et al., 2015; Tong et al., 2015; Esposito et al., 2017), and few studies are also associated with the use of RPAS for discontinuity characterization of rock slopes affected by mining activity. Salvini et al. (2017), for example, used RPAS to map discontinuities in a marble quarry and to subsequently build 3-D discrete fracture network models. McLeod et al. (2013) explored the feasibility of using RPAS-acquired video images to derive 3-D point clouds and to measure fracture orientations.
Digital images obtained from RPAS are commonly processed with the structure from motion (SfM) technique (Spetsakis and Aloimonos, 1991; Fonstad et al. 2013; Colomina and Molina, 2014; Westoby et al., 2012). SfM is based on sophisticated algorithms of image matching that use pseudo-random redundant images acquired from multiple viewpoints to reconstruct the 3-D geometry of an object or surface. In order to analyse rock outcrops, the use of RPAS multicopter results particularly suitable because it allows different geometric configurations for image acquisition (i.e. zenithal, frontal, oblique). Multiple images obtained from different angles help the image alignment procedure and limit non-linear deformations. Moreover, the relatively short distance to which multicopters can operate from rock faces allows acquisition of high-resolution images that can be used to produce high-quality topographic products and subsequent investigation.
In RPAS–SfM applications particular care is needed when georeferencing the
3-D model. As stated by Passalacqua et al. (2015), cameras fixed to RPAS
typically do not have onboard navigation systems with sufficient accuracy
for geodetic positioning. In fact, the global navigation satellite system
(GNSS) and inertial measurement unit (IMU) devices typically mounted on RPAS
are used for navigation and flight stabilization purposes and allow only a
rough estimation of airborne camera exterior orientation (Gonçalves and
Henriques, 2015). In order to obtain accurate and georeferenced 3-D models,
the use of ground control points (GCPs) surveyed with geodetic GNSS
receivers and a total station (TS) is generally employed (Francioni et al.,
2015). Nevertheless, the final accuracy is dependent not only on the
GCP-related accuracy, density and distribution within the surveyed area but
also on image quality and percentage of overlapping between single frames.
Therefore, careful planning of an RPAS photogrammetric survey plays a
crucial role in providing accurate results necessary for subsequent
analysis, such as determination of discontinuity measurements in terms of
attitude (dip direction and dip), spacing and persistence. The latter, for
example, is among the controlling factors that can have a great influence on
the stability condition of a block or slope and it cannot be fully
determined. The presence of intact rock bridges, which represent intervals of
intact rock between adjacent discontinuities (ISRM, 1978), can significantly
increase the stability of a rock slope, since the cohesion of the intact
rock is generally of at least 2 orders of magnitude greater than the shear
strength of a discontinuity (Park, 2005). In general, joint persistence (
The limit of the application of this method is that the discontinuity area
is practically impossible to measure deterministically in the field and, for
this reason, persistence is commonly measured as trace length on rock
outcrops. Jennings (1970) proposed the following Eq. (2) for persistence
calculation starting from trace length values on rock exposure:
Mathematically, it is possible to consider the presence of rock bridges in
terms of effective cohesion along the shear surface (Eberhardt et al., 2004)
by using the following Eq. (3):
Importantly, as recently reported by Tuckey and Stead (2016), in spite of the importance of intact rock bridges in slope stability, there are still no standard accepted methods for estimating the extent of rock bridges and incorporating rock bridges into slope stability analysis.
In this study, two RPAS-based photogrammetric surveys were carried out within an open-pit mine of the Apuan Alps marble district, Italy. These surveys aimed to obtain detailed topographic information of the area. The 3-D data were then used to perform a preliminary rockfall hazard assessment, requested in view of a potential restart of the mining operations interrupted some years ago. Indeed, the safety of the workforce represents a critical aspect for the exploitation of the marble quarries of the Apuan Alps. The geostructural conditions of the marble predispose the rock masses to different types of failures with different magnitudes. Slope stability analyses are therefore essential to improve safety conditions for personnel employed in the mines. However, a complete analysis of all the slopes characterizing an open-pit mine is often problematic, given their spatial extension and limitations of numerical models. For this reason, both geological and geomorphological information of the whole mining area are essential to detect and evaluate the most hazardous situations. RPAS-derived data were therefore integrated with those acquired in the field from a traditional engineering geological survey. The combined use of these data allowed preliminary 3-D analysis and evaluation of the stability conditions of a large rock block that was identified as posing a risk to the mining area.
The study area is located in the Apuan Alps marble district, in the province
of Lucca (Tuscany, Italy; Fig. 1). The open pit or quarry, named
Piastrone, is characterized by a
Plan view of the Piastrone open pit with indication of the two principal slope directions (dotted lines). Inset map shows the location of the study area.
Geological map of the Mt. Altissimo area. The green rectangle indicates the location of the quarry (modified from Giglia and Paiotti, 1963).
From a geological point of view (Fig. 2) the Piastrone open pit is located in the Apuan Alps metamorphic complex, precisely in the Mt. Altissimo Syncline (AS), belonging to the Apuane Unit (Meccheri et al., 2007). According to classical interpretation (Carmignani and Kligfield, 1990) AS resulted from a compressive tectonic phase which originated during the continental collision between the Sardinia–Corsica block and the Adria plate (Upper Oligocene–Lower Miocene). Successively, during the Early Miocene, a new ductile to brittle–ductile deformation caused by a post-compression tectonic uplift overprinted the earlier structures and generated a widespread network of joints and faults. In the Mt. Altissimo area, the main set of brittle deformation strikes SW–NE to W–E with sub-vertical dip, generally with negligible displacement except for a few cases where offsets of some 10 m have been observed (Meccheri et al., 2007).
AS involves the oldest lithologies of the Apuane Unit sequence, including
pre-Alpine basement rocks, Grezzoni dolostones, megalodont-bearing marbles
with metabreccia and chloritoid-rich phyllites, local lenses of dolomitic
marbles and Marbles
In order to assess and localize the slope stability hazard in the rocky
mining area, two RPAS surveys were carried out with direction of photo
acquisition in zenithal modality (perpendicular to the open-pit floor) and
in frontal modality (perpendicular to the rock faces). The surveys were
performed in December 2015 using the Aibotix™ Aibot X6 V1 multicopter,
which has six electric rotors, and equipped with a Nikon™ CoolpixA
digital camera (Table 1) and a GNSS–IMU system that allows recording of 3-D
coordinates (
Utilized RPAS and camera specifications.
Flight paths of the RPAS zenithal surveys
Example of a GCP measured during the RTK GNSS field survey
The zenithal survey was preliminarily designed in the laboratory with the Aibotix™ AiProFlight planning software and manually performed through single quasi-parallel flight lines. A total of 151 aerial images were acquired with a nominal overlap and sidelap of 80 and 60 %, respectively. Two flights were required to cover all sectors of the mining area (Fig. 3).
An average estimated distance between pixel centres measured on the ground
(i.e. ground sample distance – GSD) of 2.4 cm was calculated. During the
flight, a GNSS field survey was also carried out in order to ensure the
necessary spatial accuracy of the resultant images. Eight artificial
targets, 50
The GNSS survey was carried out in real-time kinematic (RTK), using geodetic receivers. A reference station was set up, recording continuous signals from the GNSS satellite constellation for more than 3 h. The positional information obtained by the reference station was then sent to a mobile receiver, using radio modem communication. Each GCP was occupied for at least 2 min with a recording interval equal to 1 s. The reference station coordinates determined using this technique were corrected by post-processing procedures using contemporary data recorded by three permanent GNSS stations (La Spezia, Pieve Fosciana and Pisa) allowing millimetric accuracy. Consequently, the rover data were shifted permitting centimetric accuracy for the eight artificial targets. The orthometric heights were also calculated by using Convergo, an Italian code for full coordinate conversion. The coordinates of the GCPs were collected in ETRF2000 and then converted to the Italian National Gauss Boaga system.
The frontal survey, with direction of acquisition parallel to the rock face, was carried out manually, without the use of the Aibotix™ AiProFlight planning software. Six flights were needed in order to cover all sectors of the mining area, resulting in a total of 448 overlapping images. The flights were executed in sub-parallel straight lines approximately 60 m distant from the rock face (Fig. 5), providing an average estimated GSD of 1.5 cm.
Plan view
As for the zenithal flight, a series of GCPs and check points (21 targets in total; Fig. 6) were measured using a reflectorless TS. Their location was decided considering a balance between an optimum spatial distribution (Fig. 6), both in space and elevation from the open-pit floor, and easy identification of points on the images. Due to the complex morphology of the slopes and the extent of the mining area, a large number of GCPs were used to orient the photogrammetric model. Two GNSS receivers, operating in static modality, were used to obtain the geographic coordinates of two points: the origin of the survey and its zero-azimuth direction. For this survey, GNSS data were corrected using contemporary data recorded by permanent GNSS stations and ellipsoidal heights were converted to orthometric heights.
Topographic survey with TS
The software Agisoft™ PhotoScan Professional version 1.2.5 (Agisoft, 2016) was used to process the images obtained from the two RPAS surveys (two zenithal flights plus six frontal flights). This software is capable of solving the camera interior and exterior orientation parameters and generating georeferenced spatial data such as 3-D point clouds, digital surface models and orthophotos. All the images acquired from the two surveys were processed with identical photogrammetric processing, in two distinct Agisoft™ PhotoScan projects; one for the zenithal flights and another for the frontal flights.
The first processing step consisted of image alignment, through which the interior and relative orientation parameters were solved. In order to improve the whole alignment process and to obtain low re-projection error, millions of tie points were automatically extracted without setting a point limit. Following image alignment, the second processing step involved georeferencing of the 3-D model in such a way as to solve the exterior orientation parameters by using the GCPs coordinates measured during the two GNSS–TS topographic surveys. The exterior orientation of images was necessary to measure the orientation, with respect to the north, and inclination of slopes and discontinuities, which are needed for the stability analysis (Firpo et al., 2011). For both surveys, a number of the measured points were used as check points to verify the model accuracy. Specifically, for the zenithal survey 2 of the 8 measured targets were used as check points, while for the frontal survey (with direction of photo acquisition parallel to the rock faces) 4 points out of 21 were used as check. Both natural and artificial targets were identified directly on the images, assigning a 3-D coordinate to each of them.
Subsequently, the “optimize” tool was utilized to adjust the estimated camera positions and to remove possible non-linear deformations, minimizing the errors due to re-projection and misalignment of the photos. Moreover, the optimization was improved by deleting all the tie points with a re-projection error greater than 1 pixel.
In a subsequent step, the zenithal and frontal dense 3-D point clouds were generated with medium quality and aggressive depth filtering settings. No automatic classification of clouds was necessary: no infrastructure was present as the mine was not operational and there was no vegetation within the area of interest.
Finally, a polygonal 3-D mesh model was created from the point cloud and used
to create the orthophoto of the open-pit area. The orthophoto has image
distortions removed due to camera characteristics (i.e. lens distortions),
camera tilt and topographic relief displacement. Unlike an uncorrected
aerial photograph with a perspective projection, an orthophoto is
geometrically corrected (“orthorectified”) and can be used to measure true
distances since it is “scale-corrected”. The corrected orthophoto image with
a spatial resolution of 1 cm pixel
In order to characterize the rock mass within the open-pit mine, a relatively large number of discontinuities were identified directly on the dense point cloud. The orientation of the selected discontinuities was manually calculated by creating patches that best fit the identified discontinuity planes in the point cloud and extracting their orientation using the Leica™ Cyclone 9.0 software. Discontinuity sets were then identified using stereographic representation (Schmidt equal-area method, lower hemisphere).
According to Mastrorocco et al. (2017), a manual deterministic fracture
mapping was also adopted because it increases the level of control of the
process, which is essential where the morphology of the quarry slope surfaces
is largely artificial (smooth cut surfaces). The collected point
cloud-derived data were then compared with those manually measured by
traditional engineering geological surveys. On the basis of engineering
geological data, geological strength index (GSI; Hoek, 1994) and rock mass
rating (RMR; Bieniawski, 1989) characterization were also applied, and a
kinematic stability analysis was carried out using the Markland test
(Markland, 1972). The latter testing was undertaken in order to identify
potential kinematic failure mechanisms that affect the slopes. The tests for
planar sliding, wedge sliding and direct toppling were conducted for the two
principal slopes (eastern slope dip direction / dip
Despite the importance of performing kinematic analysis to discover possible block instability, one of the principal limitations of this stereographic method is the inability to locate the block source areas on the slope being analysed. For this reason, the most critical blocks have been identified directly on the point cloud. The points representing the geometry of every single block were meshed in Leica™ Cyclone 9.0 and their volume estimated using reference planes corresponding to the discontinuities that demarcate or shape the respective blocks.
The collected orientation data were finally used for preliminary stability analyses using Rocscience™ Swedge software. Swedge is a 3-D software for evaluating the stability of surface wedges in rock slopes. It considers the intersection of discontinuities and allows the calculation of safety factors of the formed blocks. The software is based on classical limit equilibrium methods that usually have some limitations, such as the lack of consideration of in situ stress, strains and intact material failure (Stead et al., 2006). Nevertheless, Swedge does allow consideration of external forces by applying a force (vector with given orientation and intensity) to the formed blocks. From this, preliminary analysis of the impact of water content or other forces can also be performed. The software can also be used to assess the stability of wedges formed by a basal plane.
The image alignment process, described in Sect. 3.2, resulted in a re-projection error on the images of 0.41 pixel for the zenithal survey and 0.48 pixel for the frontal survey. The final root mean square error (RMSE) for the zenithal flights exterior orientation was equal to 0.042 m; RMSE for the frontal flights exterior orientation was equal to 0.043 m (Table 2).
Information related to the zenithal and frontal photogrammetric surveys and processing.
The final 3-D frontal and zenithal point clouds contain more than 18 000 000 and almost 13 000 000 points, respectively, with a mean point spacing varying from 1 to 4 cm.
A total of 154 discontinuity planes were selected on the 3-D point clouds derived from both frontal and zenithal surveys; related orientations were thus calculated taking into account geographic coordinates of the clouds and then plotted on a stereonet. Analysis of the stereonet allowed identification of four discontinuity sets, whose properties listed in Table 3 were obtained in the field through traditional engineering geological survey.
Characteristics of the discontinuity sets measured based on traditional manual engineering geological survey in the study area. The joint roughness coefficient (JRC) and joint compressive strength (JCS) are also indicated.
Stereonet plot of poles (Schmidt, equal area, lower hemisphere) of
the discontinuities manually collected through classical engineering
geological survey
By comparing the discontinuity planes obtained through traditional manual engineering geological survey (in Fig. 7a) and that identified on the point cloud (in Fig. 7b), a good level of congruence has been highlighted, confirming the quality of the taken approach.
Based on the discontinuity characteristics derived from RPAS and traditional
engineering geological surveys, the basic RMR (RMRb) and GSI were
calculated. The RMRb resulted equal to 67 (Table 4), while the GSI was
estimated to be between 60 and 65 using the modified chart proposed by Hoek
et al. (2013). In addition, application of Hoek et al. (2013) equation for
GSI quantification – GSI
Parameters used for RMRb determination. Chosen index values are bold.
Examples of kinematic stability analysis carried out using
Rocscience™ Dips 6.0 stereographic projection through the Wulff
equal-angle method (lower hemisphere):
These results agree with the authors' field observations and what is described in the actual quarry excavation plan (Lorenzoni, 2012). In view of the rock competency, potential instability is not related to the strength properties of the rock mass, but the intersection of discontinuity planes can locally form rock blocks with the potential for sliding or toppling. For this reason, a kinematic stability analysis was performed.
A discontinuity friction angle of 35
Potentially unstable discontinuity systems along the two different slope orientations.
Identification of two large blocks with potential for sliding; inset photo highlights the visible aperture of the basal plane.
Three different possible kinematic modes were identified, with K2b and K4 systems having the most influence on potential instability. The majority of the potential failures identified relate to planar sliding or wedge sliding, in agreement with field and SfM-based observations. The results highlight the potential for blocks to form that may be subject to gravity induced instability but, as previously stated, traditional kinematic analyses do not identify the location of these unstable blocks. Therefore, further analysis of the high-resolution images and the dense point clouds was performed in order to locate possible block source areas. More than 20 blocks were deterministically characterized in terms of size, shape and barycentric coordinates, varying from about a cubic metre to a few hundred cubic metres. In addition, the analysis identified two large blocks (Fig. 9), a few thousand cubic metres in size, with potential for sliding. These are formed by the intersection of two different faults and a discontinuity basal plane with 5 cm aperture, no infill, smooth surface and high persistence. The basal plane appears not to correspond with any of the identified discontinuity sets but is probably connected to planes of weakness of the marble in correspondence with a particular orientation of minerals crystallographic axes. The lateral and rear faults, however, may be associated with the K3a and K3b systems, respectively. The rear fault may also be associated with the E–W fault system that characterizes the geology of this area of the Apuan Alps complex (Fig. 2).
The geometric characteristics of the two blocks, including orientation of the intersecting discontinuities and volume of the meshed block, were obtained using Leica™ Cyclone 9.0 and are shown in Table 6.
Characteristics of identified blocks A and B.
The first block, Block A, is of particular interest because it daylights in the slope face and prevents Block B from sliding (similar to an active-passive wedge; Prandtl's prism transition zone; Kvapil and Clews, 1979). A main road access to the quarry is located at the base of this slope, increasing the potential risk for the area. Due to the particular geometrical configuration, Block A can be described as the key block as it daylights the rock face. In the actual setting Block B does not hold the potential for sliding as it does not daylight in the slope face, but it could play a significant role in terms of additional weight force. Nevertheless, the following preliminary stability analysis is focused on Block A. Further investigation would require an analysis of the effect of Block B on the potential for instability as this provides the “active” component of the active–passive wedge.
The geometry of Block A was deterministically re-created in Swedge using the
geometrical information shown in Fig. 10. Initially the discontinuities were
assumed to be fully persistent. This is a common approach in engineering
geology, since reliable values of persistence are almost impossible to
obtain from field mapping and most rock slope stability analyses assume that
the 100 % persistent joint exists on failure surface (Park, 2005). The
basal plane is a joint and the possible presence of rock bridges should be
carefully considered. In this case the basal plane did not show the presence
of segments of intact rock along its trace on the rock exposure;
consequently, application of Eq. (2) confirms a 100 % persistence of the
basal plane, which was used in the first analysis. In addition, the lateral
and rear or back surfaces correspond to geological faults and can be
therefore considered fully persistent. In this regard, the western lateral
surface has been included in the Swedge model has a discontinuity (Joint 1
in Fig. 10; dip direction / dip
Result of Rocscience™ Swedge 3-D preliminary slope stability analysis.
The adopted limit equilibrium solution for the slope stability analysis was
based on the Mohr–Coulomb shear strength model with a friction angle of
35
The result highlights a possible condition of instability for Block A which does not match with field observations, since the block under study has remained stable in this position for tens of years. In order to investigate the effect of uncertainty or variability of the input parameters, a sensitivity analysis was performed. In a sensitivity analysis (Fig. 11) specific parameters are varied across a range of values and the effect on the factor of safety is observed. This helps to identify the parameters that have the most effect on block stability. Since the geometrical inputs are well defined through the accurate 3-D model, the subsequent analysis focused on waviness angle (it accounts for the undulations of the joint surface, observed over distances on the order of 1 to 10 m; Miller, 1988), cohesion and friction angle of the basal plane, and water pressure. These are also the parameters with the higher input uncertainty.
Result of sensitivity analysis relative to the parameters basal joint waviness, friction angle, cohesion and water percent filled.
Results of the sensitivity analysis (Fig. 11) show that the cohesion is
clearly the parameter that has the most effect on block stability. For this
reason, the effect of this parameter was investigated in more detail in the
following analyses. Application of Eq. (3) makes it possible to determine an
effective cohesion dependent on the continuity of jointing. From this the
contribution of eventual rock bridges on the block stability can be
investigated starting from intact rock cohesion material value, which has
been determined to be approximately 16 MPa (Lorenzoni, 2012; Dumas, 1999;
data from Geomechanical Laboratory of the Centre of GeoTechnologies of Siena
University). Table 7 shows the results in terms of factor of safety obtained
from a parametric stability analysis performed with increased values of
effective cohesion, corresponding to 0, 0.5, 1, 2, 5 and 10 % of rock
bridges on the basal plane (total area of 510 m
Results of the parametric analysis increasing basal plane cohesion values.
The RPAS approach adopted in this case study, based on the combined use of high-resolution images from different perspective and accurate GNSS–TS topographic surveys, overcame data acquisition difficulties related to high steep quarry walls and provided high-resolution orthophotos of the site (1 cm pixel size). The application of RPAS instrumentation was extremely successful for reconstructing the complex morphology of the mine site where ground-based techniques (e.g. TLS, terrestrial photogrammetry) have limitations due to potential “shadow” effects and several inaccessible set-up zones due to safety reasons. GCPs measured using a TS and GNSS receivers permitted a high level of accuracy in the images exterior orientation, which is particularly important for subsequent discontinuity measurements. In particular, the final RMSE calculated on the check points (Table 2) was about 6 cm for the zenithal flight and 3 cm for the frontal flight. This small difference is mainly due to the fact that higher parts of the quarry were inaccessible for safety reason, and GCPs of the zenithal flight were only located at the bottom of the cliff, slightly increasing the RMSE of related model (0.065 m). Differently, in the case of the frontal flights performed perpendicularly to the rock faces, GCPs were collected at different heights using a TS, obtaining a frontal model of the cliff characterized by a lower error (0.03 m).
It should be underlined that GCPs for both zenithal and frontal flights are projected in the same reference system. In the end, analysis of the results confirms the good accuracy level of both models, widely adequate for the purpose of the work.
Possible limitations in the use of RPAS system can be related to the need for a pilot license and user experience on topographic survey and imagery processing. Indeed, the accuracy of the final 3-D model can be greatly affected by the quality of data collected (photos and GCPs), hardware and software capability and user expertise. Nevertheless, in this work the 3-D models have been obtained using a medium-low performance computer (Intel i5 CPU at 3.20 GHz with 16 GB RAM), using images obtained from a light compact digital camera that can easily be mounted on a low-cost RPAS. In addition, the use of GCPs overcame the necessity of an expensive IMU system for accurate image alignment. This confirms the reason of the widespread of RPAS for engineering geological investigation, mainly due to its low cost, speed of data acquisition and high safety.
In this work, the elaboration speed in terms of time employed to identify and characterize rock discontinuities has been partially decreased to guarantee consistency and quality in the interpretation of these elements. In fact, although several authors have demonstrated the reliability of automatic and semi-automatic processing of imagery and 3-D point clouds for fracture mapping (Mah et al., 2011; Vöge et al., 2013; Assali et al., 2014; Vasuki et al., 2014), a complete manual approach was adopted in this analysis because in most cases the flat and regular morphology of quarry walls only allow photointerpretation of discontinuity traces. Moreover, final visual inspection and validation of outputs is always required, even when using codes for the semi-automatic extraction of joints (Salvini et al., 2017).
In this work, the orientation of several discontinuity planes was calculated using Leica™ Cyclone 9.0 on the point cloud once its high positional accuracy level was demonstrated. This allowed for a more complete characterization of the rock mass than that which could be obtained through traditional engineering geological survey, due to limited safe access to the slopes within the site. This was particularly important also because discontinuity characteristics can vary at different heights of the rock mass due to stress relaxation induced by excavation activity. In this context, the possibility to inspect the mining area from different angles in high definition allowed identification of critical areas to be analysed in detail for safety purposes. Moreover, the possibility to use the point cloud for obtaining geometrical characteristics of blocks represented a major advantage because it allowed the correct geometrical reconstruction of a 3-D model to be used in specific software for slope stability analysis.
Characterization of the site highlighted a potential significant risk for the future workforce due to the presence of two major blocks with potential for sliding. With conservative assumptions the preliminary limit equilibrium analysis showed that the key Block A, in its present shape, is potentially unstable. This is mainly due to the fact that the basal plane dips out of the slope and daylights on the face, with a dip angle higher than the friction angle of the surface. The block is separated from the rock mass by a major fault, which can be simulated in the Swedge analysis as a tension crack. The fault can be clearly identified from the orthophoto obtained from the application of SfM method. An apparent offset has been identified on the back fault (Fig. 12). The presence of cataclasite with variable thickness is also an indication of potential for offset on the fault surface.
Detail of the zenithal orthophoto with indication of apparent offset of a major fault acting as back release surface for Block A of Fig. 9.
Details of a series of tight discontinuities above and below the basal plane under study.
In the preliminary limit equilibrium analysis, the major uncertainty is on
the basal plane (Fig. 12). Despite the presence of a continuous trace line
on the rock exposure, its full persistence in the rock mass is not clear. In
general, the presence of rock bridges plays an important role in stabilizing
removable rock blocks. Indeed, a rock block cannot fall or slide from a
slope until the rock bridges have failed. Rock bridge failure involves
fracturing of intact rock, which can be an order of magnitude stronger than
the rock mass (Kemeny and Donovan, 2005). From the undertaken sensitivity
analysis, cohesion of the basal plane was the parameter that had the most
significant influence on the block's factor of safety. This suggests that it
may not be completely persistent in this case, since the block has remained
stable over time. In this regard, the parametric analysis carried out
increasing the cohesion values showed how a small rock bridge, corresponding
to 1 % of the basal plane surface (5.1 m
The calculated rock bridge percentage seems in contradiction with the author's field observations and their experience in similar contexts, which suggest that a higher percentage of rock bridges may exist. Hudson and Priest (1983) identified two kinds of persistence relative to impersistent or intermittent joints that should be considered. Differently from impersistent joints, intermittent discontinuities consist of joint segments and intact rock bridges on the same plane. Mauldon (1994) claims that the formation of intermittent joints is geologically unlikely, unless a preferential direction of weakness exists within the rock mass. In this case, the cohesion of rock bridges along intermittent joints could be much lower than that of the intact rock, and consequently the percentage of rock bridges could be greater. This could be the case of Block A, and the presence of a series of discontinuities with similar dip direction and dip to the basal plane seems to confirm the hypothesis of a preferential plane of weakness due to the geomechanical characteristics of the marble material in that portion of the mining area (Fig. 13).
The progressive degradation of rock bridge elements with time could cause a progressive failure mechanism that has the potential to lead to a final rockfall event. This is particularly important in small engineered slopes such as the present one, where the rock mass may be continuously disturbed by excavation activity driving the slope to instability. Such mechanisms of progressive brittle fracturing of rock bridges are not considered in limit equilibrium approaches, where a small content of rock bridges adds significant apparent cohesion to the failure surface (Elmo et al., 2011; Tuckey and Stead, 2016). Therefore, in case of re-opening of mining activities, an in-depth engineering geological analysis together with the installation of a monitoring system for observing the behaviour of the rock mass over time should be considered.
The case study highlights the powerful use of RPAS technology for rock slope characterization and acquisition of accurate 3-D data for subsequent stability analysis. Specifically, an Aibotix™ Aibot X6 six-rotor multicopter was employed to obtain high-resolution topographic data of a blocky rock mass located within a quarry prone to discontinuity-controlled instability mechanisms. A detailed 3-D model of the area allowed accurate identification and geometrical measurement of the geological discontinuities that isolate significant volumes of rock. The stability analysis performed with Rocscience™ Swedge software showed that rock bridges can have a significant influence on stability conditions. The analysis highlighted the need for further detailed analysis and installation of suitable monitoring systems for future quarry operations.
These results confirm the reliability of the employed technologies to provide data for preliminary evaluation of the hazard affecting the study area. The RPAS allowed acquisition of high-resolution topographic data in an area characterized by a complex morphology where ground-based techniques would have significant limitations (e.g. TLS, terrestrial photogrammetry). It is worth noting that in mountainous environments, the use of RPAS has to be evaluated according to the local atmospheric and topographic conditions. The high temporal and spatial variability of the atmospheric conditions at high altitudes, as well as the presence of vegetation or steep and irregular slopes, could endanger the flight operations. This requires pilots with relevant experience and RPAS equipped with innovative systems to manage emergency conditions. Future analysis at this site will concentrate on the evaluation of the most useful countermeasures to reduce the risk conditions by monitoring the unstable slopes and undertaking further stability analysis including more complex 3-D discrete fracture network evaluation to assess the effects of rock bridges and elastoplastic numerical approaches to assess likely instability.
Data come from a research project between University of Siena and USL1 of Massa and Carrara (Mining Engineering Operative Unit, Department of Prevention) related to slope stability analyses in marble quarries. Data are in the public domain but authors do not have permission to make them publicly accessible via the internet.
The authors declare that they have no conflict of interest.
This article is part of the special issue “The use of remotely piloted aircraft systems (RPAS) in monitoring applications and management of natural hazards”. It is a result of the EGU General Assembly 2016, Vienna, Austria, 17–22 April 2016.
Part of the present study was undertaken within the framework of an agreement with USL1 of Massa and Carrara (Mining Engineering Operative Unit, Department of Prevention). The authors acknowledge Maura Pellegri and Domenico Gullì (USL1) and Vinicio Lorenzoni (Professional geologist) for their support of this research. Edited by: Yuichi S. Hayakawa Reviewed by: two anonymous referees