In environmental applications it is often of interest to predict the rates at which contaminant mass is discharged at a given cross section of streams and rivers. We present a Lagrangian methodology for evaluating tracer discharge (mass per unit time) at specified control cross sections (CCS) of streams. The key transport processes included in the analysis are advection, degradation/decay, and kinetically controlled mass transfer in storage zones and in bed sediment. The transport in the bed sediment is described as a diffusion process, where the tracer may sorb onto the sediment. We have derived a general solution for tracer discharge in the Laplace domain wherefrom temporal moments are computed. The derived solutions may account for deterministic changes in morphological characteristics along stretches of streams. The results are illustrated for zeroth and first two moments where we show the combined effect of advection, degradation, physicochemical mass transfer, and morphology. For illustrative purposes, we assume morphology to change downstream following power laws suggested by Langbein <1947> and Leopold and Maddock <1953>. The moments depend nonlinearly on the downstream distance, following power laws that reflect the power laws for the hydraulic geometry. We define two main dimensionless parameters, namely, kinetic storage parameter &agr;* and bed parameter M, that control the amount of tracer mass eventually discharged at any given CCS. For M≈≥0.3 the most dominant mechanism that controls the amount of ultimately discharged tracer mass is the exchange with the bed sediment. Once estimated from field data for stretches of specified streams, the dimensionless parameters can be used in the derived expressions for predictive purposes. ¿ 2000 American Geophysical Union