Estimation of AMSU-A Radiance Observation Impacts in an LETKF-Based Atmospheric Global Data Assimilation System: Comparison with EFSO and Observing System Experiments

Akira Yamazaki aApplication Laboratory, Japan Agency for Marine-Earth Science and Technology, Yokohama, Japan

Search for other papers by Akira Yamazaki in
Current site
Google Scholar
PubMed
Close
,
Koji Terasaki bMeteorological Research Institute, Japan Meteorological Agency, Tsukuba, Japan
cRIKEN Center for Computational Science, Kobe, Japan

Search for other papers by Koji Terasaki in
Current site
Google Scholar
PubMed
Close
,
Takemasa Miyoshi cRIKEN Center for Computational Science, Kobe, Japan
aApplication Laboratory, Japan Agency for Marine-Earth Science and Technology, Yokohama, Japan

Search for other papers by Takemasa Miyoshi in
Current site
Google Scholar
PubMed
Close
, and
Shunsuke Noguchi dDepartment of Earth and Planetary Sciences, Kyushu University, Fukuoka, Japan
eResearch Center for Environmental Modeling and Application, JAMSTEC, Yokohama, Japan

Search for other papers by Shunsuke Noguchi in
Current site
Google Scholar
PubMed
Close
Open access

Abstract

This work assesses the contribution of assimilating AMSU-A satellite-based radiance measurements to a global data assimilation system based on an atmospheric general circulation model and the local ensemble transform Kalman filter (LETKF). The radiance measurements were from three channels that are sensitive to the upper troposphere and lower stratosphere. The contribution of these measurements, or AMSU-A observation impact, was estimated both through ensemble-based forecast sensitivity to observations (EFSO) and observing system experiments (OSEs). Two streams of data-denial experiments for the AMSU-A observations were performed for about one month during winter in each hemisphere. The OSEs quantified the accumulated observation impact by cycling (repeating) data denials: including AMSU-A observations reduced the total observation impact for all observations of each data assimilation cycle. In contrast, EFSO estimated AMSU-A to increase the total observation impact. The opposing effects were attributed to the accumulated observation impact in the OSEs; the accumulation could stabilize the data assimilation cycles. In both experiments, the accumulated observation impact of AMSU-A was strongest in the upper troposphere, particularly in the austral midlatitudes where westerly jets exist and observations of other types are sparse. EFSO also assessed AMSU-A to have the most beneficial observation impact in similar locations. The AMSU-A observation impact tended to accumulate just downstream of where EFSO estimated the beneficial observation impact signals. The accumulated AMSU-A observation impact was tied to dynamic processes in the upper-tropospheric and general stratospheric circulation. Therefore, EFSO helps estimate the beneficial distributions of AMSU-A accumulated observation impact by considering their dynamical propagation.

Significance Statement

The Advanced Microwave Sounding Unit-A (AMSU-A) satellite radiance assimilation technique was successfully integrated into the Atmospheric General Circulation Model for the Earth Simulator (AFES)–LETKF data assimilation system. We conducted OSEs and used EFSO to assess the AMSU-A observation impact. The two estimation methods identified opposite observation impacts due to the cycling (repeating) OSEs of the AMSU-A observations. We interpreted the causes of the opposite estimations. However, even for the cycling OSEs, EFSO appeared to help estimate distributions of the accumulated observation impact. It is important to consider the dynamical propagation of accumulated observation impact in general circulation.

© 2023 American Meteorological Society. For information regarding reuse of this content and general copyright information, consult the AMS Copyright Policy (www.ametsoc.org/PUBSReuseLicenses).

Corresponding author: Akira Yamazaki, yzaki@jamstec.go.jp

Abstract

This work assesses the contribution of assimilating AMSU-A satellite-based radiance measurements to a global data assimilation system based on an atmospheric general circulation model and the local ensemble transform Kalman filter (LETKF). The radiance measurements were from three channels that are sensitive to the upper troposphere and lower stratosphere. The contribution of these measurements, or AMSU-A observation impact, was estimated both through ensemble-based forecast sensitivity to observations (EFSO) and observing system experiments (OSEs). Two streams of data-denial experiments for the AMSU-A observations were performed for about one month during winter in each hemisphere. The OSEs quantified the accumulated observation impact by cycling (repeating) data denials: including AMSU-A observations reduced the total observation impact for all observations of each data assimilation cycle. In contrast, EFSO estimated AMSU-A to increase the total observation impact. The opposing effects were attributed to the accumulated observation impact in the OSEs; the accumulation could stabilize the data assimilation cycles. In both experiments, the accumulated observation impact of AMSU-A was strongest in the upper troposphere, particularly in the austral midlatitudes where westerly jets exist and observations of other types are sparse. EFSO also assessed AMSU-A to have the most beneficial observation impact in similar locations. The AMSU-A observation impact tended to accumulate just downstream of where EFSO estimated the beneficial observation impact signals. The accumulated AMSU-A observation impact was tied to dynamic processes in the upper-tropospheric and general stratospheric circulation. Therefore, EFSO helps estimate the beneficial distributions of AMSU-A accumulated observation impact by considering their dynamical propagation.

Significance Statement

The Advanced Microwave Sounding Unit-A (AMSU-A) satellite radiance assimilation technique was successfully integrated into the Atmospheric General Circulation Model for the Earth Simulator (AFES)–LETKF data assimilation system. We conducted OSEs and used EFSO to assess the AMSU-A observation impact. The two estimation methods identified opposite observation impacts due to the cycling (repeating) OSEs of the AMSU-A observations. We interpreted the causes of the opposite estimations. However, even for the cycling OSEs, EFSO appeared to help estimate distributions of the accumulated observation impact. It is important to consider the dynamical propagation of accumulated observation impact in general circulation.

© 2023 American Meteorological Society. For information regarding reuse of this content and general copyright information, consult the AMS Copyright Policy (www.ametsoc.org/PUBSReuseLicenses).

Corresponding author: Akira Yamazaki, yzaki@jamstec.go.jp

1. Introduction

Global data assimilation systems collate data from various observation platforms such as radiosonde, aircraft, and satellites. Satellite-based radiance measurements are essential to improving data assimilation systems. An example is the Advanced Microwave Sounding Unit-A (AMSU-A), which has a particularly great observation impact (OI) (e.g., Gelaro et al. 2010; Ota et al. 2013; Hotta et al. 2017; Privé et al. 2021); therefore, incorporating AMSU-A observations is crucial to improving data assimilation systems.

Optimizing a global observing system requires proactive estimation of the contribution, or OI, of any newly added observation type. It would be useful to determine the extent and geographic location of any OI from a particular set of measurements, including AMSU-A observations. For this purpose, many studies have conducted observing system experiments (OSEs). Other recent studies have measured OI using the adjoint-based forecast sensitivity to observations (FSO; Langland and Baker 2004; Gelaro et al. 2010; Privé et al. 2021). FSO mathematically derives the impacts of all observation types by quantifying the extent to which each improves or degrades subsequent forecasts. Recent studies have applied the FSO formulation to the ensemble Kalman filter system in a diagnostic scheme known as EFSO (Kalnay et al. 2012; Ota et al. 2013; Hotta et al. 2017; Kotsuki et al. 2019). As the EFSO scheme has been successfully implemented in the Atmospheric General Circulation Model for the Earth Simulator (AFES)–local ensemble transform Kalman filter (LETKF) ensemble data assimilation system (ALEDAS; Yamazaki et al. 2021), we can assess AMSU-A OI without conducting an OSE. Here, the OI calculated by EFSO is denoted as OIEFSO. Several studies have sought to understand the extent to which (E)FSO methods can quantify OI derived by OSEs. It has usually been reported that (E)FSO performs well in estimating the OI obtained by OSEs in most situations (e.g., Hotta et al. 2017; Lawrence et al. 2019; Privé et al. 2021; Yamazaki et al. 2021). However, the (E)FSO formulation does not directly approximate OI when data denials/additions are repeated in time in an OSE (e.g., Lien et al. 2018).

Previous studies have compared the OIs derived by OSEs and by (E)FSO (e.g., Gelaro and Zhu 2009; Lawrence et al. 2019; Candy et al. 2021; Eyre 2021). Gelaro and Zhu (2009) found that FSO overall provides estimates consistent with OSEs for most observation types; they gave detailed comparisons of the OIs from FSO and OSEs in the sense of global and area-averaged values. On the other hand, Lien et al. (2018) reported that the OI estimated by EFSO cannot correspond to that from OSEs, as OSEs usually include cycling data assimilation experiments, as explained below. They concluded that an OSE with a long period is needed to collect a sufficiently large sample and that (E)FSO represents a practically useful alternative to OSEs. However, to our best knowledge, previous studies have not provided a full physical interpretation of the origins of the differences between OSEs and (E)FSO.

(E)FSO is known to estimate OI by noncycling data denials (Lien et al. 2018; Yamazaki et al. 2021). For noncycling OSEs, each data denial for particular observations is carried out once every analysis cycle. The OI is assessed by the forecast and analysis data from the control (CTL) experiment and the forecast from the OSE analysis. In noncycling OSEs, the OI can be approximated by (E)FSO (section 2b). However, cycling OSEs repeatedly exclude specific observations from each analysis cycle. By cycling data denials, the OSE analysis can depart from the CTL analysis owing to the accumulation of the OI. We call the OI in cycling OSEs the “accumulated OI.” It is not evident whether (E)FSO can estimate the accumulated OI. Practically, it would make sense to consider cycling OSEs when conducting near-stationary observational campaigns (e.g., Yamazaki et al. 2015; Inoue et al. 2015; Sato et al. 2018; Hattori et al. 2017), in which specific observations and/or observation types are repeatedly added or deleted. Feasibility studies by Chen and Kalnay (2019) and Chen and Kalnay (2020) explored how to adopt EFSO for online quality control (called proactive quality control) in cycling data assimilation systems. They reported that EFSO is helpful in improving data assimilation systems even in cycling situations through proactive quality control. Lien et al. (2018) applied EFSO as an alternative to OSEs to design newly added observation types with offline quality control. Those studies did not evaluate the differences in the distributions between accumulated OI by OSEs and OIEFSO. Therefore, we focused on the spatial distributions of AMSU-A OIEFSO and accumulated AMSU-A OI and the physical interpretation of the accumulation of AMSU-A OI during cycling data assimilation. Hence, in this study, we attempted a detailed comparison of the differences between OIEFSO and accumulated OI when cycling OSEs for AMSU-A observations are performed. To interpret the difference physically, the distribution and propagation of OI in the OSEs through the general circulations were also discussed. Physical interpretation can aid exploration of the applicability of (E)FSO.

In addition to the main purpose of this study, we also tested whether an AMSU-A satellite radiance assimilation scheme was successfully implemented in our data assimilation system. We recently restructured an in-house global atmospheric data assimilation system called ALEDAS (Enomoto et al. 2013). The updated system includes an increased model top to resolve the entire stratospheric circulation and the ability to assimilate AMSU-A satellite radiances (Terasaki and Miyoshi 2017). Here we study the feasibility of the latter update (AMSU-A satellite radiance assimilation).1 The updated system implements an AMSU-A assimilation scheme similar to that of Terasaki and Miyoshi (2017), which was initially developed for LETKF-based data assimilation systems (see also Miyoshi and Sato 2007; Miyoshi et al. 2010). Terasaki and Miyoshi (2017)’s scheme involves an observation operator that converts a model vertical profile to brightness temperature using RTTOV2 (version 11.3; Saunders et al. 2013) and online bias corrections of both the airmass and scan biases. We assimilate three channels of the AMSU-A radiances that are sensitive to temperature around the upper troposphere. It is a crucial step for an ensemble Kalman filter data assimilation system to assimilate satellite radiance data. Houtekamer and Zhang (2016)’s review noted that there are still few studies of the feasibility of assimilating satellite radiance observations in ensemble Kalman filter–based data assimilation systems. Therefore, our study may constitute a first step to studying the impact of AMSU-A satellite radiance observations on the tropospheric and stratospheric circulations using the same AMSU-A observations and assimilation scheme of Terasaki and Miyoshi (2017).

