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

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 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

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 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

This paper introduces a climate feedback kernel, referred to as the “energy gain kernel” (EGK). EGK allows for separating the net longwave radiative energy perturbations given by a Planck feedback matrix explicitly into thermal energy emission perturbations of individual layers and thermal radiative energy flux convergence perturbations at individual layers resulting from the coupled atmosphere–surface temperature changes in response to the unit forcing in individual layers. The former is represented by the diagonal matrix of a Planck feedback matrix and the latter by EGK. Elements of EGK are all positive, representing amplified energy perturbations at a layer where forcing is imposed and energy gained at other layers, both of which are achieved through radiative thermal coupling within an atmosphere–surface column. Applying EGK to input energy perturbations, whether external or internal due to responses of nontemperature feedback processes to external energy perturbations, such as water vapor and albedo feedbacks, yields their total energy perturbations amplified through radiative thermal coupling within an atmosphere–surface column. As the strength of EGK depends exclusively on climate mean states, it offers a solution for effectively and objectively separating control climate state information from climate perturbations for climate feedback studies. Given that an EGK comprises critical climate mean state information on mean temperature, water vapor, clouds, and surface pressure, we envision that the diversity of EGK across different climate models could provide insight into the inquiry of why, under the same anthropogenic greenhouse gas increase scenario, different models yield varying degrees of global mean surface warming.

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 positive phase of the Pacific meridional mode (PMM) is closely related to the onset of El Niño. Previous studies have indicated that positive sea surface temperature anomalies (SSTAs) in the central equatorial Pacific (CEP) during the spring and summer of positive PMM years primarily originate from the northeastern tropical Pacific (NETP) via positive wind–evaporation–SST feedback. We review the evolution of PMM and find weak evidence to support such a linkage. Coupled model experiments show that the positive PMM-regressed SSTAs in the NETP only account for ∼24% of those in the CEP from winter to spring, illustrating the principle that correlation does not necessarily mean causality. The strongest positive PMM SSTAs in the NETP and CEP increase El Niño intensity by 1.07°C, whereas that in the NETP alone increase El Niño intensity by 0.69°C. When the composite SSTAs in the NETP during positive PMM years are used, however, the El Niño intensity is increased merely by 0.17°C. The change in the subsurface temperature in the equatorial Pacific is curtailed for the NETP SSTAs to trigger El Niño, while the wind–evaporation–SST feedback plays a less important role. Our results indicate that the impact of PMM on El Niño might be overestimated by ∼55%. Moreover, a comprehensive understanding about the role of the tropical North Pacific on El Niño can be obtained only when the impact from the western North Pacific is considered.

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.