<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Archiving and Interchange DTD v1.0 20120330//EN" "JATS-archivearticle1.dtd">
<article xmlns:xlink="http://www.w3.org/1999/xlink">
  <front>
    <journal-meta />
    <article-meta>
      <title-group>
        <article-title>Unsupervised Severe Weather Detection Via Joint Representation Learning Over Textual and Weather Data</article-title>
      </title-group>
      <contrib-group>
        <contrib contrib-type="author">
          <string-name>Athanasios Davvetas</string-name>
          <email>tdavvetas@iit.demokritos.gr</email>
          <xref ref-type="aff" rid="aff0">0</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>Iraklis A. Klampanos</string-name>
          <email>iaklampanos@iit.demokritos.gr</email>
          <xref ref-type="aff" rid="aff0">0</xref>
        </contrib>
        <aff id="aff0">
          <label>0</label>
          <institution>National Centre for Scientific Research “Demokritos”, Institute of Informatics and Telecommunications</institution>
          ,
          <addr-line>Athens</addr-line>
          ,
          <country country="GR">Greece</country>
        </aff>
      </contrib-group>
      <abstract>
        <p>When observing a phenomenon, severe cases or anomalies are often characterised by deviation from the expected data distribution. However, non-deviating data samples may also implicitly lead to severe outcomes. In the case of unsupervised severe weather detection, these data samples can lead to mispredictions, since the predictors of severe weather are often not directly observed as features. We posit that incorporating external or auxiliary information, such as the outcome of an external task or an observation, can improve the decision boundaries of an unsupervised detection algorithm. In this paper, we increase the efectiveness of a clustering method to detect cases of severe weather by learning augmented and linearly separable latent representations. We evaluate our solution against three individual cases of severe weather, namely windstorms, floods and tornado outbreaks.</p>
      </abstract>
    </article-meta>
  </front>
  <body>
    <sec id="sec-1">
      <title>CCS CONCEPTS</title>
      <p>• Computing methodologies → Artificial intelligence ; Machine
learning; • Applied computing → Physical sciences and
engineering.
Severe weather detection, representation learning, deep learning</p>
    </sec>
    <sec id="sec-2">
      <title>INTRODUCTION</title>
      <p>Anomalies occur in the majority of datasets. They are fairly rare
and are often challenging to detect in an unsupervised setting. Due
to their lower frequency, the majority of normal samples introduces
implicit bias that results in biased predictions. From an
unsupervised perspective, one can assume that these rare occurrences can
be observed in the outliers of the data distribution. Yet, depending
on the application, searching for samples that deviate from the
expected data distribution may not improve the detection of an
unsupervised method.</p>
      <p>In some applications, the occurrence of anomalies might be
expected or it may not be trivial to detect deviation from the observed
data distribution. An example of such application is detecting cases
of severe weather. A heavy rain or windstorm may be considered
as normal, depending on the geographic region or the season, etc.
These otherwise normal circumstances may lead to natural
disasters, even in costly damages or fatalities. However, they can not
always be predicted by observing a physical quantity.</p>
      <p>
        To predict these types of occurrences, we need to incorporate
external or auxiliary information that can efectively augment the
observable features. In this paper, we investigate the efects of
incorporating external information in the form of an auxiliary task
outcome. We achieve this by utilising a deep learning method called
“Evidence Transfer” that incrementally manipulates the latent
representations of an autoencoder according to external categorical
evidence [
        <xref ref-type="bibr" rid="ref3">3</xref>
        ]. Evidence transfer allows for joint representation
learning, based on external categorical evidence retrieved from textual
sources and weather re-analysis data. Evidence transfer
successfully manipulates the initial learned representations, resulting in
increased efectiveness during individual severe weather detection.
2
2.1
      </p>
    </sec>
    <sec id="sec-3">
      <title>DATA AND METHODS</title>
    </sec>
    <sec id="sec-4">
      <title>Weather Re-analysis Data</title>
      <p>
        ERA-Interim [
        <xref ref-type="bibr" rid="ref4">4</xref>
        ] re-analysis data are produced with a sequential
