Parallel Data Assimilation Framework Research Articles

WRF-PDAF v1.0: implementation and application of an online localized ensemble data assimilation framework

Abstract. Data assimilation is a common technique employed to estimate the state and its associated uncertainties in numerical models. Ensemble-based methods are a prevalent choice, although they can be computationally expensive due to the required ensemble integrations. In this study, we enhance the capabilities of the Weather Research and Forecasting–Advanced Research WRF (WRF-ARW) model by coupling it with the Parallel Data Assimilation Framework (PDAF) in a fully online mode. Through minimal modifications to the WRF-ARW model code, we have developed an efficient data assimilation system. This system leverages parallelization and in-memory data transfers between the model and data assimilation processes, greatly reducing the need for file I/O and model restarts during assimilation. We detail the necessary program modifications in this study. One advantage of the resulting assimilation system is a clear separation of concerns between data assimilation method development and model application resulting from PDAF's model-agnostic structure. To evaluate the assimilation system, we conduct a twin experiment simulating an idealized tropical cyclone. Cycled data assimilation experiments focus on the impact of temperature profiles. The assimilation not only significantly enhances temperature field accuracy but also improves the initial U and V fields. The assimilation process introduces only minimal overhead in runtime when compared to the model without data assimilation and exhibits excellent parallel performance. Consequently, the online WRF-PDAF system emerges as an efficient framework for implementing high-resolution mesoscale forecasting and reanalysis.

HGS-PDAF (version 1.0): a modular data assimilation framework for an integrated surface and subsurface hydrological model

Abstract. This article describes a modular ensemble-based data assimilation (DA) system which is developed for an integrated surface–subsurface hydrological model. The software environment for DA is the Parallel Data Assimilation Framework (PDAF), which provides various assimilation algorithms like the ensemble Kalman filters, non-linear filters, 3D-Var and combinations among them. The integrated surface–subsurface hydrological model is HydroGeoSphere (HGS), a physically based modelling software for the simulation of surface and variably saturated subsurface flow, as well as heat and mass transport. The coupling and capabilities of the modular DA system are described and demonstrated using an idealised model of a geologically heterogeneous alluvial river–aquifer system with drinking water production via riverbank filtration. To demonstrate its modularity and adaptability, both single and multivariate assimilations of hydraulic head and soil moisture observations are demonstrated in combination with individual and joint updating of multiple simulated states (i.e. hydraulic heads and water saturation) and model parameters (i.e. hydraulic conductivity). With the integrated model and this modular DA framework, we have essentially developed the hydrologically and DA-wise robust toolbox for developing the basic model for operational management of coupled surface water–groundwater resources.

Paleoclimate data assimilation with CLIMBER-X: An ensemble Kalman filter for the last deglaciation.

Using the climate model CLIMBER-X, we present an efficient method for assimilating the temporal evolution of surface temperatures for the last deglaciation covering the period 22000 to 6500 years before the present. The data assimilation methodology combines the data and the underlying dynamical principles governing the climate system to provide a state estimate of the system, which is better than that which could be obtained using just the data or the model alone. In applying an ensemble Kalman filter approach, we make use of the advances in the parallel data assimilation framework (PDAF), which provides parallel data assimilation functionality with a relatively small increase in computation time. We find that the data assimilation solution depends strongly on the background evolution of the decaying ice sheets rather than the assimilated temperatures. Two different ice sheet reconstructions result in a different deglacial meltwater history, affecting the large-scale ocean circulation and, consequently, the surface temperature. We find that the influence of data assimilation is more pronounced on regional scales than on the global mean. In particular, data assimilation has a stronger effect during millennial warming and cooling phases, such as the Bølling-Allerød and Younger Dryas, especially at high latitudes with heterogeneous temperature patterns. Our approach is a step toward a comprehensive paleo-reanalysis on multi-millennial time scales, including incorporating available paleoclimate data and accounting for their uncertainties in representing regional climates.

Evapotranspiration prediction for European forest sites does not improve with assimilation of in situ soil water content data

