Do not fill out this information. 

On this particular page, the Initial Text and Image fields for this template are not used (inactive).  A Carousel Slideshow is displayed instead. 

To edit the Carousel Slideshow, Choose Edit Existing >  Carousel Slides.  Then navigate to your CZO's entry entitled "Research"

Research in the Luquillo Critical Zone Observatory

The Luquillo Critical Zone Observatory (LCZO) focuses on how water balances and mass fluxes differ in landscapes with contrasting lithology across a range of climatic and vegetation zones.  Parent material has been recognized as a state factor for landscape and soil development for over a century.  However, the influences of lithology on denudation, hydrologic routing, and the elemental cycles and flows is still poorly constrained in most studies.  The LCZO will use the natural laboratory of the Luquillo Mountains to determine how critical zone processes differ in watersheds underlain by quartz diorite (GD ) and volcaniclastic (VC) bedrock (Figure 1). Specifically, the LCZO will address how the bedrock is coupled to or decoupled from the hydrologic, geochemical and biogeochemical cycles in these watersheds.  A set of interrelated hypotheses (Figure 2, Table 1) will be tested taking advantage of the cohesive network of sampling sites (Table 2) and a unified data management system. This approach facilitates controlled comparisons across bedrock (GD, VC) types, landscape position (ridge, hillslope, riparian), depths (surface to bedrock), forest type (Tabonuco, Colorado, Cloud) and distance from the headwaters (upland to coastal).

Figure 1; Luquillo Critical Zone focus watersheds. With stream gages (water droplets). 

Short term responses (e.g. days to seasons) to be investigated include the influence of lithology on biogeochemical responses to storms and droughts, and the influence of lithology on hydrologic routing and watershed hydrologic budgets. Recent research indicates there are greater seasonal variations in the chemistry of atmospheric inputs and soils processes than previously expected (Heartsill et al 2007, Shanley 2008, Scholl in press).  These responses will be quantified using an improved network of weather stations, nested stream gages and multi-investigator event sampling campaigns.  Responses at intermediate time scales (e.g. decades to centuries) will be evaluated by quantifying stream channel dynamics and examining the conditions in the past recorded in the sediments of the floodplains and coastal zones.  The long-term influences of lithology on sediment generation, hillslope and landscape development will be evaluated in studies of bedrock weathering rates and rates of saprolite and soil formation.

Figure 2: Conceptual framework and relationships between the hypothesis and sampling nodes for the Luquillo Critical Zone Observatory.  Infrastructure is in red while hypothesis are in green.

Hypothesis 1: The rate of saprolite advance varies with regolith thickness and landscape position  and is fastest in GD valleys and slowest on VC ridges.  Over large areas, the rate of saprolite advance will equal the rate of denudation and can be predicted from bedrock chemistry, porefluids, and physical rock properties (Brantley, Buss, White, Heimsath, Willenbring,).
The rate of saprolite production increases with regolith thickness, which depends on landscape position and bedrock type.  The rates should be fastest in GD valleys and slowest on VC ridges. Over large areas, the rate of weathering will equal the rate of denudation resulting steady state landscape profiles.  The magnitude of the rates can be predicted from bedrock chemistry, porefluids composition, and physical rock properties

Hypothesis 2: In surface soils, chemical transformations of atmospheric inputs are decoupled from bedrock lithology and influenced by soil carbon, surface redox, and plant nutrient cycling.  Biotic influences on soil biogeochemistry decrease with storm intensity and soil depth and are greatest in surface soils of the VC during low intensity rainfalls (Johnson, Shanley, Silver, Scatena, UPR).

Hypothesis 3: The routing of water through the critical zone is determined by bedrock lithology. The residence time will be longer in areas underlain by GD and shorter in areas underlain by VC resulting in differences in water chemsitry.  However these differences will decrease with storm intensity and duration (Scholl, Scatena, Shanley, McDowell).

Hypothesis 4: Over seasonal time scales, iron reduction and the associated CO2 production will vary across the landscape based on differences in soil redox and organic matter content. Iron reduction will be greatest in riparian VC surface soils and lowest at depth on stable GD ridge tops.  At larger spatial scales and longer time scales deep weathering rates and iron cycling are closely linked to the frequency of low redox events and carbon availability (Silver, Brantley, Plante)

Hypothesis 5: The morphology and soil biogeochemistry of riparian and colluvial deposits varies systematically with lithology and distance from the headwaters.  In contrast, the vegetation and soil organic matter chemistry vary systematically with rainfall and temperature (McDowell, Plante, Silver, Scatena, Jerolmack).

Hypothesis 6: Surface erosion is the dominant source of sediment to the stream channel network across the entire landscape. While surface erosion is associated with rainfalls of moderate intensity in VC watersheds, in the GD watersheds sediment originates from the landslides associated with high intensity precipitation events (Jerolmack, Horton, Willenbring, Scatena, Shanley).

Hypothesis 7: The temporal resolution of the climatic disturbances and land use change recorded in coastal and fluvial sediments differs between lithologies.  The GD watersheds will have a higher resolution record than corresponding environments in the VC watersheds because of …(Horton, Willenbring, Jerolmack, Heimsath, Scatena).

