Search

Several configurations of the Canadian Precipitation Analysis system (CaPA) currently produce precipitation analyses at Environment and Climate Change Canada (ECCC). To improve CaPA's performance during the winter season, the impact of assimilating the IMERG V06 product (IMERG: Integrated Multi-satellitE Retrievals for GPM—Global Precipitation Measurement mission) into CaPA is examined in this study. Tests are conducted with CaPA's 10 km deterministic version, evaluated over Canada and the northern part of the United States (USA). Maps from a case study show that IMERG plays a contradictory role in the production of CaPA's precipitation analyses for a synoptic-scale winter storm over North America's eastern coast. While its contribution appears to be physically correct over southern portions of the meteorological system, and early in its intensification phase, IMERG displays unrealistic spatial structures over land later in the system's life cycle when it is located over northern (colder) areas. Objective evaluation of CaPA's analyses when IMERG is assimilated without any restrictions shows an overall decrease in precipitation, which has a mixed effect (positive and negative) on the bias indicators. But IMERG's influence on the Equitable Threat Score (ETS), a measure of CaPA's analyses accuracy, is clearly negative. Using IMERG's quality index (QI) to filter out areas where it is less accurate improves CaPA's objective evaluation, leading to better ETS versus the control experiment in which no IMERG data are assimilated. Several diagnostics provide insight into the nature of IMERG's contribution to CaPA. For the most successful configuration, with a QI threshold of 0.3, IMERG's impact is mostly found in the warmer parts of the domain, i.e., in northern US states and in British Columbia. Spatial means of the temporal sums of absolute differences between CaPA's analyses with and without IMERG indicate that this product also contributes meaningfully over land areas covered by snow, and areas where air temperature is below −2 °C (where precipitation is assumed to be in solid phase). [ABSTRACT FROM AUTHOR]

The Canadian Precipitation Analysis System (CaPA) is an operational system that uses a combination of weather gauge and ground-based radar measurements together with short-term forecasts from a numerical weather model to provide near-real-time estimates of 6- and 24-h precipitation amounts. During the winter season, many gauge measurements are rejected by the CaPA quality control process because of the wind-induced undercatch for solid precipitation. The goal of this study is to improve the precipitation estimates over central Canada during the winter seasons from 2019 to 2022. Two approaches were tested. First, the quality control procedure in CaPA has been relaxed to increase the number of surface observations assimilated. Second, the automatic solid precipitation measurements were adjusted using a universal transfer function to compensate for the undercatch problem. Although increasing the wind speed threshold resulted in lower amounts and worse biases in frequency, the overall precipitation estimates are improved as the equitable threat score is improved because of a substantial decrease in the false alarm ratio, which compensates the degradation of the probability of detection. The increase of solid precipitation amounts using a transfer function improves the biases in both frequency and amounts and the probability of detection for all precipitation thresholds. However, the false alarm ratio deteriorates for large thresholds. The statistics vary from year to year, but an overall improvement is demonstrated by increasing the number of stations and adjusting the solid precipitation amounts for wind speed undercatch. [ABSTRACT FROM AUTHOR]

