Abstract

This study evaluates the relationship between the Madden–Julian oscillation (MJO) and the occurrence of equatorial Pacific westerly wind bursts (WWBs). During the convective MJO phase, anomalous surface westerlies prevail in and west of the convective MJO center, providing favorable conditions for WWBs. Compared with the probability of WWBs expected under a null hypothesis that WWBs occur randomly, the convective MJO phase almost doubles the probability of a WWB occurring. Nevertheless, only 34.46% of WWBs co-occur with the convective MJO, which is much less than that reported in previous studies. We show that when the MJO and WWBs are defined using the same field with overlapping frequencies, the percentage of WWBs co-occurring with the convective MJO shows a significant increase. However, the higher percentage is simply caused by the fact that the strong WWBs during a convective MJO are more likely to be identified than those during the suppressed and neutral MJO phases. A total of 45.80% of WWBs are found to occur in the full MJO phase (both the convective and suppressed MJO phases), which is slightly higher than that expected based on randomness. Although the full MJO has statistically significant impact on the likelihood of WWBs, the influence from the full MJO on the tropical Pacific sea surface temperature anomaly is much weaker as compared to that from the WWBs. The relationships between the MJO and WWBs simulated in CMIP5 models are also assessed, and the percentage of WWBs that co-occur with the MJO simulated in models is in general less than that in observations.

Abstract

The dropsondes released during the Tropical Cyclone Intensity (TCI) field campaign provide high-resolution kinematic and thermodynamic measurements of tropical cyclones within the upper-level outflow and inner core. This study investigates the impact of these upper-level TCI dropsondes on analyses and prediction of Hurricane Patricia (2015) during its rapid intensification (RI) phase using an ensemble–variational data assimilation system. In the baseline experiment (BASE), both kinematic and thermodynamic observations of TCI dropsondes at all levels except the upper levels are assimilated. The upper-level wind and thermodynamic observations are assimilated in additional experiments to investigate their respective impacts. Compared to BASE, assimilating TCI upper-level wind observations improves the accuracy of outflow analyses verified against independent atmospheric motion vector (AMV) observations. It also strengthens the tangential and radial wind near the upper-level eyewall. The inertial stability within the upper-level eyewall is enhanced, and the maximum outflow is more aligned toward the inner core. Additionally, the analyses including the upper-level thermodynamic observations produce a warmer and drier core at high levels. Assimilating both upper-level kinematic and thermodynamic observations also improves the RI forecast. Compared to BASE, assimilating the upper-level wind induces more upright and inward-located eyewall convection, resulting in more latent heat release closer to the warm core. This process leads to stronger inner-core warming. Additionally, the initial warmer upper-level inner core produced by assimilating TCI thermodynamic observations also intensifies the convection and latent heat release within the eyewall, thus further contributing to the improved intensity forecasts.

Abstract

Although numerous studies have demonstrated that increasing model spatial resolution in free forecasts can potentially improve tropical cyclone (TC) intensity forecasts, studies on the impact of model resolution during data assimilation (DA) on TC prediction are lacking. In this study, using the ensemble-variational DA system for the Hurricane Weather Research and Forecasting (HWRF) Model, we investigated the individual impact of increasing the model resolution of first guess (FG) and background ensemble (BE) forecasts during DA on initial analyses and subsequent forecasts of Hurricane Patricia (2015). The impacts were compared between horizontal and vertical resolutions and also between the tropical storm (TS) and hurricane assimilation during Patricia. The results show that increasing the horizontal or vertical resolution in FG has a larger impact than increasing the resolution in BE on improving the analyzed TC intensity and structure for the hurricane stage. The result is reversed for the TS stage. These results are attributed to the effectiveness of increasing the FG resolution in intensifying the background vortex for the hurricane stage relative to the TS stage. Increasing the BE resolution contributes to improving the analyzed intensity through the better-resolved background correlation structure for both the hurricane and TS stages. Increasing horizontal resolution has an overall larger effect than increasing vertical resolution in improving the analysis at the hurricane stage and their effects are close for the analysis at the TS stage. Additionally, the more accurately analyzed primary circulation, secondary circulation, and warm-core structures via the increased resolution in DA lead to improved TC intensity forecasts.

Abstract

