Landslide inventory of Puerto Plata, Santiago de los Caballeros and Jarabacoa (Dominican Republic)

ABSTRACT Landslide inventory maps are the first and main information level in studies of reactivation of previous landslides and susceptibility, hazard and risk maps of future landslides. Due to its geographical location and tectonic setting, the Dominican Republic is exposed to several geological risks, among which landslides stand out. In this study, we present a landslide inventory elaborated in three selected zones with particular geographical, lithological, tectonic and climatic characteristics: Puerto Plata, Santiago de los Caballeros and Jarabacoa. A total of 6085 landslides have been recorded (56% in S, 29% in P and 15% in J), which have been classified according to their geometry, propagation mechanism, affected material, thickness and activity. The detected factors that control the distribution and type of gravitational deposits are lithology and geologic structure, as well as rainfall and active faulting along the Camú, Septentrional and Bonao-La Guácara fault zones.


Introduction
The Dominican Republic, due to its geographical location and tectonic setting on the northern boundary of the Caribbean Plate, is exposed to several geological risks, among which landslides are very frequent (Muñoz et al., 2009;2017;Roopnarine et al., 2020;Urien, 2011).Landslides are one of the most frequent natural hazards worldwide, causing numerous personal injuries and material damages (Haque et al., 2019).
Landslide inventory maps are fundamental in the management of these gravitational processes (van Westen et al., 2008).In this study, we present a landslide inventory of three selected zones of the Dominican Republic (Puerto Plata, Santiago de los Caballeros and Jarabacoa), made from previous geological and geomorphological maps (SYSMIN Project), photointerpretation of satellite images (Google Earth Pro; http://earth.google.com)and fieldwork.
The two main targets of this landslide inventory are: (i) making a map with a GIS-based database, describing their geometry, detection, propagation mechanism, affected material, thickness and activity and (ii) analysing the relationship between landslides and their destabilizing factors as lithology, structure, rainfall or active faulting.

Study zones
This study is focused in three zones: Puerto Plata, Santiago de los Caballeros and Jarabacoa.These zones have diverse lithological, tectonic and climatic attributes, combined with a high population density and high economic activity.These characteristics have special significance to compare landslide inventory maps and may be relevant for further analyses of landslide susceptibility.
The Hispaniola Island (second largest island of the Greater Antilles) is located in the northern boundary of the Caribbean Plate, which is moving towards the E and ENE, relative to the North American Plate (see The Geology of Hispaniola in the Main Plate).This displacement, predicted in the plate models, has been confirmed by GPS data (Symithe et al., 2015).The convergence between both plates produced an oblique collision between the Caribbean island-arc and the southern continental margin of North America in the Eocene (Mann et al., 1991;Pérez-Estaún et al., 2007), which evolved from the lower Miocene to a regime of transpressive tectonics (Dolan et al., 1998;Mann et al., 1995).From uppermost Pliocene to present day, these processes induced the uplift of the Cordillera Central, the sinistral transtension of the Cibao basin and the sinistral transpression of the Cordillera Septentrional (Escuder-Viruete & Pérez, 2020).
The study area is located in the northern sector of Hispaniola (Dominican Republic), which morphologically is constituted by two mountain ranges of WNW-ESE to NW-SE trend, separated by a depression.The northern and southern ranges are designated Cordillera Septentrional and Cordillera Central, respectively.Their main elevations are Pico Duarte (3098 m) in Cordillera Central (highest mountain in the Greater Antilles) and Loma Diego de Ocampo (1249 m) in Cordillera Septentrional.The intermediate depression is occupied by the Cibao sedimentary basin.
Due to the present-day transpressive and transtensive tectonic regime, the study area contains several large-scale active fault zones characterized by the frequent occurrence of earthquakes of variable magnitude.The Puerto Plata (P) zone is located on the north coast of the Cordillera Septentrional and is affected by the Camú Fault Zone (CFZ) (de Zoeten & Mann, 1991;Pindell & Draper, 1991).The Santiago (S) zone is located on the boundary between the Cordillera Septentrional and the Cibao basin, close to the trace of the Septentrional Fault Zone (SFZ) (Mann et al., 1998).The Jarabacoa (J) zone is located on the northern slope of the Cordillera Central, at the contact between the Bonao-La Guácara Fault Zone (BGFZ) and the Jarabacoa basin (Escuder-Viruete, 2011;Mann et al., 1991).Figure 1 shows panoramic views of the three zones, the Main Plate includes minor maps of these faults, and the lithological units affected.
The northern sector of Hispaniola, located at 19-20° north latitude, has a subtropical climate, controlled by marine influence, northeast trade winds and the mountain ranges.In this regard, the annual mean temperature and the annual accumulated precipitation are in Puerto Plata, at sea level on the windward side of the Cordillera Septentrional, 25,4° and 842 mm; in Santiago de los Caballeros, at 175 metres above sea level (masl) on the rain shadow of the Cordillera Septentrional, 24.2° and 760 mm; and in Jarabacoa, at 529 masl on the windward side of the Cordillera Central, 20.2° and 1723mm (Data collected between 1991 and 2021 by climate-data.org).These three zones are affected, from August to November, by tropical storms characterised by strong winds with heavy rainfall.

Puerto Plata zone
The Puerto Plata zone covers an area of 222.9 km 2 .The main urban areas are San Felipe de Puerto Plata, which is the third most populated city in the Dominican Republic and Imbert (Oficina Nacional de Estadística, 2012).Elevations range from 0 m on the shoreline to 766 m on Loma Isabel de Torres.The main fluvial courses are Camú, Bajabonico and Maimón rivers, and, shorter in length, San Marcos and Muñoz rivers.
The carbonated Villa Trina Formation occupies the main reliefs of the zone, as Loma Isabel de Torres (766 m) to the E (Figure 1A) and Loma El Brinco (575 m) to the SW, while the fine-grained siliciclastic units form the smooth reliefs and the valleys.

Santiago zone
The Santiago zone covers an area of 155.1 km 2 .This zone is located to the NNE of Santiago de los Caballeros, which is the second most populated city in the Dominican Republic.The main urban area is Tamboril (Oficina Nacional de Estadística, 2012).The main fluvial courses are Jacagua, Gurabo, Licey, Canca and Moca rivers.Elevations range from 219 m, just to the ESE of Canca river, to 1063 m, at the Licey river basin head.
The zone includes the following geological units: the sandstones and clays of La Toca Formation (Paleogene); the marls and nodular and calcareous siltites of Altamira Formation (Paleogene); the marls with mudstone and coralline limestone intercalations of Villa Trina Formation (Neogene); and the micritic and reef limestones of Los Haitises Formation (Neogene) (Ramírez Merino, 2010;Urien, 2011).
Tectonically, the SFZ is a WNW-striking oblique reverse to strike-slip fault (Figure 1C).The fault marks an S-facing scarp, which is occupied in its upper part by the carbonated rocks of Los Haitises and Villa Trina Formations and in the lower part by the fine-grained siliciclastic rocks of La Toca and Altamira Formations.

Jarabacoa zone
The Jarabacoa zone covers an area of 149.4 km 2 .In the zone, the main urban center is Jarabacoa, which is one of the 40 most populated cities in the Dominican Republic (Oficina Nacional de Estadística, 2012).The main fluvial courses are Yaque del Norte (longest river in the Greater Antilles), Baiguate and Jimenoa rivers as well as the Grande, Los Dajaos (I y II) and Cercado streams.Elevations range from 476 m, in the northern part of Jimenoa river, to 1526 m, between Baiguate river and the Cercado stream.
Two main geological units present in the zone are the massive basalts and locally microporphyritic and mafic tuffs of the Pelona-Pico Duarte Formation and the tonalites and leucotonalites of the Río and Junumucu batholits (Escuder-Viruete, 2011;Joubert, 2011).These units outcrop along a steep and freshlooking scarp conditioned by the BGFZ, which is a W to ENE-striking reverse oblique left-lateral fault (Figure 1D).

Methodology
The landslide inventory maps were developed according to the following stages: (1) Review of bibliographical documentation that includes: 6 geological maps at scale 1:50,000 in Puerto Plata zone (Hernáiz Huerta, 2010;Monthel, 2010) (2) Review of previous landslide reports (Muñoz et al., 2009;2017) and inventories.In 2013, the SGN initiated a landslide inventory (M.Calzadilla Ribeaux (deceased in 2019), personal communication, Departamento de Geología Ambiental y Aplicada, SGN).To date, the available database includes an inventory of 851 landslides, 7% of which are located in our study zone (Puerto Plata zone).
(3) Photointerpretation of satellite images Google Earth Pro (http://earth.google.com) at scale 1:5000.The comparison of different temporal images, from 2000 to 2021.(4) Fieldwork at scale 1:25,000 of Puerto Plata, Santiago and Jarabacoa zones (from 30th October to 19th November 2021).A drone DJI Mini 2 was used to obtain a photographic record of the main landslides.( 5) Integration of all the information, delineation of landslides and their storage in a Geographical Information System (ESRI ArcGIS 10.8).For the three zones, the topographic bases were an orthophoto from 2020 at 10 m resolution (https://docs.microsoft.com/en-us/bingmaps)and a Digital Terrain Model at 12 m resolution (LandInfo, Worldwide Mapping Limited Liability Company, Colorado, USA).From these DTMs, a contour model with 10 m of equidistance and a slope model have been derived.

Landslide detection
Depending on the detection method used, landslides have been differentiated into three types: (1) Photo-interpreted.Detected landslides using satellite images (Google Earth Pro; http://earth.google.com).
(2) Revised.Detected landslides, from previous landslides reports and inventories and satellite images, and confirmed through fieldwork.

Landslide classification
The next sections list and describe the different classifications used for the landslide inventory maps (abbreviations shown in brackets have been used in the database and in the Main Plate).Figure 2 shows some examples of the classified landslides.These classifications use both descriptive and interpretative criteria (Álvarez Muñoz, 2004).The descriptive or observable criteria include geometry, type, morphology or thickness of the affected materials, incision of drainage or density of vegetation.The interpretative criteria include propagation mechanism or activity of the detected landslides.These criteria are crucial for landslide recognition, interpreting their origin and prevent and mitigate the reactivation of existing landslides or the occurrence of new ones.

Geometry
The mapped landslides have two types of geometry: (1) Punctual.Landslides (less than 1000 m 2 ) mapped as a point in the Main Plate.
(2) Areal.Landslides (more than 1000 m 2 ) mapped as a polygon in the Main Plate.

Thickness
The landslides have been classified according to their thickness.It was no possible to estimate it in the fieldwork, so it was estimated attending to the contour model, enabling us to distinguish two types, which are briefly described below: (1) Shallow.Landslides not detectable in the contour model of 10 m equidistance.(2) Deep-seated.Landslides detectable in the contour model of 10 m equidistance and characterized by topographic criteria as scarps and concavity-convexity relationships.

Activity
Activity, or specifically state of activity, describes what is known about the timing of landslides (Cruden & Varnes, 1996).Internal morphology, as scarps or cracks, vegetation or dissection by drainage provides the first indications of landslide activity.Three distinct types were distinguished, which are briefly described below (in brackets are the classes corresponding to Keaton & DeGraff, 1996): (1) Active (active and dormant-young).Landslides with unvegetated or partly vegetated scarps and cracks.( 2) Inactive (dormant-mature).Landslides with different type or density of vegetation than adjacent terrain and slightly incised by drainage.
(3) Relict (relict and dormant-old).Landslides deeply incised by drainage network and same type and density of vegetation as adjacent terrain.

Landslide characteristics
The resulting inventory includes a total of 6085 landslides, 56% of which have been detected in Santiago zone, 29% in Puerto Plata zone and 15% in Jarabacoa zone (Table 1).
93% of shallow landslides do not intersect anthropogenic infrastructures so it is assumed most of landslides have little human disturbance.
The main results of the landslide analysis based on the classification criteria are discussed below.

Geometry
In the total of 6085 landslides, 96.3% are punctual and 3.7% areal.Additionally, the areal landslides (227 landslides) cover 85.1 km 2 , that is, the 16.1% of the area of the three zones (Table 2).

Propagation mechanism and affected material
From 308 classified landslides, the five main types are: rock slide flows (SFr), rock slides (Sr), earth flows (Fe), debris flows (Fd) and earth slide flows (SFe).
39.3% of the classified landslides are rock slide flows (SFr).They are associated with two types of scarps: lithological and related to faults.
Landslides associated to lithological scarps are conditioned by the contact between carbonated formations, more resistant to erosion, overlying siliciclastic formations, less resistant.In Puerto Plata zone, the contact between marls and limestones (Villa Trina Formation) overlying shales with blocks (San Marcos Mélange) or sandstones and clays (La Toca Formation) stands out as a lithological scarp controlled by gravitational processes.In Santiago zone, the contact between limestones (Los Haitises Formation) overlying marls and limestones (Villa Trina Formation), sandstones and clays (La Toca Formation) or marls and siltites (Altamira Formation) also gives place to a lithological scarp.In Jarabacoa zone, the basement is formed by igneous rocks, where there is no a clear lithological contrast that gives rise to greater resistance to erosion and forms a scarp.
Fault scarps with associated slide flows are those formed by the CFZ in Puerto Plata zone, the SFZ in Santiago zone and BGFZ in Jarabacoa zone.These scarps contain a highly fractured bedrock and therefore are likely to slide.
The other four types are described below.21.4% of the classified landslides are rock slides (Sr).These landslides are associated to high slopes and mainly carbonated formations.
18.5% of the classified landslides are earth flows (Fe).These landslides are associated to matrix-rich superficial formations.
11.4% of the classified landslides are debris flows (Fd).These landslides present source areas with high slopes, disponibility of poor-matrix superficial formations and a low relation source-area/deposit (Figure 3A).The deposits are clast-supported with angular and sub-angular clasts and centimetric (Figure 3B) to metric granulometries (Figure 3C).4.2% of the classified landslides are earth slide flows (SFe).These landslides are conditioned by the presence of shales with blocks (San Marcos Mélange) in Puerto Plata zone, sandstones and clays (La Toca Formation) in Santiago zone, alterites of igneous rocks and rich-matrix superficial formations in Jarabacoa zone.

Thickness
All the 5858 punctual landslides have been classified as shallow landslides and from a total of 227 areal landslides, 33.9% are shallow and 66.1% deep-seated.
Landslide activity can be correlated with relative age, reclassifying relict as very old and old landslides, inactive as mature landslides, and active as young and recent landslides (González-Díez et al., 1999;Keaton & DeGraff, 1996;Wieczorek, 1984).However, the high density of the tropical vegetation in the study zones can make this correlation inaccurate.Thus, the high percentage of inactive landslides can be correlated with old landslides, originated in a different climate (Bucci et al., 2021) or be recent landslides, covered by a dense tropical vegetation.

Analysis of destabilizing factors
Landslides can have several destabilizing factors, including precondition, preparatory, triggering and sustaining (Crozier, 1989).This understanding is beyond the scope of this paper.Hence, here the analysis will be confined to the main precondition and triggering factors.

Conditioning factors: lithology and structure
Lithology and structure are significant landslides precondition or conditioning factors.This section analyses the most significant lithologies and structures in the three zones.

Puerto Plata zone
In the total of 1751 landslides, 41.6% are located in marls and limestones (Villa Trina Formation) and 36.7% in shales with blocks (San Marcos Mélange) (Table 3).
Table 4 shows different parameters of the landslides affecting the marls and limestones (Villa Trina Formation) and the shales with blocks (San Marcos Mélange).
The inspection of Table 4 allows us to make the following considerations: (1) The marls and limestones (Villa Trina Formation) are associated to the steeper reliefs of the zone (Loma Isabel de Torres, 766 m and Loma El Brinco, 575 m), while the shales with blocks are associated to smooth reliefs (San Marcos Mélange).Thus, the carbonated Formation is interpreted as a lithology more resistant to erosion, conditioning deep-seated and more extensive slides (Mec, Am a ) and less gravitational activity (Act), while the remain siliciclastic Formation is interpreted as a lithology less resistant, conditioning shallow and less extensive flows (Mec, Am a ) as well as more gravitational activity (Act).(2) According to the stratigraphic sequence, the shales with blocks (San Marcos Mélange) underlie the marls and limestones (Villa Trina Formation).Thus, the less erosion resistance and stratigraphic position of the mainly siliciclastic lithology determine the gravitational instability of the mainly carbonated lithology.
(3) The number of punctual landslides in both lithologies is inversely proportional to the percentage of active landslides.In this sense, the marls and limestones (Villa Trina Formation), with more erosion resistance and less gravitational activity (Act), present a better preservation and a greater number of detected landslides (Nm P /km 2 ), while the shales with blocks (San Marcos Mélange), with less erosion resistance and more gravitational activity (Act), have a poorer preservation and a smaller number of detected landslides (Nm P /km 2 ).

Santiago zone
In the total of 3438 landslides, 41.6% are located in marls and siltites (Altamira Formation); 38.3% in sandstones and clays (La Toca Formation) and 8.4% in marls and limestones (Villa Trina Formation) (Table 5).Just 1% of the inventoried landslides are located in limestones (Los Haitises Formation).On the other hand, our fieldwork has revealed this lithology in many head-landslides affecting other formations (Villa Trina and La Toca Formations).
Table 6 shows different parameters of the landslides affecting the diverse lithologies.
From Table 6, we propose: (1) The sandstones and clays (La Toca Formation) are associated to the smooth reliefs.Thus, this unit is interpreted as a lithology less resistant to erosion and, consequently, a conditioning lithology of more gravitational activity (Act).(2) According to the stratigraphic sequence, the sandstones and clays (La Toca Formation) and the marls and limestones (Villa Trina Formation) underlie the limestones (Los Haitises Formation).Thus, the less erosion resistance and the  stratigraphic position of the carbonated-siliciclastic rocks determine the gravitational instability of the upper massive limestones.Consequently, in the contact between the carbonated and the siliciclastic lithologies occurs landslides, which are interpreted as deep-seated and extensive slides (Mec, Am a ), with frequent boulders of Los Haitises Formation in its heads.
Table 8 shows different parameters of the landslides controlled by the lithology.
Analyses of Table 8 allows us to make the following considerations: (1) The ultrabasic rocks and the basalts of the Pelona-Pico Duarte Formation are associated to the steeper reliefs of the zone, while the tonalite intrusions are associated to both steeper and smoother reliefs in the Jarabacoa basin.Thus, the ultrabasic rocks and basalts are interpreted as lithologies more resistant to erosion and, consequently, are lithologies of less gravitational activity (Act).However, the tonalites are interpreted as a lithology less resistant and, consequently, they determine shallow flows (Mec) and a more gravitational activity (Act).The lower resistance of tonalites may be related to their frequent supergene alteration (Figure 4).( 2) The basalts of the Pelona-Pico Duarte Formation present a highly fractured area characterized by shallow landslides, a steep relief without any associated concave uppermost part (Figure 1D) and outcrops of highly fractured basalts (Figure 5A), occasionally disaggregated (Figure 5B), adjacent to the Bonao-La Guacara Fault Zone.These data allow interpreting this area as affected by relatively shallow landslides but with a high susceptibility to develop deep-seated landslides.

Triggering factors: rainfall and active faulting
Intense rainfall and tropical storms are main triggering factors of landslides in the northern sector of Hispaniola.These factors were proposed by Muñoz et al. (2009Muñoz et al. ( , 2017) ) in the Carlos Díaz landslide (a deepseated and dormant-young landslide in the Santiago zone) or by Suárez Rodríguez et al. (2012) in La Vega zone (which includes the Jarabacoa study zone).Additionally, the Camú, Septentrional and Bonao-La Guácara fault zones are morphologically characterized by rectilinear scarps with associated both deepseated and extensive landslides.These three fault zones contain active faults and have shallow seismicity, therefore, the related scarps have been interpreted as active by Suárez Rodríguez et al. (2017b;2019) and Escuder-Viruete et al. (2021).
Taking into account both factors, intense rainfall and earthquake shaking, the landslide inventory here  elaborated can differentiate, in a first approach, which factor is more relevant.
In the total of 6085 landslides, 89% of the headlandslides are on the lowermost parts of the hillslopes and 10% on the uppermost parts (Table 9).
The observation of Table 9 allows us to establish the following considerations: (1) Following Densmore and Hovius (2000), landslide distribution can provide information on whether bedrock landslides in mountainous regions are triggered by storms or by earthquakes.
According to these authors, landslides closer to ridge crests are triggered by earthquakes and those concentrated on hillslope toes are triggered by storms.Thus, according to our map and data from Table 9, intense rainfall can be the triggering factor of the most of the landslides in the three zones (70% in P, 99% in S, 90% in J).

Discussion and conclusion
The compiled landslide inventory and the related database constitute an updated record of landslide activity in three selected zones of the Dominican Republic: Puerto Plata, Santiago and Jarabacoa.The combined photo-interpretation of satellite imagery and fieldwork made possible to detect 6085 landslides, 5858 punctual and 227 areal.Areal landslides affect the 16.1% of the area of the whole zones, corresponding to 20.0% in Santiago, 17.5% in Puerto Plata and 10.0% in Jarabacoa.The fieldwork campaign carried out in October-November of 2021 allowed to revise or identify 154 landslides, of which 81 are punctual and 73 areal.
The analysis of 308 classified landslides shows that shallow flows mainly affect to units less resistant to erosion.These units include the shales with blocks (San Marcos Mélange) in Puerto Plata zone; the tonalites intrusions (and related alterites) in Jarabacoa zone; and the sandstones and clays (La Toca Formation) and marls and siltites (Altamira Formation) in Santiago zone (see Table 6).
The deep-seated and extensive slides mainly affect to units that are more resistant to erosion and the lithological contact between more resistant rocks overlying less resistant rocks.Typically, these contacts are the marls and limestones (Villa Trina Formation) overlying the shales with blocks (San Marcos Mélange) in Puerto Plata zone; the limestones (Los Haitíses Formation) overlying the sandstones and clays (La Toca Formation) or marls and limestones (Villa Trina Formation) in Santiago zone; and the basalts (Pelona-Pico Duarte Formation) and the ultrabasic unit in Jarabacoa zone.
Rainfall is proposed as the first triggering factor over the active faulting, as indicate that 89% of the Table 9. Number and percentage of head-landslides on the lowermost parts of the hillslopes (Nm l ), on the uppermost parts (Nm u ) and with no data (null) in Puerto Plata (P), Santiago (S), Jarabacoa (J) zones and in all of them (T).head-landslides are concentrated on hillslope toes and the 10% closer to ridge crests.The higher percentage of landslide in Santiago zone (56.5% of the total landslides) stands out, compared to Puerto Plata zone (28.8%) and Jarabacoa zone (14.7%).The reasons for this difference may include the higher susceptibility to landslides of the marls and siltites of Altamira Formation (23.5% of the total landslides) and the sandstones and clays of La Toca Formation (21.6%).Future work should analyse the fault activity and the fluvio-gravitational evolution of the Septentrional Fault Zone in Santiago zone, comparing it with the Camú Fault Zone in Puerto Plata zone and the Bonao-La Guácara Fault Zone in Jarabacoa zone.

Software
We used Google Earth© (http://earth.google.com)for the visualisation of the landslides and ESRI ArcGIS 10.8 version to digitally organise and analyse the data, to use the Digital Terrain Model (LandInfo, USA), and to assemble the final map for publication and graphical editing.

Figure 1 .
Figure 1.Panoramic views: (A) rock slides, rock slides flows and earth flows at Loma Isabel Torres (766 m), main elevation of P zone; (B) rock slide flows at the CFZ scarp (P zone); (C) rock slide and rock slide flows at the SFZ scarp (S zone); (D) rock slide flows, dry debris flows, dissected debris flows and highly fractured area (hfa) at the BGFZ scarp (J zone).Letters refer to different propagation mechanisms and affected material (see Section 4.3).Dashed and continuous lines delineate, respectively, scarps and boundaries of landslides and continuous line with perpendicular strokes delineates the BGFZ scarp.Lower left key indicates the study zone and panoramic location quadrant (MecMat maps in the Main Plate) and lower right key the panoramic orientation.

Figure 2 .
Figure 2. Panoramic views: (A) Rock falls and rock slides (Loma Isabel de Torres, P zone); (B) Earth slide (P zone); (C) Earth slide flow (P zone); (D) Rock fall, rock slide and rock slide flow related to river erosion (S zone); (E) Rock slide flows and debris avalanche in the CFZ scarp (P zone); (F) Rock slide and rock slide flows perched above the river system (J zone); (G) Detail of the head of Carlos Díaz rock slide flow (S zone); (H) Outcrop of a rock slide flow (J zone).Letters refer to different propagation mechanisms and affected material (see Section 4.3).Dashed and continuous lines delineate, respectively, scarps and boundaries of landslides.Lower left key indicates the study zone and panoramic location quadrant (MecMat maps in the Main Plate) and lower right key the panoramic orientation.

Figure 3 .
Figure 3. (A) Panoramic views of some debris flows mapped in Jarabacoa zone.(B, C) Outcrops of some interpreted debris flows in Jarabacoa zone.Letters Fd refer to debris flow.Continuous line delineates the boundaries of the debris flow.Lower left key indicates the study zone and panoramic location quadrant (MecMat maps in the Main Plate) and lower right key the panoramic orientation.

Figure 5 .
Figure 5. Outcrops basalts in the BGFZ (Jarabacoa zone): (A) fractured basalts not disaggregated and (B) fractured basalts disaggregated.Lower left key indicates the study zone and panoramic location quadrant (MecMat map in the Main Plate).

Figure 4 .
Figure 4. Alterites and fresh rock cores in the tonalites (J zone).Lower left key indicates the study zone and panoramic location quadrant (MecMat map in the Main Plate).

Table 2 .
Area (A, in km 2 ) of Puerto Plata (P), Santiago (S), Jarabacoa (J) zones and in all of them (Total).The number of areal landslides (Nm a ), area (Am a , in km 2 ) and percentages relative to A (%) are also included.

Table 3 .
Lithology, number of landslides (Nm) and percentage (%) in Puerto Plata zone.Abbreviation C. refers to Complex.

Table 4 .
Parameters of the landslides affecting marls and limestones (Villa Trina Formation) and shales with blocks (San Marcos Mélange) in Puerto Plata zone.

Table 5 .
Lithology, number of landslides (Nm) and percentage (%) in Santiago zone.Abbreviations Fm. and Mb.refer to Formation and Member, respectively.

Table 6 .
Parameters of the landslides affecting limestones of Los Haitises Formation (Lime.),marlsandlimestones of Villa Trina Formation (Marl.),sandstonesandclays of La Toca Formation (Sand.) and marls and siltites of Altamira Formation (Silt.) in Santiago zone.Mec: most frequent propagation mechanism, landslide percentage and more frequent thickness (shallow or deep-seated).Nm p /km 2 : number of punctual landslides per km 2 .Am a : area affected by areal landslides.Act: percentage of active landslides.Mec and Act has been calculated according to classified landslides (119 landslides), Nm p /km 2 according to punctual landslides (3354 landslides) and Am a according to areal landslides (84 landslides) of Table1.

Table 8 .
Parameters of the landslides affecting the tonalites, the ultrabasic rocks and basalts of the Pelona-Pico Duarte Formation in Jarabacoa zone.Mec: most frequent propagation mechanism, landslide percentage and more frequent thickness (shallow or deep-seated).Nm p /km 2 : number of punctual landslides per km 2 .Am a : area affected by areal landslides.Act: percentage of active landslides.Mec and Act has been calculated according to classified landslides (45 landslides), Nm p /km 2 according to punctual landslides (863 landslides) and Am a according to areal landslides (33 landslides) of Table1.