{'fraction_non_alphanumeric': 0.026855948869223205, 'fraction_numerical': 0.022984267453294002, 'mean_word_length': 4.9129723837209305, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Improved Simulation of the Antarctic Stratospheric Final Warming by Modifying the Orographic Gravity Wave Parameterization in the Beijing Climate Center Atmospheric General Circulation Model

: The Antarctic stratospheric ﬁnal warming (SFW) is usually simulated with a substantial delay in climate models, and the corresponding temperatures in austral spring are lower than observations, implying insu ﬃ cient stratospheric wave drag. To investigate the role of orographic gravity wave drag (GWD) in modeling the Antarctic SFW, in this study the orographic GWD parameterization scheme is modiﬁed in the middle-atmosphere version of the Beijing Climate Center Atmospheric General Circulation Model. A pair of simulations are conducted to compare two orographic GWD schemes in simulating the breakdown of the stratospheric polar vortex over Antarctica. The control simulation with the default orographic GWD scheme exhibits delayed vortex breakdown and the cold-pole bias seen in most climate models. In the simulation with modiﬁed orographic GWD scheme, the simulated vortex breaks down earlier by 8 days, and the associated cold-pole bias is reduced by more than 2 K. The modiﬁed scheme provides stronger orographic GWD in the lower stratosphere, which drives an accelerated polar downwelling branch of the Brewer–Dobson circulation and, in turn, produces adiabatic warming. Our study suggests that modifying orographic GWD parameterizations in climate models would be a valid way of improving the SFW simulation over Antarctica.


Introduction
Climate models have evolved into the stage of "Earth System Models" (ESMs), in which various components of physics, chemistry, and biology are fully coupled to serve multiple purposes, including the simulation of air quality, stratospheric ozone, tropospheric chemistry, and global climate. As a major model hierarchy of ESMs, chemistry-climate models (CCMs) combine physical climate models with an explicit representation of atmospheric chemistry and have been used in investigating the Antarctic ozone depletion problem [1,2]. Despite significant improvements in modelling, the simulated Antarctic stratospheric vortex breaking down too late in the year is a common bias to in almost all CCMs. The breakdown of the Southern Hemisphere polar vortex, hereafter referred to as Antarctic stratospheric final warming (SFW), marks the transition of middle-atmospheric circulation from winter to summer. However, the Antarctic SFW simulated by a number of climate models is 1-2 weeks later than observations at 50 hPa [3][4][5]. As a result, the polar temperatures in climate models are Section 4 summarizes our findings and discusses the uncertainties associated with modifying GW parameterizations to improve the dynamical performance of climate models.

General Description of the Middle-Atmosphere Version of BCC-AGCM
In this study, the numerical simulations were carried out with the middle-atmosphere configuration of the medium-resolution Beijing Climate Center Atmospheric General Circulation Model version 3 (BCC-AGCM3-MR) [25]. The original BCC-AGCM3-MR, which participates in the Coupled Model Intercomparison Project Phase 6 (CMIP6), has a spectral truncation at wavenumber 106 with an associated Gaussian grid of 1.125 • resolution. The CMIP6 configuration of the BCC-AGCM3-MR contains 46 hybrid vertical levels with approximately 800-m vertical spacing throughout the midtroposphere and lower stratosphere. The model top is at about 1.98 hPa and the vertical spacing of the top model layer is about 3.8 km. The CMIP6 version of the BCC-AGCM3-MR uses two distinct GWD parameterizations, which represent GWs generated by subgrid-scale orography [22], and by convection [23], respectively. The convective GW parameterization depends on the deep convective heating properties and the background wind, and benefits a spontaneous stratospheric quasi-biennial oscillation [26]. In the present study, the number of model layers has been increased to 95, with an elevated model top up to 0.01 hPa or approximately 80 km. The typical vertical resolution of the 95-layer grid is 700 m from the midtroposphere to the lower stratosphere and better than 1 km up to the stratopause. In addition to the processes represented in the original BCC-AGCM3-MR, the middle-atmosphere version also includes a GW parameterization for the momentum flux deposition from fronts and jet imbalances [24]. Both of the convective and front-jet GW parameterizations used in this high-top model are source-related and supposed to cover all the GWs of nonorographic origins. Furthermore, the middle-atmosphere configuration employs a parameterization to estimate the turbulent mountain stress caused by unresolved orography, which has been proved to be important in simulating the sudden stratospheric warming events in boreal winter by influencing the propagation of planetary waves into the polar vortex [27]. The upward extension of BCC-AGCM3-MR into the middle atmosphere documented here is in the absence of chemistry, and acts as a base model for future chemistry climate modeling efforts. The primary purpose of this middle-atmosphere model is to gain a climatology of winds and temperatures that is suitable for investigating climate issues related to ozone. One prominent problem is the breakdown timing of the Southern Hemisphere stratospheric polar vortex.

Parameterization of Orographic GWs
In the default BCC-AGCM3-MR and its middle-atmosphere configuration, the drag from dissipating GWs forced by flow over orography is parameterized as in McFarlane (1987) [22]. The McFarlane (1987) scheme is a conventional orographic GWD scheme, in which the wave amplitude is represented in terms of the standard deviation of the subgrid-scale orographic heights (σ). Waves are launched at the lowest model level and propagate vertically according to Lindzen's wave saturation theory [28], which defines the saturation momentum flux as where ρ 0 , k, U, and N are the atmospheric density, the horizontal wave number, the background zonal wind, and the buoyancy frequency, respectively. When a wave becomes saturated and begins to break, the drag within a saturated region is determined from Atmosphere 2020, 11, 576 4 of 17 where α is an efficiency factor which is tunable and is set to 0.125 in this study. To investigate the impact of different orographic GW parameterizations on the simulated Antarctic SFWs, the orographic GWD scheme of Kim and Doyle (2005) [29] is implemented in the middle-atmosphere model used in this study. Compared with the McFarlane (1987) scheme, the Kim and Doyle (2005) scheme incorporates the effects of orographic asymmetry, convexity, and anisotropy, which are linked to the characteristics of the corresponding GWs. The orographic GWs in the Kim and Doyle (2005) scheme are emitted from the reference level, which is defined as the larger of the two values of 2σ and the boundary layer height. A model layer is considered unstable if the minimum Richardson number defined by

Experimental Design
As shown in Table 1, two sets of 35-year Atmospheric Model Intercomparison Project (AMIP)-type simulations were conducted with the same prescribed sea surface boundary conditions using the Hadley Centre sea surface temperatures and sea ice concentrations from 1978 to 2012. The simulation using the default McFarlane (1987) [22] scheme is referred to as M87, and the other employing the Kim and Doyle (2005) [29] scheme is referred to as KD05. For each integration, the first year is treated as the spinup, and the last 34 years are used for analysis. Both runs start from identical initial conditions. Historical forcing data are based on observations and mainly include the monthly updated greenhouse gas concentrations with zonal mean values, the stratospheric aerosols from volcanoes, the time-varying gridded ozone concentrations, the anthropogenic aerosol optical properties, the solar forcing, and the yearly global gridded land use. It should be noted that only one 35-year run for each experiment is performed, which is not enough to take into account the internal atmospheric variability adequately as in the ensemble simulations. In order to evaluate the simulated Antarctic SFWs, wind velocities and temperatures from the Japanese Meteorological Agency 55-year Reanalysis (JRA-55) spanning from 1979 to 2012 are used as observational references [30].

The Transformed Eulerian-Mean (TEM) Framework
In this study, we examine the simulated BDC and forcing from resolved planetary waves in the framework of the TEM equations [31], which can be expressed as where u is the zonal mean zonal wind, ϕ is latitude, a is the radius of the earth, z is log-pressure height, f is the Coriolis parameter, ρ 0 is the density, the overbar represents zonal average, and subscripts represent the derivatives with respect to the given variable. The first and second terms on the right-hand side of Equation (4) denote acceleration associated with the residual mean meridional circulation (v * , w * ), i.e., the BDC: where v and w are the Eulerian mean meridional circulation, θ is potential temperature, and prime represents the departure from zonal mean. The third term on the right-hand side of Equation (4) is the divergence of the Eliassen-Palm (E-P) flux where and F = F (ϕ) , F (z) is the E-P flux vector. The E-P flux divergence ∇ · F measures the forcing from resolved waves. The forcing from other processes including GWD is represented by the X term.

The Austral SFW in the BCC-AGCM
The cold-pole bias in the M87 simulation is significant. Figure 1 presents a comparison of Antarctic temperature averaged over the polar cap (60 • -90 • S) between the M87 simulation and JRA-55 reanalysis data. Figure 1a shows the JRA-55 daily climatology in the period of 1979-2012, and Figure 1b shows the difference between M87 and JRA-55. Stippled areas in this and similar comparisons denote differences that are statistically significant at the 95% confidence level, computed using a Student's t test. The cold bias, which descends from the upper to the lower stratosphere in Antarctic late winter to spring (August-November), ranges from about −2 K to −8 K being significant at the 95% level, associated with the delayed breakdown of the stratospheric polar vortex in the M87 simulation. Results from KD05 present a roughly similar bias ( Figure 1c). In detail, the modification of the orographic gravity wave parameterization has a remarkable influence on the temperature climatology in the Southern Hemisphere. Figure 1d shows the difference in the seasonal cycle of polar cap temperature between KD05 and M87. From August through November, much of the systematic cold bias has been reduced in KD05 by more than 2 K compared to M87. These differences between KD05 and M87 are significant at the 95% level. The remaining cold bias indicates that, despite the remarkable improvement of the temperature climatology, the breakdown of the southern polar vortex in the stratosphere still lags the observed behavior. Note that the comparisons shown in Figure 1b, 1c, and 1d are based on daily climatological averages for the period 1979-2012, and present sub-month perturbations which are closely associated with planetary wave modulation.
Some recent studies have reported the alleviation of the cold-pole bias [10,11]. For example, McLandress et al. (2012) [11] showed that the Southern Hemisphere cold-pole bias was reduced substantially by introducing an ad hoc source of additional orographic GWs centered on 60 • S in the Canadian Middle Atmosphere Model. However, in the present study we have merely replaced the existing McFarlane (1987) [22] orographic GWD scheme by the Kim and Doyle (2005) [29] scheme in the BCC-AGCM, without attempting to address specifically the missing orographic GWD near 60 • S where both schemes present the forcing gap as shown below. reaches the 30-hPa level in M87, whereas it is far below the 30-hPa level in the JRA-55 reanalysis data. Figure 2d shows the difference in the seasonal cycle of subpolar zonal wind between KD05 and JRA-55. While significant differences remain, much of the systematic high wind bias has been removed. During the transition from austral winter to summer in KD05, the high wind bias is reduced by more than 4 m s −1 compared to the M87 simulation, with these differences being significant at the 95% confidence level (Figure 2f). It can be expected that the date of the transition to easterlies is improved, indicated by the weaker winds in November in KD05.    Figure 2d shows the difference in the seasonal cycle of subpolar zonal wind between KD05 and JRA-55. While significant differences remain, much of the systematic high wind bias has been removed. During the transition from austral winter to summer in KD05, the high wind bias is reduced by more than 4 m s −1 compared to the M87 simulation, with these differences being significant at the 95% confidence level (Figure 2f). It can be expected that the date of the transition to easterlies is improved, indicated by the weaker winds in November in KD05. Figure 3 presents a complementary view of the zonal wind evolution in the austral spring, exhibiting the monthly zonal mean zonal wind for September, October, and November in JRA-55, M87, and KD05. Significant differences can be seen between M87 and JRA-55, such as a stratospheric jet that is up to 20 m s −1 stronger in M87 than in JRA-55 between 100 and 1 hPa in September and October. The right column of Figure 3 shows the December-February seasonal average of the zonal-mean zonal wind. Apart from the weaker subtropical jet in the Northern Hemisphere, M87 performs reasonably well in the austral summer. Also note that the stratospheric polar-night jet simulated in M87 is stronger and more upright than that in JRA-55, leading to a shift in the position of the simulated polar vortex in boreal winter. Inclusion of the modified orographic GW parameterization results in better agreement between KD05 and JRA-55. The austral spring stratospheric jet is much weaker in KD05, especially in September and October, when the differences between the simulated and observed winds are reduced to about 10 m s −1 . In the austral summer (December-February, Figure 3l), it is interesting to note that the Southern Hemisphere winds in KD05 are very similar to those in M87, while the Northern Hemisphere winds are improved, as shown by a stronger subtropical jet which is closer to the reanalysis data.
The improved simulation of the breakdown of the Southern Hemisphere stratospheric polar vortex is better illustrated in Figure 4, which shows the mean final warming dates at 50 hPa. Following Black and McDaniel (2007) [32], 5-day running means of daily data are used to calculate the final warming date as the final time that the zonal mean zonal wind at 60 • S and 50 hPa drops below 10 m s −1 and remains below that value through the following summer. The mean final warming date for the M87 simulation is about 24 days later than that for JRA-55, which is consistent with the results shown in  Figure 3 presents a complementary view of the zonal wind evolution in the austral spring, exhibiting the monthly zonal mean zonal wind for September, October, and November in JRA-55, M87, and KD05. Significant differences can be seen between M87 and JRA-55, such as a stratospheric jet that is up to 20 m s −1 stronger in M87 than in JRA-55 between 100 and 1 hPa in September and October. The right column of Figure 3 shows the December-February seasonal average of the   Figure 5 shows the latitude-height cross sections of orographic, nonorographic, and total GWD in the M87 and KD05 experiments for October, focusing on stratospheric levels from 100 to 0.1 hPa in the Southern Hemisphere. Overall, the distributions of the parameterized GWD in the two experiments are similar. The orographic GWD parameterized using the schemes of McFarlane (1987) [22]and Kim and Doyle (2005) [29], respectively, mainly occurs above 10 hPa in high latitudes. The gap in orographic GWD near 60° S is present in both experiments, consistent with the absence of topography in that latitude band. Compared with M87, significant orographic GWD enhancement in KD05 is found in the upper stratosphere over the Antarctic continent, along with moderate orographic GWD increase in the middle and lower stratosphere between 30° and 60° S where large mountains exist. On the other hand, the nonorographic GWD spreads out over a wider latitudinal region and generally increases with increasing height above 5 hPa in the mid to high latitudes with a latitudinal maximum around 70° S in the M87 and KD05 experiments. In contrast to the orographic GWs, the parameterized nonorographic GWs exert a strong drag on the zonal winds near 60° S. Although the nonorographic GWD in the two experiments is obtained from the same parameterization schemes, it weakens in the higher latitudes and strengthens in the lower latitudes above 1 hPa in KD05, probably because of the changes in the zonal-mean zonal winds induced by the enhanced orographic GWD. As a result, the changes in total GWD, i.e., the sum of orographic and nonorographic GWD, are dominated by orographic drag below about 1 hPa and nonorographic drag above 1 hPa.

Changes in GWD
It should be noted that the parameterized nonorographic GWD largely contributes to the total GWD in the stratosphere during the austral spring, both in the M87 and KD05 experiments. This differs from most climate models, in which the orographic GWD is much stronger than the nonorographic GWD. For example, McLandress (2012) found that the orographic GWD far outweighs the nonorographic GWD in the Canadian Middle Atmosphere Model and the Southern Hemisphere cold-pole bias was reduced when including an extra drag at 60° S in the orographic GWD scheme [11]. Uncertainties in the relative contribution of the nonorographic and orographic GWD to the total GWD come mainly from diverse constructions of nonorographic GWD parameterization in climate models. While the parameterization of GWs due to orography is relatively well understood, the GWs emitted by convection, fronts, and jets are often represented using uniform sources [33][34][35]. Source-related nonorographic GWD parameterizations may be more valid in modeling studies of the relative contribution of each GW during the southern SFW.  Figure 5 shows the latitude-height cross sections of orographic, nonorographic, and total GWD in the M87 and KD05 experiments for October, focusing on stratospheric levels from 100 to 0.1 hPa in the Southern Hemisphere. Overall, the distributions of the parameterized GWD in the two experiments are similar. The orographic GWD parameterized using the schemes of McFarlane (1987) [22] and Kim and Doyle (2005) [29], respectively, mainly occurs above 10 hPa in high latitudes. The gap in orographic GWD near 60 • S is present in both experiments, consistent with the absence of topography in that latitude band. Compared with M87, significant orographic GWD enhancement in KD05 is found in the upper stratosphere over the Antarctic continent, along with moderate orographic GWD increase in the middle and lower stratosphere between 30 • and 60 • S where large mountains exist. On the other hand, the nonorographic GWD spreads out over a wider latitudinal region and generally increases with increasing height above 5 hPa in the mid to high latitudes with a latitudinal maximum around 70 • S in the M87 and KD05 experiments. In contrast to the orographic GWs, the parameterized nonorographic GWs exert a strong drag on the zonal winds near 60 • S. Although the nonorographic GWD in the two experiments is obtained from the same parameterization schemes, it weakens in the higher latitudes and strengthens in the lower latitudes above 1 hPa in KD05, probably because of the changes in the zonal-mean zonal winds induced by the enhanced orographic GWD. As a result, the changes in total GWD, i.e., the sum of orographic and nonorographic GWD, are dominated by orographic drag below about 1 hPa and nonorographic drag above 1 hPa.

Changes in GWD
It should be noted that the parameterized nonorographic GWD largely contributes to the total GWD in the stratosphere during the austral spring, both in the M87 and KD05 experiments. This differs from most climate models, in which the orographic GWD is much stronger than the nonorographic GWD. For example, McLandress (2012) found that the orographic GWD far outweighs the nonorographic GWD in the Canadian Middle Atmosphere Model and the Southern Hemisphere cold-pole bias was reduced when including an extra drag at 60 • S in the orographic GWD scheme [11]. Uncertainties in the relative contribution of the nonorographic and orographic GWD to the total GWD come mainly from diverse constructions of nonorographic GWD parameterization in climate models. While the parameterization of GWs due to orography is relatively well understood, the GWs emitted by convection, fronts, and jets are often represented using uniform sources [33][34][35]. Source-related nonorographic GWD parameterizations may be more valid in modeling studies of the relative contribution of each GW during the southern SFW.

Resolved Wave Forcing Responses
The resolved wave flux from the troposphere was first examined, as shown by the 100-hPa meridional heat flux at Southern Hemisphere high latitudes for October ( Figure 6). It is clear that both waves 1 and 2 are improved significantly in KD05, indicating tropospheric wave forcing of the stratosphere is in closer agreement with JRA-55. Figure 7 shows latitude-height cross sections of monthly mean resolved wave forcing, i.e., divergence of the E-P flux, in M87 and KD05 for October throughout the southern stratosphere. As shown in JRA-55, the resolved wave drag is spread out over a wide latitudinal range, with a maximum around 1 hPa in the high latitudes (Figure 7a). Overall, the latitude-height distribution of the resolved wave forcing is reasonably reproduced in M87 (Figure 7b). The magnitude and the extent of the negative wave forcing in KD05 resembles that in M87 (Figures 7b and 7c). However, changes in the resolved wave drag in the stratosphere do occur due to changes in the parameterized GWD. Compared with M87, E-P flux divergence in KD05 becomes more negative in the lower stratosphere and more positive in the upper stratosphere. Furthermore, a notable feature in the difference of E-P flux divergence is that there is a tongue of Figure 5. Latitude-height cross sections in the stratosphere of (a,d) orographic, (b,e) nonorographic, and (c,f) total GWD (m s −1 day −1 ) in the M87 and KD05 experiments for October. Differences (KD05 minus M87) in gravity wave drag are shown in (g-i). Total GWD is the sum of the orographic and nonorographic GWD. Stippling denotes differences are statistically significant at the 95% confidence level.

Resolved Wave Forcing Responses
The resolved wave flux from the troposphere was first examined, as shown by the 100-hPa meridional heat flux at Southern Hemisphere high latitudes for October ( Figure 6). It is clear that both waves 1 and 2 are improved significantly in KD05, indicating tropospheric wave forcing of the stratosphere is in closer agreement with JRA-55. Figure 7 shows latitude-height cross sections of monthly mean resolved wave forcing, i.e., divergence of the E-P flux, in M87 and KD05 for October throughout the southern stratosphere. As shown in JRA-55, the resolved wave drag is spread out over a wide latitudinal range, with a maximum around 1 hPa in the high latitudes (Figure 7a). Overall, the latitude-height distribution of the resolved wave forcing is reasonably reproduced in M87 (Figure 7b). The magnitude and the extent of the negative wave forcing in KD05 resembles that in M87 (Figure 7b,c). However, changes in the resolved wave drag in the stratosphere do occur due to changes in the parameterized GWD. Compared with M87, E-P flux divergence in KD05 becomes more negative in the lower stratosphere and more positive in the upper stratosphere. Furthermore, a notable feature in the difference of E-P flux divergence is that there is a tongue of stronger wave forcing at 60 • S in KD05 (Figure 7d). These differences are generally significant at the 95% confidence level. Figure 8 shows the change in total wave forcing (GWD plus E-P flux divergence). Consistent with previous studies [36], changes in GWD are accompanied by changes in the E-P flux divergence of resolved waves such that partial compensation occurs between the two. One piece of evidence is that the negative changes in resolved wave drag above 1 hPa around 60 • S (Figure 7d) coincide with the positive changes in GWD in this region (Figure 5i). The compensation between different sources of mechanical forcing results in a much smoother latitudinal distribution of changes in total wave forcing than in GWD or E-P flux divergence individually. Figure 8 suggests that the reduced wave forcing in the upper stratosphere comes mainly from nonorographic GWD and resolved wave drag, whereas the enhanced wave forcing comes primarily from orographic GWD and resolved wave drag. A sharp decrease in total wave forcing is present around 66 • S. This corresponds to the abrupt decrease in the orographic GWD within grid cells of sharp topographic features surrounded by oceans, as seen in Figure 5g. The monthly mean annual cycle of the differences in the total wave forcing between the two experiments is presented in Figure 9. From August through November, the negative total wave forcing is strengthened in KD05, corresponding to the reduced cold bias as shown in Figure 1c.

Total Wave Forcing
Atmosphere 2020, 11, x FOR PEER REVIEW 11 of 17 stronger wave forcing at 60° S in KD05 (Figure 7d). These differences are generally significant at the 95% confidence level.   Figure 8 shows the change in total wave forcing (GWD plus E-P flux divergence). Consistent with previous studies [36], changes in GWD are accompanied by changes in the E-P flux divergence of resolved waves such that partial compensation occurs between the two. One piece of evidence is that the negative changes in resolved wave drag above 1 hPa around 60° S ( Figure 7d) coincide with the positive changes in GWD in this region (Figure 5i). The compensation between different sources of mechanical forcing results in a much smoother latitudinal distribution of changes in total wave forcing than in GWD or E-P flux divergence individually. Figure 8 suggests that the reduced wave forcing in the upper stratosphere comes mainly from nonorographic GWD and resolved wave drag, whereas the enhanced wave forcing comes primarily from orographic GWD and resolved wave drag. A sharp decrease in total wave forcing is present around 66° S. This corresponds to the abrupt decrease in the orographic GWD within grid cells of sharp topographic features surrounded by oceans, as seen in Figure 5g. The monthly mean annual cycle of the differences in the total wave forcing between the two experiments is presented in Figure 9. From August through November, the negative total wave forcing is strengthened in KD05, corresponding to the reduced cold bias as shown in Figure 1c.   . Difference (KD05 minus M87) in total wave forcing (m s −1 day −1 ) for October. Total wave drag is the sum of the resolved wave drag and total GWD. Stippling denotes differences are statistically significant at the 95% confidence level.

Polar Downwelling in the Southern Hemisphere
We consider next how the changes in wave forcing documented above affect downwelling (and hence temperature) at high latitudes of the Southern Hemisphere. Figure 10 shows the vertical component of the residual mean meridional circulation averaged over the Southern Hemisphere polar cap for M87 and KD05 simulations, and their difference. Compared to M87, the warming in early spring and the cooling in late spring (Figure 1c) are a consequence of enhanced and reduced downwelling in KD05, respectively (Figure 10d). The anomalous downwelling and upwelling are driven by the enhanced and reduced wave forcing, as seen in Figure 9. Note that anomalous

Polar Downwelling in the Southern Hemisphere
We consider next how the changes in wave forcing documented above affect downwelling (and hence temperature) at high latitudes of the Southern Hemisphere. Figure 10 shows the vertical component of the residual mean meridional circulation averaged over the Southern Hemisphere polar cap for M87 and KD05 simulations, and their difference. Compared to M87, the warming in early spring and the cooling in late spring (Figure 1c) are a consequence of enhanced and reduced downwelling in KD05, respectively (Figure 10d). The anomalous downwelling and upwelling are driven by the enhanced and reduced wave forcing, as seen in Figure 9. Note that anomalous upwelling locates at~0.5-3 hPa around June-August, which may be associated with reduced nonorographic GWD or resolved wave drag in the upper southern stratosphere.
Atmosphere 2020, 11, x FOR PEER REVIEW 14 of 17 upwelling locates at ~0.5-3 hPa around June-August, which may be associated with reduced nonorographic GWD or resolved wave drag in the upper southern stratosphere.

Summary and Conclusions
It is shown in this study that modifying the orographic GW parameterization results in an improved simulation of the SFW over Antarctica. A pair of experiments were conducted to

Summary and Conclusions
It is shown in this study that modifying the orographic GW parameterization results in an improved simulation of the SFW over Antarctica. A pair of experiments were conducted to compare two orographic GW schemes, namely, the McFarlane (1987) [22] scheme and the Kim and Doyle (2005) [29] scheme, in modeling the breakdown of the stratospheric polar vortex over Antarctica. The cold-pole bias is reduced by more than 2 K and the associated high wind bias is reduced by more than 4 m s −1 when the McFarlane (1987) [22] scheme is replaced by the Kim and Doyle (2005) [29] scheme. Correspondingly, the delayed SFW can be alleviated by 8 days. Further investigations on the wave forcing show that the improvements come from the accelerated downwelling branch of the BDC in the lower stratosphere, which in turn is mainly driven by the enhanced orographic GWD. This study suggests that modifying existing orographic GWD parameterizations would be a valid way of improving the SFW simulation over Antarctica.
Current representation of orographic GWD in climate models has been developed to a relatively complex level but still has uncertain formulation details. This study demonstrates that a more recent orographic GWD parameterization can go a long way in improving the cold-pole bias to a certain extent. However, comparing with statistics from the previous studies reveals that simply changing the parameterization schemes results in limited improvements compared to those that would be attained by adding a missing drag around 60 • S. The absent orographic GWD near 60 • S may be due to the simplified assumption that mountain waves only propagate vertically within the model grid column where they are triggered. The horizontal propagate of mountain waves is not taken into account in any traditional orographic GWD scheme. Recently, new orographic GW parameterizations considering horizontal wave propagation have been proposed and present significant impacts on the orographic GWD distribution [37,38]. Including the horizontal propagation effects of orographic GWs may lead to a much-improved cold-bias in climate models.
The general consensus to explain the delayed SFW over Antarctica is that climate models underestimate the GW forcing in the southern stratosphere, but the orographic or nonorographic origin of the insufficient GW drag is still controversial. Similar to orographic GWs, nonorographic GWs also act to drive the downwelling branch of the BDC. Observational studies have found that GWs over the ocean surrounding Antarctica are highly related to nonorographic sources [39][40][41][42]. Furthermore, the role of nonorographic GWs has been examined in modeling studies and it was found that parameterized GWs emitted by fronts and jets contribute significantly to the simulated Antarctic stratospheric vortex [10,12]. The effect of improving the representation of nonorographic GWs on the simulated southern SFW in climate models is beyond the scope of this paper and will be explored in future work.