Abstract. Global gridded maps (a.k.a. Level 3 products) of Earth system properties observed by satellites are central to understanding the spatiotemporal variability of these properties. They also typically serve either as inputs into biogeochemical models or as independent data for evaluating such models. Spatial binning is a common method for generating contiguous maps, but this approach results in a loss of information, especially when the measurement noise is low relative to the degree of spatiotemporal variability. Such "binned" fields typically also lack a quantitative measure of uncertainty.

Geostatistical mapping has previously been shown to make higher spatiotemporal resolution maps possible, and also provides a measure uncertainty associated with the gridded products. This study proposes a flexible moving window block kriging method that can be used as a tool for creating high spatiotemporal resolution maps from satellite data. It relies only on the assumption that the observed physical quantity exhibits spatial correlation that can be inferred from the observations. The method has several innovations relative to previously applied methods: (1) it provides flexibility in the spatial resolution of the contiguous maps, (2) it is applicable for physical quantities with varying spatiotemporal coverage (i.e., density of measurements) by utilizing a more general and versatile data sampling approach, and (3) it provides rigorous assessments of the uncertainty associated with the gridded products. The method is demonstrated by creating Level 3 products from observations of column-integrated carbon dioxide (XCO₂) from the GOSAT (Greenhouse Gases Observing Satellite) satellite, and solar induced fluorescence (SIF) from the GOME-2 (Global Ozone Monitoring Experiment-2) instrument.