Next Article in Journal
Influence of Local Climate and ENSO on the Growth of Cedrela odorata L. in Suriname
Next Article in Special Issue
Oscillatory Signatures in the Raindrop Motion Relative to the Air Medium with Terminal Velocity
Previous Article in Journal / Special Issue
Trends in Summer-Time Tropospheric Ozone during COVID-19 Lockdown in Indian Cities Might Forecast a Higher Future Risk
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Precipitation Extremes and Their Synoptic Models in the Northwest European Sector of the Arctic during the Cold Season

by
Alexander Kislov
1,*,
Tatiana Matveeva
2,3 and
Uliana Antipina
4
1
Department of Meteorology and Climatology, Lomonosov Moscow State University, Leninskie Gory MGU, GSP-1, 119991 Moscow, Russia
2
Institute of Geography Russian Academy of Science, Staromonetny Line, 29, 119017 Moscow, Russia
3
A.M. Obukhov Institute of Atmospheric Physics Russian Academy of Sciences, Pyzhyovskiy Lane, 3, 119017 Moscow, Russia
4
Israel Institute of Global Climate and Ecology, Glebovskaya 20B, 107258 Moscow, Russia
*
Author to whom correspondence should be addressed.
Atmosphere 2022, 13(7), 1116; https://doi.org/10.3390/atmos13071116
Submission received: 21 June 2022 / Revised: 11 July 2022 / Accepted: 13 July 2022 / Published: 15 July 2022
(This article belongs to the Special Issue Feature Papers in Atmosphere Science)

Abstract

:
Precipitation extrema over the Barents Sea and the neighbouring locations in Europe were analysed using data obtained from station observations and a highly detailed ERA5 re-analysis dataset. These data did not always spatially coincide (on average, coincidence was ~50%). Daily amounts of precipitation were typically higher in the observation data, although there may be a reverse picture. The analysis revealed that at several stations and in many of the ERA5 grids, the set of precipitation extremes exists as a mixture of two different subsets. The cumulative distribution functions (CDF) of the largest population in the context of both the re-analysis and observational data are well described by Pareto’s law. However, very rare cases exist in which the values deviate and exceed this base distribution value in regions possessing large values. These super-large anomalies do not obey the statistical law common to all other extremes. However, this does not mean that the extremes can be arbitrarily large. They do not exceed the marginal values that are typical for this type of climate and season. The analysis confirms that extreme precipitation in the western sector of the Arctic is caused by the penetration of moist air masses from the Atlantic in the circulation systems of intense cyclones. At certain times, mesoscale convective systems are embedded in atmospheric fronts and can significantly contribute to the formation of precipitation. Intensification of such cyclones corresponding to global warming should lead to a transformation of typical CDF, as modern outliers will become regular components of the Pareto law. This change in the statistics of extreme events reflects the nonstationarity of the climate state. The influence of polar lows on the formation of large daily precipitation amounts is not felt.

1. Introduction

