Integration of ground-penetrating radar and microgravimetric methods to map shallow caves

Ground-penetrating radar (GPR) and microgravimetric surveys have been conducted in the southern Jura mountains of western Switzerland in order to map subsurface karstic features. The study site, La Grande Rolaz cave, is an extensive system in which many portions have been mapped. By using small station spacing and careful processing for the geophysical data, and by modeling these data with topographic information from within the cave, accurate interpretations have been achieved.The constraints on the interpreted geologic models are better when combining the geophysical methods than when using only one of the methods, despite the general limitations of two-dimensional (2D) profiling. For example, microgravimetry can complement GPR methods for accurately delineating a shallow cave section approximately 10×10 m² in size. Conversely, GPR methods can be complementary in determining cavity depths and in verifying the presence of off-line features and numerous areas of small cavities and fractures, which may be difficult to resolve in microgravimetric data.     

Geophysical Survey and Archaeological Data at Masseria Grasso (Benevento,Italy)

The use of geophysical methods in metrology is a significant tool within the wide research topic of landscape archaeology context. Since 2011, the Ancient Appia Landscapes Project aims to recognize dynamics, shapes and layout of the ancient settlement located along the Appia road east of Benevento, and cyclical elements and human activities that influenced the choice of landscapes. The integration of geophysical data with an archaeological infra-site analysis allowed us to investigate the area of Masseria Grasso, about 6 km from Benevento (Campania region, Italy). In this framework, an archaeogeophysical approach (Geomagnetic and Ground Penetrating Radar) was adopted for detecting anomalies potentially correlated with buried archaeological evidences. The geomagnetic results have given a wide knowledge of buried features in a large survey highlighting significant anomalies associated with the presence of buildings, roads and open spaces. These geophysical results permitted us to define the first archaeological excavations and, successively, a detailed Ground Penetrating Radar approach has been provided highlighting the rooms and paved spaces. The overlap between archaeological dataset and geophysical surveys has also allowed recognizing the path of the ancient Appia road near the city of Benevento and hypothesize the settlement organization of the investigated area, which has been identified with the ancient Nuceriola.     

SIMULATION STUDY OF ROAD-CUT EFFECTS ON SLOPE STABILITY

Landslides and rock falls along the highway are common geological hazards in Southwest China. As an influencing factor on potential landslides behavior, roads or distance to roads have been successfully used in landslide susceptibility assessments in mountainous area. However, the relationship between the road-cut and the slope stability is not clear. Therefore, we performed two-dimensional slope stability calculation using the general limit equilibrium (GLE)method incorporated in the software SLOPE/W of GeoStudio for stability analysis of slopes. Our studies show that the man-made roads influence on the slope stability mainly exists in two ways:One is to create a new steep slope, which will result in rock falls and shallow landslides along the roads; the other is to influence the stability of the original slope, which will result in comparatively huge landslides. For the latter, our simulation study reveals that the road location, namely at which part of a natural slope to construct a road is important for the slope stability. For a natural slope with a potential slip surface, if a road is constructed at or near the slope toe where the potential slip surface surpasses, it will greatly degrade the slope's factor of safety (Fs) and make the slope unstable; however, if a rode-cut is near the top of the slope, it will increase the slope's Fs and make the slope more stable. The safety location is different for different slope angle, steeper slope needs a higher location for a safety road-cut in comparison with gentle slopes. Moreover, the slope stability decreases when loading a seismic force and it varies with the slope angle. Firstly, the Fs decreases when the slope angle increasing, and when the slope angle reaches 45°, the Fs then becomes greater with the slope angle increasing.     

Lithological classification assisted by the joint inversion of electrical and seismic data at a control site in northeast Mexico

In this paper, evidence is presented that the combination of geospectral images and geophysical signatures (resistivity–velocity cross-plots) is a good tool to provide a natural visualization of the distribution and variations of lithological features in a test site. This was confirmed by the correlation between the electrical resistivity and seismic velocity values obtained after cross-gradient joint inversion at two profiles and geotechnical information provided by shallow boreholes in a site located in the Earth Sciences School grounds in Linares, Northeastern Mexico. The results obtained from this study show how the cross-gradient joint inversion facilitates the analysis of hydrological estimates and assists in lithological classification of subsurface materials.     

Effects of forest roads on the hydrological response of a small‐scale mountain watershed in Greece

