Next Article in Journal
Determiners in the Consumer’s Purchase Decision Process in Ecotourism Contexts: A Portuguese Case Study
Next Article in Special Issue
Reliability-Based Design of Protection Net Fences: Influence of Rockfall Uncertainties through a Statistical Analysis
Previous Article in Journal
Assessment of Geotechnical Seismic Isolation (GSI) as a Mitigation Technique for Seismic Hazard Events
Previous Article in Special Issue
Rockfall Hazard Assessment in Volcanic Regions Based on ISVS and IRVS Geomechanical Indices
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

In Situ Block Size Distribution Aimed at the Choice of the Design Block for Rockfall Barriers Design: A Case Study along Gardesana Road

by
Gessica Umili
1,*,
Sabrina Maria Rita Bonetto
1,
Pietro Mosca
2,
Federico Vagnon
1 and
Anna Maria Ferrero
1
1
Department of Earth Sciences, Università degli Studi di Torino, via Valperga Caluso 35, 10125 Torino, Italy
2
Institute of Geosciences and Earth Resources—CNR, via Valperga Caluso 35, 10125 Torino, Italy
*
Author to whom correspondence should be addressed.
Geosciences 2020, 10(6), 223; https://doi.org/10.3390/geosciences10060223
Submission received: 18 May 2020 / Revised: 3 June 2020 / Accepted: 4 June 2020 / Published: 7 June 2020
(This article belongs to the Special Issue Rockfall Hazard)

Abstract

:
When studying rockfall phenomena, a single value of the block volume is not sufficient to take into account the natural variability of the geometrical features (orientation, spacing, persistence) of the discontinuity sets. Different approaches for obtaining cumulative distributions of potentially detachable block volumes are compared. A highly fractured rock mass outcropping along the western Lake Garda (Italy), consisting of prevailing limestone and interbedded marls, is studied in detail from geological and geostructural points of view. Then, a representative rock face has been selected and analyzed with traditional and non-contact survey methods to identify the main discontinuity sets and to collect spacing samples. Based on these data, in situ block size distributions for different combinations of sets are built following statistically-based approaches, without the use of a Discrete Fracture Network (DFN) generator. The validation of the obtained distributions is attempted based on the detached block surveyed at the foot of the slope. However, in this particular case study, the detached blocks cover only a minimal volume range compared to both theoretical values and visible rockfall scars. The fallen rock blocks have a marginal role in design block determination, since their volume depends on geological discontinuities (bedding and fractures) and could be affected by other processes after the detachment (e.g., fragmentation). The procedure here described should be standard practice in the study of rockfall events, and it should be uniform in European standards such as Eurocodes. Future developments should involve the scientific community for setting the percentiles of the probability distribution to be considered for block design definition.

1. Introduction

The rockfall phenomenon is defined as the detachment, fall, rolling, and bouncing of rock fragments, which can break during impacts [1]. Events range from small pebbles falling down a cliff to massive boulders of several tonnes rolling down a slope [2]. Rockfalls represent one of the major causes of fatalities associated with landslides due to the high energy and mobility of the involved rock blocks [3]. Rockfall modeling aims to assess the envelope of trajectories, the maximum runout distance, the distribution of kinematic parameters along a fall path, and the probability for a specified “design block” to stop at specific distances from the starting point [4].
However, a single value for the design block is not sufficient to take into account the natural variability of the geometric features (orientation, spacing, persistence) of the discontinuity sets. A cumulative distribution of block volumes that could be potentially generated by discontinuity sets in a rock mass is the correct approach.
The prediction of in situ block size distribution (IBSD), or better, its assessment, since the true IBSD can rarely be evaluated, has been one of the primary pursuits of mining and quarrying operations [5]. Size and shape determinations of in situ rock fragments are becoming increasingly important issues in materials mining, comminution and handling industries, and in the construction industry using stone products [6].
Many methodologies were developed to forecast IBSD. For instance, Wang proposed “the Equation Method”, which is based on a formula and a series of look-up tables [7], to define the cumulative curve for the in situ block sizes, considering negative exponential, log-normal, and uniform distributions of spacing values. The equation method uses a set of empirical equations to estimate the IBSD, relating it to the mean spacings and the mean orientations of the three primary sets of discontinuities. The equations are obtained using the exact dissection method solutions described in [8], based on the block theory by Goodman and Shi [9].
Elmouttie and Poropat [10] proposed a new method to predict IBSD in fractured rock masses using realistic DFN, robust polyhedral modeling, and a Monte Carlo sampling approach, to account for the stochastic variability in the fracture geometry. This approach is promising but not directly applicable to survey data without the use of a DFN generator.
Stavropoulou [11], following the methodology proposed originally by Priest and Hudson [12], introduced the closed-form expression of block volume cumulative distribution based on joint frequencies measured along drilled cores. This method initially developed for the quarry environment can hardly be transposed to natural slopes.
In 2006, Latham and coauthors [13] stated that, due to the rapid advances in digital imaging and semi-automated geological face mapping, software packages devoted to rock mechanics applications might soon enable IBSD to be rapidly determined from software, with a limited manual mapping of joints. Nowadays, non-contact survey methods have indeed become more robust, producing reliable joint orientation [14,15,16,17,18,19,20,21,22,23,24,25,26], trace length [27,28,29], and spacing data [30,31]. However, IBSD from survey data is not yet an option in rock mechanics codes for the study of rockfall phenomena, even though IBSD could provide useful tools for the design of flexible barriers.
Buyer and coauthors [32] presented a consistent and semi-automatic process for rock mass characterization, starting from photogrammetry and ending in numerical modeling. The Discontinuity Set Extractor [23], a semi-automatic method for identifying and analyzing the surfaces of rock outcrops, is used in combination with the software SMX Analyst (3GSM GmbH). Then, the joint network characteristics are used in 3DEC (The Itasca Consulting Group) to determine the IBSD and shape distribution at a very high level of detail and objectiveness. However, this level of automation is far from available to engineers who deal with the design of rockfall protection barriers.
In order to follow the requirements of the Limit State Design (LSD), as defined by Eurocode 7 (EC7) [33], survey methods should guarantee a high level of accuracy and precision of the results to reduce the uncertainties of the estimate of the design parameters. In essence, this is particularly relevant for the geometric characteristics of discontinuities, since the dimension and shape of potentially detachable blocks depend on their values [34].
Ferrero and coauthors [35] proposed a simple application of interval analysis to spacing data obtained from non-contact surveys for block volume estimation to quantify the influence of uncertainty on spacing and relative orientation. Since in rockfall phenomena the involved kinetic energy is strongly related to the block mass, the structural capacity of flexible barriers is defined based on the characteristic block volume, which at present is often subjectively decided by the designer. For example, EC7 [33] gives no thorough information concerning the choice of the design block; only national standards, such as UNI 11211-4 in Italy [36] and ONR 24810 in Austria [37], provide suggestions on its definition. In Italy, the characteristic block volume is defined based on the experience and knowledge acquired after an onsite survey. Conversely, ONR standards propose a percentile of the IBSD in relation to three consequence classes defined in EN 1990: 2002 [38]. However, the commonly assumed 95% fractile (or higher) value for IBSD would be very restrictive and improbable for many rockfall events [39].
A statistically robust IBSD represents a promising tool for a more practical choice of the design block, mainly because it covers the entire volume range. Therefore, the designer could take advantage of IBSD for comparing different scenarios concerning the technical features of the rockfall barrier, locations, and costs.
This work aims to compare different methods to infer IBSD, considering a rock mass outcrop in the Gargnano-Muslone area along the western Lake Garda (Brescia province, Italy; Figure 1). The investigated outcrops, consisting of prevailing limestone and interbedded marls, form steep rock faces along the main (Gardesana Road) and secondary roads and record continuous rockfall phenomena due to their high degree of fracturing.
This work also aims to bring attention to the possible limitations of the application of methodologies for estimating the return period of rockfall blocks according to their size (e.g., [40]). The survey of blocks accumulated at the slope foot and the available historical data could be insufficient, particularly in rock masses characterized by a complex geostructure as the one described here. Moreover, during impact with the ground, the detached rockfall mass can undergo fragmentation along pre-existing discontinuity planes, by the propagation and coalescence of discontinuities into rock bridges, or by the development of new cracks in intact rock [41]. The validation of the obtained curves is attempted, showing the limited contribution of data acquired at the foot of the slope in this particular case.

