We developed a large-area preprocessing framework for multisensor Landsat data, capable of processing large data volumes. Cloud and cloud shadow detection is performed by a modified Fmask code. Surface reflectance is inferred from Tanré's formulation of the radiative transfer, including adjacency effect correction. A precompiled MODIS water vapor database provides daily or climatological fallback estimates. Aerosol optical depth (AOD) is estimated over dark objects (DOs) that are identified in a combined database and image-based approach, where information on their temporal persistency is utilized. AOD is inferred with consideration of the actual target reflectance and background contamination effect. In case of absent DOs in bright scenes, a fallback approach with a modeled AOD climatology is used instead. Topographic normalization is performed by a modified C-correction. The data are projected into a single coordinate system and are organized in a gridded data structure for simplified pixel-based access. We based the assessment of the produced data set on an exhaustive analysis of overlapping pixels: 98.8% of the redundant overlaps are in the range of the expected ±2.5% overall radiometric algorithm accuracy. AOD is in very good agreement with Aerosol Robotic Network sunphotometer data (R ² : 0.72 to 0.79, low intercepts, and slopes near unity). The uncertainty in using the water vapor fallback climatology is approximately ±2.8% for the TM SWIR1 band in the wet season. The topographic correction was considered successful by an investigation of the nonrelationship between the illumination angle and the corrected radiance.