Modeling rock weathering in small watersheds

•Couple rock weathering with groundwater flow and landscape evolution models.•Assess plagioclase weathering rates at the spring watershed scale.•Describe lab and field rates in relation to aquifer hydraulic diffusivity and time. Many mountainous watersheds are conceived as aquifer media where multiple groundwater flow systems have developed (Tóth, 1963), and as bimodal landscapes where differential weathering of bare and soil-mantled rock has occurred (Wahrhaftig, 1965). The results of a weathering algorithm (Pacheco and Van der Weijden, 2012a, 2014), which integrates topographic, hydrologic, rock structure and chemical data to calculate weathering rates at the watershed scale, validated the conceptual models in the River Sordo basin, a small watershed located in the Marão cordillera (North of Portugal). The coupling of weathering, groundwater flow and landscape evolution analyses, as accomplished in this study, is innovative and represents a remarkable achievement towards regionalization of rock weathering at the watershed scale. The River Sordo basin occupies an area of approximately 51.2km2 and was shaped on granite and metassediment terrains between the altitudes 185–1300m. The groundwater flow system is composed of recharge areas located at elevations >700m, identified on the basis of δ18O data. Discharge cells comprehend terminations of local, intermediate and regional flow systems, identified on the basis of spring density patterns, infiltration depth estimates based on 87Sr/86Sr data, and spatial distributions of groundwater pH and natural mineralization. Intermediate and regional flow systems, defined where infiltration depths >125m, develop solely along the contact zone between granites and metassediments, because fractures in this region are profound and their density is very large. Weathering is accelerated where rocks are covered by thick soils, being five times faster relative to sectors of the basin where rocks are covered by thin soils. Differential weathering of bare and soil-mantled rock is also revealed by the spatial distribution of calculated aquifer hydraulic diffusivities and groundwater travel times.