data assimilation scheme during which prior information from
a forecast model is combined with the available observations in
order to estimate the state of the global atmosphere, allowing for
a better description of past atmospheric conditions. Weather
reanalysis data are gridded data (as shown in Figure 1) depicting
atmospheric variables in various timestamps and pressure levels
(gravity-adjusted height), leading to 4D variables. They cover a
time period of up to 40 years, with less than 1° spatial fineness and
6-hour temporal resolution for global region.
      </p>
      <p>
        In our experiments, we used ERA-interim data covering the time
period from January 1 1979 to May 31 2018 with 6-hour temporal
resolution (retrieved from the Research Data Archive of National
Center for Atmospheric Research in Boulder, Colorado1). The
spatial resolution is ≈ 0.7° × 0.7°, containing atmospheric variables
across 37 vertical pressure levels ranging from 1ℎ to 1000ℎ.
We reduce the region of gridded data from global region to a
Cartesian domain that covers Europe. In order to reduce the domain
of our data we used the pre-processor of Weather Research and
Forecast (WRF) Model [
        <xref ref-type="bibr" rid="ref8">8</xref>
        ], named WPS. The new spatial resolution
of our data is of 64 × 64 cells of 75 × 75 in the west-east and
south-north axes.
      </p>
      <p>
        In our study, the atmospheric variable of interest is the
geopotential height (GHT) which can be seen as a gravity-adjusted height.
GHT is often used for its predictive properties [
        <xref ref-type="bibr" rid="ref6 ref7 ref9">6, 7, 9</xref>
        ], as well as,
to extract weather patterns for other downstream tasks [
        <xref ref-type="bibr" rid="ref5">5</xref>
        ]. Severe
weather can be predicted via sequences of patterns in the
geopotential height (e.g. a cyclone can be observed as a circular pattern).
To highlight useful high-level features, such as circular shapes and
edges, we extract embeddings through a pre-trained VGG-16
network on ImageNet. We feed the VGG-16 network with 3 diferent
levels of GHT = 500, 700 and 900 ℎ in similar fashion to using
the RGB channels of an image. Therefore, a single data sample of
shape 3 × 64 × 64 is transformed into an embedding of 64 × 64,
resulting in a total of 4096 features.
2.2
      </p>
    </sec>
    <sec id="sec-5">
      <title>Textual Evidence</title>
      <p>We augment the weather-based embeddings by making use of
textual evidence for historic severe weather events, found in Wikipedia.
For example, to find severe heavy rain occurrences we search
for recorded floods. We extract categorical evidence from textual
sources of Wikipedia pages which associate a date to a severe
weather event.</p>
      <p>For our experiments we extract the following cases of extreme
events in Europe: (1) costly or deadly hailstorms, (2) floods, (3)
tornadoes and tornado outbreaks and (4) severe windstorms.</p>
      <p>Each of these event types is treated as a binary classification task
for predicting a specific severe weather case. The occurrence date
is used to both reference the weather re-analysis data, as well as
the individual tasks. Since the events listed in Wikipedia do not
typically supply exact times, we label the whole day of reference
as severe, therefore, the minimum span of an event is one day or
four 6-hour samples (we remind that the weather re-analysis data
are provided in 6-hour increments).</p>
      <p>For each of the aforementioned lists we extract the following
ifelds (for simplicity purposes, “Event” is used to represent each
individual case of severe weather): Event Name, Event Type, Afected
Countries, Location, Country Coordinates (Latitude), Country
Coordinates (Longitude), Event Description. The fields regarding the
event (name, type, location, description) are extracted from the
Wikipedia pages, while the coordinates are retrieved from querying</p>
      <sec id="sec-5-1">
        <title>1http://rda.ucar.edu/datasets/ds627.0/</title>
        <p>the GeoNames 2 API. For the majority of extracted events, country
names are used to reference the spatial extension of an event which
is stored in the “Afected Countries” field. More detailed spatial
information such as city names or state names are stored in the
“Location” field when they are available 3.
2.3</p>
      </sec>
    </sec>
    <sec id="sec-6">
      <title>Evidence Transfer</title>
      <p>
        Evidence transfer [
        <xref ref-type="bibr" rid="ref3">3</xref>
        ] is a deep learning method that incrementally