Abstract. Land surface models (LSMs) are an important tool for advancing our knowledge of the Earth system. LSMs are constantly improved to represent the various terrestrial processes in more detail. High-quality data, freely available from various observation networks, are being used to improve the prediction of terrestrial states and fluxes of water and energy. To optimize LSMs with observations, data assimilation methods and tools have been developed in the past decades. We apply the coupled Community Land Model version 5 (CLM5) and Parallel Data Assimilation Framework (PDAF) system (CLM5-PDAF) for 13 forest field sites throughout Europe covering different climate zones. The goal of this study is to assimilate in situ soil moisture measurements into CLM5 to improve the modeled evapotranspiration fluxes. The modeled fluxes will be evaluated using the predicted evapotranspiration fluxes with eddy covariance (EC) systems. Most of the sites use point-scale measurements from sensors placed in the ground; however, for three of the forest sites we use soil water content data from cosmic-ray neutron sensors, which have a measurement scale closer to the typical land surface model grid scale and EC footprint. Our results show that while data assimilation reduced the root-mean-square error for soil water content on average by 56 % to 64 %, the root-mean-square error for the evapotranspiration estimation is increased by 4 %. This finding indicates that only improving the soil water content (SWC) estimation of state-of-the-art LSMs such as CLM5 is not sufficient to improve evapotranspiration estimates for forest sites. To improve evapotranspiration estimates, it is also necessary to consider the representation of leaf area index (LAI) in magnitude and timing, as well as uncertainties in water uptake by roots and vegetation parameters.

The Impact of Profiles Data Assimilation on an Ideal Tropical Cyclone Case

Profile measurements play a crucial role in operational weather forecasting across diverse scales and latitudes. However, assimilating tropospheric wind and temperature profiles remains a challenging endeavor. This study assesses the influence of profile measurements on numerical weather prediction (NWP) using the weather research and forecasting (WRF) model coupled to the parallel data assimilation framework (PDAF) system. Utilizing the local error-subspace transform Kalman filter (LESTKF), observational temperature and wind profiles generated by WRF are assimilated into an idealized tropical cyclone. The coupled WRF-PDAF system is adopted to carry out the twin experiments, which employ varying profile densities and localization distances. The results reveal that high-resolution observations yield significant forecast improvements compared to coarser-resolution data. A cost-effective balance between observation density and benefit is further explored through the idealized tropical cyclone case. According to diminishing marginal utility and increasing marginal costs, the optimal observation densities for U and V are found around 26–27%. This may be useful information to the meteorological agencies and researchers.

Basin Scale Soil Moisture Estimation with Grid SWAT and LESTKF Based on WSN.

This research utilized in situ soil moisture observations in a coupled grid Soil and Water Assessment Tool (SWAT) and Parallel Data Assimilation Framework (PDAF) data assimilation system, resulting in significant enhancements in soil moisture estimation. By incorporating Wireless Sensor Network (WSN) data (WATERNET), the method captured and integrated local soil moisture characteristics, thereby improving regional model state estimations. The use of varying observation search radii with the Local Error-subspace Transform Kalman Filter (LESTKF) resulted in improved spatial and temporal assimilation performance, while also considering the impact of observation data uncertainties. The best performance (improvement of 0.006 m3/m3) of LESTKF was achieved with a 20 km observation search radii and 0.01 m3/m3 observation standard error. This study assimilated wireless sensor network data into a distributed model, presenting a departure from traditional methods. The high accuracy and resolution capabilities of WATERNET's regional soil moisture observations were crucial, and its provision of multi-layered soil temperature and moisture observations presented new opportunities for integration into the data assimilation framework, further enhancing hydrological state estimations. This study's implications are broad and relevant to regional-scale water resource research and management, particularly for freshwater resource scheduling at small basin scales.

The effects of assimilating a sub-grid-scale sea ice thickness distribution in a new Arctic sea ice data assimilation system

