the drag exerted by the earth surface on the atmospheric flow has long been recognized to play an important role in the general circulation of the atmosphere [ holton , 2004 ] .
this follows from considerations related to the angular momentum budget , which leads to a picture of angular momentum fluxes in the free troposphere balancing surface drag and explains the presence of surface easterlies in the tropics and surface westerlies in the midlatitudes . although it is clear that the magnitude of the surface stress controls the strength of the surface easterlies and westerlies , it has also , more surprisingly , been shown to affect the location of the surface westerly jet in the midlatitudes [ robinson , 1997 ; chen et al . , 2007 ] .
moreover , the drag exerted by topography is particularly important for the northern hemisphere ( nh ) winter extratropical circulation .
the stationary rossby ( or planetary ) waves forced by the flow over the himalayas and rockies were shown to explain the winter mean longitudinal distribution of 500 hpa heights in the nh midlatitudes [ charney and eliassen , 1949 ; held et al . ,
the drag exerted by orographically generated gravity waves has also been shown to affect the nh winter zonal mean flow both in the troposphere and the lower stratosphere since the early days of weather and climate prediction [ wallace et al . , 1983 ; palmer et al . , 1986 ;
these early studies emphasized the importance of parameterizing orographic gravity wave effects in general circulation models ( gcms ) in order to reproduce the observed mean patterns of the nh winter circulation .
later studies demonstrated that the orographic gravity wave drag affects the zonal mean flow by changing the resolved planetary waves [ mclandress and mcfarlane , 1993 ; shaw et al .
the parameterization of such orographic lowlevel blocking effects influences the representation of the nh winter circulation in both weather forecast and climate simulations [ lott and miller , 1997 ; zadra et al . , 2003 ; sandu et al .
albeit recognized as being important for both weather forecasts and climate simulations , the representation of drag processes associated with surface stress is a major source of uncertainty in gcms .
parameters controlling their strength are poorly constrained and often tuned to obtain the desired model behavior , i.e. , forecast skill [ sandu et al .
the uncertainty not only stems from the lack of observational constraints , but also from the fact that surface stress has received relatively little attention compared to other processes such as clouds , convection , or radiation . in the last couple of years ,
the interest in the surface momentum budget has increased , particularly in the framework of the world meteorological organization 's working group for numerical experimentation ( wgne ) and of a workshop on the  momentum budget and its role in weather and climate  organized in march 2015 at eynsham hall , united kingdom . in 2014 , wgne set up a model intercomparison exercise of subgrid ( or unresolved ) surface stress , the wgne drag project .
this exercise aims to compare the total subgrid surface stress ( tot ) and its partitioning among the parameterizations that most models use to represent it : the boundary layer scheme and the subgrid orography scheme ( section 2 will provide more details on the representation of subgrid surface stress in models ) .
monthly mean zonal and meridional subgrid surface stresses from the first 24 h of shortrange forecasts performed for january and july 2012 were gathered from 12 modeling groups [ zadra , 2015 ] .
since the models forecasts were starting from similar initial states , the angular momentum fluxes in the free atmosphere can be presumed to be similar and one would then expect similar total surface stresses .
first results of this exercise suggest that while the intermodel spread is relatively small over the oceans , tot is highly modeldependent over land [ zadra , 2015 ] .
as the subgrid , or parameterized stress characterizes the unresolved part of the flow , the intermodel spread found over land is due in part to differences in horizontal resolution .
nevertheless , the zonally averaged magnitude of tot can still differ by as much as 20% over land even for models of comparable resolution , such as the integrated forecast system ( ifs ) of the european centre for mediumrange weather forecasts ( ecmwf ) and the unified model ( um ) of the metoffice ( figure 1 ) .
these differences in the zonal average stress are mostly related to differences in stress over orography ( not shown ) . moreover , the partition of tot between the boundary layer and subgrid orography schemes is also modeldependent .
the um for example has almost twice as much subgrid orography stress and up to 50% less boundary layer stress than ifs in the nh ( figure 1 ) .
finally , the wgne drag project also revealed that the diurnal cycle of tot is quite different among the models ( not shown ) .
zonal averages of the magnitude of the total subgrid surface stress tot ( n / m ) and of its components for ifs ( black , scale on left ) and um ( blue , scale on left ) , from the results submitted for the wgne drag project for january 2012 , and relative difference between um and ifs ( red , scale on right ) . for ifs ,
the boundary layer stress ( left ) represents the sum of the stresses from the turbulence and turbulent orographic form drag schemes ; for um it represents the stress from the turbulence scheme which represents the orographic form drag implicitly through an orographic enhancement of the roughness ( see section 2 for details ) . for comparison with the results from the wgne drag project ,
the magnitude of the stress used to construct the zonal averages is computed in this plot only from monthly mean diurnal averages of the u and v components of the surface stress ( derived from the first 24 h of the daily 00utc forecasts ) .
the zonal averages include only land grid cells and were smoothed using a 5 running mean .
the ongoing wgne drag project emphasizes the importance of better constraining the representation of subgrid surface stress in models , especially in regions with orography .
it also highlights the need to better understand how the schemes contributing to surface stress , and their uncertainties , affect weather and climate predictions .
the present study looks into a particular aspect of this broader question , namely how some of these schemes impact the representation of the nh winter circulation on different time scales .
more precisely , our study aims to address two of the questions raised by the wgne drag project , namely ( i ) whether differences in tot of the order of those revealed by the wgne drag project matter for the representation of the nh winter circulation , and ( ii ) whether the partitioning of tot between the different schemes in a particular model is also important for predicting certain circulation aspects .
as most of the uncertainty in surface stress comes from regions with orography , we focus on the impacts of two of the schemes which represent orographic drag processes . to answer these questions
, we use the ifs to perform both 10 day weather forecasts and seasonal integrations , in which we mimic an increase in tot that is comparable to the difference between the ifs and the um shown in the wgne drag project .
this increase is simulated in two ways , namely by changing a poorly constrained parameter either in the boundary layer or in the subgrid orography scheme .
both changes concern the representation of subgrid orographic drag exerted on the flow at low levels .
the representation of the subgrid surface stress in the ifs is described in section 2 .
section 4 discusses the impacts of the changes to subgrid orographic drag on the representation of the nh winter circulation in short and mediumrange forecasts .
particular emphasis is put on whether these impacts are different depending on which of the parameterizations is modified .
the reason is that different schemes parameterize different dragrelated processes . in the ifs the boundary layer scheme represents stress associated with unresolved elements of the earth surface from horizontal scales less than 5 km , while the subgrid orography scheme ( sgo ) represents stress associated with horizontal scales between 5 km and the model resolution .
the 10 day weather forecasts and the seasonal integrations discussed in this study are performed at horizontal resolutions of approximately 32 , respectively , 75 km at the equator ( section 3 ) . at their turn , each of these two schemes represents several processes .
thus , the boundary layer scheme uses a turbulence parameterization ( turb hereafter ) to represent turbulent drag associated with surface elements such as land use or vegetation .
the turb scheme is based on monin obukhov similarity and uses prescribed roughness length tables to describe the characteristic height of these surface elements [ ifs documentation , 2015 , part iv ] .
the boundary layer scheme also includes a turbulent orographic form drag parameterization ( tofd hereafter ) to represent drag associated with subgrid orography elements with horizontal scales less than 5 km such as hills or individual mountains [ beljaars et al . , 2004 ] .
the tofd stress is formulated in terms of the variance of the terrain slope , which is derived from an empirical representation of the height spectrum in the scale range of 5 km to a few meters ( this spectrum is derived from a 1 km resolution elevation data set ) . the vertical distribution of the stress is prescribed on model levels according to the formulation proposed by wood et al . .
other models ( e.g. , the um of the metoffice ) do not have a special parameterization for turbulent orographic form drag , but represent it by artificially enhancing the surface roughness length over orography , so that the orographic form drag is implicitly represented by the turbulence parameterization .
the impact of either the tofd parameterization or the orographic enhancement of roughness is poorly documented in the literature [ wallace et al . , 1983 ] , although it has been shown to affect the nh wintertime circulation to some extent [ sandu et al .
the sgo scheme [ lott and miller , 1997 ] also represents different types of drag .
if the airflow has enough energy to pass over the socalled effective subgrid mountain height , orographic gravity waves are generated .
the subgrid stress associated with the propagation and breaking of these waves is represented by the gravity wave drag component of the sgo scheme ( gwd hereafter ) . at low levels , where the flow does not have enough energy
the stress induced by this process is described by the lowlevel blocking component of the sgo scheme ( block hereafter ) .
the depth over which the flow is blocked , and goes around the mountain , controls to a large extent the effect of the block scheme .
this blocking depth is equal to the subgrid mountain height minus the depth of the layer where the flow goes over the mountain ( and generates gravity waves ) .
the depth of this latter layer is defined by a critical froude number which makes both the gravity wave amplitude and the blocking depth stability dependent . for strong flow and/or weak stability ,
the gravity wave amplitude is large and the blocking depth is small , whereas for weak flow and very stable stratification the gravity wave amplitude is small and the blocking height large .
other models use similar schemes to represent gravity wave drag and lowlevel blocking effects [ e.g. , mcfarlane , 1987 ; webster et al . , 2003 ; brown and webster , 2004 ] , and some include an  orographic lift  component , which describes a force perpendicular to the flow that depends on the wind direction relative to the orientation of the subgrid orography [ lott , 1999 ] .
the contributions of the gwd and block schemes to the subgrid surface stress can not be disentangled in ifs , so they are presented combined hereafter as the sgo contribution to tot .
special diagnostics were added for the purpose of this study within the boundary layer scheme in order to separate the contributions of the turb and tofd schemes to tot .
figures 2 and 3 show how much each of these schemes contributes to tot in shortrange forecasts performed with the ifs for february 2014 ( the experimental setup is described in more detail in section 3 ) . in zonal average ,
the tofd and sgo schemes contribute fairly equally over land to the magnitude of tot ( figure 2 ) .
their contribution is comparable to that of the turb scheme at latitudes with a fair amount of orography ( e.g. , mid to high latitudes in the nh ) .
however , the maps shown in figure 3 indicate that the schemes act in different regions .
the tofd and sgo act as expected in regions with orography , while the turb scheme acts predominantly in regions with increased synoptic activity ( i.e. , storm tracks ) and its contribution to tot over orography is small .
these figures only illustrate 24 h averages of the monthly mean magnitude of the surface stress , ignoring thus differences in the diurnal cycle .
however , as it will be shown in the next section , the stress produced by the various schemes has also different diurnal cycles .
zonal averages of the magnitude of tot and of the contributions of the turb , tofd , and sgo schemes to tot ( n / m ) for the ctl experiment .
these zonal averages are computed from monthly means of daily surface stress magnitude ( derived from the first 24 h of the daily 00utc forecasts performed for february 2014 ) , include only land grid cells and were smoothed using a 5 running mean .
magnitude of the surface stress due to the turb , tofd , and sgo schemes ( n / m ) for the ctl experiment .
these fields represent monthly means of the daily surface stress magnitude , derived from the first 24 h of the daily 00utc forecasts performed for february 2014 .
the wgne drag project demonstrated that tot and its partition among the different schemes is largely modeldependent , particularly in regions with orography [ zadra , 2015 ] .
moreover , individual model results heavily depend on the assumptions and choices of parameter values made within the schemes contributing to tot . given that it is currently not possible to constrain either tot or its orographic component with observations on global , and even on regional scale , the assumptions and parameter choices often rely on a rather pragmatic optimization approach .
in previous work [ sandu et al . , 2013a , 2013b ] , we showed that the mediumrange forecast skill for the geopotential height during the nh winter can be improved by combining changes to the turb scheme concerning the degree of turbulent diffusion in stable layers with changes to certain parameters of the tofd and block schemes .
in the present study , we look in more detail at the individual impacts of the tofd and block schemes on the nh winter circulation .
these two schemes are likely the most poorly constrained and less documented in terms of impacts on the largescale circulation . given that they represent similar processes contributing to the surface stress ( orographic drag at low levels ) , these schemes also constitute good candidates for exploring the role of the partitioning of tot between different schemes .
the forecasts start from the ecmwf operational analysis and are performed with a vertical resolution of 137 levels at a linear spectral truncation of t639 corresponding to approximately 32 km at the equator .
second , seasonal integrations are performed with the same number of vertical levels , but at a lower spectral resolution , i.e. , t255 corresponding to approximately 75 km at the equator .
these runs are initialized with the ecmwf interim reanalysis erai [ dee et al . , 2011 ] .
they start on 1 august for every year between 1984 and 2014 and are integrated over a period of 7 months .
this setup allows to examine the climate of the model during the winter season , decemberjanuaryfebruary ( djf ) , and its response to changes in subgrid orographic drag .
the most recent operational ifs cycle ( cy41r1 operational since may 2015 ) is used for all experiments . for each of these two types of runs ,
we perform a control experiment with the reference model , an experiment with a larger magnitude of the stress resulting from the tofd scheme and an experiment with a larger magnitude of the stress resulting from the block scheme .
the 10 day forecast experiments will be denoted as ctl , high tofd ( htofd ) , and high block ( hblock ) hereafter , while the seasonal integrations will be denoted as ctlcl , htofdcl , and hblockcl .
the magnitude of the stress resulting from each scheme is enhanced by increasing the value of one of its poorly constrained constant factors , tofd for the tofd scheme and cd in the block scheme .
these factors scale the stress (  ) in the two schemes as follows : 
 ( 1)tofdtofdu|u|ftofd  blockcdu|u|fblock , where  is the air density , 
 u the horizontal wind vector , ftofd is a function of the standard deviation of the filtered subgrid orography , and fblock is function of stability , the characteristics of the subgrid orography ( i.e. , standard deviation , slope , and anisotropy ) and the direction of the flow with respect to the mountain ridges ( a detailed description of the two schemes is included in ifs documentation [ 2015 , chap . 3.4 and 4 , part iv ] ) .
the changes are such that they lead to a similar increase in the zonally averaged magnitude of tot ( defined as in the caption of figure 2 ) in the short range of the htofd and hblock forecasts .
furthermore , this increase is similar to the difference between the ifs and the um ( figure 1 ) . thus , the increase in zonal mean tot in the htofd and hblock experiments is of approximately 1015% in the nh midlatitudes ( figure 4 ) , while the difference between the ifs and um in the wgne drag project was of about 1520% ( figure 1 ) . also note that an increase in the stress from any of the schemes leads to a reduction in near surface wind , which further results in a small reduction of the stress produced by the other schemes ( figure 4 ) . in the htofd experiment for example ,
the tofd stress increases by approximately 60% , while the turb and sgo stress each decrease by 10% .
zonal averages of the magnitude of tot ( right ) and of the contributions of the turb , tofd , and sgo schemes to tot ( from left to right ) ( n / m ) for the ctl experiment ( black , scale on left ) .
the green and red curves , scale on right , show the relative change of tot and of its components in the experiments htofd and hblock with respect to the ctl .
the zonal averages are computed as in figure 2 . by design the zonally averaged magnitude of tot increases by a similar amount in the nh midlatitudes in the htofd and hblock experiments . however
, a zoom on the tot increase in the himalaya region reveals that the tofd scheme enhances the stress over a larger area than the block scheme , which seems to act more near the edges of the mountain chain ( figure 5 ) .
furthermore , figure 6 reveals that the diurnal cycle of the increase in tot is of opposite phase in the htofd and hblock experiments .
thus , the increase in tot is maximum for htofd and minimum for hblock between 6 and 12 utc ( about 1218 local time ) .
the boundary layer is mostly unstable ( or less stable ) and the nearsurface wind has a daily maximum due to strong turbulent mixing ( not shown ) .
since the tofd scheme responds to wind speed only , it will reflect the diurnal cycle of the surface wind .
however , the block scheme responds not only to wind speed , but is also strongly modulated by stability mainly through the depth of the blocked layer . when the wind speed is high and stability is weak ,
the flow the layer where the flow goes over the mountain is deep and the remaining layer where the flow is blocked is shallow .
monthly mean difference in the daily magnitude of tot ( n / m ) between the htofd ( left ) , respectively , hblock ( right ) and the ctl experiment .
time evolution of the monthly mean differences in the magnitude of tot ( n / m ) between the htofd ( green ) , respectively , hblock ( red ) and the ctl experiment during the first 5 days of the forecasts performed for february 2014 , integrated over the himalaya area showed in figure 5 .
full lines indicated 6 h averages , and dashed lines indicate diurnal averages of the difference in the magnitude of tot .
we examine the impacts of these two seemingly similar , yet different , ways of increasing tot in a particular model .
we seek to understand the consequences of these changes on the nh winter circulation , first on daily time scales and then at longer time ranges .
changes in surface pressure have traditionally been used to illustrate orographyrelated effects on the nh winter circulation since the first studies describing the parameterization of orographic gravity wave drag in gcms [ wallace et al . , 1983 ; palmer et al . , 1986 ;
the impact of enhancing the orographic drag at low levels in the htofd and hblock 10 day forecast experiments is therefore first assessed in terms of surface pressure changes at short and medium ranges ( 15 days ) .
the impacts on the representation of the nh circulation are then discussed in terms of zonal wind , temperature , and geopotential height .
the enhancement of tot , simulated by increasing either the tofd or the block stress , leads to changes in the mean surface pressure in a matter of hours ( figure 7 ) .
thus , pressure gradient changes across the regions with the largest mountain chains ( himalayas and the rocky mountains ) , as well as an increase in surface pressure over europe , can already be observed within the first 6 h of the forecasts . during
the first 24 h these changes amplify and extend to larger spatial scales in both experiments , although they are always stronger for the htofd forecasts .
mean change in surface pressure ( hpa ) in the experiments htofd ( left ) and hblock ( right ) with respect to the ctl , at lead time 6 ( top ) , 24 ( middle ) , and 120 ( bottom ) hours of the daily 00utc forecasts performed for february 2014 .
the shortrange response in surface pressure over europe is particularly intriguing because it does not clearly resemble the pattern of the pressure changes over north america and east asia .
this raises the question whether this is a local effect due to the alps or greenland or it is due to nonlocal effects caused by the changes in orographic drag over the rockies or the himalayas . more generally , the question is whether these changes in surface pressure in the short range are purely local orographic effects or are nonlocal effects due to impacts of the orographic drag on the zonal flow . to answer this question we performed additional daily 10 day forecasts for february 2014 .
the changes made in the original htofd and hblock experiments were this time applied only in selected regions of the globe . for each of the two schemes
, three new experiments were thus run in which the stress was enhanced over regions covering north america , greenland and europe , and east asia , illustrated by the green boxes in figure 8 .
mean change in surface pressure ( hpa ) in the experiments in which the tofd , respectively , the block stress , are enhanced only in the regions demarked by the dark green lines , with respect to the ctl , at lead time 24 of the daily 00utc forecasts performed for february 2014 .
the dark green hatches ( from 1 to 0 m / s , by 0.25 m / s ) indicate the deceleration of surface westerlies on the west side of the mountain ranges .
the changes in surface pressure observed after 24 h over north america and east asia in the new experiments ( figure 8 , top , bottom ) are virtually identical with those obtained in the original htofd and hblock experiments ( figure 7 ) .
the changes observed in the new experiments are also mainly due to a locally induced meridional pressure gradient ( figure 8 , middle ) .
the decrease in surface pressure seen east / south of the alps in the experiment where the tofd stress is enhanced over greenland and europe only ( figure 8 , middle ) is outweighed however in the original htofd experiment ( figure 7 ) by the increase in pressure induced by the changes in stress over east asia ( figure 8 , bottom ) . as suggested by zadra et al . , the local response in surface pressure over the largest mountain chains can be understood in terms of geostrophic balance .
the induced meridional pressure gradient across the rockies , alps , and , respectively , himalayas , corresponds to a westward geostrophic zonal wind anomaly , that is a deceleration of the midlatitude surface westerlies caused by the enhancement of the stress over orography .
a decrease in the zonal surface wind of 1 to 0.2 m s can be indeed seen on the west side of the mountain chains for all experiments ( green hatches in figure 8) . that changing the drag exerted by orography on the flow at low levels leads to a rapid response in surface pressure is hardly a surprise , given our current understanding ( zadra et al . , 2003 ) .
what is interesting instead is that although the changes made to the two schemes lead to a similar zonally averaged change in tot , the magnitude of the response in surface pressure is rather different ( figure 7 ) .
one possibility is that differences in the spatial distribution and the diurnal cycle of tot ( figures 5 and 6 ) lead to a different response in surface pressure .
the other possibility is that this response is damped to some extent in the hblock experiment due to the interactions between the gwd and block parameterizations .
indeed , the lowlevel blocking and gravity wave drag are both represented in the sgo scheme .
therefore , it is impossible to tell a priori whether a change of the stress at low levels , such as that simulated in the hblock experiment , does not also affect the orographic gravity waves .
if this was the case , the changes in surface pressure obtained in the hblock experiment could be due at least in part to changes in the propagation and breaking of the gravity waves . to investigate this possibility
, we repeated the ctl and hblock experiments while turning off the gwd part of the sgo scheme .
the difference in mean surface pressure between these two new experiments matches closely that between the original hblock and ctl experiments shown in figure 7 .
this suggests that the response obtained in the short range of the hblock forecasts is actually due to lowlevel blocking effects .
after 5 days , the changes to the tofd scheme lead to an overall increase of the surface pressure at high latitudes and a decrease at low latitudes .
the impacts on the zonal mean wind and temperature are also qualitatively similar , yet much stronger in the htofd experiment ( not shown ) .
these results corroborate the findings of zadra et al . , who discussed in detail the impacts of introducing the lowlevel blocking parameterization of lott and miller   in the gem model used at the canadian meteorological centre .
this study demonstrated that the introduction of the blocking scheme leads to a weakening of the midlatitude westerlies in the nh during wintertime , the flow being affected up to approximately 150 hpa by day 5 of the forecasts ( see their figure 6 ) .
it also induces a change in zonal mean temperature at latitudes higher than 40n characterized by a lowlevel warming and upper level cooling that also reaches 100 hpa by day 5 ( see their figure 9 ) .
relative difference of root mean square error ( rmse ) of 500 and 1000 hpa geopotential height , between the htofd ( black ) , respectively , hblock ( red ) , and the ctl forecasts , as a function of lead time .
when error bars do not cross the zero line , the performance of the respective experiment is significantly worse / better ( 95% interval ) than the ctl .
for both the experiments and the ctl forecasts the rmse was computed with respect to the analysis from which the forecasts started .
a positive difference in rmse indicates a deterioration of the model performance in the experiment with respect to the ctl .
the effects obtained when increasing the tofd stress are stronger in all circulation aspects examined so far .
the impact on the forecast performance is hence bigger in the htofd experiment ( figure 9 ) .
what is surprising however is that the change made to the tofd stress considerably degrades the ability of the model to reproduce the nh circulation in the entire troposphere , while that made to the block stress leads to a slight improvement of the model performance .
as illustrated by the increase , respectively , the reduction of the root mean square error of the geopotential height at 1000 and 500 hpa in figure 9 , these impacts are significant up to day 5 for the htofd , and day 10 for the hblock experiment .
the enhancement of tot , simulated by increasing either the tofd or the block stress , leads to changes in the mean surface pressure in a matter of hours ( figure 7 ) .
thus , pressure gradient changes across the regions with the largest mountain chains ( himalayas and the rocky mountains ) , as well as an increase in surface pressure over europe , can already be observed within the first 6 h of the forecasts . during
the first 24 h these changes amplify and extend to larger spatial scales in both experiments , although they are always stronger for the htofd forecasts .
mean change in surface pressure ( hpa ) in the experiments htofd ( left ) and hblock ( right ) with respect to the ctl , at lead time 6 ( top ) , 24 ( middle ) , and 120 ( bottom ) hours of the daily 00utc forecasts performed for february 2014 .
the shortrange response in surface pressure over europe is particularly intriguing because it does not clearly resemble the pattern of the pressure changes over north america and east asia .
this raises the question whether this is a local effect due to the alps or greenland or it is due to nonlocal effects caused by the changes in orographic drag over the rockies or the himalayas . more generally , the question is whether these changes in surface pressure in the short range are purely local orographic effects or are nonlocal effects due to impacts of the orographic drag on the zonal flow . to answer this question we performed additional daily 10 day forecasts for february 2014 .
the changes made in the original htofd and hblock experiments were this time applied only in selected regions of the globe . for each of the two schemes
, three new experiments were thus run in which the stress was enhanced over regions covering north america , greenland and europe , and east asia , illustrated by the green boxes in figure 8 .
mean change in surface pressure ( hpa ) in the experiments in which the tofd , respectively , the block stress , are enhanced only in the regions demarked by the dark green lines , with respect to the ctl , at lead time 24 of the daily 00utc forecasts performed for february 2014 .
the dark green hatches ( from 1 to 0 m / s , by 0.25 m / s ) indicate the deceleration of surface westerlies on the west side of the mountain ranges .
the changes in surface pressure observed after 24 h over north america and east asia in the new experiments ( figure 8 , top , bottom ) are virtually identical with those obtained in the original htofd and hblock experiments ( figure 7 ) .
the changes observed in the new experiments are also mainly due to a locally induced meridional pressure gradient ( figure 8 , middle ) . the decrease in surface pressure seen east / south of the alps in the experiment where the tofd stress is enhanced over greenland and europe only ( figure 8 , middle ) is outweighed however in the original htofd experiment ( figure 7 ) by the increase in pressure induced by the changes in stress over east asia ( figure 8 , bottom ) .
as suggested by zadra et al . , the local response in surface pressure over the largest mountain chains can be understood in terms of geostrophic balance .
the induced meridional pressure gradient across the rockies , alps , and , respectively , himalayas , corresponds to a westward geostrophic zonal wind anomaly , that is a deceleration of the midlatitude surface westerlies caused by the enhancement of the stress over orography .
a decrease in the zonal surface wind of 1 to 0.2 m s can be indeed seen on the west side of the mountain chains for all experiments ( green hatches in figure 8) . that changing the drag exerted by orography on the flow at low levels leads to a rapid response in surface pressure is hardly a surprise , given our current understanding ( zadra et al . ,
what is interesting instead is that although the changes made to the two schemes lead to a similar zonally averaged change in tot , the magnitude of the response in surface pressure is rather different ( figure 7 ) .
one possibility is that differences in the spatial distribution and the diurnal cycle of tot ( figures 5 and 6 ) lead to a different response in surface pressure .
the other possibility is that this response is damped to some extent in the hblock experiment due to the interactions between the gwd and block parameterizations .
indeed , the lowlevel blocking and gravity wave drag are both represented in the sgo scheme .
therefore , it is impossible to tell a priori whether a change of the stress at low levels , such as that simulated in the hblock experiment , does not also affect the orographic gravity waves . if this was the case , the changes in surface pressure obtained in the hblock experiment could be due at least in part to changes in the propagation and breaking of the gravity waves . to investigate this possibility
, we repeated the ctl and hblock experiments while turning off the gwd part of the sgo scheme .
the difference in mean surface pressure between these two new experiments matches closely that between the original hblock and ctl experiments shown in figure 7 .
this suggests that the response obtained in the short range of the hblock forecasts is actually due to lowlevel blocking effects .
after 5 days , the changes to the tofd scheme lead to an overall increase of the surface pressure at high latitudes and a decrease at low latitudes .
the impacts on the zonal mean wind and temperature are also qualitatively similar , yet much stronger in the htofd experiment ( not shown ) .
these results corroborate the findings of zadra et al . , who discussed in detail the impacts of introducing the lowlevel blocking parameterization of lott and miller   in the gem model used at the canadian meteorological centre .
this study demonstrated that the introduction of the blocking scheme leads to a weakening of the midlatitude westerlies in the nh during wintertime , the flow being affected up to approximately 150 hpa by day 5 of the forecasts ( see their figure 6 ) .
it also induces a change in zonal mean temperature at latitudes higher than 40n characterized by a lowlevel warming and upper level cooling that also reaches 100 hpa by day 5 ( see their figure 9 ) .
relative difference of root mean square error ( rmse ) of 500 and 1000 hpa geopotential height , between the htofd ( black ) , respectively , hblock ( red ) , and the ctl forecasts , as a function of lead time .
when error bars do not cross the zero line , the performance of the respective experiment is significantly worse / better ( 95% interval ) than the ctl .
for both the experiments and the ctl forecasts the rmse was computed with respect to the analysis from which the forecasts started .
a positive difference in rmse indicates a deterioration of the model performance in the experiment with respect to the ctl .
the effects obtained when increasing the tofd stress are stronger in all circulation aspects examined so far .
the impact on the forecast performance is hence bigger in the htofd experiment ( figure 9 ) .
what is surprising however is that the change made to the tofd stress considerably degrades the ability of the model to reproduce the nh circulation in the entire troposphere , while that made to the block stress leads to a slight improvement of the model performance .
as illustrated by the increase , respectively , the reduction of the root mean square error of the geopotential height at 1000 and 500 hpa in figure 9 , these impacts are significant up to day 5 for the htofd , and day 10 for the hblock experiment .
in the previous section we showed that the partitioning of tot between the tofd and block schemes matters for the representation of the nh winter circulation in 10 day weather forecasts . in this section we investigate if this partition is also important at longer time ranges .
we do that by analyzing the mean changes in circulation over the 30 winter ( djf ) seasons simulated in the seasonal integrations described in section 3 . given that the runs start on 1 august
, the examined lead time is between 4 and 7 months , which is similar to lead times of relevance for seasonal forecast integrations . in both the htofdcl and hblockcl experiments ,
the mean changes in surface pressure over these 30 djf seasons ( figure 10 , top ) resemble those after 5 days of integration in the htofd experiment ( figure 7 ) in that they show an increase in surface pressure at high latitudes and a decrease in the midlatitudes . however , this change in the meridional pressure gradient is particularly strong over the north atlantic , where it resembles a negative phase of the north atlantic oscillation .
although these changes seem at first sight stronger in the htofdcl experiment , the regions in which they are statistically significant ( according to a 95% significance level computed with a regular paired t test ) are less wide spread than in the hblockcl experiment .
mean change in surface pressure ( hpa ) ( top row ) and zonal wind at 200 hpa ( m / s ) ( bottom row ) in htofdcl and hblockcl with respect to the ctlcl experiment over the 30 djf seasons ( 19842014 ) .
the horizontal white pattern ( gridded areas ) indicate 90 ( 95)% statistical significance according to a paired t test .
the effects of orographic lowlevel drag on the upper tropospheric zonal winds are illustrated in figure 10 ( bottom ) .
the flow changes suggest that the increase in tot affects the subtropical jets in both experiments , especially over the north atlantic and north pacific , where the southward displacement of the storm tracks is most pronounced .
interestingly , there are important differences between the responses over europe , where climate models can show large biases in storm track location [ zappa et al . ,
changes in the zonal mean zonal wind are also seen in the stratosphere in both experiments ( figure 11 ) .
albeit weaker , the impacts on stratospheric winds are similar to those obtained by sigmond and scinocca   when enhancing the orographic gravity wave momentum flux ( see their figure 3 ) .
their study showed that increased orographic gravity wave momentum flux produces a basic state with weaker midlatitude winds in the wintertime lower stratosphere .
furthermore , it demonstrated that the deceleration of the lower stratospheric winds is related to changes in the resolved waves .
thus , an increased gravity wave drag strengthens the barrier to the equatorward propagation of planetary waves occurring in the midlatitude lower stratosphere .
consequently , the equatorward eliassenpalm ( ep ) flux is reduced and the planetary wave driving of the mid to highlatitude stratosphere is increased , leading to a deceleration of the polar stratospheric winds . by downward
control arguments changes in the stratospheric circulation eventually affect the tropospheric circulation [ haynes et al . , 1991 ] .
mean change in zonal mean zonal wind ( m / s , contours ) in htofdcl ( left ) and hblockcl ( right ) with respect to the ctlcl experiment over the 30 djf seasons ( 19842014 ) .
the horizontal white pattern ( gridded areas ) indicate 90 ( 95)% statistical significance according to a paired t test .
our experiments show a deceleration of the zonal winds on the poleward flank of the stratospheric polar vortex and a slight acceleration on its equatorward flank ( figure 11 ) .
these changes however appear to be statistically significant only in the hblockcl experiment . in the htofdcl experiment
, the only significant changes are found in the lower troposphere and midlatitude stratosphere . in order to understand
if these changes in stratospheric winds can also be explained by changes in the resolved wave driving , we performed an ep flux analysis similar to that of sigmond and scinocca .
a transformed eulerian mean ( tem ) momentum balance was used to analyze the forcing of the zonally averaged flow by transient plus stationary eddies , following andrews et al . ; and the components of the ep flux , which is a direct measure of the total forcing of the zonalmean state by the eddies [ edmon et al . , 1980 ] ,
the resulting ep flux arrows for the ctlcl experiment are displayed in figure 12 ( left ) .
( left ) epflux ( arrows ) and zonal mean zonal wind ( contours , m / s ) in the ctlcl experiment for the 30 djf seasons ( 19842014 ) .
( center / right ) change in epflux ( arrows ) in htofdcl , respectively , hblockcl with respect to the ctlcl experiment over the 30 djf seasons .
for clearer visualization the ep flux arrows have been scaled by the relative ranges of the two axes of the plot ( 3.14 radians by 10 pa ) .
the visibility of vectors in the stratosphere has been further enhanced by using a scaling factor proportional to the squareroot of pressure and by applying a magnifying factor ( of 3 ) above 100 hpa .
a budget for resolved wave driving is presented for the stratospheric box between 100 and 10 hpa and between 45n and 90n .
red numbers across the box represent differences of ep fluxes integrated over the box boundaries between the htofdcl / hblockcl and the ctlcl experiment , while the red numbers in the box represent the differences in time and areamean momentum deposition , or ep flux divergence , associated with the resolved waves ( all wave numbers ) within the box ( 
 104  kg  ms4 ) .
we then integrated the ep flux following kushner and polvani [ 2004 , equation ( 7 ) ] over a box bounded by 100 and 10 hpa and 4590n ( green lines in figure 12 , middle , right ) .
the red numbers on the sides of the box represent the difference in the integrated ep flux through the respective boundary , between each experiment and the ctlcl .
the values in figure 12 suggest that for both experiments the amount of upward ep flux from the troposphere into the stratosphere across the 100 hpa boundary increases , and so does the equatorward ep flux across the 45n boundary . in the hblockcl experiment the upward ep flux through the 10 hpa surface decreases slightly , while in the htofdcl it increases by an amount comparable to the increase across the 45n boundary .
overall , the hblockcl experiment shows thus an increased momentum deposition , or ep flux convergence , in the polar stratosphere , while in the htofdcl experiment this resolved wave driving changes little compared to the ctlcl experiment .
this larger momentum deposition explains why the deceleration of winds at latitudes north of 60 n is stronger ( and statistically more robust ) in the hblockcl experiment than in the htofdcl one ( figure 11 ) .
the differences in ep flux also explain the deceleration of the zonal winds above 30 hpa observed around 30n in both experiments .
the stronger deceleration observed in the htofdcl experiment appears to be associated with a stronger ep flux convergence at 2030 hpa between 30n and 40n .
the analysis of the differences in various aspects of the nh winter circulation over the 30 djf seasons suggests that the partitioning of the subgrid orographic lowlevel drag between the tofd and block schemes is also important at longer time scales . while both schemes affect the circulation in the entire atmosphere , the details of these effects depend on which scheme represents the drag .
the block scheme leads to effects somewhat similar , albeit weaker , than those of the gwd scheme discussed in sigmond et al .   and sigmond and scinocca .
in contrast to those studies , the increased highlatitude stratospheric wave drag in hblockcl is caused by more upward wave activity entering the region ( figure 12 ) , rather than by less equatorward wave activity leaving the region [ see sigmond and scinocca , 2010 , figure 6 ] ; in other words , the changes are induced from below rather than from within the stratosphere .
the turbulent orographic form drag on the other hand affects more strongly the lower troposphere and only partly the lower stratosphere .
finally , note that the main contribution to the ep flux in the stratosphere is due to wave numbers 15 , while the main contribution to the tropospheric ep flux is by wave numbers greater than 5 ( not shown ) .
therefore , the orographic form drag mainly affects the synoptic eddies while the lowlevel blocking affects the largescale waves .
the representation of surface stress is important for accurately simulating the largescale circulation in both weather forecasts and climate predictions .
the first results of the ongoing wgne drag project revealed that stateoftheart weather forecast models differ significantly in terms of the magnitude of the total subgrid surface stress , especially over land and more particularly over orography .
furthermore , models also differ in the partitioning of this stress among the different parameterization schemes . in the present study we tried to answer two questions raised by this intercomparison .
we investigated the impact of changing tot in a particular model by an amount comparable to the intermodel spread found in the wgne drag project .
the discussion focuses on changes in the representation of the nh wintertime circulation in 10 day weather forecasts and seasonal integrations .
furthermore , we examined the impact of varying the partitioning of this increase in tot between two schemes representing subgrid orographic drag , the turbulent orographic form drag and the lowlevel blocking schemes .
these schemes have been chosen as the representation of orographic drag processes is still uncertain , and their impacts have been only poorly documented in the literature . the experiments designed to address these questions consisted in increasing the zonal mean tot by 1015% in the nh midlatitudes , via an increase in either the tofd or the block contributions . such an increase in tot was shown to significantly affect various aspects of the nh winter circulation both in 10 day forecasts and in longer integrations , irrespectively of which scheme is modified .
thus , changes in surface pressure are seen in all experiments and at all examined time scales , ranging from a few hours to a few months .
a meridional surface pressure gradient across the large mountain chains is visible after a few hours of integration .
these changes are induced through geostrophic balance as a response to the deceleration of the midlatitude surface westerlies caused by the enhancement of orographic stress . at longer time ranges
the enhanced orographic stress also results in changes in the zonal wind , temperature , and geopotential at all time scales .
the amplitude and robustness of the changes in surface pressure , winds and temperature depend however on the particular scheme which is modified .
interestingly , in 10 day forecasts the increase in tofd stress has a much stronger impact on the circulation than that in block stress , while the impact is of comparable magnitude for the seasonal integrations .
more surprisingly , although both simulate an increase in stress over orography the changes to the tofd and block schemes have opposite impacts in terms of mediumrange forecast skill .
the tofd changes lead to a considerable degradation of the geopotential height skill in the troposphere ( up to 200 hpa ) up to day 5 , while the block changes slightly improve the performance throughout the entire forecast range .
the seasonal integrations suggest that the two schemes have also different effects in terms of the model climate for the djf season .
the increase in block stress results in a significant deceleration of the zonal wind on the poleward flank of the stratospheric polar vortex by changing the resolved wave driving of the mid to highlatitude stratosphere .
this was confirmed by an ep flux analysis similar to that performed by sigmond and scinocca .
the changes to the tofd scheme on the other hand mainly affect the lower troposphere and the midlatitude stratosphere . in summary , these results suggest that the representation of the nh winter circulation is sensitive not only to changes in the total subgrid ( orographic ) stress , but also to the partitioning of this stress among different schemes .
the reasons for the different circulation response to the partitioning of the total surface stress remain to be understood .
the different impacts of the two schemes discussed here could be related to inherent differences in their formulations , i.e. , stability versus non stability dependence , the diurnal cycle and the horizontal scales on which they act .
 abstracta recent intercomparison exercise proposed by the working group for numerical experimentation ( wgne ) revealed that the parameterized , or unresolved , surface stress in weather forecast models is highly modeldependent , especially over orography . 
 models of comparable resolution differ over land by as much as 20% in zonal mean total subgrid surface stress ( tot ) . 
 the way tot is partitioned between the different parameterizations is also modeldependent . in this study 
 , we simulated in a particular model an increase in tot comparable with the spread found in the wgne intercomparison . 
 this increase was simulated in two ways , namely by increasing independently the contributions to tot of the turbulent orographic form drag scheme ( tofd ) and of the orographic lowlevel blocking scheme ( block ) . 
 increasing the parameterized orographic drag leads to significant changes in surface pressure , zonal wind and temperature in the northern hemisphere during winter both in 10 day weather forecasts and in seasonal integrations . 
 however , the magnitude of these changes in circulation strongly depends on which scheme is modified . in 10 day forecasts , 
 stronger changes are found when the tofd stress is increased , while on seasonal time scales the effects are of comparable magnitude , although different in detail . at these time scales , 
 the block scheme affects the lower stratosphere winds through changes in the resolved planetary waves which are associated with surface impacts , while the tofd effects are mostly limited to the lower troposphere . 
 the partitioning of tot between the two schemes appears to play an important role at all time scales . 
Introduction
Representation of Subgrid Surface Stress in the IFS
Experimentation
Impacts in Short and MediumRange Weather Forecasts
ShortRange Response in Surface Pressure
Impacts on the NH Circulation in MediumRange Forecasts
Impacts at Longer Time Range
Conclusions