manipulates the latent representations of an autoencoder
according to external categorical evidence. In the context of evidence
transfer, any categorical variable can be utilised as evidence. The
most straight forward case of evidence is using the outcome of
an auxiliary task. Evidence transfer has been developed with the
notion that in practice the availability of external data is either not
guaranteed, or we may observe the outcome of external processes
without having explicit access to the corresponding dataset. It is
a generic method for combining external evidence in the process
of representation learning. It makes no assumptions regarding the
nature or source of external evidence. It is efective when
introduced with meaningful evidence, robust against non-corresponding
evidence and modular due to its transfer learning nature.
      </p>
      <p>Evidence transfer is a two step method. During the initialisation
step, an autoencoder is trained to reconstruct the input data of the
primary task. To ensure robustness, an intermediate step is required.
During the intermediate step a small biased evidence autoencoder
is trained to reconstruct each categorical evidence source. Evidence
autoencoder is called “biased”, due to introduced limitation in the
amount of iterations. Meaningful evidence is able to converge for
small amount of iterations, leading to a latent projection of the</p>
      <sec id="sec-6-1">
        <title>2https://www.geonames.org</title>
        <p>3Dataset available at: https://github.com/davidath/severe-weather-dataset
evidence, however non-corresponding evidence is not able to
generalise and therefore produce a uniform-like distribution. During
severe weather case detection, we avoid this step, since we know
that textual evidence is retrieved from meaningful sources.</p>
        <p>During the transfer step, the initial latent representation are
manipulated according to external evidence through the joint
optimisation of reconstructing the input, as well as, reducing the cross
entropy between an extended softmax layer of the latent space
and the external evidence. The loss function of the initialisation
step is shown in Equation 1. In Equation 2 we show the evidence
transfer step loss, where  is the set of categorical evidence sources
and  are the extended softmax layers. Structural Similarity Index
(SSIM) is used as reconstruction loss function in order to retain the
structural information of the data.</p>
        <p>1 ∑︁
 =1
ℓ = L ( ,  ′) =
  ( () ,  ′() )
1 ∑︁ 
ℓ   = ℓ +  ∗  =1  ( ,   )
(1)
(2)
2.4</p>
      </sec>
    </sec>
    <sec id="sec-7">
      <title>Class Balancing</title>
      <p>In our experiments, the original data consist of 57584 weather
reanalysis samples in 6-hour increments, while the total amount
of severe weather samples without duplicate dates are only 3136
(less than 6% of the samples). To deal with imbalanced learning,
we experiment with three diferent sampling strategies: (1)
oversampling the minority class, (2) under-sampling the majority class,
(3) combination of over-sampling and under-sampling.</p>
      <p>
        To over-sample the minority class, we use SMOTE [
        <xref ref-type="bibr" rid="ref2">2</xref>
        ]. SMOTE
generates minority class samples by joining the line segments of
k-nearest neighbors. To under-sample the majority class, we
perform random under-sample, although more sophisticated
undersampling methods such as the ENN [
        <xref ref-type="bibr" rid="ref10">10</xref>
        ] (removes data samples that
deviate from the majority of k-nearest neighbors) can also be used.
A combination of both strategies can be achieved by combining
the over-sampling with under-sampling, such as the SMOTEENN
method [
        <xref ref-type="bibr" rid="ref1">1</xref>
        ].
      </p>
      <p>In order to test the efectiveness of each sampling strategy we
experiment with using the primary task of learning representations
to detect severe weather samples by combining all severe cases into
a single class. We manipulate the initial learned space by
incorporating the ground-truth labels (i.e the binary task labels of predicting
severe from non-severe weather samples). Incorporating evidence
that exactly replicates the outcome of the primary task is not
realistic, however we use this scenario in order to investigate the best
choice of sampling strategy without introducing implicit
uncertainty from the choice of external categorical evidence. However,
to test its generalisation, we split the ground-truth into train and
test and only use the evidence labels during training with evidence
transfer.</p>
      <p>Quantitative evaluation with the micro average of precision,
