ABSTRACT We present a detailed survey of the ongoing destabilization process of the Mosul dam. The dam is located on the Tigris river and is the biggest hydraulic structure in Iraq. From a

geological point of view the dam foundation is poor due to a site geology formed by alternating strata of highly soluble materials including gypsum, anhydrite, marl and limestone. Here we

present the first multi-sensor cumulative deformation map for the dam generated from space-based interferometric synthetic aperture radar measurements from the Italian constellation

COSMO-SkyMed and the European sensor Sentinel-1a over the period 2014–2016 that we compare to an older dataset spanning 2004–2010 acquired with the European Envisat satellite. We found that

geodetic datasets suggests that the volume dissolution rate remains constant when the equivalent volume of total concrete injected during re-grouting operations is included in the

calculations. SIMILAR CONTENT BEING VIEWED BY OTHERS ADVANCED ANALYSIS OF SATELLITE DATA REVEALS GROUND DEFORMATION PRECURSORS TO THE BRUMADINHO TAILINGS DAM COLLAPSE Article Open access 04

January 2021 PRE-COLLAPSE SPACEBORNE DEFORMATION MONITORING OF THE KAKHOVKA DAM, UKRAINE, FROM 2017 TO 2023 Article Open access 03 April 2024 INSAR DATA FOR DETECTION AND MODELLING OF

OVEREXPLOITATION-INDUCED SUBSIDENCE: APPLICATION IN THE INDUSTRIAL AREA OF PRATO (ITALY) Article Open access 02 August 2024 INTRODUCTION There is a growing concern that the Mosul dam, Iraq,

is at risk of failure and consequent catastrophic flooding that could affect about 1.5 million people living near the Tigris river. Maintenance grouting to close pathways opened by water

infiltrating the underlying soluble evaporite deposits has been applied since the construction of the structure in 19841. After the Iraq war, a reconstruction plan by the United States

government and the Iraq Ministry of Water Resources was to provide for both short- and long-term solutions to this problem1. A US$27 million reconstruction project to establish a more

effective grouting system was attempted by the US Army Corps of engineers in 2007 with modest results toward maintaining the structural health of the dam. The Iraq Ministry of Water

Resources decided to keep the dam’s water level at below 319 m above sea level, thus limiting the hydrostatic pressure on the dam foundation1. With the rise of the self proclaimed Islamic

State (IS) the dam was captured by IS in August 2014 for a short period of time. Maintenance and cement grouting operations have ceased since then and the spillways remain blocked, raising

concerns about a possible dam failure. Here we present the first comprehensive multi-sensor cumulative deformation map for the dam generated from space-based synthetic aperture radar (SAR)

measurements (Fig. 1), which reveals that parts of the dam are undergoing rapid subsidence. Deformation was rapid during 2004–2010, it slowed in 2012–2014, and has increased since August

2014 when grouting operations stopped. We used data from multiple SAR satellites2 to measure the deformation at the dam in great spatial and temporal detail and shed new light on the

dynamics of the ongoing destabilization. RESULTS We process 62 images acquired in an ascending track geometry from the Italian Space Agency (ASI) COSMO-SkyMed (CSK) constellation of four

satellites to image ground deformation during 2012–2015 and 32 images acquired in a descending geometry by the European Space Agency (ESA) Sentinel-1A sensor starting October 2014 and

spanning 18 months. Six Interferograms from ESA’s Envisat satellite, from both ascending and descending geometries, were computed for monitoring the destabilization process during the period

2004–2010. A multi-temporal interferometric processing technique is applied to the CSK and Sentinel-1A phase measurements to compute two separate time-series of deformation exploiting

point-like scatterers3,4. We estimate time series of the water level at the dam wall using SAR amplitude images to measure the horizontal distance between the water and the top of the dam

and calibrate the results with satellite altimetry data. Data from the different spacecraft show that Mosul dam subsided throughout the observation period at varying rates and reveal a shift

in the location of the maximum deformation magnitude. The early Envisat dataset provides a good constraint on the first-order deformation pattern (Fig. 1a,b, Figs 2 and 3a), showing

subsidence with only small horizontal motions (Fig. 2). Figure 1c and d show the recent displacement rate along the LOS direction of the CSK and Sentinel-1A radar respectively over 8 months

