Search

Several data assimilation (DA) approaches exist to generate consistent and continuous precipitation fields valuable for hydrometeorological applications and land data assimilation. Usually, DA is based on either static or dynamic approaches. Static methods rely on deterministic forecasts to estimate background error covariance matrices, whereas dynamic approaches use ensemble forecasts. Associating the two methods is known as hybrid DA, and it has proven beneficial for different applications as it combines the advantages of both approaches. The present study intends to explore hybrid DA for the 6 h Canadian Precipitation Analysis (CaPA). Based on optimal interpolation (OI), CaPA blends forecasts and observations from surface stations and ground-based radar datasets to provide precipitation fields over the North American domain. The application of hybrid DA to CaPA consisted of finding the optimal linear combination between (i) an OI based on the Regional Deterministic Prediction System (RDPS) and (ii) an ensemble Kalman filter (EnKF) based on the 20-member Regional Ensemble Prediction System (REPS). The results confirmed the known effectiveness of the hybrid approach when low-density observation networks are assimilated. Indeed, the experiments conducted for the summer without radar datasets and for the winter (characterized by very few observations in CaPA) showed that attributing a relatively high weight to the EnKF (50 % and 70 % for summer and winter, respectively) resulted in better analysis skill and a reduction in false alarms compared with the OI method. A deterioration in the moderate- to high-intensity precipitation bias was, however, observed during summer. Reducing the weight attributed to the EnKF to 30 % alleviated the bias deterioration while improving skill compared with the OI-based CaPA.

Environment and Climate Change Canada has initiated the production of a 1980–2018, 10 km, North American precipitation and surface reanalysis. ERA-Interim is used to initialize the Global Deterministic Reforecast System (GDRS) at a 39 km resolution. Its output is then dynamically downscaled to 10 km by the Regional Deterministic Reforecast System (RDRS). Coupled with the RDRS, the Canadian Land Data Assimilation System (CaLDAS) and Precipitation Analysis (CaPA) are used to produce surface and precipitation analyses. All systems used are close to operational model versions and configurations. In this study, a 7-year sample of the reanalysis (2011–2017) is evaluated. Verification results show that the skill of the RDRS is stable over time and equivalent to that of the current operational system. The impact of the coupling between RDRS and CaLDAS is explored using an early version of the reanalysis system which was run at 15 km resolution for the period 2010–2014, with and without the use of CaLDAS. Significant improvements are observed with CaLDAS in the lower troposphere and surface layer, especially for the 850 hPa dew point and absolute temperatures in summer. Precipitation is further improved through an offline precipitation analysis which allows the assimilation of additional observations of 24 h precipitation totals. The final dataset should be of particular interest for hydrological applications focusing on transboundary and northern watersheds, where existing products often show discontinuities at the border and assimilate very few – if any – precipitation observations.