recall and F1 score metrics for the full dataset (train and test) are
presented in Table 2. Our experiments indicate that under-sampling
the majority class is the most fitting for our case. By reducing
the redundancy in the majority class, evidence transfer can more
efectively manipulate the initial representations. It additionally
allows the linear separation into two classes, while during
oversampling and combination, the implicit bias overcomes the latent
space by resulting in a single inseparable cluster.
2.5</p>
    </sec>
    <sec id="sec-8">
      <title>Method overview</title>
      <p>For all of our experiments, we follow the training procedure of
evidence transfer. First, we train a denoising stacked autoencoder
to reconstruct the primary task dataset, i.e. the weather re-analysis
data. The initialisation step is completely unsupervised, no labels
are used during this step. We consider an initial solution to our
primary task, the “baseline” solution, during which we perform
an unsupervised detection method on the initially retrieved
latent representations. We perform the same unsupervised detection
method on incrementally manipulated latent representations from
evidence transfer in order to compare its efectiveness. We supply
the additional evidence sources based on the textual severe weather
dataset.</p>
      <p>During experimental investigation of the best sampling strategy,
one class SVM was used as an unsupervised detection method.
For the cases of detecting individual severe weather cases we use
-means clustering with =2 (prediction of severe or non-severe
weather) as an unsupervised detection method, except a single case
where agglomerative clustering was used instead (ground-truth:
windstorm, evidence: tornado).
3</p>
    </sec>
    <sec id="sec-9">
      <title>EXPERIMENTAL EVALUATION</title>
      <p>We experiment with individually detecting windstorms, floods and
tornado outbreaks 4. We avoid using the hail events due to
limited amount of samples. We rotate between the diferent severe
cases by selecting one case as the ground truth and alternate
between using the rest as external evidence. For example, we select
windstorm weather samples and a portion of non-severe samples
as our primary task – ground truth, while another case, e.g flood,
is selected as the auxiliary task – external evidence. We further
under-sample the remaining non-severe weather cases in order to
match the number of severe weather samples.</p>
      <p>In Table 1, we report experimental results in terms of
precision, recall and F1-score for the anomalous class. Introducing
external evidence leads to linearly separable representations that
increases the efectiveness of clustering, and therefore detecting
the severe weather samples. Even though evidence transfer is a
scalable method that can use multiple sources of evidence, in this
case, it is not as efective, due to ground truth and external evidence
contradicting each other for some portion of the data samples.</p>
      <p>In our experiments, the final dataset consists of non-severe
samples (≈500 after under-sampling to balance the individual severe
class), one severe class as the primary task or ground truth and
one as the external evidence. As an example, consider the task of
predicting windstorm samples as ground truth and the task of
predicting flood samples as external evidence. For the task of predicting
windstorms, non-severe samples and flood samples are labelled as
“normal”. However, for the task of predicting floods, non-severe
samples and windstorm samples are labelled as “normal”. Therefore,
external evidence contradicts the ground truth during non-severe</p>
      <sec id="sec-9-1">
        <title>4Code available at: https://github.com/davidath/severe-weather-detect</title>
        <p>samples. Introducing more sources of external evidence increases
the contradiction for non-severe samples, leading to increased
uncertainty during clustering.</p>
        <p>However, both quantitatively, as shown in Table 1, as well as
qualitatively (ground-truth: windstorm, evidence: flood, depicted
in Figure 3) introducing a single source of evidence improves the
outcome of clustering method by pushing the latent
representations to become linearly separable and therefore improving the
efectiveness for both -means and agglomerative clustering.
4</p>
      </sec>
    </sec>
    <sec id="sec-10">
      <title>FUTURE WORK AND CONCLUSIONS</title>
      <p>In this paper, we investigated using evidence transfer to improve
a primary task of detecting individual cases of severe weather. By
incorporating auxiliary tasks extracted from textual sources, we
efectively manipulated the latent space of an autoencoder using
evidence transfer, in order to increase the efectiveness of severe
weather detection. Making latent representations incrementally
linearly separable resulted in improving the efectiveness of -means
and agglomerative clustering. Additionally, we investigated the best
sampling method for our imbalanced class of detecting severe cases
with non-observable predictors, by evaluating the efectiveness of
evidence transfer in one class SVM (with linear kernel) prediction.</p>
      <p>Future work is directed towards utilising the temporal aspect of