The local particle filter (LPF) and the local nonlinear ensemble transform filter (LNETF) are two moment-matching nonlinear filters to approximate the classical particle filter (PF). They adopt different strategies to alleviate filter degeneracy. LPF and LNETF localize observational impact but use different localization functions. They assimilate observations in a partially sequential and a simultaneous manner, respectively. In addition, LPF applies the resampling step, whereas LNETF applies the deterministic square root transformation to update particles. Both methods preserve the posterior mean and variance of the PF. LNETF additionally preserves the posterior correlation of the PF for state variables within a local volume. These differences lead to their differing performance in filter stability and posterior moment estimation. LPF and LNETF are systematically compared and analyzed here through a set of experiments with a Lorenz model. Strategies to improve the LNETF are proposed. The original LNETF is inferior to the original LPF in filter stability and analysis accuracy, particularly for small particle numbers. This is attributed to both the localization function and particle update differences. The LNETF localization function imposes a stronger observation impact than the LPF for remote grids and thus is more susceptible to filter degeneracy. The LNETF update causes an overall narrower range of posteriors that excludes true states more frequently. After applying the same localization function as the LPF and additional posterior inflation to the LNETF, the two filters reach similar filter stability and analysis accuracy for all particle numbers. The improved LNETF shows more accurate posterior probability distribution but slightly worse spatial correlation of posteriors than the LPF.

Abstract

The relatively weak sea surface temperature bias in the tropical Indian Ocean (TIO) simulated in the coupled general circulation model (CGCM) from the recently released CMIP6 has been found to be important in model simulations of regional and global climate. However, the cause of the bias is debated because the bias is strongly model dependent and shows marked seasonality. In this study, we separate the bias in CGCMs into bias arising from oceanic GCMs (OGCMs) and bias that is independent of OGCMs using a set of CMIP6 and OMIP6 models. We found that OGCMs contribute little to mixed layer bias in the CGCMs. The OGCM-independent bias exhibits a large-scale cold mixed layer bias in the TIO throughout the year, with an unexpectedly high degree of model consistency. By conducting a set of OGCM experiments, we show that the OGCM-independent mixed layer bias is caused mainly by surface wind bias in the utilized CGCMs. About 89% of surface wind bias in the CGCMs is due to the inability of atmospheric GCMs (AGCMs), whereas atmosphere–ocean coupling in the CGCMs has only a minor influence on surface wind bias. The bias in surface wind is also found to be the cause of subsurface temperature bias besides the ocean dynamics such as vertical mixing and vertical shear in currents. Our results indicate that correcting TIO mixed layer bias in CGCMs requires improvement in the capability of AGCM in simulating the climatological surface winds.

Abstract

Instabilities play a critical role in understanding atmospheric predictability and improving weather forecasting. The bred vectors (BVs) are dynamically evolved and flow-dependent nonlinear perturbations, indicating the most unstable modes of the underlying flow. Especially over smaller areas, however, BVs with different initial seeds may to some extent be constrained to a small subspace, missing potential forecast error growth along other unstable perturbation directions.

In this paper, the authors study the nonlinear local Lyapunov vectors (NLLVs) that are designed to capture an orthogonal basis spanning the most unstable perturbation subspace, thus potentially ameliorating the limitation of BVs. The NLLVs are theoretically related to the linear Lyapunov vectors (LVs), which also form an orthogonal basis. Like BVs, NLLVs are generated by dynamically evolving perturbations with a full nonlinear model. In simulated forecast experiments, a set of mutually orthogonal NLLVs show an advantage in predicting the structure of forecast error growth when compared to using a set of BVs that are not fully independent. NLLVs are also found to have a higher local dimension, enabling them to better capture localized instabilities, leading to increased ensemble spread.

Abstract

Nonlinear local Lyapunov vectors (NLLVs) are developed to indicate orthogonal directions in phase space with different perturbation growth rates. In particular, the first few NLLVs are considered to be an appropriate orthogonal basis for the fast-growing subspace. In this paper, the NLLV method is used to generate initial perturbations and implement ensemble forecasts in simple nonlinear models (the Lorenz63 and Lorenz96 models) to explore the validity of the NLLV method.