Abstract. In the past decade groundbreaking new satellite observations of the Arctic sea ice cover have been made, allowing researchers to understand the state of the Arctic sea ice system in greater detail than before. The derived estimates of sea ice thickness are useful but limited in time and space. In this study the first results of a new sea ice data assimilation system are presented. Observations assimilated (in various combinations) are monthly mean sea ice thickness and monthly mean sea ice thickness distribution from CryoSat-2 and NASA daily Bootstrap sea ice concentration. This system couples the Centre for Polar Observation and Modelling's (CPOM) version of the Los Alamos Sea Ice Model (CICE) to the localised ensemble transform Kalman filter (LETKF) from the Parallel Data Assimilation Framework (PDAF) library. The impact of assimilating a sub-grid-scale sea ice thickness distribution is of particular novelty. The sub-grid-scale sea ice thickness distribution is a fundamental component of sea ice models, playing a vital role in the dynamical and thermodynamical processes, yet very little is known of its true state in the Arctic. This study finds that assimilating CryoSat-2 products for the mean thickness and the sub-grid-scale thickness distribution can have significant consequences for the modelled distribution of the ice thickness across the Arctic and particularly in regions of thick multi-year ice. The assimilation of sea ice concentration, mean sea ice thickness and sub-grid-scale sea ice thickness distribution together performed best when compared to a subset of CryoSat-2 observations held back for validation. Regional model biases are reduced: the thickness of the thickest ice in the Canadian Arctic Archipelago (CAA) is decreased, but the thickness of the ice in the central Arctic is increased. When comparing the assimilation of mean thickness with the assimilation of sub-grid-scale thickness distribution, it is found that the latter leads to a significant change in the volume of ice in each category. Estimates of the thickest ice improve significantly with the assimilation of sub-grid-scale thickness distribution alongside mean thickness.

Improving the estimation of thermospheric neutral density via two-step assimilation of in situ neutral density into a numerical model

AbstractNeutral thermospheric density is an essential quantity required for precise orbit determination of satellites, collision avoidance of satellites, re-entry prediction of satellites or space debris, and satellite lifetime assessments. Empirical models of the thermosphere fail to provide sufficient estimates of neutral thermospheric density along the orbits of satellites by reason of approximations, assumptions and a limited temporal resolution. At high solar activity these estimates can be off by 70% when comparing to observations at 12-hourly averages. In recent decades, neutral density is regularly observed with satellite accelerometers on board of low Earth orbiting satellites like CHAMP, GOCE, GRACE, GRACE-FO, or Swarm. When assimilating such along-track information into global models of thermosphere–ionosphere dynamics, it has been often observed that only a very local sub-domain of the model grid around the satellite’s position is updated. To extend the impact to the entire model domain we suggest a new two-step approach: we use accelerometer-derived neutral densities from the CHAMP mission in a first step to calibrate an empirical thermosphere density model (NRLMSIS 2.0). In a second step, we assimilate—for the first time—densities predicted for a regular three-dimensional grid into the TIE-GCM (Thermosphere Ionosphere Electrodynamics General Circulation Model). Data assimilation is performed using the Local Error-Subspace Transform Kalman Filter provided by the Parallel Data Assimilation Framework (PDAF). We test the new approach using a 2-week-long period containing the 5 April 2010 Geomagnetic storm. Accelerometer-derived neutral densities from the GRACE mission are used for additional evaluation. We demonstrate that the two-step approach globally improves the simulation of thermospheric density. We could significantly improve the density prediction for CHAMP and GRACE. In fact, the offset between the accelerometer-derived densities and the model prediction is reduced by 45% for CHAMP and 20% for GRACE when applying the two-step approach. The implication is that our approach allows one to much better ’transplant’ the precise CHAMP thermospheric density measurements to satellites flying at a similar altitude. Graphical Abstract

Uncertainties in ocean biogeochemical simulations: Application of ensemble data assimilation to a one-dimensional model

Marine biogeochemical (BGC) models are highly uncertain in their parameterization. The value of the BGC parameters are poorly known and lead to large uncertainties in the model outputs. This study focuses on the uncertainty quantification of model fields and parameters within a one-dimensional (1-D) ocean BGC model applying ensemble data assimilation. We applied an ensemble Kalman filter provided by the Parallel Data Assimilation Framework (PDAF) into a 1-D vertical configuration of the BGC model Regulated Ecosystem Model 2 (REcoM2) at two BGC time-series stations: the Bermuda Atlantic Time-series Study (BATS) and the Dynamique des Flux Atmosphériques en Méditerranée (DYFAMED). We assimilated 5-day satellite chlorophyll-a (chl-a) concentration and monthly in situ net primary production (NPP) data for 3 years to jointly estimate 10 preselected key BGC parameters and the model state. The estimated set of parameters resulted in improvements in the model prediction up to 66% for the surface chl-a and 56% for NPP. Results show that assimilating satellite chl-a concentration data alone degraded the prediction of NPP. Simultaneous assimilation of the satellite chl-a data and in situ NPP data improved both surface chl-a and NPP simulations. We found that correlations between parameters preclude estimating parameters independently. Co-dependencies between parameters also indicate that there is not a unique set of optimal parameters. Incorporation of proper uncertainty estimation in BGC predictions, therefore, requires ensemble simulations with varying parameter values.

An aerosol vertical data assimilation system (NAQPMS-PDAF v1.0): development and application

Abstract. Aerosol vertical stratification is important for global climate and planetary boundary layer (PBL) stability, and no single method can obtain spatiotemporally continuous vertical profiles. This paper develops an online data assimilation (DA) framework for the Eulerian atmospheric chemistry-transport model (CTM) Nested Air Quality Prediction Model System (NAQPMS) with the Parallel Data Assimilation Framework (PDAF) as the NAQPMS-PDAF for the first time. Online coupling occurs based on a memory-based way with two-level parallelization, and the arrangement of state vectors during the filter is specifically designed. Scaling tests demonstrate that the NAQPMS-PDAF can make efficient use of parallel computational resources for up to 25 000 processors with a weak scaling efficiency of up to 0.7. The 1-month long aerosol extinction coefficient profiles measured by the ground-based lidar and the concurrent hourly surface PM2.5 are solely and simultaneously assimilated to investigate the performance and application of the DA system. The hourly analysis and subsequent 1 h simulation are validated through lidar and surface PM2.5 measurements assimilated and not assimilated. The results show that lidar DA can significantly improve the underestimation of aerosol loading, especially at a height of approximately 400 m in the free-running (FR) experiment, with the mean bias (BIAS) changing from −0.20 (−0.14) km−1 to −0.02 (−0.01) km−1 and correlation coefficients increasing from 0.33 (0.28) to 0.91 (0.53) averaged over sites with measurements assimilated (not assimilated). Compared with the FR experiment, simultaneously assimilating PM2.5 and lidar can have a more consistent pattern of aerosol vertical profiles with a combination of surface PM2.5 and lidar, independent extinction coefficients from the Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP), and aerosol optical depth (AOD) from the Aerosol Robotic Network (AERONET). Lidar DA has a larger temporal impact than that in PM2.5 DA but has deficiencies in subsequent quantification on the surface PM2.5. The proposed NAQPMS-PDAF has great potential for further research on the impact of aerosol vertical distribution.

Data assimilation of volcanic aerosol observations using FALL3D+PDAF

Abstract. Modelling atmospheric dispersal of volcanic ash and aerosols is becoming increasingly valuable for assessing the potential impacts of explosive volcanic eruptions on buildings, air quality, and aviation. Management of volcanic risk and reduction of aviation impacts can strongly benefit from quantitative forecasting of volcanic ash. However, an accurate prediction of volcanic aerosol concentrations using numerical modelling relies on proper estimations of multiple model parameters which are prone to errors. Uncertainties in key parameters such as eruption column height and physical properties of particles or meteorological fields represent a major source of error affecting the forecast quality. The availability of near-real-time geostationary satellite observations with high spatial and temporal resolutions provides the opportunity to improve forecasts in an operational context by incorporating observations into numerical models. Specifically, ensemble-based filters aim at converting a prior ensemble of system states into an analysis ensemble by assimilating a set of noisy observations. Previous studies dealing with volcanic ash transport have demonstrated that a significant improvement of forecast skill can be achieved by this approach. In this work, we present a new implementation of an ensemble-based data assimilation (DA) method coupling the FALL3D dispersal model and the Parallel Data Assimilation Framework (PDAF). The FALL3D+PDAF system runs in parallel, supports online-coupled DA, and can be efficiently integrated into operational workflows by exploiting high-performance computing (HPC) resources. Two numerical experiments are considered: (i) a twin experiment using an incomplete dataset of synthetic observations of volcanic ash and (ii) an experiment based on the 2019 Raikoke eruption using real observations of SO2 mass loading. An ensemble-based Kalman filtering technique based on the local ensemble transform Kalman filter (LETKF) is used to assimilate satellite-retrieved data of column mass loading. We show that this procedure may lead to nonphysical solutions and, consequently, conclude that LETKF is not the best approach for the assimilation of volcanic aerosols. However, we find that a truncated state constructed from the LETKF solution approaches the real solution after a few assimilation cycles, yielding a dramatic improvement of forecast quality when compared to simulations without assimilation.