We conducted AMSU-A OSEs (data-denial experiments) to assess the effects of accumulated AMSU-A OI on the general circulation in a global ensemble Kalman filter data assimilation system (ALEDAS). Furthermore, we discussed how EFSO can help estimate the accumulated AMSU-A OI.

The remainder of this paper is organized as follows. Section 2 introduces the data assimilation system, the definition of accumulated OI, and the experimental designs of the OSEs. Results of the OSEs and the physical interpretation of the difference between OIEFSO and accumulated OI are given in section 3. Section 4 discusses how accumulated AMSU-A OI propagates in the general circulation reproduced in ALEDAS. Finally, conclusions and final remarks are provided in section 5.

2. Data, definition, and experimental designs

a. Development and performance of ALEDAS

We updated ALEDAS to resolve stratospheric circulation and enable assimilation of AMSU-A satellite radiances. ALEDAS is briefly described here. An atmospheric general circulation model called AFES is used for the model part. The horizontal resolution is T119 (∼1° × 1°), and there are 56 vertical levels with the model top of ∼0.1 hPa. The cumulus convection and radiation schemes follow those of Emanuel (1991) and Sekiguchi and Nakajima (2008), respectively. The nonorographic gravity wave drag scheme of Scinocca (2003) is newly implemented (Orr et al. 2010; Baba 2019; Yamazaki and Noguchi 2023, manuscript submitted to Mon. Wea. Rev., hereafter YN23). Detailed AFES information is given by Enomoto et al. (2008, 2013), and Kodama et al. (2019). The ocean boundary condition is given by the daily 1/4° Optimum Interpolation Sea Surface Temperature (Reynolds et al. 2007). Data assimilation applies to the LETKF with an ensemble size of 63. The localization functions are 400 km horizontally and 0.4lnp vertically. The spatially uniform 10% multiplicative inflation method is adopted. These settings follow those of Enomoto et al. (2013) and Yamazaki et al. (2021). Other observations derived from the PrepBUFR3 datasets include conventional types, aircraft observations, and satellite-derived winds in addition to AMSU-A satellite radiances. ALEDAS is used to generate an experimental atmospheric global ensemble reanalysis (ALERA).

As already noted, there are two main ALEDAS updates. One is in its model part: the vertical resolution has been increased from 48 vertical levels (up to ∼3 hPa) to 56 (up to ∼0.1 hPa). YN23 describe the benefits of this update. The other is in the data assimilation part, and is particularly significant for ALEDAS’s handling of AMSU-A assimilation: the scheme is similar to that of Terasaki and Miyoshi (2017), which implements the following:

Terasaki et al. (2019) used this scheme to create multiyear global analysis data. The supplemental material (section S1) summarizes the implementation of the scheme in ALEDAS.

Implementing Terasaki and Miyoshi (2017)’s scheme, ALEDAS assimilates AMSU-A channels 6, 7, and 8, which are sensitive to temperature in the midtroposphere to lower stratosphere (Fig. 4 in Terasaki and Miyoshi 2017); the sensors are onboard polar-orbiting satellites NOAA-15, NOAA-18, and NOAA-19, and MetOp-A (Table S1). We follow Terasaki and Miyoshi’s (2017) choice of assimilated channels and satellites for AMSU-A. Sections 3 and 4 show that AMSU-A is expected to have large areas of beneficial OI over the Southern Hemisphere (SH) oceans. This is due to the strict quality control of the AMSU-A observations, which results in channels 6 and 7 being disregarded over land (see section S1.3), and to there being few conventional observations over oceans. Using the same limited AMSU-A channels as Terasaki and Miyoshi (2017) is the first step to assimilating satellite radiance observations in ALEDAS. Improving the use of satellite radiance observations would require further steady developments of previous methods.

ALEDAS can perform adequately in global atmospheric data assimilation systems, and ALERA is somewhat comparable to recent global reanalysis families (Fujiwara et al. 2017, 2022). Figure 1 shows the root-mean-square difference in the 850- and 250-hPa temperature fields, respectively, of ALERA and JRA-55 (Kobayashi et al. 2015). For comparison, the differences of ERA5 (Hersbach et al. 2020) and NCEP-NCAR (Kalnay et al. 1996) reanalyses against JRA-55 are also shown. ALERA departs from JRA-55 by no more than 2 times the difference between ERA5 and JRA-55 in the upper troposphere and 1.5 times in the lower troposphere, and it has a comparable difference between the NCEP–NCAR reanalysis and JRA-55. ALEDAS and ALERA may also be used as a global reanalysis family, because the NCEP–NCAR reanalysis has been used in enormous studies. The 500-hPa geopotential height fields from ALERA and JRA-55 look almost identical over the extratropics (Figs. 1c,f). Moreover, ALERA reproduces two stratospheric sudden warming events in the Northern Hemisphere and the Southern Hemisphere (YN23). Note that similar results were obtained using ERA5 as the reference instead of JRA-55 (see Fig. S4). In light of this, ALERA can faithfully simulate true synoptic–large-scale atmospheric circulations in the troposphere and stratosphere.

Fig. 1.
Fig. 1.

(a),(b),(d),(e) Time series of the globally averaged root-mean-square difference of ALERA (thick black lines), NCEP–NCAR reanalysis (gray lines), and ERA5 (thin black lines) relative to JRA-55 for (top) 300-hPa zonal wind (m s−1) and (middle) 850-hPa temperature (K) for the periods of December 2018–January 2019 in (a) and (b) and August–September 2019 in (d) and (e). (bottom) Snapshots of 500-hPa geopotential height fields in ALERA (red solid lines) and JRA-55 (black dashed lines) at (c) 0000 UTC 21 Dec 2018 and (f) 0000 UTC 23 Aug 2019.

Citation: Weather and Forecasting 38, 6; 10.1175/WAF-D-22-0159.1

b. Definition of accumulated OI

As mentioned in section 1, (E)FSO estimates the OI obtained by noncycling data assimilation. Yamazaki et al. (2021) physically interpreted the OI and the OIEFSO. Here, we define the accumulated OI obtained by cycling data assimilation, which can be compared with OIEFSO.

The OI of all observations in one data assimilation cycle is defined as the total OI (Yamazaki et al. 2021):
ΔE2(ef)TCef(eg)TCeg.
The errors ef and eg are given as follows:
ef=x¯f(t)x¯a(t)andeg=x¯g(t)x¯a(t),
where x(t) indicates the state vectors at time t, and the superscripts a, g, and f indicate those of the analysis, the te-h forecast from the analysis, and the te-h forecast from the background states, respectively; x¯(t) is the ensemble mean of x(t); and C is the matrix defining the norm and global averaging. In ALEDAS the evaluation time is set to te = 6 h, the assimilation window Δt is 6 h, and the referenced state is set to x¯a(t) (Yamazaki et al. 2021). Here, te is a tuning parameter set to 6 h. This is because a comparison of different evaluation times revealed that EFSO at te = 6 and 12 h can best estimate the OI obtained by noncycling OSEs (Yamazaki et al. 2021). In addition, 6-h EFSO has reasonable computational cost (Yamazaki et al. 2021). Further study is needed to elucidate whether it could be useful to use EFSO with longer te to estimate the accumulated OI. Hence, total OI or ΔE2 is equivalent to the difference in the errors between the 12- and 6-h forecasts. Note that EFSO can estimate ΔE2, as total OIEFSO (the sum of OIEFSO by all the observations of one data assimilation cycle) approximates ΔE2 (Yamazaki et al. 2021). The EFSO used here applies the same settings (parameters) as Yamazaki et al. (2021).
Total OI in ALEDAS with and without assimilating AMSU-A observations (CTL and OSE) is defined as ΔECTL2 and ΔEOSE2, respectively. Definitions of ΔECTL2 and ΔEOSE2 are the same as in Eq. (1), but the errors eCTL and eOSE are defined as follows:
eCTLf=x¯CTLf(t)x¯CTLa(t)andeCTLg=x¯CTLg(t)x¯CTLa(t),and
eOSEf=x¯OSEf(t)x¯OSEa(t)andeOSEg=x¯OSEg(t)x¯OSEa(t),
respectively, where x¯CTL(t) and x¯OSE(t) represent the state vectors in the CTL and the OSE, respectively.

Similar to Eq. (1), the OI of the OSE is evaluated by considering the difference between the two forecast errors from the analysis with AMSU-A observations (CTL) and from that without (OSE). Due to the accumulation of OI, cycling AMSU-A data denials separate x¯OSEa(t) from x¯CTLa(t). As x¯OSEa(t) relates to the initial values of x¯OSEg(t) and x¯OSEf(t), these OSE forecasts likely diverge from the corresponding CTL values [x¯CTLg(t) and x¯CTLf(t)]. The OI by the OSE thus depends on x¯CTLa(t),x¯OSEa(t),x¯CTLf(t),x¯OSEf(t),x¯CTLg(t), and x¯OSEg(t). In a noncycling OSE, on the other hand, the forecast from the background state [x¯f(t)] and the analysis as the reference state [x¯a(t)] are the same as those in the CTL, respectively: i.e., x¯OSEf(t)=x¯CTLf(t) and x¯OSEa(t)=x¯CTLa(t). The EFSO formulation starts from Eq. (1) of ΔECTL2, and EFSO can separate the relative contribution of AMSU-A observations, which constitute the total observations in CTL and are excluded from the OSE. Here, the relative contribution of AMSU-A is a constituent part of ΔECTL2. As shown in Eq. (5) of Yamazaki et al. (2021), OIEFSO approximates OI by the noncycling OSE, which is described by x¯CTLa(t),x¯CTLf(t),x¯CTLg(t), and x¯OSEg(t).

We define accumulated OI as the difference between the total OI of CTL and that of the OSE:
IAMSU-AΔECTL2ΔEOSE2.
The term IAMSU-A simply quantifies the contribution of AMSU-A observations by directly comparing the total OI of CTL and the OSE. If an AMSU-A data denial is noncycling [x¯OSEf(t)=x¯CTLf(t)] and x¯OSEa(t) is sufficiently close to x¯CTLa(t) [i.e., x¯OSEa(t)x¯CTLa(t)], then IAMSU-A is equivalent to AMSU-A OIEFSO [Eq. (5) of Yamazaki et al. (2021)]. As AMSU-A OI accumulates through cycling data denials, x¯OSEa(t) will significantly diverge from x¯CTLa(t), and then IAMSU-A could become different from AMSU-A OIEFSO.
The definition of accumulated OI [Eq. (5)] is similar but different from that of Gelaro and Zhu (2009). In the context of Gelaro and Zhu (2009), accumulated OI can be defined as follows:
iAMSU-A(eOSEg)TCeOSEg(eCTLg)TCeCTLg(eCTLg)TCeCTLg.
This definition specifically measures OI regarding forecasts from analyzed states. iAMSU-A can be directly compared with the ratio of AMSU-A OIEFSO to total OIEFSO (Gelaro and Zhu 2009). Therefore, iAMSU-A can measure the contribution rates of the AMSU-A observations. Conversely, IAMSU-A measures OI in a straightforward manner and can be compared with AMSU-A OIEFSO with the dimension of energy. In addition, IAMSU-A includes the information of eCTLf and eOSEf, which are not included in iAMSU-A. Overall, IAMSU-A is particularly specialized for quantifying OI against OIEFSO.
OI is commonly expressed as the energy norm (J kg−1). We describe all OI and errors in terms of moist total energy (Ehrendorfer et al. 1999). The moist total energy can be separated into three parts (Ota et al. 2013; Kotsuki et al. 2019):
eKE=121psps0(u2+υ2)dp,
ePE=12(1psps0cpTrT2dp+RdTrpr2ps2),and
eME=121psps0L2cpTrq2dp,
where p is the vertical pressure coordinate; ps is the surface pressure; u′, υ′, T′, q′, and ps are the errors of the zonal wind, meridional wind, temperature, specific humidity, and surface pressure, respectively, as defined by Eq. (2); cp is the specific heat capacity of air; Rd is the gas constant of dry air; L is the latent heat of condensation per unit mass; and Tr and pr are the reference temperature and surface pressure, respectively. The terms eKE, ePE, and eME, respectively, denote kinetic, potential, and moist energy, and the sum of eKE and ePE is the dry total energy. Hence, we can evaluate to which variables the AMSU-A observations contribute (see section 3e). In addition, we can collect the distributions of IAMSU-A and OIEFSO before operating global averaging in the matrix C (see Yamazaki et al. 2021). Hereafter, all the scalar values show the global averages of vertically integrated values, unless stated otherwise.

