Spatial prediction of precipitation with high resolution is a challenging task in regions with strong climate variability and scarce monitoring. For this purpose, the quasi-continuous supply of information from satellite imagery is commonly used to complement in situ data. However, satellite images of precipitation are available at coarse resolutions, and require adequate methods for spatial downscaling and calibration. The objective of this paper is to introduce and evaluate a 2-step spatial downscaling approach for monthly precipitation applied to TRMM 3B43 (from 0.25 • ≈ 27 km to 5 km resolution), resulting in 5 downscaled products for the period 01-2001/12-2011. The methodology was evaluated in 3 contrasting climatic regions of Ecuador. In step 1, bilinear resampling was applied over TRMM, and used as a reference product. The second step introduces further variability, and consists of four alternative gauge-satellite merging methods: (1) regression with in situ stations, (2) regression kriging with in situ stations, (3) regression with in situ stations and auxiliary variables, and (4) regression kriging with in situ stations and auxiliary variables. The first 2 methods only use the resampled TRMM data set as an independent variable. The last 2 methods enrich these models with auxiliary environmental factors, incorporating atmospheric and land variables. The results showed that no product outperforms the others in every region. In general, the methods with residual kriging correction outperformed the regression models. Regression kriging with situ data provided the best representation in the Coast, while regression kriging with in situ and auxiliary data generated the best results in the Andes. In the Amazon, no product outperformed the resampled TRMM images, probably due to the low density of in situ stations. These results are relevant to enhance satellite precipitation, depending on the availability of in situ data, auxiliary satellite variables and the particularities of the climatic regions.