(December 2014- July 2015). The observed subsidence is likely caused by dissolution of the dam substrate by water from the reservoir penetrating into “soluble sediments consisting of

eastward dipping beds and pods of anhydrite/gypsum, interbedded with layers of sand, silt and soft weathered marls, brecciated and clay-like marls, and fresh-to-weathered limestone”5. The

Envisat measurements from 2004 to 2010 indicate a maximum line-of-sight (LOS) velocity of −11.5 mm/year ± 6 mm/year (−12.47 mm/year vertical, negative rates indicate LOS increase or

subsidence, Fig. 3a). The fine spatial resolution (3 m) and dense temporal sampling of the CSK constellation allow us to better constrain the transient deformation during 2012–2015. A

regression analysis suggests that the dam is deforming at a maximum LOS rate of −12.1 ± 1.8 mm/year (−14.9 mm/year if considered vertical displacement) (Fig. 3b). Subsidence rates from

December 2012 to December 2013 did not exceed −2.3 ± 1.8 mm/yr (Fig. 4a). Due to the polar orbit geometry of all the SAR satellites, we do not have constraints on N-S displacements that

might be associated with fluctuating water levels and associated pressure changes on the dam. However, the CSK LOS time-series do not correlate with SAR-based water level measurements,

suggesting only a modest contribution by elastic dam deformation due to the changing water load. The most recent Sentinel-1A data confirm the subsidence rate observed with CSK (−14.6 ± 4.3 

mm/year) (Figs 1d and 4a) extending our observations until March 2016. The higher standard deviation is mainly due to the reduced number of scenes (25) increasing the noise level in the

non-linear time-series estimation (see methods section for major details). Using ascending and descending viewing geometries we can infer the horizontal (east-west component) and vertical

displacement rates using images acquired during the same time interval (Fig. 2). We oversample the results on a common 3 m grid of the CSK Stripmap-himage products. The results confirm that

the vertical velocities greatly exceed the horizontal contribution (Figs 2 and 3). Comparison of the 2004–2010 and 2014–2015 subsidence-rate profiles indicates a 300 m eastward shift of the

peak subsidence toward the dam’s main spillway (Fig. 3). The main difference between the spatial patterns of deformation during these two periods is related to the widening of the

CSK/Sentinel-1a vertical profile (Fig. 3b). The Envisat subsidence rate profile (Fig. 3a) is 5 mm/year slower compared to the CSK/Sentinel-1A profile. DISCUSSION The Mosul dam is subsiding

at a linear rate of ~15 mm/year compared to 12.5 mm/year subsidence rate in 2004–2010. The subsidence restarted at the end of 2013 after re-grouting operations stopped. The causes of the

observed linear subsidence process can be found in the human activities that have promoted the evaporite–subsidence development, primarily in gypsum deposits and may enable, in case of

continuous regrouting stop, unsaturated water to flow through or against evaporites deposits, allowing the development of small to large dissolution cavities. Large vertical movements have

resulted from the dissolution of extensive gypsum strata previously mapped beneath the Mosul dam5. Increased subsidence rate has been due to the absence of regrouting underlying the dam

basement. The subsidence seems to have a linear behavior but we can not exclude a future acceleration due to the increased dissolution speed. Although a detailed continuum mechanical

modeling of the observed subsidence is not the primary goal of this paper, the simplified model (supplementary information) does show the ability of the observed deformation to provide

first-order constraints on the primary processes involved in the dam’s basement. In Particular, simple models of the pattern of deformation put a rough constraint on the source depth,

assuming a horizontal tensile dislocation in an elastic half-space6) (See auxiliary material). Even if the adopted model does not take into consideration the plasticity conditions of the

rock succession, it is consistent with _in situ_ measurements of water temperature made by the US military using independent measurements5. The slowdown of the deformation can be interpreted

as being due to the reduced water level during the December 2012 to June 2013 dry season (Fig. 3b) and consequent low water pressure at the dam basement, which we suggest results in reduced

dissolution rates of the evaporites underlying the dam foundation. The inferred lateral propagation of the deformation source seems to be compatible with a hydrogeologic model5 constrained