c. Experimental design

We conducted two OSEs in which the entire AMSU-A observations were excluded from ALEDAS; i.e., ALEDAS generates ALERA or CTL. The OSEs were separated into two streams: one with a data assimilation cycle from 0000 UTC 21 December 2018 to the end of January 2019 (NH-winter experiment) and the other from 0000 UTC 23 August 2019 to 1800 UTC 3 October 2019 (SH-winter experiment), both initialized with ALERA. To verify our results, the CTL and the OSEs were compared using different data periods. Spinup periods of about 10 days after starting each OSE were discarded; the adjusted periods were from 0000 UTC 1 January to 1800 UTC 31 January 2019 for the NH-winter and from 0000 UTC 3 September to 1800 UTC 3 October 2019 for the SH-winter experiments. As a result, the same conclusions were obtained for the adjusted periods. The supplemental material gives further details (Figs. S6–S13 and Tables S3 and S4).

During the two OSEs, the AMSU-A observations were excluded for each data assimilation cycle; thus, these are the cycling OSEs. We try to evaluate accumulated OI by cycling data denials. Any OI localized in limited areas or altitudes can propagate and accumulate during the cycling data assimilation, affecting the general circulation.

3. Results

a. Comparison of AMSU-A accumulated OI and OIEFSO

We assess the effects of the accumulation of AMSU-A OI. Table 1 lists the OI of the OSEs and EFSO averaged over the NH-winter and the SH-winter experiments. In both OSEs, AMSU-A OIEFSO become positive (i.e., beneficial), indicating the successful implementation of Terasaki and Miyoshi (2017)’s AMSU-A assimilation scheme in ALEDAS. Furthermore, IAMSU-A is negative in both experiments, or ΔECTL2 becomes smaller than ΔEOSE2. This demonstrates that the addition of AMSU-A observations reduces total OI, which is opposite to the expectation from AMSU-A OIEFSO. Note that EFSO estimates OI well in both the CTL and the OSEs, because total OIEFSO performs well in quantifying the total OI in both cases, and their correlation coefficients are very high (>0.9). Therefore, the negative IAMSU-A results from the accumulation of AMSU-A OI.

Table 1

Total OI, total OIEFSO of the CTL and the OSEs, AMSU-A OIEFSO, and IAMSU-A for the NH-winter and SH-winter experiments. All values are averaged for each experimental period. The rightmost column shows the correlation between total OI and total OIEFSO. Units are J kg−1 except for correlation coefficients.

Table 1

The magnitudes (absolute values) of IAMSU-A are greater than those of AMSU-A OIEFSO by factors of 1.2–1.6, possibly due to the effect of the accumulation of AMSU-A OI.

The IAMSU-A and AMSU-A OIEFSO show larger values in the SH-winter experiments than in the NH-winter experiments, suggesting that AMSU-A OI is more significant and/or widespread in Southern Hemisphere winters than summers because of dynamical propagation in the extratropics (as discussed in section 4).

b. Physical interpretation of negative AMSU-A accumulated OI against the positive AMSU-A OIEFSO

We interpret here the reason for IAMSU-A being negative in contrast to the positive beneficial AMSU-A OIEFSO. First considered are the time series of ΔECTL2 and ΔEOSE2 in the NH-winter and SH-winter experiments (Figs. 2a,d): the latter is almost always larger than the former except very near the start of each experiment (21 December 2018 for the NH-winter experiment and 23 August 2019 for the SH-winter experiment). Second, we separate ΔE2 into two parts: (ef)TCef and (eg)TCeg [Eq. (1)]. The former is equivalent to the 12-h forecast error (forecast from the background) and the latter corresponds to the 6-h error (forecast from the analysis). As total OI is defined as the difference between these errors, we examine each error’s time series individually and the reason for ΔECTL2 becoming smaller than ΔEOSE2. From Figs. 2b and 2e, the errors in the CTL generally decrease relative to those in the OSEs. The above findings imply that in the CTL, or by accumulating AMSU-A OI, the short-range forecasts become closer to the analysis than those in the OSEs.

Fig. 2.
Fig. 2.

Time series of (a),(d) total OI in the CTL (black lines) and OSE (blue dashed lines); (b),(e) differences between 6-h (thin lines) and 12-h (thick lines) forecasts and their analyses (J kg−1) in the CTL (black lines) and OSE (blue dashed lines); and (c),(f) differences between the CTL and OSE analyses (J kg−1; black lines) and the analysis ensemble spreads (J kg−1) in the CTL (red upward-pointing triangles) and the OSEs (blue downward-pointing triangles). (left) NH-winter experiment: 21 Dec 2018–31 Jan 2019 and (right) SH-winter experiment: 23 Aug–3 Oct 2019.

Citation: Weather and Forecasting 38, 6; 10.1175/WAF-D-22-0159.1

Last, we consider the effects of the accumulation of AMSU-A OI. Figures 2c and 2f shows time series of the moist-total-energy differences of the analyses between the CTL and the OSEs. The accumulation of OI causes the analyses of the OSEs to depart substantially from the analysis of the CTL. The analyses of these experiments converge on different states after 10–14-day data assimilation cycles of the OSE periods. The differences converged to 7–9 (J kg−1), about 100 times the values of AMSU-A OIEFSO and 2–3 larger than the total OI values (Table 1). As the CTL and OSE forecasts start from different analysis states, the total OIs in the OSEs are larger than that in the CTL, and thus IAMSU-A can be negative against a positive AMSU-A OIEFSO. Interestingly, the disparities between the CTL and OSE analyses have similar values to the OSE analysis ensemble spreads. The analysis ensemble spreads become smaller in the CTL than those in the OSEs.

We consider the reason for the incorporation of AMSU-A observations resulting in a reduction in the total OI or IAMSU-A becoming negative. From Table 1 and Fig. 2, we sketch the effect of the accumulation of AMSU-A OI on a data assimilation system (Fig. 3). The figure shows time sequences of the analysis and background (first guess) fields in the CTL (lines C and C′, respectively) and the OSE (lines O and O′, respectively). Suppose we start cycling data denials at t = t0. Line O will depart from line C, and will reach a different equilibrium state (Figs. 2c,f): the analysis difference between lines O and C overall corresponds to the ensemble spread of the OSE analysis. Comparison of the CTL and the OSE indicates that the total OI in the CTL becomes smaller than that in the OSE (Figs. 2a,d), despite the CTL having more total observations than the OSE. The reason for this result is that the difference between the analysis and the short-range (6–12 h) forecast (i.e., the error) in CTL is smaller than the difference in the OSE (Figs. 2b,e). This further shows that the CTL increment (the difference between lines C′ and C) becomes smaller than the OSE increment (the difference between lines O′ and O). Decreasing the increment in a forecast–analysis cycle brings the short-range forecast closer to its analysis. Therefore, there is less shock in the data assimilation (forecast–analysis) cycle in the CTL, and the addition of AMSU-A observations stabilizes the data assimilation system. We can thus interpret the negative IAMSU-A as being beneficial in a cycling data assimilation system, if we define positive OIEFSO as beneficial [Eq. (1)].

Fig. 3.
Fig. 3.

Schematic representation of the accumulation of OI in a cycling OSE (data-denial experiment) and OI in the CTL. The cycling OSE starts from t0. Lines O and C show the analysis states of the OSE and of CTL, respectively; O′ and C′ are the corresponding backgrounds, or first guesses. Time is plotted on the horizontal axis and the error on the vertical axis. Thin black lines show te-h forecasts from the background state xb(t) and from the analysis xa(t) in the CTL and the OSE, where the data assimilation window is Δt. The difference between the thin lines at t + te is the total OI of the CTL and the OSE: i.e., ΔECTL2 for line C and ΔEOSE2 for line O. The accumulated OI is the difference between ΔECTL2 and ΔEOSE2.

Citation: Weather and Forecasting 38, 6; 10.1175/WAF-D-22-0159.1

c. Changes of OIEFSO in the AMSU-A OSEs

Figures 4a and 4d show the global averaged values and distributions of AMSU-A OIEFSO and those by the other observation types during the NH-winter and SH-winter experiments. Relative rankings of OIEFSO by observation type show the AMSU-A observations are estimated to have the seventh most beneficial impact among the nine types in both experiments. The AMSU-A OIEFSO distributions in Fig. 5 indicate that the AMSU-A observations have large beneficial impacts in the middle and upper troposphere in the Southern Hemisphere midlatitudes. As only channels 6–8 are assimilated, the AMSU-A OIEFSO distributions are confined to the upper troposphere. However, the distributions show that the Southern Hemisphere storm track areas benefit from the AMSU-A observations. Therefore, as storm tracks are important driving forces of the general circulation (Edmon et al. 1980; Thorncroft et al. 1993), the observations may have a generally positive influence. Regarding the other observation types, those with the three largest OIEFSO are satellite winds (SATWND), ocean scatterometer (ASCATW), and in situ land (ADPSFC) in order. OIEFSO overlapping with the large AMSU-A OIEFSO distributions are radiosondes (APDUPA), aircraft (AIRCFT), and satellite winds (SATWND).

Fig. 4.
Fig. 4.

Global averages of OIEFSO for each observation type per data assimilation cycle averaged for the periods of (left) the NH-winter experiment and (right) the SH-winter experiment. (a),(d) OIEFSO (J kg−1); (b),(e) differences between OIEFSOCTL and (minus) OIEFSOOSE (J kg−1); and (c),(f) ratios of OIEFSOOSE to OIEFSOCTL (%) for each observation type.

Citation: Weather and Forecasting 38, 6; 10.1175/WAF-D-22-0159.1

Fig. 5.
Fig. 5.

AMSU-A OIEFSO distributions for each observation type per data assimilation cycle averaged for the periods of (a) the NH-winter experiment and (b) the SH-winter experiment. Latitude–height cross sections of the AMSU-A OIEFSO distributions (J kg−1) (large plots on the left) and the OIEFSO distributions (J kg−1) of the other observation types (small plots).

Citation: Weather and Forecasting 38, 6; 10.1175/WAF-D-22-0159.1

We anticipate that the addition of AMSU-A observations will alter the relative contributions of each observation type’s OI to the total OI. Comparison of OIEFSO in the CTL and in the OSEs shows which observation types have their OI reduced by the addition of AMSU-A observations. Here, we distinguish the OIEFSO derived from the CTL as OIEFSOCTL and that from the OSEs as OIEFSOOSE. Gelaro and Zhu (2009) previously used the same approach with FSO for different observation types.

Figures 4b and 4e show the differences of OIEFSOCTL and OIEFSOOSE. The AMSU-A observations contribute to the total OIEFSOCTL by reducing the OIEFSOCTL of SATWND, ASCATW, ADPSFC, ADPUPA, AIRCAR (aircraft observations other than AIRCFT), and AIRCFT. These six observation types have the six greatest OIEFSOCTL (Figs. 4a,d). Moreover, in the OSEs, SATWND and ASCATW observations show the greatest reductions of their OIEFSOCTL. These two observation types have the two largest OIEFSOCTL (Figs. 4a,d). Figures 4c and 4f shows the ratios of OIEFSOOSE to OIEFSOCTL for each observation type. Similar to the finding for the differences, the same six observation types show the greatest reductions of OIEFSOCTL relative to OIEFSOOSE. Nonetheless, the OIEFSOCTL of ADPUPA, AIRCAR, AIRCFT, and SATWND decrease compared with the OIEFSOCTL of ASCATW and ADPSFC; the former types are observed above the upper troposphere, whereas the latter are measured near the surface (Fig. 5). Among the observation types of ADPUPA, AIRCAR, AIRCFT, and SATWND, the ratios of AIRCAR are relatively small (Figs. 4c,f).

Figure 6 shows the distributions of the differences between OIEFSOCTL and OIEFSOOSE for six observation types: ADPUPA, AIRCAR, AIRCFT, SATWND, ADPSFC, and ASCATW. As anticipated from Fig. 5, the OIEFSOCTL values (except for AIRCAR) generally decrease over the Southern Hemisphere midlatitudes (30°–60°S), which is the latitudinal band of peak beneficial AMSU-A OIEFSO. The AMSU-A observations reduce the OIEFSOCTL of these observation types. Regions of distinctly reduced OIEFSOCTL occur where both AMSU-A OIEFSO and those OIEFSOCTL have substantial beneficial values; i.e., from the lower troposphere to the upper stratosphere. AIRCAR OIEFSOCTL has beneficial peaks in the upper troposphere, similar to AMSU-A OIEFSO, although AIRCAR primarily observes the Northern Hemisphere midlatitudes. Overall, the OIs of various observation types that directly measure the Southern Hemisphere midlatitudes and/or the upper troposphere are largely substituted by AMSU-A OI.

Fig. 6.
Fig. 6.

As in Fig. 5, but with shading showing the difference between OIEFSOCTL and OIEFSOOSE (J kg−1) for the (a) NH-winter and (b) SH-winter experiments. Only results for observation types ADPUPA, AIRCAR, AIRCFT, SATWND, ADPSFC, and ASCATW are shown. Contours indicate the OIEFSOCTL (J kg−1) of each observation type.

Citation: Weather and Forecasting 38, 6; 10.1175/WAF-D-22-0159.1

It is difficult to explain the tendency of AIRCAR observations to reduce OIEFSOCTL, as the beneficial OI peaks of AIRCAR and AMSU-A observations show little overlap with each other. This may be related to the AIRCAR observed variables including temperature. Further research should assess which variables’ observations tend to be largely substituted by AMSU-A observations of brightness temperature.

d. Distributions of AMSU-A accumulated OI and OIEFSO

We now assess the spatial correspondence between beneficial IAMSU-A and AMSU-A OIEFSO. Figure 7 shows latitude–height sections of zonal mean IAMSU-A. We find distinct beneficial IAMSU-A signals in the upper troposphere at Southern Hemisphere midlatitudes in the NH-winter and SH-winter experiments. These beneficial signals largely overlap with those of beneficial AMSU-A OIEFSO in Fig. 5: the amplitudes (absolute values) of the IAMSU-A and AMSU-A OIEFSO signals are overall of the same order. The results indicate that EFSO is helpful in estimating beneficial IAMSU-A distributions. However, the IAMSU-A signals are distributed more broadly than the AMSU-A OIEFSO signals, likely because of the propagation of accumulated AMSU-A OI. Beneficial IAMSU-A is widely distributed in the troposphere and the Southern Hemisphere stratosphere: the accumulation of AMSU-A OI brings beneficial impacts globally. These findings are consistent with Table 1.

Fig. 7.
Fig. 7.

IAMSU-A (J kg−1; shading) averaged during the (a) NH-winter (21 Dec 2018–31 Jan 2019) and (b) SH-winter (23 Aug–3 Oct 2019) experiments. Latitude–height cross sections of zonal mean values are shown. Contours indicate the zonal mean wind fields (m s−1) of the CTL analysis averaged over each experimental period.

Citation: Weather and Forecasting 38, 6; 10.1175/WAF-D-22-0159.1

Similar to Fig. 7, Fig. 8 compares horizontal distributions of AMSU-A OIEFSO and IAMSU-A. Beneficial AMSU-A OIEFSO is distributed mainly over the southeastern Pacific and the Southern Ocean, and is weakly apparent in the northeastern Pacific (Figs. 8a,c). These areas seem particularly beneficial because conventional observations are sparse, but AMSU-A measurements from polar-orbit satellites are distributed uniformly worldwide. In particular, we find distinct beneficial AMSU-A OIEFSO signals over the southeastern Pacific off the coast of Peru, and Figs. 8b and 8d show the IAMSU-A signals located just downstream of the midlatitude area of peak AMSU-A OIEFSO in both experiments. These horizontal maps again suggest the propagation of accumulated AMSU-A OI in the Southern Hemisphere midlatitudes.

Fig. 8.
Fig. 8.

Horizontal distributions of vertically integrated (a),(c) AMSU-A OIEFSO and (b),(d) IAMSU-A (J kg−1) averaged for the periods of the (left) NH-winter and (right) SH-winter experiments.

Citation: Weather and Forecasting 38, 6; 10.1175/WAF-D-22-0159.1

In the Southern Hemisphere stratosphere, beneficial IAMSU-A appears in both OSEs despite the estimation that AMSU-A observations do not directly affect these regions (Fig. 5). The stratospheric IAMSU-A is attributed to the effects of the accumulation of AMSU-A OI. Small and broad beneficial IAMSU-A in the NH-winter and SH-winter experiments would stem from the propagation of accumulated OI via radiation processes, as the distribution of beneficial IAMSU-A is independent of the zonal wind fields (Fig. 7). A distinct beneficial IAMSU-A signal appears in the polar Southern Hemisphere stratosphere below 10 hPa in the SH-winter experiment, because this experiment’s AMSU-A OI accumulation is greater than that in the NH-winter experiment (Table 1 and Fig. 2). The stratospheric signal is attributed to the propagation of accumulated AMSU-A OI through dynamical processes. This is explained in section 4.

Finally, we check whether the beneficial IAMSU-A signals actually involve regions where the CTL analysis fields are improved relative those in the OSEs. Figure 9 shows the analysis fields’ root-mean-square difference distributions against JRA-55. This figure shows the difference of root-mean-square difference of the CTL minus the OSE analysis fields in moist total energy (J kg−1): more negative values indicate the CTL analysis fields to be superior to the OSE analysis fields. In both the NH-winter and SH-winter experiments, the improved (negative root-mean-square difference) regions show a good correspondence to the beneficial IAMSU-A signals in Fig. 7. Here, similar results were obtained using ERA5 as for the reference fields (Fig. S5). Therefore, regions of beneficial IAMSU-A involve improving the analysis fields. Regions of beneficial IAMSU-A signals also emerge where AMSU-A observations stabilize the forecast–analysis (i.e., data assimilation) cycles. EFSO can estimate such regions overall if we consider the OI distributions.

Fig. 9.
Fig. 9.

As in Fig. 7, but for the difference of the root-mean-square differences between the CTL analysis and the OSE analysis against JRA-55 (J kg−1); i.e., the error of CTL referenced against JRA-55 minus that of OSE. Negative values indicate that the CTL analysis is closer to JRA-55 than the OSE.

Citation: Weather and Forecasting 38, 6; 10.1175/WAF-D-22-0159.1

e. Relative contributions of AMSU-A OI to kinetic, potential, and moist energy

As OI has been normalized by the moist total energy, we can decompose the OI of AMSU-A OIEFSO or IAMSU-A to the kinetic (eKE), potential (ePE), and moist (eME) energy components [Eqs. (7)(9)]. To decompose and evaluate the relative contributions of those energy components, we can discuss which variables AMSU-A observations will contribute in ALEDAS.

Table 2 lists the relative contributions of eKE, ePE, and eME to the total OI of CTL (ΔECTL2) and total OIEFSOCTL, AMSU-A OIEFSO, and IAMSU-A. Initially, total observations at each cycle (for ΔECTL2) show contributions to the moist total energy from eKE of more than 60%, from eME of ∼23%, and from ePE of ∼13%. EFSO (for total OIEFSOCTL) performs well in estimating the energy components of ΔECTL2. Compared with ΔECTL2 or total OIEFSOCTL, the AMSU-A observations estimated by AMSU-A OIEFSO have contributions from ePE of 20%–30%, which are about double those of ePE for the total observations; in contrast, the contributions from eME are at most ∼5%, much smaller than those in ΔECTL2. The relative contributions of AMSU-A accumulated OI (IAMSU-A) show similar ratios to those in AMSU-A OIEFSO, except that the contributions to IAMSU-A from eKE and ePE become larger and smaller than those in AMSU-A OIEFSO, respectively, in the SH-winter experiment.

Table 2

Contributions (%) of kinetic (eKE), potential (ePE), and moist (eME) energy to the moist total energy (=eKE + ePE + eME) in the global averages of total OI of CTL ΔECTL2, total OIEFSOCTL (in parentheses), AMSU-A OIEFSO, and IAMSU-A for the NH-winter and SH-winter experiments. All values are averaged for the period of each experiment.

Table 2

Therefore, AMSU-A observations tend to improve the variables related to potential energy (i.e., temperature and surface pressure) among the global (total) observations. However, as the AMSU-A observations have a greater impact on the Southern Hemisphere midlatitudes than anywhere else (section 3d), kinetic (i.e., wind) and potential energy of IAMSU-A are transferred interchangeably owing to extratropical geostrophic adjustment through the accumulation (cycling data assimilation) periods, which would be more active during Southern Hemisphere winters.

Figure 10 shows AMSU-A OIEFSO and IAMSU-A distributions decomposed into their eKE, ePE, and eME components. The AMSU-A OIEFSO distributions (Figs. 10a,b) show large eKE and ePE contributions in the middle and upper troposphere in the Southern Hemisphere midlatitudes during both experiments. Their eKE and ePE distributions show similar patterns to each other, although the eKE distributions are rather broad toward the tropics. In contrast, the eME distributions of AMSU-A OIEFSO are small and mostly confined to the tropical middle troposphere, far from the major beneficial AMSU-A OIEFSO region in the Southern Hemisphere midlatitudes (Figs. 5 and 8a,c).

Fig. 10.
Fig. 10.

Latitude–height cross sections of (top) eKE, (middle) ePE, and (bottom) eME of (a),(b) AMSU-A OIEFSO and (c),(d) IAMSU-A distributions (J kg−1) averaged for the periods of the NH-winter experiment in (a) and (c) and the SH-winter experiment in (b) and (d).

Citation: Weather and Forecasting 38, 6; 10.1175/WAF-D-22-0159.1

For the IAMSU-A distributions (Figs. 10c,d), the eKE, ePE, and eME contributions show similar but broader patterns to those of the AMSU-A OIEFSO. As the eKE and ePE of the AMSU-A OIEFSO are largely distributed in the Southern Hemisphere midlatitudes, those energy components of beneficial IAMSU-A are distributed in the Southern Hemisphere extratropical stratosphere, possibly because of the propagation of accumulated OI (section 4).

In summary, AMSU-A observations have a relatively large impact on the potential energy of OI, which includes temperature improvements, compared with all the observation types (total observations). This is consistent with the AMSU-A channels used here being sensitive to atmospheric temperature. For accumulated AMSU-A OI, the potential energy and kinetic energy are exchanged with each other in the Southern Hemisphere midlatitudes through the cycling OSEs; much potential energy was transferred to kinetic energy during the SH-winter experiment. In contrast, AMSU-A observations contribute less to the moist energy of OI, because this energy is mainly distributed in the tropics, far from the Southern Hemisphere midlatitudes.

4. Discussion: Propagation of accumulated OI

Section 3 shows that the accumulation of AMSU-A OI causes IAMSU-A to become generally negative (beneficial), particularly over the upper troposphere of the Southern Hemisphere midlatitudes and over the Southern Hemisphere polar stratosphere in the Southern Hemisphere winter (Fig. 7). The distributions of AMSU-A OIEFSO overall overlap with those of beneficial IAMSU-A, but the latter distributions are broader (Figs. 5, 7, 9, and 10). This result is associated with the dynamical propagation of accumulated OI. Consequently, this section contrasts the distributions of the tropospheric AMSU-A OIEFSO and IAMSU-A signals in the upper troposphere. We will discuss the cause of the stratospheric IAMSU-A signal in the SH-winter experiment.

Figure 11 shows horizontal maps of OIEFSO and IAMSU-A at 300 hPa in regions of the Southern Hemisphere midlatitudes with notable tropospheric IAMSU-A and AMSU-A OIEFSO signals (Figs. 5, 7, 8, and 9). The AMSU-A OIEFSO distributions indicate the peak beneficial area over the southeastern Pacific off the coast of Peru in both the NH-winter and SH-winter experiments (Figs. 11a,d). The local sparsity of global observations and the Southern Hemisphere midlatitude storm tracks both contribute to the concentration of the AMSU-A OIEFSO signal in the peak area. The former is apparent from Figs. 11b and 11e. These panels show the sum of OIEFSOCTL for observation types other than AMSU-A. In greater areas of the sum of OIEFSOCTL, more observation information is obtained; i.e., the areas correspond to regions where observations are sufficiently dense even if AMSU-A observations are not available, and do not overlap with the peak area. In addition, the peak area coincides with the location of concentrated Southern Hemisphere storm-track activity (Figs. 11a,d), which is visualized by absolute values of meridional wind (e.g., Nishii and Nakamura 2010).

Fig. 11.
Fig. 11.

Shading shows 300-hPa (a),(d) AMSU-A OIEFSO; (b),(e) OIEFSO by observation types other than AMSU-A; and (c),(f) IAMSU-A (J kg−1) averaged during the NH-winter [(left) 21 Dec 2018–31 Jan 2019] and the SH-winter [(right) 23 Aug–3 Oct 2019] experiments. Contours in (a),(d) and (c),(f) indicate the absolute values of meridional wind (m s−1) and wind magnitude (m s−1) in CTL at 300 hPa, respectively, averaged over each experimental period.

Citation: Weather and Forecasting 38, 6; 10.1175/WAF-D-22-0159.1

In contrast, the beneficial IAMSU-A signals extend just downstream of the beneficial peak area of AMSU-A OIEFSO: from South America via the South Atlantic to the south Indian Ocean in both experiments. As EFSO estimates OI well in the short-range forecast period (Yamazaki et al. 2021), we speculate that AMSU-A OI tends to accumulate downstream via the cycling data assimilation, which includes short-range (6-h) forecast cycles. The beneficial IAMSU-A distributions are concentrated along the westerly jet (Figs. 11c,f), so that the accumulation and upper tropospheric dynamical propagation are related via advection and the movement of Rossby waves along the jet. Such upper-tropospheric dynamical propagation has been found in a cycling OSE (Yamazaki et al. 2015) and in short- to medium-range forecast experiments (e.g., Yamazaki et al. 2021; Sato et al. 2020). The beneficial IAMSU-A distribution in the SH-winter experiment is somewhat broader than that in the NH-winter experiment (Figs. 11c,f). This seems to be associated with the double jet structure in the Southern Hemisphere winter (Fig. 11f). In fact, Fig. 7b shows two peaks in the tropospheric IAMSU-A signal over the latitudinal bands of the tropospheric jet around 30°S and of the stratospheric jet around 60°S.

We next consider the beneficial IAMSU-A signal in the Southern Hemisphere polar stratosphere in the SH-winter experiment. In terms of the general circulation, the distinct beneficial IAMSU-A signal mostly appears at the dominant upward Eliassen–Palm flux region from the upper troposphere to the stratosphere (McIntyre 1982; Andrews et al. 1987). During the SH-winter experiment, a stratospheric sudden warming occurred in the Southern Hemisphere around 7 September (e.g., Noguchi et al. 2020; Lim et al. 2021). During the warming event, the Eliassen–Palm fluxes became poleward (YN23), which may be associated with the stratospheric IAMSU-A signal over the Southern Hemisphere polar region.

Figure 12 shows time–height sections of the area-averaged IAMSU-A over the Southern Hemisphere polar cap area (60°–90°S) at 300 hPa (upper troposphere) to 30 hPa (middle stratosphere). A peak beneficial IAMSU-A signal occurs in the upper troposphere (300 hPa) on 1–5 September. This tropospheric peak signal (1–5 September) propagates upward to the stratosphere (6–10 September). Figures 13a and 13b show longitude–height cross sections for 1–5 and 6–10 September. The tropospheric beneficial IAMSU-A signal around 90°W extends toward the stratosphere during 1–10 September. The geopotential height field tilted westward during this period, causing planetary waves to propagate higher. On 1–5 September, a blocking anticyclone formed around 90°W in the upper troposphere (Fig. 13A), demonstrating that this blocking was the source of the propagation of the tropical OI signal toward the stratosphere. The above dynamical processes thus contribute to both the tropospheric and stratospheric IAMSU-A signals in Fig. 7b.

Fig. 12.
Fig. 12.

Time series of area-averaged, 5-day running mean IAMSU-A (J kg−1) over the 60°–90°S polar cap area at (from top to bottom) 30, 50, 70, 100, and 300 hPa. Letters indicate the corresponding periods in Fig. 13.

Citation: Weather and Forecasting 38, 6; 10.1175/WAF-D-22-0159.1

Fig. 13.
Fig. 13.

(a)–(d) Longitude–height cross sections of IAMSU-A (J kg−1; shading) and geopotential height (m) in the CTL (contours) averaged over the (a),(b) 40°–70°S and (c),(d) 60°–90°S meridional bands for the SH-winter experiment. The 5-day mean fields are shown for 1–5 Sep 2019 in (a), 6–10 Sep 2019 in (b), 21–25 Sep 2019 in (c), and 26–30 Sep 2019 in (d). The contours show the 5-day mean geopotential height deviations (m) from their zonal averages, averaged over the meridional band. (A) As in (a), but for 300 hPa, the 5-day mean IAMSU-A (J kg−1; shading), and 300-hPa geopotential height fields (m) in the CTL (contours).

Citation: Weather and Forecasting 38, 6; 10.1175/WAF-D-22-0159.1

In addition to the upward propagation (Fig. 12), there is a beneficial IAMSU-A signal in the stratosphere (above 100 hPa) with an equivalent barotropic structure on 21–25 September. The longitude–height cross section of that stratospheric signal on 21–25 September (Fig. 13c) shows the beneficial IAMSU-A signal in the stratosphere at about 60°W, which is also equivalent to the barotropic structure from 30 to 100 hPa. During 21–25 September, the geopotential height field also has an equivalent barotropic or slightly eastward-tilted structure, as is often seen during stratospheric sudden warming, with planetary waves from the troposphere being absorbed or reflected (Kodera et al. 2008; Mukougawa et al. 2017). On 26–30 September (Fig. 13d), the stratospheric signal seems to propagate slightly lower, around 100 hPa, and eastward (downstream) around 120°E. With the stratospheric sudden warming, the stratospheric OI signal appears to build up and intensify in the stratosphere. It may even transmit to some degree downward toward the upper troposphere. The above stratospheric interior processes contribute to at least the stratospheric signal in Fig. 7b.

In summary, beneficial AMSU-A OIEFSO signals appear in regions of low observation density and active storm tracks in the SH, whereas beneficial IAMSU-A signals tend to accumulate in regions dynamically downstream of the beneficial OIEFSO regions. In the SH-winter experiment, AMSU-A OI was input originally in the upper troposphere, propagated along the upper tropospheric westerly jet, and was amplified by troposphere–stratosphere coupling processes, possibly due to both blocking and a stratospheric sudden warming, which could contribute to the large accumulated AMSU-A OI in the Southern Hemisphere winter hemisphere. The findings in this section provide evidence that (E)FSO can effectively evaluate the distribution of accumulated OI by considering dynamical propagation.

5. Conclusions

Using the updated AFES-LETKF data assimilation system (ALEDAS), we evaluate accumulated AMSU-A OI in cycling (repeating) OSEs. The updated ALEDAS implements the AMSU-A satellite radiance assimilation scheme initially developed by Terasaki and Miyoshi (2017). In the cycling situation a subset of observations or a particular observation type is continually denied or added, as commonly occurs during an observation campaign that continues for a certain period or incorporates a new type of satellite observation. We used the ensemble-based forecast sensitivity to observations (EFSO) technique to estimate the accumulated OI. The accumulated OI due to AMSU-A and the OI estimated by EFSO (OIEFSO) were compared through AMSU-A observing system (data denial) experiments.

First, we assessed the reproducibility of ALEDAS, finding that it enables realistic reproduction of tropospheric circulations comparable to the NCEP–NCAR reanalysis in recent global reanalysis families. Second, we defined accumulated OI slightly differently from previous studies. Our definition is based on total OI by all the observations of each data assimilation cycle. The accumulated OI of AMSU-A is defined as the difference in total OI with (CTL) and without (OSE) AMSU-A observations. Last, two cycling AMSU-A OSEs were performed for a Northern Hemisphere and a Southern Hemisphere winter to assess the accumulated AMSU-A OI. We then tested whether the accumulated OI can be estimated by OIEFSO.

AMSU-A OI had negative values that were almost of the same order as the positive (i.e., beneficial) values of AMSU-A OIEFSO; the accumulation effects caused the total OI in the CTL to become smaller than that in the OSEs. This reduction that resulted from including AMSU-A observations was explained by the CTL’s short-range forecast errors being lower than those of the OSEs. The OSE analysis states diverged from the CTL analysis states and converged to their equilibrium states within about 2 weeks; the order of their differences was nearly identical to the spread of the OSE analysis ensemble. The results also indicate that AMSU-A assimilation stabilized ALEDAS via the accumulation of AMSU-A OI. Comparisons of OIEFSO in the CTL and the OSEs show that AMSU-A OIEFSO primarily supports OIEFSO by other observation types that directly measure the Southern Hemisphere midlatitudes and/or near the upper troposphere in the global observing system.

The negative accumulated OI of AMSU-A was greatest in the Southern Hemisphere midlatitudes in both the NH-winter and SH-winter experiments, particularly in the upper troposphere where westerly jets or active storm tracks were present and observations other than AMSU-A were scarce. In the same region, the CTL analysis fields were closer to JRA-55 than the OSE fields; again, the negative accumulated OI was beneficial. The estimated AMSU-A OIEFSO showed its most beneficial peak in the same region, implying that EFSO is helpful in estimating the distributions of beneficial accumulated OI. Beneficial accumulated AMSU-A OI had broader distributions, and thus larger absolute values (amplitudes), than AMSU-A OIEFSO, possibly due to the accumulation during the cycling data assimilation. We then discussed the horizontal distribution and time evolution of the accumulated OI. In the upper troposphere of Southern Hemisphere midlatitudes, AMSU-A OI tended to accumulate downstream of the southeastern Pacific area where the OIEFSO peak signals were estimated. Accumulated AMSU-A OI had a second beneficial peak in the Southern Hemisphere polar stratosphere during the SH-winter experiment. The stratospheric peak might be attributable to the upward planetary wave propagation of a beneficial OI signal that originated from blocking in the upper troposphere during the cycling data assimilation. A stratospheric sudden warming event in the Southern Hemisphere may have contributed a stratospheric beneficial OI signal with an equivalent barotropic structure in the mid- to lower stratosphere. Therefore, stratosphere–troposphere coupling processes could contribute to amplifying the accumulation of AMSU-A OI in the Southern Hemisphere winter. In conclusion, if we consider the dynamical propagation of OI through the cycling data assimilation, both the beneficial OI peaks are connected to dynamical propagation, and EFSO can help estimate the distributions of beneficial accumulated AMSU-A OI.