Figure 1; Luquillo Critical Zone focus watersheds. With stream gages (water droplets). 

Figure 2: Conceptual framework and relationships between the hypothesis and sampling nodes for the Luquillo Critical Zone Observatory.  Infrastructure is in red while hypothesis are in green.

Research News


CZOs at AGU 2018

19 Nov 2018 - The 2018 AGU Fall Meeting will be held December 10-14 in Washington, D.C.


NSF Discovery articles focus on the CZOs.

10 May 2018 - The Discoveries section of the National Science Foundation's website on Critical Zone Observatories (CZOs).


NSF Discovery: Can an ancient ocean shoreline set the stage for a tropical forest of today?

13 Jul 2017 - Researchers at NSF Critical Zone Observatory and Long-Term Ecological Research sites are finding out.


Extreme Rainfall Doesn’t Always Mean Extreme Erosion, Penn Study Finds

10 May 2016 - A new study suggests suggest that the self-organization of bedrock river channels blunts the impact of extreme rainfall events.


Critical Zone Profile - JAIVIME EVARISTO (ecohydrologist, PhD student)

23 Mar 2016 - PhD student, Global Institute for Water Security, University of Saskatchewan


Critical Zone Profiles - Meet the people doing CZO science (Luquillo CZO)

15 Nov 2015 - Get a sense of the people and the work. Several members of the Luquillo CZO are profiled here, including students and professors.

Luquillo CZO SAVI International Scholars Program recipients anounced

10 Apr 2018 - Two Luquillo CZO researchers have been named SAVI International Scholars Program recipients,   Emma Harrison from the University of...

By river, ocean, or wind, rocks round the same way

04 Apr 2018 - The milky, smooth texture of beach glass evokes a history of turbulent transport, rough edges abraded away to produce curves. The same smooth...

Discussion paper open till April 5th in Biogeosciences

06 Mar 2018 - A discussion is open for the manuscript, "Elevating the biogeosciences within environmental research networks".

Teachers Data Jam Workshop

22 Jan 2018 - On Saturday, January 13, we successfully held our third Teacher Data Jam Workshop of the Luquillo LTER Schoolyard Program. On this occasion,...

White paper of 2017 CZ Science Meeting in Arlington, VA

01 Jan 2018 - New Opportunities for Critical Zone Science Following the June 2017 Arlington Meeting for Critical Zone Science (hosted by CZO), a white booklet...

More News >

Example Publications


Mercury sourcing and sequestration in weathering profiles at six Critical Zone Observatories. Richardson, Justin B., Arnulfo A. Aguirre, Heather L. Buss, A. Toby O'Geen, Xin Gu, Daniella M. Rempe, and Daniel deB. Richter (2018): Global Biogeochemical Cycles, doi: 10.1029/2018GB005974 Cross-CZO National


Ideas and perspectives: Strengthening the biogeosciences in environmental research networks. Richter, D.D., S.A. Billings, P.M. Groffman, E.F. Kelly, K.A. Lohse, W.H. McDowell, T.S. White, S. Anderson, D.D. Baldocchi, S. Banwart, S. Brantley, J.J. Braun, Z.S. Brecheisen, C.W. Cook, H.E. Hartnett, S.E. Hobbie, J. Gaillardet, E. Jobbagy, H.F. Jungkunst, C.E. Kazanski, J. Krishnaswamy, D. Markewitz, K. O’Neill, C.S. Riebe, P. Schroeder, C. Siebe, W.L. Silver, A. Thompson, A. Verhoef, G. Zhang (2018): Biogeosciences 15: 4815-4832 Cross-CZO National


A net ecosystem carbon budget for snow dominated forested headwater catchments: linking water and carbon fluxes to critical zone carbon storage. Perdrial J., Brooks P.D., Swetnam T., Lohse K.A., Rasmussen C., Litvak M., Harpold A.A., Zapata-Rios X., Broxton P., Mitra B., Meixner M., Condon K., Huckle D., Stielstra C., Vázquez-Ortega A., Lybrand R., Holleran M., Orem C., Pelletier J., Chorover J. (2018): Biogeochemistry 138{3): 225–243 Cross-CZO


Beyond clay: towards an improved set of variables for predicting soil organic matter content. Rasmussen C., Heckman K., Wieder W.R., Keiluweit M., Lawrence C.R., Berhe A.A., Blankinship J.C., Crow S.E., Druhan J.L., Hicks Pries C.E., Marin-Spiotta E., Plante A.F., Schädel C., Schimel J.P., Sierra C.A., Thompson A., Wagai R. (2018): Biogeochemistry 137(3): 297–306 Cross-CZO


Concentration-Discharge Relations in the Critical Zone: Implications for Resolving Critical Zone Structure, Function and Evolution. Chorover, J., Derry, L. A., McDowell, W. H. (2017): Water Resources Research 53(11): 8654–8659 Cross-CZO National


Growing new generations of critical zone scientists. Wymore, Adam S., Nicole R. West, Kate Maher, Pamela L. Sullivan, Adrian Harpold, Diana Karwan, Jill A. Marshall, Julia Perdrial, Daniella M. Rempe and Lin Ma (2017): Earth Surface Processes and Landforms 42 (14): 2498-2502 Cross-CZO National


