In a global numerical weather prediction (NWP) modeling framework we study the implementation of Gaussian uncertainty of individual particles into the assimilation step of a localized adaptive particle filter (LAPF). We obtain a local representation of the prior distribution as a mixture of basis functions. In the assimilation step, the filter calculates the individual weight coefficients and new particle locations. It can be viewed as a combination of the LAPF and a localized version of a Gaussian mixture filter, i.e., a Localized Mixture Coefficients Particle Filter (LMCPF).
Here, we investigate the feasibility of the LMCPF within a global operational framework and evaluate the relationship between prior and posterior distributions and observations. Our simulations are carried out in a standard pre-operational experimental set-up with the full global observing system, 52 km global resolution and 106 model variables. Statistics of particle movement in the assimilation step are calculated. The mixture approach is able to deal with the discrepancy between prior distributions and observation location in a real-world framework and to pull the particles towards the observations in a much better way than the pure LAPF. This shows that using Gaussian uncertainty can be an important tool to improve the analysis and forecast quality in a particle filter framework.
This paper describes a particle filter in the global NWP at Deutscher Wetterdienst (DWD). A particle filter (PF) in the global NWP at DWD is proposed and evaluated its skills in comparison with the operational system. To alleviate the degeneration, which is the largest issue in PFs with
high-dimensional systems, several approaches are effectively
incorporated such as localization, Gaussian mixture approximation in the
prior distribution, adaptive resampling, and so on (See Section 2.3). Since comprehensive formulations in this system are described, the readers can totally understand its theoretical aspects.
A new operational seasonal forecast system, Japan Meteorological Agency (JMA)/Meteorological Research Institute (MRI) Coupled Prediction System (CPS) version 3 (JMA/MRI–CPS3), has been developed. This system represents a major upgrade of the former system, CPS2. CPS3 comprises atmosphere, land, ocean, and sea ice forecast models and the necessary initialization systems for these models. For historical reforecasts, the atmospheric reanalysis dataset JRA-3Q provides initial conditions for the atmosphere and the external forcings for land, ocean, and sea ice analysis. In the operational forecast, JMA's operational atmospheric analysis is used in conjunction with JRA-3Q to initialize the system in near-real time. The land surface model is initialized using an uncoupled free simulation, forced by the atmospheric analysis. The ocean and sea ice models are initialized with the global ocean data assimilation system MOVE-G3, which incorporates a newly developed four-dimensional variational method for temperature, salinity, and sea surface height and a three-dimensional method for sea ice concentration. Compared with the previous system, the CPS3 forecast model components have approximately 2-4 times higher resolution: the atmosphere and land models are configured with ∼ 55 km horizontal resolution, with 100 vertical atmosphere layers; and the ocean and sea ice models have a resolution of 0.25° × 0.25°, with 60 vertical ocean layers. The physical processes of the atmosphere are greatly refined in CPS3 relative to CPS2, resulting in improved representation of sub-seasonal to seasonal scale variability, including the eastward propagation of the Madden–Julian Oscillation, winter blocking highs in the North Atlantic, and coupled atmosphere–ocean variability during El Niño–Southern Oscillation events. Our historical reforecast experiment for 1991-2020 suggests that CPS3 has greater forecast skill than CPS2. The usability of the model output has been improved in CPS3 by reorganizing the operation schedule to provide daily updates of five-member ensemble forecasts.
This paper describe a newly developed operational seasonal forecast system, JMA/MRI-CPS3. Ocean 4D-Var and sea ice 3D-Var data assimilation methods are newly introduced. The errors in the ocean analysis are now represented in the initial perturbations. Updated physical processes and increased resolution of the atmospheric model contribute to the improved climate reproducibility of the MJO and North Atlantic blocking highs. The introduction of a 0.25-degree-resolution ocean model provides a realistic representation of tropical instability waves and contributes to improved ENSO pattern.
In drylands, the dry vegetation coverage affects dust occurrence by modulating threshold friction velocity (or wind speed) for dust emission. However, there has been little research into quantifying the effect of dry vegetation coverage on dust occurrence. This study investigated spatial and temporal variations of dust occurrence and three definitions of strong wind frequency over the Gobi Desert and surrounding regions in March and April, months when dust occurrence is frequent, during 2001-2021. We evaluated the effects of variations in dry vegetation on dust occurrence by using the threat scores of forecasted dust occurrences for each strong wind definition. Our results indicate that dry vegetation, which was derived from the MODIS Soil Tillage Index, affects dust occurrence more remarkably in April than in March. In March, land surface parameters such as soil freeze-thaw and snow cover, in addition to dry vegetation coverage, should be considered to explain dust variations in that month. However, use of the threshold wind speed estimated from dry vegetation coverage improved the prediction accuracy of dust occurrence in April. Therefore, we propose that the dry vegetation coverage is a key factor controlling dust occurrence variations in April. The findings imply that estimation of dry vegetation coverage should be applied to dust models.
This paper proposed a machine learning method as an observation
operator for satellite radiances within a data assimilation system.
Model forecast and satellite microwave radiance observations
are used to train machine learning models to obtain the observation
operator for satellite data assimilation.Data assimilation experiments using the machine learning-based
observation operator show promising results without a separate bias
correction procedure.The machine learning-based observation operator can potentially
accelerate the development of using new satellite observations in
numerical weather prediction.
The observation operator (OO) is essential in data assimilation (DA) to derive the model equivalent of observations from the model variables. In the satellite DA, the OO for satellite microwave brightness temperature (BT) is usually based on the radiative transfer model (RTM) with a bias correction procedure. To explore the possibility to obtain OO without using physically based RTM, this study applied machine learning (ML) as OO (ML-OO) to assimilate BT from Advanced Microwave Sounding Unit-A (AMSU-A) channels 6 and 7 over oceans and channel 8 over both land and oceans under clear-sky conditions. We used a reference system, consisting of the nonhydrostatic icosahedral atmospheric model (NICAM) and the local ensemble transform Kalman filter (LETKF). The radiative transfer for TOVS (RTTOV) was implemented in the system as OO, combined with a separate bias correction procedure (RTTOV-OO). The DA experiment was performed for one month to assimilate conventional observations and BT using the reference system. Model forecasts from the experiment were paired with observations for training the ML models to obtain ML-OO. In addition, three DA experiments were conducted, which revealed that DA of the conventional observations and BT using ML-OO was slightly inferior, compared to that of RTTOV-OO, but it was better than the assimilation based on only conventional observations. Moreover, ML-OO treated bias internally, thereby simplifying the overall system framework. The proposed ML-OO has limitations due to (1) the inability to treat bias realistically when a significant change is present in the satellite characteristics, (2) inapplicability for many channels, (3) deteriorated performance, compared with that of RTTOV-OO in terms of accuracy and computational speed, and (4) physically based RTM is still used to train the ML-OO. Future studies can alleviate these drawbacks, thereby improving the proposed ML-OO.
This paper proposed a new method calculating the threshold wind speed for dust occurrence.
A new method to obtain threshold wind speed that takes account
of the interannual variations of dry vegetation cover is proposed. Dry vegetation coverage is a key factor determining interannual variations in the April dust occurrence. Other land surface factors such as soil freeze-thaw and snow
cover should be considered to explain dust occurrence variations in
March.
In the discretization of the primitive equations for numerical calculations, a formulation of a three-dimensional spectral model that uses the spectral method not only in the horizontal direction but also in the vertical direction is proposed. In this formulation, the Legendre polynomial expansion is used for the vertical discretization. It is shown that semi-implicit time integration can be efficiently done under this formulation. Then, a numerical model based on this formulation is developed and several benchmark numerical calculations proposed in previous studies are performed to show that this implementation of the primitive equations can give accurate numerical solutions with a relatively small degrees of freedom in the vertical discretization. It is also shown that, by performing several calculations with different vertical degrees of freedom, a characteristic property of the spectral method is observed in which the error of the numerical solution decreases rapidly when the number of vertical degrees of freedom is increased. It is also noted that an alternative to the sponge layer can be devised to suppress the reflected waves under this formulation, and that a “toy” model can be derived as an application of this formulation, in which the vertical degrees of freedom are reduced to the minimum.
Ishioka et al., (2022): The above paper was chosen as a JMSJ Editor's Highlight. (7 Mar. 2022)
Influence of the Global Warming on Tropical Cyclone Climatology: An Experiment with the JMA Global Model
Released on J-STAGE: July 03, 2002 | Volume 80 Issue 2 Pages 249-272
Masato SUGI, Akira NODA, Nobuo SATO
Views: 374
Structure and Evolution of Precipitation Cores in an Isolated Convective Storm Observed by Phased Array Weather Radar
Released on J-STAGE: June 14, 2021 | Volume 99 Issue 3 Pages 765-784
Yukie MORODA, Kazuhisa TSUBOKI, Shinsuke SATOH, Katsuhiro NAKAGAWA, Tomoo USHIO, Shingo SHIMIZU
Views: 240
The JRA-55 Reanalysis: General Specifications and Basic Characteristics
Released on J-STAGE: March 18, 2015 | Volume 93 Issue 1 Pages 5-48
Shinya KOBAYASHI, Yukinari OTA, Yayoi HARADA, Ayataka EBITA, Masami MORIYA, Hirokatsu ONODA, Kazutoshi ONOGI, Hirotaka KAMAHORI, Chiaki KOBAYASHI, Hirokazu ENDO, Kengo MIYAOKA, Kiyotoshi TAKAHASHI
Views: 202
Development of Synoptic-Scale Disturbances over the Tropical Western North Pacific Associated with the Boreal Summer Intraseasonal Oscillation and the Interannual Pacific–Japan Pattern
Released on J-STAGE: March 07, 2023 | Volume 101 Issue 2 Pages 103-123
Ayako SEIKI, Yu KOSAKA, Satoru YOKOI
Views: 169
Japan Meteorological Agency/Meteorological Research Institute Coupled Prediction System Version 3 (JMA/MRI–CPS3)
Released on J-STAGE: March 14, 2023 | Volume 101 Issue 2 Pages 149-169
Shoji HIRAHARA, Yutaro KUBO, Takuma YOSHIDA, Takuya KOMORI, Jotaro CHIBA, Toshinari TAKAKURA, Takafumi KANEHAMA, Ryohei SEKIGUCHI, Kenta OCHI, Hiroyuki SUGIMOTO, Yukimasa ADACHI, Ichiro ISHIKAWA, Yosuke FUJII
Views: 163