2. Geological Setting

The Lake Garda region belongs to the Southern Alps (Figure 1A), a structural domain of the Alps, and the western side of the lake consists of Mesozoic lithologies of the Lombardian Basin.
The Southern Alps identify a south-vergent fold-and-thrust belt separated from the north-vergent collisional belt by the Periadriatic Lineament and preserve part of the Mesozoic continental margin of the Adriatic plate (e.g., [42,43]). The Lombardian Basin developed as a complex subsiding area during Early Jurassic on the previously weakened crust and was separated to the east from the Venetian Platform by the extensional Ballino-Garda fault [44]. During the Alpine orogenesis, the Lombardian Basin succession was deformed by the southward propagation of the Southern Alps belt [45].
In the present structural framework of the Southern Alps, the western Lake Garda lies between the Tremosine–Tignale thrust system and the inverted Ballino–Garda fault (Figure 1B). The exposed sedimentary succession is part of the Lombardian Basin, which is tectonically overlain by Norian platform carbonates to the west [46,47,48]. Massive platform limestones of the Corna Formation (Rhaetian–Early Jurassic p.p.) are upsection followed by Lower Jurassic pelagic carbonates (mainly cherty limestone and marly limestones) of the Medolo Formation. The Middle Jurassic Concesio Formation (marlstone and marly limestone) is overlain by siliceous sediments of the Selcifero Lombardo Formation of Middle-Late Jurassic age, which is in turn covered by thinly bedded limestones of the Maiolica Formation (Early Cretaceous in age). The upper stratigraphic levels consist of very thinly stratified calcareous and argillaceous marls of the Scaglia Lombarda Formation (Late Cretaceous—Early Eocene).
Along the western Lake Garda, the whole Lombardian Basin succession is folded by km-scale anticlines and synclines (Figure 1B), and in the investigated Gargnano-Muslone area, the core of the main syncline consists of the Scaglia Lombarda Formation [48]. In this area, open to tight asymmetrical folds (Figure 2A) have mainly NE–SW trending axes and axial planes variably dipping toward NW.
In this geological setting, the analyzed rock mass is part of the Medolo Formation (Figure 2B,C), defining the southern limb of the Gargnano-Muslone syncline. This formation is mainly made of marly limestone and cherty limestone forming 10–40 cm thick beds and thin-bedded gray marls characterized by a thickness usually ranging from a few centimeters to 30–40 cm but locally significantly thicker. The bedding of the investigated outcrops dips toward NW and W at a moderate angle. In the area, mesoscopic compressional structures (top-to-SE) subparallel to the bedding have been observed.

3. Methods

This section focuses on the block volume assessment. It aims to highlight the transition from the calculation of a deterministic value to the construction of different IBSDs based on the statistical processing of spacing samples. Survey methods are also included to complete the description.

3.1. Block Volume Assessment

Miles [50] proposed analytical equations to calculate the mean and variance of volume V of the blocks, which are created either by three discontinuity sets with a negative exponential spacing probability distribution (regular subdivision) or by a random space partition with planes placed by a Poisson process with uniform density (random subdivision). The comparison between the two subdivisions shows that the random one generates volume that is on average bigger than those created by the regular subdivision (ratio of about 2).
Palmstrøm [51] provided guidelines for estimating the block volume from various types of joint density measurements. His well-known equation in the case of 3D measurements of three discontinuity sets generating a block is:
V B = S 1 S 2 S 3 sin γ 12 sin γ 23 sin γ 13 ,
where S1, S2, and S3 are the spacing of the three sets of discontinuities; γ12 is the angle between set 1 and set 2, and similarly for γ23 and γ13.
These methods have some limitations: first of all, input values (frequency or spacing of discontinuity sets) are intended as representative values, and this often results in a rough calculation of the average, without considering the actual statistical distribution of the sampled data; besides, the output is a deterministic value, and this is far from representative of a complex medium as a rock mass.

3.2. IBSD Assessment

3.2.1. Method Based on Empirical Coefficients

IBSD constitutes an evolution of block volume calculation: it represents the cumulative curve of the potentially detachable in situ blocks, and its construction considers the frequency distributions of spacing values. A particularly simple method for IBSD assessment was proposed by Lu and Latham [5]. They revised the equation method [7] and added some refinements, including the fractal spacing distribution and the influence of persistence. They obtained the following equation, describing the volume of blocks generated by three sets of discontinuities:
V i , p = C i , p ( F i m p ) q S p m 1 S p m 2 S p m 3 c o s θ c o s ϕ c o s α i = 10 ,   20 ,   , 100 ,
where Spm1, Spm2, and Spm3 are the principal mean spacing values, namely the average true spacing values of the three sets; θ, ϕ and α are the angles between the mean orientations of the three sets; and q is a constant less than 1, whose value can be assumed between 1/5 and 1/2 [52]. One should note that since the angle between a pair of sets can theoretically vary from 0° to 180°, its cosine assumes a negative value for an angle >90°: however, since cos(π − x) = −cos(x), it is sufficient to consider always the absolute value of the cosine.
The relative “impersistence factor” Fimp [52] is defined as follows:
F i m p = { S D S r S D < S r 1    S D S r ,
where SD is the mean discontinuity size. It can be estimated under different assumptions about the discontinuity shape (circular, elliptical) and using different methods. In particular, methods for estimating the mean discontinuity diameter E(D) based on the mean trace length E(l) provide formulations that can be applied to trace lengths sampled in circular windows (a review of these methods is in [53]). Sr represents the characteristic size of the rock mass under consideration, i.e., the cube root of the volume of the rock mass of interest.
Ci,p is an empirical coefficient: values for discontinuity sets with negative exponential and uniform distributions, as well as with a log-normal distribution spacing law [54], are summarized in Table 1. Lu and Latham [5] added the values for a rock mass for which the three sets of discontinuities have fractal spacing distributions (Table 1). The range is the standard deviation multiplied by a constant of 1.64. The 90% confidence interval for Ci,p can be derived from Ci,p ± Range.
In [13], a practical step-by-step methodology for IBSD assessment is reported: it includes approaches that do not rely on specialized computer software. It suggests using actual discontinuity data from traditional and non-contact surveys to apply Wang’s equation method.
  • Determine the principal spacing for the three different sets.
  • Determine the three characteristic angles α, β, and ϕ.
  • Plot the principle spacing data as histograms and assign the appropriate principal spacing distribution. Hence, determine the principal mean spacing.
Apply Equation (2) using values from Table 1 for assessing the IBSD.

3.2.2. Method Based on the Cumulative Distribution Function

As discussed in [11], the spacing of a discontinuity set can be considered a continuous random variable; therefore, its cumulative distribution function (CDF), denoted as F(x) = P{X ≤ x}, defines the probability that a given spacing value X is less than x.
The spacing data of each discontinuity set were processed using a Matlab Statistics Toolbox: gamma, negative exponential, log-normal, and Weibull distribution functions were fitted to the data. Their goodness of fit was evaluated through the one-sample Kolmogorov–Smirnov (K–S) test.
The K–S test is a nonparametric test comparing the population CDF of the data and the hypothesized CDF. The null hypothesis is that the data in vector x come from a specific distribution. The test accepts the null hypothesis if the p-value is higher than the significance level, here, it is assumed to be equal to 1%.
Once the best fitting distributions were identified for each discontinuity set, the obtained CDFs were input in Equation (1) instead of single spacing values. Regarding the angles between the considered discontinuity sets, their average values were assumed as the characteristic values, since their influence on volume variability was proven to be limited, compared to that of spacing [35]. The result is the block volume CDF, which in the following will be called Sample CDF (SCDF). Since the equation considers a group of three sets generating a block, different terms were analyzed.

3.2.3. Method Based on Monte Carlo Simulation

Based on the best fitting distribution previously obtained for a set, 100,000 input spacing values were randomly generated from that distribution, and the corresponding CDF was obtained. The same was done for all the discontinuity sets. As for the measured spacing samples, CDFs were input into Equation (1) for calculating the block volume CDF. In the following, this will be called Monte Carlo CDF (MCDF).

3.3. Surveys for Block Volume Parameters Definition

Due to the morphology of the slope, most of the outcrops were not directly accessible. Therefore, non-contact techniques and traditional geomechanical surveys were applied to define the main parameters for block volume assessment, as described in the previous sections.

3.3.1. Non-Contact Survey

Non-contact survey techniques operate on the representation of the area of study through a 3D model made of one or more dense point clouds. In this study, point clouds were obtained by both laser scanner and photogrammetric surveys, which were performed by IMAGEO S.r.l.
The laser scanner Teledyne Optech Polaris LR was used first, setting an average point density of 1 point/cm. Due to the morphology of the rock slope, the laser scanner acquisition suffered from the physical restrictions in the instrument positioning. Therefore, a photogrammetric survey from a drone was also performed. The digital images were acquired with the integrated camera of the system DJI Phantom 4 (focal length 3.61 mm, resolution 4000 × 3000 pixels).
A total of 545 images were shot, with a ground resolution of 6.6 cm/pix, covering an area of about 0.6 km2. The shooting points were aligned horizontally and vertically in routes parallel to the slope; the camera was shot with nadiral and oblique inclination. The minimal overlapping was about 80% among images of the same longitudinal route and about 60% transversally among the routes.
Laser scans and digital images were georeferenced based on Ground Control Points (GCP). GCPs were acquired by GPS differential technique using GNSS Receivers Topcon Hiper V, referring to the GNSS permanent station of Rovereto (TN). Software Agisoft Metashape was used for performing point cloud reconstruction based on digital images. Then, point clouds were merged to create a Digital Surface Model (DSM) of the entire slope.
Software Rockscan [18] was used to manually delimit discontinuity planes on the images and estimate their dip and dip direction. This method allowed surveying 1580 discontinuity planes from which the main sets were identified.
Software CloudCompare was used to measure spacing among planes belonging to the main sets. Tools included in Compass plugin were used; the user manually drew a line perpendicular to the discontinuity traces belonging to the same set and measured their relative distances, avoiding the need for corrections of the Terzaghi’s bias.

3.3.2. In Situ Traditional Survey

In the accessible area at the base of the slope, a geomechanical survey campaign was carried out according to the International Society for Rock Mechanics (ISRM) suggested methods [55]. Traditional surveys along scanlines were performed to cense all the discontinuities intersecting the scanlines and describe geometric and kinematic features, such as orientation, spacing, persistence, roughness, wall strength, aperture, degree of weathering, seepage, presence, and nature of filling.
This method made it possible to conduct the survey of 138 planes. The information on the type of discontinuity, number, and orientation of the main sets helped in the validation of the non-contact survey data.

3.4. Blocks Volume Survey at the Foot of the Slope

The Rockfall Block Size Distribution (RBSD) introduced by Ruiz-Carulla and coauthors [56] refers to the rockfall deposit. Obtaining the RBSD of a large rockfall deposit may become a challenge due to the high number of blocks to be measured. The methodology was to lay a tape measure along the blocks deposit at the foot of the slope. By walking along the tape measure, a sequence of images of the deposit was shot: in each image, a portion of the tape measure was visible. It was used to set the scale of the image content utilizing ImageJ, which is an open-source image processing program for multidimensional image data, with a focus on scientific imaging. Then, in each image, we selected blocks with three visible edges and measured them: sizes of 394 blocks were acquired. Due to the heterogeneity of block shapes, two geometrical solids were assumed to calculate block volume: parallelepiped (Equation (4)) and ellipsoid (Equation (5)).
V p a r = L 1 L 2 L 3 ,
V e l l = 4 3 π d 12 2 d 13 2 d 23 2 = π d 12 d 13 d 23 6 ,
where L1,2,3 are the lengths of three different edges of the block, and d1,2,3 are the diagonals of three different faces, which are calculated as
d i j = L i 2 + L j 2 .
Surveyed volumes will be referred to as SVpar and SVell to distinguish between the two assumed shapes.

4. Application to the Case Study

Traditional and non-contact geomechanical surveys were carried out on a representative slope in the Medolo Formation (see the yellow dot in Figure 1) interested by rockfall events. The results are reported in the following sections.

4.1. Traditional Survey

Data collected along three scanlines arranged in correspondence of accessible areas at the foot of the slope indicate the presence of six main orientations for the discontinuities, which are schematically represented in the lower hemisphere plot of Figure 3.
Considering the features observed on the field and the geological knowledge of the area, they were grouped into five main sets (labeled as K) as follows (Figure 4).
K1 is the bedding. It dips mainly toward NW at a medium-high angle and can be southward-dipping on the short limb of asymmetric folds. It is persistent and regularly presents itself along the scanlines with centimetric to decimetric spacing.
Discontinuities dipping SE at high angles consist of fractures indicated as K2. They have been identified in all scanlines and have decimetric spacing.
Fractures E and W-dipping at a high angle were observed and associated with the same set; they are named K3a and K3b, respectively. At the outcrop scale, they are regularly present and show low persistence and centimetric to metric spacing.
Besides, two sets of fractures were occasionally surveyed, dipping toward NE and SW at a medium-high angles, named K4 and K5, respectively. They show high persistence.
Table 2 reports the average orientation and spacing range of the main discontinuities sets obtained from traditional surveys.

4.2. Non-Contact Survey

Due to the large extension of the rock face, the non-contact survey allowed for the acquisition of a large number of data, assuring higher representativeness of the measurements at the slope scale.
The rock face was divided into different windows, respectively named 3, 3A, 3B, 3C, 3D, 3E, 3F, and 3G (Figure 5). The DSM of each of them was input in Rockscan, along with the corresponding images. A high number of discontinuities were recognized and described in terms of orientation (dip/dip direction) and grouped into sets (Figure 6 and Table 3).
Table 4 reports the average orientation and spacing range of the main discontinuities sets obtained from non-contact surveys.
All sets identified by the non-contact surveys are coherent with the ones from the traditional scanline surveys. SE-dipping bedding surfaces (reported as K1b in Figure 6) due to macrofolds have been identified in sectors of the outcrops inaccessible by traditional field surveys.
The number of discontinuities of each set is extremely varied. Data belonging to K2 are very few, probably because of its orientation subparallel to the slope and/or to folded bedding. K5, despite its geometrical relationship with the face, results as the second less frequent set observed.
Spacing values among discontinuities belonging to the same set were then measured directly on the DSM, using Cloud Compare. Spacing values obtained from traditional surveys were finally added to the ones collected on the DSM and processed altogether to find the best-fitting frequency distribution for data of the entire slope.
Table 5 reports, for each set, the p-value of the distributions that were accepted by the K–S test applied to the spacing sample. In the case of more than one accepted distribution, the best performance corresponds to the highest p-value. Log-normal distribution was found to provide the best fitting of four out of five spacing samples. The number of data influences the test strongly: it is evident from the results that small samples lead to uncertainty in the identification of the best fitting distributions because the test is not able to reject the null hypothesis.
Table 6 reports sample arithmetic mean (SM) and distribution mean (DM) values of the total spacing sample of each set. Figure 7 shows the spacing CDF based on the best-fitting distribution.
Among the five sets, the four more realistic combinations of three sets generating a block were considered, namely the ones that consider the strong conditioning of the bedding, the persistence of the sets K2, K4, and K5, the frequency of the K3 and the mutual geometrical relationship with the rock face. They are K1-K2-K3, K1-K2-K4, K1-K2-K5, K1-K3-K5. For the sake of simplicity, they will be referred to as terns T123, T124, T125, and T135, respectively.
Figure 8 shows IBSD calculated according to Equation (2), using coefficients Ci,p of log-normal distribution (Table 1), and assuming fully persistent joints (Fimp = 1). In our study, spacing values from traditional scanlines have been corrected based on scanline orientation. At the same time, measurements performed directly on the DSM do not require corrections: therefore, DM values (Table 6) match principal mean spacing ones and are input in Equation (2).
Figure 9 and Figure 10 show the SCDF and MCDF, respectively, of the considered terns that create a block. Since spacing CDF of K2 and K5 are truncated at 90% (Figure 7), and at least one among them is involved in the considered terns, SCDFs are all truncated at 90% cumulative frequency. Instead, Monte Carlo simulations allow to infer the entire spacing CDFs. Therefore, MCDFs are complete.
In Table 7, volumes calculated with Equation (1) considering the two different mean spacing types in Table 6 are compared: SM values produce Vsm, while DM values produce Vdm. We also show their ratios on volume correspondent to 50% cumulative frequency in MCDF (VMCDF(50%)), to demonstrate that Equation (1), which is generally used considering mean spacing values, does not produce a mean volume, but a value that in this case is from 2.2 to 3.6 times bigger.
Figure 11, Figure 12, Figure 13 and Figure 14 show the comparison of block volume distributions (IBSD, SCDF, MCDF) obtained for the four considered terns. By definition, MCDF fits SCDF: they are both based on Equation (1) and the best-fitting distributions of spacing samples. The main utility of MCDF is to allow one to associate cumulative frequencies to volumes bigger than those covered by SDCF. Instead, IBSD nature is entirely different. In Equation (2), a unique value of spacing (DM) is used for each set of the considered tern. Therefore, a sort of “base volume” is calculated, and the coefficients of Table 1 are multiplied to the same “base volume” to build the IBSD. It can be observed that IBSD intersects SCDF and/or MCDF in all the four cases. The agreement among the curves is at its maximum for T123. For the other terns, IBSD tends to produce bigger volumes with the same cumulative frequency, until the intersection is reached.
Figure 15A represents one of the images used for measuring blocks at the foot of the slope with ImageJ methodology. Figure 15B shows the frequency of detached blocks (SV) for the two assumed shapes (Equations (4) and (5)).
Table 8 reports, for the considered terms, the ratios between the volume correspondent to 100% cumulative frequency in IBSD and the maximum surveyed volume (SVmax), and between the 95% cumulative frequency in MCDF (VMCDF(95%)) and SVmax, considering the two assumed shapes (parallelepiped and ellipsoid). The ratio is an indicator of the representativity of surveyed volumes with respect to the distribution obtained from spacing values. It is evident that in this case, blocks at the foot of the slope are some order of magnitude smaller than a plausible theoretical maximum value.

5. Discussion

A few main considerations derive from the application of traditional and non-contact methods to the same rockface.
The ranges of spacing and persistence values estimated for the same set with traditional and non-contact survey techniques are different (Table 2 and Table 4). This difference depends on the geological features of the analyzed rock succession and the limit of detection of the non-contact survey.
The bedding K1 identifies the most pervasive discontinuity, and its features depend on the depositional history of the sedimentary succession, leading to changes in composition, texture, and thickness from strata to strata. The K1 spacing defined directly on the outcrops by the traditional approach is generally lower than the values obtained by the application of the non-contact method. The field study of the lithological multilayer estimated average spacing values in the order of few centimeters, with marls usually from thin to very thin bedded and limestone beds only a few centimeters thick (Figure 16A).
In comparison, the non-contact surveys have provided an average value of about 20 cm. Thus, rocky volumes comprised between two K1 surfaces drawn by non-contact surveys can be actually more bedded. Besides, K1 spacing could be decreased by possible structural discontinuities subparallel to the bedding due to the regional tectonics.
Moreover, field observations show that pervasive K3 discontinuities, crossing the bedding surfaces, are significantly branched by the lithological multilayer consisting of limestone and marls with contrasting competence (Figure 2B and Figure 16B).
The spacing of the high angle discontinuity sets (K4 and K5), being of metric and plurimetric order, was better measured using CloudCompare directly on the DSM of the slope. In fact, in the traditional surveys, due to the spacing values that characterize them, they are poorly represented or not at all. They appear either as random planes or with too few measurements to be able to derive reliable mean values statistically.
K4 and K5 discontinuities are crucial in the rock face stability: vast portions of the rock mass have been observed to be mainly isolated by the intersection of these discontinuities with the other surfaces (Figure 16C). Thus, it is likely to assume that such large volumes tend to subdivide into smaller blocks already at the time of detachment (due to the above reasons), and then, they may further reduce in size as a result of impacts during the fall.
This fact would justify the small size (values below 0.25 m3) of the blocks regularly measured at the foot of the slope, which is smaller than the ones obtained by data from traditional and non-contact surveys. Nevertheless, the volumes of detached blocks are huge in comparison to the size of single blocks due to the influence and mutual interaction of persistent and high-spaced sets (K2-K4-K5).
The volumes of the blocks on the ground could correspond then to small unitary rocky volumes, initially forming large masses isolated by the discontinuities with high spacings. Probably, they were further fragmented by the impacts during the fall, especially given the marly component that reduces the strength of the intact rock, and by the occurrence of bedding surfaces, also at a minor scale.
In light of the above observations, a statistically robust block volume distribution assumes a fundamental role in covering the entire possible volume range and associating to each value a probability of not being exceeded. Each spacing sample associated with a discontinuity set should contain a statistically sufficient number of measurements and be representative of its variability in the considered rock mass. The larger the maximum spacing, the greater the considered area must be to allow enough measurements to be made. The use of non-contact survey methods is crucial in this sense. In any way, traditional surveys cannot be neglected: their sensibility in the identification of the peculiar features of discontinuities and the associated structural elements is irreplaceable.
A robust fitting process and goodness-of-fit evaluation are also fundamental for selecting the most suitable distributions for describing spacing samples and build the proper CDFs.
The comparison of the tested block volume cumulative distributions obtained for the four considered terns shows a good agreement. MCDF is an extension of SCDF: the more complete the SCDF, the more accurate the MCDF. On this basis, MCDF allows to associate cumulative frequencies to volumes bigger than those covered by SDCF. In general, with the same cumulative frequency, we obtained IBSD volumes bigger than MCDF ones until the intersection with MCDF is reached.
Block volumes surveyed at the foot of the slope, in this case, yield data regarding only the smallest possible blocks. Visual inspection of rockfall scars (Figure 16C) confirms this fact. Since just the first part of the IBSD and MCDF overlap SV (Table 8), it is impossible to perform a real validation of the obtained distributions. Moreover, it is not possible to indicate which is the contribution of each considered tern in the total block volume distribution. It is the opinion of the authors that the ongoing research regarding the creation of a connection between the IBSD and the RBSD, considering the fragmentation process [57,58], is of fundamental importance and requires contributions from the scientific community.

6. Conclusions

This paper deals with the definition of ISBD for a highly fractured rock mass outcropping along the western Lake Garda (Italy). Different methodologies for the definition of design block volumes were discussed, and their results were compared, highlighting strengths and limitations.
Data from in situ traditional and non-contact surveys were analyzed for identifying the main discontinuity sets and their characteristics in terms of orientation and spacing. Statistical analyses on spacing values were performed for evaluating the frequency distribution of each discontinuity set. Log-normal and Gamma distributions demonstrated capability of simulating frequency distribution for discontinuity set spacing. However, the procedure described in this paper is subordinate to data availability (in terms of quality and quantity of data): the combination of non-contact and traditional surveys allows for this type of analysis, and consequently, the statistical results can be considered robust and reliable.
One of the main findings of this paper is the huge discrepancy between RBSD and ISBD results: the volume of fallen blocks distributed along the base of the rock face is far smaller than that estimated from spacing analysis on the rock mass. Consequently, if an RBSD value was used for barrier design, it would lead to rough mistakes and an underestimation of impacting energy.
In highly fractured rock masses, such as those here described, the fallen rock blocks have a marginal role in design block determination, since their volume could be affected by other processes after the detachment (e.g., fragmentation). Conversely, stratigraphy and structure play a fundamental role, since they mirror the degree of fracturing of the rock face. However, an accurate field observation, i.e., stratigraphy and structures of the rock face, is imperative for supporting any approaches to study rockfall phenomena.
The procedure here described should be standard practice in the study of rockfall events, and it should be uniform in European standards such as Eurocodes. Future developments should involve the scientific community for setting the percentiles of the probability distribution to be considered for the design block definition.

Author Contributions

Conceptualization, G.U. and S.M.R.B.; Data curation, G.U., S.M.R.B., and P.M.; Formal analysis, G.U. and S.M.R.B.; Funding acquisition, A.M.F.; Investigation, P.M.; Methodology, G.U., S.M.R.B., and P.M.; Project administration, S.M.R.B. and A.M.F.; Resources, G.U., S.M.R.B. and P.M.; Software, G.U.; Supervision, A.M.F.; Validation, G.U., S.M.R.B., and P.M.; Visualization, G.U., S.M.R.B., and P.M.; Writing—original draft, G.U., S.M.R.B., P.M., and F.V.; Writing—review and editing, G.U. and F.V. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

The authors thank IMAGEO S.r.l, who performed the laser scanner and photogrammetric surveys and provided images and DSM used in this study. The authors also thank Manuel Bertholin, who contributed to this study as part of his Master’s Degree thesis.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hungr, O.; Leroueil, S.; Picarelli, L. The Varnes classification of landslide types, an update. Landslides 2014, 11, 167–194. [Google Scholar] [CrossRef]
  2. Lambert, C.; Thoeni, K.; Giacomini, A.; Casagrande, D.; Sloan, S. Rockfall hazard analysis from discrete fracture network modelling with finite persistence discontinuities. Rock Mech. Rock Eng. 2012, 45, 871–884. [Google Scholar] [CrossRef]
  3. Ferrero, A.M.; Migliazza, M.R.; Pirulli, M.; Umili, G. Some open issues on rockfall hazard Analysis in fractured rock mass: Problems and Prospects. Rock Mech. Rock Eng. 2016, 49, 3615–3629. [Google Scholar] [CrossRef]
  4. Agliardi, F.; Crosta, G.B. High resolution three-dimensional numerical modelling of rockfalls. Int. J. Rock Mech. Min. Sci. 2003, 40, 455–471. [Google Scholar] [CrossRef]
  5. Lu, P.; Latham, J.-P. Developments in the assessment of In-Situ block size distributions of rock masses. Rock Mech. Rock Eng. 1999, 32, 29–49. [Google Scholar] [CrossRef]
  6. Wang, L.G.; Yamashita, S.; Sugimoto, F.; Pan, C.; Tan, G. A methodology for predicting the In Situ size and shape distribution of rock blocks. Rock Mech. Rock Eng. 2003, 36, 121–142. [Google Scholar] [CrossRef]
  7. Wang, H.; Latham, J.-P.; Poole, A. In-Situ block size assessment from discontinuity spacing data. Int. J. Rock Mech. Min. Sci. Geomech. Abstr. 1993, 30, A106. [Google Scholar] [CrossRef]
  8. Wang, H.; Latham, J.-P.; Poole, A.B. Predictions of block size distribution for quarrying. Q. J. Eng. Geol. Hydrogeol. 1991, 24, 91–99. [Google Scholar] [CrossRef]
  9. Goodman, R.E.; Shi, G.H. Block Theory and Its Application to Rock Engineering; Prentice Hall Inc.: Englewood Cliffs, NJ, USA, 1985; ISBN 0-t3-078189-4. [Google Scholar]
  10. Elmouttie, M.K.; Poropat, G.V. A method to estimate In Situ block size distribution. Rock Mech. Rock Eng. 2012, 45, 401–407. [Google Scholar] [CrossRef]
  11. Stavropoulou, M. Discontinuity frequency and block volume distribution in rock masses. Int. J. Rock Mech. Min. Sci. 2014, 65, 62–74. [Google Scholar] [CrossRef]
  12. Priest, S.D.; Hudson, J.A. Estimation of discontinuity spacing and trace length using scanline surveys. Int. J. Rock Mech. Min. Sci. Geomech. Abstr. 1981, 18, 183–197. [Google Scholar] [CrossRef]
  13. Latham, J.-P.; Van Meulen, J.; Dupray, S. Prediction of In-Situ block size distributions with reference to armourstone for breakwaters. Eng. Geol. 2006, 86, 18–36. [Google Scholar] [CrossRef]
  14. Kemeny, J.; Donovan, J. Rock mass characterisation using LIDAR and automated point cloud processing. Gr. Eng. 2005, 38, 26–29. [Google Scholar]
  15. Trinks, I.; Clegg, P.; McCaffrey, K.; Jones, R.; Hobbs, R.; Holdsworth, B.; Holliman, N.; Imber, J.; Waggott, S.; Wilson, R. Mapping and analysing virtual outcrops. Vis. Geosci. 2005, 10, 13–19. [Google Scholar] [CrossRef]
  16. Slob, S.; van Knapen, B.; Hack, R.; Turner, K.; Kemeny, J. Method for automated discontinuity analysis of rock slopes with three-dimensional laser scanning. Transp. Res. Rec. J. Transp. Res. Board 2005, 1913, 187–194. [Google Scholar] [CrossRef]
  17. Haneberg, W.C. Using close range terrestrial digital photogrammetry for 3-D rock slope modeling and discontinuity mapping in the United States. Bull. Eng. Geol. Environ. 2008, 67, 457–469. [Google Scholar] [CrossRef]
  18. Ferrero, A.M.; Forlani, G.; Roncella, R.; Voyat, H.I. Advanced geostructural survey methods applied to rock mass characterization. Rock Mech. Rock Eng. 2009, 42, 631–665. [Google Scholar] [CrossRef]
  19. Sturzenegger, M.; Stead, D. Close-range terrestrial digital photogrammetry and terrestrial laser scanning for discontinuity characterization on rock cuts. Eng. Geol. 2009, 106, 163–182. [Google Scholar] [CrossRef]
  20. Gigli, G.; Casagli, N. Semi-automatic extraction of rock mass structural data from high resolution LIDAR point clouds. Int. J. Rock Mech. Min. Sci. 2011, 48, 187–198. [Google Scholar] [CrossRef]
  21. Lato, M.J.; Vöge, M. Automated mapping of rock discontinuities in 3D lidar and photogrammetry models. Int. J. Rock Mech. Min. Sci. 2012, 54, 150–158. [Google Scholar] [CrossRef]
  22. Vöge, M.; Lato, M.J.; Diederichs, M.S. Automated rockmass discontinuity mapping from 3-dimensional surface data. Eng. Geol. 2013, 164, 155–162. [Google Scholar] [CrossRef]
  23. Riquelme, A.J.; Abellán, A.; Tomás, R.; Jaboyedoff, M. A new approach for semi-automatic rock mass joints recognition from 3D point clouds. Comput. Geosci. 2014, 68, 38–52. [Google Scholar] [CrossRef] [Green Version]
  24. Gomes, R.K.; De Oliveira, L.P.L.; Gonzaga, L.; Tognoli, F.M.W.; Veronez, M.R.; De Souza, M.K. An algorithm for automatic detection and orientation estimation of planar structures in LiDAR-scanned outcrops. Comput. Geosci. 2016, 90, 170–178. [Google Scholar] [CrossRef]
  25. Guo, J.; Liu, S.; Zhang, P.; Wu, L.; Zhou, W.; Yu, Y. Towards semi-automatic rock mass discontinuity orientation and set analysis from 3D point clouds. Comput. Geosci. 2017, 103, 164–172. [Google Scholar] [CrossRef]
  26. Caselle, C.; Umili, G.; Bonetto, S.; Ferrero, A.M. Application of DIC analysis method to the study of failure initiation in gypsum rocks. Géotech. Lett. 2019, 9, 35–45. [Google Scholar] [CrossRef]
  27. Sturzenegger, M.; Stead, D.; Elmo, D. Terrestrial remote sensing-based estimation of mean trace length, trace intensity and block size/shape. Eng. Geol. 2011, 119, 96–111. [Google Scholar] [CrossRef]
  28. Umili, G.; Ferrero, A.; Einstein, H.H. A new method for automatic discontinuity traces sampling on rock mass 3D model. Comput. Geosci. 2013, 51, 182–192. [Google Scholar] [CrossRef]
  29. Li, X.; Chen, J.; Zhu, H. A new method for automated discontinuity trace mapping on rock mass 3D surface model. Comput. Geosci. 2016, 89, 118–131. [Google Scholar] [CrossRef]
  30. Riquelme, A.J.; Abellán, A.; Tomás, R. Discontinuity spacing analysis in rock masses using 3D point clouds. Eng. Geol. 2015, 195, 185–195. [Google Scholar] [CrossRef] [Green Version]
  31. Buyer, A.; Schubert, W. Calculation the spacing of discontinuities from 3D point clouds. Procedia Eng. 2017, 191, 270–278. [Google Scholar] [CrossRef]
  32. Buyer, A.; Aichinger, S.; Schubert, W. Applying photogrammetry and semi-automated joint mapping for rock mass characterization. Eng. Geol. 2020, 264, 105332. [Google Scholar] [CrossRef]
  33. Eurocode 7: Geotechnical Design—Part 1: General Rules; EN 1997-1; European Committee for standardization: Brussels, Belgium, 2004.
  34. Umili, G.; Bonetto, S.; Ferrero, A.M. An integrated multiscale approach for characterization of rock masses subjected to tunnel excavation. J. Rock Mech. Geotech. Eng. 2018, 10, 513–522. [Google Scholar] [CrossRef]
  35. Ferrero, A.M.; Umili, G.; Migliazza, M.R. Some open issues on the design of protection barriers against rockfall. In Proceedings of the 49th US Rock Mechanics/Geomechanics Symposium 2015, San Francisco, CA, USA, 1–28 July 2015; pp. 3104–3109. [Google Scholar]
  36. UNI 11211–4. Opere di Difesa Dalla Caduta Massi—Parte 4: Progetto Definitivo ed Esecutivo; UNI: Rome, Italy, 2018. [Google Scholar]
  37. Technical Protection Against Rockfall—Terms and Definitions, Effects of Actions, Design, Monitoring and Maintenance; ONR 24810; Austrian Standards Institute: Vienna, Austria, 2017.
  38. Eurocode—Basis of Structural Design; EN 1990:2002+A1; European Committee for standardization: Brussels, Belgium, 2005.
  39. Vagnon, F.; Ferrero, A.M.; Umili, G.; Segalini, A. A factor strength approach for the design of rock fall and debris flow barriers. Geotech. Geol. Eng. 2017, 35, 2663–2675. [Google Scholar] [CrossRef]
  40. De Biagi, V.; Lia Napoli, M.; Barbero, M.; Peila, D. Estimation of the return period of rockfall blocks according to their size. Nat. Hazards Earth Syst. Sci. 2017, 17, 103–113. [Google Scholar] [CrossRef] [Green Version]
  41. Mavrouli, O.; Corominas, J.; Jaboyedoff, M. Size distribution for potentially unstable rock masses and In Situ rock blocks using LIDAR-generated digital elevation models. Rock Mech. Rock Eng. 2015, 48, 1589–1604. [Google Scholar] [CrossRef] [Green Version]
  42. Dal Piaz, G.V.; Bistacchi, A.; Massironi, M. Geological outline of the Alps. Episodes 2003, 26, 175–180. [Google Scholar] [CrossRef] [Green Version]
  43. Dal Piaz, G.V. The Italian Alps: A journey across two centuries of Alpine geology. J. Virtual Explor. 2010, 36, 8. [Google Scholar] [CrossRef]
  44. Bertotti, G.; Picotti, V.; Bernoulli, D.; Castellarin, A. From rifting to drifting: Tectonic evolution of the South-Alpine upper crust from the Triassic to the Early Cretaceous. Sediment. Geol. 1993, 86, 53–57. [Google Scholar] [CrossRef]
  45. Castellarin, A.; Vai, G.B.; Cantelli, L. The Alpine evolution of the Southern Alps around the Giudicarie faults: A Late Cretaceous to Early Eocene transfer zone. Tectonophysics 2006, 414, 203–223. [Google Scholar] [CrossRef]
  46. Regione Lombardia. Carta Geologica 1:250,000. Available online: http://www.w3.org/1999/xlink" xlink:href="https://www.cartografia.servizirl.it/viewer32/index.jsp?parameters={‘srsWkid’:32632,‘serviceLMOperator’:‘include’,‘widgetVisible’:‘Gestisci%20contenuto’,‘servicesLM’:[{‘wkid’:32632,‘queryAndZoom’:null,‘servicename’:’’,‘servicehost’:’’,‘type’:‘ESRI:AGSD’,‘label’:‘Carta%20Geologica%20250.000’,‘layerDefinitions’:[],‘visible’:‘true’,‘url’:‘http://www.cartografia.servizirl.it/expo/rest/services/gpt/cartageo_250/MapServer’,‘docuuid’:‘{018208BD-AD82-4D2A-B195-548D6F3432B4}’,‘layerId’:0,‘alpha’:0.7}]} (accessed on 6 May 2020).
  47. Servizio Geologico d’Italia, 1948–Foglio 35 Riva del Garda della Carta Geologica d’Italia, Scala 1:100,000. Available online: http://193.206.192.231/carta_geologica_italia/tavoletta.php?foglio=35 (accessed on 6 May 2020).
  48. Regione Lombardia. Gargnano, Piano di Governo del Territorio, Carta Geologica Est. Available online: https://www.multiplan.servizirl.it/pgtweb/pub/pgtweb.jsp (accessed on 6 May 2020).
  49. Bigi, G.; Cosentino, D.; Parotto, M.; Sartori, R.; Scandone, P. Scientific Coordination and Editing; Structural model of Italy, scale 1:500,000; CNR: Trento, Italy, 1983. [Google Scholar]
  50. Miles, R.E. The random division of space. Adv. Appl. Probab. 1972, 4, 243–266. [Google Scholar] [CrossRef]
  51. Palmstrøm, A. Characterizing rock masses by the RMi for use in practical rock engineering. Tunn. Undergr. Sp. Technol. 1996, 11, 175–188. [Google Scholar] [CrossRef]
  52. Lu, P. The Characterisation and Analysis of In-Situ and Blasted Block Size Distributions and the Blastability of Rock Masses. Ph.D. Thesis, Queen Mary and Westfield College, University of London, London, UK, 1997. [Google Scholar]
  53. Ferrero, A.M.; Umili, G. Comparison of methods for estimating fracture size and intensity applied to Aiguille Marbrée (Mont Blanc). Int. J. Rock Mech. Min. Sci. 2011, 48, 1262–1270. [Google Scholar] [CrossRef]
  54. Wang, H. Predictions of In-Situ and Blastpile Block Size Distributions of Rock Masses, with Special Reference to Coastal Requirements. Ph.D. Thesis, Queen Mary and Westfield College, London University, London, UK, 1992. [Google Scholar]
  55. International Society for Rock Mechanics. ISRM commission on standardization of laboratory and field tests: Suggested methods for the quantitative description of discontinuities in rock masses. Int. J. Rock. Mech. Min. 1978, 15, 319–368. [Google Scholar] [CrossRef]
  56. Ruiz-Carulla, R.; Corominas, J.; Mavrouli, O. A methodology to obtain the block size distribution of fragmental rockfall deposits. Landslides 2015, 12, 815–825. [Google Scholar] [CrossRef] [Green Version]
  57. Ruiz-Carulla, R.; Corominas, J.; Mavrouli, O. A fractal fragmentation model for rockfalls. Landslides 2017, 14, 875–889. [Google Scholar] [CrossRef]
  58. Marchelli, M.; De Biagi, V. Optimization methods for the evaluation of the parameters of a rockfall fractal fragmentation model. Landslides 2019, 16, 1385–1396. [Google Scholar] [CrossRef]
Figure 1. (A) Tectonic map of the central part of the Alps (modified after [49]). North-vergent collisional belt (North of Periadriatic Lineament): Austroalpine (AU), Helvetic–Dauphinois (HE), Penninic (PE). Southern Alps (SA) rest South of Periadriatic Lineament. Periadriatic lineament (PL) and fronts of buried thrusts in the Po Plain are shown. (B) Structural setting of the Lake Garda region: (TT) Tremosine–Tignale and (BG) Ballino–Garda fault systems (modified after [49]). The yellow dots in (A,B) indicate the study area.
Figure 1. (A) Tectonic map of the central part of the Alps (modified after [49]). North-vergent collisional belt (North of Periadriatic Lineament): Austroalpine (AU), Helvetic–Dauphinois (HE), Penninic (PE). Southern Alps (SA) rest South of Periadriatic Lineament. Periadriatic lineament (PL) and fronts of buried thrusts in the Po Plain are shown. (B) Structural setting of the Lake Garda region: (TT) Tremosine–Tignale and (BG) Ballino–Garda fault systems (modified after [49]). The yellow dots in (A,B) indicate the study area.
Geosciences 10 00223 g001
Figure 2. (A) Folds deforming the Jurassic succession. (B,C) Details of the marly-limestone multilayer in the Medolo Formation.
Figure 2. (A) Folds deforming the Jurassic succession. (B,C) Details of the marly-limestone multilayer in the Medolo Formation.
Geosciences 10 00223 g002
Figure 3. Plot reporting the main discontinuity sets obtained by processing scanlines data.
Figure 3. Plot reporting the main discontinuity sets obtained by processing scanlines data.
Geosciences 10 00223 g003
Figure 4. Poles plot of two different surveys along scanlines.
Figure 4. Poles plot of two different surveys along scanlines.
Geosciences 10 00223 g004
Figure 5. Rock face and survey windows (red polygons) analyzed in detail with non-contact surveys.
Figure 5. Rock face and survey windows (red polygons) analyzed in detail with non-contact surveys.
Geosciences 10 00223 g005
Figure 6. Example of a plot of a single window and the relative sets.
Figure 6. Example of a plot of a single window and the relative sets.
Geosciences 10 00223 g006
Figure 7. Spacing cumulative distribution function (CDF) obtained for each of the main discontinuity sets.
Figure 7. Spacing cumulative distribution function (CDF) obtained for each of the main discontinuity sets.
Geosciences 10 00223 g007
Figure 8. In situ block size distribution (IBSD) calculated according to Equation (2).
Figure 8. In situ block size distribution (IBSD) calculated according to Equation (2).
Geosciences 10 00223 g008
Figure 9. Sample CDF (SCDF) of the considered terns.
Figure 9. Sample CDF (SCDF) of the considered terns.
Geosciences 10 00223 g009
Figure 10. Monte Carlo CDF (MCDF) of the considered terns.
Figure 10. Monte Carlo CDF (MCDF) of the considered terns.
Geosciences 10 00223 g010
Figure 11. Comparison of block volume distributions obtained for term T123.
Figure 11. Comparison of block volume distributions obtained for term T123.
Geosciences 10 00223 g011
Figure 12. Comparison of block volume distributions obtained for term T124.
Figure 12. Comparison of block volume distributions obtained for term T124.
Geosciences 10 00223 g012
Figure 13. Comparison of block volume distributions obtained for term T125.
Figure 13. Comparison of block volume distributions obtained for term T125.
Geosciences 10 00223 g013
Figure 14. Comparison of block volume distributions obtained for term T135.
Figure 14. Comparison of block volume distributions obtained for term T135.
Geosciences 10 00223 g014
Figure 15. (A) One of the images processed with ImageJ for surveying blocks at the foot of the slope. (B) Frequency of detached blocks (SV), considering the two assumed shapes.
Figure 15. (A) One of the images processed with ImageJ for surveying blocks at the foot of the slope. (B) Frequency of detached blocks (SV), considering the two assumed shapes.
Geosciences 10 00223 g015
Figure 16. (A) Marly layers and bedding planes in limestones often not detectable by non-contact methods. (B) High-angle discontinuities fracturing bedded limestones and contributing to the small size of blocks observed onsite. (C) A big volume of rock mass isolated by the intersection of the most spaced discontinuity sets with the bedding.
Figure 16. (A) Marly layers and bedding planes in limestones often not detectable by non-contact methods. (B) High-angle discontinuities fracturing bedded limestones and contributing to the small size of blocks observed onsite. (C) A big volume of rock mass isolated by the intersection of the most spaced discontinuity sets with the bedding.
Geosciences 10 00223 g016
Table 1. Coefficients Ci,p with 90% confidence intervals for the relationship in Equation (2), considering different distributions of the product of principal mean spacing values of discontinuities (modified from [5]).
Table 1. Coefficients Ci,p with 90% confidence intervals for the relationship in Equation (2), considering different distributions of the product of principal mean spacing values of discontinuities (modified from [5]).
PassingUniformNegative Exp.Log-NormalFractal
%Ci,pRangeCi,pRangeCi,pRangeCi,pRange
100.3750.1570.3320.1310.4690.0990.46490.0128
200.7000.2920.7100.2490.9650.2071.16850.0409
301.0520.4351.2070.4231.5130.3342.16060.0748
401.4600.6071.8520.6452.2200.5423.54580.1612
501.9390.7872.7080.9843.0990.7315.31650.192
602.5481.0363.9801.5504.2871.0298.09030.3098
703.3431.3845.8672.5965.9561.50113.39200.8398
804.4951.8028.9484.5818.4972.24322.60702.2649
906.6232.69115.3329.53213.3774.22739.66605.0295
10017.7729.34838.99223.73438.27717.569108.97009.6708
Table 2. Average orientation and spacing range of the main discontinuities sets obtained from traditional surveys.
Table 2. Average orientation and spacing range of the main discontinuities sets obtained from traditional surveys.
SetDipDip DirectionSpacing Range (m)
K1503440.02–1.55
K2741580.22–3.38
K3a600900.02–2.88
K3b802760.07–3.44
K4640500.18–1.86
K5592140.15–0.96
Table 3. Number of discontinuities surveyed in each window.
Table 3. Number of discontinuities surveyed in each window.
WindowNumber of Discontinuities
3246
3A282
3B192
3C176
3D210
3E146
3F184
3G144
TOTAL1580
Table 4. Average orientation and spacing of the main discontinuity sets obtained from non-contact surveys.
Table 4. Average orientation and spacing of the main discontinuity sets obtained from non-contact surveys.
SetDipDip DirectionSpacing Range (m)
K1503100.07–2.87
K270170-
K3a700800.08–1.60
K3b80260
K4600300.30–19.21
K5402402.72–9.43
Table 5. The goodness of fit evaluated through the Kolmogorov–Smirnov (K–S) test (significance level α = 0.01). Only accepted distributions are reported.
Table 5. The goodness of fit evaluated through the Kolmogorov–Smirnov (K–S) test (significance level α = 0.01). Only accepted distributions are reported.
SetNumber of DataDistributionp-ValueBest Performance
K1266Gamma0.061Log-normal
Log-normal0.148
K212Gamma0.480Log-normal
exponential0.349
Log-normal0.718
Weibull0.488
K3241Log-normal0.263Log-normal
K4121Weibull0.024Log-normal
Log-normal0.078
K526Gamma0.166Gamma
exponential0.010
Log-normal0.099
Weibull0.128
Table 6. Sample arithmetic mean (SM) and distribution mean (DM) of the total spacing sample of each set.
Table 6. Sample arithmetic mean (SM) and distribution mean (DM) of the total spacing sample of each set.
SetSM [m]DM [m]
K10.2830.293
K20.4630.472
K30.4290.434
K42.0202.183
K53.2493.611
Table 7. Comparison of ratios between volumes calculated considering different mean spacing types.
Table 7. Comparison of ratios between volumes calculated considering different mean spacing types.
TernVsm [m3]Ratio Vsm/VMCDF(50%)Vdm [m3]Ratio Vdm/VMCDF(50%)
T1230.282.170.292.32
T1240.783.170.883.57
T1251.212.911.413.40
T1351.352.921.553.37
Table 8. Ratios calculated for the considered terms.
Table 8. Ratios calculated for the considered terms.
TernRatio VIBSD(100%)/SVmax [-]Ratio VMCDF(95%)/SVmax [-]
ParallelepipedEllipsoidParallelepipedEllipsoid
T12317.711.233.521.1
T124243.3153.4134.484.7
T125549.5346.4223.4140.8
T135204.5128.9242.9153.1

Share and Cite

MDPI and ACS Style

Umili, G.; Bonetto, S.M.R.; Mosca, P.; Vagnon, F.; Ferrero, A.M. In Situ Block Size Distribution Aimed at the Choice of the Design Block for Rockfall Barriers Design: A Case Study along Gardesana Road. Geosciences 2020, 10, 223. https://doi.org/10.3390/geosciences10060223

AMA Style

Umili G, Bonetto SMR, Mosca P, Vagnon F, Ferrero AM. In Situ Block Size Distribution Aimed at the Choice of the Design Block for Rockfall Barriers Design: A Case Study along Gardesana Road. Geosciences. 2020; 10(6):223. https://doi.org/10.3390/geosciences10060223

Chicago/Turabian Style

Umili, Gessica, Sabrina Maria Rita Bonetto, Pietro Mosca, Federico Vagnon, and Anna Maria Ferrero. 2020. "In Situ Block Size Distribution Aimed at the Choice of the Design Block for Rockfall Barriers Design: A Case Study along Gardesana Road" Geosciences 10, no. 6: 223. https://doi.org/10.3390/geosciences10060223

APA Style

Umili, G., Bonetto, S. M. R., Mosca, P., Vagnon, F., & Ferrero, A. M. (2020). In Situ Block Size Distribution Aimed at the Choice of the Design Block for Rockfall Barriers Design: A Case Study along Gardesana Road. Geosciences, 10(6), 223. https://doi.org/10.3390/geosciences10060223

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop