Photogrammetric Prediction of Rock Fracture Properties and Validation with Metric Shear Tests

: An accurate understanding of jointed rock mass behavior is important in many applications ranging from deep geological disposal of nuclear waste, to deep mining, and to urban geoengineering projects. The roughness of rock fractures and the matching of the fracture surfaces are the key contributors to the shear strength of rock fractures. In this research, push shear tests with three normal stress levels of 3.6, 6.0, and 8.5 kPa were conducted on two granite samples with artiﬁcially induced well-matching tensile fractures with sizes of 500 mm × 250 mm and 1000 mm × 500 mm. The large sample reached on average a − 60% weaker peak shear stress than the medium-sized sample, and a strong negative scale effect was observed in the peak shear strength. The roughness of the surfaces was measured using a proﬁlometer and photogrammetry. The scale-corrected proﬁlometer-based method (joint roughness coefﬁcient, JRC) underestimates the peak friction angle for the medium-sized slabs by − 27% for the medium sample and − 9% for the large sample. The photogrammetry-based ( Z (cid:48) 2 ) method produces an estimate with − 7% (medium) and + 12% (large) errors. The photogrammetry-based Z (cid:48) 2 is an objective method that consistently produces usable estimates for the JRC and peak friction angle. Author Contributions: Conceptualization, L.U.; methodology, E.C.H., L.U., M.T.; literature M.T.; software, L.U., M.T., E.C.H.; validation, L.U. and M.T.; formal analysis, M.T.; investigation, E.C.H.; resources, L.U.; data curation, E.C.H.; E.C.H., writing—review


Introduction
The roughness of rock joints and the matedness of the joint surfaces are key contributors to the shear strength of rock joints. For a jointed rock mass, in addition to the intact material strength, the strength of the rock mass is influenced by the orientation, geometry, roughness, and matedness of the rock joints. An accurate understanding of jointed rock mass behavior is important in many applications ranging from deep geological disposal of nuclear waste to deep mining to geoengineering projects below urban environments.
The scale effect is significant, as it affects shear tests and the measurement resolution of joint roughness coefficient (JRC). The JRC is used to show the degree of roughness [5,13]. In most cases, there seems to be a lack of unanimous agreement on the effects of the sample size on shear strength. As Bandis et al. [3] affirm, the majority of the cases fail to yield a concrete conclusion as to whether this effect on shear strength is negative, positive, or neutral.

Overview of the Research Methods Used
This work was conducted as follows (1) Preparation of rock samples: One large-scale slab (2.000 m × 1.000 m × 0.130 m) and one medium-scale slab (0.500 m × 0.250 m × 0.045 m), with horizontal artificial fractures. (2) Evaluation of fracture surface roughness using a profilometer and Barton and Choubey profiles. (3) Preparation of 3D photogrammetry-based models of the rock samples. (4) Calculation of JRC values by statistical methods from the 3D models. (5) Shear testing of the rock samples. (6) Evaluation and comparison of the results. (7) Comparison of the results from the tests.

Manufacturing of the Slab Pairs
The rock samples originate from intact Finnish Kuru grey granite blocks [63]. The rock fracture was induced by drilling short holes in the perimeter of the rock block and by splitting the slab by applying tensile stress using plugs, the so-called plug and feather technique. Defects of the rock block were detected based on the mine crew's experience to select the best direction for drilling. Holes placed 15 cm away from each other were drilled along a line representing a projected crack. Afterward, 10 cm plug and feather sets were embedded in the holes in such a way that the direction of the feathers' ears was perpendicular to the projected cracks. Subsequently, the plugs were gradually hammered into the holes with a rubber hammer. After observing an initial crack, the fracture direction indicates which plug should be hammered to next. After the crack has propagated throughout the whole rock block, the two parts of the rocks were split with a wrench. After splitting, the rock blocks were cut to the desired dimensions (2000 mm × 1000 mm × 130 mm and 500 mm × 250 mm × 45 mm).

Surface Roughness Measurements
The roughness of the fracture surfaces was measured in two ways. In the contact method, the JRC values were approximated by comparing the profiles defined by the profilometer with Barton and Choubey's standard profiles [5]. The second approach utilized the profiles obtained through digitalization of the fracture surfaces. The JRC values (JRC i ) obtained by both methods for each section (i) were averaged to calculate the final JRC value for each fracture surface (Equation (1))

Roughness Measurements Using the Profilometer
As shown in Figure 1, three lines in the planned shearing direction along Y-axis were drawn to measure the roughness of the slabs' surfaces. For the large sample, the lines are spaced at 250 mm intervals along the X-axis, while for the medium sample at 62.5 mm intervals. The photographed profilometer profiles were assigned every 10 cm, the JRC value was assessed and compared with the Barton and Choubey standard profiles [5].
Mapping scratched points from shear tests in different stages was accomplished through visual observation and with a ruler. The criterion used to define damage was shear or scratch marks on the surface, such as dust or debris of granite.

Roughness Measurements Using Photogrammetry
The rock fracture surfaces were scanned by taking a number of high-resolution photographs with two dip angles of 30° and 60° to reconstruct a 3D model of the fracture surfaces using the structure-from-motion photogrammetric method (Figure 2a). The camera was mounted on a camera slider, and the sample was placed in a fixed position ( Figure  2b). A photographing procedure was developed to rationalize and simplify the process. This procedure consists of scanning the entire slab via the combination of several highresolution images taken from two different tilt angles and in the positive/negative directions along X and Y axes. Each point on the slab is seen through this procedure at least eight times, in four different directions (positive and negative of X and Y directions) and from two different tilt angles (30 and 60 degrees). Two selected tilt angles are optimal to maximize point matching accuracy (Figure 2) [64].

Roughness Measurements Using Photogrammetry
The rock fracture surfaces were scanned by taking a number of high-resolution photographs with two dip angles of 30 • and 60 • to reconstruct a 3D model of the fracture surfaces using the structure-from-motion photogrammetric method (Figure 2a). The camera was mounted on a camera slider, and the sample was placed in a fixed position (Figure 2b). A photographing procedure was developed to rationalize and simplify the process. This procedure consists of scanning the entire slab via the combination of several high-resolution images taken from two different tilt angles and in the positive/negative directions along X and Y axes. Each point on the slab is seen through this procedure at least eight times, in four different directions (positive and negative of X and Y directions) and from two different tilt angles (30 and 60 degrees). Two selected tilt angles are optimal to maximize point matching accuracy (Figure 2) [64].
The photogrammetry of different parts of the slab was conducted with different procedures. The origin corner was locked with 13 images taken manually in a circular motion with the camera pointing towards the corners (three-quarters of the corner's periphery) (Figure 3a). The blue arrows in Figure 3a illustrate the position and orientation of the camera at each of the 13 steps, as well as the direction followed. Each corner was photographed in the same way while moving along the edges. The corners were shot only once per angle for a total of 144 images (72 per angle). The sides of the slab were scanned by continuous image-taking along the sides. The camera lens was pointed inwards at a 45 • angle from the edge being photographed. The path followed is represented in Figure 3a by the black dashed arrows and the numbers represent the sequence followed. The arrows show the position and orientation of the camera. Figure 3b, on the other hand, illustrates the backward direction (the purple dashed arrows) and the new position of the camera. After the perimeter had been photographed, the inner area of the slab was shot following similar steps as in the shooting of the edges. Defined paths were used to cover the central part of the slab while providing enough reference points by means of the sequence of the shooting, performing the same backward-and-forwards movement. Parallel to each axis, the lanes were spaced at every 25 cm for the large sample and at every 7.5 cm for the medium sample (Figure 3c,d).
2b). A photographing procedure was developed to rationalize and simplify the process. This procedure consists of scanning the entire slab via the combination of several highresolution images taken from two different tilt angles and in the positive/negative directions along X and Y axes. Each point on the slab is seen through this procedure at least eight times, in four different directions (positive and negative of X and Y directions) and from two different tilt angles (30 and 60 degrees). Two selected tilt angles are optimal to maximize point matching accuracy (Figure 2) [64]. The photogrammetry of different parts of the slab was conducted with different procedures. The origin corner was locked with 13 images taken manually in a circular motion with the camera pointing towards the corners (three-quarters of the corner's periphery) (Figure 3a). The blue arrows in Figure 3a illustrate the position and orientation of the camera at each of the 13 steps, as well as the direction followed. Each corner was photographed in the same way while moving along the edges. The corners were shot only once per angle for a total of 144 images (72 per angle). The sides of the slab were scanned by continuous image-taking along the sides. The camera lens was pointed inwards at a 45° angle from the edge being photographed. The path followed is represented in Figure 3a by the black dashed arrows and the numbers represent the sequence followed. The arrows show the position and orientation of the camera. Figure 3b, on the other hand, illustrates the backward direction (the purple dashed arrows) and the new position of the camera. After the perimeter had been photographed, the inner area of the slab was shot following similar steps as in the shooting of the edges. Defined paths were used to cover the central part of the slab while providing enough reference points by means of the sequence of the shooting, performing the same backward-and-forwards movement. Parallel to each axis, the lanes were spaced at every 25 cm for the large sample and at every 7.5 cm for the medium sample (Figure 3c,d). The photographs were taken using a Canon 5DS R DSLR camera and a Canon 35 mm f/1.4L II USM lens. The ISO sensitivity was set to 400, and the aperture was set to f/11. The methodology used for the acquisition of the pictures follows the basic guidelines developed in the previous stages of the KARMO-Mechanical Properties of Rock Joints project [46,65]. The set of photos obtained from this process was used as the data input to generate point clouds with the VisualSfM open-source software [66]. The point clouds obtained were then denoised and triangulated with the CloudCompare open-source software to create a mesh replicating the rocks' surfaces [67]. JRC values and damaged areas can be determined by the evaluation of the topographies of the fractures.
The comparison between the pre-test and post-test surface point clouds was con- The photographs were taken using a Canon 5DS R DSLR camera and a Canon 35 mm f/1.4L II USM lens. The ISO sensitivity was set to 400, and the aperture was set to f/11. The methodology used for the acquisition of the pictures follows the basic guidelines developed in the previous stages of the KARMO-Mechanical Properties of Rock Joints project [46,65]. The set of photos obtained from this process was used as the data input to generate point clouds with the VisualSfM open-source software [66]. The point clouds obtained were then denoised and triangulated with the CloudCompare open-source software to create a mesh replicating the rocks' surfaces [67]. JRC values and damaged areas can be determined by the evaluation of the topographies of the fractures.
The comparison between the pre-test and post-test surface point clouds was conducted in CloudCompare. It can help identify the areas on the sample surfaces that experienced the most damage during the mated stages of the shear tests. For each slab, pre-test and post-test point clouds were imported into CloudCompare. The pre-test cloud was used as a reference cloud for this step, and the post-test cloud was aligned to it. Next cloud-to-cloud absolute distance was calculated as a scalar field. As it is an absolute distance, the software does not indicate whether the deviation is positive, indicating a gap or negative, indicating an alignment error where the sample occupies the same volume.

Photogrammetric Joint Roughness Coefficient (JRC) Calculation from a 2D Line
To evaluate the joint roughness coefficient (JRC) in the photogrammetry method, statistical solutions were used, giving an estimate of the roughness amplitude and surface texture. To accomplish this objective, 2D cross-sectional lines were extracted by drawing them on the point clouds, similarly to the lines in Figure 1. In this study, root mean square (RMS) characterization of local slopes [18] and modified root mean square [21] were applied to assess the roughness. The slope measurement proposed in [18] is the RMS estimate from the local slopes of the profile with intervals between measured data points. The relationship with JRC and RMS can be represented thus (Equation (2)): where Z 2 represents the RMS (Equation (3)): where z is the height of the profile above the line, N the number of measurements, and ds the incremental width in horizontal direction. The two-dimensional sectioning profiles were normalized with a sampling interval of 0.5 mm to maintain optimal behavior of the statistical JRC modelling, as this interval was originally used for the distance between measurements (ds). The mean inclination angle of fracture asperities was merely taken into account in this method. Neglecting the scale effect and the influence of amplitudes are drawbacks of this method. To overcome these disadvantages, a general logistic function is formulated as Equation (4) [21]: where b represents a constant value when it equals 20 for suggested JRC values and λ signifies a roughness index obtained from Equation (5): where h is the mean height of the surface asperities, defined by Equation (6): where L is the length of the profile, x i and z i are the x-coordinate and the z-coordinate points along the profile, respectively, and Z 2 denotes the Modified Root Mean Square written as Equation (7): The average of the JRC values of the lines was used to determine the JRC value of each surface.

Push-Shear Test Setup
The method selected to obtain the peak and residual shear strengths of the samples was a multi-stage shear test consisting of pushing the top rock slab over the rough surface of its fixed counterpart, according to the recommendations of the ISRM revised method for determining the shear strength of rock fractures in the laboratory [68]. All the stages, their descriptions, and the normal stresses for the large and medium samples are presented in Table 1. Two shearing rates were used, 0.1 mm/min until peak shear strength, followed by a shearing rate of 0.5 mm/min until the moving half had been displaced by 2.5% (50 mm) of the length of the large sample (2000 mm) and 10% (50 mm) of the length of the medium sample (500 mm). Because of the large size of the rock samples, hydraulic cylinders to apply normal stress were ruled out. The displacement of the rock slab would have inclined the cylinders, causing changes in the angle of the applied pressure. Instead, a dead weight pressure was used on top of the slab to maintain the normal component (Table 2). Each bag was weighed and labeled with its actual weight. As a result, the normal stresses used for the shearing test are low compared to in situ rock conditions.  Figure 4 illustrates the test setup before shearing the sample. The top slab rests on the bottom slab with freedom of movement in the horizontal direction, i.e., along the shearing direction, Y-axis. For the large sample, a couple of 2.5 m steel beams are screwed onto the sides of the slab to handle it. Additional shorter steel beams are attached below the steel bars to prevent the rotation of the top slab as a result of the uneven sliding surface. These shorter beams make contact with the bottom slab through four ball-bearing plates to allow horizontal movement and reduce the friction between the elements. U-shaped steel beams with UNP100 profile are used to transmit the load and those below the frame to restrain rotation. The weight of the steel elements is considered as part of the normal stress. In the fifth stage of the shearing test, the frame was removed from the large slab. To have same amounts of normal stress for the large and medium samples, the additional weight for the medium sample was reduced.  The bottom slab rests on a wooden pallet that provides the necessary height to locate the centroid of the top slab at the same elevation as the cylinders. A steel block serves as a barrier to the bottom slab, restricting its movement. The barrier is fixed in place by two heavy-duty floor anchors.
The shearing force is applied by a pair of hydraulic cylinders parallel to the surface of the fracture and the direction of pushing is considered horizontal and transmitted to the top sample through a perpendicular steel beam with an H-shaped cross-section (HEM or H-beam) and a couple of wooden blocks fastened to the steel beams of the frame. A coupling system composed of a pair of steel plates and a ball-bearing plate is used as a connection point between the H-beam and the wooden blocks. The H-beam is restrained in its lateral movement by a pair of ball bearing plates allowing movement along the Zaxis (pushing direction) and Y-axis (vertical direction). Two HSS sections keep the guide cylinders in position. The pushing cylinders connect with the H-beam through a pair of wooden plates to center the cylinders and keep them in place while the force is applied ( Figure 5 and Appendix A).
The displacement in all directions (x-, y-, and z-) is tracked with 12 and 11 linear variable differential transformers (LVDTs) in the large and medium samples, respectively ( Figure 5). Additionally, each piston has a force sensor and an LDVT that measures the longitudinal displacement of the cylinder, and it is also used to measure the shearing dis- The bottom slab rests on a wooden pallet that provides the necessary height to locate the centroid of the top slab at the same elevation as the cylinders. A steel block serves as a barrier to the bottom slab, restricting its movement. The barrier is fixed in place by two heavy-duty floor anchors.
The shearing force is applied by a pair of hydraulic cylinders parallel to the surface of the fracture and the direction of pushing is considered horizontal and transmitted to the top sample through a perpendicular steel beam with an H-shaped cross-section (HEM or H-beam) and a couple of wooden blocks fastened to the steel beams of the frame. A coupling system composed of a pair of steel plates and a ball-bearing plate is used as a connection point between the H-beam and the wooden blocks. The H-beam is restrained in its lateral movement by a pair of ball bearing plates allowing movement along the Z-axis (pushing direction) and Y-axis (vertical direction). Two HSS sections keep the guide cylinders in position. The pushing cylinders connect with the H-beam through a pair of wooden plates to center the cylinders and keep them in place while the force is applied ( Figure 5 and Appendix A).
( Figure 5 and Appendix A).
The displacement in all directions (x-, y-, and z-) is tracked with 12 and 11 linear variable differential transformers (LVDTs) in the large and medium samples, respectively ( Figure 5). Additionally, each piston has a force sensor and an LDVT that measures the longitudinal displacement of the cylinder, and it is also used to measure the shearing displacement rate.

Estimation of Peak Friction Angle
An empirical equation (Equation (8)) was developed by Barton [15] to estimate the peak friction angle (φ ) as follows: where φ is the basic friction angle, is the joint compressive strength, σ is the normal stress on the failure plane, and is the joint roughness coefficient. For the Kuru granite, φ was measured with the tilt table test as 33° [47], for the of 216 MPa was used, and is according to Table 2. The JRC values were estimated by the aforementioned methods: the profilometer and statistical methods. The calculated values of the peak friction angle obtained in this way were compared to the values obtained from the back-calculation of the direct shear test results with Equation (9)   The displacement in all directions (x-, y-, and z-) is tracked with 12 and 11 linear variable differential transformers (LVDTs) in the large and medium samples, respectively ( Figure 5). Additionally, each piston has a force sensor and an LDVT that measures the longitudinal displacement of the cylinder, and it is also used to measure the shearing displacement rate.

Estimation of Peak Friction Angle
An empirical equation (Equation (8)) was developed by Barton [15] to estimate the peak friction angle (ϕ p ) as follows: where ϕ b is the basic friction angle, JCS is the joint compressive strength, σ n is the normal stress on the failure plane, and JRC is the joint roughness coefficient.  (9) developed by Bandis et al. [5]: where τ is the peak shear strength of the rock joints.  The photogrammetry data provide an opportunity to detect and inspect fracture surfaces in detail. The initial and the sheared surfaces are illustrated in Figures 7 and 8. By comparing the topography of a fracture surface before and after the shearing test. As detected by the photogrammetry method, damaged areas are highlighted, and the damaged areas are spread almost all over the fracture surfaces. The photogrammetry data provide an opportunity to detect and inspect fracture surfaces in detail. The initial and the sheared surfaces are illustrated in Figures 7 and 8. By comparing the topography of a fracture surface before and after the shearing test. As detected by the photogrammetry method, damaged areas are highlighted, and the damaged areas are spread almost all over the fracture surfaces.

JRC Estimated by the Profilometer
The visual comparison method was adopted for evaluating the JRC values of the fracture surfaces with the profilometer. Each line in Figure 1 was divided into 10 cm sections and the profilometer was positioned over each section. JRC values were estimated for each section by comparing the measured profiles with the 10 standard profiles. Figure 9 shows the manually defined JRC results for the rock slab surfaces before the shear tests. These values were averaged by Equation (1)   The photogrammetry data provide an opportunity to detect and inspect fracture surfaces in detail. The initial and the sheared surfaces are illustrated in Figures 7 and 8. By comparing the topography of a fracture surface before and after the shearing test. As detected by the photogrammetry method, damaged areas are highlighted, and the damaged areas are spread almost all over the fracture surfaces.     The visual comparison method was adopted for evaluating the JRC values of the fracture surfaces with the profilometer. Each line in Figure 1 was divided into 10 cm sections and the profilometer was positioned over each section. JRC values were estimated for each section by comparing the measured profiles with the 10 standard profiles. Figure  9 shows the manually defined JRC results for the rock slab surfaces before the shear tests. These values were averaged by Equation (1)

Photogrammetric JRC Values
The 2D cross-sectional YZ coordinates of each profile were extracted from the photogrammetric models. The JRC values were estimated for the divided 10 cm profiles, similarly to the profilometer method. The JRC values were then calculated for each section using Equation (2) (Figure 10). The distribution of the JRC values is between 6-8 for the large sample and 12-14 for the medium sample.
To consider the effect of scale on the JRC values, Equations (4) to (7) were employed with the extracted photogrammetric data. The calculated JRC values are shown in Figure 11. A comparison among the JRC, Z 2 , Z 2 methods shown in Table 5 reveals that the modified RMS (Z 2 ) resulted in significantly lower JRC values for the large sample from approximately 10 to approximately 7 (Figures 9a, 10a and 11a). There is no remarkable difference between the JRC from the profilometer (9.3) and Z 2 (8.9) methods for the medium sample. However, the JRC values obtained from Z 2 (10.9), or the medium sample, are greater by almost 2 JRC units or 17% of the JRC value.

Photogrammetric JRC Values
The 2D cross-sectional YZ coordinates of each profile were extracted from the photogrammetric models. The JRC values were estimated for the divided 10 cm profiles, similarly to the profilometer method. The JRC values were then calculated for each section using Equation (2) (Figure 10). The distribution of the JRC values is between 6-8 for the large sample and 12-14 for the medium sample.

11.
A comparison among the JRC, Z′2, Z2 methods shown in Table 5 reveals that the modified RMS (Z′2) resulted in significantly lower JRC values for the large sample from approximately 10 to approximately 7 (Figures 9a−11a). There is no remarkable difference between the JRC from the profilometer (9.3) and Z′2 (8.9) methods for the medium sample. However, the JRC values obtained from Z2 (10.9), or the medium sample, are greater by almost 2 JRC units or 17% of the JRC value.

Shear Test Diagrams
Figures 12-18 depict the results for shearing tests. The same methodology was used in both cases. Figure 12a shows the shear results of the large sample for the five stages of shearing (Tables 1 and 2). Figure 12b depicts the shear strengths before the peaks (Table 3).
(b) Pre-shearing of the medium sample (500 mm × 250 mm) Figure 11. The photogrammetry method and Z′2 to define the JRC roughness of the slabs (a) preshearing of the large sample (2000 mm × 1000 mm), (b) pre-shearing of the medium sample (500 mm × 250 mm).

Shear Test Diagrams
Figures 12-20 depict the results for shearing tests. The same methodology was used in both cases. Figure 12a shows the shear results of the large sample for the five stages of shearing (Tables 1 and 2). Figure 12b depicts the shear strengths before the peaks (Table  3). The shear displacement of the sample and the resultant shear strength are represented in the positive portion of the x-axis. The peak and residual envelopes are represented on the negative portion of the x-axis according to the normal stress applied. A nonlinear relationship between the normal stress and the shear strength is seen for both envelopes. This non-linearity is more prominent in the first part of the curves. Figure 12a  residual envelopes, the results present a similar linear relationship as for the large sample (Figures 14 and 15). The results are consistent between both tests; with increasing normal stress, the resulting shear stress increases. Nonetheless, it is not the case of the shear displacement required to achieve peak shear strength, as shown in Table 3. In both cases, the shear displacement to achieve the peak shear strength increased from stage 1 to stage 2, while it decreased during stage 3. The normal stresses used in both large and medium sample tests were low compared to the expected normal stresses hundreds of meters deep underground. The comparison between the results for both samples is illustrated in Figures 14 and 15. The peak and residual strengths from the large test show a clear linear correlation between the shear strength and the normal stress applied. The coefficient of determination R 2 for the peak strength is 0.95 and for the residual strength 0.98 ( Figure 14).
Furthermore, in the medium test, the linear correlation seen in the large sample is consistent in the resultant residual strength at the different stages with a determination coefficient of 0.99. The peak shear strength for the medium sample does not show such a high degree of determination (0.84). The peak shear strength for the fourth stage of the  As shown in Table 4, the scale effect does not exist in the residual friction angle. In contrast, the peak friction angle for the large test (53°) is reduced to almost 75% of the medium test result (71°) and the peak cohesion of the large test (1.3 kPa) is reduced to 36% of the medium test result (3.5 kPa).  As shown in Table 4, the scale effect does not exist in the residual friction angle. In contrast, the peak friction angle for the large test (53°) is reduced to almost 75% of the medium test result (71°) and the peak cohesion of the large test (1.3 kPa) is reduced to 36% of the medium test result (3.5 kPa). a drop in the height of the curves after each stage, but especially in the gap between the curves of stages 1 and 4 under the same normal stress. In stage 5, the dilation presents a negative value for a large part of the experiment, corresponding to the low or non-existent degree of matedness between the surfaces and the open character of the fracture. This phenomenon was also reported in [58]. They concluded that fractures with large apertures have the lowest shear strength since the matedness of such fractures is poorer, i.e., fewer resistance contact points. In contrast to the normal stress-dilation relationship of the large sample, the medium sample showed an odd response ( Figure 17). The stage 1 normal displacement was partially lost due to incorrect sensor reset and maxing out of S7, S8, S9, S10. At the beginning of the shearing, and consistently for stages 1 to 4, but most noticeably in stages 1 and 2, the dilation recorded for the top slab went up quickly with a sudden drop and fast axial displacement, corresponding to the moment when the peak shear strength was overcome. After this point, the dilation continued with an increasing trend as in the large case. While the dilation decreased with increasing normal stress and shearing stage in the large test, this statement did not become a tendency for the medium sample. Figure 17 shows there was no correlation between the normal stress applied and the dilation. With increasing  asperities are overrun by the shear force, resulting in the vertical displacement of the top slab. From the results observed in both tests, stage 1 presented the largest shear displacement before transition from negative to positive normal displacement of the four stages. This could correspond to the fact that with every further shearing stage, the micro-scale surface asperities become crushed, and the contact surface between the rock samples evens out.

Method Comparison for Peak Friction Angles
The shear strength tests results were compared with the roughness analysis made with the profilometer and the photogrammetry measurements ( Table 5). The JRC values of the direct shear tests were back-calculated with the outcomes from the first stage of the shear tests and Equation (9). The JRC values of the direct shear tests were influenced by the sample size, decreasing from 9.54 for the medium sample to 6.17 for the large-scale one by approximately 35%. The scaled JRC significantly underestimates the JRC by −28%, the scaled Z2 method underestimates the JRC by −19%, and the Z′2 method by −7%. The peak friction angles estimated with Equation (8) were close to the results of the shear test between 75° to 78°, except the peak friction angle predicted by Z2, that results in 70°.
Regarding the large-scale sample, the scaled JRC underestimates the JRC by −9% and the scaled Z2 method gives a result similar to the scaled JRC with −9% underestimation. The modified RMS (Z′2) produces a result (6.95, + 13%) comparable with the experimental result (6.17). The estimated JRC value by the modified RMS (Z′2) method and the JRC from  The shear displacement of the sample and the resultant shear strength are represented in the positive portion of the x-axis. The peak and residual envelopes are represented on the negative portion of the x-axis according to the normal stress applied. A non-linear relationship between the normal stress and the shear strength is seen for both envelopes. This non-linearity is more prominent in the first part of the curves. Figure 12a presents some downward spikes regarding the stick-slip phenomenon after the peak shear strength was achieved; this effect was also previously reported for the pulling test [49,54]. During the loading, the pistons apply incremental force until the friction opposed by the surface's roughness is overcome, resulting in shear displacement and a sudden drop in the force registered by the sensors, then building up until again reaching the limit of frictional resistance. This process was repeated several times until 2.5% of the sample length (5 cm for the large sample and 1.25 cm for the medium sample) was reached. From the results obtained and shown in Figure 12a, increasing normal stress the peak shear strength requires larger elastic deformation and displacement to reach the peak ( Table 3). The shear strength profile for stage 4 presents a lower peak compared to stage 1 with the same normal stress. Additionally, the resulting curve seems more homogeneous than the previous stages, possibly meaning that the bigger asperities on the surface opposing the displacement in the previous stages have been crushed, and the sliding is smoother.
The peak shear strength was achieved before 1 mm of shear displacement for every loading stage. As the normal stress was increased, so was the shear displacement necessary to reach the peak strength. Once the extra pressure was taken away from the samples, the displacement required to reach the peak dropped again, as shown in Table 3.
Similarly, Figure 13a illustrates the results of the shearing test on the medium sample. Figure 13b shows the shear strengths before the peaks (Table 3). In terms of the peak and residual envelopes, the results present a similar linear relationship as for the large sample (Figures 14 and 15). The results are consistent between both tests; with increasing normal stress, the resulting shear stress increases. Nonetheless, it is not the case of the shear displacement required to achieve peak shear strength, as shown in Table 3. In both cases, the shear displacement to achieve the peak shear strength increased from stage 1 to stage 2, while it decreased during stage 3.
The normal stresses used in both large and medium sample tests were low compared to the expected normal stresses hundreds of meters deep underground. The comparison between the results for both samples is illustrated in Figures 14 and 15. The peak and residual strengths from the large test show a clear linear correlation between the shear strength and the normal stress applied. The coefficient of determination R 2 for the peak strength is 0.95 and for the residual strength 0.98 ( Figure 14).
Furthermore, in the medium test, the linear correlation seen in the large sample is consistent in the resultant residual strength at the different stages with a determination coefficient of 0.99. The peak shear strength for the medium sample does not show such a high degree of determination (0.84). The peak shear strength for the fourth stage of the medium test represented 58% of the shear strength achieved during the first stage. In contrast, in the large sample, this same relation was only 25%. The linear Mohr-Coulomb failure criteria were fitted to all the results obtained from the large and medium samples and both peak shear and residual shear strengths ( Figure 14). The approximated characteristics of the mechanical properties of the rock fractures, cohesion, and friction angle are summarized in Table 4. As shown in Table 4, the scale effect does not exist in the residual friction angle. In contrast, the peak friction angle for the large test (53 • ) is reduced to almost 75% of the medium test result (71 • ) and the peak cohesion of the large test (1.3 kPa) is reduced to 36% of the medium test result (3.5 kPa).
The peak strengths and residual strengths are illustrated in Figure 15. The difference in the values of the residual strength is not significant. The peak shear strength values registered for the large sample represent 39%, 45%, 44%, and 50% of the medium sample peak shear strength for stages 1 to 4, respectively (Table 3). On the other hand, the residual strength showed a contrary behavior, with a proportion (large/medium) for residual shear strength of 95%, 100%, 100%, and 114% for stages 1 to 4, respectively (Table 3). With increasing sample size, the shear strength is reduced. Nevertheless, the relationship between the residual shear strength and the sample scale does not follow the same trend.
Suppose the degree of matedness for fresh, unfilled, and rough rock fractures is decreased and the number of contact points reduced. In that case, a decrease in the mean peak shear strength is expected due to a reduction of the larger asperities opposing the shearing [14]. This effect can be seen with the reduced peak shear strength of stages 1 and 4 under the same normal stress. Additionally, the matedness effect can be observed in the dilation curves of the different stages, with a special remark on the fifth stage of both samples. For the fifth stage, the slab was turned through 180 degrees, with the aim being to reach the lowest matedness coefficient, and then pushed along the long side as in the previous stages. As the top slab was rotated 180 degrees, the matedness of the surfaces was practically zero. During the large test, the lateral restraints of the shearing slab (top sample) had to be removed to allow transversal movement or a yaw in-plane rotation effect of the top slab while sliding over its bottom counterpart. The absence of the restraints resulted in a slightly (−0.2 kPa) lower normal stress during stage 5 and some yaw effect during the test.
Nevertheless, both the variation in the normal stress (−0.2 kPa) and the yaw effect were minimal; the top sample remained relatively straight even without the restraints. Finally, during the fifth stage, the slab was pushed until it reached 100 mm (10% of the length of the sample) for a more reliable zero matedness effect on the test. Because of the length of the test, an additional shearing rate to the normal testing rates was used for stage 5; 0.1 mm/min up to a 3.0 mm displacement, 0.5 mm/min until a 50 mm displacement, and 1.0 mm/min for the remaining 50 mm of shearing. Due to lack of matedness, the stick-slip phenomena was less prominent in the 5th stage.
For the fifth stage with the medium sample, the normal stress was set to match the 3.4 kPa used on the large sample. No further changes were made to the shear box setup for this scale, which reduced the probability of the yaw effect as lateral restraints were in place. The average normal displacement of the large sample during the shearing of the top slab was 9.86 mm after the 50 mm movement in the shearing direction. Since the pushing occurred in the uphill direction, the front sensors (S 5 and S 6 in Figure 5b) registered readings 0.7-1.0 mm higher than the rear sensors (S 3 and S 4 in Figure 5b). For the large sample, the elevation difference was 14.9 mm uphill over 2000 mm or 0.43 • , and for the medium sample the elevation difference was 2.3 mm uphill over 500 mm or 0.26 • . Considering that almost every single natural joint set encloses certain inclination degrees of its shearing plane (joint surfaces) regarding the shearing direction [69], a correction of the dilation data for the large test was performed to account for the influence of the slope effect. The same correction was applied in evaluating the dilation results of the medium sample, which also sheared uphill. Figures 16 and 17 present the adjusted curves of dilation phenomena for the large and medium slabs, respectively, once the uphill effect had been corrected.
During the large test, the dilation curve of the rock was as expected. The normal displacement curves kept growing at a reducing pace until almost reaching the point at which the curves start to flatten, at 50 mm horizontal displacement. The observed average dilation decreased with every shearing stage as the asperities opposing the greater resistance to displacement were becoming crushed ( Figure 16). The additional weight placed on top of the slab also contributed to the crushing of the asperities, as observed by a drop in the height of the curves after each stage, but especially in the gap between the curves of stages 1 and 4 under the same normal stress. In stage 5, the dilation presents a negative value for a large part of the experiment, corresponding to the low or non-existent degree of matedness between the surfaces and the open character of the fracture. This phenomenon was also reported in [58]. They concluded that fractures with large apertures have the lowest shear strength since the matedness of such fractures is poorer, i.e., fewer resistance contact points.
In contrast to the normal stress-dilation relationship of the large sample, the medium sample showed an odd response ( Figure 17). The stage 1 normal displacement was partially lost due to incorrect sensor reset and maxing out of S 7 , S 8 , S 9 , S 10 . At the beginning of the shearing, and consistently for stages 1 to 4, but most noticeably in stages 1 and 2, the dilation recorded for the top slab went up quickly with a sudden drop and fast axial displacement, corresponding to the moment when the peak shear strength was overcome. After this point, the dilation continued with an increasing trend as in the large case. While the dilation decreased with increasing normal stress and shearing stage in the large test, this statement did not become a tendency for the medium sample. Figure 17 shows there was no correlation between the normal stress applied and the dilation. With increasing normal stress and shearing stages, the bigger asperities on the fracture surface should be crushed, reducing their height and the displacement-opposing surface area. This is consistent with the quasi-identical curves of stages 2 and 4, where the results are similar to the large case.
In contrast to the trend seen in the fifth stage in the large test, the fifth stage in the medium test shows increasing dilation with the highest degree of matedness. This result can be explained by the differences in the shear apparatus setup. Since the medium slab fitted onto a shear box, the surrounding frame provided the lateral restraint of the top slab, preventing the yaw effect to a maximum degree.
Finally, as the same normal stress was applied on the rock surface with a smaller area opposing the shear displacement, the shear strength resulted in higher values that could have crushed the bigger asperities in the samples' surfaces sooner, resulting in smaller dilation values. The medium test shows about half of the dilation observed in the end phase of the test compared to the large test. In addition, it is noteworthy that every rock fracture is different. Therefore, the surface profile of both samples is not the same and could vary in a high degree between them, i.e., the height, extent, and distribution of the asperities in the surface of the slabs could be quite different, between sample, even if extracted from the same rock domain. Figure 18 illustrates the transition point from negative to positive values in the dilation experienced by the large top slab at each stage first to fourth of the shear test. The transition at each stage occurs before the slab has been displaced by 0.15 mm. As a technical note, the reversed shearing of 2nd and 3rd tests in Figure 18 are most likely some effects of the loading configuration and the sample support system and should be not interpreted as a part of the measured signal but more as loading start artefacts. The top block was well placed and matched over the bottom sample at each stage since the transition from negative to positive dilation took very little shear displacement to occur. The contrary effect was explained in [48], where the pull test resulted in a significant negative-to-positive transition in dilation at greater sheared distances, leading to a special observation about the probability of an existing mismatch at the beginning of the stages because of bad positioning. Therefore, the top part of the sample makes micro-adjustments as the shearing force is increased until the matching position, represented by the negative part of the dilation curve, is found and the surfaces are locked in position. Subsequently, the asperities are overrun by the shear force, resulting in the vertical displacement of the top slab. From the results observed in both tests, stage 1 presented the largest shear displacement before transition from negative to positive normal displacement of the four stages. This could correspond to the fact that with every further shearing stage, the micro-scale surface asperities become crushed, and the contact surface between the rock samples evens out.

Method Comparison for Peak Friction Angles
The shear strength tests results were compared with the roughness analysis made with the profilometer and the photogrammetry measurements ( Table 5). The JRC values of the direct shear tests were back-calculated with the outcomes from the first stage of the shear tests and Equation (9). The JRC values of the direct shear tests were influenced by the sample size, decreasing from 9.54 for the medium sample to 6.17 for the large-scale one by approximately 35%. The scaled JRC significantly underestimates the JRC by −28%, the scaled Z 2 method underestimates the JRC by −19%, and the Z 2 method by −7%. The peak friction angles estimated with Equation (8) were close to the results of the shear test between 75 • to 78 • , except the peak friction angle predicted by Z 2 , that results in 70 • . Regarding the large-scale sample, the scaled JRC underestimates the JRC by −9% and the scaled Z 2 method gives a result similar to the scaled JRC with −9% underestimation. The modified RMS (Z 2 ) produces a result (6.95, +13%) comparable with the experimental result (6.17). The estimated JRC value by the modified RMS (Z 2 ) method and the JRC from the shear test are similar (6.95 and 6.17), and the estimated peak friction angle and the experimental peak friction angles are 66.22 • (+6% error) and 62.51 • , respectively.
The scale correction for JRC was taken into account by Bandis et al. [3], Equation (10): where JRC scaled represents the scaled JRC values with different sizes rather than the 10 cm laboratory scale, JRC 0 signifies the average of JRC values of each 10 cm section, L scaled refers to the length of samples, and L 0 is 10 cm.

Conclusions
In this research, push shear tests were conducted on two granite samples with artificially induced well-matching tensile fractures. The length and width of the fractures of the large sample were 2000 mm × 1000 mm and 500 mm × 250 mm for the medium sample. A push test was conducted at three normal stress levels of 3.6, 6.0, and 8.5 kPa for both large and medium samples. The displacement was reset for each stress level. The larger sample reached on average a −60% weaker peak shear stress than the medium sample with the same normal stress. In other words, a strong negative scale effect was observed in the peak shear strength. No scale effect was observed in the residual shear stress values, and the residual shear stress remained similar for both fracture sizes. In conclusion, neither conceptual model by Barton and Choubey [5] with a −33% reduction of the peak strength or Johansson and Stille [58], with a −0% predicted reduction of the peak strength, predicted the results correctly.
The roughness of the surfaces was measured using a profilometer and photogrammetry. The scale-corrected profilometer-based method underestimates the joint roughness coefficient (JRC) and the peak friction angle for the medium-sized slabs, with a −27% error. The photogrammetry-based Z 2 estimate resulted in an error of −7%. For the large sample, the scale-corrected profilometer-based method leads to an underestimation of both the JRC and peak friction angle, with a −9% error. Here, the photogrammetry-based method Z 2 produces an estimate with +12% overestimation. In conclusion, the photogrammetry-based method Z 2 is an objective method that consistently produces usable estimates for the JRC and peak friction angle.