Karst Geology of Aguijan and Tinian, CNMI Cave Inventory and Structural Analysis of Development
Kevin W. Stafford, John E. Mylroie & John W. Jenson
Department of Geosciences Mississippi State University Mississippi State, MS 39762
Water and Environmental Research Institute of the Western Pacific University of Guam Mangilao, Guam 96923
Technical Report No. 106
Water & Environmental Research Institute of the Western Pacific University of Guam
September, 2004
Karst Geology of Aguijan and Tinian, CNMI Cave Inventory and Structural Analysis of Development
Kevin W. Stafford, John E. Mylroie & John W. Jenson
Department of Geosciences Mississippi State University Mississippi State, MS 39762
Water and Environmental Research Institute of the Western Pacific University of Guam Mangilao, Guam 96923
Technical Report No. 106
Water & Environmental Research Institute of the Western Pacific University of Guam
September, 2004
THE WATER RESOURCES RESEARCH INSTITUTE PROGRAM OF THE US GEOLOGICAL SURVEY, AWARD NO. 01HQPA0010, SUPPORTED THE WORK REPORTED HERE. THE CONTENT OF THIS REPORT DOES NOT NECESSARILY REFLECT THE VIEWS AND POLICIES OF THE DEPPARTMENT OF THE INTERIOR, NOR DOES THE MENTION OF TRADE NAMES OR COMMERCIAL PRODUCTS CONSTITUTE THEIR ENDORSEMENT BY THE UNITED STATES GOVERNMENT.
2 ABSTRACT
Tinian and Aguijan, Commonwealth of the Northern Mariana Islands (CNMI), are volcanic, back arc islands in the western Pacific formed by Pacific Plate subduction under the Philippine Plate. The islands are composed of Eocene volcanic cores mantled by Plio-Pleistocene carbonate facies and raised Holocene beach and reef deposits. The entire sequence has been tectonically uplifted and contains high-angle normal faults, while isostatic subsidence and scarp failures overprint tectonic brittle failure features. A cave and karst inventory on Tinian and Aguijan surveyed 114 features and is believed to adequately represent the megaporosity (cave) development. Two distinct cave classes were identified: mixing zone caves (flank margin caves and banana holes) and fissure caves. Most mixing zone caves were located in or near scarps and coastlines, often at similar elevations to nearby caves. Fissure caves were located in regions of brittle failure, forming linear features with narrow widths. Three previous sea-level positions were identified based on horizons of mixing zone caves. Seventeen freshwater discharge sites and four allogenic recharge sites were identified on Tinian. Kolmogorov-Smirnov statistical analyses and rose diagram comparisons of orientation trends found significant similarities between megaporosity and geologic structure (brittle failure) on Tinian. Analyses of small regions showed distinct relations between brittle deformation and megaporosity, while at larger scales similarities became less obvious due to the complex geologic history and physiography of the island. Based on similarities in populations of orientation trends, fissure cave development is primarily controlled by brittle failure deformation with development along faults, fractures, and joints, while mixing zone cave development is primarily controlled by fresh-water lens position but significantly influenced by brittle failure deformation. Tinian and Aguijan do not fit neatly into one classification of the Carbonate Island Karst Model. Regions of Tinian best fit the Simple, Carbonate-Cover and Composite Island Karst Models, but none easily fit the entire island. Aguijan must be classified as a Simple Carbonate Island because no geologic data has proved the presence of non-carbonate rocks interfering with the fresh-water lens, however it is probable that Aguijan does contain basement rocks that extend above sea-level as on other carbonate islands in the Marianas.
i PREFACE
WERI (Water and Environmental Institute of the Western Pacific) Tech Report #96 was published as a preliminary assessment of the karst development and water resources of Tinian and Rota, CNMI (Commonwealth of the Northern Mariana Islands), based on a field reconnaissance conducted on the islands in June 2002. Since that time, work in the CNMI has focused on the Islands of Aguijan, Rota and Tinian with a primary emphasis on the cave and karst inventories of these islands, which can be used to evaluate the hydrolgeologic evolution of the islands. The work presented in this study provides a detailed cave and karst inventory of the Islands of Aguijan and Tinian, CNMI with an analysis of structural controls on cave and karst development on the island of Tinian. The authors conducted this work as thesis research at Mississippi State University in conjunction with the University of Guam, while the US Geological Survey, through the National Institutes for Water Resources Research program, award no. 01HQGR0134, funded fieldwork. The data and findings in this report were originally published as a master’s thesis by Kevin Stafford, titled: Structural Controls on Megaporosity in Eogenetic Carbonate Rocks: Tinian, CNMI. Ultimately, the work presented herein was made possible through the support of the Municipality of Tinian and Aguijan through: Tinian Mayors Office, CNMI Department of Land and Natural Resources, CNMI Department of Historical Preservation, Commonwealth Utilities Corporation, and Northern Marianas College. Due to the large volume of data in this report, the appendices (Appendix A: Color Figures; Appendix B: Cave and Karst Inventory; Appendix C: Orientation Data; and Appendix D: Statistical Comparisons) are published as .pdf files in the attached CD or can be downloaded from the WERI website (www.uog.edu/weri). These appendices include color images of significant features, detailed maps of the inventoried karst features on Aguijan and Tinian, diagrams representing the structural geology of Tinian, and data matrices of statistical analyses. Additional research continues in the Mariana Islands, in order to further answer questions about the karst geology and water resources of the islands. A similar study to the work reported herein is currently underway by the authors and master’s thesis research by T. Montgomery Keel at Mississippi State University. Ultimately this study and future studies will be used to evaluate the fresh-water lens development on carbonate islands, which can be used to better manage the water resources of these communities.
ii
CONTENTS
ABSTRACT ...... i
PREFACE ...... ii
CONTENTS...... iii
LIST OF TABLES ...... v
LIST OF FIGURES...... vi
INTRODUCTION...... 1
LITERATURE REVIEW...... 2
Geographical and Geological Setting ...... 2 Historical Setting ...... 5 Carbonate Island Karst...... 6 Karst Features ...... 10 Epikarst ...... 10 Closed Depression...... 11 Caves...... 12 Discharge Features...... 15 Water Resources ...... 15
STUDY METHODOLOGY...... 17
Initial Site Investigation...... 17 Data Collection ...... 18 Data Reduction ...... 19 Statistical Comparison of Data ...... 22 Small-Scale Test Site Evaluation...... 23
STUDY RESULTS ...... 23
Cave and Karst Inventory ...... 23 Cave Orientations...... 24 Brittle Deformation...... 26 Scarps and Coastlines ...... 26 Rose Diagrams...... 27 Statistical Comparison ...... 28 Tinian Composite...... 28 Central Plateau ...... 29 Northern Lowland...... 29 North-Central Highland ...... 30 Median Valley...... 30 Southeastern Ridge...... 30 Small-Scale Test Sites...... 31 Carolinas Limestone Forest...... 31 Puntan Diapblo...... 31 Unai Dangkolo ...... 31
iii
CONTENTS (cont.)
DISCUSSION AND CONCLUSIONS...... 32
Tinian Cave and Karst Inventory...... 32 Aguijan Cave and Karst Inventory ...... 35 Controls on Cave and Karst Development...... 35 Island Scale Comparisons...... 38 Province Scale Comparisons...... 38 Small-Scale Test Site Comparisons...... 39 Structural Control of Caves ...... 40 Karst Development on Aguijan and Tinian ...... 40
SUMMARY ...... 41
REFERENCES CITED ...... 44
APPENDIX (attached CD)
A. Color Figures...... 48 B. Cave and Karst Inventory: Maps and Descriptions...... 63 C. Orientation Data...... 162 D. Statistical Comparison ...... 248
iv
LIST OF TABLES
Table Page
1 Cave and karst features surveyed on Tinian, CNMI: UTM location, physiography, cave type and geology ...... 65 2 Cave and karst features surveyed on Aguijan, CNMI: UTM location, physiography, and cave type...... 67 3 Fissure cave primary orientations ...... 162 4 Fissure cave segment orientations with segment length ...... 164 5 Mixing zone cave primary orientations...... 168 6 Mixing zone cave segment orientations with segment length...... 170 7 Mixing zone cave, entrance width segment orientations with segment length..... 178 8 Mixing zone cave, penetration segment orientations with segment length...... 182 9 Mixing zone cave, maximum width segment orientations with segment length...... 186 10 Fault orientations...... 190 11 Joint orientations...... 195 12 Fracture orientations measured during fieldwork ...... 197 13 Inland scarp segment orientations with segment length...... 200 14 Coastal scarp segment orientations with segment length...... 205 15 Coastline segment orientations with segment length ...... 213 16 Legend for column and row headings used in statistical comparison data matrices ...... 324 17 Tinian Composite statistical comparison data matrix ...... 325 18 Central Plateau statistical comparison data matrix ...... 326 19 Median Valley statistical comparison data matrix ...... 327 20 Northern Lowland statistical comparison data matrix ...... 328 21 North-Central Highland statistical comparison data matrix...... 329 22 Southeastern Ridge statistical comparison data matrix...... 330 23 Carolina's Limestone Forest statistical comparison data matrix ...... 331 24 Puntan Diapblo statistical comparison data matrix...... 332 25 Unai Dangkolo statistical comparison data matrix ...... 333
v
LIST OF FIGURES
Figure Page
1 Location map of the Mariana Islands...... 2 2 Physiographic map of Tinian, CNMI...... 3 3 Physiographic map of Aguijan...... 4 4 Geology section...... 4 5 Carbonate Island Karst Model ...... 7 6 Schematic diagram illustrating the relationship between island size, basement relationships to sea level and surface, and relative sea level within the Carbonate Island Model ...... 9 7 Conceptual model of cave types that form in eogenetic rocks on carbonate islands...... 13 8 Model showing the fresh-water lens morphology showing the location of basal and parabasal waters and the Ghyben-Herzberg principle...... 16 9 Example of primary and segment orientation trends measured using the apparent trend method for a typical fissure cave...... 20 10 Example of primary and segment orientation trends measured using the apparent trend method for a typical mixing zone cave ...... 20 11 Example of entrance, maximum width, and penetration trends measured using the entrance width trend method for a typical mixing zone cave ...... 20 12 Flank margin caves develop more complicated morphologies as they grow in size ...... 32 13 Conceptual model for the growth of flank margin caves ...... 33 14 Horizons of flank margin cave development show previous fresh-water lens positions ...... 34 15 Diagram showing the relationship between cave widths and lengths, which represent two distinct populations for fissure caves and mixing zone caves ...... 36 16 Diagram showing the relationship cave entrance width and cave maximum width...... 36 17 Coastline orientations for Tinian show a wide range of trends because of the elliptical shape of the shape of the island...... 37 18 Map of Unai Dangkolo region showing the close proximity of flank margin cave and fissure cave development ...... 43 19 Inland and coastal scarps on Tinian ...... 48 20 Areas of potential allogenic recharge...... 49 21 AMCS standard cave symbology...... 50 22 Location of test site areas...... 51 23 Location of cave and karst features surveyed on Aguijan and Tinian ...... 52 24 Fresh-water discharge sites located on Tinian ...... 53 25 Bamboo growing in the North Lemmai Recharge Feature with vines coating the scarp that forms the non-carbonate / carbonate contact ...... 54 26 Ponded water in the South Lemmai Recharge ...... 54 27 Location of closed depression on Tinian...... 55 28 Active quarry on Tinian ...... 56 29 Fresh-water at the land surface at Hagoi in the Northern Lowland ...... 56 30 Faults and non-carbonate rock outcrops...... 57 31 South Unai Dangkolo represents a typical cove...... 58 32 Hidden Beach Cave demonstrates well the transition from flank margin cave to cove resulting from coastal erosion...... 58
vi
LIST OF FIGURES (cont.)
Figure Page
33 Typical flank margin cave morphology ...... 59 34 Multiple levels of mixing zone dissolution exist on Tinian with at least three identified in the North-Central Highland near Mount Lasu...... 59 35 Fissure caves form narrow, linear passages that appear to be developed along zones of brittle failure ...... 60 36 Flowstone deposits on ceilings and walls ...... 60 37 Insect Bat Cave on Aguijan represents a paleo-discharge feature ...... 61 38 Scallops on the ceiling and walls of Liyang Atkiya, Aguijan...... 61 39 Diagram showing the location of modern carbonate beach deposits and primary brittle failure types...... 62 40 Map of "600 Meter" Fracture System ...... 68 41 Map of Almost Cave ...... 69 42 Map of Andyland Cave ...... 69 43 Map of Anvil Cave...... 70 44 Map of Barcinas East Cave...... 71 45 Map of Barcinas West Cave...... 72 46 Map of Barely Cave ...... 72 47 Map of Bee Hooch Cave...... 73 48 Map of Biting Mosquitoes Cave ...... 74 49 Map of Body Repel Cave...... 74 50 Map of Boonie Bee Sink...... 75 51 Map of Broken Stal Cave...... 76 52 Map of Cabrito Cave...... 76 53 Map of Cannon Cave ...... 77 54 Map of Carolinas Fracture Cave ...... 78 55 Map of Cave Without a Cave...... 79 56 Map of Cave Without a Roof...... 79 57 Map of Cavelet Cave...... 80 58 Map of Central Mendiola Cave Complex ...... 81 59 Map of Cetacean Cave ...... 82 60 Map of Chiget Fracture ...... 83 61 Map of Cobble Cave ...... 84 62 Map of Coconut Trap Cave...... 84 63 Map of Command Post Cave Complex...... 85 64 Map of Cowrie Cave ...... 86 65 Map of CUC Cave...... 86 66 Map of Danko's Misery...... 87 67 Map of Death Fracture Complex...... 88 68 Map of Diamond Cave...... 89 69 Map of Dos Cenotes Cave...... 90 70 Map of Dos Sakis Cave Complex ...... 91 71 Map of Dove Cave ...... 91 72 Map of Dripping Tree Fracture Cave...... 93 73 Map of Dump Coke Cave ...... 94 74 Map of Dynasty Cave...... 95 75 Map of East Suicide Cliff Cave ...... 96 76 Map of Edwin's Ranch Cave...... 96
vii
LIST OF FIGURES (cont.) Figure Page
77 Map of Elevator Cave ...... 96 78 Map of False Floor Cave...... 97 79 Map of Five Bee Cave Complex...... 98 80 Map of Flamingo Tail Caves...... 99 81 Map of Fleming Point Cave ...... 100 82 Map of Full Bottle Cave...... 101 83 Map of Gecko Cave ...... 101 84 Map of Goat Cave ...... 102 85 Map of Goat Fracture Cave...... 103 86 Map of Half-Dozen Cave...... 103 87 Map of Headless Tourist Pit...... 104 88 Map of Hermit Crab Cave...... 105 89 Map of Hidden Beach Cave ...... 106 90 Map of Hollow Column Cave ...... 106 91 Map of Insect Bat Cave...... 107 92 Map of Isotope Cave ...... 108 93 Map of John's Small Cave...... 109 94 Map of Lasu Recharge Cave...... 110 95 Map of Leprosy Caves ...... 111 96 Map of Leprosy Discharge Feature...... 111 97 Map of "600 Meter" Fracture System ...... 113 98 Map of Liyang Barangka ...... 114 99 Map of Liyang Dangkolo...... 115 100 Map of Liyang Diapblo...... 116 101 Map of Liyang Gntot...... 117 102 Map of Liyang Lomuk ...... 118 103 Map of Liyang Mohlang ...... 119 104 Map of Liyang Popporput ...... 120 105 Map of Liyang Sampapa ...... 120 106 Map of Liyang Umumu...... 121 107 Map of Lizard Cave ...... 121 108 Map of Lower Suicide Cliff Cave Complex ...... 123 109 Map of Masalok Fracture Cave Complex ...... 124 110 Map of Mendiola Arch Cave Complex...... 125 111 Map of Metal Door Cave ...... 126 112 Map of Metal Spike Cave Complex...... 127 113 Map of Metal Stretcher Cave ...... 128 114 Map of Modified Cave...... 129 115 Map of Monica Wants to be Like Kevin Cave...... 129 116 Map of Natural Arch Cave...... 130 117 Map of North Unai Dangkolo ...... 132 118 Map of Northern Playground Cave ...... 133 119 Map of Nuestra Señora de Santo Lourdes Cave Complex...... 134 120 Map of Orange Cave ...... 135 121 Map of Orphan Kids Cave Complex...... 136 122 Map of Pebble Cave...... 137 123 Map of Pepper Cave...... 137 124 Map of Piña Cave Complex ...... 138
viii
LIST OF FIGURES (cont.)
Figure Page
125 Map of Playground Cave ...... 138 126 Map of Plunder Cave ...... 139 127 Map of Radio Inactive Cave ...... 140 128 Map of Red Snapper Cave ...... 141 129 Map of Rock Hammer Cave ...... 141 130 Map of Rogue Cave ...... 143 131 Map of Rootcicle Cave...... 143 132 Map of Scorpion Cave ...... 144 133 Map of Screaming Bat Cave ...... 144 134 Map of Skip Jack Cave ...... 145 135 Map of Skull Cave Complex...... 146 136 Map of Skylight Cave ...... 147 137 Map of Solitary Cave ...... 147 138 Map of South Mendiola Cave ...... 148 139 Map of South Unai Dangkolo ...... 149 140 Map of Spider Cave ...... 149 141 Map of Swarming Termites Cave ...... 150 142 Map of Swiftlet Cave ...... 151 143 Map of Swimming Hole Cave Complex...... 152 144 Map of Toppled Column Cave...... 153 145 Map of Tridactid Cave Complex...... 152 146 Map of Twin Ascent Caves...... 155 147 Map of Unai Chiget...... 156 148 Map of Unai Lamlam...... 156 149 Map of Unai Masalok...... 157 150 Map of Water Cave ...... 158 151 Map of Waypoint Cave ...... 159 152 Map of West Lasu Depression ...... 160 153 Map of West Suicide Cliff Cave Complex...... 161 154 Rose diagrams of fissure cave primary orientations ...... 163 155 Rose diagrams of fissure cave, five-meter segment orientations ...... 166 156 Rose diagrams of fissure cave, ten-meter segment orientations ...... 167 157 Rose diagrams of mixing zone cave, primary orientations ...... 169 158 Rose diagrams of mixing zone cave, five-meter segment orientations...... 173 159 Rose diagrams of mixing zone cave, ten-meter segment orientations ...... 174 160 Rose diagrams of composite cave, primary orientations...... 175 161 Rose diagrams of composite cave, five-meter segment orientations ...... 176 162 Rose diagrams of composite cave, ten-meter segment orientations...... 177 163 Rose diagrams of mixing zone cave entrance width orientations ...... 179 164 Rose diagrams of mixing zone cave entrance width, five-meter segment orientations...... 180 165 Rose diagrams of mixing zone cave entrance width, ten-meter segment orientations...... 181 166 Rose diagrams of mixing zone cave penetration orientations...... 183 167 Rose diagrams of mixing zone cave penetration, five-meter segment orientations...... 184
ix
LIST OF FIGURES (cont.)
Figure Page
168 Rose diagrams of mixing zone cave penetration, ten-meter segment orientations...... 185 169 Rose diagrams of mixing zone cave maximum width orientations...... 187 170 Rose diagrams of mixing zone cave maximum width, five-meter segment orientations...... 188 171 Rose diagrams of mixing zone cave maximum width, ten-meter segment orientations...... 189 172 Rose diagrams of fault orientations...... 192 173 Rose diagrams of fifty-meter, fault segment orientations ...... 193 174 Rose diagrams of one hundred-meter, fault segment orientations ...... 194 175 Rose diagrams of joint orientations...... 196 176 Rose diagrams of orientations of fractures measured during fieldwork ...... 199 177 Rose diagrams of inland scarp orientations...... 202 178 Rose diagrams of inland scarp, fifty-meter segment orientations...... 203 179 Rose diagrams of inland scarp, one hundred-meter segment orientations...... 204 180 Rose diagrams of coastal scarp orientations ...... 207 181 Rose diagrams of coastal scarp, fifty-meter segment orientations...... 208 182 Rose diagrams of coastal scarp, one hundred-meter segment orientations...... 209 183 Rose diagrams of composite scarp orientations ...... 210 184 Rose diagrams of composite scarp, fifty-meter segment orientations...... 211 185 Rose diagrams of composite scarp, one hundred-meter segment orientations...... 212 186 Rose diagrams of coastline orientations...... 215 187 Rose diagrams of coastline, fifty-meter segment orientations ...... 216 188 Rose diagrams of coastline, one hundred-meter segment orientations...... 217
x
KARST GEOLOGY OF AGUIJAN AND TINIAN, CNMI Cave Inventory and Structural Analysis of Development
INTRODUCTION
This study has three primary objectives: eogenetic limestone (Jenson et al., 2002; 1) conduct a cave and karst inventory on Mylroie and Jenson, 2002; Mylroie and Tinian and Aguijan, because no inventory Jenson, 2001). It provides a fundamental existed, 2) investigate the cave and karst and systematic framework for describing development on Tinian to determine if hydrogeologic karst evolution in young influences on mega-porosity development carbonate rocks, by incorporating the effects are associated with brittle deformation in of mixing zone dissolution, glacio-eustacy, eogenetic rocks, and 3) evaluate of the tectonism, and lithologic variations. In this islands of Tinian and Aguijan in relation to model distinctive cave morphologies the Carbonate Island Karst Model (Jenson et develop as eogenetic karst. These al., 2002; Mylroie and Jenson, 2002), in distinctive morphologies include flank order to further advance the understanding margin caves, banana holes, pit caves, and of eogenetic karst development on carbonate stream caves. Traditionally, the affects of islands. structural and lithologic controls on The Tinian municipality of the eogenetic karst development have been Commonwealth of the Northern Mariana greatly overlooked because the original Islands (CNMI) governs Tinian and model was developed to describe features Aguijan. The Mariana Islands are a observed in the structurally and volcanic, back arc island chain in the lithologically simple islands of the Bahamas western Pacific formed by Pacific Plate (Mylroie and Carew, 1995a). subduction under the Philippine Plate along Structural control of karst development the Mariana Trench. Aguijan covers 7.2 is common in continental settings, where the square kilometers and Tinian covers 102 rock exhibits low porosity and is highly square kilometers. A geologic study was fractured (Klimchouk and Ford, 2000). conducted on Tinian in 1960 (Doan et al., Eogenetic karst is generally associated with 1960), which described the island as an caves formed by mixing zone dissolution, Eocene volcanic core mantled by Plio- primarily independent of structural controls, Pleistocene coralliferous and algal carbonate in rocks that have never been buried beyond facies and raised Holocene beach and reef the range of meteoric diagenesis and retain deposits. The entire sequence has been high primary and syngenetic porosities. uplifted and contains high-angle normal Recently, Jenson and coworkers (2002) faults produced from arc tectonism. Tinian recognized the importance of structural and has experienced greater than 100 meters of lithologic controls on carbonate island karst uplift since the Pleistocene, with an and modified the Carbonate Island Karst estimated 1.8 meters of uplift in the Model to account for the effects of Holocene (Dickenson, 1999). Aguijan has deformation and lithology on dissolution only been studied briefly (Tayama, 1936) that occurs on carbonate islands that are or and no geologic map has been produced for were tectonically active and/or have intricate the island; however, it is presumed to have a interfingering of carbonate and non- depositional and tectonic history similar to carbonate rocks. Tinian based on its proximity (~9 kilometers The cave and karst inventory of this southwest of Tinian). study was conducted in two intensive field The Carbonate Island Karst Model seasons based on information gathered predicts cave and karst development in during a reconnaissance survey in June 2002
1
(Stafford et al., 2002). During fieldwork, test). Using similar techniques, Nelson caves and karst features were documented, (1988) and Barlow and Ogden (1992) have surveyed and classified by cave type. Sites shown that continental, telogenetic karst of fresh-water discharge and allogenic development is related to the orientations of recharge were located. Orientations of regional faults, fractures, and lineaments, brittle failure structures were measured. implying that cave development is primarily After fieldwork was complete, maps were controlled by brittle deformation. produced of the cave and karst features that Evaluation of the cave and karst features were surveyed. Maps of features on Tinian on Tinian and Aguijan allows them to be were analyzed to determine cave orientation. evaluated with the Carbonate Island Karst The cave orientation data was compared to Model. However, the lack of previous fault orientations, scarp orientations, and geologic studies on Aguijan and the coastline orientations for similarity. Data complex tectonic history of Tinian do not comparisons were evaluated as independent allow for simple classification of these data pairs both visually (rose diagrams) and islands. statistically (Kolmogorov-Smirnov 2-sample
LITERATURE REVIEW
GEOGRAPHICAL AND maximum elevation of 157 meters (Figure GEOLOGICAL SETTING 3). Both islands have wet-dry tropical The Mariana Islands are located in the climates with a distinct rainy season (July- western Pacific Ocean and comprise a total September) and dry season (February- of 17 islands (Figure 1). Guam is the largest March). Annual rainfall averages 200 and southern-most island in the Marianas centimeters and temperature ranges from 20 and the only one that is not politically o to 32o Celsius (Gingerich and Yeatts, 2000; affiliated with the Commonwealth of the Butler, 1992; Tracey et al., 1964; Doan et Northern Mariana Islands (CNMI). The al., 1960; Cloud et al., 1956). relative position of the islands north of Guam in order of increasing distance is: Rota, Aguijan, Tinian, Saipan and Medinilla, which are carbonate islands. The remainder of the island chain is volcanic (Cloud et al., 1956). Tinian is located approximately 3000 kilometers east of Asia, approximately 180 kilometers north, northwest of Guam, and approximately 10 kilometers south of Saipan, while Aguijan is located approximately 9 kilometers southwest of Tinian. Tinian (Latitude: 15.01oN, Longitude: 145.62oE) has a surface area of 102 square kilometers with 51.2 kilometers of coastline and a maximum elevation of 187 meters (Figure 2). Aguijan (Latitude: 14.85oN, Longitude: 145.57oE) has a surface area of 7.2 square kilometers Figure 1: Location map of carbonate with 12.4 kilometers of coastline and a islands in the Marianas.
2
The Mariana Ridge, on which the convergence of the Pacific Plate to the east islands of Tinian and Aguijan are located, is and the Philippine Plate to the west. The formed along a volcanic arc located Mariana Islands are situated on an older approximately 160 kilometers west of the island arc that is separated by the Mariana Mariana Trench, which is the world's Trough from the younger Mariana West deepest trench with a maximum depth of Ridge, approximately 300 kilometers to the 11,035 meters (Gross, 1982). This west, which is developing in the Mariana subduction zone is created at the back-arc basin. This shift in ridge development was initiated by a change in Pacific Plate subduction geometry approximately 43 million years ago (mya) (Dickinson, 1999; Mink and Vacher, 1997; Reagan and Miejer, 1984). Doan and coworkers (1960) used topography and spatial relationships to divide Tinian into five physiographic regions: Northern Lowland, North- Central Highland, Central Plateau, Median Valley and Southeastern Ridge (Figure 2). The Northern Lowland comprises the broad, flat, nearly horizontal surface that slopes gently upward from the west coast to Sabanettan Chiget. Located above the Central plateau and midway between the east and west coasts is the North-Central Highland, which contains the highest point (162 meters) in northern Tinian at Mount Lasu. The Central Plateau includes the central portion of the island and is isolated by steep slopes and bounding scarps associated with north-south faults. In the south and east-central regions, the Median Valley expresses little relief, but forms a broad depression Figure 2: Physiographic map of Tinian, CNMI, with important bounded by faults. The features and locations identified (adapted from Doan et al., 1960). Southeastern Ridge, which
3
includes Kastiyu, the highest point on Tinian (187 meters), is developed on two principal fault blocks (Doan et al., 1960). Because little geologic work has been performed on Aguijan (Butler, 1992; Tayama, 1936), the following physiographic provinces are proposed based on the classification system used on Rota, which exhibits a similar terraced topography (Sugiwara, 1934): Upper Terrace, defined by elevations greater than 100 meters, which form a broad, relatively flat plateau and reaches a maximum elevation of 157 meters; Middle Terrace, defined by elevations between 50 and 100 meters, which is nearly absent on the northern side of the island and best developed on the southeastern side; and Lower Terrace, defined by elevations less Figure 4: Geology section (adapted from Doan et al., 1960). than 50 meters, which includes the steep cliffs that form the coastline and the lowest pyroclastics were probably ejected from a bench. This classification is proposed here, submerged vent. The limestone units are in order to establish a distinction between subdivided into two formations: Tagpochau karst development located in different Limestone and Mariana Limestone. The regions on Aguijan (Figure 3). Tagpochau Limestone covers 16% of surface area of Tinian and is early Miocene in age. The formation is composed of three contemporaneously deposited facies: detrital (Tt), argillaceous (Tta), and sandy (Tts). The detrital facies comprises the majority of the formation and is composed primarily of fragments of biogenic calcium carbonate with calcite cement. The other two facies comprise only a small portion of the Tagpochau Limestone. The Mariana Limestone covers 83% of the surface area of Tinian and was deposited in Pliocene and Figure 3: Physiographic map of Aguijan. Pleistocene time. The formation was subdivided based upon the presence of The geology of Tinian was described by constructional or detrital compositions Doan and coworkers (1960) and remains the into seven facies: constructional most detailed geologic study for the island. coralliferous facies (Qtmcc), constructional Tinian is composed of volcanic tuffs and algal facies (QTmca), detrital coralliferous breccias covered with coralline and algal facies (QTmc), detrital shelly facies limestone (Figure 4). The igneous rocks (QTms), detrital Halimeda facies (QTmh), (Tinian pyroclastic) comprise less than three detrital argillaceous facies (QTmu), and square kilometers and retain only relict detrital undifferentiated facies (QTmu). structures and textures because of extensive Overlying these deposits in coastal regions weathering. Based on the presence of are Holocene limestones, developing sands foraminifera within the sediments, the and gravels, and reefs (Siegrest, 1988; Doan
4
et al., 1960; Burke, 1953). Recent fieldwork control lasted until World War One (Hunt indicates that late Pleistocene limestones and Wheeler, 2000). may also overlie the Mariana Limestone, In 1914 Japan took control of the because some limestone present in coastal Mariana Islands at the start of World War areas appears to correspond with oxygen One. After the war, Japan's control as isotope stage 5e, suggesting that Mariana administrator was recognized by the League Limestone deposition does not extend to the of Nations and remained so until World War Pleistocene/Holocene boundary (Stafford et Two. During the Japanese occupation, vast al., 2002). tracts of land were cleared of coconut trees While no detailed study of the geology and tropical forest, so that the extensive of Aguijan has been conducted, the same sugar cane plantations could be developed classification for rock units used on Tinian covering 58% (Bormann, 1992) to 90% will be applied, based on its close proximity, (McClure, 1977) of the land area. In 1944, until future work produces a more detailed United States military forces secured the geologic map of Aguijan. In addition, it is islands in some of the bloodiest battles of presumed that Aguijan is similar to other World War Two. Tinian became islands in the southern Marianas, in that a instrumental in America’s Pacific war non-carbonate core of pyroclastic rocks campaign. The northern third of the island exists beneath the exposed carbonate rocks was developed into a large airbase named that crop out, while the development of three North Field, consisting of four major distinct terrace levels on Aguijan suggests a runways. During full operation, this was the similar tectonic history (Stafford et al., in largest and busiest airfield in the world with press; Doan et al., 1960; Sugiwara, 1934). two B-29’s taking off simultaneously every 45 seconds on mission days. From North HISTORICAL SETTING Field, the two atomic bombs, Fat Man and The Mariana Islands were first settled Little Boy, were assembled and loaded into by the Chamorro people around 1500 B.C.E. the B-29 bombers Enola Gay and Bock’s (before common error). In 1521, Ferdinand Car for drops on Hiroshima and Nagasaki, Magellan, discovered the Marianas and respectively (McClure, 1977). named them Islas de los Ladrones (Islands In 1947, the United Nations gave the of the Thieves), which were then renamed United States trusteeship over the Mariana Islas de las Marianas (Islands of the Islands, which continued until the 1970’s. Marianas) in the early 1600’s after Maria In 1970 discussion began between the Ana of Austria, the widow of the Spain’s Marianas and the United States over the king, Phillip the IV. The islands remained termination of the trustee agreement, which in Spanish control for almost four centuries, led to a plebiscite in 1975, which negotiated during which time much migration by a covenant with the Unites States. The indigenous people from the Caroline Islands Covenant to Establish a Commonwealth of occurred. In 1899, after Spain’s loss of the the Northern Marianas Islands in Political Spanish-American war, the Spanish Union with the United States of America controlled islands of Micronesia (Mariana negotiated an agreement where the citizens Islands, Caroline Islands and Palau) where of the Northern Mariana Islands became sold to Germany. During German control, self-governing as a Commonwealth, while the northern Marianas greatly developed Guam remained a United States territory. agricultural and fishing economies by The covenant enabled the Commonwealth of immigrating additional Carolinian people the Northern Mariana Islands (CNMI) to from Chuuk and Yap, as well as Japanese retain the benefits of U.S. citizenship, people from nearby islands like Okinawa, excluding the right to vote in U.S. who developed a lucrative agricultural trade presidential elections, which ensured market of copra and sugar cane. German continued economic support. In exchange,
5
the U.S. was given lease of 75 square 1. The fresh water/salt water boundary kilometers of land in the Northern Mariana creates mixing dissolution, and Islands, including over two thirds of Tinian produces organic-trapping horizons at (Hunt and Wheeler, 2000). both the upper and lower boundaries Today, Tinian and Aguijan are governed of the fresh-water lens. by a single municipality as part of the 2. Glacio-eustacy has moved the fresh- Commonwealth of the Northern Mariana water lens up and down through a Islands (CNMI). Aguijan is uninhabited and vertical range of over 100 m in the Tinian hosts a population of approximately Quaternary. 2000 individuals. While the majority of the 3. Local tectonics can overprint the Tinian remains under U.S. military control, glacio-eustatic sea level events, small farms and tourism, including scuba adding complexity to the record. diving and casinos, provide commerce 4. Carbonate islands can be divided into (Bormann, 1992). four categories based on basement/sea level relationships: CARBONATE ISLAND KARST i. Simple carbonate islands (no non- This study enables the studies in karst carbonate rocks). geology and hydrology performed on Guam ii. Carbonate cover islands (non- (Mylroie et al., 2001) and in the Caribbean carbonate rocks beneath a carbonate and Atlantic (Frank et al., 1998; Mylroie and veneer). iii. Composite islands (carbonate and Carew, 1995a; Mylroie et al., 1995a) to be non-carbonate rocks exposed on the applied and evaluated to different settings. surface). A primary objective of investigation on iv. Complex islands (faulting and facies Tinian is to advance the understanding of interfingering create complex the karst hydrology of carbonate islands, carbonate/non-carbonate while refining the general Carbonate Island relationships). Karst Model (CIKM). The Carbonate Island 5. The karst is eogenetic, i.e., it has Karst Model has been designed as the developed in carbonate rocks that are definitive model for the hydrologic young and have never been buried development of island karst. Modern below the range of meteoric carbonate islands are unique due to diagenesis. extensive interaction between fresh and saline groundwater within young, porous Carbonate islands described by the rock, which produces a unique geologic and Carbonate Island Karst Model are composed hydrologic history that is different than that of young limestones and are heavily in continental settings (Vacher and Mylroie, influenced by meteoric waters and mixing- 2002; Mylroie and Jenson, 2002; Mylroie et zone dissolution, creating an environment al., 2001; Mylroie and Vacher, 1999, and for eogenetic karst to develop because of references therein). their close proximity to the site of Karst forming in marine conditions on deposition. Vacher and Mylroie (2002, p. carbonate coasts and islands can be 183) define eogenetic karst as “the land explained by the Carbonate Island Karst surface evolving on, and the pore system Model (Figure 5) (Stafford et al., 2003; developing in, rocks undergoing eogenetic, Jenson et al., 2002; Mylroie and Jenson, meteoric diagenesis.” Eogenetic carbonate 2002; Mylroie and Jenson, 2001). Its main rocks are diagenetically young and have not aspects are: undergone extensive cementation and compaction, as opposed to telogenetic carbonate rocks that are diagenetically mature, have been buried below the range of meteoric diagenesis, and are currently
6
Figure 5: Carbonate Island Karst Model (CIKM) (adapted from Mylroie et al., 2001).
7
exposed at the earth surface by tectonic on islands like the Bahamas, where the rock uplift and erosion of overlying strata is diagenetically young, fresh-water interacts (Klimchouk and Ford, 2000). Typically, extensively with the salt-water, and the small carbonate islands are the sites of fresh-water lens has been significantly eogenetic karst undergoing meteoric affected by glacio-eustacy. Therefore, diagenesis, however, it is not possible to islands that have telogenetic karst, are state that carbonate islands develop solely removed from the affects of glacio-eustacy, eogenetic karst and continental settings and/or do not exhibit mixing zone develop solely telogenetic karst, because dissolution should have karst landforms and fresh-water / sea-water interaction and a features reported as karst on islands, while fresh-water lens affected by glacio-eustacy islands that exhibit eogenetic karst, are are also required for the formation of affected by glacio-eustacy, and exhibit carbonate island karst, as defined by the mixing zone dissolution should have karst Carbonate Island Karst Model (Jenson et al., landforms and features reported as island 2002; Mylroie and Jenson, 2002). karst (Vacher and Mylroie, 2002). Eogenetic karst development occurs on In eogenetic rocks that have not islands such as Bermuda, the Bahamas, and undergone compaction and cementation, the the Marianas, as well as in continental rocks tend to initially exhibit high matrix settings such as the Biscayne aquifer of porosity and moderate permeability, but can Florida, where the carbonate rocks are develop secondary vuggy porosity as a result young. Telogenetic karst development of meteoric and fresh-water diagenesis. The occurs on islands such as Gotland (Sweden) matrix porosity will generally decrease with and Kephallenia (Greece), as well as in age as secondary cementation infills the pore continental settings such as Kentucky space, while permeability increases as (United States) and England (Great Britain), preferential flow develops extensive where the rocks are of appreciable age and horizontal flow routes. Over time, the bulk have been buried beyond the range of porosity remains the same or decreases, but meteoric diagenesis (Vacher and Mylroie, the vertical hydraulic conductivity of the 2002). Because eogenetic karst and island rock decreases while the horizontal karst are not necessarily synonymous, a hydraulic conductivity increases (Vacher distinction should be made between island and Mylroie, 2002). Due to the proximity of karst and karst on islands. Island karst eogenetic karst to marine waters in island occurs in specific environments that must settings, it is highly susceptible to changes include three basic parameters: 1) in sea level, which results in the migration interaction of fresh-water and sea-water, of the fresh-water lens and flow routes in which produces mixing zone dissolution, 2) response to glacioeustatic and tectonic a fresh-water lens that was affected by changes. glacioeustatic changes in the Quaternary, The four conceptual model and 3) karst that is eogenetic. Karst on classifications of carbonate islands are based islands may exhibit some of the on island composition (Figure 5). Simple characteristics of island karst, but these carbonate islands are composed of only islands are not true island karst unless all carbonate rocks at the surface and to a depth three basic parameters are present. Islands below the base of the fresh-water lens like Jamaica and Puerto Rico have karst (Figure 5A). The Bahama Islands are a landforms (e.g. cockpits and mogotes) in good example. Carbonate cover islands are their interior that are developed in rocks of composed of only carbonate rocks at the appreciable age and that are not influenced surface, but have non-carbonate rocks that by glacio-eustacy or fresh-water / salt-water interact with the fresh-water lens without interaction. This makes the karst landforms being exposed at the surface (Figure 5B). karst on islands. True island karst occurs Bermuda is a good example. Carbonate
8
cover and simple carbonate islands can In traditional continental settings, easily shift from one form to the other as a structural and lithologic controls have been result of changes in relative sea level (Figure recognized as having a significant influence 6). Composite islands are composed of both on karst development. Klimchouk and Ford non-carbonate and carbonate rocks on the (2000, p. 57) state: “Bedding planes, joints, surface (Figure 5C). Barbados is a good and faults are planar breaks that serve as the example. These three island types represent principal structural guides for groundwater a classification scheme where the two end flow in almost all karstified rocks.” In member environments are islands island karst settings the significance of completely composed of carbonate rocks at structural and lithologic controls have often one end and islands completely devoid of been overlooked because of the dominant carbonate rocks at the other end, with role of mixing-zone dissolution in the intermediate compositions described by the hydrogeologic system and because the model (Figure 6). The fourth type are original models for island karst were complex islands, which are characterized as developed on simple carbonate islands in the having complex geologies as a result of Bahamas. Jenson and coworkers (2002) and faulting, the interfingering of different Mylroie and coworkers (2001) have recently lithologies, and the presence of both recognized the importance of structural and carbonate and non-carbonate rocks (Figure lithologic controls on carbonate island karst 5D). Guam and Saipan are excellent by modifying the CIKM with the addition of examples (Jenson et al., 2002; Mylroie et al., the fourth type, the complex island. This 2001). addition enables the effects of deformation
Figure 6: Schematic diagram illustrating the relationship between island size, basement relationships to sea level and surface, and relative sea level for simple, carbonate-cover, and composite islands within the Carbonate Island Karst Model (CIKM) (Mylroie and Jenson, 2001, Fig.2, p. 54).
9
and lithology to be accounted for on 2002). Faults are surfaces or narrow zones carbonate islands that are or were that show relative displacement similar to tectonically active and/or have intricate fractures, but over larger regions and interfingering of carbonate and non- generally associated with regional tectonics carbonate rocks. (Twiss and Moores, 1992). Faulting is Lithologic controls include variations in reported throughout Tinian, including the rock composition between rock formations boundaries between the Southeastern Ridge and within the same rock formation, where and the Median Valley and between the beds subdivide the unit. Variations may Central Plateau and the Northern Lowland, include grain size, grain origin, grain where fault blocks moved independently sorting, and chemical composition. These during island emergence (Doan et al., 1960). variations provide routes for preferential Joints, fractures, and faults generally dissolution within beds and formations of enhance fluid movement, but may also favorable composition, while less favorable restrict it. Enhanced flow can occur along beds and formations may restrict dissolution planar features that provide a surface for and fluid movement (Klimchouk and Ford, fluid movement both laterally and vertically 2000; Sasowsky and White, 1994; Palmer, across lithologic boundaries, creating greater 1991; White, 1988). connectivity within the subsurface. Structural controls include both brittle Restricted flow can result if recrystallization and ductile deformation. Ductile (slickensides) or secondary infilling (caliche deformation results in the folding of rock dikes) develop along the planar surface, bodies, which adds to the complexity of creating a barrier for fluid movement. Flow lithologic controls. Brittle deformation restriction may also occur if faulting results includes joints, fractures, and faults, which in the offsetting of different rock units, are often found in tectonically active, which produces contacts between carbonate carbonate islands where uplift and and non-carbonate rocks that did not subsidence create complex geology. previously exist (Klimchouk and Ford, Fractures are surfaces where the rock has 2000; Sasowsky and White, 1994; Palmer, been broken, occur over a range of 1991; White, 1988). centimeters to meters and fall into several categories, including: extension – relative KARST FEATURES rock motion is perpendicular to the surface; Karst features on Tinian can be shear – relative motion is parallel to the classified into four broad categories; surface; and oblique – relative motion is a epikarst, closed depressions, caves, and both parallel and perpendicular to the discharge features. Identification, surface. Joints are specific types of classification, and spatial distribution of fractures, which show only minor extension different karst morphologies provide a basis (Twiss and Moores, 1992). Fractures for understanding the hydrology of the reported on Tinian and Aguijan are region. primarily associated with gravity slides resulting from bank margin or cliff margin Epikarst failure along steep scarps. Joints reported Epikarst is the zone of dissolutional occur inland of the bank margin/cliff margin sculpturing (karren) that is present on the fractures, possibly as a result of rock surface and upper few meters of bedrock in expansion from a decrease in lateral pressure carbonate regions. Karren has been from the scarp failure or as unloading described as minor solutional forms, which structures associated with isostatic range from millimeters to meters in scale subsidence, and perpendicular to coastlines, (White, 1988). In general, epikarst is which may be associated with coastal independent of environmental setting; erosion or regional faulting (Stafford et al., however, in coastal regions, where salt spray
10
is in active contact with carbonate rocks, an the deep vadose zone via fractures and pits environment is created for the production of and supply water directly to the phreatic biokarst (Viles, 1988) or phytokarst (Folk et zone (Mylroie et al., 1999). Epikarst al., 1973); however, the biological affects development and extent is a primary are often overstated (Mylroie and Carew, controlling factor on the quantity of water 1995b). The surface of the karren in that enters phreatic storage via vadose paths eogenetic karst is generally extremely rough and it is possible that it may serve as a and may not support extensive soil profiles location of significant water storage (Jocson because soil-forming material is in short et al., 2002; Jocson et al., 1999; and Jenson supply. Taborosi and coworkers (in press) et al., 1997). Epikarst on Tinian appears have recently reported additional karren identical to that seen on Guam and Saipan morphologies on eogenetic karst on Guam. (Taborosi, 2000). In the Mariana Limestone on Guam most Individuals involved in the development weathered surfaces are extremely jagged and of land in karst regions dominated by reminiscent of biokarst or phytokarst and is epikarst should be aware that any attributed to a polygenetic origin that modifications to the land surface and includes mixing dissolution from rainwater epikarst would alter the drainage dynamics mixed with salt spray, salt weathering, of the area. Ponding basins may actually dissolution by meteoric waters, and exhibit lower infiltration rates because of biological weathering (Taborosi et al., in high sediment loads. Runoff events and soil press). However, other karren morphologies erosion may also increase as a result of the are also present, which include more modification of the natural landscape rounded karren in inland regions and (Mylroie and Carew, 1997). dissolution along brittle failure planes producing enlarged vertical and horizontal Closed Depressions joints as wells as solution pans (kaminetzas) In carbonate island environments, closed (Taborosi et al., in press). Carbonate rocks depressions can be classified into three on most islands generally have little general categories: dissolution, natural insoluble material to produce soils, construction, and human modification. suggesting that the insoluble material is These different types can be extremely hard exclusively of eolian origin. The soils that to differentiate based on appearance and it is are present are generally piped downward possible that the features may have been any into voids and dissolutional cavities within or all of the three classification types at the karst system. Eogenetic karst regions some point in their development. In have specific vegetation types because of the carbonate islands, dissolutional depressions limited development of soils, which affects are generally small to moderate in size as a infiltration rates of water entering the karst result of their young age and the nature of system based on soil thickness, composition autogenic recharge. However, in areas and its presence, absence or modification where non-carbonate rocks are exposed at (human development). Beneath karren, the surface, streams can develop that dissolutional bedrock debris, and soil, which provide allogenic recharge to the karst comprise the surficial epikarst, the system, possibly forming large depressions remainder of the epikarst zone is composed as are seen on Guam (Mylroie et al., 2001; of solutional fissures, holes, and shallow Taborosi, 2000). Natural construction small cavities in the bedrock (Mylroie et al., depressions are those that formed at the time 2001). Epikarst contributes water to the the rocks were deposited or are the result of deep vadose zone as diffuse flow and subsequent deformational processes. through the integration of flow paths created Mylroie and coworkers (2001) report that by cavities and fissures in the lower epikarst. natural construction depressions are the This integrated flow can effectively bypass most common form on simple carbonate
11
islands, and anthropogenic depressions (e.g. borrow pits, refuse disposal sites, artificial quarries, landfills, artificial drainage ponds, drainage ponds in residential areas, and and storage ponds) are often constructed in features associated with World War Two these pre-existing features (Mylroie et al., (bomb pits, defensive positions, etc.). Initial 1999). field investigations have revealed that Tinian exhibits areas where dissolution- significant natural depressions are found in type closed depressions are formed. association with volcanic outcrop exposures. Dissolution depressions have been seen at The large closed depressions seen elsewhere the contacts between exposed volcanic are probably natural construction features or outcrops and carbonate outcrops (Stafford et features produced by human modification, al., 2002). On Tinian, four exposures of possibly to pre-existing natural depressions, volcanic rocks are recorded near Sabanettan of the land surface (Stafford et al., 2002). Mangpang, Bañaderon Lemmai, and Laderan Apaka. At the three northern Caves locations (two at Sabanettan Mangpang and Caves are natural openings in the earth one at Bañaderon Lemmai), there are closed that can be characterized based on their size, depressions in the limestone outcrops, which shape, length and overall geometry (White, show allogenic streams descending into 1988). Solution caves, which have been them from the volcanic outcrops. An initial formed by the dissolution of bedrock by field investigation of Bañaderon Lemmai circulating groundwater, are present has shown that these features are similar to throughout the islands of Tinian and those seen on Guam and are providing point Aguijan. Solution caves on carbonate source recharge into the karst system islands can be grouped into various (Stafford et al., 2002). In the south-central categories. In the Mariana Islands, five region near Laderan Apaka, volcanic distinct cave types have been documented: outcrops are exposed on a north-facing cliff. banana hole, flank margin cave, fissure To the north of the outcrops, a large closed cave, pit cave, and stream cave (Figure 7). depression (Sisonyan Makpo) exists in In some areas, it is difficult to discern the which the Municipal Wells are located exact origin and original extent of some (United States Department of the Interior caves, because of the extensive modification Geological Survey, 1983), which may be a of some features for Japanese military natural construction feature formed by purposes during World War Two (Taborosi complex faulting on the island, based on its and Jenson, 2002). In addition to human location between the two prominent ridges modification, horizontal notches cut into (Kastiyu/Carolinas and Piña) of the cliff faces create an additional classification Southeastern Ridge and the low-lying problem. Traditionally these have been Median Valley. Weathering of volcanic identified as bioerosion notches, but may rocks and talus accumulation may be also form by lateral corrosion and cliff partially armoring the slopes of the cliff and retreat, or they may be the remnants of flank the closed depression, increasing recharge to margin caves (Mylroie et al., 1999). the closed depression and possibly leading Accurate classification of karst features is to lateral corrosion of the closed depression. integral to interpreting the hydrogeology of Throughout Tinian, human-modified island karst. depressions are also commonly reported Banana holes (Figure 7) are shallow, (United States Department of the Interior small chambers that formed at the top of the Geological Survey, 1983) and may represent fresh-water lens. They are isolated features modified constructional depressions similar and exhibit a morphology with a width to to those reported on other carbonate islands depth ratio greater than one (Harris et al., (Mylroie et al., 1999). These include 1995). On Tinian no feature has been modern and ancient quarries (latte stone), currently identified as a definite banana
12
Figure 7: Conceptual model of cave types that form in eogenetic rocks on carbonate islands. hole. A feature located several hundred complicating the morphology and meters to the northwest of the Lasso Shrine identification of flank margin caves are appears to have a complex history that could bioerosion notches formed by wave erosion be explained as the stacking of two or more and invertebrate borings at sea level during banana hole features vertically on one stillstands (Mylroie and Carew, 1991). another as sea level fluctuated. Later Flank margin caves have been reported collapse by upward stoping could then form at various locations. The most significant the collapsed areas seen in the lower levels flank margin cave development is in the of this feature (Stafford et al., 2002). Suicide Cliffs on the southern end of Tinian, Flank margin caves (Figure 7) are were numerous entrances can be seen inland formed in the distal margin of the freshwater from the road on the southern coastal lens where thinning of the lens and mixing terrace. These caves are located of fresh and saline waters creates approximately 50 meters upslope from the dissolutionally aggressive phreatic waters lower terrace and are developed along a (Mylroie et al., 1995a) and are the most consistent horizon that is believed to reflect common morphology type that has been a past sea-level stillstand. Along the reported on Tinian (Stafford et al., 2003, northern edge of the Southeastern Ridge, 2002). Flank margin caves exhibit globular near Laderan Apaka, a series of flank morphologies with a wide range of sizes that margin caves are located, which are less may be connected, but often remain as extensive due to more complete erosion and isolated chambers. Flank margin caves can slope retreat. Near Lasso Shrine on northern be used as indicators of previous sea-level Tinian, a series of modified flank margin stillstands and have potential as tools for caves are located in the eastern cliffs. Along evaluating differential rates of uplift in the east coast, flank margin caves that are at tectonically active carbonate islands (Carew various stages of erosion due to cliff retreat and Mylroie, 1995; Mylroie et al., 1995a,b). and coastal erosion are present near Unai Flank margin caves may collapse in coastal Masalok and Unai Dangkolo. At Unai regions to form coves and caletas, which Dangkolo, the largest flank margin cave retain only remnants of the original cave documented on Tinian is located morphology (Back et al., 1984). Further approximately 200 meters inland. This cave
13
is breached on the surface by ceiling descending meteoric waters and have a collapse and can only be entered from above depth to width ratio greater than one. Pit via a 10-meter descent. The cave consists of caves often develop as a series of shafts with several large chambers that intersect, connecting lateral sections, act as vadose creating a complex system that is laterally fast-flow routes for water entering the extensive covering and area greater than subsurface, and effectively drain the epikarst 1300 square meters with cave passages (Mylroie and Carew, 1995a). Although developed around a central chamber common in certain settings (Harris et al., approximately 15 meters tall and 35 meters 1995; Mylroie and Carew, 1995a), only one in diameter (Stafford et al., 2002). pit cave has been reported on Tinian. It is Fissure caves (Figure 7) are developed located 10 meters from the cliff edge on the along joints, fractures, or faults, in which south coast, where it connects to a low-level preferential flow along the planar surface bench just above sea level. It is has resulted in enhanced dissolution rates approximately 25 meters deep with an (White, 1988). Three basic types of fissure entrance width of 3 meters, but due to its caves have been identified, based on their location and near-direct connection to the morphology and spatial relationship to the ocean, it has little effect on aquifer recharge. island. At the coastal edge, development However, its presence indicates the perpendicular to the coastline has been possibility of other similar features in the reported, which in several cases has resulted area, which may be covered at the surface in caves that often discharge freshwater and by collapse and at shallow depths could be penetrate inland up to 30 meters as tubular open, acting as vadose fast flow routes passages with distinct joints located in their (Stafford et al., 2002). ceilings. The second type has formed as Stream caves (recharge caves, Figure 7), fractures parallel to the coastline or scarps as hydrologically active caves that are fed by a result of cliff-margin or bank-margin allogenic water, have been documented in failure. These fractures may reach the water regions were contacts between carbonate table where they can affect the flow and non-carbonate rocks occur. Meteoric dynamics of the lens by locally distorting water is unable to infiltrate efficiently into the lens or by providing vadose routes that the volcanic rocks, and forms surface intercept the normal diffuse flow in the streams, which channel water to the outcrop phreatic zone (Mylroie et al., 1995c; Aby, periphery where it descends into the 1994). Dissolutionally enhanced fractures carbonate rocks as point source allogenic have produced caves over 40 meters deep on recharge (Mylroie et al., 2001). Tinian. The third type occur where Tinian volcanic rocks outcrop in four dissolution along a fault plane, in regions and closed depressions are conjunction with collapse, has formed caves, associated with the periphery of each of which tend to develop at moderate to steep these outcrops. Due to time and logistics, angles along the dip of the fault, and the only one of these features, located to the east caves extend laterally along the strike of the of the largest outcrop at Bañaderon Lemmai, fault. All three types of fissure caves are was investigated in June 2002 (Stafford et hydrologically important because they al., 2002). This site contains a large provide vadose fast flow routes for water sinkhole that appears to be armored with through the subsurface, either as recharge or volcanic sediments that restrict infiltration, discharge features and may distort the local but on the eastern edge, a small cave is lens morphology and flow dynamics present that shows evidence of allogenic (Stafford et al., in press; Stafford et al., recharge during rain events (Stafford et al., 2002). 2002). It is expected that other closed Pit caves (Figure 7) are vertical shafts depressions associated with these volcanic that have developed by the dissolution of outcrops will exhibit similar features formed
14
by dissolution at the non- WATER RESOURCES carbonate/carbonate contact, with variations Fresh-water is partitioned in a lens dependant on sediment infilling and the above salt-water and is affected by the volume of allogenic recharge that is presence of non-carbonate rocks interacting received. with the fresh-water lens. Basal water occurs where the fresh-water lens is directly Discharge Features underlain by sea-water, while parabasal Discharge volume and types vary water occurs where the fresh-water lens is greatly on carbonate islands and spatially directly underlain by non-carbonate, within the same island. Three general basement rocks (Figure 8; Mink and Vacher, categories of discharge features can be 1997). Basal waters develop a lens defined: seeps, springs, and submarine thickness defined by the Ghyben-Herzberg freshwater vents. Of these types, seeps Principle: represent diffuse discharge, while the other types represent focused discharge points. Z = [(ρf)/(ρs - ρf)]h, Seeps occur extensively in coastal areas where calcareous sand covers the bedrock where Z is the depth of the fresh-water lens and disperses emerging water over a large below sea level, h is the height of the lens area. Springs are locations in which water above sea level, and ρf and ρs are the emerges from the bedrock along preferential densities of fresh-water and salt-water flow paths generally defined by bedding respectively (Figure 8; White, 1988; Raeisi planes or fractures along coastlines. Springs and Mylroie, 1995). The fresh-water/salt- may be laterally extensive if controlled by water lens thickens as the fresh-water head bedding planes or they may be restricted to a depresses the interface below sea level, single point discharge if controlled by relative to the density difference in the two fractures. Submarine freshwater vents are waters: areas where freshwater discharge occurs below tidal level along the island periphery Z = αh, (Jocson et al., 2002; Mylroie et al., 1999). Seeps and springs have been where α = (ρf)/(ρs - ρf). With ρf = 1:00 3 3 documented in coastal regions on Tinian. g/cm and ρs = 1.025 g/cm , α = 40; Minor seeps are present at Unai Dangkolo therefore, for every one meter or foot of and Unai Masalok on Tinian, while springs fresh-water head, the fresh-water/salt- have been documented in various coastal water interface is depressed 40 meters areas. On Tinian, fracture caves were or feet (1:40 ratio) (White, 1988; Raeisi reported along the west coast near Unai and Mylroie, 1995). Masalok where freshwater could be Island environments have limited water observed mixing with salt water and resources and specific problems because of discharging from solutionally widened the morphology of the fresh-water lens bedding planes that extend approximately 30 (water with a chloride concentration <250 meters inland (Stafford et al, 2002). mg/L (Gingerich and Yeatts, 2000)). Submarine freshwater vents are not reported Because of the limited extent of island for Tinian and Aguijan, but features are aquifers and the characteristics of eogenetic expected to exist that are similar to those karst, island aquifers are extremely that have been reported from Guam susceptible to contamination. (Mylroie et al., 2001; Jenson et al., 1997), Contamination from the surface can include: where freshwater discharges below sea human and animal wastes, fertilizers, level. detergents, pesticides, herbicides, petroleum spills, and solvent spills. Subsurface contamination can occur from salt-water
15
Figure 8: Model showing the fresh-water lens morphology showing the location of basal and parabasal waters and the Ghyben-Herzberg principle (adapted from Mink and Vacher, 1997 and Raeisi and Mylroie, 1995). intrusion, where over pumping produces a for military purposes, while other regions cone-of-depression in the lens, which then are less accessible due to terrain. Of the up-cones salt-water into the fresh-water lens inhabited region, 60% is an undeveloped, at a one to forty ratio. Because of the ratio public, rural area. The remaining 40% of of up-coned water to lens thickness above this region is composed of residential and sea level, a small decrease in the lens commercial lots that include a casino resort, thickness above sea level will create a large small businesses, farming, grazing and cone of depression at depth in basal waters housing. The residential and commercial (Gingerich and Yeatts, 2000; Mylroie and lots provide the largest concerns for Carew, 1997). When non-carbonate rocks contamination from surface spills and extend into or through the freshwater lens, biological waste disposal. Tinian currently basement rocks restrict upconing of has no sewer facility, instead human waste is parabasal water because salt-water does not disposed of through septic and seepage directly underlie fresh-water. Therefore, it tanks, leaching fields or holding tanks, is important for island communities to which may not be adequate for preventing effectively manage their water resources in groundwater contamination, because of the order to continue to produce potable water, thin soil profile and rapid infiltration rates by limiting water extraction at specific sites associated with eogenetic, carbonate rocks. and by utilizing zones of parabasal Other sources of contamination include the groundwater (Mink and Vacher, 1997). airport, several quarries and a solid waste In the 1995 census, 2,631 people were disposal dump, which provide paths for listed as living on Tinian primarily in the direct recharge into the aquifer by bypassing Median Valley and parts of the adjacent the soil surface and epikarst, which have Central Plateau, which comprises 25% of been removed (Gingerich and Yeatts, 2000; the total island surface area. This limited Bormann, 1992). region of population is partially the result of On Tinian, USGS investigations show the United States control of the northern that the maximum lens thickness is third of the island, the Northern Lowland, approximately 12 meters (40 feet) in the
16
center of the Median Valley, with slight With the stable fresh-water lens thinning near the Municipal Well and Marpi morphology on Tinian, the location and Marsh (Gingerich and Yeatts, 2000). The identification of karst features becomes Municipal Well (a Maui-type, infiltration important for preventing groundwater well) supplies the majority of the islands contamination. Understanding the spatial water needs at a rate of approximately 4.5 x distribution and extent of fissure caves, pit 106 liters/day (~1.2 Mgal/day). Today water caves, and recharge features, which can produced from the Municipal well has a transport contaminants rapidly to the lens, chloride concentration of 180 mg/L, which may enable government planners to regulate is 100 mg/L greater than when the well was activities near sensitive areas. Similarly, constructed in 1945. The USGS monitored identified discharge features can be used as the aquifer thickness in relation to rain sampling points for monitoring possible events in the wet season of 1993 and groundwater contamination. Identification reported an aquifer thickening of 90-150 cm of allogenic recharge may be used to predict (3-5 ft), while reports in the dry season of regions where parabasal waters may be 1994 showed an aquifer thinning of 30-60 thicker due to increased recharge, such that cm (1-2 ft). These results revealed that salt-water intrusion risks from water annual fluctuations are minor and lens extraction may be reduced. Spatial thickness is dependent more upon long-term distributions of flank margin and banana rainfall patterns instead of annual hole caves may provide insight into the fluctuations, for maintaining overall diffuse flow characteristics of previous morphology (Gingerich and Yeatts, 2000). fresh-water lenses, making it possible to This stable lens morphology indicates that better evaluate current lens morphology. If contamination problems have the potential positive correlations do exist between brittle to produce long term effects due to long deformation and megaporosity, then residence time for water within the aquifer. regional structure may provide insights into groundwater behavior.
STUDY METHODOLOGY
This study had three objectives. It was the study were used to evaluate the islands initially developed to inventory, survey, and of Tinian and Aguijan in relation to the classify the cave and karst features on Carbonate Island Karst Model. Tinian and Aguijan, because no such database existed for the islands. In addition INITIAL SITE INVESTIGATION to the inventory, this study was developed to The initial site investigation included a evaluate whether or not a statistical reconnaissance of Tinian and analysis of the comparison could be made between physiography and geology of the island. megaporosity and zones of brittle failure in The reconnaissance was conducted in June eogenetic rocks on carbonate islands, using 2002 (Stafford et al., 2002) and provided data collected from the island of Tinian basic information about the cultural, where basic geologic studies have been physical and logistical aspects of Tinian, conducted (Doan et al., 1960). The study while establishing relationships with local was conducted in five major phases, in order government bodies that would be crucial for to reach the two objectives: 1) initial site in-depth studies on the island. During this investigation, 2) data collection, 3) data reconnaissance areas that were reported by reduction, 4) statistical comparison of data, island residents as having significant cave and 5) small-scale test site evaluation. development were visited, including, but not Results from the two primary objectives of limited to: Suicide Cliffs, Unai Dangkolo,
17
and Mount Lasu. This initial site relation to outcrops of igneous rock. The investigation demonstrated that most known closed depressions that were found to be caves, reported by local residents and proximal to the igneous outcrops were hunters, occur predominantly along scarps, identified as locations of possible allogenic coastlines, and closed depressions, which recharge developed by dissolution, while the agreed with previous investigations on closed depressions that were distal to the carbonate islands (Stafford et al., 2002; igneous outcrops probably had non- Mylroie et al., 1999). dissolutional (i.e. constructional or human After the initial reconnaissance, the modified) origin (Figure 20, Appendix A). islands geography and geology were The reconnaissance and GIS analyzed with a geographical information investigation of Tinian provided a system (GIS) that was produced for the framework for fieldwork and data collection island using ArcView 3.2 (ESRI, 2000). A on Tinian. Because there is little published digital elevation model (DEM) was created geology on Aguijan and a reconnaissance for the island using spatial data transfer visit was not possible, no detailed fieldwork standard (SDTS) compliant raster data with plan could be developed. However, this 10-meter postings produced by the National project extended the cave and karst Mapping Program of the United States inventory to Aguijan when the opportunity Geological Survey (USGS, 2001a,b,c,d,e), for access to the island was available. which provides greater resolution than the 1:25,000 topographic map produced by the DATA COLLECTION USGS (USGS, 1983). The DEM was then Two types of data were collected during overlain with a scanned geology map fieldwork: 1) cave and karst surveys, and 2) produced by Doan and coworkers (1960). It structural orientations of zones of brittle was geo-referenced using Image Analyst failure. These two datasets were collected (ESRI, 2000) in order to scale and align the over the course of two intensive field geology map with DEM, then all igneous seasons (December 10, 2002 to January 07, outcrops and faults were manually digitized 2003 and May 04, 2003 to June 1, 2003) and in order to create shape files of these provided the database for analysis and features. The DEM and digitized geologic correlation in this study. During the course features were then used as the basic dataset of fieldwork, coastal and scarp needed to delineate prominent scarps and investigations were limited to regions that closed depressions on the island. could be accessed within an acceptable risk Scarps were defined as any change in level that would not greatly endanger safety. slope greater than twenty degrees, which Cave and karst surveys were conducted enabled the identification of all major scarps on Aguijan and Tinian based on the initial in the island interior and coastline (Figure site investigation and additional reports 19, Appendix A). Because of the DEM cell provided by local residents during the course size, smaller scarps that might contain cave of fieldwork. This fieldwork focused on entrances were excluded; however, it did not known caves, coastlines, scarps, and closed eliminate investigations of smaller coastal depressions. Caves were surveyed in scarps because coastlines were identified as accordance with current international sites of cave development during the initial standards for cave cartography and mapping reconnaissance. Closed depressions were established by the National Speleological defined as regions that were lower than the Society (NSS) (Dasher, 1994) and the surrounding topography on all sides. Only Association for Mexican Cave Studies closed depressions greater than 10 meters in (AMCS) (Sprouse and Russell, 1980). diameter could be identified based on the Individual surveys were conducted using a limitations of the DEM grid size. The Suunto compass, Suunto inclinometer, and closed depressions where analyzed in fiberglass tape, in association with a field
18
sketch recorded by experienced project on Aguijan and Tinian, 2) analysis of cave sketchers. Caves and other karst features, maps and delineation of primary cave including discharge and recharge features, orientations and cave segment orientations which did not warrant survey, were photo- for features on Tinian, 3) analysis of faults documented and recorded. Discharge and joints reported by Doan and coworkers volumes on discharge features were (1960), 4) analysis of scarp and coastline estimated (minimal discharge and significant orientations on Tinian, and 5) orientation discharge), because coastal conditions did data reduction and production of rose not allow for measurements of salinity and diagrams for structural and cave orientations temperature to be taken at most sites. on Tinian. During fieldwork, features were Data from surveyed cave and karst classified by cave type based on their features was reduced using the software appearance in accordance with the package WALLS (McKenzie, 2002), which Carbonate Island Karst Model. The feature enables compass, inclinometer and tape types included: banana hole, discharge measurements taken during surveys to be cave/feature, fissure cave, flank margin plotted with corrections for minor loop cave, recharge cave/feature, and pit cave. closure errors and regional magnetic When satellite coverage permitted, declination. The corrected line plots for Universal Transverse Mercator System each cave were then used as a basis for (UTM) coordinates and elevation were drafting spatially correct final maps of the recorded with the global positioning system features using Corel Xara 2.0 (Xara, 1997). (GPS), in order to establish accurate location Final cartographic products were created by information for the feature. When satellite overlaying field sketch notes and corrected coverage was not possible, because of line plots. The field sketch notes were vegetation or topography, cave locations scaled and rubber-banded to match the were identified on 1:25,000 USGS corrected line plots, then features recorded topographic maps (USGS, 1983), which on the field sketch notes were manually were then used to determine UTM digitized using standard cave symbology coordinates and elevation for features. established by the NSS and AMCS (Figure Structural orientations of planes of 21, Appendix A; Sprouse and Russell, 1980) brittle failure were measured using a Suunto in order to create accurate maps of the cave compass during fieldwork on Tinian. and karst features inventoried. Aguijan was excluded from this phase of These maps were used to delineate data collection, because no published primary cave orientations and segment geologic map for the island exists and orientations for features surveyed on Tinian, available time on the island was limited. using two methods: 1) apparent trends and Orientations were taken in areas where a 2) entrance oriented trends. Cave joint or fracture could be observed in the orientations have been used in previous bedrock and was not obscured by karren or studies in continental settings to correlate phytokarst development. Orientations were cave and karst development with regional only taken where the joint or fracture structure (Nelson, 1988; Barlow and Ogden, continued over a distance of several meters, 1982), where the orientations of individual cut through more than one bedding plane or cave segments that have a consistent trend several meters of bedrock, and in areas are measured and compared against regional where the bedrock appeared to be in situ. brittle failure features to determine if a correlation between the two populations DATA REDUCTION exist. Data reduction included five phases: 1) Using the apparent trend method, all survey data reduction and production of Tinian cave maps were analyzed maps of the cave and karst features surveyed individually and a primary axis was defined
19
through the cave based on the maximum 1988; Barlow and Ogden, 1982). Two length of the cave and the location of the segment lengths were used to form two sets breached entrance. Cave segments were of data that would reduce the subjectivity of then defined by the orientation of individual the parameters that were used to define chambers, passages, and wall characteristics segment length, such that smaller cave such as large pockets and chamber alcoves sections would be included in the 5-meter (Figures 9 and 10). This technique is highly segment data, while only larger cave subjective, but was used because of its sections would be included in the 10-meter similarity to studies in caves where passages segment data. tend to be linear (Nelson, 1988; Barlow and The entrance width trend method of Ogden, 1982). The primary and segment delineating cave orientations was developed cave trends were then measured and and performed in an attempt to reduce the compiled. The segment cave trends were subjectivity in data reduction for flank length-weighted into 5 and 10-meter margin and banana hole type caves, which increments, such that each 5 or 10-meter tend to form globular or elliptical chambers section of cave segment was counted as an instead of the more linear passages seen in individual orientation measurement in order fissure caves. Using the entrance orientation to give greater significance to longer method, up to three orientations were segment trends during data analysis (Nelson, measured for flank margin and banana hole type cave maps on Tinian (Figure 11). In this technique, the orientation of the entrances to caves that were entered horizontally was measured. Next a maximum penetration measurement was calculated near perpendicular to the entrance (90o ± 15o from the entrance orientation). In caves that had been breached by ceiling Figure 9: Example of primary and segment collapse and were entered vertically, no orientation trends measured using the entrance orientation was measured, but the apparent trend method for a typical penetration measurement was defined as the fissure cave. longest dimension of the cave. Based on the
Figure 10: Example of primary and segment Figure 11: Example of entrance, maximum width, orientation trends measured using and penetration trends measured using the apparent trend method for a the entrance width trend method for a typical mixing zone cave. typical mixing zone cave.
20
penetration measurement, a maximum width package GEOrient 9.2 (Holcombe, 2003), was measured near perpendicular to the using five-degree orientation sectors in order penetration orientation (90o ± 15o from the to be able to visually compare the penetration orientation). In cases where the orientation patterns of caves, zones of brittle cave entrance was the maximum width, the failure, scarps, and coastlines. Rose same measurement was reported for both diagrams were produced for each type of entrance width and the maximum width. orientation data in each of the five These measurements were then compiled as physiographic provinces of Tinian (Central primary entrance width, penetration length, Plateau, Median Valley, Northern Lowland, and maximum width orientations and North-Central Plateau, and Southeastern corresponding segment orientations that Ridge) and for the entire island of Tinian were length-weighted in 5 and 10-meter (Tinian Composite). These rose diagrams segments similar to the cave segments include the following types: defined in the apparent method. Faults and joints reported by Doan and 1. fault orientations coworkers (1960) were measured for length 2. fault, 5-meter segment orientations and orientation. The data for faults was 3. fault, 10-meter segment orientations length-weighted by 50 and 100-meter 4. joint orientations segments in order to apply greater 5. orientations of fractures measured during fieldwork significance to faults that extended over 6. inland scarp orientations greater distances. Two segment lengths 7. inland scarp, 5-meter segment orientations were used in order to reduce the subjectivity 8. inland scarp, 10-meter segment in selecting segment lengths, such that orientations shorter lengths would be included in the 50- 9. coastal scarp orientations meter segment data and only longer 10. coastal scarp, 5-meter segment segments would be included in the 100- orientations meter segment data. These two segment 11. coastal scarp,10-meter segment lengths were defined because they are orientations proportional to the segment definition used 12. all scarps orientations 13. all scarps, 5-meter segment orientations for cave passages, but measured at a scale 14. all scarps, 10-meter segment orientations that is one order of magnitude greater, 15. coastline orientations because 5 and 10-meter lengths could not be 16. coastline, 5-meter segment orientations accurately measured from the geology 17. coastline, 10-meter segment orientations reported by Doan and coworkers (1960). 18. fissure cave orientations Similarly, all scarps that were identified 19. fissure cave, 5-meter segment orientations during the initial site investigation and all 20. fissure cave, 10-meter segment coastlines were divided into linear segments orientations in order to eliminate the effects of minor 21. mixing zone cave orientations variations in coastal erosion and mass 22. mixing zone cave, 5-meter segment orientations wasting. The linear segments were 23. mixing zone cave, 10-meter segment measured for orientation and length; with orientations length weighted, 50 and 100-meter segments 24. all cave types orientations applied to both scarp and coastline 25. all cave types, 5-meter segment orientations as was defined for fault orientations segments, because the DEM cell size does 26. all cave types, 10-meter segment not allow for accurate measurements less orientations than 10 meters. 27. mixing zone cave penetration orientations In the final phase of data reduction, all 28. mixing zone cave penetration, 5-meter orientation measurements were reduced to segment orientations 29. mixing zone cave penetration, 10-meter produce rose diagrams using the software segment orientations
21
30. mixing zone cave entrance width 7. inland scarp, 5-meter segment orientations orientations 31. mixing zone cave entrance width, 5-meter 8. inland scarp, 10-meter segment segment orientations orientations 32. mixing zone cave entrance width, 10- 9. coastal scarp orientations meter segment orientations 10. coastal scarp, 5-meter segment 33. mixing zone cave maximum width orientations orientations 11. coastal scarp, 10-meter segment 34. mixing zone cave maximum width, 5- orientations meter segment orientations 12. all scarps orientations 35. mixing zone cave maximum width, 10- 13. all scarps, 5-meter segment orientations meter segment orientations 14. all scarps, 10-meter segment orientations STATISTICAL COMPARSION 15. coastline orientations OF DATA 16. coastline, 5-meter segment orientations Statistical comparisons were calculated 17. coastline, 10-meter segment in order to determine if different populations orientations of orientation data were similar within each These correlations were performed in of the five physiographic provinces of order to test three different null hypotheses Tinian and the entire island of Tinian. The (H ): comparisons were calculated in order to 0 determine if there is a significant H (1): Regional brittle deformation and relationship between: 1) brittle deformation 0 karst development represent of eogenetic rocks and the development of significantly different populations eogenetic karst; 2) scarp orientations and the or population distributions and development of eogenetic karst; and 3) develop independently. coastline orientations and the development of eogenetic karst. H (2): Regional scarp positions and karst Statistical comparisons were conducted 0 development represent using the Kolmogorov-Smirnov 2-sample significantly different populations test (Burt and Barber, 1996; Till, 1974; or population distributions and Miller and Kahn, 1962), which performs a develop independently. non-parametric test on two independent sample populations to determine if they H (3): Regional coastline positions and represent similar populations or population 0 karst development represent distributions. During the statistical significantly different populations comparison phase, several populations of or population distributions and orientations were evaluated, which develop independently. correspond to the 35 orientation types that were used to produce rose diagrams during The first null hypothesis (H (1)) was the data reduction phase. Each of the 35 0 formulated in order to determine if there is a orientation types was evaluated for relationship between brittle deformation and correlation with: karst development on Tinian, while the second and third null hypotheses (H0 (2 and 1. fault orientations 3)) were formulated to determine if there is a 2. fault, 5-meter segment orientations relationship between karst development and 3. fault, 10-meter segment orientations coastlines or scarps, where the edge of the 4. joint orientations paleo fresh-water lens would be expected 5. orientations of fractures measured based on the Carbonate Island Karst Model during fieldwork (Mylroie and Vacher, 1999). Because of the 6. inland scarp orientations wide variation in orientations of the data in
22
this study, a high significance level (P ≤ province scale, and 3) site scale. Because of 0.01) was used to ensure that the populations the wide distribution of cave and karst or population distributions that were features, one square kilometer test sites were compared represented significantly similar chosen for additional analyses. In order to datasets. If the first null hypothesis is compare the results of analyses from test rejected, then the regional populations of sites with the larger regions of Tinian, five joints, fractures, and faults are similar to the parameters where required to exist within regional karst development and a positive the one square kilometer boundary, relationship between the two populations including: 1) three or more surveyed can be inferred. If the second null features, 2) both fissure and mixing zone hypothesis is rejected, then the regional caves, 3) coastline, 4) faults reported by populations of scarp orientations are similar Doan and coworkers (1960), and 5) fracture to the regional karst development and a orientations measured during fieldwork. positive relationship between the two The requirements for choosing small- populations can be inferred. If the third null scale test sites required the completion of hypothesis is rejected, then the regional fieldwork and the cave and karst inventory populations of coastline orientations are of this study. Once the inventory stage was similar to the regional karst development complete, three test sites were chosen for and a positive relationship between the two analysis (Figure 22, Appendix A): Carolinas populations can be inferred. Limestone Forest, Puntan Diapblo, and Unai Dangkolo. At each of these sites, statistical SMALL-SCALE TEST comparisons were evaluated using the SITE EVALUATION Kolmogorov-Smirnov 2-sample test, with In addition to the analyses performed for the same 35 parameters evaluated for each the five physiographic provinces and the of the physiographic provinces and the island of Tinian, smaller-scale test sites were entire island of Tinian. analyzed to determine if any relationships existed at various scales: 1) island scale, 2)
STUDY RESULTS
CAVE AND KARST INVENTORY zone type caves (banana hole and flank The mapping portion of cave and karst margin type caves) were the most prominent inventory surveyed 114 caves or cave cave types found during the inventory (81% complexes (green areas, Figure 23, on Aguijan; 77% on Tinian). While Tinian Appendix A; Appendix B): 26 on Aguijan exhibited a larger diversity and quantity of and 88 on Tinian. The features were caves, the ratio of fissure caves to mixing classified by morphology type and grouped zone caves on both islands was 0.24 (0.235 by physiographic province they were found for Aguijan; 0.238 for Tinian), with fissure in (Table 1, Appendix B; Table 2, Appendix caves accounting for almost 20% of the B). On Aguijan, the cave survey produced caves recorded. maps of 1 banana hole, 5 fissure caves, and A total of seventeen locations were 20 flank margin caves. On Tinian, the cave found on Tinian, which exhibited fresh- survey produced maps of 3 banana holes, 5 water discharge (Figure 24, Appendix A). discharge features, 12 fissure caves, 65 flank Most of these features were identified by margin caves, 1 pit cave, and 2 recharge schlieren mixing at the coastline, but only features; however, 4 of the discharge five of them were associated with cave features were also classified as fissure caves, development and surveyed. Fresh-water making a total of 16 fissure caves. Mixing discharge observed at Unai Dangkolo, Unai
23
Masalok and Taga Beach appeared as fresh-water at the time of survey (Figure 29, diffuse discharge through carbonate sand Appendix A), while the largest closed beach deposits, while the discharge at other depression on the island (Sisonyan Makpo) sites appeared as focused discharge along is less than 3 meters above mean sea-level at fractures and bedding planes. The two its lowest elevation and is the site of the discharge features that exhibited the greatest islands primary municipal well (Makpo discharge were located on the east and west Wells). The three unconfirmed recharge coasts of the island at Gecko Cave and features are located north of Mount Lasu Barcinas Cove, respectively. near small outcrops of non-carbonate rocks, Investigation of closed depressions but due to their small size and dense identified four definite regions of allogenic vegetation in the area, a positive recharge in the North-Central Highland, confirmation could not be made. with small caves associated with two of them (Lasu Recharge Cave and West Lasu CAVE ORIENTATIONS Depression Cave). The presence of The two methods of cave orientation accumulated detritus and the lack of analysis (apparent trend and entrance width sediment coatings on the feature walls is trend) were applied to orientation data for evidence that recharge is rapid and that mapped caves on Tinian. Aguijan was water does not pond at these features during excluded from this analysis because no recharge events. The other two recharge geologic map has been published for the features (North and South Lemmai Recharge area, which would compliment the cave Features) covered larger areas, contained orientation data when comparing cave and contacts between carbonate and non- karst development to brittle failure features. carbonate rocks, and showed vegetative Each cave was analyzed separately, evidence of water ponding within the closed including those that occurred on maps where depressions prior to entering the subsurface a complex or series of caves were surveyed (Figure 25, Appendix A). Ponded water was together as one map, making it possible to present in portions of South Lemmai have several datasets of orientation Recharge Feature when investigated (Figure measurements for some cave maps. Using 26, Appendix A). the apparent trend method all fissure caves The other sixteen closed depressions and mixing zone caves (flank margin and identified during the initial site investigation banana hole type caves) were analyzed for that were not confirmed as allogenic primary orientation and segment recharge features were also investigated orientations. The single pit cave and two (Figure 27, Appendix A), either physically recharge caves were excluded from all or through communication with local datasets except composite cave types (all residents. Five closed depressions were cave types) because of their small sample identified as quarries or borrow pits, eight size. Using the entrance trend method, the were identified as natural constructional mixing zone caves were analyzed for features, and three appear to be recharge entrance width, penetration length and features. Four of the modified closed maximum width, while fissure caves, the depressions (quarries and borrow pits) single pit cave and the two recharge caves showed evidence of excavation in the past, were excluded. The fissure caves were while the quarry located near Barcinas Cove excluded from the entrance width trend is being actively excavated today (Figure 28, method because they are linear features with Appendix A). The eight natural distinct segment orientations similar to construction features showed no evidence of telogenetic caves (Nelson, 1988; Barlow and excavation or allogenic recharge; however, Ogden, 1982), while the single pit cave and Hagoi in the Northern Lowland is less than 2 two recharge caves were excluded again meters above mean sea-level and contained because of their small sample size.
24
The apparent trend method was applied Southeastern Ridge: 115 orientations, 213 because of its similarity to previous studies five-meter segments and 128 ten-meter (Nelson, 1988; Barlow and Ogden, 1982). Apparent trend analysis of fissure caves The entrance width trend method was analysis yielded 18 primary cave performed in order to reduce the subjectivity orientations (Table 3, Appendix C): 5 in the of orientation measurements that exist in the Central Plateau; 5 in the Median Valley; 0 in apparent trend method. Measurements used the North-Central Highland; 1 in the in the entrance width trend method included Northern Lowland; and 7 in the entrance width, penetration length, and Southeastern Ridge. Fissure cave segment maximum width, with 10 caves excluded analysis using the apparent trend method from the entrance width analysis because yielded 147 orientations with 297 five-meter their entrances were entered vertically as a segments and 135 ten-meter segments result of ceiling collapse. Entrance width (Table 4, Appendix C): analysis using the entrance trend method yielded 107 entrance widths with 324 five- Central Plateau: 15 orientations, 24 five- meter segments and 163 ten-meter segments meter segments and 12 ten-meter (Table 7, Appendix C): segments Median Valley: 69 orientations, 100 five- Central Plateau: 40 orientations, 105 five- meter segments and 44 ten-meter meter segments and 59 ten-meter segments segments North-Central Highland: 0 orientations Median Valley: 20 orientations, 89 five- and segments meter segments and 43 ten-meter Northern Lowland: 2 orientations, 2 five- segments meter segments and 2 ten-meter segments North-Central Highland: 5 orientations, Southeastern Ridge: 61 orientations, 171 11 five-meter segments and 6 ten-meter five-meter segments and 77 ten-meter segments segments Northern Lowland: 1 orientation, 22 five- meter segments and 11 ten-meter Apparent trend analysis of mixing zone segments caves yielded 128 primary orientations Southeastern Ridge: 42 orientations, 97 (Table 5, Appendix C): 46 in the Central five-meter segments and 44 ten-meter Plateau; 27 in the Median Valley; 8 in the segments North-Central Highland; 1 in the Northern Lowland; and 46 in the Southeastern Ridge. Penetration length analysis using the Mixing zone cave segment analysis yielded entrance trend method produced 118 388 orientations with 980 five-meter orientations with 369 five-meter segments segments and 480 ten-meter segments and 190 ten-meter segments (Table 8, (Table 6, Appendix C): Appendix C):