Abstract: Travel time distributions (TTDs) provide an effective way to describe the transport and mixing processes of water parcels in a subsurface hydrological system. A major challenge in characterizing catchment TTD is quantifying the travel times in deep groundwater and its contribution to the streamflow TTD. Here, we develop and test a novel modeling framework for an integrated assessment of catchment scale TTDs through explicit representation of 3D-groundwater dynamics. The proposed framework is based on the linkage between a flux tracking scheme with the surface hydrologic model (mHM) for the soil-water compartment and a particle tracking scheme with the 3D-groundwater model OpenGeoSys (OGS) for the groundwater compartment. This linkage provides us with the ability to simulate the spatial and temporal dynamics of TTDs in these different hydrological compartments from grid scale to regional scale. We apply this framework in the Nägelstedt catchment in central Germany. Simulation results reveal that both shape and scale of grid-scale groundwater TTDs are spatially heterogeneous, which are strongly dependent on the topography and aquifer structure. The component-wise analysis of catchment TTD shows a time-dependent sensitivity of transport processes in soil zone and groundwater to driving meteorological forcing. Catchment TTD exhibits a power-law shape and fractal behavior. The predictive uncertainty in catchment mean travel time is dominated by the uncertainty in the deep groundwater rather than that in the soil zone. Catchment mean travel time is severely biased by a marginal error in groundwater characterization. Accordingly, we recommend to use multiple summary statistics to minimize the predictive uncertainty introduced by the tailing behavior of catchment TTD.