The effects of land use changes on the ecology and hydrology of natural watersheds have long been debated. However, less attention has been given to the hydrological effects of forest roads. Although less studied, several researchers have claimed that streamflow changes related to forest roads can cause a persistent and pervasive effect on hillslope hydrology and the functioning of the channel system. The main potential direct effects of forest roads on natural watersheds hydrologic response are runoff production on roads surfaces due to reduced infiltration rates, interruption of subsurface flow by road cutslopes and rapid transfer of the produced runoff to the stream network through roadside ditches. The aforementioned effects may significantly modify the total volume and timing of the hillslope flow to the stream network. This study uses detailed field data, spatial data, hydro‐meteorological records, as well as numerical simulation to investigate the effects of forest roads on the hydrological response of a small‐scale mountain experimental watershed, which is situated in the east side of Penteli Mountain, Attica, Greece. The results of this study highlight the possible effects of forest roads on the watersheds hydrological response that may significantly influence direct runoff depths and peak flow rates. It is demonstrated that these effects can be very important in permeable watersheds and that more emphasis should be given on the impact of roads on the watersheds hydrological response. Copyright © 2014 John Wiley & Sons, Ltd.     

Geophysical exploration for geothermal low enthalpy resources in Lipari Island, Italy

Integrated geophysical surveys were performed in two sites, Fossa di Fuardo and Terme di San Calogero in Lipari Island, Southern Italy with the intent of the exploration of low-enthalpy geothermal fluids. Both sites show strong geochemical and geologic evidences of hydrothermal activity. The geophysical methods consist of two microgravimetric surveys, two 2D geoelectric profiles, a seismic reflection profile and a five seismic refraction profiles. The seismic methods allowed us to locate the main subsurface seismic discontinuities and to evaluate their geometrical relationships. The gravity field was used to constraint the seismic discontinuities, while the electric prospecting let discriminate more conductive areas, which could correspond to an increase in thermal fluid circulation in the investigated sites.The results obtained by the different geophysical methods are in good agreement and permit the definition of a reliable geo-structural model of the subsurface setting of the two investigated areas. A low-enthalpy geothermal reservoir constituted by a permeable pyroclastic and lava sequence underlying two shallow impermeable formations was found at Fossa del Fuardo. The reservoir is intersected by some sub-vertical faults/fractures that probably play an important role in convoying the thermal water up to the surface. At the other site, Terme di S. Calogero, the geophysical surveys showed that an intense circulation of fluids affects the subsurface of the area. This circulation concentrates along a ENE-trending fault located at a little distance from the thermal resort. The hot fluids may upraise along the fault if the width of the ascent area is smaller than 20 m.     

Detection of sinkholes using 2D electrical resistivity imaging

Sinkholes in dolomitic areas are notoriously difficult geophysical targets, and selecting an appropriate geophysical solution is not straightforward. Electrical resistivity imaging, or tomography (RESTOM) is well suited to mapping sinkholes because of the ability of the technique for detecting resistive features and discriminating subtle resistivity variations. RESTOM surveys were conducted at two sinkhole sites near Pretoria, South Africa. The survey areas are located in the dolomites of the Lyttelton Formation, which forms part of the Malamani Subgroup and Chuniespoort Group of the Transvaal Supergroup. The survey results suggest that RESTOM is an ideal geophysical tool to aid in the detection and monitoring of sinkholes and other subsurface cavities.     

Electrical resistivity and seismic refraction tomography to detect buried cavities

Near‐surface cavities can pose serious hazards to human safety, especially in highly urbanized town centres. The location of subsurface voids, the estimation of their size and the evaluation of the overburden thickness are necessary to assess the risk of collapse. In this study, electrical resistivity tomography (ERT) and seismic refraction tomography data are integrated in a joint interpretation process for cavity location in the city of Rome. ERT is a well established and widely employed method for cavity detection. However, additional information provided by seismic refraction tomography is capable of eliminating some potential pitfalls in resistivity data interpretation. We propose that the structure of the cavities defined by ERT can be used as a base to optimize seismic refraction tomography investigations within the framework of a joint interpretation process. Data integration and the insertion of a priori information are key issues for reducing the uncertainties associated with the inversion process and for optimizing both acquisition procedures and computation time. Herein, the two geophysical methods are tested on both synthetic and real data and the integration of the results is found to be successful in detecting isolated cavities and in assessing their geometrical characteristics. The cavity location inferred by geophysical non‐invasive methods has been subsequently confirmed by direct inspection.     