Recently, interest has increased in regard to the study of the meteorological conditions of the Arctic due to the economic development of this region. This region is an area of intensive shipping, fishery, and natural gas production. From a climate perspective, the region is of great interest due to the observation that this area has heated more quickly than has any other region on Earth. This is the so-called Arctic amplification [1] that has been occurring and that results from dynamic processes in the upper ocean layers due to increasing upward heat flux (so-called “atlantisation process” [2]), and in the atmosphere due to intensification of poleward heat advection [3,4]. These effects are particularly significant in the Barents Sea. The geographical position of this region causes the Barents Sea to exist as a transitional zone where the climate is primarily influenced by the North Atlantic Ocean. Notably, relatively little research has been published regarding the climate of this region.
The purpose of this article is to study extreme precipitation that occurs on the marine surface and the neighbouring portion of land consisting of the Scandinavian and Cola peninsulas during the cold season. We focus on a small portion of the Pan-Arctic domain to more clearly delineate its regional characteristics. We have taken into account that in the sub-Arctic region, the winter season lasts throughout the entire interval from November to March. At this time of year, particularly large amounts of precipitation are associated with cyclones that arise from the west and southwest and sometimes originate over this region. The extreme values of precipitation during the warm season were studied in our recent work [5].
A large portion of the Barents Sea is covered with ice during the cold season. From the perspectives of temperature, roughness, and other characteristics, the sea surface merges with land. The last two decades have been characterized by significant climate heating, a reduction in the surface of the old and first-year sea ice in the Arctic, and the appearance of a significantly larger ice-free marine surface compared to that which existed earlier.
The Arctic region is characterised by a sparse in situ observational network. Re-analysis data provide an important alternative for filling this gap, as they provide global coverage and combine weather forecast models and the assimilation of observations from a wide variety of sources. For this purpose, a horizontally detailed re-analysis of ERA5 was performed. This product (see below) was developed by the European Centre for Medium-Range Weather Forecasts (ECMWF).
Extreme precipitation events can be characterised using various indicators. Here, we used daily amounts. For statistical description of their extrema, the Weibull, Frechet, and Gumbel distributions (and their relatives, such as beta and gamma distributions and their mixes) and also energetic distributions (such as the Pareto distribution) have traditionally been applied (see [6,7,8,9,10,11]). We continue to work in this tradition. Extreme value theory assumes that the data selected for analysis must be independent and identically distributed. However, extreme value theory can be extended for dependent time series [12,13] or for data that have previously been transformed into independent series [14,15]. Regarding the identical distribution, our analysis utilised different extreme values (wind speed [15], precipitation [5], and wave climate [16] and revealed that the set of extremes obtained from observations or from modelling is a mixture of two different subsets. The probability distribution function of the largest population volume aptly describes one of the above-mentioned lows. This approximation was used as the base distribution. The largest events belonging to this base distribution are sometimes referred to as “black swans” based on the terminology introduced by N. Taleb [17]. However, there are very rare cases in which the values deviate and exceed the base distribution in the large-value area. These unique events are termed “dragons” based on the terminology introduced by D. Sornette [18].
Statistical analysis that could be used to differentiate extreme values is substantial not only due to the need for practical goals but also because they permit, in some cases, the detection of the origins of extrema, as various statistical distributions sometimes suggest different originating mechanisms. This hypothesis motivated our interest in detecting, analysing, and understanding such different extrema and their nature. It is important to understand an extension of the ability of atmospheric re-analysis to reproduce extrema. Hence, we plan to investigate utility of the ERA5 data (meteorological information with a fine spatial resolution) in regard to simulating the peculiarities of precipitation extremes.
The next step of our analysis is to detect the synoptic processes responsible for such events. We understand that precipitation extrema in high latitudes connect with extratropical cyclones (see, for example, [19]); however, this mechanism has not been fully explored. Therefore, we focus on case studies that aid us in gaining a detailed understanding of the processes involved in extreme precipitation occurring along frontal structures, in extreme cyclones, and in polar lows [20].
In the following section, we describe the study area and briefly summarise the information regarding the datasets. The next sections provide evidence of a Pareto distribution in the station observation data. Subsequently, we describe the re-analysis data and compare the fit to the observations. In the next section, we describe the synoptic conditions that lead to precipitation extremes. The Section 5 concludes the paper.

2. Data

The study was performed in the sub-Arctic realm of the Barents Sea and included both the Atlantic Northeast and inland European territory. A dataset of daily precipitation amounts (covering the period 1966–2018) was used from stations located within the domain (Figure 1).
The ERA5 re-analysis developed by the ECMWF was used as a source of re-analysis data. Its horizontal resolution was increased to 0.25 × 0.25. Compared to the values in a previous successful ERA-Interim re-analysis, the number of vertical levels was increased to 137 pressure levels from 1 to 1000 hPa, the temporal resolution was changed to hourly, and the list of output parameters was extended. Furthermore, the number of assimilated observations was enhanced (approximately five-fold greater than that of ERA-Interim) [21,22].
For the analysis of synoptic situations, we used 6 h synoptic maps that were archived by the Deutscher Wetterdienst (https://www1.wetter3.de/archiv_dwd_dt.html, accessed on 1 February 2022). Information regarding polar lows that was required to assess their roles in precipitation extrema was provided by the STARS data set image database (https://projects.met.no/polarlow/stars-dat/, accessed on 1 February 2022) and PANGAEA (https://doi.org/10.1594/PANGAEA.903058, accessed on 1 February 2022). These data were analysed in previous studies [23,24,25,26].

3. Methodology

To apply statistical approaches, we composed our data based on the independence condition. Practically, this indicates that the data sample was required to include only independent extreme values. We demonstrated (both for the station and for the ERA5 grid cells) that independence is achieved if the time series consists of the daily precipitation amount for each day. This dead-time interval (1 day) was obtained via autocorrelation function analysis as a period for the disappearance of the correlation between consequent fluctuations.
During the winter, the Barents Sea can be either open water or covered with ice at various concentrations. In our analysis, we only used open water samples from the western portion of the sea. Over the ice-free sea, precipitation can occur in the form of either snow or rain. During the cold season, precipitation over the cold land surface falls mainly in the form of snow.
Our analysis revealed that the daily precipitation time series cannot be satisfactorily described by a low-parameter formula from the known probability distribution laws. However, an approximation can be achieved by utilizing a mixture of distributions. For example, a mixture of several (at least five) beta distributions reproduced all details of the cumulative distribution function (CDF). However, this CDF is described by 30 parameters that must be calculated from the samples. It is difficult to recognise this solution as satisfactory, as such a parameterisation is unstable and spatially non-universal (because at each grid point or station, the CDF exhibits a unique number of parameters).
Therefore, an alternative approach was utilised in which the total sample size was not analysed, but only its largest values that exceeded a certain threshold. The statistics for extreme precipitation are described by the Pareto distribution for realising the peaks-over-threshold method. The CDF for this approach is:
W(p) = 1 − ((ΣPth)/(ΣP))γ
where p is daily precipitation amount and Pth is threshold value.
We can stress that the advantage of a multiparameter formula (based on a mixture of beta distributions) over a simpler Pareto formula is small. Taking into account this circumstance and also the knowledge that we study the statistical regime primarily to identify the physical patterns of the extreme formation, the Pareto formula was used.
Moreover, our analysis allowed us to obtain a single threshold value (9 mm) for the entire region. This provides an important advantage, as it is now possible to concentrate our attention on the analysis of γ as a single parameter. This threshold value removed ~98% of the sample, and approximately 150–200 values remained for analysis. The Pareto parameters were estimated using the maximum likelihood method. To this end, expression (1) can be replaced by:
ln(1 − W(p)) = γ(ln(ΣPth) − ln(ΣP))
In Figure 2, several Pareto plots are presented as calculated based on the station data; a straight line was recovered if the sample exhibits a Pareto distribution. We can visually observe the quality of the description and also quantify it based on the coefficient of determination (R2). Which provides a measure of approximation success. In a mathematical sense, the use of R2 is related to the application of the Cramer–Mises–Smirnov statistical criterion. At all stations, we observed that the majority of the points of the empirical CDF exhibited close approximations of a Pareto distribution. This is mentioned above the base distribution denoting the community of black swans. The application of the Kolmogorov–Smirnov test also indicated that there was no reason not to trust the Pareto distribution.
However, all peculiarities are again not described. There are several points of the CDF that depict predominantly high values that do not belong to the base distribution. We have not tried to find common formulae but conclude that, as was mentioned previously, representatives of such populations were marked as “dragons” (in comparison to “black swans” representing the largest events belonging to base distribution).
We will continue to follow this terminology by designating the appearance of these anomalies as the “dragon effect” to indicate the differences in samples belonging to various groups.

4. Results and Discussions

4.1. Statistical Features of the Observed Daily Precipitation Sums

The same statistical properties observed in the statistics of extremes could be considered to be results of the same organisation principle, and this suggests a common originating mechanism for each representative of this population. This idea allows us to understand that a large extreme is not distinguished from its small siblings, with the exception of its large power.
The outliers from the base distribution can be detected based on obvious breaks in the tails of the precipitation distributions (Figure 2a). Here, there is a group of several points arranged along a line on the coordinate area of the Pareto distribution. This result falls under the classification discussed in the Section 1, where the sample data of the same item refer to different distribution functions. In other cases, there were only a small number of such extreme values that did not fall into the base distribution (Figure 2c), and finally, there were no such specific events in some stations (Figure 2b,d). A very important observation is that we can easily diagnose events adhering to the base distribution (there are several examples demonstrating that such a diagnosis is not simple and requires different methods adapted to the specific problem [27]).
Note that in statistics, there is the test of a null hypothesis, which states that two samples are derived from the same population, against an alternative hypothesis, which states that a specific population tends to possess larger values than does the other. This is the Mann–Whitney U test (the Kruskal–Wallis test extends the Mann–Whitney U test when there are more than two sets). However, to compare two (or more) sets of observations, they must be selected a priori. In the case of precipitation extremes, such a method is not necessary, as all data refer to the same nomenclature.
Consider the property of the base distribution in more detail. The parameter γ of the Pareto distribution calculated for all the stations is presented in Table 1. Most of the data are in the range of 3 to 5, but sometimes, γ is much higher. They refer to stations located in the interior of the continent (such as Sodankila) or that are blocked by mountains from the influence of the sea (such as Alta or Kvikkjokku). Concerning the Kolguyev Island Northern station, the sample size of anomalies exceeding the threshold was very small, and this may have led to a distortion of the result. It should be noted that in conditions where the threshold value is approximately the same, a significantly larger γ indicates that large anomalies occur less frequently.

4.2. The Pareto Distributions in Data of the ERA5

The next step of the analysis was to investigate the extent to which the above-mentioned peculiarities of precipitation extremes are reproduced by ERA5 re-analysis. The establishment of the correspondence between simulation products and near-surface observations could help us to assess the quality of the modelling products and their ability to reproduce the precipitation extremes. Additionally, it was important to extend our analysis to the marine area, as this cannot be implemented based on station data only.
In Figure 3, we plotted several CDFs on the bases of the ERA5 data. Again, as observed in the station data, practically all points of the CDFs from the EPA5 exhibit close approximations of a Pareto distribution. Only a small number of points do not consist of this base distribution. In some cases, these points exit the general dependence and were arranged along a line. They belong again in the Pareto distribution; however, there are much larger values for γ. In other cases, there are several or only one such extreme value, and finally, there are no such specific events. These results fall under the situation discussed above, where the sample data for the same item refer to different distribution functions.
It is important to complement these similar properties of the CDFs by comparing information regarding rare events (black swans and dragons) obtained from observations and by using the ERA5 for the same day at the same grid point. For such a comparison, 10 dates were selected (based on the CDFs) in which the largest anomalies were observed. Next, we considered the extrema at the grid points close to the station and studied if the ERA5 anomalies exhibit similar sizes and origins on the same dates.
First, it was observed that the dates did not always coincide. Coincidence occurred in 10–80% of cases, and on average, coincidence was ≈50% (Figure 4). The absolute values of the daily amounts of precipitation were higher in the observational data, although there may be a reverse picture (Figure 4).
For the selected extrema, we analysed if they belonged to the community of black swans (B) or dragons (D). Species coincidence (D-D or B-B) was noted in approximately 60% of cases (Figure 5). In 35% of cases, event D (according to observations) was denoted as B according to the ERA5 data. The opposite was true only in 5% of the cases. This result reflects the observation that the absolute values of the extrema according to observations typically exceed those yielded by reanalysis.
Thus, the extremes of daily precipitation are adequate (from the point of view of belonging to a particular family) and are reflected in the re-analysis data in only approximately 25–30% of cases. In another 15–20% of cases, D in the observational data is converted to B in ERA5.
It should be noted that despite the observation that the extrema in the series of observational data and in the re-analysis data do not always appear synchronously, the exponents γ (see expression (1)) extracted from both the observation and ERA5 data are approximately the same. All ERA5 data over the land territory of the selected region lie between the mean values of the observed γ and its standard deviations. Moreover, ERA5 data were approximately the same over land and sea areas.
Another important effect detected among the samples belonging to the base distribution is the gradual increase in the parameter γ from west to east (Figure 6), which indicates that the frequency of large anomalies decreases in this direction. As mentioned, in conditions where the threshold value is approximately the same, a significantly larger γ indicates that large anomalies occur less frequently. Considering that large anomalies in precipitation are maximally frequent with a cyclonic type of circulation (e.g., [19]; see the next section), we can conclude that this effect reflects a decrease in cyclonic activity in the west–east direction, and cyclones are typically destroyed in the Yamal–Taimyr region.
Therefore, the ERA5 generates extrema that exhibit statistical properties that are similar to observations; however, in both cases, these are random processes, and their extrema do not always coincide.

4.3. Precipitation Extrema and Extratropical Cyclones

As mentioned previously, one of the important features of the CDFs of precipitation extrema is the observation that most powerful cases did not fall into the base distribution (Figure 2 and Figure 3). Such “super-large” anomalies are unique, and therefore, they are presented in the sample as extremely rare events. Additionally, as mentioned, their return periods exceed the return periods of the same anomalies prescribed by base distribution. Typically, it is difficult to generalise the features of the dragon effect in space from both station data and ERA5 data, as the parameters of “super-large” anomalies vary at neighbouring stations or at neighbouring grid points. We attempted to use synoptic analysis to determine the circulation conditions that were responsible for the “superlarge” precipitation anomalies, and these efforts are coordinated with the dragon events as determined from the CDFs (Figure 2 and Figure 3). Priority was given to events that were simultaneously recorded as anomalous at several stations (or grid nodes).
Thus, in many locations in the Barents Sea region, anomalies of the highest ranks were noted on 10 January 2019 (see Table 2). At that time, the region was under the influence of a cyclone in the mature stage, and the centre of this cyclone was located in the Svalbard region (Figure 7).
Air from the temperate latitudes of the Atlantic Ocean was transferred by cyclonic circulation to the Barents Sea region. As a result of occluding, the warm and moist air spread into the middle and upper troposphere covered not only the warm sector of the cyclone but also filled the central and rear areas of the vortex. The analysis of the thermodynamic diagrams (Atmospheric Soundings (uwyo.edu)) proves that the layered rain/snow clouds (a complex of Ns-As-Cs) occupied the entire troposphere. Therefore, significant precipitation fell over a large area, and the anomaly was designated as the greatest at many points within the entire region. The largest daily amount was almost 30 mm. However, for the Arctic (particularly in winter) this is a significant anomaly.
Consider the following episode (Figure 8): The situation appears similar to that discussed earlier (Figure 7). However, atmospheric sounding data show that in the lower and middle troposphere (between 800 and 500 hPa), the temperature is close to the moist adiabatic lapse rate, thereby manifesting the development of convective clouds. The shear of wind in the lower troposphere confirms the ascending air movements generated by the frontal system. Moist convection conditions and large wind shear, together with a large moisture content, contributed to the development of a pseudo-adiabatic process to initiate much precipitation.

4.4. Precipitation Extrema and Role of the Polar Lows

It is interesting to compare extreme precipitation to polar lows. For example, an interesting feature of the situation on 6–7 January was the presence of a polar low that was located near a pixel (72 N, 28 E) [26]. It is clear that it could not be responsible for the observed wide zone of precipitation anomalies that were particularly distinguished along 73 N (see Table 2).
Consider the next situation on 27 November 2019 when the entire region was again under the influence of a polar cyclone. Along its western periphery, the polar low moved outside the baroclinic zone to the south. The following maximum of precipitation that appeared in a few days (29–30 November 2019) over a large area (Table 2) was not related to this polar low. The extreme precipitation was generated by a synoptic-scale vortex.
In general, our analysis of several hundred cases of polar lows that are collected in different archives [23,24,25,26] revealed that their role in the formation of precipitation extremes of daily sums is insignificant. Only in rare episodes (for example, on 19 January 2000) can we note a weak local extreme of precipitation occupying one or two pixels near the location of the polar low (Table 2); however, the anomaly was far from the highest ranks. Small polar lows, fast movement speeds, and short lifetimes do not contribute to the formation of extended long-term precipitation anomalies. Additionally, it is established that at least some portion of polar lows are formed within the outbreaks of cold air from ice fields to the warm sea surface. The moisture reserve in cold air is small, and based on this, there is no reason to expect a large amount of precipitation.

5. Conclusions

  • It was demonstrated that the Pareto distributions are close for the extremes that were extracted from both the re-analysis data and observational data. However, the strongest episodes did not fall into this distribution. For them, the correspondence between the absolute value of the anomaly and its probability is completely lost, thereby signifying that any anomaly can occur. This means that super-large anomalies do not obey the statistical law common to all other extremes. It is no wonder, according to this sign, they are labelled as dragons. However, the statement “any anomalies” does not mean that the extremes can be arbitrarily large. They do not exceed the marginal values that are typical for this type of climate and season. This effect appears to be one of saturation when the limit is approached.
  • Let us focus on the widely used practice of utilizing quantile values (k = 0.95 or 0.99) for the analysis of precipitation anomalies. Clearly, this procedure is correct if the CDF is adequate for the entire domain and if the distribution of events exceeding the selected quantile value is given by the expression 1-W(pk). However, in the cases discussed in this article, this leads to distortions, as the structure of the “tip of the tail” of the distribution is completely different from that of the regular base function. The procedure of “pulling” the base Pareto law into the region of “superlarge” precipitation significantly underestimates the return periods of the observed precipitation extrema. This circumstance requires careful use of quantile analysis. For this, it is necessary to consider the situation in which the samples belong to different distribution functions.
  • Extreme precipitation in the Barents region of the Arctic during the cold period was caused by the advection of moist air masses from the Atlantic Ocean. The transfer occurs in the circulation system of the intense polar front cyclones. Sometimes mesoscale convective systems are embedded in atmospheric fronts, thereby making a significant contribution to the formation of precipitation. The influence of polar lows on the formation of large amounts of precipitation was practically non-existent.
  • Intensive cyclones penetrate the Arctic more often during climate warming. In winter, this happens due to the modification of the regional circulation, the increase in meridional processes during blocking [28], and due to the increase in the open water zone (due to the reduction in the area of sea ice), in which it is easier for cyclones to move. Cyclones work as a positive feedback factor. This should lead to a transformation of typical CDF, as modern “chaotic” dragons will become “regular” black swans. This change in the statistics of extreme events reflects the nonstationarity of the climate state.

Author Contributions

Conceptualization, A.K. and T.M.; methodology, A.K.; software, T.M.; validation, T.M. and U.A.; formal analysis, A.K.; investigation, T.M.; resources, U.A.; data curation, U.A.; writing—original draft preparation, U.A.; writing—review and editing, U.A.; visualization, T.M.; supervision, A.K.; project administration, A.K.; funding acquisition, U.A. All authors have read and agreed to the published version of the manuscript.

Funding

The research of A.Kislov was funded by Russian Ministry of Science and Higher Education contract number 075-15-2021-574. The research of T. Matveeva was funded by the Institute of Geography, Russian Academy of Sciences, the State Assignment number FMGE-2019-0009, as well the Russian Science Foundation (RSF) grant number 19-17-00242.

Data Availability Statement

General data used in this study is archived in the repository (Datasets Generated: “Mendeley Data”, https://data.mendeley.com/datasets/vdy2nksk4h/1, accessed on 1 February 2022).

Conflicts of Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

References

  1. Serreze, M.C.; Barry, R.G. Processes and Impacts of Arctic Amplification: A Research Synthesis. Glob. Planet. Chang. 2011, 77, 85–96. [Google Scholar] [CrossRef]
  2. Polyakov, I.V.; Pnyushkov, A.V.; Alkire, M.B.; Ashik, I.M.; Baumann, T.M.; Carmack, E.C.; Goszczko, I.; Guthrie, J.; Ivanov, V.V.; Kanzow, T.; et al. Greater role for Atlantic inflows on sea-ice loss in the Eurasian Basin of the Arctic Ocean. Science 2017, 356, 285–291. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Cohen, J.; Screen, J.A.; Furtado, J.C.; Barlow, M.; Whittleston, D.; Coumou, D.; Francis, J.; Dethloff, K.; Entekhabi, D.; Overland, J.; et al. Recent Arctic Amplification and Extreme Mid-Latitude Weather. Nat. Geosci. 2014, 7, 627–637. [Google Scholar] [CrossRef] [Green Version]
  4. Ye, K.; Wu, R.; Liu, Y. Interdecadal Changes of Eurasian Snow, Surface Temperature and Atmospheric Circulation in the Late 1980s. J. Geophys. Res. Atmos. 2015, 120, 2738–2753. [Google Scholar] [CrossRef]
  5. Kislov, A.; Antipina, U.I.; Korneva, I.A. Precipitation extrema over the European sector of the Arctic during the summer time: Statistics and synoptic models. Russ. Meteorol. Hydrol. 2012, 46, 434–443. [Google Scholar] [CrossRef]
  6. Zolina, O.; Kapala, A.; Simmer, C.; Gulev, S.K. Analysis of extreme precipitation over Europe from different reanalyses: A comparative assessment. Glob. Planet. Chang. 2004, 44, 129–161. [Google Scholar] [CrossRef]
  7. Ben Alaya, M.A.; Zwiers, F.W.; Zhang, X. A bivariate approach to estimating the probability of very extreme precipitation events. Weather Clim. Extrem. 2020, 30, 100290. [Google Scholar] [CrossRef]
  8. Katz, R.W. Extreme value theory for precipitation: Sensitivity analysis for climate change. Adv. Water Resour. 1999, 23, 133–139. [Google Scholar] [CrossRef]
  9. Kharin, V.V.; Zwiers, F.W.; Zhang, X.B.; Hegerl, G.C. Changes in temperature and precipitation extremes in the IPCC ensemble of global coupled model simulations. J. Clim. 2007, 20, 1419–1444. [Google Scholar] [CrossRef] [Green Version]
  10. Semenov, V.A.; Bengtsson, L. Secular trends in daily precipitation characteristics: Greenhouse gas simulation with a coupled AOGCM. Clim. Dyn. 2002, 19, 123–140. [Google Scholar] [CrossRef] [Green Version]
  11. Maraun, D.; Wetterhall, F.; Ireson, A.M.; Chandler, R.E. Precipitation downscaling under climate change. Recent developments to bridge the gap between dynamical models and the end user. Rev. Geophys. 2010, 48, RG3003. [Google Scholar] [CrossRef]
  12. Beirlant, J.; Goegebeur, Y.; Segers, J.; Teugels, J.; De Waal, D.; Ferro, C. Statistics of Extremes: Theory and Applications; Willey Series in Probability and Statistics; John Wiley & Sons Ltd.: Chichester, UK, 2004. [Google Scholar] [CrossRef]
  13. Coles, S. An Introduction to Statistical Modeling of Extreme Values; Springer Series in Statistics; Springer: Berlin/Heidelberg, Germany, 2001. [Google Scholar]
  14. Coles, S.; Walshaw, D. Directional Modelling of Extreme Wind Speeds. J. R. Stat. Soc. Ser. C (Appl. Stat.) 1994, 43, 139–157. [Google Scholar] [CrossRef]
  15. Kislov, A.; Matveeva, T. An Extreme Value Analysis of Wind Speed over the European and Siberian Parts of Arctic Region. Atmos. Clim. Sci. 2016, 6, 205–223. [Google Scholar] [CrossRef] [Green Version]
  16. Myslenkov, S.; Platonov, V.; Kislov, A.; Silvestrova, K.; Medvedev, I. Wave climate and storm activity in the Kara sea. Thirty-Nine-Year Wave Hindcast, Storm Activity, and Probability Analysis of Storm Waves in the Kara Sea, Russia. Water 2021, 13, 648. [Google Scholar] [CrossRef]
  17. Taleb, N.N. The Black Swan: The Impact of the Highly Improbable, 2nd ed.; Penguin: New York, NY, USA, 2010. [Google Scholar]
  18. Sornette, D. Dragon-Kings, Black Swans and the Prediction of Crises. Int. J. Terraspace Sci. Eng. 2009, 2, 1–18. [Google Scholar] [CrossRef]
  19. Kislov, A.; Matveeva, T. The Monsoon over the Barents Sea and Kara Sea. Atmos. Clim. Sci. 2020, 10, 339–356. [Google Scholar] [CrossRef]
  20. Rasmussen, E.A.; Turner, J. Polar Lows: Mesoscale Weather Systems in the Polar Regions; Cambridge University Press: Cambridge, UK, 2003. [Google Scholar]
  21. Hersbach, H.; Dee, D. ERA-5 Reanalysis Is in Production. ECMWF Newsl. 2016, 147, 5–6. [Google Scholar]
  22. Hersbach, H.; Bell, B.; Berrisford, P.; Horanyi, A.; Sabater, J.M.; Nicolas, J.; Radu, R.; Schepers, D.; Simmons, A.; Soci, C.; et al. Global Reanalysis: Goodbye ERAInterim, Hello ERA5. ECMWF Newsl. 2019, 159, 17–24. [Google Scholar]
  23. Rojo, M.; Claud, C.; Noer, G.; Carleton, A. In situ measurements of surface winds, waves, and sea state in polar lows over the North Atlantic. J. Geophys. Res. Atmos. 2019, 124, 700–718. [Google Scholar] [CrossRef] [Green Version]
  24. Smirnova, J.E.; Golubkin, P.A.; Bobylev, L.P.; Zabolotskikh, E.V.; Chapron, B. Polar low climatology over the Nordic and Barents seas based on satellite passive microwave data. Geophys. Res. Lett. 2015, 42, 5603–5609. [Google Scholar] [CrossRef] [Green Version]
  25. Revokatova, A.; Nikitin, M.; Rivin, G.; Rozinkina, I.; Nikitin, A.; Tatarinovich, E. High-Resolution Simulation of Polar Lows over Norwegian and Barents Seas Using the COSMO-CLM and ICON Models for the 2019–2020 Cold Season. Atmosphere 2021, 12, 137. [Google Scholar] [CrossRef]
  26. Golubkin, P.; Smirnova, J.; Bobylev, L. Satellite-Derived Spatio-Temporal Distribution and Parameters of North Atlantic Polar Lows for 2015–2017. Atmosphere 2021, 12, 224. [Google Scholar] [CrossRef]
  27. Sornette, D.; Ouillon, G. Dragon-Kings: Mechanisms, Statistical Methods and Empirical Evidence. Eur. Phys. J. Spec. Top. 2012, 205, 1–26. [Google Scholar] [CrossRef] [Green Version]
  28. Kislov, A.; Sokolikhina, N.; Semenov, E.; Tudriy, K. Blocking Anticyclone in the Atlantic Sector of the Arctic as an Example of an Individual Atmospheric Vortex. Atmos. Clim. Sci. 2017, 7, 323–336. [Google Scholar] [CrossRef] [Green Version]
Figure 1. The area of interest and locations of observation stations.
Figure 1. The area of interest and locations of observation stations.
Atmosphere 13 01116 g001
Figure 2. Empirical CDFs of daily precipitation amounts (station observations) straightened on the coordinate axis of the Pareto distribution, and a linear regression line corresponding to the Pareto function. (a) Fruholmen Fyr, (b) Kandalaksha, (c) Teriberka, (d) Arkhangelsk. In all cases, R2 > 0.96 for basic distribution (blue dots). The yellow dots are deviations from the Pareto function (the so-called “dragons”).
Figure 2. Empirical CDFs of daily precipitation amounts (station observations) straightened on the coordinate axis of the Pareto distribution, and a linear regression line corresponding to the Pareto function. (a) Fruholmen Fyr, (b) Kandalaksha, (c) Teriberka, (d) Arkhangelsk. In all cases, R2 > 0.96 for basic distribution (blue dots). The yellow dots are deviations from the Pareto function (the so-called “dragons”).
Atmosphere 13 01116 g002
Figure 3. The CDFs of daily precipitation amounts (ERA5 data) straightened on the coordinate axis of the Pareto distribution, and linear regression line corresponding to the Pareto function: (a) 75 N, 24 E (the lower points on the chart correspond 10 January 2019 and 30 November 2019); (b) 74 N, 33 E (the lower point on the chart corresponds 30 November 2019); (c) 74 N, 41 E (the lower point on the chart corresponds to 30 November 2019, 5th bottom point on the chart corresponds to 19 January 2000; (d) 73 N, 31 E (the lower point on the chart corresponds 6 January 2009). In all cases R2 > 0.98 for basic distribution (blue dots). The yellow dots are deviations from the Pareto function (the so-called “dragons”). CDFs discussed below (see text).
Figure 3. The CDFs of daily precipitation amounts (ERA5 data) straightened on the coordinate axis of the Pareto distribution, and linear regression line corresponding to the Pareto function: (a) 75 N, 24 E (the lower points on the chart correspond 10 January 2019 and 30 November 2019); (b) 74 N, 33 E (the lower point on the chart corresponds 30 November 2019); (c) 74 N, 41 E (the lower point on the chart corresponds to 30 November 2019, 5th bottom point on the chart corresponds to 19 January 2000; (d) 73 N, 31 E (the lower point on the chart corresponds 6 January 2009). In all cases R2 > 0.98 for basic distribution (blue dots). The yellow dots are deviations from the Pareto function (the so-called “dragons”). CDFs discussed below (see text).
Atmosphere 13 01116 g003
Figure 4. The proportion of coincidences in the 10 observed (at all stations—see the Figure 1) largest precipitation extremes with similar extrema according to ERA5 data at these grid points (abscissa axis) and the ratio of observed precipitation (Pobs) to the total of ERA5 data (Pera5) (ordinate axis).
Figure 4. The proportion of coincidences in the 10 observed (at all stations—see the Figure 1) largest precipitation extremes with similar extrema according to ERA5 data at these grid points (abscissa axis) and the ratio of observed precipitation (Pobs) to the total of ERA5 data (Pera5) (ordinate axis).
Atmosphere 13 01116 g004
Figure 5. Combinations of coincidences or inconsistencies in the nomenclature of events D and B according to observations and according to re-analysis data.
Figure 5. Combinations of coincidences or inconsistencies in the nomenclature of events D and B according to observations and according to re-analysis data.
Atmosphere 13 01116 g005
Figure 6. Distribution of the parameter γ (describing the samples belonging to the base distribution) in the west–east direction in the framework 67.5–75 N based on ERA5 data.
Figure 6. Distribution of the parameter γ (describing the samples belonging to the base distribution) in the west–east direction in the framework 67.5–75 N based on ERA5 data.
Atmosphere 13 01116 g006
Figure 7. Weather map valid 18 UTC 10 January 2019.
Figure 7. Weather map valid 18 UTC 10 January 2019.
Atmosphere 13 01116 g007
Figure 8. Weather map valid 00 UTC 7 January 2009.
Figure 8. Weather map valid 00 UTC 7 January 2009.
Atmosphere 13 01116 g008
Table 1. The Pareto distribution parameter γ calculated for each station.
Table 1. The Pareto distribution parameter γ calculated for each station.
Name of Observation Stationsγ
Bardufoss4.0
Tromso/langnes3.0
Tromsoe3.2
Alta6.8
Fruholmen Fyr3.8
Kirkenes lufthavn5.2
Sodankylae7.0
Vayda-Guba3.8
Polyarnoe3.9
Umba3.7
Cape Kanin Nos4.2
Cape Svyatoy Nos4.0
Kandalaksha4.7
Kanevka4.8
Cape Konstantinovsky2.9
Indiga4.6
Kovdor3.5
Kolguyev Island Northern5.8
Krasnoshchelye3.8
Mezen4.2
Arkhangelsk5.0
Teriberka3.9
Vardo5.4
Murmansk4.6
Shoina3.8
Sosnovets island3.8
Janiskoski4.0
Slettnes Fyr3.8
Mean4.4
Mean ± Std3.4–5.4
Table 2. Characteristics of precipitation extrema, discussed in the text.
Table 2. Characteristics of precipitation extrema, discussed in the text.
DataMaximum of Precipitation, mmRegion of AppearancesPrecipitable Water *, kg·m−2Rank
10 January 20192875 N 16–28 E51–3
1874 N 16–29 E<10
1873 N 19–24 E<10
2072 N 16–27 E<10
1871 N 20–21 E<20
1570 N 16–17 E<20
6 January 2009 – 7 January 20093073 N 26–35 E430
30 November 20192775 N 17–37 E71–3
2974 N 26–42 E1–3
1973 N 37–41 E<10
1572 N 16–27 E<20
19 January 20001674 N 40–42 E3<10
1473 N 37–39 E<10
* It was calculated as averaged value based on atmospheric sounding data at stations surrounding the Barents Sea (Murmansk 69.0 N, 33,1 E; Shoina 67.5 N, 44.2 E; Malye Karmakuly 72.4 N, 52.7 E; Polargmo Im. Krenkelja 80.6 N, 58.0 E).
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kislov, A.; Matveeva, T.; Antipina, U. Precipitation Extremes and Their Synoptic Models in the Northwest European Sector of the Arctic during the Cold Season. Atmosphere 2022, 13, 1116. https://doi.org/10.3390/atmos13071116

AMA Style

Kislov A, Matveeva T, Antipina U. Precipitation Extremes and Their Synoptic Models in the Northwest European Sector of the Arctic during the Cold Season. Atmosphere. 2022; 13(7):1116. https://doi.org/10.3390/atmos13071116

Chicago/Turabian Style

Kislov, Alexander, Tatiana Matveeva, and Uliana Antipina. 2022. "Precipitation Extremes and Their Synoptic Models in the Northwest European Sector of the Arctic during the Cold Season" Atmosphere 13, no. 7: 1116. https://doi.org/10.3390/atmos13071116

APA Style

Kislov, A., Matveeva, T., & Antipina, U. (2022). Precipitation Extremes and Their Synoptic Models in the Northwest European Sector of the Arctic during the Cold Season. Atmosphere, 13(7), 1116. https://doi.org/10.3390/atmos13071116

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop