Ensemble Data Assimilation System Research Articles

Investigating the mechanisms of an intense coastal rainfall event during TAHOPE/PRECIP-IOP3 using a multiscale radar ensemble data assimilation system

Abstract The joint Taiwan-Area Heavy Rain Observation and Prediction Experiment (TAHOPE)/Prediction of Rainfall Extremes Campaign In the Pacific (PRECIP) field campaign between Taiwan and the United States took place from late May to mid-August in 2022. The field campaign aimed to understand the dynamics, thermodynamics, and predictability of heavy rainfall events in the Taiwan area. This study investigated the mechanisms of a heavy rainfall event that occurred on 6–7 June during the intensive observation period-3 (IOP3) of the field campaign. Heavy rainfall occurs on Taiwan’s western coast when a Meiyu front hovers in northern Taiwan. A multiscale radar ensemble data assimilation system based on the successive covariance localization (SCL) method is used to derive a high-resolution analysis for forecasts. Two numerical experiments are conducted with the use of convective-scale (RDA) or multiscale (MRDA) corrections in the assimilation of the radial velocity from operational radars at Chigu and Wufen, and the additional S-Pol radar deployed at Hsinchu during the field campaign. Compared with RDA, MRDA results in large-area wind corrections, which help reshape and relocate a low-level mesoscale vortex, a key element of this heavy rainfall event, offshore of western central Taiwan and enhances the front intensity offshore of northwestern Taiwan. Consequently, MRDA improves the 6-h heavy rainfall prediction over the coast of western Taiwan and better represents the elongated rainband in northern Taiwan during the 3- to 6-h forecast. Sensitivity experiments demonstrate the importance of assimilating winds from Chigu and S-Pol radar in establishing low-level mesoscale vortex and convergence zones.

Insights from very‐large‐ensemble data assimilation experiments with a high‐resolution general circulation model of the Red Sea

AbstractEnsemble Kalman Filters (EnKFs), which assimilate observations based on statistics derived from an ensemble of samples of ocean states, have become the norm for ocean data assimilation (DA) and forecasting. These schemes are commonly implemented with inflation and localization techniques to increase their ensemble spread and to filter out spurious long‐range correlations resulting from the limited‐size ensembles imposed by computational burden constraints. Such ad‐hoc methods were found to be not necessary in ensemble DA experiments with simplified ocean/atmospheric models and large ensembles. Here, we conduct a series of one‐year‐long ensemble experiments with a fully realistic EnKF‐DA system in the Red Sea using tens ‐to thousands of ensemble members. The system assimilates satellite and in‐situ observations and accounts for model uncertainties by integrating a 4‐km‐resolution ocean model with European Center for Medium Range Weather Forecast (ECMWF) atmospheric ensemble fields, perturbed internal physics and initial conditions for forecasting. OceanOur results indicate that accounting for model uncertainties is more beneficial than simply increasing the ensemble size, with the improvements due to large ensembles leveling off at about 250 members. Besides, and in contrast to what is commonly observed with simplified models, the investigated ensemble DA system still required localization even when implemented with thousands of members. These findings are explained by: (i) amplified spurious long‐range correlations produced by the low‐rank nature of the ECMWF atmospheric forcing ensemble; and (ii) non‐Gaussianity generated by the perturbed internal physical parameterization schemes. Large‐ensemble forcing fields and non‐Gaussian DA methods might be needed to get full benefits from large ensembles in ocean DA.

РАЗВИТИЕ МЕТОДОВ И ТЕХНОЛОГИЙ ОПЕРАТИВНОГО УСВОЕНИЯ ДАННЫХ МЕТЕОРОЛОГИЧЕСКИХ НАБЛЮДЕНИЙ В ГИДРОМЕТЦЕНТРЕ РОССИИ

A brief description of the system functioning at the Hydrometcenter of Russia for operational assimilation of meteorological observations and results of the work on its development are given. The main direction of its development is creating an original ensemble variational data assimilation system using neural network modeling techniques. Originality of the system consists in a new multiscale convolutional analysis. Results of testing the new data assimilation technique for a two-dimensional case are presented. In numerical experiments, the convolutional analysis has demonstrated higher accuracy than the traditional ensemble variational approach. A new multiscale technique is described, and results of its testing are given. The second direction of the development is increasing efficiency of satellite data assimilation. Results of retrieving sea ice concentration using data of AMSR2 and MTVZA-GYa microwave radiometers using machine learning methods are presented. A technique of assimilation of delayed observations is proposed and tested.

Potential impact of Aeroclipper observations targeting tropical cyclone in the Western Pacific

AbstractThe Aeroclipper is a new balloon device that can be attracted and captured by tropical cyclones (TC) and perform continuous in situ measurements at the air–sea interfaces. To estimate the potential effect of Aeroclipper observations on the analysis of TCs, virtual Aeroclipper observations targeting TC Haima (October 2016) were synthesized using an idealized surface pressure distribution and best track data and were assimilated using an ensemble data assimilation system. Results show that the assimilation of Aeroclipper measurements may provide a more accurate representation of the TC pressure, wind, and temperature in analyses. This also leads to improved precipitation around the Philippines. The ensemble spread shows that the Aeroclipper measurement assimilation has an impact on the analyses that extends into the tropics from the early stages of TC development. These impact signals propagate westward with easterly waves and eastward with large‐scale convective disturbances. Although the underlying mechanisms need to be further examined and tested using real Aeroclipper measurements, the present study shows that these balloons could provide valuable observations to improve the precision of analyses in presence of a TC. This is a first step toward a study of the impact of the Aeroclipper measurement on TC forecast.

Comparison of 4-dimensional variational and ensemble optimal interpolation data assimilation systems using a Regional Ocean Modeling System (v3.4) configuration of the eddy-dominated East Australian Current system

Abstract. Ocean models must be regularly updated through the assimilation of observations (data assimilation) in order to correctly represent the timing and locations of eddies. Since initial conditions play an important role in the quality of short-term ocean forecasts, an effective data assimilation scheme to produce accurate state estimates is key to improving prediction. Western boundary current regions, such as the East Australia Current system, are highly variable regions, making them particularly challenging to model and predict. This study assesses the performance of two ocean data assimilation systems in the East Australian Current system over a 2-year period. We compare the time-dependent 4-dimensional variational (4D-Var) data assimilation system with the more computationally efficient, time-independent ensemble optimal interpolation (EnOI) system, across a common modelling and observational framework. Both systems assimilate the same observations: satellite-derived sea surface height, sea surface temperature, vertical profiles of temperature and salinity (from Argo floats), and temperature profiles from expendable bathythermographs. We analyse both systems' performance against independent data that are withheld, allowing a thorough analysis of system performance. The 4D-Var system is 25 times more expensive but outperforms the EnOI system against both assimilated and independent observations at the surface and subsurface. For forecast horizons of 5 d, root-mean-squared forecast errors are 20 %–60 % higher for the EnOI system compared to the 4D-Var system. The 4D-Var system, which assimilates observations over 5 d windows, provides a smoother transition from the end of the forecast to the subsequent analysis field. The EnOI system displays elevated low-frequency (>1 d) surface-intensified variability in temperature and elevated kinetic energy at length scales less than 100 km at the beginning of the forecast windows. The 4D-Var system displays elevated energy in the near-inertial range throughout the water column, with the wavenumber kinetic energy spectra remaining unchanged upon assimilation. Overall, this comparison shows quantitatively that the 4D-Var system results in improved predictability as the analysis provides a smoother and more dynamically balanced fit between the observations and the model's time-evolving flow. This advocates the use of advanced, time-dependent data assimilation methods, particularly for highly variable oceanic regions, and motivates future work into further improving data assimilation schemes.

Assessment and Ensemble-Based Analysis of the Landfalling Typhoon Muifa (2022)

By considering the uncertainties in the initial field, model physical processes, and lateral boundary conditions, the Shanghai Weather And Risk Model System-Ensemble Prediction System (SWARMS-EN) is constructed. According to the prediction results of typhoon Muifa (2022), the daily track error of SWARMS-EN within 5 days is 70.6 km, 142.2 km, 129.1 km, 174.5 km, and 203.5 km, respectively. When compared with the Typhoon Ensemble Data Assimilation and Prediction System (TEDAPS) and the Global Ensemble Forecast System (GEFS) of the National Centers for Environmental Prediction (NCEP) in homogeneous conditions, SWARMS-EN performs better than TEDAPS within 72 h and better than GEFS beyond 72 h in track forecasting. This indicates an improvement in forecasting accuracy. The ensemble spread within two days is less than the root mean square error (RMSE), according to an analysis of the relationship between ensemble RMSE and spread, which shows that SWARMS-EN has no apparent systematic bias overall. The system has improved the ensemble RMSE and spread, indicating that it can better represent the uncertainty of the forecast and produce more reliable forecasts. Additionally, SWARMS-EN provides the landfall forecast five days in advance. The ensemble-based analysis suggests that the large-scale circulation is the primary factor contributing to the forecast differences among members, and the strong steering flow provides an indication of the landfalling forecast. The analysis of the ensemble characteristics of the initial field indicates that the initial perturbation between the wind field and the temperature field in the dynamically unstable region (such as near a tropical cyclone) exhibits flow dependence, and the small perturbation shows continuity throughout the entire troposphere. The distribution of ensemble spread and disturbance energy exhibited a reasonable growth stage as the forecast lead time increased. Disturbance internal energy dominated the lower troposphere, while the upper troposphere was mainly characterized by disturbance kinetic energy. Disturbance kinetic energy played a leading role in the evolution process. This conclusion further confirms the importance of paying attention to the initial small perturbations near TC in order to optimize the initial perturbation.

Improving Afternoon Thunderstorm Prediction over Complex Terrain with the Assimilation of Dense Ground-Based Observations: Four Cases in the Taipei Basin

Abstract In this study, we investigate the impact of assimilating densely distributed Global Navigation Satellite System (GNSS) zenith total delay (ZTD) and surface station (SFC) data on the prediction of very short-term heavy rainfall associated with afternoon thunderstorm (AT) events in the Taipei basin. Under weak synoptic-scale conditions, four cases characterized by different rainfall features are chosen for investigation. Experiments are conducted with a 3-h assimilation period, followed by 3-h forecasts. Also, various experiments are performed to explore the sensitivity of AT initialization. Data assimilation experiments are conducted with a convective-scale Weather Research and Forecasting–local ensemble transform Kalman filter (WRF-LETKF) system. The results show that ZTD assimilation can provide effective moisture corrections. Assimilating SFC wind and temperature data could additionally improve the near-surface convergence and cold bias, further increasing the impact of ZTD assimilation. Frequently assimilating SFC data every 10 min provides the best forecast performance especially for rainfall intensity predictions. Such a benefit could still be identified in the earlier forecast initialized 2 h before the start of the event. Detailed analysis of a case on 22 July 2019 reveals that frequent assimilation provides initial conditions that can lead to fast vertical expansion of the convection and trigger an intense AT. This study proposes a new metric using the fraction skill score to construct an informative diagram to evaluate the location and intensity of heavy rainfall forecast and display a clear characteristic of different cases. Issues of how assimilation strategies affect the impact of ground-based observations in a convective ensemble data assimilation system and AT development are also discussed. Significance Statement In this study, we investigate the impact of frequently assimilating densely distributed ground-based observations on predicting four afternoon thunderstorm events in the Taipei basin. While assimilating GNSS-ZTD data can improve the moisture fields for initializing convection, assimilating surface station data improves the prediction of rainfall location and intensity, particularly when surface data are assimilated at a very high frequency of 10 min.

Bias Correction of Tropical Cyclone Intensity for Ensemble Forecasts Using the XGBoost Method

Abstract In this study, the extreme gradient boosting (XGBoost) algorithm is used to correct tropical cyclone (TC) intensity in ensemble forecast data from the Typhoon Ensemble Data Assimilation and Prediction System (TEDAPS) at the Shanghai Typhoon Institute (STI), China Meteorological Administration (CMA). Results show that the forecast accuracy of TC intensity may be improved substantially using the XGBoost algorithm, especially when compared with a simple ensemble average of all members in the ensemble forecast [as depicted by the ensemble average (EnsAve) algorithm in this study]. The forecast errors for maximum wind speed (MWS) and minimum sea level pressure (MSLP) have been reduced by a significant margin, ranging from 6.3% to 18.4% for MWS and from 4% to 14.9% for MSLP, respectively. The performance of the XGBoost algorithm is overall better than that of the EnsAve algorithm, although there are a few samples when it is worse. The bias analysis shows that TEDAPS underpredicts the MWS and overpredicts the MSLP, meaning that the TEDAPS underestimates TC intensity. However, the XGBoost algorithm can reduce the bias to improve the forecast accuracy of TC intensity. Specifically, it achieves a reduction of over 20% in forecast errors for both the MWS and MSLP of typhoons compared to the EnsAve algorithm, indicating the XGBoost algorithm’s particular advantage in forecasting intense TCs. These results indicate that the TC intensity forecast can be substantially improved using the XGBoost algorithm, relative to the EnsAve algorithm.

Ingesting GOES-16 fire radiative power retrievals into Warn-on-Forecast System for Smoke (WoFS-Smoke)

Background The record number of wildfires in the United States in recent years has led to an increased focus on developing tools to accurately forecast their impacts at high spatial and temporal resolutions. Aims The Warn-on-Forecast System for Smoke (WoFS-Smoke) was developed to improve these forecasts using wildfire properties retrieved from satellites to generate smoke plumes in the system. Methods The WoFS is a regional domain ensemble data assimilation and forecasting system built around the concept of creating short-term (0–6 h) forecasts of high impact weather. This work extends WoFS-Smoke by ingesting data from the GOES-16 satellite at 15-min intervals to sample the rapidly changing conditions associated with wildfires. Key results Comparison of experiments with and without GOES-16 data show that ingesting high temporal frequency data allows for wildfires to be initiated in the model earlier, leading to improved smoke forecasts during their early phases. Decreasing smoke plume intensity associated with weakening fires was also better forecast. Conclusions The results were consistent for a large fire near Boulder, Colorado and a multi-fire event in Texas, Oklahoma, and Arkansas, indicating a broad applicability of this system. Implications The development of WoFS-Smoke using geostationary satellite data allows for a significant advancement in smoke forecasting and its downstream impacts such as reductions in air quality, visibility, and potentially properties of severe convection.

Precursory Analysis Ensemble Spread Signals That Foreshadow Stratospheric Sudden Warmings

Abstract This study conducts a thorough investigation into the behaviors of analysis ensemble spreads linked to stratospheric sudden warming (SSW) events. A stratosphere-resolving ensemble data assimilation system is used here to document the evolution of analysis spread leading up to a pair of warming events. Precursory signals of the increased ensemble spreads were found a few days prior to two SSW events that occurred during December 2018 and August–September 2019 in the Northern and Southern Hemispheres, respectively. The signals appeared in the upper and middle stratosphere and did not appear at lower heights. When the signals appeared, it was found that both tendency by forecast and analysis increment in a forecast-analysis (data assimilation) cycle simultaneously became large. An empirical orthogonal function analysis showed that the dominant structures of the precursory signals were equivalent barotropic and were 90° out-of-phase with the analysis ensemble-mean field. Over the same period, the upper and middle stratosphere became more susceptible to barotropic instability than in their previous states. We conclude that the differing growth of barotropically unstable modes across ensemble members can amplify spread during the lead-up to SSW events. Significance Statement Winds in the winter stratospheric polar vortex are typically westerly. Occasionally, however, warming over the pole leads to a reversal of the flow through a process known as stratospheric sudden warming. These events are difficult to predict even in state-of-the-art analysis and forecasting systems. In this study, we identify a precursor signal in the form of increased ensemble spread that appears to originate from differing realizations of growing barotropic modes across the ensemble. This signal could serve as a useful forecasting tool by enhancing situational awareness in the lead-up to potential stratospheric sudden warming events.

Performance of a hybrid gain ensemble data assimilation system based on the GRAPES_Meso model

In this paper, a hybrid gain data assimilation (HGDA) method was introduced into the mesoscale version of the Global/Regional Assimilation and Prediction System (GRAPES_Meso) model. To further improve the forecasting performance of the GRAPES_Meso model, a scheme combining the HGDA and multiscale incremental analysis update (IAU) methods was first proposed in this work. To evaluate the performance of the HGDA method in the GRAPES_Meso model, different initialization schemes, including the GRAPES three-dimensional variational (3DVAR), HGDA and HGDA multiscale IAU schemes, were compared for a tropical cyclone (TC) case (Lupit, 2021). The results showed that the HGDA scheme outperformed the GRAPES-3DVAR scheme in forecasting the TC position and intensity. Using the optimal relaxation times, the HGDA multiscale IAU scheme attained the best performance in forecasting the TC position, intensity and precipitation. In addition, 3-month quasioperational comparative experiments were conducted with different assimilation schemes. Rainfall validation showed that the HGDA multiscale IAU method yielded a higher threat score (TS) and equitable threat score (ETS) and a lower bias score than did the 3DVAR multiscale IAU method. Especially in the first 18 h of the forecast, the HGDA scheme could increase the TS and ETS values by 11% and 6%, respectively, compared to the control and 3DVAR experiments. The results demonstrated that the HGDA method resulted in better prediction performance of the GRAPES_Meso model than did the GRAPES-3DVAR method.

Evaluating the Multiscale Implementation of Valid Time Shifting within a Real-Time EnVar Data Assimilation and Forecast System for the 2022 HWT Spring Forecasting Experiment

Abstract Multiscale valid time shifting (VTS) was explored for a real-time convection-allowing ensemble (CAE) data assimilation (DA) system featuring hourly assimilation of conventional in situ and radar reflectivity observations, developed by the Multiscale Data Assimilation and Predictability Laboratory. VTS triples the base ensemble size using two subensembles containing member forecast output before and after the analysis time. Three configurations were tested with 108-member VTS-expanded ensembles: VTS for individual mesoscale conventional DA (ConVTS) or storm-scale radar DA (RadVTS) and VTS integrated to both DA components (BothVTS). Systematic verification demonstrated that BothVTS matched the DA spread and accuracy of the best-performing individual component VTS. The 10-member forecasts showed BothVTS performs similarly to ConVTS, with RadVTS having better skill in 1-h precipitation at forecast hours 1–6, while Both/ConVTS had better skill at later hours 7–15. An objective splitting of cases by 2-m temperature cold bias revealed RadVTS was more skillful than Both/ConVTS out to hour 10 for cold-biased cases, while BothVTS performed best at most hours for less-biased cases. A sensitivity experiment demonstrated improved performance of BothVTS when reducing the underlying model cold bias. Diagnostics revealed enhanced spurious convection of BothVTS for cold-biased cases was tied to larger analysis increments in temperature than moisture, resulting in erroneously high convective instability. This study is the first to examine the benefits of a multiscale VTS implementation, showing that BothVTS can be utilized to improve the overall performance of a multiscale CAE system. Further, these results underscore the need to limit biases within a DA and forecast system to best take advantage of VTS analysis benefits.

Global Scale Inversions from MOPITT CO and MODIS AOD

Top-down observational constraints on emissions flux estimates from satellite observations of chemical composition are subject to biases and errors stemming from transport, chemistry and prior emissions estimates. In this context, we developed an ensemble data assimilation system to optimize the initial conditions for carbon monoxide (CO) and aerosols, while also quantifying the respective emission fluxes with a distinct attribution of anthropogenic and wildfire sources. We present the separate assimilation of CO profile v9 retrievals from the Measurements of Pollution in the Troposphere (MOPITT) instrument and Aerosol Optical Depth (AOD), collection 6.1, from the Moderate Resolution Imaging Spectroradiometer (MODIS) instruments. This assimilation system is built on the Data Assimilation Research Testbed (DART) and includes a meteorological ensemble to assimilate weather observations within the online Community Atmosphere Model with Chemistry (CAM-chem). Inversions indicate an underestimation of CO emissions in CAMS-GLOB-ANT_v5.1 in China for 2015 and an overestimation of CO emissions in the Fire INventory from NCAR (FINN) version 2.2, especially in the tropics. These emissions increments are consistent between the MODIS AOD and the MOPITT CO-based inversions. Additional simulations and comparison with in situ observations from the NASA Atmospheric Tomography Mission (ATom) show that biases in hydroxyl radical (OH) chemistry dominate the CO errors.

Making the Case for High-Resolution Regional Ocean Reanalyses: An Example with the Red Sea

Abstract Coastal oceans host 40% of the world population and amount to $1.5 trillion of the global economy. Studying, managing, and developing the coastal regions require decades-long information about their environment. Long-term ocean measurements are, however, lacking for most coastal regions and often global reanalyses are used instead. These are, however, coarse in nature and tuned for the global circulations. The Red Sea (RS) is a narrow basin connected to the Indian Ocean through the Bab-al-Mandab Strait. Despite being the busiest commercial crossroad and hosting the world’s third largest coral reef system, the RS lacks long-term observations. A recent increase in population and an unprecedented acceleration in governmental and industrial developments further emphasized the need for long-term datasets to support its development and the sustainability of its habitats, and to understand its response to a changing climate. Toward this end, we have generated a 20-yr high-resolution reanalysis for the RS (RSRA) using a state-of-the-art ensemble data assimilation system incorporating available observations. Compared to global reanalyses, RSRA provides a markedly better description of the RS general and mesoscale circulation features, their variability, and trends. In particular, RSRA accurately captures the three-layer summer transport through the Bab-al-Mandab, simulated as two-layer transport by some global reanalyses. It further reproduces the seasonal anomalies, whereas global reanalyses misidentify some seasons as anomalous. Global reanalyses further overestimate the interannual variations in salinity, misrepresent the trend in temperature, and underestimate the trend in sea level. Our study clearly emphasizes the importance of generating dedicated high-resolution regional ocean reanalyses.

Development of a Polar Mesocyclone and Associated Environmental Characteristics During the Heavy Snowfall Event in Sapporo, Japan, in Early February 2022

AbstractIn early February 2022, the Sapporo region, north Japan, experienced record‐breaking snowfall. During this period, a mesoscale cyclonic rotation corresponding to a polar mesocyclone (PMC) was observed over the Sea of Japan by a geostationary satellite. This study investigates important processes involved in this heavy snowfall event with attention to the role of the PMC using a 100‐member regional ensemble forecast. The initial conditions for the ensemble forecast are obtained by a regional‐scale ensemble data assimilation system that assimilates conventional observations only. Several ensemble members successfully predict heavy snowfall in the Sapporo region with an 18 hr lead time. Composite analyses reveal that these better ensemble members predict that a southward‐moving PMC with a precipitation system develops from small vortices over the ocean and approaches to the Sapporo region, whereas worse members instead predict precipitation systems shifted to the north of the Sapporo region. Further composite analyses demonstrate that the better ensemble members have environmental characteristics such as a deeper and colder upper‐level trough and large low‐level baroclinicity, which are favorable for the development of PMCs. In particular, the latter is associated with reverse shear in which PMCs tend to move southward. Therefore, these environmental factors are key to predicting the heavy snowfall event in Sapporo associated with the PMC.

C-Coupler3.0: an integrated coupler infrastructure for Earth system modelling

Abstract. The community coupler (C-Coupler) for Earth system modelling is a coupler family that was developed in China in 2010. C-Coupler3.0, the latest version, is fully compatible with the previous version, C-Coupler2, and is an integrated infrastructure with new features, i.e. a series of parallel-optimization technologies for accelerating coupling initialization and reducing memory usage, a common halo-exchange library for developing a parallel version of a model, a common module-integration framework for integrating a software module (e.g. a flux algorithm, a parameterization scheme, and a data assimilation method), a common framework for conveniently developing a weakly coupled ensemble data assimilation system, and a common framework for flexibly inputting and outputting fields in parallel. Specifically, C-Coupler3.0 is able to handle coupling under much finer resolutions (e.g. more than 100 million horizontal grid cells) with fast coupling initialization and successful generation of remapping-weight files.

Arctic sea ice data assimilation combining an ensemble Kalman filter with a novel Lagrangian sea ice model for the winter 2019–2020

Abstract. Advanced data assimilation (DA) methods, widely used in geophysical and climate studies to merge observations with numerical models, can improve state estimates and consequent forecasts. We interface the deterministic ensemble Kalman filter (DEnKF) to the Lagrangian neXt generation Sea Ice Model, neXtSIM. The ensemble is generated by perturbing the atmospheric and oceanic forcing throughout the simulations and randomly initialized ice cohesion. Our ensemble–DA system assimilates sea ice concentration (SIC) from the Ocean and Sea Ice Satellite Application Facility (OSI-SAF) and sea ice thickness (SIT) from the merged CryoSat-2 and SMOS datasets (CS2SMOS). Because neXtSIM is computationally solved on a time-dependent evolving mesh, it is a challenging application for ensemble–DA. As a solution, we perform the DEnKF analysis on a fixed and regular reference mesh, on which model variables are interpolated before the DA and then back to each member's mesh after the DA. We evaluate the impact of assimilating different types of sea ice observations on the model's forecast skills of the Arctic sea ice by comparing satellite observations and a free-run ensemble in an Arctic winter period, 2019–2020. Significant improvements in modeled SIT indicate the importance of assimilating weekly CS2SMOS SIT, while the improvements of SIC and ice extent are moderate but benefit from daily ingestion of the OSI-SAF SIC. For most of the winter, the correlation between SIT and SIC is weaker, which results in little cross-inference between the two variables in the assimilation step. Overall, the ensemble–DA system based on the stand-alone sea ice model demonstrates the feasibility of winter Arctic sea ice prediction with good computational efficiency. These results open the path toward operational implementation and the extension to multi-year assimilation.

Performance of a Convective-Scale Ensemble Prediction System on 2017 Warm-Season Afternoon Thunderstorms over Taiwan

Abstract A 16-member convective-scale ensemble prediction system (CEPS) developed at the Central Weather Bureau (CWB) of Taiwan is evaluated for probability forecasts of convective precipitation. To address the issues of limited predictability of convective systems, the CEPS provides short-range forecasts using initial conditions from a rapid-updated ensemble data assimilation system. This study aims to identify the behavior of the CEPS forecasts, especially the impact of different ensemble configurations and forecast lead times. Warm-season afternoon thunderstorms (ATs) from 30 June to 4 July 2017 are selected. Since ATs usually occur between 1300 and 2000 LST, this study compares deterministic and probabilistic quantitative precipitation forecasts (QPFs) launched at 0500, 0800, and 1100 LST. This study demonstrates that initial and boundary perturbations (IBP) are crucial to ensure good spread–skill consistency over the 18-h forecasts. On top of IBP, additional model perturbations have insignificant impacts on upper-air and precipitation forecasts. The deterministic QPFs launched at 1100 LST outperform those launched at 0500 and 0800 LST, likely because the most-recent data assimilation analyses enhance the practical predictability. However, it cannot improve the probabilistic QPFs launched at 1100 LST due to inadequate ensemble spreads resulting from limited error growth time. This study points out the importance of sufficient initial condition uncertainty on short-range probabilistic forecasts to exploit the benefits of rapid-update data assimilation analyses. Significance Statement This study aims to understand the behavior of convective-scale short-range probabilistic forecasts in Taiwan and the surrounding area. Taiwan is influenced by diverse weather systems, including typhoons, mei-yu fronts, and local thunderstorms. During the past decade, there has been promising improvement in predicting mesoscale weather systems (e.g., typhoons and mei-yu fronts). However, it is still challenging to provide timely and accurate forecasts for rapid-evolving high-impact convection. This study provides a reference for the designation of convective-scale ensemble prediction systems; in particular, those with a goal to provide short-range probabilistic forecasts. While the findings cannot be extrapolated to all ensemble prediction systems, this study demonstrates that initial and boundary perturbations are the most important factors, while the model perturbation has an insignificant effect. This study suggests that in-depth studies are required to improve the convective-scale initial condition accuracy and uncertainty to provide reliable probabilistic forecasts within short lead times.

Understanding the impact of assimilating FORMOSAT‐7/COSMIC‐2 radio occultation refractivity on tropical cyclone genesis: Observing system simulation experiments using Hurricane Gordon (2006) as a case study

AbstractStudies have shown that assimilating the radio occultation (RO) observations, including those from the FORMOSAT‐3/COSMIC (constellation observing systems for meteorology, ionosphere, and climate) (FS3‐C), provides positive impacts on tropical cyclone (TC) forecasts. The FS3‐C's successor, the FORMOSAT‐7/COSMIC‐2 (FS7‐C2), provides denser spatial data coverage over the Tropics and Subtropics, where severe weather systems often occur. This study investigates the impact of FS7‐C2 refractivity profiles on the prediction of TC genesis. A quick observing system simulation experiment is conducted for the period when Hurricanes Helene and Gordon (2006) occurred over the North Atlantic Ocean using a regional ensemble data assimilation system. Though assimilating FS3‐C or FS7‐C2 ROs successfully reproduces Helene's development, assimilating FS7‐C2 ROs better captures the genesis and development of Gordon with abundant moisture and vorticity in Gordon's core region, providing conditions favorable for the development of deep convection. A minimum area‐mean total precipitable water vapor of 54 mm, as well as the existence of mid‐level cyclonic vorticity (e.g., 500 hPa), at the storm core region in the initial condition is required for forecasting Gordon's genesis. Also, the assimilation of FS7‐C2 ROs in our experiments reduces the 500 hPa geopotential error by 22% and improves probabilistic quantitative precipitation forecast compared with assimilating FS3‐C ROs. Two sensitivity tests are conducted to evaluate the impact of low‐level negatively biased FS7‐C2 RO profiles and the removal of FS7‐C2 data below 3 km on Gordon's genesis. The former test does not degrade Gordon's genesis forecast skills due to a dipole error correlation between the background ROs and the moisture field over an observed RO profile near Gordon. The latter test does degrade Gordon's forecast skills but is still better than the assimilation of FS3‐C ROs since the features of low‐level moisture and mid‐level vorticity are preserved to some extent.

Accounting for Correlated AMV Satellite Observation Errors in the Ensemble Data Assimilation System

Accounting for Correlated AMV Satellite Observation Errors in the Ensemble Data Assimilation System