Coupling the Community Land Model version 5.0 to the parallel data assimilation framework PDAF: description and applications

Abstract. Land surface models are important for improving our understanding of the Earth system. They are continuously improving and becoming better in representing the different land surface processes, e.g., the Community Land Model version 5 (CLM5). Similarly, observational networks and remote sensing operations are increasingly providing more data, e.g., from new satellite products and new in situ measurement sites, with increasingly higher quality for a range of important variables of the Earth system. For the optimal combination of land surface models and observation data, data assimilation techniques have been developed in recent decades that incorporate observations to update modeled states and parameters. The Parallel Data Assimilation Framework (PDAF) is a software environment that enables ensemble data assimilation and simplifies the implementation of data assimilation systems in numerical models. In this study, we present the development of the new interface between PDAF and CLM5. This newly implemented coupling integrates the PDAF functionality into CLM5 by modifying the CLM5 ensemble mode to keep changes to the pre-existing parallel communication infrastructure to a minimum. Soil water content observations from an extensive in situ measurement network in the Wüstebach catchment in Germany are used to illustrate the application of the coupled CLM5-PDAF system. The results show overall reductions in root mean square error of soil water content from 7 % up to 35 % compared to simulations without data assimilation. We expect the coupled CLM5-PDAF system to provide a basis for improved regional to global land surface modeling by enabling the assimilation of globally available observational data.

Strongly Coupled Data Assimilation of Ocean Observations Into an Ocean‐Atmosphere Model

AbstractWe compare strongly coupled data assimilation (SCDA) and weakly coupled data assimilation (WCDA) by analyzing the assimilation effect on the estimation of the ocean and the atmosphere variables. The AWI climate model (AWI‐CM‐1.1) is coupled with the parallel data assimilation framework (PDAF). Only satellite sea surface temperature data are assimilated. For WCDA, only the ocean variables are directly updated by the assimilation. For SCDA, both the ocean and the atmosphere variables are directly updated by the assimilation. Both WCDA and SCDA improve ocean state and yield similar errors. In the atmosphere, WCDA gives slightly smaller errors for the near‐surface temperature and wind velocity than SCDA. In the free atmosphere, SCDA yields smaller errors for the temperature, wind velocity, and specific humidity than WCDA in the Arctic region, while in the tropical region, the errors are generally larger.

Efficient ensemble data assimilation for coupled models with the Parallel Data Assimilation Framework: example of AWI-CM (AWI-CM-PDAF 1.0)

Abstract. Data assimilation integrates information from observational measurements with numerical models. When used with coupled models of Earth system compartments, e.g., the atmosphere and the ocean, consistent joint states can be estimated. A common approach for data assimilation is ensemble-based methods which utilize an ensemble of state realizations to estimate the state and its uncertainty. These methods are far more costly to compute than a single coupled model because of the required integration of the ensemble. However, with uncoupled models, the ensemble methods also have been shown to exhibit a particularly good scaling behavior. This study discusses an approach to augment a coupled model with data assimilation functionality provided by the Parallel Data Assimilation Framework (PDAF). Using only minimal changes in the codes of the different compartment models, a particularly efficient data assimilation system is generated that utilizes parallelization and in-memory data transfers between the models and the data assimilation functions and hence avoids most of the file reading and writing, as well as model restarts during the data assimilation process. This study explains the required modifications to the programs with the example of the coupled atmosphere–sea-ice–ocean model AWI-CM (AWI Climate Model). Using the case of the assimilation of oceanic observations shows that the data assimilation leads only to small overheads in computing time of about 15 % compared to the model without data assimilation and a very good parallel scalability. The model-agnostic structure of the assimilation software ensures a separation of concerns in which the development of data assimilation methods can be separated from the model application.

Improving the ocean and atmosphere in a coupled ocean–atmosphere model by assimilating satellite sea‐surface temperature and subsurface profile data