by ground measurements in which the dam substrate is characterized by an increased hydraulic gradient toward the east abutment: promoting fluid pathways connecting the reservoir to the

downstream rocks in the subsurface. The model identifies the reservoir-induced pressure as the main cause of an increase in the dissolution rates. Moreover the abundant presence of readily

soluble gypsum below the eastern dam foundation close to the main spillway seems compatible with the observed subsidence shift and the eastward deepening of the dissolution front. These

observations agree with the previously mapped geometry of the sedimentary strata5, characterized by eastward dipping gypsum beds. SAR data analysis together with the inferred lateral

propagation from knowledge of the geology and hydraulic connectivity suggest that the locus of foundation dissolution has shifted eastward toward the dam’s principal spillway. Model

Inversions for the locus of volume loss in 2004–2010 infer a depth of 70 meters below the surface and a deepening in 2014–2015 of 160 meters below the surface (see auxiliary material). The

annual increase of water levels in the future, especially during the summer months, could lead to an increase of water pressure in the dam foundation, which combined with the absence of

regrouting could speed up the dissolution of the dam substrate and promote its destabilization. The Markov chain Monte Carlo (MCMC) approach5 enables us to calculate the dissolved volume

loss rates at different times (34.1 · 102 ± 9 m3/year during 2004–2010 vs. 44.6 · 102 ± 9 m3/year in 2014–2015), indicating a speedup of the dissolution that was likely caused by the lack of

maintenance and regrouting. To support this hypothesis, we calculate the sum of the injected cement volume rate1 and dissolved volume-loss rate in 2004–2010 (15 · 102 m3/year) and compare

it to the total dissolved volume in 2014–2015, when regrouting was stopped. We find that the estimates of total dissolution rates for the two observation periods agree within 10% (See

auxiliary material). This volume change analysis supports the hypothesis that a roughly steady dissolution rate holds on a yearly basis. However, the joint time-series analysis of water

levels and deformation on a weekly timescale points to the reduced water-level at the dam wall as a possible cause for the subsidence slow down during the period spanning December

2012–August 2014. We cannot confirm this short-term dynamic effect during the period 2004–2010 because of the lack of sufficiently frequent Envisat acquisitions. We can summarize our

conclusions in the following statements: * a This study highlights how the availability of new constellations of SAR sensors supported by historical SAR databases are enabling dam

deformation monitoring in which InSAR measurements are accurate enough to measure sub-kilometer scale geologic and anthropogenic processes, providing critical information fundamental for

mitigating anthropogenic hazards. * b Our observational strategy together with simple elastic models is in agreement with independent studies confirming both the depth of the past

dissolution source and the eastward shift of the dissolution front that occurred after the regrouting stopped. METHODS We used an MT-InSAR approach on the COSMO-SkyMed stripmap and Senintel

TOPS data based on a persistent scatterers (PS) interferometry processing chain. This stochastic approach3 has been widely adopted for monitoring deformation at dams7,8,9 when a copious

number of acquisitions (>20) is available. The ENVISAT dataset is characterized by only 6 interfrograms, hence the PS technique is not applicable and a standard DInSAR stacking technique

has been adopted (i.e. Unwrapped interferograms have been averaged and weighted with their temporal baseline). The PS technique assumes a set of N + 1 coregistered SLC SAR images. We

generate N single look differential interferograms with respect to a single master image. The master image is chosen maximizing the expected coherence of the generated interferograms

following: Where , T, _F__D_ and , _T__C_, _F__DC_ are the perpendicular, temporal and doppler baselines and their critical values respectively. A first rough identification of PS is

performed considering all the pixel with an inverse amplitude stability of 0.8 and choosing a pixel R relatively far from the dam as a reference point, the interferometric phase difference

between two neighbors pixels _x__r_ and _x_ _s_ in the interferogram i can be written as the sum of 5 terms: Where _ΔC__topo_ and _ΔC__V_ are proportional to the DEM errors and mean velocity

difference between pixels _x__r_ and _x_ _s_ respectively. _Δw__r,s_ contains the sum of three contributions related to atmospheric phase screen (APS) noise and non linear motion. We assume