This study assesses the potential of satellite imagery of vertically integrated columns of dry-air mole fractions of CO2 (XCO2) to constrain the emissions from cities and power plants (called emission clumps) over the whole globe during 1 year. The imagery is simulated for one imager of the Copernicus mission on Anthropogenic Carbon Dioxide Monitoring (CO2M) planned by the European Space Agency and the European Commission. The width of the swath of the CO2M instruments is about 300 km and the ground horizontal resolution is about 2 km resolution. A Plume Monitoring Inversion Framework (PMIF) is developed, relying on a Gaussian plume model to simulate the XCO2 plumes of each emission clump and on a combination of overlapping assimilation windows to solve for the inversion problem. The inversion solves for the 3 h mean emissions (during 08:30–11:30 local time) before satellite overpasses and for the mean emissions during other hours of the day (over the aggregation between 00:00–08:30 and 11:30–00:00) for each clump and for the 366 d of the year. Our analysis focuses on the derivation of the uncertainty in the inversion estimates (the “posterior uncertainty”) of the clump emissions. A comparison of the results obtained with PMIF and those from a previous study using a complex 3-D Eulerian transport model for a single city (Paris) shows that the PMIF system provides the correct order of magnitude for the uncertainty reduction of emission estimates (i.e., the relative difference between the prior and posterior uncertainties). Beyond the one city or few large cities studied by previous studies, our results provide, for the first time, the global statistics of the uncertainty reduction of emissions for the full range of global clumps (differing in emission rate and spread, and distance from other major clumps) and meteorological conditions. We show that only the clumps with an annual emission budget higher than 2 MtC yr−1 can potentially have their emissions between 08:30 and 11:30 constrained with a posterior uncertainty smaller than 20 % for more than 10 times within 1 year (ignoring the potential to cross or extrapolate information between 08:30–11:30 time windows on different days). The PMIF inversion results are also aggregated in time to investigate the potential of CO2M observations to constrain daily and annual emissions, relying on the extrapolation of information obtained for 08:30–11:30 time windows during days when clouds and aerosols do not mask the plumes, based on various assumptions regarding the temporal auto-correlations of the uncertainties in the emission estimates that are used as a prior knowledge in the Bayesian framework of PMIF. We show that the posterior uncertainties of daily and annual emissions are highly dependent on these temporal auto-correlations, stressing the need for systematic assessment of the sources of uncertainty in the spatiotemporally resolved emission inventories used as prior estimates in the inversions. We highlight the difficulty in constraining the total budget of CO2 emissions from all the cities and power plants within a country or over the globe with satellite XCO2 measurements only, and calls for integrated inversion systems that exploit multiple types of measurements.

A large fraction of fossil fuel CO2 emissions emanate from “hotspots”, such as cities (where direct CO2 emissions related to fossil fuel combustion in transport, residential, commercial sectors, etc., excluding emissions from electricity-producing power plants, occur), isolated power plants, and manufacturing facilities, which cover a small fraction of the land surface. The coverage of all high-emitting cities and point sources across the globe by bottom-up inventories is far from complete, and for most of those covered, the uncertainties in CO2 emission estimates in bottom-up inventories are too large to allow continuous and rigorous assessment of emission changes (Gurney et al., 2019). Space-borne imagery of atmospheric CO2 has the potential to provide independent estimates of CO2 emissions from hotspots. But first, what a hotspot is needs to be defined for the purpose of satellite observations. The proposed space-borne imagers with global coverage planned for the coming decade have a pixel size on the order of a few square kilometers and a XCO2 accuracy and precision of  ppm for individual measurements of vertically integrated columns of dry-air mole fractions of CO2 (XCO2). This resolution and precision is insufficient to provide a cartography of emissions for each individual pixel. Rather, the integrated emission of diffuse emitting areas and intense point sources is sought. In this study, we characterize area and point fossil fuel CO2 emitting sources which generate coherent XCO2 plumes that may be observed from space. We characterize these emitting sources around the globe and they are referred to as “emission clumps” hereafter. An algorithm is proposed to identify emission clumps worldwide, based on the ODIAC global high-resolution 1 km fossil fuel emission data product. The clump algorithm selects the major urban areas from a GIS (geographic information system) file and two emission thresholds. The selected urban areas and a high emission threshold are used to identify clump cores such as inner city areas or large power plants. A low threshold and a random walker (RW) scheme are then used to aggregate all grid cells contiguous to cores in order to define a single clump. With our definition of the thresholds, which are appropriate for a space imagery with 0.5 ppm precision for a single XCO2 measurement, a total of 11 314 individual clumps, with 5088 area clumps, and 6226 point-source clumps (power plants) are identified. These clumps contribute 72 % of the global fossil fuel CO2 emissions according to the ODIAC inventory. The emission clumps is a new tool for comparing fossil fuel CO2 emissions from different inventories and objectively identifying emitting areas that have a potential to be detected by future global satellite imagery of XCO2. The emission clump data product is distributed from https://doi.org/10.6084/m9.figshare.7217726.v1.