Considering the development of data assimilation systems, the results obtained here indicate that Terasaki and Miyoshi’s (2017) scheme of assimilating AMSU-A satellite radiance data has been successfully implemented in ALEDAS, a global ensemble Kalman filter–based system. This development is significant for ALEDAS, as it previously could assimilate only conventional, aircraft, and satellite-wind observations.

As discussed in several previous studies, designing actual observing systems requires assessment of the differences between (E)FSO estimations and OSE outputs (e.g., Gelaro and Zhu 2009; Lawrence et al. 2019; Candy et al. 2021; Eyre 2021). As ALEDAS can moderately reproduce the actual general circulation, the results obtained here could be helpful in the above research streams. This study found that EFSO is useful for estimating distributions of the analysis improvement due to the accumulated effects of AMSU-A OI (section 3d). In the next step, we will investigate whether EFSO is useful for estimating the impact of accumulated AMSU-A OI on the forecast fields.

Assimilating AMSU-A satellite radiances is an essential part for our OSE studies that have been done using the previous ALEDAS. The latest ALEDAS is anticipated to outperform the previous one by successfully assimilating AMSU-A observations. Furthermore, the current ALEDAS could be used to assimilate satellite radiance observations that observation campaigns might specifically obtain. Further updates to ALEDAS to make it more useful include using Terasaki and Miyoshi (2017)’s scheme to assimilate more AMSU-A channels that are sensitive to the stratosphere (Miyoshi and Sato 2007; Miyoshi et al. 2010) and incorporating other satellite radiance observations. Exchanging several model schemes that have recently been successfully implemented into AFES (e.g., Baba 2020; Baba and Ogata 2022) could also further improve the performance of ALEDAS.

1

Regarding a benefit of the former update, we report elsewhere that the data assimilation system can efficiently reproduce stratospheric sudden warming events (YN23).

2

The abbreviation of Radiative Transfer for Television and Infrared Observation Satellite (TIROS) Operational Vertical Sounder (TOVS).

3

The abbreviation of “Prepared” Binary Universal Form for the Representation of meteorological data.

Acknowledgments.

We appreciate the constructive comments of two anonymous reviewers. The Earth Simulator at JAMSTEC was used for ALEDAS and to create ALERA. Dr. Yuya Baba has implemented the non-orographic gravity wave drag scheme in AFES (Baba 2019). This work was supported by JSPS KAKENHI (20H01976, 18K13617, and 19H05702) and the Arctic Challenge for Sustainability II (ArCS II).

Data availability statement.

Datasets of the analysis and background ensemble mean and spread fields of ALERA are available from https://www.jamstec.go.jp/esc/fes/dods/alera3/plev. NCEP PrepBUFR and AMSU-A brightness temperature data are also available online (https://doi.org/10.5065/Z83F-N512 and https://doi.org/10.5065/DWYZ-Q852).

REFERENCES

  • Andrews, D. G., J. R. Holton, and C. B. Leovy, 1987: Middle Atmosphere Dynamics. Academic Press, 489 pp.

  • Baba, Y., 2019: Spectral cumulus parameterization based on cloud-resolving model. Climate Dyn., 52, 309334, https://doi.org/10.1007/s00382-018-4137-z.

    • Search Google Scholar
    • Export Citation
  • Baba, Y., 2020: Shallow convective closure in a spectral cumulus parameterization. Atmos. Res., 233, 104707, https://doi.org/10.1016/j.atmosres.2019.104707.

    • Search Google Scholar
    • Export Citation
  • Baba, Y., and T. Ogata, 2022: Resolution dependence of tropical cyclones simulated by a spectral cumulus parameterization. Dyn. Atmos. Ocean, 97, 101283, https://doi.org/10.1016/j.dynatmoce.2022.101283.

    • Search Google Scholar
    • Export Citation
  • Candy, B., J. Cotton, and J. Eyre, 2021: Recent results of observation data denial experiments. Met Office Weather Science Tech. Rep. 641, 35 pp., https://www.metoffice.gov.uk/binaries/content/assets/metofficegovuk/pdf/research/weather-science/frtr_641_2021p.pdf.

  • Chen, T.-C., and E. Kalnay, 2019: Proactive quality control: Observing system simulation experiments with the Lorenz’96 model. Mon. Wea. Rev., 147, 5367, https://doi.org/10.1175/MWR-D-18-0138.1.

    • Search Google Scholar
    • Export Citation
  • Chen, T.-C., and E. Kalnay, 2020: Proactive quality control: Observing system experiments using the NCEP Global Forecast System. Mon. Wea. Rev., 148, 39113931, https://doi.org/10.1175/MWR-D-20-0001.1.

    • Search Google Scholar
    • Export Citation
  • Edmon, H. J., Jr., B. J. Hoskins, and M. E. McIntyre, 1980: Eliassen–Palm cross-sections for the troposphere. J. Atmos. Sci., 37, 26002616, https://doi.org/10.1175/1520-0469(1980)037<2600:EPCSFT>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Ehrendorfer, M., R. M. Errico, and K. D. Raeder, 1999: Singular-vector perturbation growth in a primitive equation model with moist physics. J. Atmos. Sci., 56, 16271648, https://doi.org/10.1175/1520-0469(1999)056<1627:SVPGIA>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Emanuel, K. A., 1991: A scheme for representing cumulus convection in large-scale models. J. Atmos. Sci., 48, 23132329, https://doi.org/10.1175/1520-0469(1991)048<2313:ASFRCC>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Enomoto, T., A. Kuwano-Yoshida, N. Komori, and W. Ohfuchi, 2008: Description of AFES 2: Improvements for high-resolution and coupled simulations. High Resolution Numerical Modelling of the Atmosphere and Ocean, K. Hamilton and W. Ohfuchi, Eds., Springer, 77–97, https://doi.org/10.1007/978-0-387-49791-4_5.

  • Enomoto, T., T. Miyoshi, Q. Moteki, J. Inoue, M. Hattori, A. Kuwano-Yoshida, N. Komori, and S. Yamane, 2013: Observing-system research and ensemble data assimilation at JAMSTEC. Data Assimilation for Atmospheric, Oceanic and Hydrologic Applications (Vol. II), S. Park and L. Xu, Eds., Springer, 509–526, https://doi.org/10.1007/978-3-642-35088-7_21.

  • Eyre, J. R., 2021: Observation impact metrics in NWP: A theoretical study. Part I: Optimal systems. Quart. J. Roy. Meteor. Soc., 147, 31803200, https://doi.org/10.1002/qj.4123.

    • Search Google Scholar
    • Export Citation
  • Fujiwara, M., and Coauthors, 2017: Introduction to the SPARC Reanalysis Intercomparison Project (S-RIP) and overview of the reanalysis systems. Atmos. Chem. Phys., 17, 14171452, https://doi.org/10.5194/acp-17-1417-2017.

    • Search Google Scholar
    • Export Citation
  • Fujiwara, M., G. L. Manney, L. J. Gray, and J. S. Wright, Eds., 2022: SPARC Reanalysis Intercomparison Project (S-RIP) final report. SPARC Rep. 10, 635 pp., WCRP-17/2020, https://doi.org/10.17874/800dee57d13.

  • Gelaro, R., and Y. Zhu, 2009: Examination of observation impacts derived from observing system experiments (OSEs) and adjoint models. Tellus, 61A, 179193, https://doi.org/10.1111/j.1600-0870.2008.00388.x.

    • Search Google Scholar
    • Export Citation
  • Gelaro, R., R. H. Langland, S. Pellerin, and R. Todling, 2010: The THORPEX observation impact intercomparison experiment. Mon. Wea. Rev., 138, 40094025, https://doi.org/10.1175/2010MWR3393.1.

    • Search Google Scholar
    • Export Citation
  • Hattori, M., A. Yamazaki, S.-Y. Ogino, P. Wu, and J. Matsumoto, 2017: Impact of the radiosonde observations of cold surge over the Philippine Sea on the tropical region and the Southern Hemisphere in December 2012. SOLA, 13, 1924, https://doi.org/10.2151/sola.2017-004.

    • Search Google Scholar
    • Export Citation
  • Hersbach, H., and Coauthors, 2020: The ERA5 global reanalysis. Quart. J. Roy. Meteor. Soc., 146, 19992049, https://doi.org/10.1002/qj.3803.

    • Search Google Scholar
    • Export Citation
  • Hotta, D., T.-C. Chen, E. Kalnay, Y. Ota, and T. Miyoshi, 2017: Proactive QC: A fully flow-dependent quality control scheme based on EFSO. Mon. Wea. Rev., 145, 33313354, https://doi.org/10.1175/MWR-D-16-0290.1.

    • Search Google Scholar
    • Export Citation
  • Houtekamer, P. L., and F. Zhang, 2016: Review of the ensemble Kalman filter for atmospheric data assimilation. Mon. Wea. Rev., 144, 44894532, https://doi.org/10.1175/MWR-D-15-0440.1.

    • Search Google Scholar
    • Export Citation
  • Inoue, J., A. Yamazaki, J. Ono, K. Dethloff, M. Maturilli, R. Neuber, P. Edwards, and H. Yamaguchi, 2015: Additional Arctic observations improve weather and sea-ice forecasts for the northern sea route. Sci. Rep., 5, 16868, https://doi.org/10.1038/srep16868.

    • Search Google Scholar
    • Export Citation
  • Kalnay, E., and Coauthors, 1996: The NCEP/NCAR 40-Year Reanalysis Project. Bull. Amer. Meteor. Soc., 77, 437472, https://doi.org/10.1175/1520-0477(1996)077<0437:TNYRP>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Kalnay, E., Y. Ota, T. Miyoshi, and J. Liu, 2012: A simpler formulation of forecast sensitivity to observations: Application to ensemble Kalman filters. Tellus, 64A, 18462, https://doi.org/10.3402/tellusa.v64i0.18462.

    • Search Google Scholar
    • Export Citation
  • Kobayashi, S., and Coauthors, 2015: The JRA-55 reanalysis: General specifications and basic characteristics. J. Meteor. Soc. Japan, 93, 548, https://doi.org/10.2151/jmsj.2015-001.

    • Search Google Scholar
    • Export Citation
  • Kodama, C., A. Kuwano-Yoshida, S. Watanabe, T. Doi, H. Kashimura, and T. Nasuno, 2019: JAMSTEC model intercomparision project (JMIP). JAMSTEC Rep. Res. Dev., 28, 534, https://doi.org/10.5918/jamstecr.28.5.

    • Search Google Scholar
    • Export Citation
  • Kodera, K., H. Mukougawa, and S. Itoh, 2008: Tropospheric impact of reflected planetary waves from the stratosphere. Geophys. Res. Lett., 35, L16806, https://doi.org/10.1029/2008GL034575.

    • Search Google Scholar
    • Export Citation
  • Kotsuki, S., K. Kurosawa, and T. Miyoshi, 2019: On the properties of ensemble forecast sensitivity to observations. Quart. J. Roy. Meteor. Soc., 145, 18971914, https://doi.org/10.1002/qj.3534.

    • Search Google Scholar
    • Export Citation
  • Langland, R. H., and N. L. Baker, 2004: Estimation of observation impact using the NRL atmospheric variational data assimilation adjoint system. Tellus, 56A, 189201, https://doi.org/10.3402/tellusa.v56i3.14413.

    • Search Google Scholar
    • Export Citation
  • Lawrence, H., N. Bormann, I. Sandu, J. Day, J. Farnan, and P. Bauer, 2019: Use and impact of Arctic observations in the ECMWF numerical weather prediction system. Quart. J. Roy. Meteor. Soc., 145, 34323454, https://doi.org/10.1002/qj.3628.

    • Search Google Scholar
    • Export Citation
  • Lien, G.-Y., D. Hotta, E. Kalnay, T. Miyoshi, and T.-C. Chen, 2018: Accelerating assimilation development for new observing systems using EFSO. Nonlinear Processes Geophys., 25, 129143, https://doi.org/10.5194/npg-25-129-2018.

    • Search Google Scholar
    • Export Citation
  • Lim, E.-P., and Coauthors, 2021: The 2019 Southern Hemisphere stratospheric polar vortex weakening and its impacts. Bull. Amer. Meteor. Soc., 102, E1150E1171, https://doi.org/10.1175/BAMS-D-20-0112.1.

    • Search Google Scholar
    • Export Citation
  • McIntyre, M. E., 1982: How well do we understand the dynamics of stratospheric warmings? J. Meteor. Soc. Japan, 60, 3765, https://doi.org/10.2151/jmsj1965.60.1_37.

    • Search Google Scholar
    • Export Citation
  • Miyoshi, T., and Y. Sato, 2007: Assimilating satellite radiances with a Local Ensemble Transform Kalman Filter (LETKF) applied to the JMA global model (GSM). SOLA, 3, 3740, https://doi.org/10.2151/sola.2007-010.

    • Search Google Scholar
    • Export Citation
  • Miyoshi, T., Y. Sato, and T. Kadowaki, 2010: Ensemble Kalman filter and 4D-Var intercomparison with the Japanese operational global analysis and prediction system. Mon. Wea. Rev., 138, 28462866, https://doi.org/10.1175/2010MWR3209.1.

    • Search Google Scholar
    • Export Citation
  • Mukougawa, H., S. Noguchi, Y. Kuroda, R. Mizuta, and K. Kodera, 2017: Dynamics and predictability of downward-propagating stratospheric planetary waves observed in March 2007. J. Atmos. Sci., 74, 35333550, https://doi.org/10.1175/JAS-D-16-0330.1.

    • Search Google Scholar
    • Export Citation
  • Nishii, K., and H. Nakamura, 2010: Three-dimensional evolution of ensemble forecast spread during the onset of a stratospheric sudden warming event in January 2006. Quart. J. Roy. Meteor. Soc., 136, 894905, https://doi.org/10.1002/qj.607.

    • Search Google Scholar
    • Export Citation
  • Noguchi, S., Y. Kuroda, K. Kodera, and S. Watanabe, 2020: Robust enhancement of tropical convective activity by the 2019 Antarctic sudden stratospheric warming. Geophys. Res. Lett., 47, e2020GL088743, https://doi.org/10.1029/2020GL088743.

    • Search Google Scholar
    • Export Citation
  • Orr, A., P. Bechtold, J. Scinocca, M. Ern, and M. Janisková, 2010: Improved middle atmosphere climate and forecasts in the ECMWF model through a non-orographic gravity wave drag parameterization. J. Climate, 23, 59055926, https://doi.org/10.1175/2010JCLI3490.1.

    • Search Google Scholar
    • Export Citation
  • Ota, Y., J. C. Derber, T. Miyoshi, and E. Kalnay, 2013: Ensemble-based observation impact estimates using the NCEP GFS. Tellus, 65A, 20038, https://doi.org/10.3402/tellusa.v65i0.20038.

    • Search Google Scholar
    • Export Citation
  • Privé, N. C., R. M. Errico, R. Todling, and A. E. Akkraoui, 2021: Evaluation of adjoint-based observation impacts as a function of forecast length using an observing system simulation experiment. Quart. J. Roy. Meteor. Soc., 147, 121138, https://doi.org/10.1002/qj.3909.

    • Search Google Scholar
    • Export Citation
  • Reynolds, R. W., T. M. Smith, C. Liu, D. B. Chelton, K. S. Casey, and M. G. Schlax, 2007: Daily high-resolution-blended analyses for sea surface temperature. J. Climate, 20, 54735496, https://doi.org/10.1175/2007JCLI1824.1.

    • Search Google Scholar
    • Export Citation
  • Sato, K., J. Inoue, A. Yamazaki, J.-H. Kim, A. Makshtas, V. Kustov, M. Maturilli, and K. Dethloff, 2018: Impact on predictability of tropical and mid-latitude cyclones by extra Arctic observations. Sci. Rep., 8, 12104, https://doi.org/10.1038/s41598-018-30594-4.

    • Search Google Scholar
    • Export Citation
  • Sato, K., J. Inoue, and A. Yamazaki, 2020: Performance of forecasts of hurricanes with and without upper-level troughs over the mid-latitudes. Atmosphere, 11, 702, https://doi.org/10.3390/atmos11070702.

    • Search Google Scholar
    • Export Citation
  • Saunders, R., and Coauthors, 2013: RTTOV-11 science and validation report. EUMETSAT Tech. Rep. NWPSAF-MO-TV-032, 62 pp., https://www.nwpsaf.eu/site/download/documentation/rtm/docs_rttov11/rttov11_svr.pdf.

  • Scinocca, J. F., 2003: An accurate spectral nonorographic gravity wave drag parameterization for general circulation models. J. Atmos. Sci., 60, 667682, https://doi.org/10.1175/1520-0469(2003)060<0667:AASNGW>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Sekiguchi, M., and T. Nakajima, 2008: A k-distribution-based radiation code and its computational optimization for an atmospheric general circulation model. J. Quant. Spectrosc. Radiat. Transfer, 109, 27792793, https://doi.org/10.1016/j.jqsrt.2008.07.013.

    • Search Google Scholar
    • Export Citation
  • Terasaki, K., and T. Miyoshi, 2017: Assimilating AMSU-A radiances with the NICAM-LETKF. J. Meteor. Soc. Japan, 95, 433446, https://doi.org/10.2151/jmsj.2017-028.

    • Search Google Scholar
    • Export Citation
  • Terasaki, K., S. Kotsuki, and T. Miyoshi, 2019: Multi-year analysis using the NICAM-LETKF data assimilation system. SOLA, 15, 4146, https://doi.org/10.2151/sola.2019-009.

    • Search Google Scholar
    • Export Citation
  • Thorncroft, C. D., B. J. Hoskins, and M. E. McIntyre, 1993: Two paradigms of baroclinic-wave life-cycle behaviour. Quart. J. Roy. Meteor. Soc., 119, 1755, https://doi.org/10.1002/qj.49711950903.

    • Search Google Scholar
    • Export Citation
  • Yamazaki, A., J. Inoue, K. Dethloff, M. Maturilli, and G. König-Langlo, 2015: Impact of radiosonde observations on forecasting summertime Arctic cyclone formation. J. Geophys. Res. Atmos., 120, 32493273, https://doi.org/10.1002/2014JD022925.

    • Search Google Scholar
    • Export Citation
  • Yamazaki, A., T. Miyoshi, J. Inoue, T. Enomoto, and N. Komori, 2021: EFSO at different geographical locations verified with observing system experiments. Wea. Forecasting, 36, 12191236, https://doi.org/10.1175/WAF-D-20-0152.1.

    • Search Google Scholar
    • Export Citation