AbstractAn ensemble‐based data assimilation framework for a coupled ocean–atmosphere model is applied to investigate the influence of assimilating different types of ocean observations on the ocean and atmosphere simulation. The data assimilation is performed with the parallel data assimilation framework (PDAF) for the climate model AWI‐CM. Observations of the ocean, namely satellite sea‐surface temperature (SST) and temperature and salinity profiles, are assimilated into the ocean component. The atmospheric state is only influenced by the model dynamics. Different assimilation scenarios were carried out with different combinations of observations to investigate to what extent the assimilation into the coupled model leads to a better estimation of the state of the ocean as well as the atmosphere. The influence of the data assimilation is assessed by comparing the ocean prediction with dependent and independent ocean observations. For the atmosphere, the assimilation result is compared with the ERA‐Interim atmospheric reanalysis data. The ocean temperature and salinity are improved by all the assimilation scenarios in the coupled system. The assimilation leads to a response of the atmosphere throughout the troposphere and impacts the global atmospheric circulation. Globally the temperature and wind speed are improved in the atmosphere on average.

An offline framework for high-dimensional ensemble Kalman filters to reduce the time to solution

Abstract. The high computational resources and the time-consuming IO (input/output) are major issues in offline ensemble-based high-dimensional data assimilation systems. Bearing these in mind, this study proposes a sophisticated dynamically running job scheme as well as an innovative parallel IO algorithm to reduce the time to solution of an offline framework for high-dimensional ensemble Kalman filters. The dynamically running job scheme runs as many tasks as possible within a single job to reduce the queuing time and minimize the overhead of starting and/or ending a job. The parallel IO algorithm reads or writes non-overlapping segments of multiple files with an identical structure to reduce the IO times by minimizing the IO competitions and maximizing the overlapping of the MPI (Message Passing Interface) communications with the IO operations. Results based on sensitive experiments show that the proposed parallel IO algorithm can significantly reduce the IO times and have a very good scalability, too. Based on these two advanced techniques, the offline and online modes of ensemble Kalman filters are built based on PDAF (Parallel Data Assimilation Framework) to comprehensively assess their efficiencies. It can be seen from the comparisons between the offline and online modes that the IO time only accounts for a small fraction of the total time with the proposed parallel IO algorithm. The queuing time might be less than the running time in a low-loaded supercomputer such as in an operational context, but the offline mode can be nearly as fast as, if not faster than, the online mode in terms of time to solution. However, the queuing time is dominant and several times larger than the running time in a high-loaded supercomputer. Thus, the offline mode is substantially faster than the online mode in terms of time to solution, especially for large-scale assimilation problems. From this point of view, results suggest that an offline ensemble Kalman filter with an efficient implementation and a high-performance parallel file system should be preferred over its online counterpart for intermittent data assimilation in many situations.

Seasonal Arctic Sea Ice Prediction Using a Newly Developed Fully Coupled Regional Model With the Assimilation of Satellite Sea Ice Observations

AbstractTo increase our capability to predict Arctic sea ice and climate, we have developed a coupled atmosphere‐sea ice‐ocean model configured for the pan‐Arctic with sufficient flexibility. The Los Alamos Sea Ice Model is coupled with the Weather Research and Forecasting Model and the Regional Ocean Modeling System in the Coupled Ocean‐Atmosphere‐Wave‐Sediment Transport modeling system. It is well known that dynamic models used to predict Arctic sea ice at short‐term periods strongly depend on model initial conditions. Parallel Data Assimilation Framework is implemented into the new modeling system to assimilate sea ice observations and generate skillful model initialization, which aid in the prediction procedures. The Special Sensor Microwave Imager/Sounder sea ice concentration, the CyroSat‐2, and Soil Moisture and Ocean Salinity sea ice thickness are assimilated with the localized error subspace transform ensemble Kalman filter. We conduct Arctic sea ice prediction for the melting seasons of 2017 and 2018. Predictions with improved initial sea ice conditions show reasonable sea ice evolution and small biases in the minimum sea ice extent, although the ice refreezing is delayed. Our prediction experiments suggest that the use of appropriate uncertainty for the observed sea ice thickness can lead to improved spatial distribution of the initial ice thickness and thus the predicted sea ice distribution. Our new modeling system initialized by the output of the National Centers for Environmental Prediction Climate Forecast System seasonal forecasts with data assimilation can significantly increase the sea ice prediction skills in sea ice extent for the entire Arctic as well as in the Northern Sea Route compared with the predictions by the National Centers for Environmental Prediction Climate Forecast System.

Improving soil moisture and runoff simulations at 3 km over Europe using land surface data assimilation

