Next Article in Journal
A Fisheye Image Matching Method Boosted by Recursive Search Space for Close Range Photogrammetry
Previous Article in Journal
Semantic Segmentation of SLAR Imagery with Convolutional LSTM Selectional AutoEncoders
 
 
Comment published on 8 December 2020, see Remote Sens. 2020, 12(24), 4011.
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Pre-Collapse Space Geodetic Observations of Critical Infrastructure: The Morandi Bridge, Genoa, Italy

1
NASA Jet Propulsion Laboratory, California Institute of Technology, Pasadena, CA 91109, USA
2
Department of Architecture and Civil Engineering, University of Bath, Bath BA2 7AY, UK
3
Raser limited, 33050 Ruda, UD, Italy
4
Italian Space Agency, 75100 Matera (MT), Italy
5
Geoportale Nazionale, Italian Ministry of the Environment, 00147 Rome, Italy
*
Author to whom correspondence should be addressed.
Remote Sens. 2019, 11(12), 1403; https://0-doi-org.brum.beds.ac.uk/10.3390/rs11121403
Submission received: 30 April 2019 / Revised: 6 June 2019 / Accepted: 7 June 2019 / Published: 12 June 2019
(This article belongs to the Section Urban Remote Sensing)

Abstract

:
We present a methodology for the assessment of possible pre-failure bridge deformations, based on Synthetic Aperture Radar (SAR) observations. We apply this methodology to obtain a detailed 15-year survey of the Morandi bridge (Polcevera Viaduct) in the form of relative displacements across the structure prior to its collapse on August 14th 2018. We generated a displacement map for the structure from space-based SAR measurements acquired by the Italian constellation COSMO-SkyMed and the European constellation Sentinel-1A/B over the period 2009–2018. Historical satellite datasets include Envisat data spanning 2003–2011. The map reveals that the bridge was undergoing an increased magnitude of deformations over time prior to its collapse. This technique shows that the deck next to the collapsed pier was characterized since 2015 by increasing relative displacements. The COSMO-SkyMed dataset reveals the increased deformation magnitude over time of several points located near the strands of this deck between 12th March 2017 and August 2018.

Graphical Abstract

1. Introduction

Over the past 20 years Interferometric Synthetic Aperture Radar (InSAR) and Persistent Scatterer Interferometry (PSI) have established themselves as a methodology capable of monitoring surface deformations in urban environments, thanks to the presence and density of coherent radar targets. The imaging capability of the current generation of SAR constellations (e.g., COSMO-SkyMed/SAOCOM, TerraSAR-X/PAZ and Senintel-1A/B) has been significantly improved, to the extent that SAR images are now provided with meter to sub-meter spatial resolution, and are acquired with a revisit time of a few days, in most cases according to a regular acquisition schedule which over the years can result in hundreds of images per individual site [1].
There is a wealth of scientific literature demonstrating the benefits of InSAR to monitor a wide spectrum of urban processes, such as subsidence [2,3,4,5], foundation settlement [1,6,7], deformation of new buildings [8], surface effects due to tunneling [9,10,11], ground motions associated to faults and local tectonics [12,13,14]. The implementation on real case studies proves the increasing use of these techniques on a routine basis and how they can influence the management of infrastructure assets and preservation of the built environment [1,15,16].
However, when the InSAR technique is applied to urban areas, the interpretation of phenomena occurring at a local scale or at the level of a single structure remains challenging. While a stand-alone InSAR multi-temporal analysis can be used as a straightforward indicator of structural instability [17,18], only a synergistic approach allows to make sense of information inferred from different satellite datasets [19]. Moreover, probabilistic procedures are required to infer the three-dimensional field of displacement enabling the structural evaluation of individual assets [1].
Here we present a methodology for the assessment of single structure deformations, based on Synthetic Aperture Radar (SAR) observations and Markov Chain Monte Carlo (MCMC) approach. We applied this approach to the Morandi bridge producing a comprehensive multi-sensor time-series deformation map generated from space-based SAR measurements. First, we show the average displacements of the main structural elements of the bridge and the local displacements of points which are exhibiting an increased rate of deformation. Then, we model the deformation inferred from multi-geometry observations using a MCMC approach to solve for three-dimensional deformation of the single decks and their base. The evaluation of the uncertainties associated with the method is presented in the methodology section. Finally, we discuss the dataset records in light of their compatibility with possible collapse mechanisms, showing the potential operational use of InSAR millimeter accuracy time-series for early warning, hindcasting and forensic engineering applications.

1.1. The Morandi Bridge

The Polcevera viaduct (also known as Morandi bridge) is a prestressed concrete cable-stayed bridge composed of 11 spans (Figure 1) for a total length of about 1300 m [20]. It connected the motorways Milano-Genova and Genova-Savona, in northern Italy, crossing the Polcevera river on an area composed of marine and alluvial sand. Starting from the southern side, the bridge is supported by two single piers (1 to 2), six V-shaped piers (3 to 8) and three A-shaped frames (9 to 11) (Figure 1). Piers 3 to 8 consist of two inclined piers connected by a double cantilever deck of variable length [21], while piers 9 to 11 were designed as independent “balanced systems” made of a three-span continuous deck resting on four supports. Each of the three-span decks is supported externally by the anchorage of two prestressed stay-cables passing over a 90 m high mast. The anchorages are connected by transverses located at a distance of 10 m from the cantilever tips. The two internal supports are the top end of two pairs of diverging inclined legs. Each balanced system is supported by a concrete ribbed foundation raft on concrete piles. On August 14th 2018, the system number 9, including an approximately 240 m long section of the adjacent deck suddenly collapsed, causing the death of 43 people.
The bridge was built between 1960 and 1966. In 1993, a detailed inspection of the balanced system number 11 reported oxidation of the metallic membrane protecting the strands together with humidity, absence of grout, slacked or sheared strands, and signs of tensile rupture due to steel cross-section reduction [22]. The 1993 report informed the structure management authorities of the ongoing material degradation and micro-cracking phenomena occurring within the prestressed concrete [23]. The balanced system 11 was subject to retrofitting and many other maintenance activities were carried over during the years [24]. The strands replacement for the balanced system 9 was planned to start two months after the collapse [23].