The performance of the NLLV method is compared comprehensively and systematically with other methods such as the bred vector (BV) and the random perturbation (Monte Carlo) methods. In experiments using the Lorenz63 model, the leading NLLV (LNLLV) captured a more precise direction, and with a faster growth rate, than any individual bred vector. It may be the larger projection on fastest-growing analysis errors that causes the improved performance of the new method. Regarding the Lorenz96 model, two practical measures—namely the spread–skill relationship and the Brier score—were used to assess the reliability and resolution of these ensemble schemes. Overall, the ensemble spread of NLLVs is more consistent with the errors of the ensemble mean, which indicates the better performance of NLLVs in simulating the evolution of analysis errors. In addition, the NLLVs perform significantly better than the BVs in terms of reliability and the random perturbations in resolution.

Abstract

The observation accuracy of the surface air temperature less than 0.1 K is a requirement, stated by the meteorological and climatological community. However, the accuracy of a temperature sensor inside a shield is affected by a number of factors including solar radiation, wind speed, upwelling longwave radiation, air density, sun elevation angle, sun azimuth angle, underlying surface, precipitation, moisture, structure, and coating of the radiation shield. Due to these factors, the temperature error of the temperature sensor may be much larger than 1 K under adverse conditions. To improve the observation accuracy, this paper proposed a spherical temperature sensor array. A series of analytical calculations based on a computational fluid dynamics (CFD) method is performed to verify the design principle of this sensor array. The calculation results show that the temperature error ratio can be assumed as a constant. To verify the accuracy of this sensor array, simulations and observation experiments are conducted. The simulation results show that the mean difference between the temperature provided by this sensor array and the reference air temperature is 0.072 K. The field experiment results show that a root-mean-square error (RMSE) and a mean absolute error (MAE) between the temperature provided by this sensor array and the reference air temperature are 0.173 and 0.153 K, respectively.

Abstract

Ensemble prediction is a widely used tool in weather forecasting. In particular, the arithmetic mean (AM) of ensemble members is used to filter out unpredictable features from a forecast. AM is a pointwise statistical concept, providing the best sample-based estimate of the expected value of any single variable. The atmosphere, however, is a multivariate system with spatially coherent features characterized with strong correlations. Disregarding such correlations, the AM of an ensemble of forecasts removes not only unpredictable noise but also flattens features whose presence is still predictable, albeit with somewhat uncertain location. As a consequence, AM destroys the structure, and reduces the amplitude and variability associated with partially predictable features. Here we explore the use of an alternative concept of central tendency for the estimation of the expected feature (instead of single values) in atmospheric systems. Features that are coherent across ensemble members are first collocated to their mean position, before the AM of the aligned members is taken. Unlike earlier definitions based on complex variational minimization (field coalescence of Ravela and generalized ensemble mean of Purser), the proposed feature-oriented mean (FM) uses simple and computationally efficient vector operations. Though FM is still not a dynamically realizable state, a preliminary evaluation of ensemble geopotential height forecasts indicates that it retains more variance than AM, without a noticeable drop in skill. Beyond ensemble forecasting, possible future applications include a wide array of climate studies where the collocation of larger-scale features of interest may yield enhanced compositing results.

Abstract

Whether the state-of-the-art CMIP5 models have different El Niño types and how the degree of modeled El Niño diversity would be impacted by the future global warming are still heavily debated. In this study, cluster analysis is used to investigate El Niño diversity in 30 CMIP5 models. As the method does not rely on any prior knowledge of the patterns of El Niño seen in observations, it provides a practical way to identify the degree of El Niño diversity in models. Under the historical scenario, most models show a poor degree of El Niño diversity in their own model world, primarily due to the lopsided numbers of events belonging to the two modeled El Niño types and the weak compactness of events in each cluster. Four models are found showing significant El Niño diversity, yet none of them captures the longitudinal distributions of the warming centers of the two El Niño types seen in the observations. Heat budget analysis of the sea surface temperature (SST) anomaly suggests that the degree of modeled El Niño diversity is highly related to the climatological zonal SST gradient over the western-central equatorial Pacific in models. As the gradient is weakened in most models under the future high-emission scenario, the degree of modeled El Niño diversity is further reduced in the future. The results indicate that a better simulation of the SST gradient over the western-central equatorial Pacific might allow a more reliable simulation/projection of El Niño diversity in most CMIP5 models.