Supplementary Materials

Save
  • Andrews, D. G., J. R. Holton, and C. B. Leovy, 1987: Middle Atmosphere Dynamics. Academic Press, 489 pp.

  • Baba, Y., 2019: Spectral cumulus parameterization based on cloud-resolving model. Climate Dyn., 52, 309334, https://doi.org/10.1007/s00382-018-4137-z.

    • Search Google Scholar
    • Export Citation
  • Baba, Y., 2020: Shallow convective closure in a spectral cumulus parameterization. Atmos. Res., 233, 104707, https://doi.org/10.1016/j.atmosres.2019.104707.

    • Search Google Scholar
    • Export Citation
  • Baba, Y., and T. Ogata, 2022: Resolution dependence of tropical cyclones simulated by a spectral cumulus parameterization. Dyn. Atmos. Ocean, 97, 101283, https://doi.org/10.1016/j.dynatmoce.2022.101283.

    • Search Google Scholar
    • Export Citation
  • Candy, B., J. Cotton, and J. Eyre, 2021: Recent results of observation data denial experiments. Met Office Weather Science Tech. Rep. 641, 35 pp., https://www.metoffice.gov.uk/binaries/content/assets/metofficegovuk/pdf/research/weather-science/frtr_641_2021p.pdf.

  • Chen, T.-C., and E. Kalnay, 2019: Proactive quality control: Observing system simulation experiments with the Lorenz’96 model. Mon. Wea. Rev., 147, 5367, https://doi.org/10.1175/MWR-D-18-0138.1.

    • Search Google Scholar
    • Export Citation
  • Chen, T.-C., and E. Kalnay, 2020: Proactive quality control: Observing system experiments using the NCEP Global Forecast System. Mon. Wea. Rev., 148, 39113931, https://doi.org/10.1175/MWR-D-20-0001.1.

    • Search Google Scholar
    • Export Citation
  • Edmon, H. J., Jr., B. J. Hoskins, and M. E. McIntyre, 1980: Eliassen–Palm cross-sections for the troposphere. J. Atmos. Sci., 37, 26002616, https://doi.org/10.1175/1520-0469(1980)037<2600:EPCSFT>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Ehrendorfer, M., R. M. Errico, and K. D. Raeder, 1999: Singular-vector perturbation growth in a primitive equation model with moist physics. J. Atmos. Sci., 56, 16271648, https://doi.org/10.1175/1520-0469(1999)056<1627:SVPGIA>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Emanuel, K. A., 1991: A scheme for representing cumulus convection in large-scale models. J. Atmos. Sci., 48, 23132329, https://doi.org/10.1175/1520-0469(1991)048<2313:ASFRCC>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Enomoto, T., A. Kuwano-Yoshida, N. Komori, and W. Ohfuchi, 2008: Description of AFES 2: Improvements for high-resolution and coupled simulations. High Resolution Numerical Modelling of the Atmosphere and Ocean, K. Hamilton and W. Ohfuchi, Eds., Springer, 77–97, https://doi.org/10.1007/978-0-387-49791-4_5.

  • Enomoto, T., T. Miyoshi, Q. Moteki, J. Inoue, M. Hattori, A. Kuwano-Yoshida, N. Komori, and S. Yamane, 2013: Observing-system research and ensemble data assimilation at JAMSTEC. Data Assimilation for Atmospheric, Oceanic and Hydrologic Applications (Vol. II), S. Park and L. Xu, Eds., Springer, 509–526, https://doi.org/10.1007/978-3-642-35088-7_21.

  • Eyre, J. R., 2021: Observation impact metrics in NWP: A theoretical study. Part I: Optimal systems. Quart. J. Roy. Meteor. Soc., 147, 31803200, https://doi.org/10.1002/qj.4123.

    • Search Google Scholar
    • Export Citation
  • Fujiwara, M., and Coauthors, 2017: Introduction to the SPARC Reanalysis Intercomparison Project (S-RIP) and overview of the reanalysis systems. Atmos. Chem. Phys., 17, 14171452, https://doi.org/10.5194/acp-17-1417-2017.

    • Search Google Scholar
    • Export Citation
  • Fujiwara, M., G. L. Manney, L. J. Gray, and J. S. Wright, Eds., 2022: SPARC Reanalysis Intercomparison Project (S-RIP) final report. SPARC Rep. 10, 635 pp., WCRP-17/2020, https://doi.org/10.17874/800dee57d13.

  • Gelaro, R., and Y. Zhu, 2009: Examination of observation impacts derived from observing system experiments (OSEs) and adjoint models. Tellus, 61A, 179193, https://doi.org/10.1111/j.1600-0870.2008.00388.x.

    • Search Google Scholar
    • Export Citation
  • Gelaro, R., R. H. Langland, S. Pellerin, and R. Todling, 2010: The THORPEX observation impact intercomparison experiment. Mon. Wea. Rev., 138, 40094025, https://doi.org/10.1175/2010MWR3393.1.

    • Search Google Scholar
    • Export Citation
  • Hattori, M., A. Yamazaki, S.-Y. Ogino, P. Wu, and J. Matsumoto, 2017: Impact of the radiosonde observations of cold surge over the Philippine Sea on the tropical region and the Southern Hemisphere in December 2012. SOLA, 13, 1924, https://doi.org/10.2151/sola.2017-004.

    • Search Google Scholar
    • Export Citation
  • Hersbach, H., and Coauthors, 2020: The ERA5 global reanalysis. Quart. J. Roy. Meteor. Soc., 146, 19992049, https://doi.org/10.1002/qj.3803.

    • Search Google Scholar
    • Export Citation
  • Hotta, D., T.-C. Chen, E. Kalnay, Y. Ota, and T. Miyoshi, 2017: Proactive QC: A fully flow-dependent quality control scheme based on EFSO. Mon. Wea. Rev., 145, 33313354, https://doi.org/10.1175/MWR-D-16-0290.1.

    • Search Google Scholar
    • Export Citation
  • Houtekamer, P. L., and F. Zhang, 2016: Review of the ensemble Kalman filter for atmospheric data assimilation. Mon. Wea. Rev., 144, 44894532, https://doi.org/10.1175/MWR-D-15-0440.1.

    • Search Google Scholar
    • Export Citation
  • Inoue, J., A. Yamazaki, J. Ono, K. Dethloff, M. Maturilli, R. Neuber, P. Edwards, and H. Yamaguchi, 2015: Additional Arctic observations improve weather and sea-ice forecasts for the northern sea route. Sci. Rep., 5, 16868, https://doi.org/10.1038/srep16868.

    • Search Google Scholar
    • Export Citation
  • Kalnay, E., and Coauthors, 1996: The NCEP/NCAR 40-Year Reanalysis Project. Bull. Amer. Meteor. Soc., 77, 437472, https://doi.org/10.1175/1520-0477(1996)077<0437:TNYRP>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Kalnay, E., Y. Ota, T. Miyoshi, and J. Liu, 2012: A simpler formulation of forecast sensitivity to observations: Application to ensemble Kalman filters. Tellus, 64A, 18462, https://doi.org/10.3402/tellusa.v64i0.18462.

    • Search Google Scholar
    • Export Citation
  • Kobayashi, S., and Coauthors, 2015: The JRA-55 reanalysis: General specifications and basic characteristics. J. Meteor. Soc. Japan, 93, 548, https://doi.org/10.2151/jmsj.2015-001.

    • Search Google Scholar
    • Export Citation
  • Kodama, C., A. Kuwano-Yoshida, S. Watanabe, T. Doi, H. Kashimura, and T. Nasuno, 2019: JAMSTEC model intercomparision project (JMIP). JAMSTEC Rep. Res. Dev., 28, 534, https://doi.org/10.5918/jamstecr.28.5.

    • Search Google Scholar
    • Export Citation
  • Kodera, K., H. Mukougawa, and S. Itoh, 2008: Tropospheric impact of reflected planetary waves from the stratosphere. Geophys. Res. Lett., 35, L16806, https://doi.org/10.1029/2008GL034575.

    • Search Google Scholar
    • Export Citation
  • Kotsuki, S., K. Kurosawa, and T. Miyoshi, 2019: On the properties of ensemble forecast sensitivity to observations. Quart. J. Roy. Meteor. Soc., 145, 18971914, https://doi.org/10.1002/qj.3534.

    • Search Google Scholar
    • Export Citation
  • Langland, R. H., and N. L. Baker, 2004: Estimation of observation impact using the NRL atmospheric variational data assimilation adjoint system. Tellus, 56A, 189201, https://doi.org/10.3402/tellusa.v56i3.14413.

    • Search Google Scholar
    • Export Citation
  • Lawrence, H., N. Bormann, I. Sandu, J. Day, J. Farnan, and P. Bauer, 2019: Use and impact of Arctic observations in the ECMWF numerical weather prediction system. Quart. J. Roy. Meteor. Soc., 145, 34323454, https://doi.org/10.1002/qj.3628.

    • Search Google Scholar
    • Export Citation
  • Lien, G.-Y., D. Hotta, E. Kalnay, T. Miyoshi, and T.-C. Chen, 2018: Accelerating assimilation development for new observing systems using EFSO. Nonlinear Processes Geophys., 25, 129143, https://doi.org/10.5194/npg-25-129-2018.

    • Search Google Scholar
    • Export Citation
  • Lim, E.-P., and Coauthors, 2021: The 2019 Southern Hemisphere stratospheric polar vortex weakening and its impacts. Bull. Amer. Meteor. Soc., 102, E1150E1171, https://doi.org/10.1175/BAMS-D-20-0112.1.

    • Search Google Scholar
    • Export Citation
  • McIntyre, M. E., 1982: How well do we understand the dynamics of stratospheric warmings? J. Meteor. Soc. Japan, 60, 3765, https://doi.org/10.2151/jmsj1965.60.1_37.

    • Search Google Scholar
    • Export Citation
  • Miyoshi, T., and Y. Sato, 2007: Assimilating satellite radiances with a Local Ensemble Transform Kalman Filter (LETKF) applied to the JMA global model (GSM). SOLA, 3, 3740, https://doi.org/10.2151/sola.2007-010.

    • Search Google Scholar
    • Export Citation
  • Miyoshi, T., Y. Sato, and T. Kadowaki, 2010: Ensemble Kalman filter and 4D-Var intercomparison with the Japanese operational global analysis and prediction system. Mon. Wea. Rev., 138, 28462866, https://doi.org/10.1175/2010MWR3209.1.

    • Search Google Scholar
    • Export Citation
  • Mukougawa, H., S. Noguchi, Y. Kuroda, R. Mizuta, and K. Kodera, 2017: Dynamics and predictability of downward-propagating stratospheric planetary waves observed in March 2007. J. Atmos. Sci., 74, 35333550, https://doi.org/10.1175/JAS-D-16-0330.1.

    • Search Google Scholar
    • Export Citation
  • Nishii, K., and H. Nakamura, 2010: Three-dimensional evolution of ensemble forecast spread during the onset of a stratospheric sudden warming event in January 2006. Quart. J. Roy. Meteor. Soc., 136, 894905, https://doi.org/10.1002/qj.607.

    • Search Google Scholar
    • Export Citation
  • Noguchi, S., Y. Kuroda, K. Kodera, and S. Watanabe, 2020: Robust enhancement of tropical convective activity by the 2019 Antarctic sudden stratospheric warming. Geophys. Res. Lett., 47, e2020GL088743, https://doi.org/10.1029/2020GL088743.

    • Search Google Scholar
    • Export Citation
  • Orr, A., P. Bechtold, J. Scinocca, M. Ern, and M. Janisková, 2010: Improved middle atmosphere climate and forecasts in the ECMWF model through a non-orographic gravity wave drag parameterization. J. Climate, 23, 59055926, https://doi.org/10.1175/2010JCLI3490.1.

    • Search Google Scholar
    • Export Citation
  • Ota, Y., J. C. Derber, T. Miyoshi, and E. Kalnay, 2013: Ensemble-based observation impact estimates using the NCEP GFS. Tellus, 65A, 20038, https://doi.org/10.3402/tellusa.v65i0.20038.

    • Search Google Scholar
    • Export Citation
  • Privé, N. C., R. M. Errico, R. Todling, and A. E. Akkraoui, 2021: Evaluation of adjoint-based observation impacts as a function of forecast length using an observing system simulation experiment. Quart. J. Roy. Meteor. Soc., 147, 121138, https://doi.org/10.1002/qj.3909.

    • Search Google Scholar
    • Export Citation
  • Reynolds, R. W., T. M. Smith, C. Liu, D. B. Chelton, K. S. Casey, and M. G. Schlax, 2007: Daily high-resolution-blended analyses for sea surface temperature. J. Climate, 20, 54735496, https://doi.org/10.1175/2007JCLI1824.1.

    • Search Google Scholar
    • Export Citation
  • Sato, K., J. Inoue, A. Yamazaki, J.-H. Kim, A. Makshtas, V. Kustov, M. Maturilli, and K. Dethloff, 2018: Impact on predictability of tropical and mid-latitude cyclones by extra Arctic observations. Sci. Rep., 8, 12104, https://doi.org/10.1038/s41598-018-30594-4.

    • Search Google Scholar
    • Export Citation
  • Sato, K., J. Inoue, and A. Yamazaki, 2020: Performance of forecasts of hurricanes with and without upper-level troughs over the mid-latitudes. Atmosphere, 11, 702, https://doi.org/10.3390/atmos11070702.

    • Search Google Scholar
    • Export Citation
  • Saunders, R., and Coauthors, 2013: RTTOV-11 science and validation report. EUMETSAT Tech. Rep. NWPSAF-MO-TV-032, 62 pp., https://www.nwpsaf.eu/site/download/documentation/rtm/docs_rttov11/rttov11_svr.pdf.

  • Scinocca, J. F., 2003: An accurate spectral nonorographic gravity wave drag parameterization for general circulation models. J. Atmos. Sci., 60, 667682, https://doi.org/10.1175/1520-0469(2003)060<0667:AASNGW>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Sekiguchi, M., and T. Nakajima, 2008: A k-distribution-based radiation code and its computational optimization for an atmospheric general circulation model. J. Quant. Spectrosc. Radiat. Transfer, 109, 27792793, https://doi.org/10.1016/j.jqsrt.2008.07.013.

    • Search Google Scholar
    • Export Citation
  • Terasaki, K., and T. Miyoshi, 2017: Assimilating AMSU-A radiances with the NICAM-LETKF. J. Meteor. Soc. Japan, 95, 433446, https://doi.org/10.2151/jmsj.2017-028.

    • Search Google Scholar
    • Export Citation
  • Terasaki, K., S. Kotsuki, and T. Miyoshi, 2019: Multi-year analysis using the NICAM-LETKF data assimilation system. SOLA, 15, 4146, https://doi.org/10.2151/sola.2019-009.

    • Search Google Scholar
    • Export Citation
  • Thorncroft, C. D., B. J. Hoskins, and M. E. McIntyre, 1993: Two paradigms of baroclinic-wave life-cycle behaviour. Quart. J. Roy. Meteor. Soc., 119, 1755, https://doi.org/10.1002/qj.49711950903.

    • Search Google Scholar
    • Export Citation
  • Yamazaki, A., J. Inoue, K. Dethloff, M. Maturilli, and G. König-Langlo, 2015: Impact of radiosonde observations on forecasting summertime Arctic cyclone formation. J. Geophys. Res. Atmos., 120, 32493273, https://doi.org/10.1002/2014JD022925.

    • Search Google Scholar
    • Export Citation
  • Yamazaki, A., T. Miyoshi, J. Inoue, T. Enomoto, and N. Komori, 2021: EFSO at different geographical locations verified with observing system experiments. Wea. Forecasting, 36, 12191236, https://doi.org/10.1175/WAF-D-20-0152.1.

    • Search Google Scholar
    • Export Citation
  • Fig. 1.

    (a),(b),(d),(e) Time series of the globally averaged root-mean-square difference of ALERA (thick black lines), NCEP–NCAR reanalysis (gray lines), and ERA5 (thin black lines) relative to JRA-55 for (top) 300-hPa zonal wind (m s−1) and (middle) 850-hPa temperature (K) for the periods of December 2018–January 2019 in (a) and (b) and August–September 2019 in (d) and (e). (bottom) Snapshots of 500-hPa geopotential height fields in ALERA (red solid lines) and JRA-55 (black dashed lines) at (c) 0000 UTC 21 Dec 2018 and (f) 0000 UTC 23 Aug 2019.

  • Fig. 2.

    Time series of (a),(d) total OI in the CTL (black lines) and OSE (blue dashed lines); (b),(e) differences between 6-h (thin lines) and 12-h (thick lines) forecasts and their analyses (J kg−1) in the CTL (black lines) and OSE (blue dashed lines); and (c),(f) differences between the CTL and OSE analyses (J kg−1; black lines) and the analysis ensemble spreads (J kg−1) in the CTL (red upward-pointing triangles) and the OSEs (blue downward-pointing triangles). (left) NH-winter experiment: 21 Dec 2018–31 Jan 2019 and (right) SH-winter experiment: 23 Aug–3 Oct 2019.