Next Article in Journal
Marine Sediments from a Contaminated Site: Geotechnical Properties and Chemo-Mechanical Coupling Processes
Next Article in Special Issue
Catalogue of the Geological Effects of Earthquakes in Spain Based on the ESI-07 Macroseismic Scale: A New Database for Seismic Hazard Analysis
Previous Article in Journal
Ecosystem Services Assessment and Valuation of Atmospheric Magnesium Deposition
Previous Article in Special Issue
Paleoliquefaction Studies and the Evaluation of Seismic Hazard
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Sequential Earthquake Damage Assessment Incorporating Optimized sUAV Remote Sensing at Pescara del Tronto

1
Department of Civil and Environmental Engineering, Brigham Young University, Provo, UT 84602, USA
2
Department of Chemical Engineering, Brigham Young University, Provo, UT 84602, USA
3
Department of Architecture and Design (DAD), Politecnico di Torino, Viale Mattioli, 39, 10125 Torino (TO), Italy
*
Author to whom correspondence should be addressed.
Geosciences 2019, 9(8), 332; https://doi.org/10.3390/geosciences9080332
Submission received: 13 June 2019 / Revised: 9 July 2019 / Accepted: 23 July 2019 / Published: 29 July 2019

Abstract

:
A sequence of large earthquakes in central Italy ranging in moment magnitudes (Mw) from 4.2 to 6.5 caused significant damage to many small towns in the area. After each earthquake in 2016 (24 August and 26 October), automated small unmanned aerial vehicles (sUAV) acquired valuable imagery data for post-hazard reconnaissance in the mountain village of Pescara del Tronto, and were applied to 3D reconstruction using Structure-from-Motion (SfM). In July 2018, the site was again monitored to obtain additional imagery data capturing changes since the last visit following the 30 October 2016 Earthquake. A genetic-based mission-planning algorithm that delivers optimal viewpoints and path planning was field tested and reduced the required photos for 3D reconstruction by 9.1%. The optimized 3D model provides a better understanding of the current conditions of the village, when compared to the nadir models, by containing fewer holes on angled surfaces, including an additional 17% surface area, and with a comparable ground-sampling distance (GSD) of ≈2.4 cm/px (≈1.5 cm/px when adjusted for camera pixel density). The resulting three time-lapse models provide valuable metrics for ground motion, progression of damage, resilience of the village, and the recovery progress over a span of two years.

1. Introduction

