scholarly journals The unusually long cold spell and the snowstorm Filomena in Spain in January 2021

2022 ◽  
Author(s):  
Philipp Zschenderlein ◽  
Heini Wernli

Abstract. In early January 2021, Spain was affected by two extreme events – an unusually long cold spell and a heavy snowfall event associated with extratropical cyclone Filomena. For example, up to 50 cm of snow fell in Madrid and the surrounding areas in 4 days. Already during 9 days prior to the snowfall event, anomalously cold temperatures at 850 hPa and night frosts prevailed over large parts of Spain. During this period, anomalously cold and dry air was transported towards Spain from central Europe and even from the Barents Sea. The storm Filomena, which was responsible for major parts of the snowfall event, developed from a precursor low-pressure system over the central North Atlantic. Filomena intensified due to interaction with an upper-level potential vorticity (PV) trough, which was the result of anticyclonic wave breaking over Europe. In turn, this wave breaking was related to an intense surface anticyclone and upper-level ridge, whose formation was strongly influenced by a warm conveyor belt outflow of a cyclone off the coast of Newfoundland. The most intense snowfall occurred on 09 January and was associated with a sharp air mass boundary with an equivalent potential temperature difference at 850 hPa across Spain exceeding 20 K. Overall, the combination of pre-existing cold surface temperatures, the optimal position of the air mass boundary, and the dynamical forcing for ascent induced by Filomena and its associated upper-level trough were all essential – and in parts physically independent – ingredients for this extreme snowfall event to occur.

2021 ◽  
Author(s):  
Philipp Zschenderlein ◽  
Heini Wernli

<p>In January 2021, large parts of Spain were affected by an unusually long cold spell and exceptional snowfall associated with the winter storm Filomena. According to the Spanish weather service AEMET, snow heights of nearly 50 cm were registered in and around Madrid. During the days after Filomena, record-breaking low temperatures were measured at many stations.</p><p>Already during the days before the arrival of storm Filomena, anomalously cold temperatures at 850 hPa and night frosts at the surface prevailed over large parts of Spain. During these days in early January, the air flow towards Spain was predominantly northeasterly and advected cold air masses from Central Europe, as revealed by backward trajectories that were initialised near the surface over Spain. The land surface progressively cooled down during the days prior to the heavy snowfall, which then prevented the snow from melting when reaching the surface. Therefore, this cold spell preconditioning seems to be very important for the extreme consequences of the snowfall event.</p><p>The storm Filomena affected Spain between 8 and 10 January. It developed from a precursor low-pressure system between the Azores and Madeira. The precursor low-pressure system itself developed on 2 January 2021 between the northeastern US and Nova Scotia, rapidly intensified along a potential vorticity (PV) streamer and propagated southeastwards. Between 4 and 6 January, the cyclone, now located near the Azores, was associated with a PV cut-off and eventually decayed into multiple centres. Out of this decaying low-pressure system, Filomena developed and reached Spain on 8 January.</p><p>The most intense snowfall occurred on 9 January and affected large parts of Spain, except for southwestern Spain, where temperatures were too high and all precipitation fell as rain. Filomena was associated with an intense air mass boundary, with dry and cold air in the north and warm and humid air in the south. Equivalent potential temperature differences at 850 hPa across Spain exceeded 20 K. Along the warm frontal part of this air mass boundary, strong ascending airstreams, intensified by the dynamics of Filomena, led to cloud formation. Due to the unusually cold lowermost troposphere, snow was not melting before reaching the land surface, and the surface snow layer could therefore easily grow.</p><p>Overall, the combination of the already cold temperatures near the surface, the optimal position of the air mass boundary, and the dynamical forcing for ascent at this intense baroclinic zone associated with Filomena were essential ingredients for this extreme snow fall event to occur.</p>


2021 ◽  
Author(s):  
Antoine Hochet ◽  
Rémi Tailleux ◽  
Till Kuhlbrodt ◽  
David Ferreira