weather re-analysis data. For our experiments, we mostly focused
(a) Baseline of “Windstorm - Flood” Combination
(b) Evidence transfer combination of “Windstorm - Flood”
on using embeddings extracted from an image recognition task.
However, retrieving temporally-aware embeddings from raw data,
e.g. via a recurrent autoencoder, could improve the individual
detection of severe weather cases by exploiting the temporal aspect of
the data. Additionally, since the under-sampling strategy appears
to perform better for this problem, it would be beneficial to increase
the total amount of severe weather samples from additional sources.</p>
    </sec>
    <sec id="sec-11">
      <title>ACKNOWLEDGMENTS</title>
      <p>This work has been supported by the Industrial Scholarships
program of Stavros Niarchos Foundation.</p>
    </sec>
  </body>
  <back>
    <ref-list>
      <ref id="ref1">
        <mixed-citation>
          [1]
          <string-name>
            <surname>Gustavo</surname>
            <given-names>E. A. P. A.</given-names>
          </string-name>
          <string-name>
            <surname>Batista</surname>
          </string-name>
          ,
          <string-name>
            <surname>Ronaldo C. Prati</surname>
          </string-name>
          , and Maria Carolina Monard.
          <year>2004</year>
          .
          <article-title>A Study of the Behavior of Several Methods for Balancing Machine Learning Training Data</article-title>
          .
          <source>SIGKDD Explor. Newsl. 6</source>
          ,
          <issue>1</issue>
          (
          <year>June 2004</year>
          ),
          <fpage>20</fpage>
          -
          <lpage>29</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref2">
        <mixed-citation>
          [2]
          <string-name>
            <surname>Nitesh</surname>
            <given-names>V.</given-names>
          </string-name>
          <string-name>
            <surname>Chawla</surname>
            , Kevin W. Bowyer, Lawrence O. Hall, and
            <given-names>W. Philip</given-names>
          </string-name>
          <string-name>
            <surname>Kegelmeyer</surname>
          </string-name>
          .
          <year>2002</year>
          .
          <article-title>SMOTE: Synthetic Minority over-Sampling Technique</article-title>
          .
          <source>J. Artif. Int. Res</source>
          .
          <volume>16</volume>
          ,
          <issue>1</issue>
          (
          <year>June 2002</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref3">
        <mixed-citation>
          [3]
          <string-name>
            <given-names>Athanasios</given-names>
            <surname>Davvetas</surname>
          </string-name>
          , Iraklis A.
          <string-name>
            <surname>Klampanos</surname>
            , and
            <given-names>Vangelis</given-names>
          </string-name>
          <string-name>
            <surname>Karkaletsis</surname>
          </string-name>
          .
          <year>2019</year>
          .
          <article-title>Evidence Transfer for Improving Clustering Tasks Using External Categorical Evidence</article-title>
          .
          <source>In The International Joint Conference on Neural Networks (IJCNN)</source>
          .
          <source>IEEE</source>
          , 1-
          <fpage>8</fpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref4">
        <mixed-citation>
          [4]
          <string-name>
            <given-names>D. P.</given-names>
            <surname>Dee</surname>
          </string-name>
          ,
          <string-name>
            <given-names>S. M.</given-names>
            <surname>Uppala</surname>
          </string-name>
          ,
          <string-name>
            <given-names>A. J.</given-names>
            <surname>Simmons</surname>
          </string-name>
          ,
          <string-name>
            <given-names>P.</given-names>
            <surname>Berrisford</surname>
          </string-name>
          ,
          <string-name>
            <given-names>P.</given-names>
            <surname>Poli</surname>
          </string-name>
          ,
          <string-name>
            <given-names>S.</given-names>
            <surname>Kobayashi</surname>
          </string-name>
          ,
          <string-name>
            <given-names>U.</given-names>
            <surname>Andrae</surname>
          </string-name>
          ,
          <string-name>
            <given-names>M. A.</given-names>
            <surname>Balmaseda</surname>
          </string-name>
          , G. Balsamo,
          <string-name>
            <given-names>P.</given-names>
            <surname>Bauer</surname>
          </string-name>
          ,
          <string-name>
            <given-names>P.</given-names>
            <surname>Bechtold</surname>
          </string-name>
          ,
          <string-name>
            <surname>A. C. M. Beljaars</surname>
            , L. van de Berg,
            <given-names>J.</given-names>
          </string-name>
          <string-name>
            <surname>Bidlot</surname>
            ,
            <given-names>N.</given-names>
          </string-name>
          <string-name>
            <surname>Bormann</surname>
            ,
            <given-names>C.</given-names>
          </string-name>
          <string-name>
            <surname>Delsol</surname>
            ,
            <given-names>R.</given-names>
          </string-name>
          <string-name>
            <surname>Dragani</surname>
            ,
            <given-names>M.</given-names>
          </string-name>
          <string-name>
            <surname>Fuentes</surname>
            ,
            <given-names>A. J.</given-names>
          </string-name>
          <string-name>
            <surname>Geer</surname>
            ,
            <given-names>L.</given-names>
          </string-name>
          <string-name>
            <surname>Haimberger</surname>
            ,
            <given-names>S. B.</given-names>
          </string-name>
          <string-name>
            <surname>Healy</surname>
            ,
            <given-names>H.</given-names>
          </string-name>
          <string-name>
            <surname>Hersbach</surname>
            ,
            <given-names>E. V.</given-names>
          </string-name>
          <string-name>
            <surname>Hólm</surname>
            ,
            <given-names>L.</given-names>
          </string-name>
          <string-name>
            <surname>Isaksen</surname>
            ,
            <given-names>P.</given-names>
          </string-name>
          <string-name>
            <surname>Kållberg</surname>
            ,
            <given-names>M.</given-names>
          </string-name>
          <string-name>
            <surname>Köhler</surname>
            ,
            <given-names>M.</given-names>
          </string-name>
          <string-name>
            <surname>Matricardi</surname>
            ,
            <given-names>A. P.</given-names>
          </string-name>
          <string-name>
            <surname>McNally</surname>
            ,
            <given-names>B. M.</given-names>
          </string-name>
          <string-name>
            <surname>Monge-Sanz</surname>
            ,
            <given-names>J.-J.</given-names>
          </string-name>
          <string-name>
            <surname>Morcrette</surname>
            ,
            <given-names>B.-K.</given-names>
          </string-name>
          <string-name>
            <surname>Park</surname>
            ,
            <given-names>C.</given-names>
          </string-name>
          <string-name>
            <surname>Peubey</surname>
            , P. de Rosnay,
            <given-names>C.</given-names>
          </string-name>
          <string-name>
            <surname>Tavolato</surname>
            ,
            <given-names>J.-N.</given-names>
          </string-name>
          <string-name>
            <surname>Thépaut</surname>
            , and
            <given-names>F.</given-names>
          </string-name>
          <string-name>
            <surname>Vitart</surname>
          </string-name>
          .
          <year>2011</year>
          .
          <article-title>The ERA-Interim reanalysis: configuration and performance of the data assimilation system</article-title>
          .
          <source>Quarterly Journal of the Royal Meteorological Society</source>
          <volume>137</volume>
          ,
          <issue>656</issue>
          (
          <year>2011</year>
          ),
          <fpage>553</fpage>
          -
          <lpage>597</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref5">
        <mixed-citation>
          [5]
          <string-name>
            <surname>Iraklis</surname>
            <given-names>A.</given-names>
          </string-name>
          <string-name>
            <surname>Klampanos</surname>
            , Athanasios Davvetas, Spyros Andronopoulos, Charalambos Pappas, Andreas Ikonomopoulos, and
            <given-names>Vangelis</given-names>
          </string-name>
          <string-name>
            <surname>Karkaletsis</surname>
          </string-name>
          .
          <year>2018</year>
          .
          <article-title>AutoencoderDriven Weather Clustering for Source Estimation during Nuclear Events</article-title>
          .
          <source>Environmental Modelling &amp; Software</source>
          <volume>102</volume>
          (
          <year>April 2018</year>
          ),
          <fpage>84</fpage>
          -
          <lpage>93</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref6">
        <mixed-citation>
          [6]
          <string-name>
            <surname>Mikhail</surname>
            <given-names>A.</given-names>
          </string-name>
          <string-name>
            <surname>Krinitskiy</surname>
          </string-name>
          , Yulia A.
          <string-name>
            <surname>Zyulyaeva</surname>
          </string-name>
          , and
          <string-name>
            <surname>Sergey</surname>
            <given-names>K.</given-names>
          </string-name>
          <string-name>
            <surname>Gulev</surname>
          </string-name>
          .
          <year>2019</year>
          . (9
          <year>2019</year>
          ). https://doi.org/10.6084/m9.figshare.
          <volume>9851099</volume>
          .v1
        </mixed-citation>
      </ref>
      <ref id="ref7">
        <mixed-citation>
          [7]
          <string-name>
            <given-names>T. N.</given-names>
            <surname>Krishnamurti</surname>
          </string-name>
          ,
          <string-name>
            <given-names>K.</given-names>
            <surname>Rajendran</surname>
          </string-name>
          ,
          <string-name>
            <given-names>T. S. V.</given-names>
            <surname>Vijaya</surname>
          </string-name>
          <string-name>
            <surname>Kumar</surname>
          </string-name>
          , Stephen Lord, Zoltan Toth, Xiaolei Zou, Steven Cocke,
          <string-name>
            <given-names>Jon E.</given-names>
            <surname>Ahlquist</surname>
          </string-name>
          ,
          <string-name>
            <given-names>and I. Michael</given-names>
            <surname>Navon</surname>
          </string-name>
          .
          <year>2003</year>
          .
          <article-title>Improved Skill for the Anomaly Correlation of Geopotential Heights at 500 hPa</article-title>
          .
          <source>Monthly Weather Review</source>
          <volume>131</volume>
          ,
          <issue>6</issue>
          (
          <year>2003</year>
          ),
          <fpage>1082</fpage>
          -
          <lpage>1102</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref8">
        <mixed-citation>
          [8]
          <string-name>
            <given-names>John</given-names>
            <surname>Michalakes</surname>
          </string-name>
          , Jimy Dudhia,
          <string-name>
            <given-names>D.</given-names>
            <surname>Gill</surname>
          </string-name>
          , Tom Henderson,
          <string-name>
            <given-names>J.</given-names>
            <surname>Klemp</surname>
          </string-name>
          ,
          <string-name>
            <given-names>W.</given-names>
            <surname>Skamarock</surname>
          </string-name>
          , and
          <string-name>
            <given-names>Wei</given-names>
            <surname>Wang</surname>
          </string-name>
          .
          <year>2004</year>
          .
          <article-title>The Weather Reseach and Forecast Model: Software Architecture and Performance</article-title>
          .
          <source>11th ECMWF Workshop on the Use of High Performance Computing In Meteorology.</source>
        </mixed-citation>
      </ref>
      <ref id="ref9">
        <mixed-citation>
          [9]
          <string-name>
            <given-names>Murat</given-names>
            <surname>Türkeş</surname>
          </string-name>
          , UM Sümer, and
          <string-name>
            <given-names>G</given-names>
            <surname>Kiliç</surname>
          </string-name>
          .
          <year>2002</year>
          .
          <article-title>Persistence and periodicity in the precipitation series of Turkey and associations with 500 hPa geopotential heights</article-title>
          .
          <source>Climate Research - CLIMATE RES 21 (05</source>
          <year>2002</year>
          ),
          <fpage>59</fpage>
          -
          <lpage>81</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref10">
        <mixed-citation>
          [10]
          <string-name>
            <given-names>D. L.</given-names>
            <surname>Wilson</surname>
          </string-name>
          .
          <year>1972</year>
          .
          <article-title>Asymptotic Properties of Nearest Neighbor Rules Using Edited Data</article-title>
          .
          <source>IEEE Transactions on Systems, Man, and Cybernetics SMC-2</source>
          ,
          <issue>3</issue>
          (
          <year>1972</year>
          ),
          <fpage>408</fpage>
          -
          <lpage>421</lpage>
          .
        </mixed-citation>
      </ref>
    </ref-list>
  </back>
</article>