Consistent and continuous fields provided by precipitation analyses are valuable for hydrometeorological applications and land data assimilation modeling, among others. Providing uncertainty estimates is a logical step in the analysis development, and a consistent approach to reach this objective is the production of an ensemble analysis. In the present study, a 6-h High-Resolution Ensemble Precipitation Analysis (HREPA) was developed for the domain covering Canada and the northern part of the contiguous United States. The data assimilation system is the same as the Canadian Precipitation Analysis (CaPA) and is based on optimal interpolation (OI). Precipitation from the Canadian national 2.5-km atmospheric prediction system constitutes the background field of the analysis, while at-site records and radar quantitative precipitation estimates (QPE) compose the observation datasets. By using stochastic perturbations, multiple observations and background field random realizations were generated to subsequently feed the data assimilation system and provide 24 HREPA members plus one control run. Based on one summer and one winter experiment, HREPA capabilities in terms of bias and skill were verified against at-site observations for different climatic regions. The results indicated HREPA's reliability and skill for almost all types of precipitation events in winter, and for precipitation of medium intensity in summer. For both seasons, HREPA displayed resolution and sharpness. The overall good performance of HREPA and the lack of ensemble precipitation analysis (PA) at such spatiotemporal resolution in the literature motivate further investigations on transitional seasons and more advanced perturbation approaches. [ABSTRACT FROM AUTHOR]

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 one 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 8:30-11:30 local time) before satellite overpasses and for the mean emissions during other hours of the day (over the aggregation between 0:00-8:30 and 11:30-0:00) for each clump and for the 366 days 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 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 per year can potentially have their emissions between 8:30 and 11:30 constrained with a posterior uncertainty smaller than 20 % for more than 10 times within one year (ignoring the potential to cross or extrapolate information between 8: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 8: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 of systematic assessment of the sources of uncertainty in the spatiotemporally-resolved emission inventories used as prior estimates in the inversions. We highlight the difficulty to constrain global and national fossil fuel CO2 emissions with satellite XCO2 measurements only, and calls for integrated inversion systems that exploit multiple types of measurements. [ABSTRACT FROM AUTHOR]

This study presents an evaluation of the performance of the dynamically dimensioned search (DDS) algorithm when calibrating the hydrological component of the Visualizing Ecosystems for Land Management Assessments (VELMA) ecohydrological model. Two calibration strategies were tested for the initial parameter values: (1) a 'high-cost strategy', where 100 sets of initial parameter values were randomly chosen within the overall parameter space, and (2) a 'low-cost strategy', where a unique set of initial parameter values was derived from the available field data. Both strategies were tested for six different values of the maximum number of model evaluations ranging between 100 and 10,000. Results revealed that DDS is able to converge rapidly to a good parameter calibration solution of the VELMA hydrological component regardless of the parameter initialization strategy used. The accuracy and convergence efficiency of the DDS algorithm were, however, slightly better for the low-cost strategy. This study suggests that initializing the parameter values of complex physically based models using information on the watershed characteristics can increase the efficiency of the automatic calibration procedures. [ABSTRACT FROM AUTHOR]

A large fraction of fossil fuel CO2 emissions occur within "hotspots", such as cities and power plants, which cover a very small fraction of the land surface. Although some of these emission hotspots are monitored closely, there is no detailed emission inventory for most of them. Space-borne imagery of atmospheric CO2 has the potential to provide independent estimates of CO2 emissions from hotspots. But first, what is a hotspot 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 < 1 ppm for individual pixels. This resolution and precision is insufficient to provide a cartography of emissions for each individual pixel. Rather, the integrated emission of the diffuse emitting area and the intense point sources are sought. In this study, we address the question of the global characterization of area and point fossil fuel CO2 emitting sources (those hotspots are called emission clumps hereafter) that may cause coherent XCO2 plumes in space-borne CO2 images. 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 5,088 area clumps and 6,226 point-source clumps (power plants), are identified. These clumps contribute 72 % of the global fossil fuel CO2 emissions according to the ODIAS 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. [ABSTRACT FROM AUTHOR]

Background: Satellite imagery will offer unparalleled global spatial coverage at high-resolution for long term cost-effective monitoring of CO2 concentration plumes generated by emission hotspots. CO2 emissions can then be estimated from the magnitude of these plumes. In this paper, we assimilate pseudo-observations in a global atmospheric inversion system to assess the performance of a constellation of one to four sun-synchronous low-Earth orbit (LEO) imagers to monitor anthropogenic CO2 emissions. The constellation of imagers follows the specifications from the European Spatial Agency (ESA) for the Copernicus Anthropogenic Carbon Dioxide Monitoring (CO2M) concept for a future operational mission dedicated to the monitoring of anthropogenic CO2 emissions. This study assesses the uncertainties in the inversion estimates of emissions ("posterior uncertainties"). Results: The posterior uncertainties of emissions for individual cities and power plants are estimated for the 3 h before satellite overpasses, and extrapolated at annual scale assuming temporal auto-correlations in the uncertainties in the emission products that are used as a prior knowledge on the emissions by the Bayesian framework of the inversion. The results indicate that (i) the number of satellites has a proportional impact on the number of 3 h time windows for which emissions are constrained to better than 20%, but it has a small impact on the posterior uncertainties in annual emissions; (ii) having one satellite with wide swath would provide full images of the XCO2 plumes, and is more beneficial than having two satellites with half the width of reference swath; and (iii) an increase in the precision of XCO2 retrievals from 0.7 ppm to 0.35 ppm has a marginal impact on the emission monitoring performance. Conclusions: For all constellation configurations, only the cities and power plants with an annual emission higher than 0.5 MtC per year can have at least one 8:30–11:30 time window during one year when the emissions can be constrained to better than 20%. The potential of satellite imagers to constrain annual emissions not only depend on the design of the imagers, but also strongly depend on the temporal error structure in the prior uncertainties, which is needed to be objectively assessed in the bottom-up emission maps. [ABSTRACT FROM AUTHOR]