A sequence of large earthquakes caused considerable damage to multiple cities surrounding the faults in the Central Italian Apennine Mountain Range [1,2,3,4,5,6,7]. Analysis of high-magnitude earthquake sites assists post-disaster recovery and preventative efforts; an example in Iran is shown by Karimzadeh et al. [8]. Images of the affected area were collected by small unmanned aerial vehicles (sUAV) and 3D digital models of the small mountain village of Pescara del Tronto were made using Structure-from-Motion (SfM). Three models of Pescara del Tronto represent the post-disaster conditions after the August 2016 earthquake, after the October 2016 earthquake, and in July 2018, two years after the series of earthquakes; these models narrate the progression of damage and recovery. Volumetric measurements, movement detection, post-hazard reconnaissance, exterior safety inspections of buildings, documented recovery developments, and more can be gleaned from the sequential 3D reconstructed models. Pescara del Tronto presents a case study for comparing sequential changes of the city, comparing automated grid and optimized camera location methods, and an introduction to additional analysis (statistical and geotechnical).
The 2016 Italy earthquakes devastated communities along the Apennine Mountain Range along a series of normal faults with rake angles of 80 to 100 . The August 2016 earthquake caused significant damage in the following villages: Arquata del Tronto, Accumoli, Amatrice, and Pescara del Tronto [4,5]. The October events caused significant damage primarily in the villages of Visso, Ussita, and Norcia [2,6], and exacerbated the damage in Pescara del Tronto. Other analyses of the 2016 earthquake sequence include geohazard forensics of Norcia in Barone and Di Maggio [9] as well as numerical modeling of ground deformation in Valerio et al. [10], the mechanisms of the earthquakes in Zhong et al. [11], and Amatrice building damage assessment using synthetic aperture radar (SAR) in Karimzadeh and Matsuoka [12]. The earthquake events appear in a gap between two recent M6.1 events—the 1997 Umbria-Marche earthquake to the north-west and the 2009 L’Aquila earthquake to the south-east. Figure 1 from Stewart et al. [2] depicts the three major earthquakes in relation to the Umbria-Marche and L’Aquila earthquake events. This gap was determined to be an area of elevated risk prior to the 2016 earthquake events [4,5].
In total, a sequence of 17 major earthquake and aftershocks with moment magnitudes greater than 4.2 Mw were recorded in central Italy between August and November 2016. These events were initially recorded by the Italian National Seismic Network (Rete Sismica Nazionale, RSN) owned by the Italian Institute of Geophysics (Instituto Nazionale di Geofisica e Vulcanologia, INGV) near the villages of Amatrice, Accumoli, Norcia, and Visso [13,14]. Later, the number of recorded events increased to 49 by the United States Geological Survey (USGS) as recorded on the USGS web site (https://earthquake.usgs.gov/earthquakes/). Three main events identify the major earthquakes: M6.2 (24 August 2016), M6.1 (26 October 2016), and M6.6 (30 October 2016) [1,2,3,4,5] with the remainder classified as aftershocks.
Due to the temporal and geographical proximity of the series of earthquakes, the city of Pescara del Tronto was devastated by these events. The city lies within approximately 29 km of the all three main earthquake events and 8 km from the M6.2 24 August earthquake event (see Figure 2). The peak ground acceleration (PGA) in city for the M6.2 24 August, M6.1 26 October, and M6.6 30 October events were 0.48 g, 0.10 g, and 0.34 g, respectively [4,5]. The closest ground motion station to the events’ epicenters was located at 42.962 N 13.050 E at a distance of 29.6 km from Pescara del Tronto, as recorded on the USGS web site for each event (https://earthquake.usgs.gov/earthquakes/). The M6.2 24 August and M6.6 30 October events had larger PGA in Pescara del Tronto than the M6.1 26 October events due to proximity. The M6.2 24 August and M6.6 30 October events had the most devastating consequences on Pescara del Tronto, leaving the city uninhabitable until cleanup efforts conclude.
To the authors’ knowledge, the referenced Geotechnical Extreme Events Reconnaissance (GEER) reports follow international guidelines for post-disaster reconnaissance, but the scope of this paper is neither Rapid Mapping (hours to days post-disaster) nor Risk and Recovery Mapping (weeks to a couple months after the event) as described by Copernicus Emergency Management (https://emergency.copernicus.eu/) [1,2,3,4,5,15]. As the sequential data is as far out as 2 years post-disaster, this paper is an additional ad hoc representation of the earthquake recovery at Pescara del Tronto.

2. Method of Preliminary sUAV Reconnaissance and Modeling in 2016

The use of sUAVs in disaster reconnaissance prevents researchers and surveyors from having to enter many high-risk areas and reduces the time spent collecting surveyed data at every desired point. A recent example of high resolution sUAV data acquisition is Pavelka et al. [16] of the Pista Geoglyph in Peru, and expertise with sUAVs (and LiDAR) as shown in Pellicani et al. [17] demonstrates that sUAVs support and compliment geomorphic interpretations of disaster sites. Creating useful and accurate models with SfM software requires a series of photos at various normal and oblique angles that eventually tie into ground control points (GCPs) in a 3D reconstructed model. Three methods were used to collect the required imagery set of Pescara del Tronto during the two 2016 reconnaissance missions: basic nadir grids planned using the eMotion software to control the eBee system, multiple manual flights, and autonomous flights through the GSPro app.
During the 2016 reconnaissance trips, multiple sUAVs collected the images for reconstructing 3D models of Pescara del Tronto. Due to widespread use and an easily programmable software developer kit (SDK), the main method of collecting the images at Pescara del Tronto was with DJI products. During the 2016 reconnaissance trips, the platform used for collecting images was a DJI Phantom 4 ProTM quadrotor sUAV, equipped with a 4 K video camera with a 1/2.3” CMOS sensor, 94 field of view, 12.4 MP images, and a focal length set to infinity [3,18] was used. The Phantom 4 Pro system weighs 1.388 kg, has a maximum flight time of 28 min. In addition to multiple manual flights, GSPro was used as the autopilot to fly the drone and capture the desired nadir grid images.
The group from the Politecnico di Torino Geomatics provided a SenseFlyTM eBee. The eBee is a small fixed-wing sUAV manufactured by SenseFly at about 0.6 kg of platform weight, and can collect up to 8 km2 of surface area and 40 min of flight time depending on payload for the reconnaissance mission. The equipped payload included a digital camera Canon Power Shot S100, 1/1.7” Canon CMOS sensor, focal length of 5.2 mm, and 12 MP images [6]. While the platform can be navigated manually by a ground pilot, a piece of mission-planning software called eMotion [19] pre-programmed the system with GPS aerial waypoints. In addition to the models developed from the DJI Phantom 4 Pro, the eBee’s high altitude images provide an expanded view of the site and the surrounding terrain.
To georeference the models, GCPs were placed throughout the city and matched to multiple photos within the SfM software. Pescara del Tronto contained a total of 23 points with five 50 cm2 fabricated markers and 18 natural GCPs. The markers were surveyed by a Real Time Kinematic (RTK) configuration using the GNSS Geomax Zenith 35TM (www.geomax.com) [3]. Images showing the layout of the GCP markers and layout in Pescara del Tronto are shown below (see Figure 3 and Figure 4 from the GEER report) [1].
After the images and ground control points were obtained during the 2016 reconnaissance projects, 3D models developed and recorded current conditions and detected deformations that took place between the two events [6] (see Figure 5 and Figure 6 for the models). The models were made through BentleyTM ContextCapture modeling software which determines 3D point locations within a scene using GPS and camera metadata. The data is stored within the photos and used to calculate interspatial intersections—molding/stitching the photos together by identifying matching figures in the overlapping photos as in Koenderink and van Doorn [20]. SfM assigns point x, y, or z positions to similar features in each picture and spatially stitches each point together into a cohesive object [21]. However, the quality of the final product is generally determined by the skill and experience of the pilot constructing a grid overlap and obtaining additional images through manual flights.

3. Explanation of the Algorithms Implemented in 2018

Due to the mountainous terrain surrounding Pescara del Tronto, standard grid pattern photo gathering is inefficient and makes the collection of the required images very difficult [22]. Many flight planning applications have assisted autopilot software; however, flight path algorithms rarely consider occlusions, topographical features outside of elevation change, and do not change or optimize the six degrees of freedom data for camera positions in addition to the desired flight path. In post-hazard environments, prompt assessments and responses are essential. The time required to create a model is dependent upon the number of photos evaluated and number of tie points that are matched in a set of photos. Grid overlap flight paths not only take redundant photos, but require a pilot to estimate occluded areas, oblique positions, and take manual photos before generating an optimal model [20,23,24]. Moreover, locations similar to Pescara del Tronto have large hills, steep slopes, and cliffs that leave missing model sections with nadir grid flight paths. The nadir flight paths may go through a high unforeseen hill or may not collect the required images that are normal to a cliff or complex 3D surfaces. The optimized view and path planning algorithm developed at Brigham Young University (BYU) chooses optimal camera locations and oblique angles to capture the required images at various locations and orientations, thus, reducing the redundant photos and required time for 3D reconstruction. Genetic principles are incorporated into view-planning and field tested as part of the sequential earthquake reconnaissance case study. The basics of the discussed algorithms (assuming that the reader already has a general knowledge of combinatorial optimization) are followed by an explanation of the sUAV SfM work by Martin et al. [23] and is implemented at Pescara del Tronto.

3.1. Optimized View and Path Planning Algorithms

To overcome the downsides of the nadir grid approach and obtain better orthogonal views of Pescara del Tronto, a view-planning algorithm delineates optimized photo locations and angles through greedy heuristics [25]. View-planning identifies the best sensor locations for observing an object or site by planning oblique images for optimized SfM [26,27]. As shown in Martin et al. [23], when comparing models that use optimized view-planning with ones that use standard grid patterns, models using additional view-planning yield a greater accuracy of up to 43 % . These results led to the use of algorithmic view-planning (specifically model-based view-planning) at Pescara del Tronto.
The sUAV SfM used in Italy builds off previous work by Palmer et al. [28]. Various computational algorithms optimize what the sUAV needs to “see” as well as the best route to do so [23,29,30,31,32,33]. The path planning algorithm, follows a weighted Traveling Salesman Problem (TSP). Combining greedy and TSP algorithms effectively addresses non-deterministic polynomial-time hard (NP-hard) problems [34]. These computational tools organize the data points in a practical manner that produces photogrammetric results that look down vertically as well as at various oblique angles. The optimized approach is safer to people and equipment by taking into account the 3D terrain and produces more complete photogrammetric models. View and path planning algorithmic execution smooths out the sUAV mission and results (see Figure 7 for a basic graphic of the algorithms).
The view-planning algorithm follows the basic structure of a greedy algorithm [35]. While relatively simplistic, it quickly decides from a histogram where an sUAV can take a picture of the environment at the few locations that are worth taking a picture. The greedy algorithm follows the structure shown in the left of Figure 7 demonstrating a sub-optimal (highlighted by the red circles) but still effective solution by taking the smaller number in each binary option. While quick and generally effective, the downside to greedy heuristics is that no attempt to remember the previous choices is made, so a sub-optimal scenario could be a chosen result.
The path planning algorithm solves a weighted TSP (see the right portion of Figure 7) [36]. The Christofides algorithm is used to solve the TSP to a certain guaranteed degree of optimization as explained in Hoppe [37] and is used despite further research into computational solutions to the TSP as in Bożejko et al. [38]. The TSP path planning replaces the nadir grid approach in a manner that optimizes resources.

3.2. Evolutionary View-Planning

The optimized method applied at Pescara del Tronto deals with evolutionary model-based view-planning and closely follows Martin et al. [23] using a priori terrain data, a fitness function, inheritance, diversity and elitism, and convergence.

3.2.1. A Priori Terrain Data

A priori data is obtained from iterative flights and/or rough elevation models from open sources such as Google Earth and the USGS [39,40]. From the data, camera angles (azimuth and image elevation angles while neglecting roll (analogous to pitch, yaw, and roll with 1 less degree of freedom)), and X, Y, and Z coordinates are determined from an outlined polygon of the area of interest. The optimized approach directs the sUAV and obtains oblique angles by using the camera gimbal [41]. Figure 8 from Martin et al. [23] graphically shows a tiny data set (histogram) of possible camera locations and angles. At Pescara del Tronto, the similar histogram contains thousands of potential views with more noticeable camera angles.

3.2.2. Fitness Function

The histogram data is organized following fitness function f (see Equation (1)) as an objective where C is coverage (total visible points) and F c is functional coverage (a percentage that describes what the sUAV needs to “see” on its mission) [23]. For coverage C (see Equation (2)), v is visible points from the images, and T is the points describing the total sUAV mission area. For functional coverage F c (see Equation (3)), P is points that will be “seen” by the sUAV, and T is the same as in Equation (2). In other words: P is the points that are captured by at least 3 camera shots; F c is a normalized subset of C as P is divided by the points representing the entire terrain space T. The objective seeks to choose a minimum subset of the histogram data set, from the reduced histogram we generate sUAV camera locations and poses to maximize model quality (in accordance with similar procedures [42]).
f = C + F c
C = v T
F c = P T

3.2.3. Inheritance

To achieve the optimal sUAV camera locations and orientations, the genetic algorithm uses inheritance as in Martin et al. [23]. From the histogram, two random image solutions become mother and father, undertake a blend crossover with 0.7 probability, and the weighted average becomes the new image solution. Equations (4) and (5) further describe image crossover for each image variable ( 1 5 where each variable has 0.5 crossover probability) and r is randomly chosen between ( 0 , 1 ) . The inheritance continues until there are two new image solutions sets.
n e w i m a g e A e l e v a t i o n = r ( m o t h e r i m a g e e l e v a t i o n ) + ( 1 r ) ( f a t h e r i m a g e e l e v a t i o n )
n e w i m a g e B e l e v a t i o n = ( 1 r ) ( m o t h e r i m a g e e l e v a t i o n ) + ( r ) ( f a t h e r i m a g e e l e v a t i o n )

3.2.4. Diversity and Elitism

To ensure diversity in the population (image sets) simple random selection prevents premature solutions (that occur from a seeded tournament selection). A 5 % chance to mutate occurs during crossover on each variable (akin to a chromosome mutating [43]). A dynamic mutation parameter of α = 1 keeps initial values in variable bounds. An intermediate variable a accompanies the dynamic mutation (see Equation (6)) where g is current generation number and G is total generation number. Using the same example as in Martin et al. [23] using Equation (7), E is current image easting, E n e w is new image easting, E m a x is maximum allowed northing, E m i n is minimum allowed northing, and r e is a random number between ( E m i n , E m a x ) .
a = 1 g 1 G α
E n e w = E n e w = E m i n + ( r e E m i n ) a ( E E m i n ) 1 a if E E n e w = E m a x ( E m a x r e ) a ( E m a x E ) 1 a otherwise
The diversity of the image solution set is culled by elitism towards optimal results. Traits of better parents are weighted higher than lesser parents such that the better 25 % of solutions replace the lesser 25 % of solutions in each child solution set (e.g., an angle of a camera pose that views more points is more likely to propagate across generations). Preserving better potential solutions leads to convergence.

3.2.5. Convergence

Every generation of the algorithm is tested for convergence using Equation (8) from Martin et al. [23] (with appropriate parameter tuning). In the convergence equation, a c averages the current generation score and m is the mean of the average scores from the previous ten generations.
| a c m | a c < 0.00005
Upon convergence, the solution removes obsolete images that cannot “see” enough terrain points and grooms the image set by removing image locations that pass through or are under terrain. Cleaning the solution image set increases safety and avoids potential collisions that may have been generated. Once optimal, the matrix of remaining solutions is ready for the TSP as in Section 3.1, and, once completed, the data set is ready for a sUAV mission.

3.2.6. Why Evolutionary View-Planning Matters

A priori terrain data, fitness function, inheritance, diversity and elitism, and convergence produce a superior image set from genetic principles. The main downside of the optimized algorithms is that the use of evolutionary view-planning in the field is a bit slower than a manual approach (but is faster than a comparable grid approach). This is because of the required time at each camera location is longer due to stopping to stabilize while positioning the gimbal each time. However, even though the optimized approach requires fewer photos, the sUAV “sees” more than the grid approach that has the camera pointed in a single direction at all times. The optimized method captures pictures in areas that need extra detail rather than taking pictures at every location (which fetches excess data).

4. Application of the Algorithm to Pescara del Tronto in 2018

Pescara del Tronto serves as a prime location for testing the optimized path and view-planning algorithms. The 2018 visit to Pescara del Tronto yielded three models: one representing the data collected using an industry standard automated nadir flight planning application (a lawnmower approach with the camera pointing down); another model using the view-planning algorithm [44], and the final model created by combining all collected images, including those from manual flights, into a single model—the third model is not included in further quantitative discussion as the resolution is a skewed average of the other two models that use many more resources, but the combined model is used in Section 5 and Section 6 for qualitative comparisons. A sequential comparison of the 2016 and 2018 models is discussed in later sections of this paper.
During the July 2018 trip, an additional sUAV was used to collect images for reconstructing 3D models of Pescara del Tronto. A DJI Inspire 2, equipped with a Zenmuse X4S 4K video camera with a 1” CMOS sensor, 84 field of view, 20 MP images, and a focal length of 24 mm to collect the images for the nadir grid model. The Inspire 2 system weighs 3.44 kg, and has a maximum flight time of 27 min with the gimbal and Zenmuse X4S payload. Originally, the plan was for the Inspire 2 to fly the optimized view and flight path due to the Zenmuse X4S being significantly better than the Phantom’s camera. However, due to unforeseen events, the Optimized View and Flight Path algorithms were used on the Phantom 4 Pro with the same payload as the 2016 flights (all specifications from the cameras are kept the same except for pixel density).
Multiple ground control points were surveyed over the span of a week at Pescara del Tronto. During the first few days, images collected by using the nadir approach and manual flights were obtained in addition to 22 GCPs. Later in the week, the site was flown with the optimized view and path planning algorithms and was surveyed with 12 new GCPs (see Figure 9). The different number of GCPs was due to different groups performing the surveys because the GCPs were removed before the optimized view and path planning algorithms were used in the later part of the week.
The surface used to generate the cameras positions for Pescara del Tronto combines a Google Earth keyhole markup language (.kml) file and a point cloud (.ply) file of the area. The .ply file was obtained from the previous models of the area in 2016. By using the point cloud in addition to the .kml file, the view-planning algorithm identifies better camera locations for the flight compared to generated optimal locations based solely on the .kml file.
In 2018, the Phantom 4 Pro carried out the optimized mission and the DJI Inspire 2 flew the nadir grid mission. The Inspire 2 had a Zenmuse X4S camera attached that rates up to 4K like the Phantom 4 Pro camera but with 20 MP instead of 12 MP [45]. Figure 10 and Figure 11 demonstrate some of the camera locations for the nadir grid and optimized mission plans, respectively. Once the drone captured all the images at the desired locations, each image with the georeferenced metadata and additional GPS data was processed by Bentley’s ContextCapture 3D modeling software [46]. Table 1 and Table 2 summarize quantitative results using the GSPro autopilot (a typical industry nadir grid) and the optimized approach (targeted automated flight).
Table 1 reviews GCP, photo count, Surface Area, and ground-sampling distance (GSD). GSD is a measurement of model resolution derived from pixel density and camera distance from the target. The model created by the optimized view-planning algorithm, uses 30 % fewer photos per square kilometer while providing a more detailed model. The 3D model made from the optimized view locations had significantly fewer holes and a greater surface area by 17.1 % . Even though the nadir 3D model’s GSD is slightly better than the optimized equivalent, it is because of camera quality. This makes it difficult to compare the two models from the GSD. As the Inspire 2 is 167 % pixels compared to the Phantom 4 Pro, the adjusted GSD (as if both flights were done with the Inspire 2) based on MP differences would be 1.5 cm/px for the optimized approach and 2.0 cm/px for the nadir grid mission. The differences in camera qualities explain the discrepancy of the nadir-based results, and if projected would have higher average GSD than the optimized flight. However, even though both models have similar GSD, an analysis that does not include the GSD follows (see Table 2).
The quantitative comparison in Table 2 combines a four by five grid over Pescara del Tronto. At intersections of the grid, the photo count is found. Areas further up the slope have fewer photos in the nadir model, and areas of overlapping flights paths over the steeper slopes had a maximum photo count of 67. The analysis reveals that the model built from the optimized approach delivers a lower standard deviation of 6.78 when compared to the nadir approach that yields a standard deviation of 17.47 (258% the spread of the optimized approach). On average, each point delivers a mean photo count of 23.2 when optimized and a photo count of 28.13 from the nadir approach. These results correspond to the total photo counts of 716 and 912 respectively; this confirms that sloped regions of high interest and flux are better represented by the optimized approach.

5. Model Observations of Temporal Changes Two Years after the Earthquake

With the 2018 reconnaissance mission, the new techniques involved in collecting the images at desired locations with viewpoints greatly reduce the holes in the models and provide a more complete model of the city with its many slope failures near heavily forested regions (see Figure 12). The textured models, derived from the dense point-clouds, aid in post-reconnaissance visual inspections. The sequential models create “snapshots” in time that are revisited and compared with other “snapshots” to highlight patterns, differences, and continual new discoveries.
The 2016 and 2018 models of Pescara del Tronto allow changes in the environment over time to be observed. The optimized flight path for the sUAVs improves observations and makes further study of the area possible. The optimized approach analyzes drop offs, rockslides, and structures at high spatial resolution and is qualitatively clearer than the grid approach. Going to the site to collect data, fly the optimized path, and later create a model economizes time.
The models demonstrate the structural rebuilding process, geotechnical changes, and environmental regrowth that allows changes to be quantified (see major regrowth and environmental changes at Pescara del Tronto in Figure 12). Temporal changes have damaged infrastructure within the city. Reconstruction efforts focus mainly on stabilizing slopes and cleaning the rubble off roads that provide access to bridges. The following images focus on particular points of interest: specifically, the retaining wall and lower landslide (see Figure 13).
As recorded in the models, areas throughout the city are being reclaimed by vegetation. Areas where this is most evident are the upper slopes and lower retaining walls (see Figure 12). New trees and vegetation have begun to reclaim the damaged areas and provide natural protection and stabilization to the failed slopes (see Figure 13). Moreover, while rubble covers the majority of the 2016 models, the 2018 models record vast numbers of trees and vegetation that continue to reclaim most of the damaged areas. Because the soil is loose near the retaining wall, it is steep and thick vegetation has yet to develop. These small details are difficult to study when on site but are clearly recorded in the models—smaller details are clearer from optimized than grid reconstructions. For further research and questions, each location can easily be “revisited” as well as other places of interest that may have been missed while on site through means of the sequential models.
The central part of the city against a steep slope and retaining wall noticeably changed from the 2016 to 2018 models (see Figure 13). Over the span of two years, the steeper sloped area of the lower landslide (Section A) did not progress in its failure mode. Rather, in 2018, the slope was overgrown with foliage. The foliage is beneficial to the slope stability because the slope is reinforced by vines, bushes, and young trees. The main section of the slope failure (Section B), begins to be reclaimed by the foliage. However, with the consecutive quakes and further erosion, the slope continues to fail and slide onto the lower road. The third area of note is the large retaining wall failure (Section C). After the retaining wall fell during the October quakes, wind and rain continue to erode the shear sections of the cliff faster than foliage can reclaim and strengthen the slope. At the toe of the failure, foliage begins to grow on the rubble and likely will grow up the slopes in coming years (the optimized results reconstruct vegetation better than grid results). See http://prismweb.groups.et.byu.net/gallery2/ for access to the 3D models themselves [44].

6. Results of Observed Damage over Time

The value of these models for earthquake studies is greatly enhanced because of their sequential nature. From the three models of Pescara del Tronto, temporal “snapshots” are available between the August and October 2016 earthquakes, just after the October 2016 earthquake, and then 21 months after that same October earthquake (July 2018). This sequential data is rare, and allows for the opportunity to study changes over time in structures, slope failures, and other earthquake induced dangers and phenomena. The results of observed damage are similar to post-disaster situation maps and grading as used by Copernicus Emergency Management, but due to a longer scope, is an ad hoc addition for sequential earthquake recovery analysis at Pescara del Tronto [15]. The two older grid-based models provide needed data, but the third optimized model increases clarity at points of interest. Nadir-based models are used for the 2016 models and the combined nadir/optimized-based model is used for the 2018 model for qualitative comparisons. Due to the prohibitive size of the digital surface model (DSM) point-clouds without decimation, the primary focus of this section is qualitative visual analysis of the sequential damage over time (see Section 4 for quantitative comments. The study captures moments in time that give key visual information about the behavior of infrastructure developments and recovery after earthquakes.

6.1. Damage Observations between the August and October 2016 Events and the July 2018 Follow-Up Visit

Observations and data from this section are largely referenced from the GEER 2017 report of these events [1]. Significant damage was sustained in Pescara del Tronto after each of these earthquake events. The August events ( M a x P G A = 0.48 g) left nearly half of the structures collapsed. After the October event ( M a x P G A = 0.34 g) structures with damage levels of D4 or D5 increased from 65 % to 70 % .
Between the October event in 2016 and the 2018 July visit, only about 20 % of the shown building/structural areas (see Figure 14 and Figure 15) have been completely cleared away with about 9 % partially cleared. Given that a span of over 18 months had passed between the earthquake and the return visit, the time lapse shows how slowly disaster areas can recover.
For analysis and detailed inspection of the damage, Pescara del Tronto has been divided up the into seven focus areas (see Figure 16). The priority areas are grouped by similar damage characteristics. The most damaged areas lie within the dense city center; area one studies the structural failure of densely spaced houses. Areas two and three comprise the geotechnical failure of a retaining wall and failed slopes. Areas four and five are slope stability failures of nearby residences. Area six records the exposed pipeline near a gravel pit on the north side of the city. Area seven shows signs of a slope stability failure next to a bridge abutment.

6.1.1. Structural Damage Area (Location 1)

Location 1 (see Figure 16) is characterized by talus deposits with dense structural damage and rubble [1]. After the first earthquake in August 2016, this area still had many damaged or compromised structures standing and was scattered with rubble. After the October 2016 earthquake, almost no structures were left standing and the whole area was covered in dense rubble (see Figure 17). The rubble is often dense enough that individual building locations become hard to distinguish. Specific damage to structures is shown in Figure 14 and cleanup is graphically analyzed in Figure 15.
In the 2018 visit, much of the area was left as found following the October 2016 event; however, a portion was cleared out towards the north-west. Data from 2018 shows just under 70 % of Pescara del Tronto has yet to be cleared and recovered (see Figure 15). As previously stated, ≈20% of the area has been completely cleared and ≈9% has partially been cleared of debris. The process of clearing out the rubble and improvements for future reconstruction are documented two years since the cataclysmic series of earthquakes in 2016.

6.1.2. Landslides Impacting Strada Statale (Location 2)

A set of two large landslides above Strada Statale were observed after the August event in Area 2 (see Location 2 in Figure 16). The landslide as of September 2016 measured 71 m wide and 34 m tall (vertical measurement). Significant movement of this landslide is visible when comparing the September and December 2016 models (see Figure 18). Much of the road and retaining wall, above Strada Statale No. 4 of the northern landslide (Area 2a), which was still standing in the September model was observed to have fallen and collapsed in the December 2016 model. The collapse increased the landslide vertical height from 34 to 39 m. Large limestone boulders that are visible in the December 2016 and July 2018 models have dislodged and fallen from beneath the road at the top of the landslide. No significant changes occurred during the interim of the December 2016 and July 2018 models, though vegetation was beginning to cover much of the landslide area, and especially the area that did not fail between the September and December 2018 models.
The smaller landslide to the south (Area 2b) as of September 2016 was 18 m wide and 27 m tall. The December 2016 model shows increased width from 18 to 24 m while maintaining the same 27 m height. No significant changes were observed between the December 2016 and July 2018 models although overgrowth in the 2018 model prevents close observations; however, the reconstructed vegetation is crisper in the optimized model when compared to the grid-based models.

6.1.3. Failed Retaining Wall Surrounding Village (Location 3)

A 24 m tall and 85 m long masonry retaining wall supported a large portion of the village. This wall sustained partial failure after the August 2016 earthquake with two sections collapsing that were 5 m wide and respectively 13 m and 8 m tall. The retaining wall was destroyed after the October 2016 event leaving the valley full of debris. Though a few houses collapsed between the September 2016 and December 2016 models and a small area of debris was cleared above the retaining wall in the July 2018 model, a cloud comparison revealed no significant changes to the slope after reconnaissance in November 2016 (see Figure 19 and Figure 20 that were previously discussed in Figure 13). The angle of repose of the native talus behind the retaining wall was measured as 53 . Again, note that the sharp angles and vegetation of the optimized model are cleaner than previous grid-based models.

6.1.4. Localized Landslides along Gully Wall (Location 4)

Note Figure 21, the first two examples (a) and (b) appear as part of the GEER report in Stewart et al. [3] and the third (c) is from a combined nadir/optimized model generated July 2018. After the first earthquake, the gully wall depressed slightly; however, after the second earthquake in 2016, the 9 m thick and 20 m wide slope catastrophically fell leaving the 12 m scarp at an angle of nearly 52 , and the residential structure after August 2016 has crumbled into the gully as of November 2016 [3]. As of July 2018, the slope continues to erode and fail, but there have been some repairs to the road and slope. The stairs in the upper right of the shown portion for all three models, are more accurate in the optimized model than in the older models.

6.1.5. Landslide on Northern Slope of Village (Location 5)

Towards the northern end of Pescara del Tronto rests a shallow landslide in the desolated cityscape addressed in Stewart et al. [3] and is paraphrased and expounded upon (see Figure 22, the first two (a) and (b) appear as part of the GEER report in Stewart et al. [3] and the third (c) is from the newer 2018 model. The landslide first occurred as part of the August 2016 earthquake with lesser effects on infrastructure. The October 2016 earthquake did not increase the depth of the landslide; however, the width of the landslide nearly doubled. Note how different seasons effect the amount of vegetation in each photo, the November 2016 SfM allows for better visualization due to less vegetation than the Summer 2016 and Summer 2018 models. The July 2018 visual demonstrates the regrowth of vegetation (seasonal) as well as continued expansion of the landslide with build-up of structural debris down slope. The November 2016 and July 2018 sections are comparable despite the differences due to vegetation and the further spread of debris although Figure 22c reconstructs vegetation more clearly.

6.1.6. Haul Road and Exposed Pipeline near the Gravel Pit (Location 6)

As discussed in Stewart et al. [3] and shown in Figure 23 the pipelines off the haul road and near the gravel pit were affected by the earthquakes and over time ((a) and (b) are referenced in the GEER report). Additionally, Figure 23c shows the same location as of July 2018. From August 2016 to November 2016 about 8 m of the larger pipeline became exposed by weathering; the lower and smaller pipeline became visible in the same time frame. From November 2016 to July 2018 below the haul road, a hole in-between the larger pipeline and the haul road developed. The lower pipeline became more prominent with continued erosion. The pipelines sag and are more prominently exposed with the passage of time; the exposure is especially evident because the optimized approach forms a more detailed reconstruction of the pipes than the grid approach. The haul road has weathered overtime and has not been repaired since the slope of the hill degraded after the earthquake events until July 2018.

6.1.7. Large Slope Failure North of Village (Location 7)

Figure 24, sections (a) (August 2016), (b) (November 2016), and (c) (July 2018), depicts the large slope failure just north of the village. From August 2016 to November 2016 the landslide grew and devastated additional buildings. Although the landslide progressed, fixing the battered infrastructure and moving the displaced rocks and soils was not done. Looking at the same slope failure in July 2018, shows that the landslide progressed further, and the hillside eroded significantly. Much of the rubble in this area has been cleared away, and newly constructed terraced retaining walls help protect the area from future damage. The darker portions of the 2018 model are mats designed to hold the soil in place and prevent further damage. Looking at the copse of trees at the bottom middle of each model reveals that the optimized reconstruction better captures vegetation. Preventative efforts seek to mitigate further damage to the area.

7. Conclusions and Future Work

The sequence of earthquakes caused considerable damage throughout Pescara del Tronto. The widespread progression of damage in focus areas becomes difficult to track and manage. However, by the aid of sUAVs, digital imagery sensors collected data and provided crucial details after the August and October 2016 earthquakes. In addition, the July 2018 post-recovery reconnaissance sUAV mission, two years later, assessed the progression of damage. These images generate dense point-clouds and textured mesh models that are visually and analytically compared to obtain information about the construction of the city, geotechnical background of the city’s foundation, and the surrounding countryside.
During the 2018 reconnaissance trip, optimized view and path planning algorithms, increase the completeness of the 3D models by providing needed oblique angles (despite increased vegetation growth). When compared to the other model created from data obtained by a typical flight path developed by GSPro, optimized results (using genetic principles) cover 17 % more surface area, better represent objects, contain fewer holes on the slopes and sites of interest, requires 30 % fewer photos, has a higher resolution adjusted GSD of 1.5 cm/px, and yields an overall GSD of 2.4 cm/px, and is a case study of the applied algorithms for post-earthquake sUAV remote sensing. Optimized sUAV missions as a resource in post-hazard reconnaissance strengthens the ability to delineate and rate areas of focus affected by earthquakes. The constructed models provide a detailed reference source for the tracking of recovery efforts. Seven areas of interest were identified (see Figure 16). The focus areas depict two main causes of damage: geotechnical slope failures; and structural failure of old masonry buildings. Sequential changes to each area of interest are compared using models from August 2016, October 2016, and July 2018. A specifically promising avenue of future work is to compare complete DSMs of sequential changes with the quantitative tools available from CloudCompare. The progression of damage and recovery are assessed and show many areas being reclaimed by vegetation. While complete event details are difficult to analyze in person, event specifics and recovery are recorded in the sequential models. Further study of each location can easily be “revisited” (as well as other places of interest) and research of focus areas, throughout the models, promotes better understanding of the recovery that takes place after a series of earthquakes devastates a city.

Author Contributions

This paper was written, reviewed and edited by M.F., C.V., B.B., N.H., J.D. and J.P. with M.F. and N.H. as project administrators. M.F. led the data curation for the 2018 reconnaissance trip to Pescara del Tronto and the conceptualization of this research paper. B.B. created the models and performed the formal data analysis. M.F. validated the results of the analysis and the models. Visualization was accomplished by M.F., C.V., B.B. and N.H. Additional reviews and edits were done by B.H., J.H. (Joshua Hammond), J.J. and C.F., J.H. (John Hedengren) and K.F. provided funding acquisition and supervision for this work along with C.F.

Funding

This work has been funded by the Center for Unmanned Aircraft Systems (C-UAS), a National Science Foundation Industry/University Cooperative Research Center (I/UCRC) under NSF Award No. IIP-1650547 along with significant contributions from C-UAS industry members.

Acknowledgments

The authors would like to acknowledge the assistance of the Politecnico di Torino Geomatics Group, INGV, RSN, GEER, and Process Research and Intelligent Systems Modeling (PRISM) Group (Abe Martin) and the Research in Optimized Modeling (ROAM) Group. The collaborative efforts between the groups Brigham Young University and Politecnico di Torino Geomatic Group resulted in a wealth of data for creating detailed models of the affected areas in Pescara del Tronto. Models from the various sUAV missions in Italy can be viewed (see http://prismweb.groups.et.byu.net/gallery2/) [44].

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
BYUBrigham Young University
C-UASCenter for Unmanned Aircraft Systems
DSMDigital Surface Model
GCPGround Control Point
GSDGround-Sampling Distance
GEERGeotechnical Extreme Events Reconnaissance
GPSGlobal Positioning System
INGVInstituto Nazionale di Geofisica e Vulcanologia
I/UCRCIndustry/University Cooperative Research Center
PGAPeak Ground Acceleration
PRISMProcess Research and Intelligent Systems Modeling
NSFNational Science Foundation
ROAMResearch in Optimized Aerial Modeling
RMSRoot Mean Square
RSNRete Sismica Nazionale
RTKReal Time Kinematic
SDKSoftware Developer Kit
SfMStructure-from-Motion
sUASSmall Unmanned Aerial Systems
sUAVSmall Unmanned Aerial Vehicle
TSPTraveling Salesman Problem
USGSUnited States Geological Survey

References

  1. Stewart, J.P. Engineering Reconnaissance Following the 2016 M 6.0 Central Italy Earthquake: Ver 1; GeotechnicalExtreme Events Reconnaissance Association: Berkeley, CA, USA, 2016. [Google Scholar]
  2. Stewart, J.; Lanzo, G.; Ausilio, E.; Cairo, R.; Bozzoni, F.; Capatti, M.C.; della Pasqua, F.; Dezi, F.; Di Sarno, L.; Durante, M.; et al. Engineering Reconnaissance Following the October 2016 Central Italy Earthquakes; GeotechnicalExtreme Events Reconnaissance Association: Berkeley, CA, USA, 2017. [Google Scholar]
  3. Stewart, J.; Lanzo, G.; Ausilio, E.; Cairo, R.; Bozzoni, F.; Capatti, M.C.; della Pasqua, F.; Dezi, F.; Di Sarno, L.; Simonelli, A.L.; et al. Engineering Reconnaissance Following the October 2016 Central Italy Earthquakes—Version 2; GeotechnicalExtreme Events Reconnaissance Association: Berkeley, CA, USA, 2017. [Google Scholar]
  4. Stewart, J.; Lanzo, G.; Aversa, S.; Bozzoni, F.; Dashti, S.; Di Sarno, L.; Durante, M.; Simonelli, A.L.; Penna, A.; Foti, S.; et al. Engineering Reconnaissance of the 24 August 2016 Central Italy Earthquake. Version 2; GeotechnicalExtreme Events Reconnaissance Association: Berkeley, CA, USA, 2016. [Google Scholar]
  5. Stewart, J.; Zimmaro, P.; Lanzo, G.; Mazzoni, S.; Ausilio, E.; Aversa, S.; Bozzoni, F.; Cairo, R.; Chiara, M.; Castiglia, M.; et al. Reconnaissance of 2016 Central Italy Earthquake Sequence. Earthq. Spectra 2018, 34, 1547–1555. [Google Scholar] [CrossRef] [Green Version]
  6. Franke, K.W.; Lingwall, B.N.; Zimmaro, P.; Kayen, R.E.; Tommasi, P.; Chiabrando, F.; Santo, A. Phased Reconnaissance Approach to Documenting Landslides Following the 2016 Central Italy Earthquakes. Earthq. Spectra 2018, 34, 1693–1719. [Google Scholar] [CrossRef]
  7. Lanzo, G.; Tommasi, P.; Ausilio, E.; Aversa, S.; Bozzoni, F.; Cairo, R.; d’Onofrio, A.; Durante, M.G.; Foti, S.; Giallini, S.; et al. Reconnaissance of geotechnical aspects of the 2016 Central Italy earthquakes. Bull. Earthq. Eng. 2018, 1, 1–38. [Google Scholar] [CrossRef]
  8. Karimzadeh, S.; Feizizadeh, B.; Matsuoka, M. From a GIS-based hybrid site condition map to an earthquake damage assessment in Iran: Methods and trends. Int. J. Disaster Risk Reduct. 2017, 22, 23–36. [Google Scholar] [CrossRef]
  9. Barone, P.M.; Di Maggio, R.M. Forensic Investigations of Geohazards: The Norcia 2016 Earthquake. Geosciences 2018, 8, 316. [Google Scholar] [CrossRef]
  10. Valerio, E.; Tizzani, P.; Carminati, E.; Doglioni, C.; Pepe, S.; Petricca, P.; De Luca, C.; Bignami, C.; Solaro, G.; Castaldo, R.; et al. Ground Deformation and Source Geometry of the 30 October 2016 Mw 6.5 Norcia Earthquake (Central Italy) Investigated through Seismological Data, DInSAR Measurements, and Numerical Modelling. Remote Sens. 2018, 10, 1901. [Google Scholar] [CrossRef]
  11. Zhong, S.; Xu, C.; Yi, L.; Li, Y. Focal Mechanisms of the 2016 Central Italy Earthquake Sequence Inferred from High-Rate GPS and Broadband Seismic Waveforms. Remote Sens. 2018, 10, 512. [Google Scholar] [CrossRef]
  12. Karimzadeh, S.; Matsuoka, M. Building Damage Assessment Using Multisensor Dual-Polarized Synthetic Aperture Radar Data for the 2016 M 6.2 Amatrice Earthquake, Italy. Remote Sens. 2017, 9, 330. [Google Scholar] [CrossRef]
  13. Istituto Nazionale di Geofisica e Vulcanologia. INGV Centro Nazionale Terremoti: Seismic Network IV; Istituto Nazionale di Geofisica e Vulcanologia: Rome, Italy, 2017. [Google Scholar] [CrossRef]
  14. INGV Sede Irpinia. RETE Sismica Nazionale. 2019. Available online: http://www.gm.ingv.it/ (accessed on 15 March 2019).
  15. Copernicus Emergency Management Service. 2019. Available online: https://emergency.copernicus.eu/ (accessed on 2 April 2019).
  16. Pavelka, K.; Šedina, J.; Matoušková, E. High Resolution Drone Surveying of the Pista Geoglyph in Palpa, Peru. Geosciences 2018, 8, 479. [Google Scholar] [CrossRef]
  17. Pellicani, R.; Argentiero, I.; Manzari, P.; Spilotro, G.; Marzo, C.; Ermini, R.; Apollonio, C. UAV and Airborne LiDAR Data for Interpreting Kinematic Evolution of Landslide Movements: The Case Study of the Montescaglioso Landslide (Southern Italy). Geosciences 2019, 9, 248. [Google Scholar] [CrossRef]
  18. DJI. Phantom 4 Pro. 2019. Available online: https://www.dji.com/phantom-4-pro/info (accessed on 20 April 2019).
  19. Group, P. eMotion. 2019. Available online: https://www.sensefly.com/software/emotion/ (accessed on 20 April 2019).
  20. Koenderink, J.J.; van Doorn, A.J. Affine structure from motion. J. Opt. Soc. Am. A 1991, 8, 377. [Google Scholar] [CrossRef] [PubMed]
  21. Barrett, B.J. Field Validation of an Advanced Autonomous Method of Exterior Dam Inspection Using Unmanned Aerial Vehicles. Ph.D. Thesis, Brigham Young University, Provo, UT, USA, 2018. [Google Scholar]
  22. Scaioni, M.; Crippa, J.; Corti, M.; Barazzetti, L.; Fugazza, D.; Azzoni, R.; Cernuschi, M.; Diolaiuti, G.A. Technical Aspects Related to the Application of SFM Photogrammetry in High Mountain. ISPRS Int. Arch. Photogramm. Remote Sens. Spat. Inf. Sci. 2018, XLII-2, 1029–1036. [Google Scholar] [CrossRef]
  23. Martin, R.; Rojas, I.; Franke, K.; Hedengren, J. Evolutionary View Planning for Optimized UAV Terrain Modeling in a Simulated Environment. Remote Sens. 2015, 8, 26. [Google Scholar] [CrossRef]
  24. ASPRS. The Imaging & Geospatial Information Society; ASPRS: Bethesda, MD, USA, 2018. [Google Scholar]
  25. Jung, D.; Dong, Y.; Frisk, E.; Krysander, M.; Biswas, G. Sensor selection for fault diagnosis in uncertain systems. Int. J. Control 2018, 91, 1–11. [Google Scholar] [CrossRef]
  26. Henriques, M.; Roque, D. Unmanned Aerial Vehicles (UAV) as a Support to Visual Inspections of Concrete Dams; LNEC: Lisbon, Portugal, 2015. [Google Scholar]
  27. Zhao, J.; Cheung, S.C. Optimal visual sensor planning. In Proceedings of the 2009 IEEE International Symposium on Circuits and Systems, Taipei, Taiwan, 24–27 May 2009; pp. 165–168. [Google Scholar] [CrossRef]
  28. Palmer, L.M.; Franke, K.W.; Abraham Martin, R.; Sines, B.E.; Rollins, K.M.; Hedengren, J.D. Application and Accuracy of Structure from Motion Computer Vision Models with Full-Scale Geotechnical Field Tests. In IFCEE 2015; American Society of Civil Engineers: Reston, VA, USA, 2015; pp. 2432–2441. [Google Scholar] [CrossRef]
  29. Liu, J.; Sridharan, S.; Fookes, C. Recent Advances in Camera Planning for Large Area Surveillance. ACM Comput. Surv. 2016, 49, 6. [Google Scholar] [CrossRef]
  30. Rahimian, P.; Kearney, J.K. Optimal Camera Placement for Motion Capture Systems. IEEE Trans. Vis. Comput. Graph. 2016, 23, 1209–1221. [Google Scholar] [CrossRef] [PubMed]
  31. Arabnia, H.R.; Deligiannidis, L. Emerging Trends in Image Processing, Computer Vision and Pattern Recognition; Elsevier: Amsterdam, The Netherlands, 2015; pp. 1–609. [Google Scholar] [CrossRef]
  32. Martin, R.; Blackburn, L.; Pulsipher, J.; Franke, K.; Hedengren, J. Potential Benefits of Combining Anomaly Detection with View Planning for UAV Infrastructure Modeling. Remote Sens. 2017, 9, 434. [Google Scholar] [CrossRef]
  33. Lu, E.H.-C.; Yang, Y.W. A hybrid route planning approach for logistics with pickup and delivery. Expert Syst. Appl. 2019, 118, 482–492. [Google Scholar] [CrossRef]
  34. Tauer, G.; Date, K.; Nagi, R.; Sudit, M. An incremental graph-partitioning algorithm for entity resolution. Inf. Fusion 2019, 46, 171–183. [Google Scholar] [CrossRef]
  35. Cormen, T.H.; Leiserson, C.E.; Rivest, R.L.; Stein, C. Introduction to Algorithms, 3rd ed.; MIT Press: Cambridge, MA, USA, 2009; pp. 357–478. [Google Scholar]
  36. Gass, S.I.; Fu, M.C. (Eds.) Encyclopedia of Operations Research and Management Science: Traveling Salesman Problem; Springer: Boston, MA, USA, 2013. [Google Scholar] [CrossRef]
  37. Hoppe, C.; Wendel, A.; Zollmann, S.; Pirker, K.; Irschara, A.; Bischof, H.; Kluckner, S. Photogrammetric camera network design for micro aerial vehicles. In Proceedings of the Computer Vision Winter Workshop (CVWW), Mala Nedelja, Slovenia, 1–3 February 2012; Volume 8, pp. 1–3. [Google Scholar]
  38. Bożejko, W.; Gnatowski, A.; Niżyński, T.; Affenzeller, M.; Beham, A. Local Optima Networks in Solving Algorithm Selection Problem for TSP. In Proceedings of the DepCoS-RELCOMEX 2018, Brunow, Poland, 2–6 July 2018; pp. 83–93. [Google Scholar] [CrossRef]
  39. Dollison, R. The National Map: New Viewer, Services, and Data Download: U.S. Geological Survey Fact Sheet 2010–3055; United States Geological Survey: Reston, VA, USA, 2016. [CrossRef]
  40. Whitten, A. Decentralized Planning for Autonomous Agents Cooperating in Complex Missions. Master’s Thesis, Massachusetts Institute of Technology, Cambridge, MA, USA, 2010. [Google Scholar]
  41. Dobrokhodov, V.N.; Yakimenko, O.A.; Jones, K.D.; Kaminer, I.; Bourakov, E.; Kitsios, I.; Lizarraga, M. New generation of rapid flight test prototyping system for small unmanned air vehicles. In Proceedings of the AIAA Modeling and Simulation Technologies Conference, Hilton Head, SC, USA, 20–23 August 2007; pp. 1–26. [Google Scholar]
  42. Li, R.; Fulp, E.W. Evolutionary Approaches for Resilient Surveillance Management. In Proceedings of the 2017 IEEE Security and Privacy Workshops (SPW), San Jose, CA, USA, 25–25 May 2017. [Google Scholar] [CrossRef]
  43. Li, R.; Fulp, E. Using Evolutionary Approaches to Manage Surveillance Cameras in Dynamic Environments. In Proceedings of the 17th IEEE International Conference on Trust, Security and Privacy in Computing and Communications and 12th IEEE International Conference on Big Data Science and Engineering, New York, NY, USA, 1–3 August 2018; pp. 191–196. [Google Scholar] [CrossRef]
  44. PRISM, B. BYU PRISM 3D Model Gallery. 2019. Available online: http://prismweb.groups.et.byu.net/gallery2/ (accessed on 18 April 2019).
  45. DJI. Inspire 2. 2019. Available online: https://www.dji.com/inspire-2/info (accessed on 28 March 2019).
  46. Bentley Systems. Reality Modeling Software. 2019. Available online: https://www.bentley.com/ (accessed on 18 April 2019).
Figure 1. Map of central Italy showing moment tensors of major earthquakes since 1997 and the intermediate gap areas [2].
Figure 1. Map of central Italy showing moment tensors of major earthquakes since 1997 and the intermediate gap areas [2].
Geosciences 09 00332 g001
Figure 2. Main earthquake events in relation to Pescara del Tronto.
Figure 2. Main earthquake events in relation to Pescara del Tronto.
Geosciences 09 00332 g002
Figure 3. Fabricated GCP markers (a), surveyed natural formations (b) [1].
Figure 3. Fabricated GCP markers (a), surveyed natural formations (b) [1].
Geosciences 09 00332 g003
Figure 4. Locations of the 23 GCPs throughout Pescara del Tronto (red denote fabricated GCP markers, and yellow markers denote natural formation GCP markers [1].
Figure 4. Locations of the 23 GCPs throughout Pescara del Tronto (red denote fabricated GCP markers, and yellow markers denote natural formation GCP markers [1].
Geosciences 09 00332 g004
Figure 5. Model of Pescara del Tronto after the August 2016 Earthquake events.
Figure 5. Model of Pescara del Tronto after the August 2016 Earthquake events.
Geosciences 09 00332 g005
Figure 6. Model of Pescara del Tronto after the October 2016 Earthquake events.
Figure 6. Model of Pescara del Tronto after the October 2016 Earthquake events.
Geosciences 09 00332 g006
Figure 7. (Left) basic example of a greedy algorithm; (Right) demonstration of path planning.
Figure 7. (Left) basic example of a greedy algorithm; (Right) demonstration of path planning.
Geosciences 09 00332 g007
Figure 8. Graphical representation of possible camera locations and poses from a histogram [23].
Figure 8. Graphical representation of possible camera locations and poses from a histogram [23].
Geosciences 09 00332 g008
Figure 9. Locations of the 12 GCPs in Pescara Del Tronto used in the 2018 Optimized model (red represent surveyed fabricated GCP markers, green indicate natural GCP markers, and orange is the base station location).
Figure 9. Locations of the 12 GCPs in Pescara Del Tronto used in the 2018 Optimized model (red represent surveyed fabricated GCP markers, green indicate natural GCP markers, and orange is the base station location).
Geosciences 09 00332 g009
Figure 10. Industry standard mission plan: camera locations and angles represented as triangular prisms.
Figure 10. Industry standard mission plan: camera locations and angles represented as triangular prisms.
Geosciences 09 00332 g010
Figure 11. Optimized mission plan: camera locations and angles represented as triangular prisms.
Figure 11. Optimized mission plan: camera locations and angles represented as triangular prisms.
Geosciences 09 00332 g011
Figure 12. Model results of the temporal progress at Pescara del Tronto two years after the earthquake.
Figure 12. Model results of the temporal progress at Pescara del Tronto two years after the earthquake.
Geosciences 09 00332 g012
Figure 13. Lower retaining wall and slope failure progression over two years. Section A: Slope failure reclaimed by foliage; Section B: Further failure of the slope; Section C: Complete failure of retaining wall.
Figure 13. Lower retaining wall and slope failure progression over two years. Section A: Slope failure reclaimed by foliage; Section B: Further failure of the slope; Section C: Complete failure of retaining wall.
Geosciences 09 00332 g013
Figure 14. Map of building damage and cleared areas in Pescara Del Toronto following the October and August 2016 earthquakes and 2018 follow-up visit. (a) After the August 2016 earthquake; (b) After the October 2016 earthquake; (c) 2018 follow-up visit [1].
Figure 14. Map of building damage and cleared areas in Pescara Del Toronto following the October and August 2016 earthquakes and 2018 follow-up visit. (a) After the August 2016 earthquake; (b) After the October 2016 earthquake; (c) 2018 follow-up visit [1].
Geosciences 09 00332 g014
Figure 15. Bar graph showing the level and percentage of buildings damaged between the October and August 2016 earthquakes and the percentage of buildings cleared as of July 2018. Total building areas counted were 143. (a) After the August 2016 earthquake; (b) After the October 2016 earthquake; (c) 2018 follow-up visit. Modified from the GEER 2016 report [1].
Figure 15. Bar graph showing the level and percentage of buildings damaged between the October and August 2016 earthquakes and the percentage of buildings cleared as of July 2018. Total building areas counted were 143. (a) After the August 2016 earthquake; (b) After the October 2016 earthquake; (c) 2018 follow-up visit. Modified from the GEER 2016 report [1].
Geosciences 09 00332 g015
Figure 16. Pescara del Tronto Locality Map (similar to GEER map) [1]. Areas of Focus: 1—Structural Damage of Residential Area, 2—Landslides Impacting the Strada Statale, 3—Failed Retaining Wall Surrounding the Village, 4—Localized Landslides along the gully wall, 5—Landslide on Northern Slope of the Village, 6—Haul Road and Exposed Pipeline, 7—Large Slope Failure North of the Village.
Figure 16. Pescara del Tronto Locality Map (similar to GEER map) [1]. Areas of Focus: 1—Structural Damage of Residential Area, 2—Landslides Impacting the Strada Statale, 3—Failed Retaining Wall Surrounding the Village, 4—Localized Landslides along the gully wall, 5—Landslide on Northern Slope of the Village, 6—Haul Road and Exposed Pipeline, 7—Large Slope Failure North of the Village.
Geosciences 09 00332 g016
Figure 17. Location 1: Progression of structural damage and cleanup in the densely populated area of Pescara del Tronto; (a) September 2016, (b) December 2016, and (c) July 2018.
Figure 17. Location 1: Progression of structural damage and cleanup in the densely populated area of Pescara del Tronto; (a) September 2016, (b) December 2016, and (c) July 2018.
Geosciences 09 00332 g017
Figure 18. Location 2: Progression of landslides above Strada Statale; (a) September 2016, (b) December 2016, and (c) July 2018.
Figure 18. Location 2: Progression of landslides above Strada Statale; (a) September 2016, (b) December 2016, and (c) July 2018.
Geosciences 09 00332 g018
Figure 19. Location 3: Progression of failed retaining wall with location of cross section shown; (a) September 2016, (b) December 2016, and (c) July 2018.
Figure 19. Location 3: Progression of failed retaining wall with location of cross section shown; (a) September 2016, (b) December 2016, and (c) July 2018.
Geosciences 09 00332 g019
Figure 20. Location 3: Slope cross section comparison between the three models: September 2016, December 2016, and July 2018.
Figure 20. Location 3: Slope cross section comparison between the three models: September 2016, December 2016, and July 2018.
Geosciences 09 00332 g020
Figure 21. Location 4: Time progression of localized landslides along the gully wall in; (a) August 2016, (b) November 2016, and (c) July 2018.
Figure 21. Location 4: Time progression of localized landslides along the gully wall in; (a) August 2016, (b) November 2016, and (c) July 2018.
Geosciences 09 00332 g021
Figure 22. Location 5: Time progression of landslides on the northern slope of the village; (a) August 2016, (b) November 2016, and (c) July 2018.
Figure 22. Location 5: Time progression of landslides on the northern slope of the village; (a) August 2016, (b) November 2016, and (c) July 2018.
Geosciences 09 00332 g022
Figure 23. Location 6: haul road and exposed pipelines near the gravel pit (a) August 2016, (b) November 2016, and (c) July 2018.
Figure 23. Location 6: haul road and exposed pipelines near the gravel pit (a) August 2016, (b) November 2016, and (c) July 2018.
Geosciences 09 00332 g023
Figure 24. Location 7: large slope failure north of village (a) August 2016, (b) November 2016, and (c) July 2018.
Figure 24. Location 7: large slope failure north of village (a) August 2016, (b) November 2016, and (c) July 2018.
Geosciences 09 00332 g024
Table 1. Data from Nadir Grid-Based and Optimized Models.
Table 1. Data from Nadir Grid-Based and Optimized Models.
GCPTotal PhotosSurface Area ( km 2 )GSD (cm/px)Adjusted GSD (cm/px)
Nadir Grid229120.352.02.0
Optimized127160.412.41.5
Table 2. Quantitative Comparison of Photo Count per Point Between the Nadir Grid-Based and Optimized Models.
Table 2. Quantitative Comparison of Photo Count per Point Between the Nadir Grid-Based and Optimized Models.
Points TestedMeanStandard DeviationMinimum NumberMaximum Number
Nadir Grid3028.1317.47667
Optimized3023.206.781038

Share and Cite

MDPI and ACS Style

Freeman, M.; Vernon, C.; Berrett, B.; Hastings, N.; Derricott, J.; Pace, J.; Horne, B.; Hammond, J.; Janson, J.; Chiabrando, F.; et al. Sequential Earthquake Damage Assessment Incorporating Optimized sUAV Remote Sensing at Pescara del Tronto. Geosciences 2019, 9, 332. https://doi.org/10.3390/geosciences9080332

AMA Style

Freeman M, Vernon C, Berrett B, Hastings N, Derricott J, Pace J, Horne B, Hammond J, Janson J, Chiabrando F, et al. Sequential Earthquake Damage Assessment Incorporating Optimized sUAV Remote Sensing at Pescara del Tronto. Geosciences. 2019; 9(8):332. https://doi.org/10.3390/geosciences9080332

Chicago/Turabian Style

Freeman, Michael, Cory Vernon, Bryce Berrett, Nicole Hastings, Jeff Derricott, Jenessa Pace, Benjamin Horne, Joshua Hammond, Joseph Janson, Filiberto Chiabrando, and et al. 2019. "Sequential Earthquake Damage Assessment Incorporating Optimized sUAV Remote Sensing at Pescara del Tronto" Geosciences 9, no. 8: 332. https://doi.org/10.3390/geosciences9080332

APA Style

Freeman, M., Vernon, C., Berrett, B., Hastings, N., Derricott, J., Pace, J., Horne, B., Hammond, J., Janson, J., Chiabrando, F., Hedengren, J., & Franke, K. (2019). Sequential Earthquake Damage Assessment Incorporating Optimized sUAV Remote Sensing at Pescara del Tronto. Geosciences, 9(8), 332. https://doi.org/10.3390/geosciences9080332

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