scholarly journals High-Resolution (2.5 km) Ensemble Precipitation Analysis across Canada

2020 ◽  
Vol 21 (9) ◽  
pp. 2023-2039
Author(s):  
Dikra Khedhaouiria ◽  
Stéphane Bélair ◽  
Vincent Fortin ◽  
Guy Roy ◽  
Franck Lespinas

AbstractConsistent 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.

2016 ◽  
Author(s):  
G. J. Schürmann ◽  
T. Kaminski ◽  
C. Köstler ◽  
N. Carvalhais ◽  
M. Voßbeck ◽  
...  

Abstract. We describe the Max Planck Institute Carbon Cycle Data Assimilation System (MPI-CCDAS) built around the tangent-linear version of the land surface scheme of the MPI-Earth System Model v1 (JSBACH). The simulated terrestrial biosphere processes (phenology and carbon balance) were constrained by observations of the fraction of photosynthetically active radiation (TIP-FAPAR product) and by observations of atmospheric CO2 at a global set of monitoring stations for the years 2005–2009. The system successfully, and computationally efficiently, improved average foliar area and northern extra-tropical seasonality of foliar area when constrained by TIP-FAPAR. Global net and gross carbon fluxes were improved when constrained by atmospheric CO2, although the system tended to underestimate tropical productivity. Assimilating both data streams jointly allowed the MPI-CCDAS to match both observations (TIP-FAPAR and atmospheric CO2) equally well as the single data stream assimilation cases, therefore overall increasing the appropriateness of the resultant parameter values and biosphere dynamics. Our study thus highlights the role of the TIP-FAPAR product in stabilising the underdetermined atmospheric inversion problem and demonstrates the value of multiple-data stream assimilation for the simulation of terrestrial biosphere dynamics. The constraint on regional gross and net CO2 flux patterns is limited through the parametrisation of the biosphere model. We expect improvement on that aspect through a refined initialisation strategy and inclusion of further biosphere observations as constraints.


2017 ◽  
Author(s):  
Wei He ◽  
Ivar R. van der Velde ◽  
Arlyn E. Andrews ◽  
Colm Sweeney ◽  
John Miller ◽  
...  

Abstract. We have implemented a regional carbon dioxide data assimilation system based on the CarbonTracker Data Assimilation Shell (CTDAS) and a high-resolution Lagrangian transport model, the Stochastic Time-Inverted Lagrangian Transport model driven by the Weather Forecast and Research meteorological fields (WRF-STILT). With this system, named as CTDAS‑Lagrange, we simultaneously optimize terrestrial biosphere fluxes and four parameters that adjust the lateral boundary conditions (BCs) against CO2 observations from the NOAA ESRL North America tall tower and aircraft Programmable Flask Packages (PFPs) sampling program. Least-squares optimization is performed with a time-stepping ensemble Kalman smoother, over a time window of 10 days and assimilating sequentially a time series of observations. Because the WRF-STILT footprints are pre-computed, it is computationally efficient to run the CTDAS-Lagrange system. To estimate the uncertainties of the optimized fluxes from the system, we performed sensitivity tests with various a priori biosphere fluxes (SiBCASA, SiB3, CT2013B) and BCs (optimized mole fraction fields from CT2013B and CTE2014, and an empirical data set derived from aircraft observations), as well as with a variety of choices on the ways that fluxes are adjusted (additive or multiplicative), covariance length scales, biosphere flux covariances, BC parameter uncertainties, and model-data mismatches. In pseudo-data experiments, we show that in our implementation the additive flux adjustment method is more flexible in optimizing NEE than the multiplicative flux adjustment method, and that the CTDAS-Lagrange system has the ability to correct for the potential biases in the lateral boundary conditions and to resolve large biases in the prior biosphere fluxes. Using real observations, we have derived a range of estimates for the optimized carbon fluxes from a series of sensitivity tests, which places the North American carbon sink for the year 2010 in a range from −0.92 to −1.26 PgC/yr. This is comparable to the TM5-based estimates of CarbonTracker (version CT2016, −0.91 ± 1.10 PgC/yr) and CarbonTracker Europe (version CTE2016, −0.91 ± 0.31 PgC/yr). We conclude that CTDAS-Lagrange can offer a versatile and computationally attractive alternative to these global systems for regional estimates of carbon fluxes, which can take advantage of high-resolution Lagrangian footprints that are increasingly easy to obtain.


2020 ◽  
Author(s):  
Bertrand Cluzet ◽  
Matthieu Lafaysse ◽  
Marie Dumont ◽  
Emmanuel Cosme ◽  
Clément Albergel

<p>In mountainous areas, detailed snowpack models are essential to capture the high spatio-temporal variability of the snowpack. This task is highly challenging, and models suffer from large simulation errors. In these regions, in-situ observations are scarce, while remote sensing observations are generally patchy owing to complex physiographic features (steep slopes, forests, shadows,...) and weather conditions (clouds). This point is stressing the need for a spatially coherent data assimilation system able to propagate the informations into unobserved locations.</p><p>In this study, we present CRAMPON (CRocus with AssiMilation of snowPack ObservatioNs), an ensemble data assimilation system ingesting snowpack observations in a spatialized context. CRAMPON quantifies snowpack modelling uncertainties with an ensemble and reduces them using a Particle Filter. Stochastic perturbations of meteorological forcings and the multi-physical version of Crocus snowpack model (ESCROC) are used to build the ensemble. Two variants of the Sequential Importance Resampling Particle Filter (PF) were implemented to tackle the common PF degeneracy issue that arises when assimilating a large number of observations. In a first approach (so-called global approach), the observations information is spread across topographic conditions by looking for a global analysis. Degeneracy is mitigated by inflating the observation error covariance matrix, with the side effect of reducing the impact of the assimilation. In a second approach (klocal), we propagate the information and mitigate degeneracy by a localisation of the PF based on background correlation patterns between topographic conditions.</p><p>Here, we investigate the ability of CRAMPON to globally benefit from partial observations in a conceptual semi-distributed domain which accounts for the main features of topographic-induced snowpack variability. We compare simulations without assimilation with experiments assimilating synthetic observations of the Height of Snow and VIS/NIR reflectance. This setup demonstrates the ability of CRAMPON to spread the information of various snow observations into unobserved locations.</p>


2015 ◽  
Vol 47 (5) ◽  
pp. 051401
Author(s):  
Yoichi Ishikawa ◽  
Teiji In ◽  
Satoshi Nakada ◽  
Kei Nishina ◽  
Hiromichi Igarashi ◽  
...  

2009 ◽  
Vol 78 (2) ◽  
pp. 237-248 ◽  
Author(s):  
Yoichi Ishikawa ◽  
Toshiyuki Awaji ◽  
Takahiro Toyoda ◽  
Teiji In ◽  
Kei Nishina ◽  
...  

2016 ◽  
Vol 17 (11) ◽  
pp. 2743-2761 ◽  
Author(s):  
Patrick D. Broxton ◽  
Xubin Zeng ◽  
Nicholas Dawson

Abstract There is a large uncertainty of snow water equivalent (SWE) in reanalyses and the Global Land Data Assimilation System (GLDAS), but the primary reason for this uncertainty remains unclear. Here several reanalysis products and GLDAS with different land models are evaluated and the primary reason for their deficiencies are identified using two high-resolution SWE datasets, including the Snow Data Assimilation System product and a new dataset for SWE and snowfall for the conterminous United States (CONUS) that is based on PRISM precipitation and temperature data and constrained with thousands of point snow observations of snowfall and snow thickness. The reanalyses and GLDAS products substantially underestimate SWE in the CONUS compared to the high-resolution SWE data. This occurs irrespective of biases in atmospheric forcing information or differences in model resolution. Furthermore, reanalysis and GLDAS products that predict more snow ablation at near-freezing temperatures have larger underestimates of SWE. Since many of the products do not assimilate information about SWE and snow thickness, this indicates a problem with the implementation of land models and pinpoints the need to improve the treatment of snow ablation in these systems, especially at near-freezing temperatures.


Atmosphere ◽  
2020 ◽  
Vol 11 (4) ◽  
pp. 365
Author(s):  
Yuanbing Wang ◽  
Yaodeng Chen ◽  
Jinzhong Min

An efficient regional hybrid ensemble-variational (EnVar) data assimilation method using the global-ensemble-model-augmented error covariance is proposed and preliminarily tested in this study. This method uses the global ensemble error covariance as the complementary low-resolution regional ensemble error covariance. The high-resolution dynamic ensemble mean is used as the first guess in hybrid EnVar and then re-centered to the updated high-resolution dynamic ensemble perturbations after minimization analysis. In this study, the proposed method is implemented into the Weather Research and Forecasting Model’s (WRF) data assimilation system coupled with the ensemble transform Kalman filter (ETKF) and preliminarily tested for numerical weather prediction during the Mei-Yu season over eastern China. It is found that the experiment containing fewer regional dynamic ensemble members but augmented with global ensemble error covariance obtains similar results to the experiment containing many more regional dynamic ensemble members. However, the former experiment only takes up one third of the latter experiment’s computational cost. The method proposed in this study also outperforms the 3DVar, hybrid EnVar using the pure global ensemble error covariance, as well as the hybrid EnVar using regional ETKF ensemble with a smaller size. The method proposed in this paper effectively combines the contributions of the ensemble error covariance from both the global and the regional models to produce better initial conditions for the regional WRF data assimilation system.


Sign in / Sign up

Export Citation Format

Share Document