2. Materials and Methods

The methodology adopted in this paper followed two main steps including SAR multi-temporal interferometric analysis (MT-InSAR) and Markov Chain Monte Carlo three-dimensional inversion.

2.1. Multi-Temporal InSAR

Multi-temporal InSAR combines several acquisitions highlighting the temporal behavior of the relative displacement between each acquisition. In this paper surface deformation measurements have been generated using the SARPROZ software extending the standard linear PS technique [25,26] to solve for non-linear motion with no a priori information. SRTM has been used as a supporting DEM for the PS processing. COSMO-SkyMed (CSK) and Sentinel-1 (SNT) data have been specifically processed for this study while the Envisat time-series are part of a nationwide monitoring program [27]. Table 1 shows information related to the dataset processed and used in this study. We consider a set of N+1 coregistered single look complex (SLC) SAR images and generates a set of N single look differential interferograms (for both SNT and CSK) with respect to a master imaged chosen maximizing the following coherence formula:
γ t o t a l = γ t e m p o r a l γ s p a t i a l γ d o p p l e r [ 1 f ( T T C ) ] [ 1 f ( B B C ) ] [ 1 f ( F D F D C ) ]
f ( x ) = { x ,   for   x 1 1 ,   for   x 1
where, B, T, FD and BC, TC, FDC are the perpendicular, temporal and doppler baselines and their critical values, respectively. The selected master images for CSK ascending, descending and SNT descending were acquired on 1st January 2015, 11th June 2015 and 20th July 2017, respectively. After selecting a master image, we applied the proposed algorithm only on pixels showing an inverse amplitude stability (i.e., ratio between the mean calibrated image amplitudes and its standard deviation) higher than 0.65, choosing 4 reference points as explained in the results section of this manuscript. The advantage of this approach is to consider only persistent scatterers (PS) that are unaffected by geometrical decorrelation, i.e., PS in interferograms characterized by perpendicular baseline larger than the critical baseline still contain useful phase information [25,26]. Specifically, our output consisted of four different time-series.
The interferometric phase difference between two neighbor pixels xr and xs in the interferogram i can be written as the sum of five terms:
Δ ϕ i ( x r , x s , T ) =   Δ C t o p o · B + Δ C V · T + Δ w r , s i
where Δ C t o p o and Δ C V are proportional to the DEM errors and mean velocity difference between pixels xr and xs, respectively. Δ w r , s contains the sum of three contributions related to non-linear motion and atmospheric phase screen (APS) noise. Assuming | Δ w r , s | < π for neighbor pixels it is possible to estimate Δ C t o p o and Δ C V by maximizing the absolute value of the temporal coherence | ξ P S | for neighbor pixels:
| ξ P S | = | i = 1 N e j ( Δ C t o p o · B + Δ C V · T Δ ϕ i ( x r , x s , T ) ) N |
The following steps consisted in unwrapping every interferogram and estimating the phase term Δ w r , s i . Once the data were unwrapped we estimated the low pass component using a three sample temporal base-line weighted moving average and assumed this residual phase term as an estimation of the non-linear motion contribution. These steps have been performed considering only neighbor pixels. The same processing was then applied considering only differences between the reference point and all persistent scatterers including the estimated APS in Equation (4).
The accuracy of the estimated velocities and DEM errors can be derived by a regression analysis considering the irregularly sampled data (perpendicular baseline and temporal baseline in the case of the CSK constellation. As in [28] we define the standard deviation of scatterer height ( σ Δ h 2 ) and velocity ( σ Δ v 2 ) :
σ Δ h 2 ( λ R s i n θ 4 π ) 2 σ ϕ 2 N σ B 2
σ Δ v 2 ( λ 4 π ) 2 φ σ ϕ 2 N σ B t 2
where R is the target-sensor distance, θ is the incidence angle, λ is the sensor wavelength, σ ϕ 2 is the phase noise variance supposed independent from the acquisition and σ B 2 , σ B t 2 are the variance of the perpendicular and temporal baselines respectively. For further details on the technique we refer the reader to [25,26,28].

2.2. Markov Chain Monte Carlo

We used a Markov chain Monte Carlo (MCMC) approach [1] (Algorithm 1) to solve for three-dimensional deformation and their relative errors. We estimated the East North and Vertical velocity components (in mm/yr, Figure 2) for each pillar on the northern and southern side of the bridge, at the deck level and at each pillar ground level. Using Equation (7), we estimated for each incidence angle θ and heading angle φ (Table 2), the probability density functions (PDF) from histograms of retained solutions (Figure 2). We ran the MCMC algorithm iteratively until reaching three million kept solutions. The most probable parameters were estimated from the PDF together with 95% confidence intervals. Figure 2 shows the inverted parameters for the CSK and Sentinel-1A data. When too many parameters are inverted at the same time, the PDFs for several of the parameters do not form smooth, Gaussian-like distributions, but instead are multi-modal. We found a stable solution only when three independent line of sight were available (i.e., CSK ascending, descending and SNT descending). The input LOS parameters are listed in Table 1 while the obtained results are shown in Figure 2 and Figure 3.
V L O S = [ v A S C C S K v D S C C S K v D S C S N T ] = [ E sin ( θ A S C C S K ) cos ( φ A S C C S K ) + N sin ( θ A S C C S K ) sin ( φ A S C C S K ) + U   cos ( θ A S C C S K ) E sin ( θ D S C C S K ) cos ( φ D S C C S K ) + N sin ( θ D S C C S K ) sin ( φ D S C C S K ) + U   cos ( θ D S C C S K ) E sin ( θ D S C S N T ) cos ( φ D S C S N T ) + N sin ( θ D S C S N T ) sin ( φ D S C S N T ) + U   cos ( θ D S C S N T ) ]
Algorithm 1. Markov Chain Monte Carlo Algorithm steps
 1. Initialize E, N, U initial value Calculate CSK, ASC/DSC and SNT LOS projections δ L O S
 2. Calculate error χ 1 = ( V L O S δ L O S ) 2  
 3. Random increment in (E,N,U)
 4. Recalculate δ L O S and χ 2 = ( V L O S δ L O S ) 2
If uniformly distributed random number in the interval (0,1) < exp [ ( χ 1   χ 2 2 ) ]
 5. Then keep solution and iterate from points 1. to 5. until reaching 3 MLN of kept solutions
Else reiterate from points 3 to 5.

3. Results

We processed 130 and 148 images acquired in an ascending and descending geometry, respectively, from the Italian COSMO-SkyMed (CSK) constellation to infer ground deformation during 2009–2018. A total of 134 and 136 images in ascending and descending geometry, respectively, have been acquired by the European Sentinel-1A/B sensors starting in September 2015 and spanning 3 years. A total of 34 ascending and 24 descending images have been acquired by the Envisat sensor starting in 2003 and spanning 7 years, showing no relevant deformation of the bridge pillars during this time-span (Figure S1). A multi-temporal interferometric processing technique (MT-InSAR) was applied to the CSK and Sentinel-1A/B phase measurements to compute line of sight time-series of deformation exploiting point-like scatterers (PS) which are points located on the Earth surface and anthropogenic structures characterized by a stable electromagnetic signature in time [1]. While both the SNT and CSK descending datasets captured the temporal behavior of points located on the bridge deck, the SNT ascending dataset only recorded points located on the ground and for this reason has been discarded from our analysis (Figure 4 and Figure S2). Because of the differential nature of the MT-InSAR measurements, the time series of cumulative deformation were calculated relative to a reference area. Different reference areas were used, in order to highlight different phenomena characterizing the bridge deck and its base. In this analysis we used four reference areas (Figure 4, Panel A, Panel F): the first reference area was located on the west side of the bridge, close to pier 5 and 6, at the deck level (REF1), while the other three were positioned on the north east side of each balanced system 9–11 (REF2–4) (Figure 1A). REF1 was used as a reference area for analyzing the behavior of the entire bridge and the interaction of each balanced system with the ground, while REF2–4 were used for characterizing the three-dimensional displacement field of the deck, independently for each system. Time series of deformation using REF1 show that average deck deformation recorded from ascending and descending data were characterized by linear trends during the entire satellite acquisition campaign (grey pattern in Movies 1–5, see supplementary material). Movie 1 and Figure 5 show the time-series deformation for the descending line of sight (LOS) of the CSK satellites. The deformation refers to the northern side of the bridge for the three balanced systems at the deck level. In Movie 1 all the PS were located on the north side of the bridge (Figure 4, Figure 5). We found one PS (red dot in Movie 1 and Point CSK DSC 1 in Figure 1), located at the junction between the northwest deck level and the northwest strand, which was characterized by an abrupt acceleration starting March 12th 2017. The LOS deformation rate of this point increased by seven times from 10 ± 2.8 mm/yr to 70 ± 2.8 mm/yr away from the satellite, as shown in Movie 1 and Figure 1. The descending SNT dataset (Movie 2) recorded similar average trends (9.4 ± 2.8 mm/yr) for all the scatterers located on the northern side of the bridge, although there are no persistent scatterers in the SNT dataset covering the area of the bridge characterized by an anomalous speed increase. However, SNT descending scatterers located on the southern side of the balanced system 9 highlighted an area close to the junction between the southern west deck level and the southern west strand that was characterized by a motion toward the satellite at a rate of about 10.3 mm/yr and two deformation transient steps of 1 and 2 cm occurring on January 21st 2017 and April 28th 2018, respectively (Movie 3). These steps were incompatible with phase unwrapping errors only occurring at multiples of half wavelength (2.8 cm for C-band 5.6 cm Sentinel-1A/B wavelength). The balanced system 9 moved away from the satellite at a constant rate of 6.7 ± 2.8 mm/yr. PS recorded by the CSK ascending were located on the southern side of the balanced systems (Movie 3) and moved by about 6.6 ± 5.8 mm/yr away from the satellite, confirming the SNT average deformation trend of the deck southern side of the balanced system 9. As for CSK descending, we found two PS (red and blue dots in Movie 4) located at the junction between the south-west (Point CSK ASC 1 in Figure 1) and east deck level (Point CSK ASC 2 in Figure 1) and the south-west and east strand, respectively, which were characterized by an abrupt acceleration starting March 12th 2017. The LOS deformation rate of these points increases by 10 times from –3.3 ± 2.8 mm/yr to 33.3 ± 2.8 mm/yr away from the satellite for point CSK ASC1, and from –1.6 ± 2.8 mm/yr to 13.3 ± 2.8 mm/yr for point CSK ASC2 (Figure 1). Movies 6–8 show time-series of cumulative deformation at the ground level (using REF1) for CSK descending, SNT descending and CSK ascending, respectively. Given that neither GPS nor other ground truths were available, it is not possible to establish if REF1 was affected by displacement. This made the single balanced system analysis difficult to interpret when using REF1. We interpret this behavior later on when inferring the three-dimensional deformation field. We found that for the deck 9, PS located on the north and southern side moved by 1.2 ± 2.8 mm/yr and 2.8 ± 2.8 mm/yr towards the sensor, respectively, while PS located at the ground level moved by 1.5 ± 2.8 mm/yr away from the satellite. The points located on the north side of both decks 10 and 11 were characterized by a 0.6 ± 2.8 mm/yr motion toward the sensor along the LOS. An insufficient number of PS was found on the southern side of the deck 9 and 10, making the trend analysis not feasible in this area. The ground under pier 10 was moving at 0.5 ± 2.8 mm/yr away from the sensor, similarly to the ground under pier 9. A comprehensive interpretation of LOS displacements for each dataset remains a difficult task. For this reason, it is fundamental to infer three-dimensional displacement field and reproject our data into East, North and Vertical deformation trends. We analyzed LOS measurements from three independent datasets, using the reference areas REF2–4 to characterize the three-dimensional deformation of each deck prior to collapse and filter out deformation potentially characterizing REF 1. We modeled the multi-geometry inferred deformation using a Markov chain Monte Carlo (MCMC) approach to solve for three-dimensional deformation of the bridge deck. Each pier was designed as an independent balanced system made of a 171 m continuous deck connected to two simply-supported Gerber beam systems of 36 m. It is therefore reasonable to assume that each deck and beam can be approximated as an independent rigid block. Under these assumptions we infer three-dimensional deformation trends of each side (south and north) of the deck and calculate its deformation along three main independent axes. We use an MCMC approach to infer a three-dimensional displacement field at each side of the deck using available LOS directions and Equation (7) presented in the methodology section. The MCMC approach indicated that the northwest side of deck 9 moved (starting February 2015) at a rate of about 3.5 ± 2.1 mm/yr west, 14.1 ± 7.5 mm/yr south, and 4.6 ± 1.1 mm/yr down compared to the northeast side (i.e., REF2). The southwest side of deck 9 moved at a rate of about 7.1 ± 2 mm/yr west, 9.8 ± 6.1 mm/yr north, and 0.9 ± 1.1 mm/yr down compared to the northeast side, while the southeast side of deck 9 moved at a rate of about 7.5 ± 3.5 mm/yr west, 6.15 ± 6 mm/yr north, and 1.4 ± 1.1 mm/yr down compared to REF2. Figure 3 shows a three-dimensional characterization of the bridge deck and Table 3 indicates the numeric values inferred using the MCMC approach (Figure 2), suggesting that the deck 9 was characterized since 2015 by increased deformation accumulation on each side of the deck compared to deck 10 and deck 11. The trend inversion of the three-dimensional displacement also showed that on deck 9 the deformation accumulated mostly in the north–south direction, and decks 9–11 were relatively stable compared to the ground underneath their corresponding pier.

4. Discussion

This study highlights how the availability of new constellations of SAR could be applied to bridge deformation monitoring. While these techniques cannot decisively distinguish between stress accumulation or material degradation processes, they are useful to detect structural distress signs. A dedicated study of the available literature also confirms how these methods proved to be crucial in characterizing three-dimensional displacements with millimetric accuracy and capable of highlighting anomalous behaviors in near-real time [17,18]. The comprehensive multi-sensor dataset available for the Morandi bridge collapse allows us to make several unique considerations regarding the overall quality of the presented data, the future use of such data for a near-real time monitoring system and bridge construction guidelines.
Our InSAR single line of sight multi temporal analysis coupled with the three-dimensional MCMC inversion scheme show an increased accumulation of relative displacements of deck 9 compared to decks 10 and 11, starting from 2015. The analysis does not show macroscopic deformations characterizing the base of the piers 9–11. The absence of such deformations is also confirmed by the IG18 report by the Commission of the Ministry of Infrastructures and Transportation [29]. Furthermore, there are several PS recorded from different dataset and localized on different areas of deck 9 showing an increased acceleration. These PS are located close to the strands on the north and south west side of deck 9 respectively (Figure 1, Figure 5, Movie 1 and Movie 3) and show an anomalous behavior consisting in an abrupt acceleration starting 12th March 2017 at the junction between the deck and the strand located on the northwest, southwest and east sides of deck 9 respectively.
At the time of writing, the causes of the Morandi bridge collapse have not been officially identified [29]. No consensus has been found yet on the failure mechanism and the conditions which led to that mechanism. Different potential failure mechanisms have been evaluated [30]. Possible collapse modes included: (a) The loss of capacity of one of the stays connecting the deck with the antenna of balanced system 9; (b) the shear failure of the cantilever beam of system 9 or the connecting simply supporting beam; (c) a local failure at the connection between the stay and the deck. These mechanisms have been considered in relation to their ability of triggering the complete failure of one of the stays and the subsequent failure of the whole system 9 [30]. With reference to the pre-existing conditions which could have led to one of these mechanisms, a general material degradation has been documented for both the stays and the deck beams [23,29]. Our results show deformation patterns compatible with both the accelerating degradation of a stay strand, leading to its elongation, and a deterioration of the deck–stay connection [31].
The IG18 analysis does not highlight any of these deformation trends, indicating that previous analyses of CSK satellite-based datasets have been reported to give inconclusive results on pre-collapse variations of the displacement field in time. The discrepancy between the two different results can be ascribable to several reasons such as different temporal coverage (i.e., the IG18 only covers the period 2016–2018 whereas our CSK analysis covers 8 years 2010–2018, different coherence threshold used for the data selection (our analysis only considers pixels with a temporal coherence higher than 0.6 whereas the coherence threshold for IG18 is not known), absence of line of sight diversity and the lack of SNT data, together with the three-dimensional inversion scheme used in this study.

5. Conclusions

The synergistic use of CSK, SNT and Envisat played a key role in our study. The SNT data have been acquired at a rate nearly three times faster than CSK (6 days repeat vs 16 days repeat for SNT and CSK respectively), while the large CSK orbital tube and better CSK spatial resolution (3 x 3 meters vs. 4 x 22 meters range and azimuth resolution for CSK and SNT respectively) allowed us to better geolocate the analyzed scatterers in the three-dimensional space by a factor of three, as confirmed by independent studies [32,33]. The SNT tight orbital tube introduces a further level of difficulty when geolocating the observed scatterers located on the northern/southern side of the bridge. In our study we used a reference ground control point for georeferencing our data and only considered scatterers along the bridge aligning with optical images. The points not aligning with the optical images are not located at the deck level and hence discarded from our analysis. However, the severe limitations introduced in urban area by a small orbital tube represent a straightforward advantage for geophysical applications in mountainous area characterized by severe slopes introducing geometrical decorrelations. While the Envisat historical data allowed disclosing the presence of any pre-existing long-term trend characterizing the bridge pillars, no data were available for the deck level, principally due to the spatial resolution and the repeat time of the Envisat satellite (25 x 5 meters ground range/azimuth resolution and an average of one acquisition every three months). The SNT dataset revealed to be fundamental for the three-dimensional deformation decomposition when merged with the CSK data. The absence of scatterers located at the deck level for CSK descending (on the southern side of the bridge) and SNT ascending leads to the suggestion of introducing artificial scatterers such as corner-reflectors or transponders embedded within the bridge design. These devices could bridge the gaps between the SNT and CSK constellations, allowing a precise characterization of infrastructures deformation and the capability to monitor specific areas of the bridge according to the system design.
Overall this study shows that either CSK and SNT InSAR measurements are now accurate enough to measure millimeter scale deformation processes, providing information, that can help mitigate deformation-related hazards to large-scale infrastructure.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2072-4292/11/12/1403/s1, Table S1. LOS trends from different satellites and line of sight recorded during the period 2015–2009. Table S2. Dataset characteristics used in the MCMC inversion. Figure S1. Envisat multi-temporal line of sight (LOS) analysis. Figure S2. Sentinel-1A/B (SNT) ascending multi-temporal Interferometric SAR analysis of the Morandi Bridge. Movie1-8. Time-series deformation for the ascending/descending line of sight (LOS) of the CSK/Sentinel satellites recorded on the north/south side of the bridge, respectively.

Author Contributions

P.M. wrote the paper, set up the methodology, processed and analyzed the CSK, SNT and Envisat data. G.G. helped P.M. writing the paper and provided the structural engineering analysis of the Morandi bridge. D.P. provided support with the SARPROZ software and the time-series analysis quality and interpretation. A.C. G.M. provided CSK data access and participated to time-series analysis data interpretation. C.T. provided access to the historical dataset and participated to time-series analysis data interpretation.

Funding

This research was funded by National Aeronautics and Space Administration.

Acknowledgments

The research was carried out at the Jet Propulsion Laboratory, California Institute of Technology, under a contract with the National Aeronautics and Space Administration. We thank the Italian Space Agency (ASI) for providing COSMO-SkyMed data for this project. Original COSMO-SkyMed product ASI Agenzia Spaziale Italiana (2009–2018). The Sentinel data were provided courtesy of the European Space Agency. Envisat time-series analysis have been provided by the Italian Ministry of the Environment through the Geoportale Nazionale database. (http://www.pcn.minambiente.it/mattm/).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Milillo, P.; Bürgmann, R.; Lundgren, P.; Salzer, J.; Perissin, D.; Fielding, E.; Biondi, F.; Milillo, G. Space geodetic monitoring of engineered structures: The ongoing destabilisation of the Mosul dam, Iraq. Sci. Rep. 2016, 6, 37408. [Google Scholar] [CrossRef]
  2. Cigna, F.; Lasaponara, R.; Masini, N.; Milillo, P.; Tapete, D. Persistent Scatterer Interferometry Processing of COSMO-SkyMed StripMap HIMAGE Time Series to Depict Deformation of the Historic Centre of Rome, Italy. Remote Sens. 2014, 6, 12593–12618. [Google Scholar] [CrossRef] [Green Version]
  3. Boni, R.; Bosino, A.; Meisina, C.; Novellino, A.; Bateson, L.; McCormack, H. A Methodology to Detect and Characterize Uplift Phenomena in Urban Areas Using Sentinel-1 Data. Remote Sens. 2018, 10, 607. [Google Scholar] [CrossRef]
  4. Chatterjee, R.; Fruneau, B.; Rudant, J.; Roy, P.; Frison, P.-L.; Lakhera, R.; Dadhwal, V.; Saha, R.; Dadhwal, V. Subsidence of Kolkata (Calcutta) City, India during the 1990s as observed from space by Differential Synthetic Aperture Radar Interferometry (D-InSAR) technique. Remote Sens. Environ. 2006, 102, 176–185. [Google Scholar] [CrossRef]
  5. Floris, M.; Fontana, A.; Tessari, G.; Mulè, M. Subsidence Zonation Through Satellite Interferometry in Coastal Plain Environments of NE Italy: A Possible Tool for Geological and Geomorphological Mapping in Urban Areas. Remote Sens. 2019, 11, 165. [Google Scholar] [CrossRef]
  6. Stramondo, S.; Bozzano, F.; Marra, F.; Wegmuller, U.; Cinti, F.; Moro, M.; Saroli, M. Subsidence induced by urbanisation in the city of Rome detected by advanced InSAR technique and geotechnical investigations. Remote Sens. Environ. 2008, 112, 3160–3172. [Google Scholar] [CrossRef]
  7. Bozzano, F.; Esposito, C.; Mazzanti, P.; Patti, M.; Scancella, S. Imaging Multi-Age Construction Settlement Behaviour by Advanced SAR Interferometry. Remote Sens. 2018, 10, 1137. [Google Scholar] [CrossRef]
  8. Ma, P.; Lin, H.; Lan, H.; Chen, F. Multi-dimensional SAR tomography for monitoring the deformation of newly built concrete buildings. ISPRS J. Photogramm. Remote Sens. 2015, 106, 118–128. [Google Scholar] [CrossRef]
  9. Milillo, P.; Giardina, G.; DeJong, M.J.; Perissin, D.; Milillo, G. Multi-Temporal InSAR Structural Damage Assessment: The London Crossrail Case Study. Remote Sens. 2018, 10, 287. [Google Scholar] [CrossRef]
  10. Giardina, G.; Milillo, P.; DeJong, M.; Perissin, D.; Milillo, G. Evaluation of InSAR monitoring data for post-tunnelling settlement damage assessment. Struct. Control Health Monit. 2019, 26, e2285. [Google Scholar] [CrossRef]
  11. Roccheggiani, M.; Piacentini, D.; Tirincanti, E.; Perissin, D.; Menichetti, M. Detection and Monitoring of Tunneling Induced Ground Movements Using Sentinel-1 SAR Interferometry. Remote Sens. 2019, 11, 639. [Google Scholar] [CrossRef]
  12. Rousset, B.; Jolivet, R.; Simons, M.; Lasserre, C.; Riel, B.; Milillo, P.; Çakir, Z.; Renard, F. An aseismic slip transient on the North Anatolian Fault. Geophys. Res. Lett. 2016, 43, 3254–3262. [Google Scholar] [CrossRef] [Green Version]
  13. Karimzadeh, S.; Matsuoka, M. A Weighted Overlay Method for Liquefaction-Related Urban Damage Detection: A Case Study of the 6 September 2018 Hokkaido Eastern Iburi Earthquake, Japan. Geosci. 2018, 8, 487. [Google Scholar] [CrossRef]
  14. Di Naccio, D.; Kastelic, V.; Carafa, M.M.C.; Esposito, C.; Milillo, P.; Di Lorenzo, C. Gravity Versus Tectonics: The Case of 2016 Amatrice and Norcia (Central Italy) Earthquakes Surface Coseismic Fractures. J. Geophys. Res. Earth Surf. 2019, 124, 994–1017. [Google Scholar] [CrossRef]
  15. Chaussard, E.; Milillo, P.; Bürgmann, R.; Perissin, D.; Fielding, E.; Baker, B. Remote sensing of ground deformation for monitoring groundwater management practices: Application to the Santa Clara Valley during the 2012–2015 California drough. J. Geophys. Res. Solid Earth 2017, 122, 8566–8582. [Google Scholar] [CrossRef]
  16. Confuorto, P.; Di Martire, D.; Infante, D.; Novellino, A.; Papa, R.; Calcaterra, D.; Ramondini, M. Monitoring of remedial works performance on landslide-affected areas through ground- and satellite-based techniques. Catena 2019, 178, 77–89. [Google Scholar] [CrossRef]
  17. Sousa, J.J.; Bastos, L. Multi-temporal SAR interferometry reveals acceleration of bridge sinking before collapse. Nat. Hazards Earth Syst. Sci. 2013, 13, 659–667. [Google Scholar] [CrossRef]
  18. Selvakumaran, S.; Plank, S.; Geiß, C.; Rossi, C.; Middleton, C. Remote monitoring to predict bridge scour failure using Interferometric Synthetic Aperture Radar (InSAR) stacking techniques. Int. J. Appl. Earth Obs. Geoinf. 2018, 73, 463–470. [Google Scholar] [CrossRef]
  19. Milillo, P.; Riel, B.; Minchew, B.; Yun, S.-H.; Simons, M.; Lundgren, P. On the Synergistic Use of SAR Constellations Data Exploitation for Earth Science and Natural Hazard Response. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2016, 9, 1095–1100. [Google Scholar] [CrossRef]
  20. Morandi, R. Concrete Cable-stayed Structures; L’Industria Italiana del Cemento: Rome, Italy, 1980; Volume 10, pp. 887–910. (In Italian) [Google Scholar]
  21. Morandi, R. Il viadotto sul Polcevera per l’autostrada Genova-Savona; L’Industria Italiana del Cemento: Rome, Italy, 1967; Volume 12, pp. 849–873. (In Italian) [Google Scholar]
  22. Gentile, C.; Cabrera, F.M.Y. Dynamic investigation of a repaired cable-stayed bridge. Earthq. Eng. Struct. Dyn. 1997, 26, 41–59. [Google Scholar] [CrossRef]
  23. Bazzucchi, F.; Restuccia, L.; Ferro, G. Considerations over the Italian road bridge infrastructure safety after the Polcevera viaduct collapse: Past errors and future perspectives. Frattura ed Integrità Strutturale 2018, 12, 400–421. [Google Scholar] [CrossRef]
  24. Malerba, P.G. Inspecting and repairing old bridges: Experiences and lessons. Struct. Infrastruct. Eng. 2014, 10, 443–470. [Google Scholar] [CrossRef]
  25. Ferretti, A.; Prati, C.; Rocca, F. Permanent scatterers in SAR interferometry. IEEE Trans. Geosci. Remote Sens. 2001, 39, 8–20. [Google Scholar] [CrossRef]
  26. Ferretti, A.; Prati, C.; Rocca, F. Nonlinear subsidence rate estimation isong permanent scatterers in differential SAR interferometry. IEEE Trans. Geosci. Remote Sens. 2000, 38, 2202–2212. [Google Scholar] [CrossRef]
  27. Costantini, M.; Ferretti, A.; Minati, F.; Falco, S.; Trillo, F.; Colombo, D.; Novali, F.; Malvarosa, F.; Mammone, C.; Vecchioli, F.; et al. Analysis of surface deformations over the whole Italian territory by interferometric processing of ERS, Envisat and COSMO-SkyMed radar data. Remote Sens. Environ. 2017, 202, 250–275. [Google Scholar] [CrossRef]
  28. Colesanti, C.; Ferretti, A.; Novali, F.; Prati, C.; Rocca, F. SAR monitoring of progressive and seasonal ground deformation using the permanent scatterers technique. IEEE Trans. Geosci. Remote Sens. 2003, 41, 1685–1701. [Google Scholar] [CrossRef]
  29. Mortellaro, A.; Ievolella, G.; Lombardo, F.; Nuti, C.; Nanzi, I. Relazione della Commissione Ispettiva Ministeriale crollo Viadotto Polcevera; Ministero delle Infrastrutture e dei Trasporti: Roma, Italy, 2018. Available online: http://www.mit.gov.it/comunicazione/news/ponte-crollo-ponte-morandi-commissione-ispettiva-genova/ponte-morandi-online-la (accessed on 19 May 2019).
  30. Calvi, G.M.; Moratti, M.; O’Reilly, G.J.; Scattarreggia, N.; Monteiro, R.; Malomo, D.; Calvi, P.M.; Pinho, R. Once upon a Time in Italy: The Tale of the Morandi Bridge. Struct. Eng. Int. 2018, 29, 198–217. [Google Scholar] [CrossRef]
  31. Wolff, M.; Starossek, U. Cable loss and progressive collapse in cable-stayed bridges. Bridg. Struct. 2009, 5, 17–28. [Google Scholar] [CrossRef]
  32. Costantini, M.; Malvarosa, F.; Minati, F.; Trillo, F.; Vecchioli, F. Complementarity of high-resolution COSMO-SkyMed and medium-resolution Sentinel-1 SAR interferometry: Quantitative analysis of real target displacement and 3D positioning measurement precision, and potential operational scenarios. In Proceedings of the 2017 IEEE International Geoscience and Remote Sensing Symposium (IGARSS) 2017, Fort Worth, TX, USA, 23–28 July 2017; pp. 4602–4605. [Google Scholar]
  33. Bovenga, F.; Belmonte, A.; Refice, A.; Pasquariello, G.; Nutricato, R.; Nitti, D.O.; Chiaradia, M.T. Performance Analysis of Satellite Missions for Multi-Temporal SAR Interferometry. Sensors 2018, 18, 1359. [Google Scholar] [CrossRef]
Figure 1. The Morandi bridge before collapse. (A) Two-dimensional structural scheme of the bridge. From West to East, the bridge includes two single piers (1 to 2), six V-shaped piers (3 to 8) and three independent balanced systems (9 to 11). Dimensions in meters. (B) Three-dimensional view of the bridge.
Figure 1. The Morandi bridge before collapse. (A) Two-dimensional structural scheme of the bridge. From West to East, the bridge includes two single piers (1 to 2), six V-shaped piers (3 to 8) and three independent balanced systems (9 to 11). Dimensions in meters. (B) Three-dimensional view of the bridge.
Remotesensing 11 01403 g001
Figure 2. Probability histograms for the Morandi bridge decks 9–11. Markov chain Monte Carlo (MCMC) histograms output inferring three-dimensional deformation fields. LOS displacement rates used as inputs are listed in Table 2.
Figure 2. Probability histograms for the Morandi bridge decks 9–11. Markov chain Monte Carlo (MCMC) histograms output inferring three-dimensional deformation fields. LOS displacement rates used as inputs are listed in Table 2.
Remotesensing 11 01403 g002
Figure 3. Three-dimensional displacement rate of deformation of the Morandi bridge between 2015 and 2018 (in mm/yr). Inferred value of deformation using the Markov chain Monte Carlo (MCMC) approach described in the methodology section. All displacement components refer to the coordinate system at the bottom right. North, West and Up axes are characterized by independent arbitrary scales for visualization purposes. Line of sight displacement rates used as inputs are listed in Table 2. MCMC outputs are listed in Table 3 and Figure 2.
Figure 3. Three-dimensional displacement rate of deformation of the Morandi bridge between 2015 and 2018 (in mm/yr). Inferred value of deformation using the Markov chain Monte Carlo (MCMC) approach described in the methodology section. All displacement components refer to the coordinate system at the bottom right. North, West and Up axes are characterized by independent arbitrary scales for visualization purposes. Line of sight displacement rates used as inputs are listed in Table 2. MCMC outputs are listed in Table 3 and Figure 2.
Remotesensing 11 01403 g003
Figure 4. COSMO-SkyMed (CSK) and Sentinel-1A/B (SNT) multi-temporal Interferometric SAR analysis of the Morandi Bridge. (A) CSK descending dataset: (top) persistent scatterers (PS) height (center) PS deck level cumulative LOS displacements spanning 2009–2018; (bottom) LOS displacements for points located on the ground spanning 2009–2018. (B) SNT descending dataset: (top) persistent scatterers height; (center) PS deck level cumulative LOS displacements spanning 2015–2018; (bottom) LOS displacements for points located on the ground spanning 2015–2018. (C) CSK ascending dataset: (top) persistent scatterers height; (center) PS deck level cumulative LOS displacements spanning 2011–2018; (bottom) LOS displacements for points located on the ground spanning 2011–2018. (D) Zoom out of the CSK descending dataset spanning 2009–2018: (top) persistent scatterers height; (bottom) Cumulative LOS displacements. (E) Zoom out of the SNT descending dataset spanning 2015–2018: (top) persistent scatterers height; (bottom) Cumulative LOS displacements. (F) Zoom out of the CSK ascending dataset spanning 2011–2018: (top) persistent scatterers height; (bottom) Cumulative LOS displacements. All PS height and cumulative displacements maps share the same scale as shown in panel B. Panels A, B and C are rotated and aligned with the bridge east–west extension. All panels are in georeferenced WGS84 coordinate systems.
Figure 4. COSMO-SkyMed (CSK) and Sentinel-1A/B (SNT) multi-temporal Interferometric SAR analysis of the Morandi Bridge. (A) CSK descending dataset: (top) persistent scatterers (PS) height (center) PS deck level cumulative LOS displacements spanning 2009–2018; (bottom) LOS displacements for points located on the ground spanning 2009–2018. (B) SNT descending dataset: (top) persistent scatterers height; (center) PS deck level cumulative LOS displacements spanning 2015–2018; (bottom) LOS displacements for points located on the ground spanning 2015–2018. (C) CSK ascending dataset: (top) persistent scatterers height; (center) PS deck level cumulative LOS displacements spanning 2011–2018; (bottom) LOS displacements for points located on the ground spanning 2011–2018. (D) Zoom out of the CSK descending dataset spanning 2009–2018: (top) persistent scatterers height; (bottom) Cumulative LOS displacements. (E) Zoom out of the SNT descending dataset spanning 2015–2018: (top) persistent scatterers height; (bottom) Cumulative LOS displacements. (F) Zoom out of the CSK ascending dataset spanning 2011–2018: (top) persistent scatterers height; (bottom) Cumulative LOS displacements. All PS height and cumulative displacements maps share the same scale as shown in panel B. Panels A, B and C are rotated and aligned with the bridge east–west extension. All panels are in georeferenced WGS84 coordinate systems.
Remotesensing 11 01403 g004
Figure 5. Time-series deformation for the ascending/descending line of sight (LOS) of the CSK/Sentinel-1A/B satellites. Animated time-series of bridge deformation are shown in Movies 1–8.
Figure 5. Time-series deformation for the ascending/descending line of sight (LOS) of the CSK/Sentinel-1A/B satellites. Animated time-series of bridge deformation are shown in Movies 1–8.
Remotesensing 11 01403 g005
Table 1. Characteristics of the COSMO-SkyMEd (CSK) and Sentinel-1A/B (SNT) dataset used in this study. Heading angles are calculated positive clockwise from north. Table S1 shows acquisition dates and perpendicular baselines of each dataset.
Table 1. Characteristics of the COSMO-SkyMEd (CSK) and Sentinel-1A/B (SNT) dataset used in this study. Heading angles are calculated positive clockwise from north. Table S1 shows acquisition dates and perpendicular baselines of each dataset.
DatasetIncidence Angle (Degree)Heading Angle (Degree)Average Acq. Sampling (Days) Resolution RNG (m) × AZ (m)
CSK Ascending33.98–10.6163 × 3
CSK Descending26.63–168.5163 × 3
SNT Descending33.90–169.265 × 20
Table 2. Line of sight (LOS) trends from different satellites and line of sight recorded during the period 2015–2018. Reference point for each pillar are displayed in Figure 4A. West/East are referred to points located at the deck level (located on the west/east side) near the strands while average refers to the average trend of the deck. “/” indicates no persistent scatterer (PS) coverage in the area. Systems 9–11 refer to the balanced system numbers as listed in Figure 1 and Figure 4. North, South and Ground indicates where the pixels are located on the northern side, southern side, or at the base of the bridge, respectively. 2018–2010 trends for CSK data are presented in Table S2.
Table 2. Line of sight (LOS) trends from different satellites and line of sight recorded during the period 2015–2018. Reference point for each pillar are displayed in Figure 4A. West/East are referred to points located at the deck level (located on the west/east side) near the strands while average refers to the average trend of the deck. “/” indicates no persistent scatterer (PS) coverage in the area. Systems 9–11 refer to the balanced system numbers as listed in Figure 1 and Figure 4. North, South and Ground indicates where the pixels are located on the northern side, southern side, or at the base of the bridge, respectively. 2018–2010 trends for CSK data are presented in Table S2.
2018–2015Sentinel DSC (mm/yr)CSK DSC (mm/yr)CSK ASC (mm/yr)
WESTAVG.EASTWESTAVG.EASTWESTAVG.EAST
System 9 North–4.8–2.90–5.3–3.70–0.5–0.60
System 10 North52.203.61.20–3.6–40
System 11 North–4.8–1.30–1.30.30–0.41.60
System 9 South–0.50.6–0.5///5.61.46.8
System 10 South//////–3–3.4–2.7
System 11 South0.4–0.70///0.1–0.20.4
System 9 Ground/–3//0.1//–0.6/
System 10 Ground/–0.6//–0.1//–3.3/
System 11 Ground/–1.5//2////
Table 3. Three-dimensional deformation trends from different satellites and line of sight recorded during the period 2018–2015. West, average and east are referred to the west and east is referred to the closest point located at the deck level near the strands, while average refers to the average displacement of the deck. “/” indicates no PS coverage in the area. Deck –11 refers to the balanced system numbers as listed in Figure 1. E,N,U refer to East, North and Up, E prefix indicates error. B indicates bottom. 2018–2010 trends for CSK data are presented in Table S2.
Table 3. Three-dimensional deformation trends from different satellites and line of sight recorded during the period 2018–2015. West, average and east are referred to the west and east is referred to the closest point located at the deck level near the strands, while average refers to the average displacement of the deck. “/” indicates no PS coverage in the area. Deck –11 refers to the balanced system numbers as listed in Figure 1. E,N,U refer to East, North and Up, E prefix indicates error. B indicates bottom. 2018–2010 trends for CSK data are presented in Table S2.
Trends (mm/yr)Accuracy(mm/yr)
ENUEEENEU
Deck 9 NW–3.5–14.1–4.62.17.51.1
Deck 9 SW–7.19.8–0.926.11
Deck 9 SE–7.56.15–1.43.561.1
Deck 9 B–1.2–5.6–225.31.5
Deck 10 NW–2.4–8.9–3.12.14.51.1
Deck 10 B2.54.8–1.462.16.51.1
Deck 11 NW–2.8–1.522.632.241.1
Deck 11 SW0.52–3.81–0.90.75.11.3
Deck 11 SE0.1–2.180.11.361.5

Share and Cite

MDPI and ACS Style

Milillo, P.; Giardina, G.; Perissin, D.; Milillo, G.; Coletta, A.; Terranova, C. Pre-Collapse Space Geodetic Observations of Critical Infrastructure: The Morandi Bridge, Genoa, Italy. Remote Sens. 2019, 11, 1403. https://0-doi-org.brum.beds.ac.uk/10.3390/rs11121403

AMA Style

Milillo P, Giardina G, Perissin D, Milillo G, Coletta A, Terranova C. Pre-Collapse Space Geodetic Observations of Critical Infrastructure: The Morandi Bridge, Genoa, Italy. Remote Sensing. 2019; 11(12):1403. https://0-doi-org.brum.beds.ac.uk/10.3390/rs11121403

Chicago/Turabian Style

Milillo, Pietro, Giorgia Giardina, Daniele Perissin, Giovanni Milillo, Alessandro Coletta, and Carlo Terranova. 2019. "Pre-Collapse Space Geodetic Observations of Critical Infrastructure: The Morandi Bridge, Genoa, Italy" Remote Sensing 11, no. 12: 1403. https://0-doi-org.brum.beds.ac.uk/10.3390/rs11121403

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

Article Metrics

Back to TopTop