Geochemical evolution of the Critical Zone across variable time scales informs concentration-discharge relationships: Jemez River Basin Critical Zone Observatory. McIntosh J.C., Schaumberg C., Perdrial J., Harpold A., Vázquez-Ortega A., Rasmussen C., Vinson D., Zapata-Rios X., Brooks P.D., Meixner T., Pelletier J., Derry L., Chorover J. (2017): Water Resources Research 53(5): 4169–4196 Cross-CZO National


Designing a network of critical zone observatories to explore the living skin of the terrestrial Earth. Brantley, S.L., McDowell, W.H., Dietrich, W.E., White, T.S., Kumar, P., Anderson, S., Chorover, J., Lohse, K.A., Bales, R.C., Richter, D., Grant, G., and Gaillardet, J. (2017): Earth Surface Dynamics, 5, 841–860 Cross-CZO National


Reviews and syntheses: on the roles trees play in building and plumbing the critical zone. Brantley, Susan L., David M. Eissenstat, Jill A. Marshall, Sarah E. Godsey, Zsuzsanna Balogh-Brunstad, Diana L. Karwan, Shirley A. Papuga, Joshua Roering, Todd E. Dawson, Jaivime Evaristo, Oliver Chadwick, Jeffrey J. McDonnell, Kathleen C. Weathers (2017): Biogeosciences, 14, 5115-5142 Cross-CZO National


Controls on deep critical zone architecture: a historical review and four testable hypotheses. Riebe, C. S., Hahm, W. J., Brantley, S. L. (2017): Earth Surface Processes and Landforms, 42 (1): 128–156 Cross-CZO National


Special issue of The Earth Scientist about the Critical Zone and the US NSF Critical Zone Observatory (CZO) program. CZO Education/Outreach team (2016): The Earth Scientist, Volume XXXII, Issue 3, Fall 2016 Cross-CZO National


Variation of organic matter quantity and quality in streams at Critical Zone Observatory watersheds. Miller, Matthew P., Boyer, Elizabeth W., McKnight, Diane M., Brown, Michael G., Gabor, Rachel S., Hunsaker, Carolyn T., Iavorivska, Lidiia, Inamdar, Shreeram, Johnson, Dale W., Kaplan, Louis A., Lin, Henry, McDowell, William H., Perdrial, Julia N. (2016): Water Resources Research, 52 (10): 8202–8216 Cross-CZO


Reconstructing the transport history of pebbles on Mars. Szabó T., Domokos, G., Grotzinger, J.P., Jerolmack, D.J. (2015): Nature Communications


A field comparison of multiple techniques to quantify groundwater–surface-water interactions. González-Pinzón, R., Ward, A., Hatch, C., Wlostowski, A., Singha, K., Gooseff, M., Haggerty, R., Harvey, J., Cirpka, O., and Brock, J. (2015): Freshwater Science, 34(1) Cross-CZO


Analytical model for flow duration curves in seasonally dry climates. Müller, M.F., Dralle, D.N., and Thompson, S.E. (2014): Water Resources Research 50(7): 5510-5531. Cross-CZO

Interdisciplinary Teaching about Earth and the Environment for a Sustainable Future Book Chapter: Implementing and Assessing InTeGrate Critical Zone Science Materials in an Undergraduate Geoscience Program. Ashlee Dere, Carol Engelmann, Timothy White, Adam Wymore, Adam Hoffman, James Washburne, Martha Conklin (2018): Springer Cross-CZO National

Fire, flood, and drought—Extreme climate events alter flow paths and stream chemistry. Murphy, S.F., McCleskey, R.B., Martin, D.A., Writer, J.H., and Ebel, B.A. (2018): Journal of Geophysical Research: Biogeosciences, 123 Cross-CZO

Potential for iron reduction increases with rainfall in montane basaltic soils of Hawaii. Hodges, C., E. King, J. Pett-Ridge, A. Thompson (2018): Soil Science Society of America Journal 82:176-185 Cross-CZO

Influence of pO2 on iron redox cycling and anaerobic organic carbon mineralization in a humid tropical forest soil. Chen, Chunmei, Christof Meile, Jared Wilmoth, Diego Barcellos, and Aaron Thompson (2018): Environmental Science & Technology 52 (14): 7709-7719 Cross-CZO

Hurricane Maria in the US Caribbean: Disturbance Forces, Variation of Effects, and Implications for Future Storms. Van Beusekom, A.E. Álvarez-Berríos N.L. Gould W.A., Quiñones M. González, G. (2018): Remote sensing

Order from disorder: do soil organic matter composition and turnover co-vary with iron phase crystallinity?. Steven, J. Hall, Asmeret A. Berhe, Aaron Thompson. (2018): Biogeochemistry

Litterfall Production Prior to and during Hurricanes Irma and Maria in Four Puerto Rican Forests. Xianbin Liu, Xiucheng Zeng, Xiaoming Zou, Grizelle González, Chao Wang, Si Yang (2018): Forests

More Publications >