AbstractThe representation of ocean heat uptake in Simple Climate Models used for policy advice on climate change mitigation strategies is often based on variants of the one-dimensional Vertical Advection/Diffusion equation (VAD) for some averaged form of potential temperature. In such models, the effective advection and turbulent diffusion are usually tuned to emulate the behaviour of a given target climate model. However, because the statistical nature of such a “behavioural” calibration usually obscures the exact dependence of the effective diffusion and advection on the actual physical processes responsible for ocean heat uptake, it is difficult to understand its limitations and how to go about improving VADs. This paper proposes a physical calibration of the VAD that aims to provide explicit traceability of effective diffusion and advection to the processes responsible for ocean heat uptake. This construction relies on the coarse-graining of the full three-dimensional advection diffusion for potential temperature using potential temperature coordinates. The main advantage of this formulation is that the temporal evolution of the reference temperature profile is entirely due to the competition between effective diffusivity that is always positive definite, and the water mass transformation taking place at the surface, as in classical water mass analyses literature. These quantities are evaluated in numerical simulations of present day climate and global warming experiments. In this framework, the heat uptake in the global warming experiment is attributed to the increase of surface heat flux at low latitudes, its decrease at high latitudes and to the redistribution of heat toward cold temperatures made by diffusive flux.


2007 ◽  
Vol 135 (6) ◽  
pp. 2168-2184 ◽  
Author(s):  
Gregory L. West ◽  
W. James Steenburgh ◽  
William Y. Y. Cheng

Abstract Spurious grid-scale precipitation (SGSP) occurs in many mesoscale numerical weather prediction models when the simulated atmosphere becomes convectively unstable and the convective parameterization fails to relieve the instability. Case studies presented in this paper illustrate that SGSP events are also found in the North American Regional Reanalysis (NARR) and are accompanied by excessive maxima in grid-scale precipitation, vertical velocity, moisture variables (e.g., relative humidity and precipitable water), mid- and upper-level equivalent potential temperature, and mid- and upper-level absolute vorticity. SGSP events in environments favorable for high-based convection can also feature low-level cold pools and sea level pressure maxima. Prior to 2003, retrospectively generated NARR analyses feature an average of approximately 370 SGSP events annually. Beginning in 2003, however, NARR analyses are generated in near–real time by the Regional Climate Data Assimilation System (R-CDAS), which is identical to the retrospective NARR analysis system except for the input precipitation and ice cover datasets. Analyses produced by the R-CDAS feature a substantially larger number of SGSP events with more than 4000 occurring in the original 2003 analyses. An oceanic precipitation data processing error, which resulted in a reprocessing of NARR analyses from 2003 to 2005, only partially explains this increase since the reprocessed analyses still produce approximately 2000 SGSP events annually. These results suggest that many NARR SGSP events are not produced by shortcomings in the underlying Eta Model, but by the specification of anomalous latent heating when there is a strong mismatch between modeled and assimilated precipitation. NARR users should ensure that they are using the reprocessed NARR analyses from 2003 to 2005 and consider the possible influence of SGSP on their findings, particularly after the transition to the R-CDAS.


2009 ◽  
Vol 66 (2) ◽  
pp. 495-507 ◽  
Author(s):  
Lawrence Coy ◽  
Stephen Eckermann ◽  
Karl Hoppel

Abstract The major stratospheric sudden warming (SSW) of January 2006 is examined using meteorological fields from Goddard Earth Observing System version 4 (GEOS-4) analyses and forecast fields from the Navy Operational Global Atmospheric Prediction System–Advanced Level Physics, High Altitude (NOGAPS-ALPHA). The study focuses on the upper tropospheric forcing that led to the major SSW and the vertical structure of the subtropic wave breaking near 10 hPa that moved low tropical values of potential vorticity (PV) to the pole. Results show that an eastward-propagating upper tropospheric ridge over the North Atlantic with its associated cold temperature perturbations (as manifested by high 360-K potential temperature surface perturbations) and large positive local values of meridional heat flux directly forced a change in the stratospheric polar vortex, leading to the stratospheric subtropical wave breaking and warming. Results also show that the anticyclonic development, initiated by the subtropical wave breaking and associated with the poleward advection of the low PV values, occurred over a limited altitude range of approximately 6–10 km. The authors also show that the poleward advection of this localized low-PV anomaly was associated with changes in the Eliassen–Palm (EP) flux from equatorward to poleward, suggesting an important role for Rossby wave reflection in the SSW of January 2006. Similar upper tropospheric forcing and subtropical wave breaking were found to occur prior to the major SSW of January 2003.


2021 ◽  
pp. 17-28
Author(s):  
A. V. Gochakov ◽  
◽  
O. Yu. Antokhina ◽  
V. N. Krupchatnikov ◽  
Yu. V. Martynova ◽  
...  