Abstract. Accurate and reliable hydrologic simulations are important for many applications such as water resources management, future water availability projections and predictions of extreme events. However, the accuracy of water balance estimates is limited by the lack of large-scale observations, model simulation uncertainties and biases related to errors in model structure and uncertain inputs (e.g., hydrologic parameters and atmospheric forcings). The availability of long-term and global remotely sensed soil moisture offers the opportunity to improve model estimates through data assimilation with complete spatiotemporal coverage. In this study, we assimilated the European Space Agency (ESA) Climate Change Initiative (CCI) derived soil moisture (SM) information to improve the estimation of continental-scale soil moisture and runoff. The assimilation experiment was conducted over a time period 2000–2006 with the Community Land Model, version 3.5 (CLM3.5), integrated with the Parallel Data Assimilation Framework (PDAF) at a spatial resolution of 0.0275∘ (∼3 km) over Europe. The model was forced with the high-resolution reanalysis COSMO-REA6 from the Hans Ertel Centre for Weather Research (HErZ). The performance of assimilation was assessed against open-loop model simulations and cross-validated with independent ESA CCI-derived soil moisture (CCI-SM) and gridded runoff observations. Our results showed improved estimates of soil moisture, particularly in the summer and autumn seasons when cross-validated with independent CCI-SM observations. The assimilation experiment results also showed overall improvements in runoff, although some regions were degraded, especially in central Europe. The results demonstrated the potential of assimilating satellite soil moisture observations to produce downscaled and improved high-resolution soil moisture and runoff simulations at the continental scale, which is useful for water resources assessment and monitoring.

Arctic‐Wide Sea Ice Thickness Estimates From Combining Satellite Remote Sensing Data and a Dynamic Ice‐Ocean Model with Data Assimilation During the CryoSat‐2 Period

AbstractExploiting the complementary character of CryoSat‐2 and Soil Moisture and Ocean Salinity satellite sea ice thickness products, daily Arctic sea ice thickness estimates from October 2010 to December 2016 are generated by an Arctic regional ice‐ocean model with satellite thickness assimilated. The assimilation is performed by a Local Error Subspace Transform Kalman filter coded in the Parallel Data Assimilation Framework. The new estimates can be generally thought of as combined model and satellite thickness (CMST). It combines the skill of satellite thickness assimilation in the freezing season with the model skill in the melting season, when neither CryoSat‐2 nor Soil Moisture and Ocean Salinity sea ice thickness is available. Comparisons with in situ observations from the Beaufort Gyre Exploration Project, Ice Mass Balance Buoys, and the NASA Operation IceBridge demonstrate that CMST reproduces most of the observed temporal and spatial variations. Results also show that CMST compares favorably to the Pan‐Arctic Ice‐Ocean Modeling and Assimilation System product and even appears to correct known thickness biases in the Pan‐Arctic Ice‐Ocean Modeling and Assimilation System. Due to imperfect parameterizations in the sea ice model and satellite thickness retrievals, CMST does not reproduce the heavily deformed and ridged sea ice along the northern coast of the Canadian Arctic Archipelago and Greenland. With the new Arctic sea ice thickness estimates sea ice volume changes in recent years can be further assessed.

Improving sea ice thickness estimates by assimilating CryoSat‐2 and SMOS sea ice thickness data simultaneously

The impact of assimilating weekly CryoSat‐2 sea ice thickness data together with daily SMOS sea ice thickness and daily SSMIS sea ice concentration data on the sea ice fields of a coupled sea ice–ocean model of the Arctic Ocean is investigated. The sea‐ice model is based on the Massachusetts Institute of Technology general circulation model (MITgcm) and the assimilation is performed by a localized Singular Evolutive Interpolated Kalman (LSEIK) filter coded in the Parallel Data Assimilation Framework (PDAF). A period of three months from 1 November 2011 to 30 January 2012 is selected to assess the skill of the assimilation system in the cold season. Compared to the unassimilated solution and a solution where only sea ice concentration is assimilated, the model–data misfits are substantially reduced in areas of both thick and thin ice. The sea ice thickness estimates agree significantly better with in situ observations in the central Arctic Ocean than the sea ice thickness obtained from assimilating SMOS data alone, while the sea ice concentration shows very small improvements. The sea ice fields obtained by the joint assimilation of SMOS and CryoSat‐2 data also have lower errors in thickness and concentration than those obtained from directly assimilating a statistically merged SMOS and CryoSat‐2 sea ice thickness product. These lower errors suggest that model dynamics play a significant role in data blending.