Water Quality and Landscape Processes of Four Watersheds in Eastern Puerto Rico Professional Paper 1789 U.S. Department of the Interior U.S. Geological Survey Cover : Main channels of four study rivers in eastern Puerto Rico. The top row shows the rivers with naturally forested watersheds (left, Río Mameyes; right, Río Icacos), and the bottom row shows the rivers with agriculturally developed watersheds (left, Río Canóvanas; right, Río Cayaguás). Rivers that drain volcanic and volcaniclastic bedrock with little to no quartz are on the left, and rivers that drain quartz-rich granitic bedrock are on the right. Whether natural or developed, channels of rivers draining bedrock with little to no quartz contain mostly cobbles and boulders, and channels of rivers draining quartz-rich bedrock contain quartz-rich sand. Water Quality and Landscape Processes of Four Watersheds in Eastern Puerto Rico Edited by Sheila F. Murphy and Robert F. Stallard Contributions by Heather L. Buss, William A. Gould,1 Matthew C. Larsen, Zhigang Liu,2 Sebastián Martinuzzi,3 Sheila F. Murphy, Robert F. Stallard, Isabel K. Parés-Ramos,4 Arthur F. White, and Xiaoming Zou2,5 1 U.S. Department of Agriculture Forest Service, International Institute of Tropical Forestry 2 Institute for Tropical Ecosystem Studies, University of Puerto Rico 3 Department of Forest and Wildlife Ecology, University of Wisconsin 4 Department of Biology, University of Puerto Rico 5 XishuangbannaTropical Botanical Garden, Chinese Academy of Sciences Volume comprises Chapters A, B, C, D, E, F, G, H, and I Professional Paper 1789 U.S. Department of the Interior U.S. Geological Survey U.S. Department of the Interior KEN SALAZAR, Secretary U.S. Geological Survey Marcia K. McNutt, Director U.S. Geological Survey, Reston, Virginia: 2012 For more information on the USGS—the Federal source for science about the Earth, its natural and living resources, natural hazards, and the environment, visit http://www.usgs.gov or call 1–888–ASK–USGS. For an overview of USGS information products, including maps, imagery, and publications, visit http://www.usgs.gov/pubprod To order this and other USGS information products, visit http://store.usgs.gov Any use of trade, firm, or product names is for descriptive purposes only and does not imply endorsement by the U.S. Government. Although this information product, for the most part, is in the public domain, it also may contain copyrighted materials as noted in the text. Permission to reproduce copyrighted items must be secured from the copyright owner. Suggested citation: Murphy, S.F., and Stallard, R.F., eds., 2012, Water quality and landscape processes of four watersheds in eastern Puerto Rico: U.S. Geological Survey Professional Paper 1789, 292 p. iii Volume Contents Foreword ...............................................................................................................................................iv By Matthew C. Larsen Executive Summary ..............................................................................................................................v By Sheila F. Murphy and Robert F. Stallard A. Physiography, Geology, and Land Cover of Four Watersheds in Eastern Puerto Rico ...............................................................................................................................1 By Sheila F. Murphy, Robert F. Stallard, Matthew C. Larsen, and William A. Gould B. Land Use, Population Dynamics, and Land-Cover Change in Eastern Puerto Rico .............................................................................................................................25 By William A. Gould, Sebastián Martinuzzi, and Isabel K. Parés-Ramos C. Hydrology and Climate of Four Watersheds in Eastern Puerto Rico ..................................43 By Sheila F. Murphy and Robert F. Stallard D. Atmospheric Inputs to Watersheds of the Luquillo Mountains in Eastern Puerto Rico .............................................................................................................................85 By Robert F. Stallard E. Water Quality and Mass Transport in Four Watersheds in Eastern Puerto Rico ...........113 By Robert F. Stallard and Sheila F. Murphy F. Landslides and Sediment Budgets in Four Watersheds in Eastern Puerto Rico ............153 By Matthew C. Larsen G. Effects of Earthworms on Slopewash, Surface Runoff, and Fine-Litter Transport on a Humid-Tropical Forested Hillslope in Eastern Puerto Rico..................................179 By Matthew C. Larsen, Zhigang Liu, and Xiaoming Zou H. Weathering, Landscape Equilibrium, and Carbon in Four Watersheds in Eastern Puerto Rico ............................................................................................................199 By Robert F. Stallard I. Weathering Processes in the Icacos and Mameyes Watersheds in Eastern Puerto Rico ...........................................................................................................................249 By Heather L. Buss and Arthur F. White Appendix 1. Data Processing and Computation to Characterize Hydrology and Compare Water Quality of Four Watersheds in Eastern Puerto Rico .........................263 Appendix 2. Methods Used to Analyze Water Quality of Four Watersheds in Eastern Puerto Rico ...........................................................................................................................289 iv Foreword By Matthew C. Larsen Humans face a future limited by the increasing effects of environmental problems and resource scarcity. The world’s population is now consuming water and energy, and using agricultural and other land resources, at per capita rates that will be challenging to sustain. Environmental degradation is intensifying conflict over resources. Local consequences are exacerbated by climate change, manifested as changes in the timing and distribution of rainfall and runoff. As our influence on Earth’s lands and waters grows more substantial, we need to understand and quantify these effects to preserve our well-being. Earth sciences are central to meeting this need because they guide the monitoring, assessment, and modeling of natural processes that help us to understand our interaction with the environment. The U.S. Geological Survey has a long history of providing authoritative earth science to inform local and national decision makers as they confront difficult choices in a changing world. The low-latitude regions of Earth, in particular, are undergoing profound, rapid landscape change. Humid-tropical watersheds of eastern Puerto Rico, at 18° North latitude, provide a natu- ral laboratory that allows researchers to quantify and better understand this change and how it impinges on the human world. In recognition of this opportunity, the National Science Founda- tion’s Long-Term Ecological Research program and the U.S. Geological Survey’s Water, Energy, and Biogeochemical Budgets (WEBB) program have operated in eastern Puerto Rico for several decades. WEBB program scientists have studied hydrologic, geologic, geomorphic, and anthro- pogenic processes in eastern Puerto Rico in an effort to generalize their results to the broader tropical world. Long time-series of streamflow, water-quality, suspended-sediment, meteorologi- cal, and land-use data, collected by Federal and by Commonwealth of Puerto Rico agencies, record the effects of landscape changes that are characteristic of environmental challenges that we all face. This publication summarizes research based on these data, and it provides insights into how landscape change may increase soil erosion, reservoir sedimentation, and landslide hazards while decreasing soil-nutrient content and surface-water quality. In the 1990s the Puerto Rico Tourism Company invited travelers to “discover the continent of Puerto Rico.” In a very different sense, Puerto Rico’s landscape and the results of our studies there provide a microcosm of the challenges faced by societies on islands and continental land- masses. What choices will we make as we confront future limitations of land and resources? What scientific information will we need to guide these choices? The science presented in this volume improves our understanding of the ramifications of landscape changes by using analy- sis of watershed-scale geologic, hydrologic, and geomorphic processes, and it can help guide policymakers and resource managers as they confront difficult choices. Executive Summary By Sheila F. Murphy and Robert F. Stallard Puerto Rico is in a state of rapid, ongoing change. Locally, agricultural lands are undergoing reforestation, while coastal areas are becoming heavily urbanized. The area is also changing because of the introduction of nonnative species, water supply projects, and the construction of roads and other infrastructure. Superimposed on these local phenomena are slower, larger scale changes, such as the deposition of airborne pollutants and natural and human-induced climate change. Owing to the island’s steep topography, low water-storage capacity, and dependence on trade-wind precipitation, Puerto Rico’s people, ecosystems, and water supply are vulnerable to extreme weather such as hurricanes, floods, and droughts. Eastern Puerto Rico offers a natural laboratory for separating geologic and land-cover influences from regional- and global-scale influences because of its various bedrock types and the changing land cover surrounding intact, mature forest of the Luquillo Experimental Forest (which is contiguous with El Yunque National Forest). Accordingly, a multiyear assessment of hydrological and biogeochemical processes was designed to develop an understanding of the effects of these differences on local climate, streamflow, water quality, and ecosystems, and to form the basis for a long-term and event- based program of climate and hydrologic monitoring. Focusing on small watersheds allows for integrated studies of hydrologic and chemical pro- cesses, owing to minimal spatial variability of geology, land use, and climate. For two decades, the U.S. Geological Survey has been evaluating the processes controlling fluxes of water, energy, and elements throughout a range of temporal and spatial scales in small watersheds at five sites in different parts of the nation. The Water, Energy, and Biogeochemical Budgets Proj- ect in eastern Puerto Rico represents a montane, humid-tropical environment, in which lie four watersheds of differing geology and land use. Two watersheds are located on coarse-grained granitic rocks (Icacos and Cayaguás), and two are located on fine-grained volcanic rocks and vol- caniclastic sedimentary rocks (Mameyes and Canóvanas). For each bedrock type, one watershed is covered with mature rainforest (Icacos and Mameyes); the other is undergoing reforestation after being used as agricultural land (Cayaguás and Canóvanas). These watersheds, like most of the rest of Puerto Rico, were subjected to intensive agriculture in the 19th and early 20th century, but they have been undergoing reforestation as a result of a shift from an agricultural economy to an industrial one and subsequent human migration to urban areas (discussed in chapters A and B in this volume). Puerto Rico lies directly in the path of the easterly trade winds, which deliver steady rainfall to the mountains and steer weather systems called tropical waves toward the island. Hurricanes and tropical storms derived from these systems typically deliver the majority of yearly rainfall and occur from May to December (chapter C). Northern cold fronts can also deliver heavy rainfall for several days at a time, typically from December through February. These storms vary greatly in frequency and intensity, contributing to substantial interannual variation in precipitation and stream discharge. The largest storms can have profound geomorphic consequences, such as landslides, debris flows, deep gullying on deforested lands, excavation and suspension of sedi- ment in stream channels, and delivery of a substantial fraction of annual stream sediment load (chapters F and G). Past deforestation and agricultural activities in the Cayaguás and Canóvanas watersheds led to profoundly accelerated erosion and soil loss, and this material continues to be remobilized during large storms. v vi Rainfall varies greatly over small distances in eastern Puerto Rico owing to differences in elevation, topographic position, aspect, and proximity to the ocean. The Icacos and Mameyes watersheds, located on the eastern side of the Luquillo Mountains, are the wettest of the four watersheds studied, and their highest elevations receive more than 4,000 millimeters of rain annually (chapter C). Precipitation increases with elevation in these watersheds. The Canóvanas and Cayaguás watersheds, located on the western side of the Luquillo Mountains, are consid- erably drier, and precipitation and elevation are not correlated. Precipitation and runoff in all watersheds show large interannual variation and are highest in years when major storms—such as Hurricanes Hortense (1996) and Georges (1998)—strike eastern Puerto Rico. These large storms typically produce similar runoff in all of the study watersheds, suggesting that higher annual runoff in the eastern, windward side of the Luquillo Mountains (which includes the Mameyes and Icacos watersheds) is caused by smaller, more-routine rain events. When one considers watershed-wide water budgets, the windward or leeward aspect of a watershed is more important than differing geology and land cover. Regional weather patterns and consequent sources of air masses influence the type and timing of atmospheric contributions to eastern Puerto Rico. Nitrogen loads in precipitation at a National Atmospheric Deposition Project site in eastern Puerto Rico have roughly doubled since measure- ments began in 1985 (chapter D). Eastern Puerto Rico also receives marine salts and Saharan Desert dust in rainfall. The proportion of material delivered by these sources varies season- ally; deposition of marine salts is greatest in January, whereas material from North America is deposited primarily in January, April, and May, and Saharan Desert dust peaks in June and July. Saharan dust typically contributes enough alkalinity in June and July to neutralize acidity in precipitation. During large storms, entrainment of ocean water can lead to highly elevated concentrations of chloride in stream waters (chapter E). Because infrequent, large storms play a major role in erosion of landscapes and can lead to a change in hydrologic flow paths, the Water, Energy and Biogeochemical Budgets Project focused on high-runoff events; we sampled 263 storms, including all major hurricanes that occurred between 1991 and 2005 (chapter C, appendix 1). Nearly 5,000 routine and event samples were analyzed for parameters that allow determination of denudation rates based on suspended and dissolved loads; 860 of these samples were analyzed for a comprehensive suite of chemical constituents. Of samples analyzed for comprehensive chemistry and for sediment, 543 were collected at runoff rates greater than 1 millimeter per hour, 256 at rates exceeding 10 millimeters per hour, and 3 at rates exceeding 90 millimeters per hour. Streams have rarely been sampled during events with such high runoff rates. The rivers studied are generally similar in water-quality characteristics. Most chemical con- stituents show similar trends in the four watersheds, which imply considerable similarity in runoff generation and flow-path structuring despite differences in geology, soils, land cover, and weathering styles (chapter E). The rivers with lowest mean-annual runoff rates and highest ratios of evapotranspiration to runoff (Cayaguás and Canóvanas) tend to have higher concentra- tions of nonbioactive constituents. These developed watersheds typically have higher concentra- tions of nutrients (potassium, nitrate, ammonium ion, phosphate), perhaps indicating additional agricultural or wastewater sources. Projecting watershed yields to a common, intermediate mean-annual runoff (1,860 millimeters per year) generally decreased or did not change the range of yields of constituents that are the primary indicators of chemical weathering, biological activ- ity on the landscape, or atmospheric contributions (dissolved bedrock, sodium, silica, chloride, dissolved organic carbon, and calcium), further indicating no dominant influence of either vii geology or land cover (chapter H). Magnesium and inorganic carbon showed a dependence on geology, possibly due to the presence of carbonates or mafic rocks. Projected yields of nutrients and particulate constituents (suspended solids and particulate organic carbon), however, were far in excess of equilibrium yields, and they were much greater for developed landscapes as compared with forested watersheds, consistent with the known effects of land clearing, agricul- tural activities, and domestic wastewater inputs. Physical and chemical weathering rates of the four watersheds studied are high. Bedrock in the Icacos (granitic rock) and Mameyes (volcaniclastic rock) watersheds have some of the highest documented rates of chemical weathering of silicate rocks in the world (chapter I). Physical denudation rates based on mass balances are higher than expected for a steady-state system; this excess is substantial in all watersheds except the Mameyes (forested watershed on volca- niclastic bedrock; chapter H). Deforestation and agriculture can explain the accelerated physical erosion in the two developed watersheds (Canóvanas and Cayaguás). Physical erosion rates of the granitic watersheds are seven-fold as great as those for the volcaniclastic watersheds, owing to greater permeability and thus higher rates of water filtration and greater susceptibility to landsliding (chapters C and F). The reason for such high rates of physical erosion in the Icacos watershed (forested watershed on granitic bedrock) is unclear but may be related to changes in forest quality or to the history of road construction. The elevated physical erosion drives an increased particulate organic carbon flux, one that is large and is important to the carbon cycle. This increased flux is also sustainable because soil-carbon replacement is rapid. It is crucial to understand long-term geomorphical, hydrological, and biogeochemical processes in tropical regions, because these regions occupy about a quarter of Earth’s land surface, yet they contribute a substantially higher fraction of the water, solutes, and sediment discharged to the world’s oceans. Nearly half of Earth’s population lives in the tropics, and therefore development stresses are intense and can potentially harm soil resources, water quality, and water supply and in addition increase landslide and flood hazards. Small watersheds in eastern Puerto Rico provide an excellent opportunity to examine these processes and their connection to climate, geology, and land cover. The 15-year Water, Energy, and Biogeochemical Budget dataset, which includes discharge, field parameters, suspended sediment, major cations and anions, and nutrients, is available from the U.S. Geological Survey’s National Water Information System (http://waterdata.usgs.gov/nwis). The dataset provides a baseline for characterizing future environmental change and will improve our understanding of the interdependencies of land, water, and biological resources and their responses to changes in climate and land use. Because eastern Puerto Rico resembles many tropical regions in terms of geology and patterns of development, implications from this study are transferable to other tropical regions facing deforestation, rapid land-use change, and climate change. Acknowledgments The work in Puerto Rico could not have been done without the assistance of numerous people. In alphabetical order, an incomplete list includes the following: George Aiken, Ellen Axtmann, Sean Baran, Susan Brantley, Mary-Margaret Coates, Paul Collar, Iris Concepción, Russ Curtis, Pedro Díaz, Marcelle Fabregas, Robert Hirsch, Laura Hubbard, Ariel Lugo, Deborah Martin, Joel Martinez, Erik Oerter, Jennifer Riggs, Manuel Rosario, Cathy Rubin, Abdul Wahab Sadeqi, F.N. Scatena, Jamie Shanley, Angel Torres-Sánchez, Heriberto Torres, Joe Troester, Ank Webbers, and viii Allen Zack. Helpful reviews of individual chapters and appendices were provided by Alex Blum, W. Troy Baisden, Brian Ebel, Allen Gellis, Jonathan Godt, G. González, Cliff R. Hupp, Andy Kurtz, Ariel Lugo, William McDowell, John Moody, Jake Peters, Fred Quiñones, Deborah Repert, Susan Riggins, F.N. Scatena, Martha Scholl, Jamie Shanley, Jonathan Stock, Gordon Tribble, Frank Wadsworth, John Walker, and Richard Webb. Two individuals deserve special acknowledgment: Ellen Axtmann helped implement and ran chemical analyses for the first half of the program, and Angel Torres-Sánchez, who was the first technician, was crucial in field-testing procedures and in helping adapt to the complexities of event sampling in a challenging environment. Physiography, Geology, and Land Cover of Four Watersheds in Eastern Puerto Rico By Sheila F. Murphy, Robert F. Stallard, Matthew C. Larsen, and William A. Gould Chapter A of Water Quality and Landscape Processes of Four Watersheds in Eastern Puerto Rico Edited by Sheila F. Murphy and Robert F. Stallard Professional Paper 1789–A U.S. Department of the Interior U.S. Geological Survey Contents Abstract ...........................................................................................................................................................5 Introduction.....................................................................................................................................................5 Site Description ..............................................................................................................................................5 Physiography .........................................................................................................................................7 Geology .................................................................................................................................................10 Geology in Watersheds Dominated by Volcaniclastic Rock ..............................................10 Geology in Watersheds Dominated by Granitic Rock ..........................................................12 Soil and Saprolite ................................................................................................................................13 Soils in Watersheds Dominated by Volcaniclastic Rock ....................................................13 Soils in Watersheds Dominated by Granitic Rock ................................................................14 Vegetation and Land Use ...................................................................................................................15 Land Cover of Forested Watersheds ......................................................................................16 Land Cover of Developed Watersheds ...................................................................................18 Summary........................................................................................................................................................18 Acknowledgments .......................................................................................................................................20 References ....................................................................................................................................................20 Figures 1. Maps of Puerto Rico and study watersheds ............................................................................6 2. Digital elevation model of elevation and relief of Puerto Rico, showing major mountain ranges and study watersheds ..................................................................................7 3. Digital elevation model of relief of eastern Puerto Rico showing study watersheds, and elevation profiles along three lines of cross section ......................................................8 4. Diagrams of percentage of watershed area with each aspect, for the Canóvanas, Cayaguás, Icacos, and Mameyes watersheds ........................................................................9 5–9. Maps of the following: 5. Generalized geology of Puerto Rico .................................................................................10 6. Detailed geology of study watersheds, Puerto Rico .....................................................11 7. Ecological life zones of Puerto Rico .................................................................................15 8. Generalized land cover of Puerto Rico in 2003 ...............................................................16 9. Detailed land cover of study watersheds in 2003 ..........................................................17 10. Diagram of population density in a barrio in Canóvanas watershed and a barrio in Cayaguás watershed .................................................................................................................20 Tables 1. Geomorphic and geographic characteristics of study watersheds ....................................9 2. Geology of study watersheds, eastern Puerto Rico .............................................................12 3. Mineral content of granitic intrusions, eastern Puerto Rico ...............................................13 4. Land cover and land-cover dynamics in study watersheds, eastern Puerto Rico ..........19 3Abbreviations Used in This Report cm centimeter km kilometer km2 square kilometer m meter mm millimeter Ma mega-annum, 1 million years old USDA U.S. Department of Agriculture USGS U.S. Geological Survey WEBB Water, Energy, and Biogeochemical Budgets Conversion Factors Multiply By To obtain Length centimeter (cm) 0.3937 inch (in.) millimeter (mm) 0.03937 inch (in.) meter (m) 3.281 foot (ft) kilometer (km) 0.6214 mile (mi) Area square kilometer (km2) 0.3861 square mile (mi2) Abstract Four watersheds with differing geology and land cover in eastern Puerto Rico have been studied on a long-term basis by the U.S. Geological Survey to evaluate water, energy, and biogeochemical budgets. These watersheds are typical of tropical, island-arc settings found in many parts of the world. Two watersheds are located on coarse-grained granitic rocks that weather to quartz- and clay-rich, sandy soils, and two are located on fine-grained volcanic rocks and volcaniclastic sedimentary rocks that weather to quartz-poor, fine-grained soils. For each bedrock type, one watershed is covered with mature forest, and the other watershed, like most of Puerto Rico, has transformed from relatively undis- turbed pre-European forest to intensive agriculture in the 19th and early 20th centuries, and further to ongoing refor- estation that began in the middle of the 20th century. The comparison of water chemistry and hydrology in these water- sheds allows an evaluation of the effects of land-use history and geology on hydrologic regimes and erosion rates. This chapter describes the physiography, geology, and land cover of the four watersheds and provides background information for the remaining chapters in this volume. Introduction Five field sites across the nation are being monitored long-term as part of the U.S. Geological Survey’s (USGS) Water, Energy, and Biogeochemical Budgets (WEBB) program (Baedecker and Friedman, 2000). The Puerto Rico WEBB site, located in eastern Puerto Rico, represents a mon- tane, humid-tropical environment. Its geology and land cover are typical of large regions of the humid tropics, particularly montane regions that are affected by hurricanes (Larsen and Stallard, 2000). The Puerto Rico WEBB project has been evaluating hydrological, chemical, and sediment budgets and processes in four watersheds of differing geology and land use in eastern Puerto Rico since 1991 (Larsen and others, 1993). Two watersheds are located on coarse-grained granitic rocks, and two are located primarily on a combination of fine-grained volcanic and volcaniclastic sedimentary rocks. For each bed- rock type, one watershed is covered with mature forest, and the other watershed has been affected by agricultural land use typical of eastern Puerto Rico. The comparison of water chem- istry and hydrology in these watersheds allows an evaluation of the effects of land-use history and geology on hydrologic regimes and erosion rates. This chapter describes and compares the physiography, geology, soils, vegetation, and land cover of the watersheds. Climate and hydrology are described in Murphy and Stallard (2012), land-cover change is discussed in Gould and others (2012), and erosional processes are described in Larsen (2012). These comparisons provide context for data presented in the remaining chapters. Site Description Puerto Rico is a 9,000-square kilometer (km2) island within the Greater Antilles, located about 1,700 km south- east of Miami, Fla., at about latitude 18°N., longitude 66°W. (fig. 1). The WEBB study watersheds discussed here and in other chapters of this professional paper are gaged subwater- sheds of larger watersheds having the same name. The Mam- eyes (17.8 km2) and Icacos (3.26 km2) watersheds are located inside El Yunque National Forest, a 113-km2 forest preserve administered by the U.S. Department of Agriculture (USDA) Forest Service. El Yunque National Forest is conterminous with the Luquillo Experimental Forest (fig. 1), which has been the site of a National Science Foundation Long-Term Ecological Research project since 1985, and it is designated an International Biosphere Reserve by the United Nations Educational, Scientific, and Cultural Organization. The Cay- aguás (26.4 km2) and Canóvanas (25.5 km2) watersheds are located within the partially urbanized, agriculturally devel- oped yet reforesting, Río Grande de Loíza watershed (fig. 1). Physiography, Geology, and Land Cover of Four Watersheds in Eastern Puerto Rico By Sheila F. Murphy,1 Robert F. Stallard,1 Matthew C. Larsen,2 and William A. Gould3 1U.S. Geological Survey, 3215 Marine St. Suite E127, Boulder, Colo. 80303, USA 2U.S. Geological Survey, 436 National Center, Reston, Va. 20192, USA 3U.S. Department of Agriculture Forest Service, International Institute of Tropical Forestry, Jardín Botánico Sur, 1201 Calle Ceiba, Río Piedras, Puerto Rico 00926 6 Water Quality and Landscape Processes of Four Watersheds in Eastern Puerto Rico Area shown in B A B FLORIDA PUERTO RICO 25° 85° 80° 75° 70° 65° 60° 0 0 200 400 KILOMETERS 200 400 MILES 20° 15° 10° BAHAMAS DOMINICAN REPUBLIC HAITI CUBA JAMAICA L e s s e r A n t i l l e s G r e a t e r A n t i l l e s C A R I B B E A N S E A A T L A N T I C O C E A N Cayaguás San Juan 18°20’ 66°10’ 66°00’ 65°50’ 65°40’ 18°10’ A T L A N T I C O C E A N 0 0 10 KILOMETERS 10 MILES 5 5 N:\Jeff\den11_hwcg00_0224_pp_murphy\chapter_A_figures\BW_figures\BWfigure_01A.ai Canóvanas C A Y E Y M O U N T A I N S C A R I B B E A N S E A Mameyes Icacos L U Q U I L L O M O U N T A I N SRío Grande de Loíza Watershed Base from United Nations Digital Soil Map of the World, 2007 1:5,000,000, Geographic coordinate system, degrees World Geodetic Survey 1984 (WGS 1984) Base from U.S. Geological Survey National Atlas, 2002 1:2,000,000, Geographic coordinate system, degrees North American Datum of 1983 (NAD 1983) EXPLANATION Luquillo Experimental Forest Watershed boundary Study watershed boundary Figure 1. Location of Puerto Rico and study watersheds, eastern Puerto Rico. Physiography, Geology, and Land Cover of Four Watersheds in Eastern Puerto Rico 7 Physiography Puerto Rico is mountainous with central highlands surrounded by flat-lying coastal plains and alluvial valleys (fig. 2). Mountain ranges generally are oriented east-west and dominate the southern two-thirds of the island. Stream channels are steep, and most river valleys are deeply incised. Accelerated weathering in the north-facing watersheds, com- pared with the slower weathering on the drier southern slopes, has shifted the divide south such that streams are longer and their slopes are gentler on the north side of the island than on the south (Ramos-Ginés, 1999). The Mameyes, Icacos, and Canóvanas watersheds drain the rugged Luquillo Mountains (fig. 2). The Mameyes water- shed is located on the eastern side of the Luquillo Mountains (fig. 3). Slopes in this watershed mostly face north, northeast, and east (fig. 4), directly into the dominant east and northeast wind directions (Calvesbert, 1970). The Canóvanas watershed is on the western side of the Luquillo Mountains, with south- west to northwest aspects, facing away from prevailing winds. Relief in the Mameyes and Canóvanas watersheds is almost 1,000 meters (m) (table 1). The smaller Icacos watershed is located at elevations between about 620 and 832 m, and slopes typically face east, southeast, and southwest (fig. 4). A ridge approaching 800 m elevation separates the Icacos and Mameyes watersheds. Five peaks are higher than 1,000 m in the Luquillo Mountains (U.S. Geological Survey, 1967); the highest, El Toro (1,075 m), is located between the Icacos and Canóvanas watersheds (fig. 3). The Cayaguás watershed, located southwest of the Luquillo Mountains, consists of low hills and alluvial valleys; the Cayey Mountains are located to the west (fig. 2). Slopes in the Cayaguás watershed show no predominant orientation (fig. 4). Hillslopes and channels are steep in all of the WEBB watersheds, and mean hillslopes range from 0.189 (10.7°) in the Cayaguás watershed to 0.365 (20.0°) in the Mameyes watershed (table 1). Figure 2. Elevation and relief of Puerto Rico, showing major mountain ranges and study watersheds (outlined). N:\Jeff\den11_hwcg00_0224_pp_murphy\chapter_A_figures\BW_figures\BWfigure_02A.ai EXPLANATION Elevation—In meters −29–0 1–200 201–400 401–600 601–800 801–1,000 Base from U.S. Geological Survey National Elevation Dataset 30-meter digital elevation model 0 0 10 20 KILOMETERS 10 20 MILES Study watershed boundary Canóvanas CayaguásC O R D I L L E R A C E N T R A L CAYEY MOUNTAINS LUQUILLO MOUNTAINS Icacos Mameyes 66°30’67°00’ 18°30’ 18°00’ 66°00’ C A R I B B E A N S E A A T L A N T I C O C E A N P U E R T O R I C O 8 Water Quality and Landscape Processes of Four Watersheds in Eastern Puerto Rico A A' B'B C C' El Cacique El Toro Pico del Oeste Pico del Este El Yunque 65°35’65°40’65°45’65°50’65°55’66°00’66°05’ 18°20’ 18°15’ 18°10’ 18°05’ Base from U.S. Geological Survey National Elevation Dataset 30-meter digital elevation model N:\Jeff\den11_hwcg00_0224_pp_murphy\chapter_A_figures\BW_figures\BWfigure_03A.ai Cayaguás Mameyes Icacos Canóvanas 0 0 4 8 KILOMETERS 8 MILES4 A'A EXPLANATION Elevation—In meters –29–0 1–200 201–400 401–600 601–800 801–1,000 Study watershed boundary Line of section Peaks higher than 1,000 meters A B C A R I B B E A N S E A P U E R T O R I C O WEST B EAST B' METERS WEST 1,000 800 600 400 400 0 METERS 1,000 800 600 400 400 0 A METERS SOUTHWEST 1,000 800 600 400 400 0 C NORTHEAST C' EAST A' VERTICAL EXAGGERATION X10 0 0 2 4 KILOMETERS 4 MILES2 C an óv an as Cayaguás Cayaguás Ic ac os Ic ac os M am ey es Mameyes Figure 3. A, Elevation and relief of eastern Puerto Rico, showing study watersheds (outlined), and three lines of cross section; B, elevation profiles along lines of cross section. Physiography, Geology, and Land Cover of Four Watersheds in Eastern Puerto Rico 9 N:\Jeff\den11_hwcg00_0224_pp_murphy\chapter_A_figures\BW_figures\BWfigure_04A.ai 0 45° 90°270° 135° 315° 225° 180° 0 45° 90°270° 135° 315° 225° 180° 0 45° 90°270° 135° 315° 225° 180° 0 45° 90°270° 135° 315° 225° 180° North North North NorthA. Canóvanas B. Cayaguás C. Icacos D. Mameyes 12 12 8 8 4 4 0 Pe rc en ta ge o f w at er sh ed a re a 12 12 8 8 4 4 0 Pe rc en ta ge o f w at er sh ed a re a Figure 4. Percentage of watershed area with each aspect, Canóvanas, Cayaguás, Icacos, and Mameyes watersheds (determined from elevation data from U.S. Geological Survey). Table 1. Geomorphic and geographic characteristics of study watersheds. [Values calculated from geographic information system data unless otherwise noted. Slope values are dimensionless] Characteristic Icacos Mameyes Canóvanas Cayaguás Area, square kilometers1 3.26 17.8 25.5 26.4 Minimum elevation, meters 620 83 70 156 Maximum elevation, meters 832 1,050 956 445 Mean elevation, meters 686 508 464 287 Mean hillslope of watershed2 0.222 0.365 0.255 0.189 Mean channel slope2 0.073 0.21 0.151 0.12 Main channel length, kilometers2 2.01 13.6 21.33 23.5 Total channel length, kilometers2 2.91 24.02 34.37 49.46 Dominant bedrock type Granitic Volcaniclastic Volcaniclastic Granitic Dominant land use history Forest Forest Agricultural Agricultural 1Díaz and others (2004). 2Larsen (1997). 10 Water Quality and Landscape Processes of Four Watersheds in Eastern Puerto Rico Geology Puerto Rico is a volcanic island-arc terrain with a geologic record spanning about 150 million years (Bawiec, 1998). It was formed by volcanism and sedimentation characteristic of tectonically active plate boundaries, and it consists of a core of volcanic and plutonic rocks surrounded by younger sedi- mentary rocks. The WEBB watersheds overlie igneous rocks typical of island-arc terrains. The Luquillo Mountains and the Loíza watershed are underlain largely by marine-deposited volcaniclastic rock with intrusions of granitic rock (figs. 5 and 6; Broedel, 1961; Seiders, 1971a,b,c; M’Gonigle, 1978, 1979; Rogers and others, 1979; Bawiec, 1998). Folding and large-scale strike-slip faulting have deformed strata, resulting in westerly dips such that strata grow progressively younger toward the northwest (Jolly and others, 1998). Geology in Watersheds Dominated by Volcaniclastic Rock The Mameyes and Canóvanas watersheds are underlain primarily by Cretaceous, marine-deposited, quartz-poor volcani- clastic rocks. The bedrock of the Mameyes watershed is primar- ily Fajardo, Tabonuco, and Hato Puerco Formations (41, 38, and 3 percent of the watershed, respectively; another 18 percent of the watershed is underlain by the Río Blanco quartz diorite intrusion, discussed in the next section). The Canóvanas water- shed is almost entirely underlain by the Hato Puerco and Lomas Formations (87 and 12 percent, respectively; fig. 6, table 2). The oldest of these formations, the Fajardo Formation, is about 105 to 100 million years old (Ma) (Jolly and others, 1998). Its lithology varies substantially and it has been divided into sev- eral members (Briggs, 1973; Briggs and Aguilar-Cortés, 1980). 0 0 10 20 KILOMETERS 10 20 MILES 66°30’67°00’ 18°30’ 18°00’ 66°00’ Canóvanas Cayaguás Icacos Mameyes EXPLANATION Geology Quaternary alluvium Tertiary through Cretaceous limey sedimentary rock Tertiary through Cretaceous nonlimey sedimentary rock Tertiary through Cretaceous granitoid rock Tertiary through Cretaceous gabbroic rock Cretaceous Subaerial volcaniclastic rock Tertiary marine volcaniclastic rock Tertiary and Cretaceous marine volcaniclastic rock Cretaceous marine volcaniclastic rock Alteration terrane age unknownSubmarine basalt and chert terrane Ultramafic and amphibolite terrane Water Geology modified from U.S. Geological Survey Open-File Report 98–38, Bawiec, 1998 Study watershed boundary P U E R T O R I C O Area shown in 6B Area shown in 6A N:\Jeff\den11_hwcg00_0224_pp_murphy\chapter_A_figures\BW_figures\BWfigure_05A.ai C A R I B B E A N S E A A T L A N T I C O C E A N Base from U.S. Geological Survey digital data, 1998 1:100,000, polyconic projection Central Meridian 67°W, Latitude of projection origin 18°N North American Datum of 1927 (NAD 1927) Figure 5. Generalized geology of Puerto Rico (from Bawiec, 1998); study watersheds outlined. Physiography, Geology, and Land Cover of Four Watersheds in Eastern Puerto Rico 11 0 0 1 2 KILOMETERS 1 2 MILES 0 0.5 0 0.5 1 KILOMETER 1 MILE A B Base from U.S. Geological Survey digital data, 1998 1:100,000, polyconic projection Central Meridian 67°W, Latitude of projection origin 18°N North American Datum of 1927 (NAD 1927) Base from U.S. Geological Survey digital data, 1998 1:100,000, polyconic projection Central Meridian 67°W, Latitude of projection origin 18°N North American Datum of 1927 (NAD 1927) Kh Kh Kh TKh TKh TKh TKh TKh Kfa Kfa Kta Kta Kta Klo Klo TKap TKha Canóvanas Icacos Mameyes Qa Kf Kf TKap TKapQa Qa Qa Qa Qa Qa Qa Ksl Ksl Ksl Kpgg Cayaguás Qa N:\Jeff\den11_hwcg00_0224_pp_murphy\chapter_A_figures\BW_figures\BWfigure_06A.ai Kpgqm 18°06’ 18°09’ 18°07’ 18°08’ 65°58’ 65°57’ 65°56’ 65°55’ 18°19’ 18°18’ 18°17’ 18°16’ 18°20’ 65°53’ 65°52’ 65°51’ 65°50’ 65°49’ 65°48’ 65°47’ 65°46’ 65°45’ Study watershed boundary Note: Geology modified from U.S. Geological Survey Open-File Report 98–38, Bawiec, 1998 EXPLANATION Alluvium Figuera Lava Lomas Formation Hato Puerco Formation Tabonuco Formation Fajardo Formation Hydrothermally altered rock Río Blanco quartz diorite Augite andesite porphyry Punta Guayanes quartz monzonite Punta Guayanes granodiorite San Lorenzo granodiorite-quartz diorite Fault Qa Kf Klo TKha Kpgqm TKh TKap Kh Kta Kpgg Ksl Kfa Geology Qa Figure 6. Detailed geology of study watersheds, eastern Puerto Rico (from Bawiec, 1998). 12 Water Quality and Landscape Processes of Four Watersheds in Eastern Puerto Rico Thick-bedded tuff breccias, with large clasts of andesite, domi- nate. They are interbedded with coarse-grained tuff, tuffaceous sandstone, and black cherty siltstone. The Fajardo Formation is largely noncalcareous, except for sporadic calcareous sandstone and siltstone beds and calcareous tuffs (Briggs, 1973). The Fajardo Formation grades into the overlying Tabo- nuco Formation, which is about 100 Ma (Jolly and others, 1998). It contains more mudstone and less very thick bedded sandstone and breccia than the Fajardo Formation (Seiders, 1971b). About 60 percent of the Tabonuco Formation consists of volcanic sandstones, which are usually light gray, calcare- ous, moderately sorted, and fine to very coarse grained, and which are composed primarily of angular volcanic rock frag- ments and subordinate plagioclase, clinopyroxene, and rare quartz grains (Seiders, 1971b). About 30 percent of the Tabo- nuco Formation consists of dark gray, thin-bedded to rarely thick-bedded, laminated mudstones, which can be slightly to highly calcareous. About 10 percent of the Tabonuco Forma- tion consists of volcanic breccia and conglomerate. The Hato Puerco Formation overlies the Tabonuco Formation and is about 97 to 90 Ma (fig. 6; Jolly and others, 1998). More than half of the formation is composed of andes- itic to basaltic volcanic sandstone, which contains volcanic lithic clasts, grains of plagioclase, clinopyroxene, and horn- blende, and sparse grains and boulders of limestone (Seiders, 1971a). About 30 to 40 percent of the formation consists of volcanic breccias that contain angular to subrounded pieces of andesite or basalt, with fragments of limestone that can make up as much as 30 percent of the breccia. Calcareous mudstone and andesitic to basaltic porphyritic lava are present in minor amounts. The Hato Puerco contains more carbonate material than the Fajardo or Tabonuco Formations, particularly in the vicinity of the Río Canóvanas stream channel, where outcrops of calcareous shales and tuffaceous limestone are common (Meyerhoff, 1931; Seiders, 1971a). Calcareous rocks require relatively small outcrop areas to have a marked effect on stream chemistry (Stallard, 1995), and the larger amount of carbonate rocks in the Canóvanas watershed compared with the Mameyes watershed is reflected in stream chemistry (Stallard and Murphy, 2012). The Lomas Formation, found in the southern part of the Canóvanas watershed, consists of gray-green, thick-bedded, poorly sorted andesitic to basaltic volcanic breccia and sandstone, with minor pillow lava locally (Seiders, 1971c). Strata commonly contain pumice, vesicular lava fragments, and grains of plagioclase and clinopyroxene, and calcite is locally abundant in the matrix (Seiders, 1971c). This formation is seen at the surface only in fault blocks and is not observed in depositional contact with other formations; on the basis of its similarities with other formations in the region, it has been estimated to be between 90 and 85 Ma (Jolly and others, 1998). The intrusion of the Río Blanco quartz diorite (dis- cussed in the next section) metamorphosed the Fajardo, Tabonuco, and Hato Puerco Formations within a 2–4 km contact zone (Bawiec, 1998), resulting in rocks that are typi- cally harder; the divide between the Mameyes and Icacos watersheds was formed on the slightly hardened rocks of the metamorphic aureole. In the metamorphosed rocks, hornblende, garnet, epidote, chlorite, iron sulfide minerals, gold, copper, and quartz, and aplitic veins have been identi- fied (Seiders, 1971a,c; Bawiec, 1998). Gold-, silver-, and copper-bearing deposits were first worked by the Spanish in the early 1600s; production was limited to several small adits within the contact zone (Wadsworth, 1970; Cardona, 1984; Bawiec, 1998). Placer gold from weathered metamorphosed rocks was mined intermittently in the Mameyes and Loíza watersheds beginning in the 1500s. Geology in Watersheds Dominated by Granitic Rock The Icacos and Cayaguás watersheds are primarily underlain by Upper Cretaceous and Tertiary granitic rocks. Almost all (99 percent) of the Icacos watershed is underlain by the Río Blanco quartz diorite (fig. 6, table 2), which intruded the surrounding volcaniclastic rocks about 49 to 42 Ma (Smith and others, 1998). The Río Blanco quartz diorite is light-gray, medium- to coarse-grained, and commonly contains plagio- clase, quartz, amphibole, and minor biotite in a fine- to medium- crystalline groundmass of quartz, amphibole, orthoclase, and plagioclase, and with accessory zircon, tourmaline, sphene, and epidote (Seiders, 1971a,b; table 3). Plagioclase ranges from fresh to sericitized. Hornblende is commonly altered to chlorite, calcite, and epidote (Seiders, 1971a). Biotite is commonly fresh, but it is also altered to chlorite (Murphy and others, 1998). In some areas the Río Blanco quartz diorite contains pyrite, molyb- denite, or disseminated chalcopyrite (Bawiec, 1998). Table 2. Geology of study watersheds, eastern Puerto Rico.1 Watershed and formation Symbol Rock type Percent Canóvanas Hato Puerco Formation Kh Volcaniclastic 87.2 Lomas Formation Klo Volcaniclastic 11.5 Río Blanco quartz diorite TKh Intrusive 0.5 Alluvium Qa Alluvium 0.4 Augite andesite porphyry Tkap Intrusive 0.3 Cayaguás San Lorenzo granodiorite- quartz diorite Ksl Intrusive 88.0 Punta Guayanes quartz monzonite Kpgqm Intrusive 11.3 Alluvium Qa Alluvium 0.6 Punta Guayanes granodiorite Kpgg Intrusive 0.02 Icacos Río Blanco quartz diorite TKh Intrusive 98.5 Fajardo Formation Kfa Volcaniclastic 1.2 Tabonuco Formation Kta Volcaniclastic 0.4 Mameyes Fajardo Formation Kfa Volcaniclastic 40.9 Tabonuco Formation Kta Volcaniclastic 37.9 Río Blanco quartz diorite TKh Intrusive 17.8 Hato Puerco Formation Kh Volcaniclastic 3.3 Alluvium Qa Alluvium 0.04 1Adapted from Bawiec (1998). Physiography, Geology, and Land Cover of Four Watersheds in Eastern Puerto Rico 13 About 88 percent of the Cayaguás watershed is under- lain by the San Lorenzo granodiorite–quartz diorite intrusion (fig. 6, table 2), which has been dated to be about 73 Ma (Rogers and others, 1979; Bawiec, 1998). This intrusion occupies a 500-km² area in southeastern Puerto Rico. It is a medium- to dark-gray, granitic-textured, medium-crystalline, generally massive granodiorite which grades into quartz dio- rite and locally into quartz monzonite. It is composed gener- ally of plagioclase, potassium feldspar, quartz, hornblende, and minor biotite and accessory magnetite, sphene, apatite, and zircon (Rogers and others, 1979) (table 3). Many plagio- clase grains are partly sericitized and argillized. Hornblende is commonly altered to chlorite and epidote, and biotite is commonly altered to chlorite. About 11 percent of the Cayaguás watershed has been mapped as Punta Guayanes quartz monzonite, which has been dated at about 66 Ma (Rogers and others, 1979). This rock type ranges in composition from quartz monzonite to granodiorite, but in general it is granitic-textured, medium crystalline, and composed of abundant plagioclase, potas- sium feldspar, and quartz, with subordinate hornblende, minor biotite, and accessory magnetite, sphene, apatite, and zircon (table 3). Most plagioclase grains are sericitized to various degrees, and hornblende is commonly altered to chlorite and epidote (Rogers and others, 1979). The southern part of the watershed was mapped to greater detail by Rogers and others (1979) than the northern part (Broedel, 1961); Broedel clas- sified all rocks simply as granodiorite. The greater amount of quartz and lesser amount of plagioclase in the quartz monzo- nite may make the monzonite more resistant to weathering than the granodiorite, and they may explain why hillslopes underlain by monzonite in the watershed have higher eleva- tion and greater relief that elsewhere in the basin (Larsen, 1997). The quartz monzonite may actually underlie a larger portion of the Cayaguás basin than the reported 11 percent, because this accentuated topography is expressed elsewhere in areas mapped as granodiorite. The effect of geology on channel morphology in the WEBB watersheds is evident in the bed material and par- ticle size (see cover photograph). Stream channels in the Icacos and Cayaguás watersheds typically have sand-bed channels and floodplains that are dominantly sandy allu- vium with occasional boulders. Median grain size (D50) is 0.6 millimeters (mm) and 0.5 mm, respectively, in the two channels (Larsen, 1997). Channels in the Mameyes and Canóvanas watersheds are generally characterized by little to no floodplain and channels lined with bedrock, boulders, or cobbles, with little bed sediment. Median grain size in these two channels is 70 mm and 110 mm, respectively (Larsen, 1997). Soil and Saprolite Intense, continuous weathering in the warm, humid- tropical climate of the Luquillo Mountains and the Loíza watershed has produced a mantle of soil overlying both gra- nitic and volcaniclastic bedrock that is as much as 24 m thick (Simon and others, 1990; White and others, 1998). Below the surface-soil horizons, which show strong bioturbation (Brown and others, 1983; Brown and others, 1995; González and others, 1996) and subsoil movement (Lewis, 1974; Larsen and others, 1999), is commonly a layer of saprolite. Saprolite is soil that is derived from bedrock that has been thoroughly decomposed through isovolumetric chemical weathering; it remains in place and retains the appearance of the bedrock structure (American Geological Institute, 1976). The com- position of the underlying bedrock, particularly the quartz content, determines the chemical and physical properties of the saprolite (Meyerhoff, 1931; Schellekens and others, 2004; Buss and White, 2012; Stallard, 2012). Soils in Watersheds Dominated by Volcaniclastic Rock The volcaniclastic rocks underlying the Mameyes and Canóvanas watersheds are generally quartz poor, and they are composed predominantly of basic feldspars and ferro- magnesian minerals (Seiders, 1971a,b,c; Briggs, 1973). In the saprolite, these minerals weather to kaolinite and other clay minerals of the 1:1 lattice type, iron and aluminum oxides, and small amounts of quartz (Huffaker, 2002). Relic grains of plagioclase and mafic minerals are recognizable in some hand specimens of the saprolite (Scatena, 1989). Geo- electric soundings indicate a gradual change from unweath- ered bedrock to less-weathered saprolite to highly weathered saprolite; active weathering occurs in the less-weathered saprolite (Schellekens and others, 2004). Clay-lined faults and joints, which are common in the saprolite, make primary sedimentary structures difficult to distinguish. Table 3. Mineral content of granitic intrusions, eastern Puerto Rico. [--, not available] Mineral Average mineral content of bedrock (weight percent) Río Blanco1 San Lorenzo2 Quartz monzonite of Punta Guayanes2 Plagioclase 55.8 51 35 Quartz 26.0 20.5 33 Amphibole 7.3 11 1.5 K-feldspar 5.0 11 28 Biotite 3.6 4.5 2 Accessory3 1.4 -- -- Apatite 0.2 -- -- Iron oxide 0.7 1.5 0.5 Total 100 99.5 100 1Seiders (1971b). 2Rogers and others (1979). 3Accessory minerals: Río Blanco—epidote, sphene, tourmaline, zircon; San Lorenzo—apatite, sphene, zircon, and minor augite, myrmekite; Punta Guayanes—apatite, sphene, zircon. The volcaniclastic saprolite is overlain by clay-rich, quartz-poor soils that are typically 0.8 to 1.0 m deep (Huf- faker, 2002). Soils developed on the volcaniclastic rocks are finer grained than those developed on the granitic rocks, and they have lower permeability and infiltration rates than granitic soils (Simon and others, 1990; Murphy and Stallard, 2012). Dense clays in the volcaniclastic soils retard deep infiltration, causing water to follow a shallow trajectory and reach streams quickly (McDowell and others, 1992). Because of the generally high soil moisture content in the Luquillo Mountains, soils are quickly saturated during storms, leading to saturation overflow (Larsen and others, 1999). As such, during periods of intense rainfall, water is commonly observed flowing over the surface in rills and ephemeral gullies. At the highest elevations in the Mameyes watershed, cor- responding with palo colorado forest (see Vegetation and Land Use section), soils overlying the volcaniclastic rock consist primarily of the Yunque clay loam and Moteado clay soils, which are Oxisols, and the Guayabota clay loam and Palm clay soils, which are Inceptisols (Huffaker, 2002). These soils range from poorly drained (on concave slopes) to moderately well drained (on convex slopes). The Inceptisols presumably reflect the greater rainfall in the highest elevations and wetter areas near rivers. At lower elevations, corresponding with tabonuco forest, soils overlying the volcaniclastic rocks consist primarily of Zarzal clay and Cristal clay loam soils, which are Oxisols, and range from somewhat poorly drained (concave slopes) to well drained (convex slopes) (Huffaker, 2002). The Canóvanas watershed contains a variety of soil types, including Oxisols (Yunque, Moteado, Zarzal, and Cristal clay or clay loams, discussed in the previous paragraph, and the deep, moderately well drained Los Guineos clay), Incepti- sols (Mucara silty clay loam and Caguabo clay loam, which are shallow to moderately deep and well drained), Ultisols (Humatas clay and Naranjito silty clay loam, which are deep to moderately deep and well drained), and minor Mollisols (Reilly silty clay loam, found in the floodplain) (Boccheciamp, 1977; Huffaker, 2002). As a result of past agricultural cultiva- tion, much of the original surface layers of several of the soils in the lower part of the Canóvanas watershed, including the Los Guineos, Humatas, Mucara, and Naranjito soils, have been removed by erosion (Boccheciamp, 1977). Soils in Watersheds Dominated by Granitic Rock The quartz-rich granitic rocks underlying the Icacos and Cayaguás watersheds weather to a saprolite composed of quartz, weathered biotite, and kaolinite, and other clays and sesquioxides (Boccheciamp, 1978; Murphy and others, 1998; White and others, 1998). Huffaker (2002) contends that the granitic rocks weather more rapidly than the volcanic rocks, and that if other conditions are equal, the regolith overlying the granitic rocks is thicker than that in volcanic areas. The Río Blanco quartz diorite, which underlies the Icacos watershed, has been the subject of extensive studies on chemical weather- ing rates (including White and Blum, 1995; Dong and others, 1998; Murphy and others, 1998; White and others, 1998; Schulz and White, 1999; Turner and others, 2003; Buss and others, 2005, 2008, 2009; Fletcher and others, 2006; Pett-Ridge and others, 2009; Minyard and others, 2011; Buss and White, 2012). The relatively young age of the Río Blanco quartz diorite, combined with high rainfall and temperature, results in the Icacos watershed having one of the highest documented chemical weathering rates of granitic rocks in the world (White and Blum, 1995; White and others, 1998). The weathering front has been estimated to proceed about 1 centimeter (cm) per 100 years on the basis of cosmogenic evidence (Brown and others, 1995). Stallard (1995) argues that minerals such as intergranular calcite and volatiles in fluid inclusions derived from the late-stage cooling of magma greatly enhance the sus- ceptibility of younger plutons, such as the Río Blanco quartz diorite, to weathering, compared with identical rock types of Precambrian age, where volatiles have migrated out and late-stage minerals have recrystallized. The transition between the quartz-diorite bedrock and saprolite consists of a 20- to 60-cm-thick zone characterized by fracture-bound concentric shells, termed rindlets (Turner and others, 2003; Fletcher and others, 2006; Buss and others, 2008; Buss and White, 2012). The majority of weathering in the watershed occurs in this tran- sition zone, as individual rindlets become progressively more altered away from joint planes. The San Lorenzo granodiorite– quartz diorite, which underlies the Cayaguás watershed, has not been well studied, but because of its similar composition and age it presumably weathers in a comparable manner. On hillslopes, saprolite is overlain by 0.5–2 m of soil, which is typically sandy and loamy, with a thin, organic-rich A horizon and a layer of translocated clays concentrated in the B horizon (Simon and others, 1990; Huffaker, 2002). Soil thick- ness on ridges in the Icacos watershed has been observed to be more than ten-fold greater (Simon and others, 1990). Soils on the granitic rocks have higher permeability and infiltration rates than those on the volcaniclastic rocks, allowing greater water infiltration (Simon and others, 1990). The coarse soil is easily eroded and is highly susceptible to landslides when subjected to intensive agricultural activity or other disturbance (Larsen, 2012). Soils overlying the Río Blanco quartz diorite in the Ica- cos watershed primarily consist of Picacho and Utuado loamy soils. These Inceptisols are very deep (greater than 150 cm) and somewhat poorly drained (Huffaker, 2002). The Icacos flood plain and low terraces at middle elevations are overlain by the Icacos loam, an Inceptisol which is very deep (greater than 150 cm) and somewhat poorly drained. Quartz persists through the entire profile, becoming more fine grained nearer the surface (Brown and others, 1995). Soils at the highest elevations in the Cayaguás water- shed consist primarily of the Lirios silty clay loam, which is a well-drained Ultisol (Boccheciamp, 1978). Soils at lower elevations primarily consist of the Pandura sandy loam, an Inceptisol. Owing to past agricultural cultivation, both of these soils have lost much of their original surface lay- ers through erosion. By using geomorphic evidence, Larsen 14 Water Quality and Landscape Processes of Four Watersheds in Eastern Puerto Rico Physiography, Geology, and Land Cover of Four Watersheds in Eastern Puerto Rico 15 and Santiago-Román (2001) estimate that the surface of the Cayaguás basin has been lowered an average of 660 mm since about 1820, comparable to the 500 mm of surface-lowering estimated by using cosmogenic Be10 (Brown and others, 1998). Vegetation and Land Use All of the natural ecological life zones of Puerto Rico (as developed by Ewel and Whitmore, 1973) are forest (fig. 7). Three life zones dominate the island: subtropical moist for- est (58 percent of the main island), subtropical wet forest (23 percent), and subtropical dry forest (18 percent) (Birdsey and Weaver, 1982). Subtropical rain forest, lower montane wet forest, and lower montane rain forest occupy a little more than 1 percent of the island at high elevations, including most of the Mameyes and Icacos watersheds. Historical accounts suggest that, like many other Caribbean islands, Puerto Rico was heavily forested (about 95 percent) prior to European settlement (Wadsworth, 1950; Birdsey and Weaver, 1982). The island’s diverse physiography and tropical location result in diverse flora and fauna and a high level of endemism (Gould and others, 2008); however, plant diversity is lower than that at mainland tropical sites (Ashton and others, 2004). In the 16th century, European settlers began clearing forests for pasture, cropland, timber, and fuel. By 1828, forest covered about 66 percent of Puerto Rico (Wadsworth, 1950). Pressure on land resources increased in the 19th and 20th centuries owing to population growth, expanding production of sugar cane and coffee, economic depression, and the effects of several severe hurricanes. By the late 1940s, forest cover of Puerto Rico had declined to 6 percent; cropland and pasture each accounted for about 42 percent of land cover; the remain- ing 10 percent contained urban areas and wasteland (Koenig, 1953; Birdsey and Weaver, 1982). During this period, Puerto Rico was one of the most severe1y deforested regions in the world. Soil erosion related to this deforestation released high loads of sediment from hillslopes to footslopes, valley floors, N:\Jeff\den11_hwcg00_0224_pp_murphy\chapter_A_figures\BW_figures\BWfigure_07A.ai Suptropical moist forest Subtropical rain forest Lower montane rain forest EXPLANATION Lower montane wet forest Subtropical dry forest Subtropical wet forest Holdridge Life Zone Map data from Ewel and Whitmore, 1973, as digitized by the Luquillo Long-Term Ecological Research program Canóvanas Cayaguás Icacos Mameyes 18°30’ 18°00’ 67°00’ 66°30’ 66°00’ 0 0 10 20 KILOMETERS 10 20 MILES C A R I B B E A N S E A A T L A N T I C O C E A N Base from USDA Forest Service digital data, 1973, 1:240,000 Puerto Rico State Plane Coordinate System of 1927 Standard parallels 18°26’N and 18°02’N, central meridian 66°26’W Puerto Rico Datum P U E R T O R I C O Figure 7. Ecological life zones of Puerto Rico (from Luquillo Long-Term Ecological Research program; digitized from Ewel and Whitmore, 1973). 16 Water Quality and Landscape Processes of Four Watersheds in Eastern Puerto Rico and streams channels (Larsen, 2012). Much of the remaining forest was in the Luquillo Mountains, because access to this area was limited by steep slopes, high annual rainfall, and its designation as a reserve by the Spanish crown in 1876 and later by the U.S. Government (Zimmerman and others, 1995). A shift in Puerto Rico’s economic base that began in the late 1940s, partly as a result of the U.S. Government’s “Opera- tion Bootstrap,” led to reforestation (Rudel and others, 2000). In 1934, about 43 percent of the island’s gross national prod- uct came from agriculture but only 7 percent from industry; by 1996, agriculture and industry contributed 1 percent and 41 percent respectively to the gross national product (López and others, 2001). As Puerto Ricans moved from rural areas to cities on the island or migrated to the United States, agricultural land cover in Puerto Rico declined by 95 percent from 1951 to 2000 (Kennaway and Helmer, 2007). Much of this abandoned land reverted to forest. By 1980, forest cover of Puerto Rico had increased to 31 percent (Franco and others, 1997). In 2003, forest, woodland, and shrubland represented about 44 percent of land area (Gould and others, 2012) (figs. 8 and 9). Meanwhile, population density of Puerto Rico has increased more than threefold during the last century, resulting in one of the highest densities in the world: more than 400 people/km2 (López and others, 2001). Puerto Rico now has the highest road density of any Caribbean island, 2.5 km of paved roads per square kilometer (Lugo, 1996). On the basis of population density and remote sensing of developed areas, about 16 percent of the island is considered urban land use, 50 percent is suburban, and 34 percent is rural (Martinuzzi and others, 2007). Eleven percent of the island has developed surface, and about 60 percent of this surface corresponds with high-density development (Mar- tinuzzi and others, 2007). In the San Juan metropolitan area, filling of the mangroves and mudflats has restricted circula- tion and impaired water quality (Webb and Gómez-Gómez, 1998). Land has been converted to urban uses mostly in lowland regions and on lesser slopes. Land Cover of Forested Watersheds The four study watersheds were selected, in part, for differing land cover. The Icacos watershed is almost entirely covered by mature montane wet evergreen forest, which N:\Jeff\den11_hwcg00_0224_pp_murphy\chapter_A_figures\BW_figures\BWfigure_08A.ai EXPLANATION Land Cover Barren Forest Nonforest vegetation Urban Woodland and shrubland Water 0 0 10 20 KILOMETERS 10 20 MILES Study watershed boundary 67°00’ 18°30’ 18°00’ 66°30’ 66°00’ C A R I B B E A N S E A A T L A N T I C O C E A N Cayaguás Canóvanas Icacos Mameyes m 01 KILOMETERS Land cover data from USDA Forest Service, International Institute of Tropical Forestry, Puerto Rico GAP Analysis Project (2003) Base from USDA Puerto Rico GAP Analysis Project digital data, 2006 Lambert Conformal Conic projection Standard parallels 18°26’N and 18°02’N, central meridian 66°26’W North American Datum of 1983, High Accuracy Reference Network (NAD 1983 HARN) P U E R T O R I C O Area shown in figure 9 Figure 8. Generalized land cover of Puerto Rico in 2003 (from Gould and others, 2008). Physiography, Geology, and Land Cover of Four Watersheds in Eastern Puerto Rico 17 N:\Jeff\den11_hwcg00_0224_pp_murphy\chapter_A_figures\BW_figures\BWfigure_09A.ai Map data modified from U.S. Department of Agriculture Forest Service, International Institute of Tropical Forestry, Puerto Rico GAP Analysis Project (2003) Cayaguás 18°20’ 18°15’ 18°10’ 18°05’ 66°00’ 66°55’ 66°50’ 66°45’ Canóvanas Icacos Mameyes EXPLANATION Land Cover Elfin woodland forest Sierra Palm forest Palo colorado forest Tabonuco forest Other wet forest Moist forest Dry forest Grassland and pasture Agriculture Urban Barren Water Study watershed boundary 0 0 6 KILOMETERS3 3 6 MILES C A R I B B E A N S E A Base from USDA Puerto Rico GAP Analysis Project digital data, 2006 Lambert Conformal Conic projection Standard parallels 18°26’N and 18°02’N, central meridian 66°26’W North American Datum of 1983, High Accuracy Reference Network (NAD 1983 HARN) Figure 9. Detailed land cover of study watersheds in 2003 (from Gould and others, 2008). 18 Water Quality and Landscape Processes of Four Watersheds in Eastern Puerto Rico includes about 85 percent palo colorado forest, 13 percent sierra palm forest, and 2 percent other montane wet forest (Gould and others, 2008) (fig. 9; table 4). Prominent trees are palo colorado (Cyrilla racemiflora), caimitillo (Micropho- lis garciniaefolia, M. chrysophylloides), and laurel sabino (Magnolia splendens) (Foster and others, 1999). The Mam- eyes watershed also consists primarily of mature montane wet evergreen forest, but it is larger and has a much greater variation in elevation than the Icacos watershed, and thus it has more varied forest. Land cover is about 57 percent tabonuco forest, 25 percent sierra palm forest, 14 percent palo colorado forest, 1 percent elfin woodland, and 2 percent other montane wet forest (Gould and others, 2008, 2012) (table 4). Above 600 m, vegetation in the Mameyes watershed is pri- marily palo colorado and sierra palm forest. Elfin woodland, characterized by short, slow-growing vegetation, develops on ridges above 900 m elevation (Weaver, 1986, 1989, 1990). Elevations below 600 m are primarily tabonuco forest. This forest is dominated by the tabonuco tree (Dacryodes excelsa), but 170 tree species have been described for this forest type (Foster and others, 1999). Forest cover in much of the Icacos and Mameyes water- sheds persisted during the nadir of forest cover in Puerto Rico in the 1940s, and those watersheds have maintained a con- tinuous forest canopy since that time (table 4). As discussed above, steep topography, harsh conditions, and governmental protection limited access for development or forest conver- sion. Although the Luquillo Mountains were less affected by development than other areas of Puerto Rico, they are not pristine; cattle production, shade coffee, timber extraction, subsistence farming, and charcoal production were practiced in the late 1800s and in the 1900s, primarily below 600 m elevation. The USDA Forest Service purchased much of the lower-elevation land in the Luquillo Mountains in the 1930s (Zimmerman and others, 1995). Nonnative bamboo was introduced from Asia in the 1930s and 1940s along roads to stabilize slopes; however, during major disturbances such as hurricanes, extensive bamboo root systems can actually help cause catastrophic slope failures (O’Connor and others, 2000). Bamboo is expanding downstream in streams within the Mameyes watershed at a rate of 8 meters per year (O’Connor and others, 2000). Large stands of bamboo are observed in the Mameyes and Icacos watersheds, especially near roads. Land Cover of Developed Watersheds In the Canóvanas watershed, land cover is currently about 71 percent forest, 26 percent grassland or pasture, and 3 percent urban (Gould and others, 2008) (table 4). Above 600 m elevation, land cover is primarily palo colorado for- est, similar to the Icacos and Mameyes watersheds (fig. 9, table 4). Between 400 and 600 m, land cover is primarily tabonuco forest and grassland or pasture. Below 400 m, land cover is primarily secondary lowland moist forest and grassland or pasture. The Canóvanas watershed has a mixed land-use history. Elevations above 600 m were forested in the middle 20th century, but most of the lower watershed was used for pasture, coffee, and fruit crops, and to a lesser degree, bananas and sugar cane (Kennaway and Helmer, 2007). Between the 1950s and the 1970s, the predominant land cover shifted from agriculture and pasture to forest, and forest cover in the watershed increased by 44 percent (table 4). These decades correspond with a low point in population density of the Cubuy barrio, which is entirely within the Canóvanas watershed and represents the majority of the land area in the watershed (fig. 10). Current population density has returned to that of the 1950s, but land use is pre- dominantly residential as opposed to agricultural (Martinuzzi and others, 2007; Parés-Ramos and others, 2008). Current land cover in the Cayaguás watershed is 57 percent pasture, 41 percent forest or shrubland (primarily lowland moist forest and montane wet forest) and 2 percent urban (fig. 9, table 4). Fifty years ago, the Cayaguás watershed was almost entirely devoted to agricultural activities, such as pasture, fruit crops, and tobacco (Kennaway and Helmer, 2007). Forest cover has increased about 20 percent since the late 1970s, primarily at the expense of pasture and grasslands (table 4). Population densities in the Quebrada Arenas barrio, which is entirely within the Cayaguás watershed and repre- sents most of the land area in the watershed, reached a low point in the 1980 census (fig. 10). Population has since been rising owing to increased residential land use (Martinuzzi and others, 2007; Parés-Ramos and others, 2008). Summary Four watersheds in eastern Puerto Rico have been studied since 1991 as part of the U.S. Geological Survey’s Water, Energy, and Biogeochemical Budget program. These watersheds differ in geology, soils, and land cover. The Icacos watershed is underlain by granitic rock, which weathers to quartz- and clay-rich, sandy soils. Land cover is primarily old-growth wet forest, and human disturbance has been minimal. The Mameyes watershed is primarily under- lain by volcaniclastic rock, which weathers to quartz-poor, fine-grained soils. This watershed also contains largely old- growth wet forest. The Canóvanas watershed is underlain by volcaniclastic rocks and soils that are similar to those in the Mameyes, but historical agricultural practices at lower elevations have led to substantial erosion of soils. This watershed has largely reforested since the 1950s. The Cay- aguás watershed is underlain by granitic rock; this water- shed was also historically used for agriculture and has been deeply eroded. It has undergone some reforestation since the 1970s. Past soil erosion in the Canóvanas and Cayaguás watersheds released high loads of sediment from hillslopes to footslopes, valley floors, and streams channels. The refor- estation of these watersheds, and of Puerto Rico, may serve as a prototype for the hydrological and chemical response of other tropical areas that are shifting from an agricultural to an industrial economy. Physiography, Geology, and Land Cover of Four Watersheds in Eastern Puerto Rico 19 Ta bl e 4. La nd c ov er a nd la nd -c ov er d yn am ic s in s tu dy w at er sh ed s, e as te rn P ue rto R ic o. [L an d- co ve r d yn am ic s a re in di ca te d as p er ce nt ag e in cr ea se o r d ec re as e fr om 1 93 6, 1 95 0 (L ug o an d ot he rs , 2 00 4) , a nd 1 97 7 (R am os a nd L ug o, 1 99 4) fo r w oo dy (f or es t, sh ru bl an d, a nd w oo dl an d) , h er ba ce ou s (g ra ss la nd , p as tu re s, ag ric ul tu re ), an d de ve lo pe d ar ea s. In 1 93 6 an d 19 50 , t he C ay ag uá s w at er sh ed w as p rim ar ily a gr ic ul tu ra l a nd n on fo re st ed . - -, no t a na ly ze d] La nd c ov er ty pe Ic ac os M am ey es Ca nó va na s Ca ya gu ás 20 03 la nd c ov er 1 pe rc en t Pe rc en t ch an ge s in ce 20 03 la nd c ov er pe rc en t Pe rc en t ch an ge s in ce 20 03 la nd c ov er pe rc en t Pe rc en t ch an ge s in ce 20 03 la nd c ov er pe rc en t Pe rc en t ch an ge s in ce 19 36 19 50 19 77 19 36 19 50 19 77 19 36 19 50 19 77 19 36 19 50 19 77 W oo dy 99 .9 0 0 0 99 .6 1 0 0 71 .3 47 44 2 41 .7 -- -- 20 Pa lo c ol or ad o fo re st 85 .4 -- -- -- 13 .9 -- -- -- 14 .2 -- -- -- 0. 00 -- -- -- Si er ra p al m fo re st 12 .7 -- -- -- 25 .3 -- -- -- 1. 60 -- -- -- 0. 00 -- -- -- E lfi n w oo dl an d 0. 00 -- -- -- 1. 07 -- -- -- 0. 38 -- -- -- 0. 00 -- -- -- Ta bo nu co fo re st 0. 00 -- -- -- 57 .0 -- -- -- 24 .5 -- -- -- 0. 00 -- -- -- O th er m on ta ne w et f or es t, sh ru bl an d, w oo dl an d 1. 84 -- -- -- 2. 36 -- -- -- 10 .4 -- -- -- 21 .5 -- -- -- Lo w la nd m oi st fo re st , s hr ub la nd , w oo dl an d 0. 00 -- -- -- 0. 00 -- -- -- 20 .2 -- -- -- 19 .8 -- -- -- A ba nd on ed a nd a ct iv e co ffe e pl an ta tio ns 0. 00 -- -- -- 0. 00 -- -- -- 0. 00 -- -- -- 0. 37 -- -- -- H er ba ce ou s 0. 15 0 0 0 0. 32 –1 0 0 25 .6 8 –4 1 –4 8 3 56 .7 -- -- –2 3 M oi st g ra ss la nd s a nd p as tu re s 0. 15 -- -- -- 0. 32 -- -- -- 25 .6 8 -- -- -- 56 .7 -- -- -- D ev el op ed 0. 00 0 0 0 0. 02 0 0 0 2. 79 3 3 –4 1. 55 -- -- 2 Lo w -d en si ty u rb an d ev el op m en t 0. 00 -- -- -- 0. 02 -- -- -- 2. 79 -- -- -- 1. 55 -- -- -- O th er 0. 00 0 0 0 0. 05 0 0 0 0. 25 0 0 0 0. 08 -- -- 0 B ar re n 0. 00 -- -- -- 0. 04 -- -- -- 0. 25 -- -- -- 0. 00 -- -- -- Fr es hw at er 0. 00 -- -- -- 0. 00 4 -- -- -- 0. 00 -- -- -- 0. 08 -- -- -- To ta l 10 0 10 0 10 0 10 0 1 F ro m G ou ld a nd o th er s, 20 08 . 20 Water Quality and Landscape Processes of Four Watersheds in Eastern Puerto Rico Acknowledgments This chapter was greatly improved by reviews by John F. Walker and Richard M.T. Webb of the U.S. Geological Survey and F.N. Scatena of the University of Pennsylvania. References American Geological Institute, 1976, Dictionary of geological terms: Garden City, N.Y., Anchor Press. 472 p. Ashton, P.S., and CTFS Working Group, 2004, Floristics and vegetation of the forest dynamics plots, in Losos, E.C., and Leigh, E.G. Jr., eds., Tropical forest diversity and dynamism—Findings from a large-scale plot network: Chicago, Ill., The University of Chicago Press, p. 90–102. Baedecker, M.J., and Friedman, L.C., 2000, Water, energy, and biogeochemical budgets—A watershed research program: U.S. Geological Survey Fact Sheet 165–99, 4 p. Bawiec, W. J., 1998, Geology, geochemistry, geophysics, mineral occurrences, and mineral resource assessment for the commonwealth of Puerto Rico: U.S. Geological Survey Open-File Report 98–38, CD-ROM. Birdsey, R.A., and Weaver, P.L., 1982, The forest resources of Puerto Rico: U.S. Department of Agriculture Forest Service, Southern Forest Experiment Station, Resources Bulletin SO–85, 65 p. Boccheciamp, R.A., 1977, Soil survey of the Humacao area of eastern Puerto Rico: U.S. Department of Agriculture Soil Conservation Service, 103 p. Boccheciamp, R.A., 1978, Soil survey of San Juan area of Puerto Rico: U.S. Department of Agriculture Soil Conserva- tion Service, 141 p. Briggs, R.P., 1973, The Lower Cretaceous Figuera Lava and Fajardo Formation in the stratigraphy of northeastern Puerto Rico: U.S. Geological Survey Bulletin 1372–G, 10 p. Briggs, R.P., and Aguilar-Cortés, E., 1980, Geologic map of the Fajardo and Cayo Icacos quadrangles, Puerto Rico: U.S. Geological Survey Miscellaneous Geologic Investigations Map I–1153, scale 1:20,000, 1 sheet. Broedel, C.H., 1961, Preliminary geologic map showing iron and copper prospects in the Juncos quadrangle, Puerto Rico: U.S. Geological Survey Miscellaneous Geologic Investiga- tions Map I–326. Brown, E.T., Stallard, R.F., Larsen, M.C., Bour, D.L., Raisbeck, G.M., and Yiou, F., 1998, Pre-anthropogenic denudation rates of a perturbed watershed (Cayaguás River, Puerto Rico) esti- mated from in-situ-produced 10Be in river-borne quartz: Earth and Planetary Science Letters, v. 160, p. 723–728. Brown, E.T., Stallard, R.F., Larsen, M.C., Raisbeck, G.M., and Yiou, F., 1995, Denudation rates based on accumula- tion of in-situ-produced 10Be compared with watershed mass balance results in the Luquillo Experimental Forest, Puerto Rico: Earth and Planetary Science Letters, v. 129, p. 193–202. Brown, Sandra, Lugo, A.E., Silander, Susan, and Liegel, Leon, 1983, Research history and opportunities in Luquillo Experimental Forest: U.S. Department of Agriculture For- est Service, Southern Forest Experimental Station, New Orleans, LA, General Technical Report SO–44, 128 p. Buss, H.L., Bruns, M.A., Schultz, M.J., Moore, J., Mathur, C.F., and Brantley, S.L., 2005, The coupling of biologi- cal iron cycling and mineral weathering during saprolite formation, Luquillo Mountains, Puerto Rico: Geobiology, v. 3, p. 247–260. Buss, H.L., Mathur, R., White, A.F., and Brantley, S.L., 2009, Phosphorous and iron cycling in deep saprolite, Luquillo Mountains, Puerto Rico: Chemical Geology, v. 269, p. 52–61. Figure 10. Population density in the Cubuy barrio (in Canóvanas watershed) and the Quebrada Arenas barrio (in Cayaguás watershed) (from U.S. Census, 2000). 1950 1960 1970 1980 1990 Po pu la tio n de ns ity , i n pe op le p er s qu ar e ki lo m et er 2000 150 125 100 75 50 25 0 Year EXPLANATION Cubuy barrio Quebrada Arenas barrio N:\Jeff\den11_hwcg00_0224_pp_murphy\chapter_A_figures\BW_figures\BWfigure_10A.ai Physiography, Geology, and Land Cover of Four Watersheds in Eastern Puerto Rico 21 Buss, H.L., Sak, P.B., Webb, S.M., and Brantley, S.L., 2008, Weathering of the Río Blanco quartz diorite, Luquillo Mountains, Puerto Rico—Coupling oxidation, dissolution, and fracturing: Geochimica et Cosmochimica Acta, v. 72, p. 4488–4507. Buss, H.L., and White, A.F., 2012, Weathering processes in the Río Icacos and Río Mameyes watersheds in eastern Puerto Rico, ch. I in Murphy, S.F., and Stallard, R.F., eds., Water quality and landscape processes of four watersheds in east- ern Puerto Rico: U.S. Geological Survey Professional Paper 1789, p. 249–262. Calvesbert, R.J., 1970, Climate of Puerto Rico and the U.S. Virgin Islands: U.S. Department of Commerce, Climatogra- phy of the United States 60–52, 29 p. Cardona, W.A., 1984, El Yunque mineral prospects, eastern Puerto Rico: Caribbean Journal of Science, v. 20, p. 79–87. Díaz, P.L., Aquino, Zaida, Figueroa-Alamo, Carlos, García, René, and Sánchez, A.V., 2004, Water resources data for Puerto Rico and the U.S. Virgin Islands, water year 2002: U.S. Geological Survey Water-Data Report PR–02–1, 652 p. Dong, Hailiang, Peacor, D.R., and Murphy, S.F., 1998, TEM study of progressive alteration of igneous biotite to kaolinite throughout a weathered soil profile: Geochimica et Cosmo- chimica Acta, v. 62, p. 1881–1887. Ewel, J.J., and Whitmore, J.L., 1973, The ecological life zones of Puerto Rico and the Virgin Islands: U.S. Department of Agriculture Forest Service, Institute of Tropical Forestry, Research Paper ITF–18, Río Piedras, Puerto Rico, 72 p. Fletcher, R.C., Buss, H.L., and Brantley, S.L., 2006, A sphe- roidal weathering model coupling porewater chemistry to soil thickness during steady-state denudation: Earth and Planetary Science Letters, v. 244, p. 444–457. Foster, D.R., Fluet, M., and Boose, E.R., 1999, Human or natural disturbance—Landscape-scale dynamics of the tropical forests of Puerto Rico: Ecological Applications, v. 9, p. 555–572. Franco, P.A., Weaver, P.L., and Eggen-McIntosh, S., 1997, Forest resources of Puerto Rico, 1990: U.S. Department of Agriculture Forest Service, Southern Research Station, Ashville, NC, Resource Bulletin SRS–22, 45 p. González, Grizelle, Zou, Xiaoming, and Borges, Sonia, 1996, Earthworm abundance and species composition in aban- doned tropical croplands—Comparisons of tree plantations and secondary forests: Pedobiologia, v. 40, p. 385–391. Gould, W.A., Alarcón, C., Fevold, B., Jiménez, M.E., Mar- tinuzzi, Sebastian, Potts, G., Quiñones M., Solórzano, M., and Ventosa, E., 2008, The Puerto Rico Gap Analysis Proj- ect, vol. 1, Land cover, vertebrate species distributions, and land stewardship: U.S. Department of Agriculture Forest Service, International Institute of Tropical Forestry, Río Pie- dras, Puerto Rico, General Technical Report IITF–GTR–39, 165 p. Gould, W.A., Martinuzzi, Sebastian, and Parés-Ramos, I.K., 2012, Land use, population dynamics, and land-cover change in eastern Puerto Rico, ch. B in Murphy, S.F., and Stallard, R.F., eds., Water quality and landscape processes of four watersheds in eastern Puerto Rico: U.S. Geological Survey Professional Paper 1789, p. 25–42. Huffaker, Loren, 2002, Soil survey of Caribbean National Forest and Luquillo Experimental Forest, Commonwealth of Puerto Rico: U.S. Department of Agriculture Natural Resources Conservation Service, Washington, D.C., 181 p. Jolly, W.T., Lidiak, E.G., Schellekens, J.H., and Santos, J., 1998, Volcanism, tectonic and stratigraphic correlations in Puerto Rico, in Lidiak, E.G., and Larue, D.K., eds., Tectonics and geochemistry of the northeastern Carib- bean: Geological Society of America Special Paper 322, p. 1–34. Kennaway, T., and Helmer, E.H., 2007, The forest types and ages cleared for land development in Puerto Rico: GIScience and Remote Sensing, v. 44, p. 356–382. Koenig, Nathan, 1953, A comprehensive agricultural program for Puerto Rico: Washington, D.C., U.S. Department of Agriculture and Commonwealth of Puerto Rico, 290 p. Larsen, M.C., 1997, Tropical geomorphology and geo- morphic work—A study of geomorphic processes and sediment and water budgets in montane humid-tropical forested and developed watersheds, Puerto Rico: Boulder, University of Colorado Geography Department, Ph.D. dissertation, 341 p. Larsen, M.C., 2012, Landslides and sediment budgets in four watersheds in eastern Puerto Rico, ch. F in Murphy, S.F., and Stallard, R.F., eds., Water quality and landscape processes of four watersheds in eastern Puerto Rico: U.S. Geological Survey Professional Paper 1789, p. 153–178. Larsen, M.C., Collar, P.D., and Stallard, R.F., 1993, Research plan for the investigation of water, energy, and biogeochem- ical budgets in the Luquillo mountains, Puerto Rico: U.S. Geological Survey Open-File Report 92–150, 19 p. Larsen, M.C., and Santiago-Román, A., 2001, Mass wasting and sediment storage in a small montane watershed—An extreme case of anthropogenic disturbance in the humid tropics, in Dorava, J.M., Palcsak, B.B., Fitzpatrick, F., and Montgomery, D., eds., Geomorphic processes and riverine habitat: American Geophysical Union Monographs, Water Science & Application Series, v. 4, p. 119–138. Larsen, M.C., and Stallard, R.F., 2000, Water, energy, and biogeochemical budgets, Luquillo Mountains, Puerto Rico: U.S. Geological Survey Fact Sheet 163–99, 4 p. 22 Water Quality and Landscape Processes of Four Watersheds in Eastern Puerto Rico Larsen, M.C., Torres-Sánchez, A.J., and Concepción, I.M., 1999, Slopewash, surface runoff, and fine-litter transport in forest and landslide scars in humid-tropical steeplands, Luquillo Experimental Forest, Puerto Rico: Earth Surface Processes and Landforms, v. 24, p. 481–502. Lewis, L.A., 1974, Slow movement of earth under tropical rain forest conditions: Geology, v. 2, p. 9–10. López, T.M., Aide, T.M., and Thomlinson, J.R., 2001, Urban expansion and the loss of prime agricultural lands in Puerto Rico: Ambio, v. 30, p. 49–54. Lugo, A.E., 1996, Caribbean island landscapes—Indicators of the effects of economic growth on the region: Environment and Development Economics, v. 1, p. 128–136. Lugo, A.E., López, T.M., Ramos González, O.M., and Velez, L.L., and others, 2004, Urbanización de los terrenos en la periferia de El Yunque: U.S. Department of Agriculture Forest Service Report WO–66, 29 p. Martinuzzi, Sebastian, Gould, W.A., and Ramos González, O.M., 2007, Land development, land use, and urban sprawl in Puerto Rico—Integrating remote sensing and population cen- sus data: Landscape and Urban Planning, v. 79, p. 288–297. McDowell, W.H., Bowden, W.B., and Asbury, C.E., 1992, Riparian nitrogen dynamics in two geomorphologically distinct tropical rain forest watersheds—Subsurface solute patterns: Biogeochemistry, v. 18, p. 53–75. Meyerhoff, H.A., 1931, The geology of the Fajardo district, Porto Rico: New York Academy of Sciences, Scientific survey of Porto Rico and the Virgin Islands, v. II, pt. 3, p. 201–360. M’Gonigle, J.W., 1978, Geologic map of the Humacao quad- rangle, Puerto Rico: U.S. Geological Survey Miscellaneous Investigations Series Map I–1070, scale 1:20,000. M’Gonigle, J.W., 1979, Geologic map of the Naguabo and part of the Punta Puerca quadrangles, Puerto Rico: U.S. Geological Survey Miscellaneous Investigations Series Map I–1099, scale 1:20,000. Minyard, M.L., Bruns, M.A., Martinez, C.E., Liermann, L.J., Buss, H.L., and Brantley, S.L., 2011, Halloysite nanotubes and bacteria at the saprolite-bedrock interface, Río Icacos watershed, Puerto Rico: Soil Science Society of America Journal, v. 75, p. 348–356. Murphy, S.F., Brantley, S.L., Blum, A.E., White, A.F., and Dong, H., 1998, Chemical weathering in a tropical watershed, Luquillo Mountains, Puerto Rico II—Rate and mechanism of biotite weathering: Geochimica et Cosmo- chimica Acta, v. 62, p. 227–243. Murphy, S.F., and Stallard, R.F., 2012, Hydrology and climate of four watersheds in eastern Puerto Rico, ch. C in Murphy, S.F., and Stallard, R.F., eds., Water quality and landscape processes of four watersheds in eastern Puerto Rico: U.S. Geological Survey Professional Paper 1789, p. 43–84. O’Connor, P.J., Covich, A.P., Scatena, F.N., and Loope, L.L., 2000, Nonindigenous bamboo along headwater streams of the Luquillo Mountains, Puerto Rico—Leaf fall, aquatic leaf decay and patterns of invasion: Journal of Tropical Ecology, v. 16, p. 499–516. Parés-Ramos, I.K., Gould, W.A., and Aide, T.M., 2008, Agri- cultural abandonment, suburban growth, and forest expan- sion in Puerto Rico between 1991 and 2000: Ecology and Society, v. 13, no. 2; p. 1. Last accessed on September 18, 2008 at http://www.ecologyandsociety.org/vol13/iss2/art1/ Pett-Ridge, J.C., Derry, L.A., and Kurtz, A.C., 2009, Sr isotopes as a tracer of weathering processes and dust inputs in a tropical granitoid watershed, Luquillo Moun- tains, Puerto Rico: Geochimica et Cosmochimica Acta, v. 73, p. 25–43. Ramos, O.M., and Lugo, A.E., 1994, Mapa de la vegetación de Puerto Rico: Acta Cientifica, v. 8, p. 63–66. Ramos-Ginés, Orlando, 1999, Estimation of magnitude and frequency of floods for streams in Puerto Rico—New empirical models: U.S. Geological Survey Water-Resources Investigations Report 99–4142, 41 p. Rogers, C.L., Cram, C.M., Pease, M.H., and Tischler, M.S., 1979, Geologic map of the Yabucoa and Punta Tuna quad- rangles, Puerto Rico: U.S. Geological Survey Miscellaneous Investigations Map I–1086, scale 1:20,000. Rudel, T.K., Perez-Lugo, M., and Zichal, H., 2000, When fields revert to forest—Development and spontaneous refor- estation in postwar Puerto Rico: Professional Geographer, v. 52, p. 386–397. Scatena, F.N., 1989, An introduction to the physiography and history of the Bisley experimental watersheds in the Luquillo Mountains of Puerto Rico: U.S. Department of Agriculture General Technical Report SO–72, 22 p. Schellekens, Jaap, Scatena, F.N., Bruijnzeel, L.A., van Dijk, A.I.J.M., Groen, M.M.A., and van Hogezand, R.J.P., 2004, Stormflow generation in a small rainforest catchment in the Luquillo Experimental Forest, Puerto Rico: Hydrological Processes, v. 18, p. 505–530. Schulz, M.S., and White, A.F., 1999, Chemical weathering in a tropical watershed, Luquillo Mountains, Puerto Rico III—Quartz dissolution rates: Geochimica et Cosmochimica Acta, v. 63, p. 337–350. Seiders, V.M., 1971a, Cretaceous and lower Tertiary stratigra- phy of the Gurabo and El Yunque quadrangles, Puerto Rico: U.S. Geological Survey Bulletin 1294–F, 58 p., 1 pl. Physiography, Geology, and Land Cover of Four Watersheds in Eastern Puerto Rico 23 Seiders, V.M., 1971b, Geologic map of the El Yunque Quad- rangle, Puerto Rico. U.S. Geological Survey Miscellaneous Geological Investigations Map I–658, scale 1:20,000. Seiders, V.M., 1971c, Geologic map of the Gurabo quad- rangle, Puerto Rico: U.S. Geological Survey Miscellaneous Geologic Investigations Map I–657, scale 1: 20,000. Simon, Andrew, Larsen, M.C., and Hupp, C.R., 1990, The role of soil processes in determining mechanisms of slope failure and hillslope development in a humid-tropical forest, eastern Puerto Rico, in Kneuper, P.L.K., and McFadden, L.D., eds., Soils and landscape evolution: Geomorphology, v. 3, p. 263–286. Smith, A.L., Schellekens, J.H., and Muriel Díaz, A.-L., 1998, Batholiths as markers of tectonic change in the northeastern Caribbean: Geological Society of America Special Paper, v. 322, p. 99–122. Stallard, R.F., 1995, Relating chemical and physical erosion, in White, A.F., and Brantley, S.L., eds., Chemical weath- ering rates of silicate minerals: Washington, D.C., Mineralogical Society of America Reviews in Mineralogy, v. 31, p. 543–564. Stallard, R.F., 2012, Weathering, landscape equilibrium, and carbon in four watersheds in eastern Puerto Rico, ch. H in Murphy, S.F., and Stallard, R.F., eds., Water quality and landscape processes of four watersheds in eastern Puerto Rico: U.S. Geological Survey Professional Paper 1789, p. 199–248. Stallard, R.F., and Murphy, S.F., 2012, Water quality and mass transport in four watersheds in eastern Puerto Rico, ch. E in Murphy, S.F., and Stallard, R.F., eds., Water quality and landscape processes of four watersheds in eastern Puerto Rico: U.S. Geological Survey Professional Paper 1789, p. 113–152. Turner, B.F., Stallard, R.F., and Brantley, S.L., 2003, Inves- tigation of in situ weathering of quartz diorite bedrock in the Río Icacos Basin, Luquillo Experimental Forest, Puerto Rico: Chemical Geology, v. 202, p. 313–341. U.S. Census Bureau, 2000, United States Census 2000: last accessed on April 2, 2011 at http://www.census.gov/main/ www/cen2000.html U.S. Geological Survey, 1967, El Yunque quadrangle, Puerto Rico: U.S. Geological Survey topographic map, 7.5 minute series, 1 sheet. Wadsworth, F.H., 1950, Notes on the climax forests of Puerto Rico and their destruction and conservation prior to 1990: Caribbean Forester, v. 11, p. 38–56. Wadsworth, F.H., 1970, Review of past research in the Luquillo Mountains of Puerto Rico, in Odum, H.T., and Pigeon, R.F., eds., A tropical rain forest: Springfield, Va., National Technical Information Service, p. B33–B46. Weaver, P.L., 1986, Growth and age of Cyrilla racemiflora L. in montane forests of Puerto Rico: Interciencia, v. 11, p. 221–227. Weaver, P.L., 1989, Rare trees in the Colorado Forest of Puerto Rico’s Luquillo Mountains: Natural Areas Journal, v. 9, p. 169–173. Weaver, P.L., 1990, Succession in the elfin woodland of the Luquillo Mountains of Puerto Rico: Biotropica, v. 22, p. 83–89. Webb, R.M.T., and Gómez-Gómez, F., 1998, Synoptic survey of water quality and bottom sediments, San Juan Bay estu- ary system, Puerto Rico, December 1994 –July 1995: U.S. Geological Survey Water-Resources Investigations Report 97–4144, 48 p. + appendix. White, A.F., and Blum, A.E., 1995, Effects of climate on chemical weathering rates in watersheds: Geochimica et Cosmochimica Acta, v. 59, p. 1729–1747. White, A.F., Blum, A.E., Schulz, M.S., Vivit, D.V., Ston- estrom, D.A., Larsen, M.C., Murphy, S.F., and Eberl, D.D., 1998, Chemical weathering in a tropical watershed, Luquillo Mountains, Puerto Rico I—Long-term versus short-term weathering fluxes: Geochimica et Cosmochimica Acta, v. 62, p. 209–226. Zimmerman, J.K., Aide, T.M., Rosario, M., Serrano, M., and Herrera, L., 1995, Effects of land management and a recent hurricane on forest structure and composition in the Luquillo Experimental Forest, Puerto Rico: Forest Ecology and Management, v. 77, p. 65–76.