for neighbors pixels. We estimate _ΔC__topo_, _ΔC__V_ maximizing the absolute value of the temporal coherence for neighbors pixels. Where: We then unwrap every interferogram estimating the

the phase term . We consider the temporal low pass and high pass component of: the low pass component is considered the non linear component of motion while for each interferogram a low pass

spatial filter is applied on the high pass component to estimate the APS. These steps have been performed considering only neighbor’s pixels. The same processing is then applied considering

only differences between the reference point R and all persistent scatterers including the estimated APS in eq. (4). The accuracy of the estimated velocities and DEM errors can be derived

by a regression analysis taking into account the irregularly sampled data (perpendicular baseline and temporal baseline in the case of the CSK constellation. As in ref. 4: where R is the

target-sensor distance, θ is the incidence angle, λ is the sensor wavelength, is the phase noise variance supposed independent of the acquisition and , are the variance of the perpendicular

and temporal baselines respectively. MT-InSAR provides ground displacements along the LOS. Given the geometry diversification our dataset it is possible to combine ascending and descending

cumulative displacement maps covering the same temporal period to highlight vertical and east-west deformations where: Where 39.2, 35.5 and 22.7 degrees are the incidence angles for

Sentinel, CSK and ENVISAT respectively. In this formulation we considered the C-band and X-band sensors penetration depth negligible when combining ascending and descending data. For what

concern the water level estimation we manually measure in a GIS environment the distance between the top of the dam and the water level in the SAR images. Our measurements measure the

pixels’ distance and estimate a scaling (slope) factor by calibrating the amplitude based water levels with the JASON-1 altimetrer. The standard errors are calculated taking into account the

pixel size of the sensors (Range + Azimuth). ADDITIONAL INFORMATION HOW TO CITE THIS ARTICLE: Milillo, P. _et al_. Space geodetic monitoring of engineered structures: The ongoing

destabilization of the Mosul dam, Iraq. _Sci. Rep._ 6, 37408; doi: 10.1038/srep37408 (2016). PUBLISHER'S NOTE: Springer Nature remains neutral with regard to jurisdictional claims in

published maps and institutional affiliations. REFERENCES * Bowen, S. W. Office of the special Inspector General for Iraq Reconstruction. Relief and reconstruction funded work at Mosul dam

Mosul, Iraq. (date of access: 16/03/2016). (http://cybercemetery.unt.edu/archive/sigir/20131001121159/http:/www.sigir.mil/files/assessments/PA-07-105.pdf) (2007). * Milillo, P. et al. On the

Synergistic Use of SAR Constellations’ Data Exploitation for Earth Science and Natural Hazard Response. IEEE JSTARS, 9(3), 1095–1100 (2016). ADS  Google Scholar  * Ferretti, A., Prati, C.

& Rocca, F. Nonlinear subsidence rate estimation using permanent scatterers in differential SAR interferometry. IEEE trans. Geo. Rem. Sens. 38(5), 2202–2212 (2000). Article  ADS  Google

Scholar  * Ferretti, A., Prati, C. & Rocca, F. (2001). Permanent scatterers in SAR interferometry. IEEE Transactions on geoscience and remote sensing, 39(1), 8−20.Okada, Y. Surface

deformation due to shear and tensile faults in a half-space. Bulletin of the seismological society of America, 75(4), 1135–1154 (1985). Google Scholar  * Kelley, J. R. et al. Geologic

setting of Mosul Dam and its engineering implications. (date of access: 16/03/2016). (http://el.erdc.usace.army.mil/elpubs/pdf/tr07-10.pdf) (2007) * Lundgren, P., Samsonov, S. V., López

Velez, C. M. & Ordoñez, M. Deep source model for Nevado del Ruiz Volcano, Colombia, constrained by interferometric synthetic aperture radar observations. Geophysical Research Letters,

42(12), 4816–4823 (2015). Article  ADS  Google Scholar  * Wang, T., Perissin, D., Rocca, F. & Liao, M. S. Three Gorges Dam stability monitoring with time-series InSAR image analysis.

Science China Earth Sciences, 54(5), 720–732 (2011). Article  Google Scholar  * Di Martire, D. et al. Comparison between differential SAR interferometry and ground measurements data in the

displacement monitoring of the earth-dam of Conza della Campania (Italy). Remote sensing of environment, 148, 58–69 (2014). Article  ADS  Google Scholar  * Milillo, P. et al. Monitoring dam

structural health from space: Insights from novel InSAR techniques and multi-parametric modeling applied to the Pertusillo dam Basilicata, Italy. International Journal of Applied Earth

Observation and Geoinformation, 52, 221–229. (2016). Article  ADS  Google Scholar  Download references ACKNOWLEDGEMENTS We thank the Italian Space Agency (ASI) for providing COSMO-SkyMed

data for this project. Original COSMO-SkyMed product ASI Agenzia Spaziale Italiana (2012–2016). The Envisat-ASAR and Sentinel data were provided courtesy of the European Space Agency. Mosul

Jason Altimeter data provided by the United States department of agriculture, foreign agricultural service. The The work of P. Milillo was sponsored by the National Aeronautics and Space

Administration Postdoctoral Program. The research was carried out at the Jet Propulsion Laboratory, California Institute of Technology, under a contract with the National Aeronautics and

Space Administration. AUTHOR INFORMATION AUTHORS AND AFFILIATIONS * NASA Jet Propulsion Laboratory, California Institute of Technology, 4800 Oak Grove Drive, Pasadena, 91109, (CA), USA

Pietro Milillo, Paul Lundgren & Eric Fielding * Dept. of Earth and Planetary Science, University of California, Berkeley, 389 McCone Hall, Berkeley, 94720, (CA), USA Roland Bürgmann *

GFZ German Research Centre for Geosciences, Physics of Earthquakes and Volcanoes, Telegrafenberg, 14473, Potsdam, Germany Jacqueline Salzer * Lyles School of Civil Engineering, Purdue

University, 550 Stadium Mall Drive, 47907, West Lafayette (IN), USA Daniele Perissin * University of L’Aquila, Italy Filippo Biondi * Italian Space Agency, Contrada Terlecchia, Matera,

75100, (MT), Italy Giovanni Milillo Authors * Pietro Milillo View author publications You can also search for this author inPubMed Google Scholar * Roland Bürgmann View author publications

You can also search for this author inPubMed Google Scholar * Paul Lundgren View author publications You can also search for this author inPubMed Google Scholar * Jacqueline Salzer View

author publications You can also search for this author inPubMed Google Scholar * Daniele Perissin View author publications You can also search for this author inPubMed Google Scholar * Eric

Fielding View author publications You can also search for this author inPubMed Google Scholar * Filippo Biondi View author publications You can also search for this author inPubMed Google

Scholar * Giovanni Milillo View author publications You can also search for this author inPubMed Google Scholar CONTRIBUTIONS P.M., wrote the main manuscript helped by R.B., P.L., D.P. and

J.S. P.M., R.B., P.L., J.S. and E.F. worked on the modeling and interpretation of the modeled results. J.S. focused the Envisat SAR data and prepared figs 2–3. G.M. and F.B. Processed the

COSMO-SkyMed data. D.P. provided the processing software and helped P.M., G.M. and F.B. Interpreting the InSAR time-series results. P.M. Processed the Sentinel-1A data and the altimeter

data. All authors reviewed the manuscript. ETHICS DECLARATIONS COMPETING INTERESTS The authors declare no competing financial interests. ELECTRONIC SUPPLEMENTARY MATERIAL SUPPLEMENTARY

INFORMATION RIGHTS AND PERMISSIONS This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included

in the article’s Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain

permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/ Reprints and permissions ABOUT THIS ARTICLE

CITE THIS ARTICLE Milillo, P., Bürgmann, R., Lundgren, P. _et al._ Space geodetic monitoring of engineered structures: The ongoing destabilization of the Mosul dam, Iraq. _Sci Rep_ 6, 37408

(2016). https://doi.org/10.1038/srep37408 Download citation * Received: 05 July 2016 * Accepted: 28 October 2016 * Published: 06 December 2016 * DOI: https://doi.org/10.1038/srep37408 SHARE

THIS ARTICLE Anyone you share the following link with will be able to read this content: Get shareable link Sorry, a shareable link is not currently available for this article. Copy to

clipboard Provided by the Springer Nature SharedIt content-sharing initiative