A watershed (also drainage basin, river basin, or catchment) is defined as "… the area that topographically appears to contribute all the water that passes through a specified cross section of a stream (the outlet)" (Dingman 2015). In this chapter, I choose to use the term "watershed" as it is a broadly used one; it should be understood more as small watersheds or catchments. Watersheds are the fundamental units that support river networks, the blood vessels at Earth's surface ultimately draining into the ocean. Watersheds are complex hydro-biogeochemical reactors. They receive water, mass, and energy, transport them to distinct compartments, and transform them into various forms (Fig.1). Hydrological processes partition precipitation to the atmosphere, to the ground surface, and to the subsurface, eventually entering streams. Similarly, plants translate sunlight, water, CO 2 , and nutrients into organic matters (leaves, stems, roots) that fall and deposit in soil. As water routes through soils, it interacts with roots, microbes, and reactive gases (i.e., CO 2 and O 2), releases solutes, and ultimately transporting them out of watersheds. The water flow (discharge) and solute concentrations measured at stream outlets therefore reflect convoluted signature of ecohydrological and biogeochemical coupling. The process interactions and feedbacks are dictated not only by hydroclimatic forcing but also the architecture of watersheds, in particular the above-ground characteristics such as land cover and surface topography, as well as the below-ground structure including soil depth, soil type, geology, and root architecture. The external forcing and internal idiosyncrasies dictate the magnitude, timing, and spatial distribution of water flow and chemistry (Chorover et al. 2011; Brooks et al. 2015), giving rise to non-linear emergent behaviors that are unique of watershed reactive transport processes. Understanding feedbacks and non-linearity requires process-based models that integrate multiple interacting processes. Such integration however does not come easily with traditional boundaries of hydrology versus biogeochemistry relevant disciplines. As will be discussed later in this chapter, relevant model development has been advancing along two separate lines: hydrology models that solve for water storage and fluxes at the watershed scale and beyond (Fatichi et al. 2016), and reactive transport models (RTMs) that center on aqueous and solid concentration changes arising from transport and multi-component biogeochemical reactions typically in "closed" groundwater systems without much interactions with "open" watersheds directly receiving precipitation and sun light (Steefel et al. 2015; Li et al. 2017b). This comes along with a history of hydrologists often trained as physicists studying fluid mechanics, and biogeochemists typically grow up as geologists, chemists, or environmental engineers. There are however considerable needs to reach beyond disciplinary boundaries and integrate the two lines to develop watershed reacti...