浅层地震资料解释陷阱(英文)

高分辨率浅层地震方法是在近地表调查中使用最为广泛的方法。然而,在许多情况下,地震资料的解释经常会出现错误。在本文中,我们介绍了三个例子,分析了造成P波,SH波,多道的面波(MASW)地震资料解释的错误原因,大都是由于在表面或地下条件约束不确当引起的。第一个例子是P波反射剖面上的一个波的特征被解释为浅层断裂带,但后来证实它是由高水平的背景噪音引起的,因为采集测线通过了一个公路交叉口。第二个例子是SH波反射地震剖面上一个波特征被解释为是逆倾向滑断层,但有针对性的钻探表明,它是一个侵入到基岩面的一个深层局部侵蚀。最后,第三个例子,MASW调查剖面上,一个陡倾特征一开始被解释为基岩谷。然而,后来的钻探表明这是一个非常软的湖泊沉积物,后者严重损坏了应用面波频段。虽然最初的解释是不正确的,但这刺激地球物理学家和地质学家之间的讨论,并强调地球物理数据采集的时候,采集之前以及采集之后需要科学家之间有意义的合作与讨论。     

Coastal and inland karst morphologies driven by sea level stands: a GIS based method for their evaluation

Sea level is the base level for groundwater circulation in coastal aquifers. The evolution of karst surface landforms and subsurface drainage systems in these aquifers has been conditioned in geological time by tectonics and glacio‐eustatic sea‐level changes. Present morpho‐structural settings and the type/distribution of karst surface and subsurface forms have developed in different carbonate formations according to differences in lithology, climate and exposure time, all driving the intensity of morphologic and karst processes. The repeated and significant changes of groundwater level linked to ‘sea‐level changes’ have had the most important role in driving the continuous evolution of karstic drainage systems, and has resulted in most cases in a multiphase karst. This study aims at defining a general method for identifying, in karst coastal settings, the elevations of flat or low topographic gradient surfaces (using morphometric analysis of Digital Elevation Models (DEMs) and geographical information systems (GISs), and their comparison with elevations of distinctive karstic levels (passages, lateral solution cavities) observed in vertical shafts and horizontal caves. Of the elevations of flat or low topographic gradient surfaces only those agreeing, within ±10 m or ±20 m, with elevation ranges marked by the high frequency of distinctive karst levels were considered as representative of the more probable past sea‐level stands. The method is applied to a regional coastal carbonate formation in southern Italy, by using a 10 m DEM and information on 140 complex caves and 85 shafts. Of the 15 elevations indicated by DEM analysis [620, 600, 470, 450, 425, 385, 355, 315, 270, 250, 205, 180, 150, 110, and 70 m above sea level (a.s.l.)], 13 match clearly those highlighted by significant frequencies of distinctive karstic levels. These elevations are validated by comparison to the elevation of terraces and karst plains indicated in the literature. Copyright © 2012 John Wiley & Sons, Ltd.     

Runoff and road erosion at the plot and road segment scales,St John,US Virgin Islands

Previous studies have identified unpaved roads as the primary source of erosion on St John in the US Virgin Islands, but these studies estimated road erosion rates only as annual averages based primarily on road rill measurements. The goal of this project was to quantify the effect of unpaved roads on runoff and sediment production on St John, and to better understand the key controlling factors. To this end runoff and sediment yields were measured from July 1996 to March 1997 from three plots on naturally vegetated hillslopes, four plots on unpaved road surfaces and two cutslope plots. Sediment yields were also measured from seven road segments with contributing areas ranging from 90 to 700 m². With respect to the vegetated plots, only the two largest storm events generated runoff and there was no measurable sediment yield. Runoff from the road surface plots generally occurred when storm precipitation exceeded 6 mm. Sediment yields from the four road surface plots ranged from 0·9 to 15 kg m⁻² a⁻¹, and sediment concentrations were typically 20–80 kg m⁻³. Differences in runoff between the two cutslope plots were consistent with the difference in upslope contributing area. A sprinkler experiment confirmed that cross‐slope roads intercept shallow subsurface stormflow and convert this into surface runoff. At the road segment scale the estimated sediment yields were 0·1 to 7·4 kg m⁻² a⁻¹. Road surface runoff was best predicted by storm precipitation, while sediment yields for at least three of the four road surface plots were significantly correlated with storm rainfall, storm intensity and storm runoff. Sediment yields at the road segment scale were best predicted by road surface area, and sediment yields per unit area were most strongly correlated with road segment slope. The one road segment subjected to heavy traffic and more frequent regrading produced more than twice as much sediment per unit area than comparable segments with no truck traffic. Particle‐size analyses indicate a preferential erosion of fine particles from the road surface and a rapid surface coarsening of new roads. Published in 2001 by John Wiley & Sons, Ltd.     

Integrated geophysical and chemical study of saline water intrusion

Surface geophysical surveys provide an effective way to image the subsurface and the ground water zone without a large number of observation wells. DC resistivity sounding generally identifies the subsurface formations-the aquifer zone as well as the formations saturated with saline/brackish water. However, the method has serious ambiguities in distinguishing the geological formations of similar resistivities such as saline sand and saline clay, or water quality such as fresh or saline, in a low resistivity formation. In order to minimize the ambiguity and ascertain the efficacy of data integration techniques in ground water and saline contamination studies, a combined geophysical survey and periodic chemical analysis of ground water were carried out employing DC resistivity profiling, resistivity sounding, and shallow seismic refraction methods. By constraining resistivity interpretation with inputs from seismic refraction and chemical analysis, the data integration study proved to be a powerful method for identification of the subsurface formations, ground water zones, the subsurface saline/brackish water zones, and the probable mode and cause of saline water intrusion in an inland aquifer. A case study presented here illustrates these principles. Resistivity sounding alone had earlier failed to identify the different formations in the saline environment. Data integration and resistivity interpretation constrained by water quality analysis led to a new concept of minimum resistivity for ground water-bearing zones, which is the optimum value of resistivity of a subsurface formation in an area below which ground water contained in it is saline/brackish and unsuitable for drinking.     

Runoff generation and sediment production on unpaved roads,footpaths and agricultural land surfaces in northern Thailand

Rainfall simulation was used to examine runoff generation and sediment transport on roads, paths and three types of agricultural fields in Pang Khum Experimental Watershed (PKEW), in mountainous northern Thailand. Because interception of subsurface flow by the road prism is rare in PKEW, work focused on Horton overland flow (HOF). Under dry antecedent soil moisture conditions, roads generated HOF in c. 1 min and have event runoff coefficients (ROCs) of 80 per cent, during 45 min, c. 105 mm h⁻¹ simulations. Runoff generation on agricultural fields required greater rainfall depths to initiate HOF; these surfaces had total ROCs ranging from 0 to 20 per cent. Footpaths are capable of generating erosion‐producing overland flow within agricultural surfaces where HOF generation is otherwise rare. Paths had saturated hydraulic conductivity (K_s) values 80–120 mm h⁻¹ lower than those of adjacent agricultural surfaces. Sediment production on roads exceeded that of footpaths and agricultural lands by more than eight times (1·23 versus < 0·15 g J⁻¹). Typically, high road runoff volumes (owing to low K_s, c. 15 mm h⁻¹) transported relatively high sediment loads. Initial road sediment concentrations exceeded 100 g l⁻¹, but decayed with time as loose surface material was removed. Compared with the loose surface layer, the compacted, underlying road surface was resistant to detachment forces. Sediment concentration values for the road simulations were slightly higher than data obtained from a 165 m road section during a comparable natural event. Initial simulation concentration values were substantially higher, but were nearly equivalent to those of the natural event after 20 min simulation time. Higher sediment concentration in the simulations was related to differences in the availability of loose surface material, which was more abundant during the dry‐season simulations than during the rainy season natural event. Sediment production on PKEW roads is sensitive to surface preparation processes affecting the supply of surface sediment, including vehicle detachment, maintenance activities, and mass wasting. The simulation data represent a foundation from which to begin parameterizing a physically based runoff/erosion model to study erosional impacts of roads in the study area. Copyright © 2000 John Wiley & Sons, Ltd.     

地震-电成像联合探测试验——以淄博市活断层探测为例

在城市浅部活断层探测中,地震和电成像是2种主要的地球物理勘探手段。它们既可独立开展工作,也可在地质构造复杂地段进行联合探测,以更合理地解释探测区域的地质构造。文中利用电阻率和地震纵波速度之间的简单关系,进行了地震-电成像联合反演的初步探索。结果表明,通过参考电成像图像的电阻率分布,可获得合理的准二维地震波叠加速度模型,最终获得探测区域合理的地震时/深剖面,深化对地质构造的认识和提供合理的解释     

Persistence of road runoff generation in a logged catchment in Peninsular Malaysia

Measurements of saturated hydraulic conductivity (K_s) and diagnostic model simulations show that all types of logging road/trail in the 14·4 ha Bukit Tarek Experimental Catchment 3 (BTEC3) generate substantial Horton overland flow (HOF) during most storms, regardless of design and level of trafficking. Near‐surface K_s(0–0·05 m) on the main logging road, skid trails and newly constructed logging terraces was less than 1, 2 and 34 mm h^?1, respectively. Near‐surface K_s on an abandoned skid trail in an adjacent basin was higher (62 mm h^?1), owing to the development of a thin organic‐rich layer on the running surface over the past 40 years. Saturated hydraulic conductivity measured at 0·25 m below the surface of all roads was not different (all <6 mm h^?1) and corresponded to the K_s of the adjacent hillslope subsoil, as most roads were excavated into the regolith more than 0·5–1 m. After 40 years, only limited recovery in near‐surface K_s occurred on the abandoned skid trail. This road generated HOF after the storage capacity of the upper near‐surface layer was exceeded during events larger than about 20 mm. Thus, excavation into low‐K_s substrate had a greater influence on the persistence of surface runoff production than did surface compaction by machinery during construction and subsequent use during logging operations. Overland flow on BTEC3 roads was also augmented by the interception of shallow subsurface flow traveling along the soil–saprolite/bedrock interface and return flow emerging from the cutbank through shallow biogenic pipes. The most feasible strategy for reducing long‐term road‐related impacts in BTEC3 is limiting the depth of excavation and designing a more efficient road network, including minimizing the length and connectivity of roads and skid trails. Copyright © 2007 John Wiley & Sons, Ltd.     

The origin,typology, spatial distribution and detrimental effects of the sinkholes developed in the alluvial evaporite karst of the Ebro River valley downstream of Zaragoza city (NE Spain)

Three types of sinkhole have been mapped in a 50 km² stretch of the Ebro River valley downstream of Zaragoza: large collapse sinkholes, large shallow subsidence depressions and small cover-collapse sinkholes. The sinkholes relate to the karstification of evaporitic bedrock that wedges out abruptly downstream, giving way to a shale substratum. Twenty-three collapse sinkholes, up to 50 m in diameter by 6 m deep, and commonly hosting saline ponds, have been identified in the floodplain. They have been attributed to the upward stoping of dissolutional cavities formed within the evaporitic bedrock by rising groundwater flows. Twenty-four large shallow subsidence depressions were mapped in the floodplain. These may reach 850 m in length and were formed by structurally controlled interstratal karstification of soluble beds (halite or glauberite? and gypsum) by rising groundwater flow and the progressive settlement of the overlying bedrock and overburden sediments. A total of 447 small cover-collapse, or dropout, sinkholes have been recognized in a perched alluvial level along the southern margin of the valley. These sinkholes result from the upward propagation of voids through the alluvial mantle caused by the downward migration of detrital sediments into dissolutional voids. The majority of these sinkholes, commonly 1·5–2 m in diameter, are induced by human activities. Over the karstic bedrock, there is a significant increase in sinkhole density downstream. This is interpreted as being a result of the evaporitic bedrock wedging out and the convergence of the groundwater flow lines in the karstic aquifer. The collapse sinkholes in this area, locally with a probability of occurrence higher than 140 sinkholes/km²/year, cause substantial damage to the linear infrastructures, buildings and agriculture, and they might eventually cause the loss of human lives. Copyright © 2006 John Wiley & Sons, Ltd.     

深部金属矿产资源地球物理勘查与应用

随着找矿勘探工作程度的不断深入,许多露天的和近地表的金属矿产资源已基本上被查明,在地表浅部(第一深度空间:0~500 m深度)找到大型或超大型金属矿床的难度将必越来越大,而金属矿产资源的短缺已在日益加剧.为了缓解当今的资源危机,必须进行深部找矿勘探(第二深度空间:500~2000 m深度).这是因为在地壳深部具有良好的成矿环境和找矿潜力,地球物理勘探技术的发展与成效已使得深部找矿成为可能.本文通过分析和讨论国内外典型金属矿产资源找矿过程的实践表明:地球物理勘探方法发挥着重要作用,因为它具有大探测深度、高精度和高分辨率的特点,可为深部金属矿勘查提供有效信息,是第二深度空间找矿勘探的有力手段.近年来,国内外的实践表明,深部矿产资源的地球物理勘探取得了尚在不断取得重要成效.     

2D magnetic resonance tomography applied to karstic conduit imaging

Karstic conduits play a crucial role for water supply in many parts of the world. However, the imaging of such targets is generally a difficult task for most geophysical methods. Magnetic Resonance Sounding (MRS) is a geophysical method designed for imaging of water bearing structures. Initially, MRS was developed for characterizing horizontally stratified aquifers. However, when applying a 1D MRS measuring setup to the imaging of 2D–3D targets, the size of which may be much smaller than the loop, the accuracy and the lateral resolution may not be sufficient. We have studied the possibility of simultaneously processing several MRS aligned along a profile to perform a Magnetic Resonance Tomography (MRT). This work emphasizes the gain of resolution for 2D–3D imagery of MRT versus the interpolation of 1D inversion results of MRS along the same profile. Numerical modelling results show that the MRT response is sensitive to the size and location of the 2D target in the subsurface. Sensitivity studies reveal that by using the coincident transmitting/receiving (TX/RX) setup and shifting the loop around the anomaly area, the depth, section and position of a single karstic conduit with a size smaller than the MRS loop size can be resolved. The accuracy of the results depends on the noise level and signal level, the latter parameter being linked to the depth and volume of the karstic conduit and the water content in the limestone matrix. It was shown that when applying MRT to the localization of 2D anomalies such as karstic conduits, the inclination of the geomagnetic field, the orientation of the MRT profile and the angle of crossover of the conduit by the MRT profile must be taken into account. Otherwise additional errors in interpretation should be expected. A 2D inversion scheme was developed and tested. Both numerical and experimental results confirm the efficiency of the developed approach.     

Gravity and magnetic investigation of maar volcanoes,Auckland volcanic field,New Zealand

Detailed gravity and aeromagnetic data over maars in the Auckland volcanic field reveal contrasting anomalies, even where surface geology is similar. Pukaki and Pukekiwiriki, almost identical maars marked by sediment-filled craters and tuff rings, have gravity and magnetic anomalies of − 6 g.u. and 20 nT, and 8 g.u. and 160 nT, respectively. The Domain and Waitomokia maars, with similar tuff rings but each with a small central scoria cone, have gravity and magnetic anomalies of 32 g.u. and 300 nT, and 21 g.u. and 310 nT, respectively. These differences in geophysical expression are attributed to varying volumes of dense, magnetic basalt in the form of shallow bowl-shaped bodies up to several hundreds of metres in diameter and up to 140 m thick beneath the maar centres. These bodies are interpreted as solidified magma that ponded into early-formed phreatomagmatic explosion craters. Where magma supply was limited relative to groundwater availability, no residual subsurface basalt occurs (as at Pukaki); continued magma supply, but limited groundwater, resulted in ponding (e.g. at Pukekiwiriki) and eventually the building of a scoria cone (as at Domain and Waitomokia). There is no evidence in these geophysical data for diatreme structures below the maars or for shallow and/or extensive feeder dykes associated with these maars. If diatreme structures do occur, their lack of geophysical signature must be a consequence of either their small geophysical contrast with host Miocene sediments and/or masking by the stronger anomalies associated with the subsurface basalt. In addition, any magma conduits appear to be confined centrally beneath the maars, at least to shallow depths (upper 100 m).     

Road Network and Damaged Buildings in Urban Areas: Short and Long-term Interaction

In this paper, a model able to analyse the seismic behaviour of road network in urban areas, considering interaction between buildings and roads is presented. Damage to buildings and short-term countermeasures, such as propping, can affect roads and even block them, reducing capacity of the road network. Two successive phases are considered. In the first, immediately after the seismic event, emergency services have to reach the relevant buildings. Here, the network topology is of interest. In the second, in the long term, network capacity is not yet completely restored and network demand has to account for displaced people due to unusable buildings. In this case, road serviceability is of interest. In order to consider uncertainties on building behaviour, a probabilistic approach is developed and the network is analysed by means of Monte Carlo simulation. The methodology is applied to the municipality of Potenza (southern Italy) evaluating in the short term, the probability of strategic buildings are not connected and, in the long term, road serviceability.