Many large-scale dynamic phenomena in the Earth’s atmosphere are associated with the processes of propagation and breaking of Rossby waves. A new method for identifying the Rossby wave breaking (RWB) is proposed. It is based on the detection of breakings centers by analyzing the shape of the contours of potential vorticity or temperature on quasimaterial surfaces: isentropic and iserthelic (surfaces of constant Ertel potential vorticity (PV)), with further RWB center clustering to larger regions. The method is applied to the set of constant PV levels (0.3 to 9.8 PVU with a step of 0.5 PVU) at the level of potential temperature of 350 K for 12:00 UTC. The ERA-Interim reanalysis data from 1979 to 2019 are used for the method development. The type of RWB (cyclonic/anticyclonic), its area and center are determined by analyzing the vortex geometry at each PV level for every day. The RWBs obtained at this stage are designated as elementary breakings. Density-Based Spatial Clustering of Applications with Noise algorithm (DBSCAN) was applied to all elementary breakings for each month. As a result, a graphic dataset describing locations and dynamics of RWBs for every month from 1979 to 2019 is formed. The RWB frequency is also evaluated for each longitude, taking into account the duration of each RWB and the number of levels involved, as well as the anomalies of these parameters.


2021 ◽  
Author(s):  
Antoine Hochet ◽  
Remi Tailleux ◽  
Till Kuhlbrodt ◽  
David Ferreira

Abstract The representation of ocean heat uptake in Simple Climate Models used for policy advice on climate change mitigation strategies is often based on variants of the one-dimensional Vertical Advection/Diffusion equation (VAD) for some averaged form of potential temperature. In such models, the effective advection and turbulent diffusion are usually tuned to emulate the behaviour of a given target climate model. However, because the statistical nature of such a \behavioural" calibration usually obscures the exact dependence of the effective diffusion and advection on the actual physical processes responsible for ocean heat uptake, it is difficult to understand its limitations and how to go about improving VAD. This paper proposes a physical calibration of the VAD that aims to provide explicit traceability of effective diffusion and advection to the processes responsible for ocean heat uptake. This construction relies on the coarse-graining of the full three-dimensional advection diffusion for potential temperature using potential temperature coordinates. The main advantage of this formulation is that the temporal evolution of the reference temperature profile is entirely due to the competition between effective diffusivity that is always positive definite, and the water mass transformation taking place at the surface, as in classical water mass analyses literature. These quantities are evaluated in numerical simulations of present day climate and global warming experiments. In this framework, the heat uptake in the global warming experiment is attributed to the increase of surface heat flux at low latitudes, its decrease at high latitudes and to the redistribution of heat toward cold temperatures made by diffusive flux.


2021 ◽  
Author(s):  
Toru Terao ◽  
Fumie Murata ◽  
Yusuke Yamane ◽  
Masashi Kiguchi ◽  
Azusa Fukushima ◽  
...  

<p>The Asian summer monsoon system is the strongest monsoon circulation on the Earth. A huge reversal of meridional temperature gradient develops over the area covering the hemispheric region due to strong diabatic heating associated with convective activities. Vigorous conventions reach the upper troposphere providing a great amount of high potential temperature airmass. This high potential temperature air mass originates from the high equivalent potential temperature airmass accumulated in the lower troposphere over the Asian monsoon region. The highest potential temperature tropospheric air mass is observed only over the Asian summer monsoon region. To get a total view of the Asian summer monsoon circulation system, we focused on the mass budget of the upper-tropospheric air mass with a potential temperature between 355K to 370K. The non-conservative change of the air mass corresponds with the diabatic heating due to the convective activities, and the diabatic cooling due to the radiative process. To analyze the radiative cooling process that takes place in the upper troposphere, we utilized hourly GSMaP pixel values to detect rain-free pixels of the ERA5 dataset. We calculated the non-conservative air mass tendency over the rain-free pixels on a daily and 0.5 degrees spatio-temporal scale. We found the radiative equilibrium amount of high potential temperature air mass and the Newtonian cooling process with a relaxation time scale of 6 to 7 days. We will show the quantitative estimates of the total convective process of the Asian summer monsoon system associated with the convective clouds and radiative processes, through the mass budget of 355K-370K potential temperature air mass. We will further show results of the evaluation of the accuracy of TRMM and GPM products using our high-resolution tipping bucket raingauge network distributed over the Northeastern Indian subcontinent.</p>


2019 ◽  
Vol 19 (9) ◽  
pp. 6535-6549 ◽  
Author(s):  
Bojan Škerlak ◽  
Stephan Pfahl ◽  
Michael Sprenger ◽  
Heini Wernli

Abstract. Upper-level fronts are often associated with the rapid transport of stratospheric air along tilted isentropes to the middle or lower troposphere, where this air leads to significantly enhanced ozone concentrations. These plumes of originally stratospheric air can only occasionally be observed at the surface because (i) stable boundary layers prevent an efficient vertical transport down to the surface, and (ii) even if boundary layer turbulence were strong enough to enable this transport, the originally stratospheric air mass can be diluted by mixing, such that only a weak stratospheric signal can be recorded at the surface. Most documented examples of stratospheric air reaching the surface occurred in mountainous regions. This study investigates two such events, using a passive stratospheric air mass tracer in a mesoscale model to explore the processes that enable the transport down to the surface. The events occurred in early May 2006 in the Rocky Mountains and in mid-June 2006 on the Tibetan Plateau. In both cases, a tropopause fold associated with an upper-level front enabled stratospheric air to enter the troposphere. In our model simulation of the North American case, the strong frontal zone reaches down to 700 hPa and leads to a fairly direct vertical transport of the stratospheric tracer along the tilted isentropes to the surface. In the Tibetan Plateau case, however, no near-surface front exists and a reservoir of high stratospheric tracer concentrations initially forms at 300–400 hPa, without further isentropic descent. However, entrainment at the top of the very deep boundary layer (reaching to 300 hPa over the Tibetan Plateau) and turbulence within the boundary layer allows for downward transport of stratospheric air to the surface. Despite the strongly differing dynamical processes, stratospheric tracer concentrations at the surface reach peak values of 10 %–20 % of the imposed stratospheric value in both cases, corroborating the potential of deep stratosphere-to-troposphere transport events to significantly influence surface ozone concentrations in these regions.


2009 ◽  
Vol 137 (6) ◽  
pp. 1972-1990 ◽  
Author(s):  
Stanley B. Trier ◽  
Robert D. Sharman

Abstract Widespread moderate turbulence was recorded on three specially equipped commercial airline flights over northern Kansas near the northern edge of the extensive cirrus anvil of a nocturnal mesoscale convective system (MCS) on 17 June 2005. A noteworthy aspect of the turbulence was its location several hundred kilometers from the active deep convection (i.e., large reflectivity) regions of the MCS. Herein, the MCS life cycle and the turbulence environment in its upper-level outflow are studied using Rapid Update Cycle (RUC) analyses and cloud-permitting simulations with the Weather Research and Forecast Model (WRF). It is demonstrated that strong vertical shear beneath the MCS outflow jet is critical to providing an environment that could support dynamic (e.g., shearing type) instabilities conducive to turbulence. Comparison of a control simulation to one in which the temperature tendency due to latent heating was eliminated indicates that strong vertical shear and corresponding reductions in the local Richardson number (Ri) to ∼0.25 at the northern edge of the anvil were almost entirely a consequence of the MCS-induced westerly outflow jet. The large vertical shear is found to decrease Ri both directly, and by contributing to reductions in static stability near the northern anvil edge through differential advection of (equivalent) potential temperature gradients, which are in turn influenced by adiabatic cooling associated with the mesoscale updraft located upstream within the anvil. On the south side of the MCS, the vertical shear associated with easterly outflow was significantly offset by environmental westerly shear, which resulted in larger Ri and less widespread model turbulent kinetic energy (TKE) than at the northern anvil edge.


2013 ◽  
Vol 26 (13) ◽  
pp. 4535-4549 ◽  
Author(s):  
Giacomo Masato ◽  
Brian J. Hoskins ◽  
Tim Woollings

Abstract This paper generalizes and applies recently developed blocking diagnostics in a two-dimensional (2D) latitude–longitude context, which takes into consideration both mid- and high-latitude blocking. These diagnostics identify characteristics of the associated wave breaking as seen in the potential temperature θ on the dynamical tropopause, particularly the cyclonic or anticyclonic direction of wave breaking (“DB index”) and the relative intensity (“RI index”) of the air masses that contribute to blocking formation. The methodology is extended to a 2D domain and a cluster technique is deployed to classify mid- and high-latitude blocking according to the wave-breaking characteristics. Midlatitude blocking is observed over Europe and Asia, where the meridional gradient of θ is generally weak, whereas high-latitude blocking is mainly present over the oceans, to the north of the jet stream, where the meridional gradient of θ is much stronger. They occur on the equatorward and poleward flank of the jet stream, respectively, where the horizontal shear ∂u/∂y is positive in the first case and negative in the second case. A regional analysis is also conducted. Warm-cyclonic blocking over the Pacific and cold-anticyclonic blocking over Europe are identified as the most persistent types and are associated with large amplitude anomalies in temperature and precipitation. Finally, the high-latitude cyclonic events seem to correlate well with low-frequency modes of variability over the Pacific and Atlantic Oceans.


Sign in / Sign up

Export Citation Format

Share Document