Ground collapses over abandoned underground mines constitute a new environmental risk in the world. The high risk associated with subsurface voids, together with lack of knowledge of the geometric and geomechanical features of mining areas, makes abandoned underground mines one of the current challenges for countries with a long mining history. In this study, a stability analysis of Montevecchia marl mine is performed in order to validate a general approach that takes into account the poor local information and the variability of the input data. The collapse risk was evaluated through a numerical approach that, starting with some simplifying assumptions, is able to provide an overview of the collapse probability. The final results is an easy-accessible-transparent summary graph that shows the collapse probability. This approach may be useful for public administrators called upon to manage this environmental risk. The approach tries to simplify this complex problem in order to achieve a roughly risk assessment, but, since it relies on just a small amount of information, any final user should be aware that a comprehensive and detailed risk scenario can be generated only through more exhaustive investigations.
Historic underground mining activity has left us with a critical legacy all over the world. The risk associated with abandoned mines requires careful consideration because these dismissed mines can generate several forms of damage to the surrounding environment: collapse of large areas, subsidence, pollution, etc. As reported by several authors [1–5], abandoned mines have serious implications for the land: a number of authors have studied the impact of mines on surface stability [6–17], while others have focused on pollution of the environment [18–24].
Abandoned mines are now considered a new environmental risk  and, as a consequence, in some countries mine closure requires long-term risk assessment, as well as preparation of closure plans aimed at protecting public health and mitigating environmental damage to the surrounding areas [13, 26, 27]. Although a census has identi˝ed more than 2400 abandoned mines in Italy , law does not currently promote any form of long-term management of mine sites, and guidelines for mine rehabilitation have not yet been defined. Abandoned mines therefore represent a continuing threat to public health in Italy because the closure of these sites will increase the risk of collapse or pollution of the surrounding area. As for many other geo-engineering issues, lack of knowledge and uncertainty are crucial points in problem-solving . In abandoned mines, uncertainty can be related to both the intrinsic heterogeneity of the geological features, and the lack of detailed knowledge of voids in terms of tunnel section, depth of the tunnels, the number of exploitation levels, and so on. This problem is still more critical because no funding is generally allocated to investigate these areas, due to the absence of economic interest once the mine site has been abandoned. Additionally, in some cases direct investigation cannot be performed inside mines due to stability problems or to the presence of water.
Because of unavailable or inaccurate information with respect to other natural risks [30–33] assessment of dismissed mines demands for approaches based on simplifying assumptions for the quantification of risk .
In this work, we propose a methodological approach to evaluate the risk of collapse of abandoned mine sites, addressing the issues related to data unavailability and data uncertainty. The final goal is to provide decision makers with useful tools  capable to set action priorities minimizing the use of public funds.
We focus on a case-study of an underground mine in Northern Italy that was abandoned in 1958 after a massive collapse involving all mining levels. The inaccessibility of this site and the di˚culty to obtain information makes this case a prototype test-bed for the definition of a general methodology for collapse risk scenario evaluation at abandoned mine sites. Section 2 brie˛y presents the case study and the proposed methodology. In Section 3 the di˙erent steps of the methodology are detailed explained thanks to their application to the case study. Section 4 discusses the results of the risk analysis through a “summary approach”. Finally, a general conclusion is reported highlighting the importance to properly consider data uncertainty even if only a rough risk analysis may be achieved.
Though the proposed methodology is here applied to a specific site, it could be easily exported to other case-studies.
2 Materials and methods
The case-study is located in Montevecchia, a small historic village standing on the ˛anks of a 500 m-high isolated hill in Northern Italy. From a geological point of view, the area features Flysch facies, called Flysch di Bergamo, which consists primarily of a sequence of shale, marl and sandstone. A huge megabed, the Missaglia Megabed , is present in the ˛yschoid complex (Figure 1a). According to the literature [37, 38], the Megabed is considered a unique deposit produced by exceptional sedimentary processes such as debris ˛ows and turbidity currents. The term Megabed is used for beds with great thickness and di˙erent composition from their host sequences. The Missaglia Megabed is 40 m thick and consists of megabreccia in the lower part (2-3 m of thickness) and a homogeneous segment of calcarenite-marl in the upper part. From a tectonic point of view, folds are the dominant structures in this area, with a sequence of anticline-syncline WNWESE trending folds. Figure 1b shows the geological section of the exploited Megabed, in which, as stated above, the Megabed is composed of a chaotic band, basal breccia, and a more homogenous calcarenite-marl unit.
Montevecchia mine (Figure 2a) was exploited from the late twenties for almost 3 decades to extract marl rock that was used to produce concrete. According to collected documents, the mining site could be roughly divided into three sections: the east section, where six tunnels were excavated, the central section, which had 7 levels and was the most heavily exploited (in terms of tunnel section), and the west section, which had still seven levels but with smaller size (Figure 2b). All levels develop eastwards from mine entrance, with a maximum length of 450 m and maximum depth of 100 m. The uppermost mine level is at an elevation of 430 a.s.l. and approximately 45 m below ground surface, and the levels were excavated with a planned distance (from base to base) of about 15 m. Nevertheless, a number of miners reported that several tunnels were enlarged (up to 20 m) and the distance between levels was reduced, regardless of the mine exploitation plans. These modifications were most likely performed to obtain greater economic advantages, to the detriment of mine safety.
Table 1 lists the most important events according to the collected historical documents. In 1956 the mine experienced first significant stability problems, with partial collapses of the roof of the lowest level, and, subsequently, between the third and the forth levels. Extraction activities continued but were limited to the first 3 shallower levels and just in the middle section of the mine. In 1958 a large sinkhole of approximately 10000 m2 occurred on top of the hill (Figure 3) because of a massive collapse involving all mining levels in the central part of the mine. Mining activity ceased, the site was abandoned, and the road running above the mine along the hill ridge had to be deviated and repaired; a parking lot that still exists was also created in the area previously a˙ected by the sinkhole. A number of investigations were recently performed in order to evaluate the safety conditions of the mine and the possible risk of other events like the 1958 sinkhole . Failures of the areas to the east and to the west of the previous collapse would directly concern human settlements.
|Year||Main facts and events|
|1928||Opening of marl extraction activities.|
|1930||A technical report states that exploitation of Montevecchia mine cannot cause problems to natural springs.|
|1931||Only the first level is exploited with a the tunnel having a section of about 3-5 m.|
|1932-35||Exploitation of the first level is completed and new deeper levels with larger sections are excavated.|
|1936||Exploitation of the third level is completed.|
|1940||A technical report confirms two important springs has vanished most likely because of mine exploitation. The most important spring in the area vanishes at the end of the year.|
|1947||Exploitation of the 4th and 5th levels is completed. The distance between tunnels is about 15 m from base to base.|
|1954||Several tunnels are exploited with increasing section (7-8 m). There are 7 mining levels and overexploitation a˙ected the central section of the mine.|
|1955||Miners report that some tunnels were enlarged (up to 20m) and the distance between levels was reduced, regardless of the mine exploitation plans|
|1956||First significant failures with partial collapses of the roof of the lowest level, and, subsequently, between the third and the forth levels. Activities are limited to extractions at the first 3 levels and just in the middle section of the mine.|
|1958||Extensive collapse of the mine at all levels and end of extraction activities.|
For a better understanding of the risk associated to the presence of abandoned mines, we propose a 3-step approach depicted in the ˛owchart in Figure 4 and described hereafter. The first step of the methodology concerns the characterization of the problem with the aim of defining a conceptual model of the abandoned mine as detailed as possible. The input data to this step include historical documents and miners’ statements, which are mainly useful to understand the geometrical features of the underground voids, along with geomechanical parameters that are often obtained by means of direct investigations in the field.
The second step of the proposed methodology involves the simplification of the problem through the definition of a limited number of reference states or scenarios, according to which the outcomes of the last step of the procedure will be grouped.
The third and final step is risk assessment and it is achieved through a numerical approach. First of all a stability analysis is performed and then, due to the lack of knowledge coupled with the intrinsic uncertainty associated with abandoned mines, a sensitivity and uncertainty analysis are carried out. Starting from the idea of a simple approach, sensitivity and uncertainty analysis are performed around a reference case: the stability analysis reported in the Paragraph 3.3.1. Finally, a summary of all the possible outcomes are discussed in order to define the final collapse probability.
3.1 Conceptual modelling
For our case-study, the analysis of mining claims, plans, cross sections and technical reports was particularly useful to determine the number of tunnels and their sections (see Paragraph 2.1). Instead, geological and geophysical investigations were carried out to characterize the features of the rock mass. Geological mapping was performed outside the mine across two outcrops in the Flysch facies and in the Megabed and results are listed in Table 2. A Schmidt hammer was also used to measure joint compressive strength (JCS) in Missaglia Megabed. Field data were analyzed in order to gather useful information for defi-nition of the geomechanical properties of marl and ˛ysch. Rock mass classification systems have been successfully applied to provide reliable input data for numerical approaches. Rock mass rating (RMR- ), Q-system (Q) and Geological Strength Index (GSI- ) are the main classification systems used in the mining industry. In this work these methods have been applied to both Flysch and Megabed outcrops (Table 3).
|FLYSCH DI BERGAMO|
|Joint set||dip direction (°)||dip (°)||Persistence (%)||Spacing (m)||aperture (m)|
|Joint set||dip direction (°)||dip (°)||Persistence (%)||Spacing (m)||aperture (m)|
|JCS 52 Mpa|
Though use of rock mass classification is a conventional procedure for defining the input data for numerical simulation, the outcome of a geophysical survey can be used to better estimate some geomechanical properties of the rock mass. Seismic is definitely the most interesting method among the geophysical techniques when the objective is assessment of the rock mass’s mechanical properties. Three techniques that appeared to be promising for the case at hand were selected: refraction seismic, multi-channel analysis of surface waves (MASW), and transillumination experiments (a sort of tomography with very limited angular coverage because of the restrictions on entering the mine). Thanks to the geophysical surveys rock layers with di˙erent properties in term of P- and S-waves were identified [39, 43], from which the rock’s elastic parameters can be estimated.
3.2 Collapse options
For a rough risk assessment a simple classification of the numerical outcomes to facilitate the comparison of the results of several simulations is required. Therefore, an easy way to quantify the risk collapse is introduced by the use of some references states . This is one important simplifying assumption in order to permit a final summary risk analysis that takes into account all the performed simulations. The proposed collapse options are: no collapse (stability), minor collapse in two tunnels at most with some movements of blocks (that means partial instability) and extensive collapse that involves all tunnels, like the event of 1958 (instability). The three risk reference levels are then coded as green, yellow and red. An example is reported in Figure 5.
3.3 Numerical modelling
The management of several natural disasters (landslides, ˛oods, hurricanes and so on) is often based upon the hypothesis of di˙erent risk scenarios [45, 46], defined through use of the available information. Despite this, risk scenarios are valuable only if believable and rigorous from a scientific point of view [47–49].
For the case study, the available data do not permit to carry out a detailed analysis. Therefore, among many different approaches, a simplified numerical method fits the aim of the paper because capable of estimating a hint of the collapse risk with the current knowledge of the site. To accomplish this goal, numerical code was used in two steps: 1) the stability analysis; 2) the multiple scenario that takes into account the variability of the case-study considering a range of possible values of the input data in numerical analysis.
For what concerns the numerical code, jointed rock masses are usually simulated by the distinct element method  or by discontinuous deformation analysis  and among these, in the present work, the distinct element method 3DEC (a three-dimensional code developed by Itasca in 2008) was chosen to provide the stability analysis of this abandoned mine site.
3.3.1 Stability analysis
The first step deals with a simple stability analysis. This analysis is then considered as the reference case around which the uncertainty and sensitivity analysis are performed. A stability analysis with the lack of knowledge of the input data is not a straightforward procedure and, often does not lead to a unique outcome. It is therefore necessary to introduce in a clear way the simplifying assumptions made to overcome the lack of information:
the Missaglia Megabed is simulated with a 110 x 160 x 2 m (width x height x thickness) model that permits to simulate the Megabed and surrounding material up to the contact between Marl and Flysch.
a distinction is made between tunnels features and geomechanical dataset because the uncertainty in tunnel features is an intrinsic problem for this specific risk and, moreover, the level of confidence of this information is poor. Therefore the variability of tunnel sections is included in the analysis with the use of four di˙erent tunnel sections. For each section an outcome is produced. The features of the first tunnel arrangement are reported below:
– horseshoe-shaped tunnel with a section of 7 x 8 m on all seven levels.
– first level located at a depth of 45 m.
– 15 m interval between levels from base to base.
In the other arrangements tunnel size is progressively enlarged, and only the first level remains fixed at 7 x 8 m, as reported in all historical documents.
he geomechanical features are defined coupling all available information . When necessary, the geomechanical properties were also integrated with values from the literature [52, 53]. All the properties of Flysch and Marl are reported in Table 4. The uncertainty of these data will be considered only in the next paragraph, with the sensitivity and uncertainty analysis. The vertical to horizontal stress ratio, k, was fixed at 1.5 and the Mohr-Coulomb constitutive criterion was applied to simulate block deformation, whereas the Coulomb slip failure criterion was adopted for discontinuities. GSI values provide a geomechanical characterization of Flysch, which is simulated as an equivalent continuum medium because of its peculiar structure.
external triggering factors are not considered.
|Numerical model geomechanical characterization|
|Rock matrix Young modulus||Modulus ratio||31||GPa|
|Rock matrix Poisson’s ratio||Hoek Brown law||0.28|
|Joint normal sti˙ness||Literature||0.8||GPa/m|
|Joint shear sti˙ness||Literature||0.4||GPa/m|
|Joint set 1||Survey||123°/86°|
|Joint set 2||Survey||205°/80°|
|Joint set 3||Survey||Sub horizontal|
|Flysch di Bergamo (GSI 40)|
|Rock matrix Young modulus||Modulus ratio||2.9||GPa|
|Rock matrix Poisson’s ratio||Hoek Brown law||0.3|
|Cohesion||Hoek Brown law||0.3||MPa|
|Tensile strength||Hoek Brown law||0.8||MPa|
Figure 6 shows the outcomes of this analysis. This result clearly displays how enlarged voids induce a state of instability in the system, especially at deeper levels. The stability analysis highlights some kinematic features of the mass movements: the deepest levels are prone to collapse and probably this mechanism consequently trigger progressive failure of the upper parts until the system finds a new equilibrium.
3.3.2 Multiple risk scenario
When uncertainty is not limited to few properties, a probabilistic analysis that takes into account the extreme variability of input data is recommended . A sensitivity and uncertainty analysis may be used to have an overview of the real risk considering the contribution of uncertainties in risk analysis .
The stability analysis is considered as the reference case around which the uncertainty and sensitivity analysis are carried out. Sensitivity analysis attempts to define the change in final outcomes due to the variability of input data. A full sensitivity analysis is performed by changing the main properties. Each change was done while the rest of dataset remains constant. In this way the most important input data that mainly contribute to the uncertainty of the risk analysis are defined (Figure 7). The properties ranges explored in the sensitivity analysis are reported in Table 5.
|Sensitivity test||Stability Analysis Value||Variation Interval|
|Joints Friction angle||30°||From 20°to 35°|
|Joints Normal sti˙ness||0.8 GPa/m||From 0.60 GPa/m to 1.00 GPa/m|
|Joints Shear sti˙ness||0.4 GPa/m||From 0.30 GPa/m to 0.50 GPa/m|
|Joints Persistence||0.7||From 0.4 to 1.0|
|Joints Cohesion||1 MPa||From 0 MPa to 3 MPa|
The sensitivity analysis can be summarized as follow:
Joint cohesion and joint friction have no in˛uence on tunnel convergence at existing confining stresses.
Joint sti˙ness and joint persistence are significantly linked to tunnel stability.
Considering the sensitivity analysis, five di˙erent sets of geomechanical properties were used in the uncertainty analysis to generate a range of possible risk scenario. As before, the set named as reference case refers to the first hypothesis of stability analysis. The other four sets are defined accordingly with the Table 6.
|Set 1||Set 2||Reference Case||Set 3||Set 4|
|Joints normal sti˙ness||0.5 GPa/m||0.5 GPa/m||0.8 GPa/m||1 GPa/m||1 GPa/m|
|Joints Friction angle||20°||30°||30°||30°||30°|
As before, four representative tunnel dimensions are considered. For the uncertainty analysis a specific attention is paid to the state of “in situ” rock stress, a condition that remains a challenge in many geo-engineering practices. Indeed, the state of rock stress [56, 57] is di˚cult to achieve even when stress measurements are performed. Numerical simulations were thus conducted for di˙erent values of k ; vertical stress was estimated considering the overburden material, while a range of horizontal stresses was taken into consideration.
The graphs in Figure 8 report the outcomes of 80 simulations (i.e., the combination of 5 di˙erent sets of input dataset with 4 tunnel sections and 4 k values); this Figure highlights how tunnel section and k ratio a˙ect the results.
E˚cient management of any kind of environmental risk requires the definition of a “risk scenario” closest to reality that should be used as the basis for decision-making. This is essential, but, for abandoned mines characterized by great uncertainty, it is di˚cult to achieve. The problem of poor knowledge can be solved in two ways: the first is to reduce the lack of knowledge through additional research and data collection (often very di˚cult for logistic and economic reasons), while the second way attempts to conduct a general analysis considering all possible risk scenarios. This second approach fits well the purpose of this work and it is particular suitable for highly unreliable outcomes. Therefore, the solution must be sought in the collection of di˙erent data and in uncertainty analysis  aimed to generate a set of risk scenarios.
On the basis of the available information and on the integration of sensitivity and uncertainty analysis, a set of scenarios was defined. The results of multiple scenario (Figure 8) were summed considering the reference states (Figure 5) in order to generate an overview of the collapse probability, directly related to tunnel sections. This analysis is reported in the summary graph (Figure 9), that shows the probability of each scenario.
More in detail, the final result may be useful for an overview of the actual residual risk for human settlements located in the western and eastern areas of the mine (Figure 1), which have not yet collapsed as well as for land use and management.
Before starting with the evaluation of the residual risk in Montevecchia mine, the 1958 event was analyzed to validate the summary graph (Figure 9). Considering that historical data reported the central part, where the sinkhole occurred, as the more exploited area with the larger tunnel section, the last case of Figure 9 was considered. As expected the 12x10 m tunnel section shows less than 20% of stable condition, confirming the critical situation that led to the collapse of the central area.
After this validation it was possible to evaluate the residual risk for the other two sections (Figure 9). According to most of the historical documents, the eastern part of the mine was not so heavily exploited and its sections should be 10x8 meters, resulting in a situation of relative stability; it has about a 55% chance of being stable, and only a 20% chance of developing significant instability. The western part has a similar condition, and exploitation was performed on only six levels. Thus this area can also be considered almost stable. Therefore new investigations of the eastern and western areas of the mine site are not currently necessary because the situation seems fairly stable, thus the residual risk in west and east areas of Montevecchia mine is not a priority for land management. Such information will help also to determine if more investigations are necessary. The summary graph in a simple way may allow decision makers to minimize the use of public funds, avoiding expensive studies which a priori are often not useful for a real improvement in the reliability of the final risk scenario.
The assumptions introduced in the previous paragraphs must be considered to evaluate the limit of the approach:
only three reference states to classify outcomes are taken into account and all the results are linked only with one of these states. Moreover only four di˙erent tunnel sections are considered. This is a strong simplification but in order to avoid too many numerical simulations a limited number of cases was chosen. Of course this directly a˙ects the results.
a rigorous probabilistic analysis was not performed. For instance a Montecarlo analysis may improve the analysis of this risk but the numerical analysis could require long computational time and it was not the purpose of this work.
none model uncertainty related to the code was considered. It could be interesting to test also other codes.
k value was considered one of the main uncertainties and several simulations were performed for each of this value. Of course a detailed knowledge of this boundary condition may lead to an improvement of the summary graph.
Considering this limitations, the risk scenario should be refined and updated over time as new data become available. Moreover, to really improve the reliability of this approach, a monitoring system capable of detecting detailed information on rock mass behavior may be used. For instance, remote sensing techniques  could be powerful tools for evaluating surface displacements in the area surrounding abandoned mines [62, 63].
Geo-hazards related to the long-term stability of abandoned mines located in urbanized areas are a challenge for all countries with a mining history. When underground mines may a˙ect human settlements, it is very important to check the stability of mines, even at sites no longer in use. Nowadays, the risks related with active mines are well known and thoroughly studied, but one of the major problems related to underground mines concerns abandoned mine sites. It may be di˚cult to obtain all the information required to quantify the risk associated with these voids after mine site closure, especially in the case of old mines without historical documentation. Moreover, neither guidelines nor procedures are currently available for assessing how to approach this problem. In keeping with these environmental problems, long-term analysis of these sites should be planned.
Starting from these considerations, a simple approach is here proposed. The purpose of this approach is to give an overview of the collapse risk of abandoned mine considering how to tackle input data uncertainty. Frequently, in risk analysis a strong interference between uncertainty and decision making is present and therefore it is crucial to generate an easy, accessible and transparent decision procedure capable to solve also complex problem characterized by poor local information. The summary graph fits well the requirements of transparency (the assumptions and the limits are declared and discussed), accessibility (data are reported in the text and, graphs are easier to understand) and ease (the method is based upon an easy and simple approach).
 Smith G.J., Rosenbaum M.S., Recent underground investigations of abandoned chalk mine workings beneath Norwich City, Norfolk. Eng. Geol., 1993, 36, 67-78.10.1016/0013-7952(93)90019-9Search in Google Scholar
 Hawkins A.B., McDonald C., Hazards associated with old stone mine workings: a case history from Combe Down, Bath. In: Skipp B.O. (Ed.) Risk and reliability in ground engineering, Thomas Telford, 1994, 215-226.Search in Google Scholar
 Wojtkowiak F., Didier C., Principles for a safe closure of old mine shafts and adits. In: Proc. 9th ISRM Congress, Balkema, Rotterdam, 1999, 23-30.Search in Google Scholar
 Bell F.G., Bullock S.E.T., Halbich T.F.J., Lindsay P., Environmental impacts associated with an abandoned mine in the Withbank Coalfield. S. Afr. Int. J. Coal. Geol., 2001, 45(2-3), 195-216.10.1016/S0166-5162(00)00033-1Search in Google Scholar
 Fasani Bianchi G., Bozzano F., Cercato M., The underground cavity network of south-eastern Rome (Italy): An evolutionary geological model oriented to hazard assessment. B. Eng. Geol. Environ., 2011, 70 (4), 533-542.10.1007/s10064-011-0360-0Search in Google Scholar
 Betournay M.C., Mitri H.S., Hassani F., Chimneying Disintegration Failure Mechanism of Hard Rock Mines. 987-996 in Proceedings 1st North American Rock Mechanics Symposium, Austin, TX, Balkema A.A. (Ed.), 1994, 987-996.Search in Google Scholar
 Betournay M.C., The Stability of Shallow Stopes of Hard Rock Mines. PhD Thesis, McGill University, 1995.Search in Google Scholar
 Bekendam R.F., Pottgens J.J., Ground movements over the coal mines of southern Limburg, The Netherlands, and their relation to rising mine waters. IAHS Publ., 1995, 234, 3-12.Search in Google Scholar
 Zhou W., The formation of sinkholes in karst mining areas in China and some methods of prevention. Environ. Geol., 1997, 31(1-2), 50-58.10.1007/s002540050163Search in Google Scholar
 Longoni L., Papini M., The hydrogeology problems of disused mines in Olgiate Molgora (LC). In: First International Conference on Underground Spaces, The New Forest, UK, 2008.10.2495/US080051Search in Google Scholar
 Papini M., Longoni L., Dell’Orto K., Emergency guidelines for two abandoned mines on Piani di Resinelli area. In: First International Conference on Underground Spaces, The New Forest, UK, 2008.10.2495/US080111Search in Google Scholar
 Betti D., Buscarnera G., Castellanza R., Nova R., Numerical analysis of the life-time of an abandoned gypsum mine. In: Proc. 12th International Conference IACMAG 2008, Goa (India), 2008, 1210-1218.Search in Google Scholar
 Bétournay M.C., Abandoned metal mine stability risk evaluation. Risk Analysis, 2009, 29(10), 1355-1370.10.1111/j.1539-6924.2009.01267.xSearch in Google Scholar
 Castellanza R., Nova R., Orlandi G., Flooded gypsum mine remedial by chamber filling, J. Geotech. Geoenviron., 2010, ASCE 136(4), 629-639.10.1061/(ASCE)GT.1943-5606.0000249Search in Google Scholar
 Vongpaisal S., Li G., Pakalnis R., Brady T., New development of expert system module for a decision-making on mine stope stability in underground blasthole mining operations. Int. J. Min. Reclamat. Environ., 2011, 25, 41-51.10.1080/17480930.2011.552279Search in Google Scholar
 Suchowerska A., Merifield R., Carter, J., Clausen J., Prediction of underground cavity roof collapse using the Hoek-Brown failure criterion. Comput. Geotech., 2012, 44, 93-103.10.1016/j.compgeo.2012.03.014Search in Google Scholar
 Helm P.R., Davie C.T., Glendinning S., Numerical modeling of shallow abandoned mine working subsidence a˙ecting transport infrastructure. Eng. Geol., 2013, 154, 6-19.10.1016/j.enggeo.2012.12.003Search in Google Scholar
 Vinkler F., Piguet J.P., Numerical analysis of the long-term stability of an abandoned mine: impact of the groundwater level rise. In: Proc. FLAC Symp.on numerical modeling in Geomech., Balkema, Rotterdam, 1999.Search in Google Scholar
 Lee J.Y. et al., Potential groundwater contamination with toxic metals in around an abandoned Zn mine, Korea. Water Air Soil Pollut., 2005, 165, 167-185.10.1007/s11270-005-4637-4Search in Google Scholar
 Gemici Ü., Evaluation of the water quality related to the acid mine drainage of an abandoned mercury mine (Ala³ehir, Turkey). Environ. Monit. Assess., 2008, 147, 93-106.10.1007/s10661-007-0101-9Search in Google Scholar
 Zobrist J. et al., Environmental and socioeconomic assessment of impacts by mining activities-a case study in the Certej River catchment, Western Carpathians, Romania. Environ. Sci. Pollut. Res. Int., 2009, 16, S14-26.10.1007/s11356-008-0068-2Search in Google Scholar
 Pauwels H., Pettenati M., Gree˛e C., The combined e˙ect of abandoned mines and agriculture on groundwater chemistry. J. Contam. Hydrol., 2010 115(1-4), 64-78.10.1016/j.jconhyd.2010.04.003Search in Google Scholar
 Jordan G., Szucs A., Geochemical landscape analysis: Development and application to the risk assessment of acid mine drainage. A Case study in central Sweden Landscape Research, 2011, 36 (2), 231-261.10.1080/01426397.2010.547569Search in Google Scholar
 Romero F.M., Núñez L., Gutiérrez M.E., Armienta M.A., Ceniceros-Gómez A.E., Evaluation of the potential of indigenous calcareous shale for neutralization and removal of arsenic and heavy metals from acid mine drainage in the Taxco mining area, Mexico. Arch. Environ. Contam. Toxicol., 2011, 60 (2), 191-203.10.1007/s00244-010-9544-zSearch in Google Scholar
 Castellanza R., Nova R., Orlandi G., Evaluation and remediation of an abandoned gypsum mine. J. Geotech. Geoenviron., 2009, 136(4), 629-639.10.1061/(ASCE)GT.1943-5606.0000249Search in Google Scholar
 Hutchinson D.J., Phillips C., Cascante G., Risk considerations for crown pillar stability assessment for mine closure planning. Geotech. Geol. Eng, 2002, 20.1 41-64.10.1023/A:1013852722768Search in Google Scholar
 Ruegsegger L.R., Manual for abandoned underground mine inventory and risk assessment. United States Department of Transportation, Federal Highway Administration, Infrastructure Division and the Ohio Department of Transportation, (FHWA IF-99-007), 2009.Search in Google Scholar
 Berry P., Bandini A., Dacquino C., Classificazione dei siti minerari sotterranei abbandonati mediante un indice di rischio statico-strutturale. (Classification of abandoned mine sites through use of a static risk index) In: Recupero e Valorizzazione delle miniere dismesse: lo stato dell’arte in Italia (Recovery and promotion of abandoned mines: the state of the art in Italy), ISPRA, Quaderni, Ambiente e società n. 3, 2011, 31-43 (in Italian).Search in Google Scholar
 Fisher B., Eberhardt E., Assessment of Parameter Uncertainty Associated with Dip Slope Stability Analyses as a Means to Improve Site Investigations, J. Geotech. Geoenviron. Eng., 2012, 138(2), 166-173.10.1061/(ASCE)GT.1943-5606.0000515Search in Google Scholar
 Park I., Lee J., Saro L., Ensemble of ground subsidence hazard maps using fuzzy logic. Open Geosci., 2014, 6(2), 207-218.10.2478/s13533-012-0175-ySearch in Google Scholar
 Apostolidis E., Koukis G., Engineering-geological conditions of the formations in the Western Thessaly basin, Greece. Open Geosci., 2013, 5(3), 407-422.10.2478/s13533-012-0200-1Search in Google Scholar
 Longoni L., Papini M., Arosio D., Zanzi L., Brambilla D., A new geological model for Spriana landslide. Bull. Eng. Geol. Environ., 2014, 73(4), 959-970.10.1007/s10064-014-0610-zSearch in Google Scholar
 Longoni L., Papini M., Brambilla D., Arosio D., Zanzi L., The role of the spatial scale and data accuracy on deep-seated gravitational slope deformation modeling: The Ronco landslide, Italy, Geomorphology, 2016, 253, 74-82.10.1016/j.geomorph.2015.09.030Search in Google Scholar
 Eckert N., Parent E., Faug T., Naaim M., Optimal design under uncertainty of a passive defense structure against snow avalanches: from a general Bayesian framework to a simple analytical model. Nat. Hazards Earth Syst. Sci., 2008, 8(5), 1067-1081.10.5194/nhess-8-1067-2008Search in Google Scholar
 Bathrellos G.D., Gaki-Papanastassiou K., Skilodimou H.D., Papanastassiou D., Chousianitis, K.G., Potential suitability for urban planning and industry development using natural hazard maps and geological-geomorphological parameters. Environ. Earth Sci., 2012, 66(2), 537-548.10.1007/s12665-011-1263-xSearch in Google Scholar
 Bernoulli D., Bischel M., Bolli H.M., Haering M.O., Hochuli P.A., KLleboth P., The Missaglia Megabed, a catastrophic deposit in the Upper Cretaceous Bergamo Flysch, Northern Italy. Eclogae Geol. Helv., 1981, Basel 74(2), 421-442.Search in Google Scholar
 Rosell J., Wieczorek J., Main features of megaturbidites in the Eocene of southern Pyrenees. Ann. Soc. Geol. Pol., 1989, 59, 3-16.Search in Google Scholar
 Bersezio R., Fornaciari M., Gelati R., Napolitano A., Valdisturlo A., The significance of the Upper Cretaceous to Miocene clastic wedges in the deformation history of the Lombardian southern Alps. Géol. Alp., 1993, (Grenoble) 69, 3-20.Search in Google Scholar
 Arosio D., Longoni L., Papini M., Zanzi L., Seismic characterization of an abandoned mine site. Acta Geophys., 2013, 61, 611-623.10.2478/s11600-012-0090-0Search in Google Scholar
 Bieniawski Z., Engineering classification of jointed rock masses. Civil Engineer in South Africa, 1973, 15(12), 335-344.Search in Google Scholar
 Barton N., Lien R., Lunde J., Engineering classification of rock masses for the design of tunnel support. Rock Mechanics, 1974, 6(4), 189-236.10.1007/BF01239496Search in Google Scholar
 Hoek E., Brown E.T., Practical estimates of rock mass strength Int. J. Rock Mech. Min. Sci., 1997, 34(8), 1165-1186.10.1016/S1365-1609(97)80069-XSearch in Google Scholar
 Arosio D., Zanzi L., Longoni L., Papini M., Shallow seismic investigations of an abandoned mine. In: Near Surface Geoscience 2012 - 18th European Meeting of Environmental and Engineering Geophysics, September 3-5, Paris, France, 2012.10.3997/2214-4609.20143316Search in Google Scholar
 Bathrellos G.D., Kalivas D.P., Skilodimou H.D., GIS-based landslide susceptibility mapping models applied to natural and urban planning in Trikala, Central Greece. Estud Geol, 2009, 65(1), 49-65.10.3989/egeol.08642.036Search in Google Scholar
 Groves D.G., Lempert R.J., A new analytic method for finding policy-relevant scenarios. Global Environ. Chang., 2006, 17(1), 73-85.10.1016/j.gloenvcha.2006.11.006Search in Google Scholar
 Ozturk U., Tarakegn Y.A., Longoni L., Brambilla D., Papini M., Jensen J., A simplified early-warning system for imminent landslide prediction based on failure index fragility curves developed through numerical analysis. Geomat. Nat. Haz. Risk, 2015 (in press).10.1080/19475705.2015.1058863Search in Google Scholar
 Ranger N., Niehörster F., Deep uncertainty in long-term hurricane risk: Scenario generation and implications for future climate experiments. Global Environ. Chang., 2012, 22(3), 703-712.10.1016/j.gloenvcha.2012.03.009Search in Google Scholar
 Bathrellos G.D., Gaki-Papanastassiou K., Skilodimou H.D., Skianis G.A., Chousianitis K.G., Assessment of rural community and agricultural development using geomorphological-geological factors and GIS in the Trikala prefecture (Central Greece). Stoch. Environ. Res. Risk. Assess., 2013, 27(2), 573-588.10.1007/s00477-012-0602-0Search in Google Scholar
 Papadopoulou-Vrynioti K., Bathrellos G.D., Skilodimou H.D., Kaviris G., Makropoulos K., Karst collapse susceptibility mapping considering peak ground acceleration in a rapidly growing urban area. Eng. Geol., 2013, 158, 77-88.10.1016/j.enggeo.2013.02.009Search in Google Scholar
 Cundall P.A., A computer model for simulating progressive large-scale movements in blocky rock system. In: Proceedings of the International Symposium Rock Fracture, ISRM. Proc 2, 1971, 129-136.Search in Google Scholar
 Shi G.H., Discontinuous Deformation Analysis - A New Numerical Model for the Statics and Dynamics of Block Systems. PhD Thesis. University of California, Berkeley, USA, 1988.Search in Google Scholar
 Eberhardt E., Stead D., Coggan J., Numerical analysis of initiation and progressive failure in natural rock slopes-the 1991 Randa rockslide. Int. J. Rock Mech. Min. Sci., 2004, 41(1), 69-87.10.1016/S1365-1609(03)00076-5Search in Google Scholar
 Yasitli N., Unver B., 3D numerical modeling of longwall mining with top-coal caving. Int. J. Rock Mech. Min. Sci., 2005, 42(2), 219-235.10.1016/j.ijrmms.2004.08.007Search in Google Scholar
 Chamberlain S., Modarres M., Compressed natural gas bus safety: a quantitative risk assessment. Risk Anal., 2005, 25(2), 377-387.10.1111/j.1539-6924.2005.00596.xSearch in Google Scholar PubMed
 Loucks D.P., van Beek E., with contributions from, Stedinger J.R., Dijkman J.P.M., Villars M.T., Water resources systems planning and management: an introduction to methods, models and applications. Princeton Hall, Englewood Cli˙s, 2005.Search in Google Scholar
 Amadei B., Stephansson O., Rock stress and its measurement. Chapman, Hall, London, 1997.10.1007/978-94-011-5346-1Search in Google Scholar
 Hudson J.A., Cornet F.H. Special issue on rock stress estimation Int. J. Rock Mech. Min. Sci., 2003, 40(7-8), 955.10.1016/j.ijrmms.2003.08.001Search in Google Scholar
 Martin C., Kaiser P., Christiansson R., Stress, instability and design of underground excavations. Int. J. Rock Mech. Min. Sci., 2003, 40(7), 1027-1047.10.1016/S1365-1609(03)00110-2Search in Google Scholar
 Ioanna I., Dimitrios R., Theodora P., Paris T., Geotechnical and mineralogical properties of weak rocks from Central Greece. Open Geosci., 2009, 1(4), 431-442.10.2478/v10085-009-0029-0Search in Google Scholar
 Kamberis E., Bathrellos G., Kokinou E., Skilodimou, H., Correlation between the structural pattern and the development of the hydrographic network in a portion of the Western Thessaly Basin (Greece). Open Geosci., 2012, 4(3), 416-424.10.2478/s13533-011-0074-7Search in Google Scholar
 Scaioni M., Longoni L., Melillo, V., Papini M., Remote sensing for landslide investigations: An overview of recent achievements and perspectives. Remote Sens., 2014, 6, 9600-9652.10.3390/rs6109600Search in Google Scholar
 Cuenca M.C., Hooper A.J., Hanssen R.F., Surface deformation induced by water influx in the abandoned coal mines in Limburg, the Netherlands observed by satellite radar interferometry. J. Appl. Geophys., 2013, 88, 1-11.10.1016/j.jappgeo.2012.10.003Search in Google Scholar
 Chang L., Hanssen R.F., Detection of cavity migration and sinkhole risk using radar interferometric time series, Remote Sens. Environ., 2014, 147, 56-64.10.1016/j.rse.2014.03.002Search in Google Scholar
 Bersezio R., Fornaciari M., Gelati R., Geological map of the Southalpine foothills between Brianza and Iseo Lake. Mem. Soc. Geol. It., 1990, 45, 107-110.Search in Google Scholar
© 2016 L. Longoni et al., published by De Gruyter Open
This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 3.0 License.