Global-Scale Synchronization in the Meteorological Data: A Vectorial Analysis That Includes Higher-Order Differences

: To examine the evidence of global warming, in recent years, there has been a growing interest in the statistical analysis of time-dependent meteorological data. In this paper, for 116 observational stations in the world, sequential variations of the monthly distributions of meteorological data are analyzed vectorially. For specific monthly data, temperatures and precipitations are chosen, both of which are averaged over three decades. Climate change can be revealed through the intersecting angle between two 33-dimensional vectors being composed with monthly mean values. Subsequently, the angle data for the entire stations are analyzed statistically and compared between the former (1931 – 1980) and the latter (1951 – 2010) periods. Irrespective of the period and the hemisphere, the variation of the angles is found to show the exponential growth as a function of their latitudes. Furthermore, consistent with other studies, this trend is shown to become stronger in the latter period, indicating that the so-called snow/ice-albedo feedback occurs. In contrast to the temperatures, for the precipitations, no significant correlation is found between the angle and the latitude. To examine the albedo effect in more detail, a regional analysis for 75 stations in Japan is carried out as well. Numerical results show that the effect is significant even for the relatively narrow latitudinal range (19%) of the hemisphere. Finally, a synchronization of the monthly patterns of temperatures is given between the northern district of Japan and both North America and Eastern Europe.


Introduction
An agreement has been reached on the signs of global warming that have been found since the latter half of the former century [1,2]. Indeed, both the atmosphere and ocean have been getting warmer, the amounts of ice and snow have been decreasing, and the sea level has continued to rise. It has been shown that these were caused by the accumulating greenhouse gases due to human activity. In 2015, the concentration of carbon dioxide in the Earth's atmosphere recorded 400 ppm [3], and at the same time, the mean temperature of the world became +1 °C higher than before the Industrial Revolution [4]. According to an estimation with the Representative Concentration Pathways (RCP) 8.5 scenario [2], it is reported that the difference attains to +2 °C in the early period of the present century, and furthermore, it attains up to +3 °C in the latter period of the century. In December of the same year, to present a long-term target for suppressing global warming, the Paris Agreement was adopted, and became effective in November of the subsequent year, when the Intergovernmental Panel on Climate Change (IPCC) had a large bearing upon its realization. Locally, warming may be accelerated at a rate higher than the estimation. Furthermore, it is predicted that so-called extreme meteorological phenomena, such as record-breaking temperatures as well as heavy rains will continue to be stronger in their scale and occur more frequently. Disasters due to climate change have started to be detected. For example, only recently has the Australian government reported the extinction of mammals similar to rats possibly because of the sea-level rising on an islet in the Great Barrier Reef. The increasing sea temperature is causing serious damage to the coral reef. Note that, along with the Marshall Islands in the east of the West-Pacific Micronesia, both the Great Barrier Reef along the east coast of Australia and the reef in the vicinity of the Okinawa Island of Japan are now in a critical condition. In the latter case, the coral bleaching has been arising from the rising sea temperature. Although the damage may not be so critical as that responsible for the extinction of an endangered species, reports of heatstroke are increasing substantially in the summer of Japan, where wet heat waves exceeding 35 °C continue sporadically for many days. To reduce the mortality risk of runners, the Tokyo 2020 Olympic Games Organizing Committee set the starting time of the marathon race at 6:00 a.m., which is substantially earlier than in the 1964 Tokyo Olympics. Furthermore, medical advisors suggested starting the race even earlier (5:30 a.m.) [5].
In recent years, not to mention the conventional fields such as geophysics, meteorology, climatology, aerometry, glaciology, and oceanography , a large number of research articles concerning climate change and/or global warming have also been published in the context of statistical physics [50][51][52][53][54][55][56][57][58][59][60][61][62][63][64][65]. The methodologies worth noting range over the method of detrended fluctuation analysis (DFA) [51], a variationally optimized Markov analysis combined with a successive over-relaxation algorithm [52], Monte Carlo simulations [53,56,57], feedforward artificial neural networks [54], convergent cross mapping (CCM) [59], wavelet transformation method [60], the method using Mahalanobis distance metrics [61], multifractal temporally weighted detrended fluctuation analysis [63], horizontal multivariate singular spectral analysis techniques [64], and the formalism of pitchfork bifurcation with the Lyapunov exponent [65]. In this paper, for 116 observational stations in the world (97 in the Northern Hemisphere and 19 in the Southern Hemisphere, as dotted on the map in Figure 1), sequential variations of the monthly distributions of meteorological data are analyzed vectorially. For specific monthly data, precipitations as well as temperatures are chosen, both of which are averaged over thirty years. According to the chronology [66][67][68] the entire eighty years are segmented into three periods: I: 1931-1960, II: 1951-1980, and III: 1981-2010. Climate change can be revealed through the cross angle between two twelve-dimensional vectors being composed with monthly mean values. In addition to the World temperature, analysis is made for the precipitation. Subsequently, in order to discuss the albedo feedback in detail, a regional analysis is carried out. The snow/ice-albedo effect can be explained by a positive feedback of enhanced light absorption due to the degenerate cover of snow and ice [1,2,69]. Specifically, 75 observational stations in Japan, which cover from Naze (28°23′ N, 129°30′ E) to Wakkanai (45°25′ N, 141°41′ E), are considered. While the longitudinal range of the spots approximately covers Darwin (130°50′ E) to Cairns (145°47′ E), latitudinally, the range is confined in Montreal (45°28′ N) to Cape Canaveral (28°23′ N) or in Lyon (45°43′ N) to Las Palmas de Gran Canaria (28°06′ N), the Islas Canarias. It is interesting to investigate whether the effect is significant even for the relatively narrow latitudinal range (19%) of the hemisphere. Finally, a synchronization of the monthly patterns of temperatures is given between selected spots in Japan and those in North America as well as Eastern Europe.

Methodology
In this paper, to enhance the precision and to avoid the production of a number of tie data (data having the same value), additional 21 components consisting of eleven first-order and ten second-order differences are included in each vector composed with monthly mean values. With this extension, one eventually obtains the 33-dimensional vector. Besides the simple linear regression [70][71][72], in order to examine the sign of the curvature, nonlinear regressions [73] are adopted. Identifying the (plus or minus) sign of curvature of the data as a function of the latitude is of particular concern for finding evidence of the warming.

Modified Vectorial Analysis
We consider two vectors p and q being composed by twelve components where <ui> and <xi> (i = 1, 2, …, 12, for January to December) represent monthly mean values of the meteorological quantities such as temperatures and precipitations (for specific numeric, see, e.g., Equations (C1) and (C2) in Appendix C of Supplementary Material); the values are averaged further over the three decades. The intersecting angle between the vectors can be obtained by where p·q indicates the scalar product between p and q. For instance, for p and q being produced with monthly averaged temperature distributions in the two subsequent periods (Period I to II or Period II to III), climate change between the periods can be evaluated by the magnitude of the angle. This procedure, however, is not free from disadvantage, because a number of quasi-tie data are generated, which might be responsible for making the subsequent statistical analysis ill-conditioned.
Here, j = 1, 2, …, 11 and k = 1, 2, …, 10. Note that <vj > and <yj> represent the rate of change whereas <wk> and <zk> imply the 'monthly variation of curvature.' In comparison between the monthly patterns bearing a close resemblance, besides the first-order difference the inclusion of the second-order difference is useful for discriminating their details. Although an unlimited extension to joining much higher-order differences is possible, no one can find the physical meaning of the numerical artifacts. In this paper, we concentrate on the relatively slowly varying patterns and do not deal with chaotic sequences as seen in the climate oscillations. In Equations (4) and (5), the independence among components is no longer ensured. It should be noted here that, as our method does not require the calculation of the eigenvector, the loss of independence is not as important as in the principal component analysis (PCA) [70], in which it is necessary to strictly preserve the orthogonalities among components of a vector. Upon the substitution of Equations (4) and (5) into Equation (3), one can improve the method substantially. Eventually the climate change between two periods as well as between two stations can be detected by calculating the rotation angle between the vectors p and q defined in the 33-dimensional co-ordinate. It seems that the present method with the rotation of the angle is useful for revealing potential synchronizations between meteorological data, along with the one using the Minkowski distance Dr between the two vectors; with a parameter r the distance is reduced to the city-block (L1, taxicab, Manhattan, or more generally Bray-Curtis, Canberra, and Lance-Williams), the Euclidean, and the Chebyshev distance, respectively, for r = 1, 2, and ∞. Only recently have the Minkowski distance functions been used to detect teleconnections between temperatures, aerosol optical paths, and precipitation time series [74].

Statistical Analysis
Subsequently the angle data θi (i = 1, 2, …, n; n being the number of observational stations) will be dealt with statistically. Below, as specific regressions on the angles versus the coordinates, three functions are chosen: Logarithmic regression: θ= a + b lnχ, where ln abbreviates the natural logarithm; χ represents the coordinate of the station; a and b are real constants to be determined with the least square fit, specifically a = mY -b mX and b = sXY/sX 2 ; (X, Y) = (χ, θ), (χ, lnθ), and (lnχ, θ) for Equation (8), (9), and (10), respectively. Here, mY (mX), sXY, and sX 2 , respectively, are the arithmetical mean of Y (X), the covariance of X and Y, and the variance of X. It should be emphasized here that the regressions can deal not only with the simple linearity (Equation (8)) but with the growth (Equation (9) with a positive curvature) as well as the saturation model (Equation (10) with a negative curvature). The validity of the present regression model can be checked by using the degree of fit, |r|, with the Pearson's formula r = sXY/ (sX sY ) (0 < |r| < 1; sY 2 represents the variance of Y.), together with the Durbin-Watson ratio, d (0 < d < 4) [70,72].
with the summation Σ for i = 1 to n -1; with the summation Σ for i = 1 to n. Here, n indicates the summation of regression points; the hat on Y represents the point on the regression line. With level alpha test being done, it can be judged that if 0 < d < dL (dL is the lower critical value) there exists a positive correlation between adjacent points on the sequence of the residual data ei (i = 1, 2, … , n) and that if dU < d < 2 (dU is the upper critical value) there is no correlation between them. Note that, for dL < d < dU, any judgment is impossible. Therefore, a null hypothesis that 'H: there is a correlation between the neighboring residual data' is rejected solely for dU < d < 2. (For d > 2, d must be replaced by 4 -d.) For typical levels the two critical values are obtainable from numerical tables available. In this paper the level 1% test, i.e., α= 0.01, will be adopted.
In what follows, the angle data generated for the entire stations on the map of Figure 1 are analyzed and are compared between the former (I/II: 1931-1980) and the latter (II/III: 1951-2010) period. If the trend obeys the exponentially growing curve with a positive curvature in the plot of the divergence versus the latitude, and at the same time its trend is shown to become stronger in the latter period, one can give hard evidence to show that the so-called snow/ice-albedo effect due to the positive feedback [1,2,69] may be responsible for the global warming.

Data
All the meteorological data employed in this paper are taken from a series of the Chronological Scientific Tables [66][67][68], which have over more than ninety years been edited by the National Astronomical Observatory, Japan. The observation points in the Tables were sampled without bias by the editors using Observing Stations, published by the World Meteorological Organization (WMO). By the Japan Meteorological Agency, both the temperature and the precipitation data are averaged over thirty years. Note that the original data were cited from those of the Global Historical Network, which were delivered by the National Climate Data Center (NCDC), USA, and from those of meteorological messages for monthly climatic variations, which had been received from meteorological organizations located all over the world. When necessary, data are supplemented with those of the National Center for Atmospheric Research (NCAR), USA. In Japan, the Tables have been winning broad support not only from researchers, engineers, and teachers, but from students and others engaged in science. Over the series of the Chronological Scientific Tables, the entire observational term is segmented into three periods: I   [66], II (1951II ( -1980 [67], and III (1981-2010) [68]. Note that because of an uncertain editorial reason the former two periods overlap for 1951-1960; data of 1921-1950 are not available. In principle a complete survey is conducted and neither random nor intentional sampling is employed. The World observational stations are listed in Tables A1 and A2 in Appendix A of Supplementary Material, along with their coordinates and heights. In the Northern Hemisphere (Table S1), the latitudes range from 01°22′ N (Singapore) to 73°30′ N (Ostrov Dikson, Russia). Of the 97 stations, only the two spots (2%), #1: Ostrov Dikson and #2: Barrow, are located within the Arctic Circle (66°33′39″ N), whereas the 19 spots (20%), #79: Kolkata to #97: Singapore, are confined in the region between the equator and the Tropic of Cancer (23°26′21″N). In the Southern Hemisphere (Table S2), the latitudes cover 01°19′ S (Nairobi, Kenya) to 54°48′ S (Ushuaia, Argentina). Of the 19 stations, 8 (42%), #12: Rio De Janeiro to #19: Nairobi, are in the region between the equator and the Tropic of Capricorn (23°26′21″ S); there is no station within the Antarctic Circle (66°33′39″ S). The longitudinal distribution is (E, W) = (69, 28) and (12,7), respectively. For the unit of temperature, the Celsius scale is adopted.

Global Analysis
For monthly average temperature data, the dependence of the rotation angle on the latitude is plotted in Figure 2 for the Northern Hemisphere and in Figure 3 for the Southern Hemisphere. The top-thirty rankings of the rotation angle are given, respectively, in Tables 1 and 2, and in Table 3 and Table 4. For monthly averaged precipitation data, the top-thirty rankings of the angle are listed in Tables     to Period II . The curve superimposed in the plots represents the optimal exponential regression (r = 0.738 with d = 1.788 > dU = 1.55 for n = 97). Note that the Urumchi's angle θ = 17.80° on the top in Table 1 is truncated. (b) Rotation from Period II  to Period III . The curve in the plots represents the optimal exponential regression (r = 0.843 with d = 1.089 < dL = 1.51 for n = 97).
ΔT indicates the increment of the mean temperature from the former to the latter period.
the best fit (r, d) = (0.819, 1.406) for c = 1.7, exhibiting a power-law growth, though the magnitude of d remains smaller than dL.) The stations that have raised ranking substantially (bold-faced in Table 2) are #38→#2: Anchorage (61°09′ N), #26→#6: St. Petersburg (59°58′ N), #34→#8: Moskva (55°50′ N), #88→#17: Muenchen (48°21′ N), #65→#20: Vancouver (49°11′ N), #80→#21: Denver (39°46′ N), and #63→#30: Bucuresti (44°30′ N), most of which are sited on the relatively high latitude, implying the enhanced snow/ice-albedo feedback due to the global warming. Autocorrelation between the rotation for the latter (Period II-III) and that for the former (Period I-II) has been found to be r = 0.764 and d = 1.710 > dU = 1.55 (n = 97) with the best logarithmic fit. (3) From the results plotted on Figure 3 and from those listed in Tables 3 and 4, it is found that although the effect due to the warming is not as strong as in the northern counterpart the conclusion for the Northern Hemisphere is also applicable to the Southern Hemisphere. Autocorrelation between the rotation for the latter (Period II-III) and that for the former (Period I-II) has been found to be r = 0.843 and d = 1.806 > dU = 1.13 (n = 19), but with the best exponential fit, showing contrast to the best logarithmic fit seen for the northern counterpart. (4) In contrast to the temperature, for the World precipitation, independently of the hemisphere as well as the period, no significant correlation has been seen in the latitudinal dependence. For the respective sequence the top-thirty rankings of the rotation angle are given in Tables   Here, one sees a series of nulls. The stations that have raised the ranking considerably (bold-faced in Table S4)

Regional Analysis
In a regional analysis, we concentrate on 75 observational stations in Japan [66][67][68], the co-ordinates of which are given in Table B1 in Appendix B of Supplementary Material. For monthly average temperature data, the dependence of the rotation angle on the latitude is plotted in Figure 4a for Period I  to Period II  and in Figure 4b for Period II (1951II ( -1980 to Period III (1981-2010). The top-twenty rankings of the rotation angle are given, respectively, in Tables 5, 6. In addition to the latitudinal dependence, for the respective sequence, the longitudinal dependence of the angle is shown in Figure 5. Moreover, to discuss the snow-albedo effect in more detail, data of snow cover are given in Figure 6 as a function of the latitude. Lastly, top-twenty Japanese stations in the rotation angle of the monthly distributions of precipitations are listed in Tables B2, B3 in Appendix B of Supplementary Material. With these results, one can conclude as follows: (1) In Figure 4a, which plots the rotation between Period I  and Period II , in the same way as its global counterparts (Figures 2a and 3a) the best fit has been found to be exponential but with r = 0.642 and d = 1.867 > dU = 1.50 (n =75). (cf. for the linear regression of Equation (8), r = 0.597 with d = 1.876.) It can be seen from Table 5 that the top-twenty ranking is occupied by the stations of the northern areas being located on Hokkaido (#2 -#8, and #14) and the Tohoku Region (#10, #13, and #15), as well as the Nagano Prefecture (#1, #9, and #18).
Moreover, it is interesting to note that besides #4: Sapporo and #15: Sendai there are two megacities, #11: Tokyo and #19: Osaka, on that ranking. -   Table 6 that the occupation by the three cold districts remains high, namely, there are nine stations in Hokkaido (#1-#7, #9, and #20), five stations in the Tohoku Region (#8, #12 -#14, and #18), and one station in the Nagano Prefecture (#10). Autocorrelation between the rotation for the latter (Period II-III) and that for the former (Period I-II) has been found to be r = 0.687 and d = 1.322 < dL = 1.45 for n = 75 with the best linear fit.     exponentially on their longitudes, with (r, d) = (0.831, 1.402). In a comparison between Figure 5a, b, it should be noted that, unlike the observation on the comparison between Figure 4a, b, no increase of the gradient of the regressive curve is seen in the longitudinal dependence of the angle, concluding that from both plots, effects due to the snow/ice-albedo feedback are not detectable longitudinally.  Figure  6a and that of Figure 6b shows that in Japan the latitudinal dependence of snow cover decreases gradually. Note that the gradient of the regressive line on Figure 6b is 96% that on Figure 6a.  Table S8 and Table S9 in Appendix B of Supplementary Material. In Table B3    ----------------------------------------

Toward Detecting Synchronization
The global change of climate may give rise to the so-called teleconnection [75][76][77][78], the unexpected synchronization of meteorological data among faraway spots on Earth. Here, in contrast to the traditional teleconnection studies in which either the phase synchronization or the synchronized chaos was of main concern, we regard the two spots as climatically synchronous if the cross angle between the two monthly mean vectors over the three decades decreases monotonically over the entire period, I to III. Of the 97C2 (= 4656) combinations possible on the Northern Hemisphere, it has been found that the monthly average temperature patterns on spots in the northern district of Japan and those in North America as well as in Eastern Europe exhibit a synchronous behavior. Main results are summarized in Tables 7, 8. Here, five Japanese stations dispersed in Table S1 Tables 7 and 8, we remark as follows: (1) First, the complete sequences of the cross angle of the monthly distributions of temperature between three North American stations (#37: Boston, #38: Chicago, and #46: Denver in Table S1) and the five Japanese stations are summarized in Table 7. It appears that synchronization arises   Table 7. The sequence of the cross angle of the monthly distributions of temperatures between three North-American and five Japanese stations.
Ｐｅｒｉｏｄ Ⅰ Ｐｅｒｉｏｄ Ⅱ Ｐｅｒｉｏｄ Ⅲ - Ｐｅｒｉｏｄ Ⅰ Ｐｅｒｉｏｄ Ⅱ Ｐｅｒｉｏｄ Ⅲ - solely at Sapporo, which shows striking contrast to the decreasing synchronization at the other three (Tokyo, Osaka, and Fukuoka). (2) Subsequently the results between three East-European stations (#20: Praha, #31: Bucuresti, and #36: Sofia in Table S1) and the five stations in Japan are given in Table 8. The tendency similar to North America (Table 7) can be seen for Eastern Europe.  Table 8. The sequence of the cross angle of the monthly distributions of temperatures between three East-European and the five Japanese stations.
Ｐｅｒｉｏｄ Ⅰ Ｐｅｒｉｏｄ Ⅱ Ｐｅｒｉｏｄ Ⅲ - Ｐｅｒｉｏｄ Ⅰ Ｐｅｒｉｏｄ Ⅱ Ｐｅｒｉｏｄ Ⅲ - . Although it appears that, at least in the six combinations, the two temperature distributions are no doubt in synchronization, the effect is found to be much weaker than that revealed in the northern counterpart.

Global Analysis
Comparison among Figures 2 and 3 indicates quantitatively that the climate change in recent years has occurred globally independently of the hemisphere, but the situation is far more apparent in the northern countries on the Northern Hemisphere than the countries on the Southern Hemisphere. It should be stressed that comparison between Figures 2a and 2b reproduces the conclusion of former studies that owing to the snow/ice-albedo feedback [63,[79][80][81][82] the regions with the higher latitudes on the Northern Hemisphere are most susceptible of the warming [83,84]. The effect can be explained by an increasing light absorption arising from degenerating cover of snow and ice that show, in comparison with the seawater and the ground, relatively higher reflectivity (at least 85 to 90% for snow whereas 65% for ice) of sunlight [1,2,69]. There are three reasons why the warming is more intense in the Northern Hemisphere (Figure 2) than the Southern Hemisphere ( Figure 3): (1) The emission of major greenhouse gases such as the carbon dioxide, methane, ozone, and chlorofluorocarbons (CFCs) is unevenly distributed toward the former hemisphere, particularly in the region over 30° N to 60° N where there are many industrialized countries. The emission from the latter hemisphere cannot bear comparison with this. (2) As is well-known the distribution of the Earth's continents is unevenly distributed toward the Northern Hemisphere where the land area amounts to 1.0 × 10 10 ha, whereas that of the Southern Hemisphere is 4.9 × 10 9 ha. Note that the land has a specific heat lower than the sea. (3) In striking contrast to the critical circumstance in the Arctic, the temperature growth on the Antarctica is not as steep as the northern counterpart. The reason is that the extremely thick ice (2.0 to 2.5 km) on the Antarctic Continent is free from the exposure of the ground, making ineffective the positive feedback due to the reduced albedo (i.e., the enhanced absorption) of the surface. Note again that the optical reflectivity for snow is at least 85% to 90%.

Regional Analysis
It should be emphasized that comparison between Figure 4a and Figure 4b reproduces the conclusion of the global analysis, i.e., owing to the positive feedback due to the decreased albedo the cold districts with relatively high latitude in Japan are more susceptive of the global warming. Again, the effect can be attributed to an increasing light absorption arising from degenerating cover of snow that shows high reflectivity of sunlight [1,2,69]. As in Europe and in North America, Japan has often been hit by the subsequent record-breaking heat waves. Anomalous weather phenomena such as a "tropical night," a night when the temperature does not fall below 25 °C, have become too frequent to make the headlines on newspapers. To declare a state of emergency on spots where the temperature exceeds 35 °C, besides a conventional "tropical day" on which the temperature is 30 °C or above, the Japan Meteorological Agency invented a "day of intense heat." In addition to the intensifying strength of the chronic heat wave in the midsummer, to date there have been reports of disasters in connection with the climate crisis: (1) A recent report of numerical simulations by the Okhotsk Sea Ice Museum of Hokkaido [85] has predicted that in the near future there will be no sea ice drifting from a north region of the Sea of Okhotsk to the northern coast of Hokkaido, the northernmost island of Japan. The sea ice is produced in the sea of northeastern area of Sakhalin, and the south drift is caused either by the East Sakhalin Current or by an Okhotsk anticyclone that produces an air mass [86,87].
According to the report [84] the temperature around the northern area of the Sea has risen two degrees over the last fifty years, the increment of which amounts to the three times larger than the world average. Eventually, the rise in temperature has resulted in the 30% shrinkage of the ice area over the entire Sea. It should be noted here that the freezing temperature of the sea ice is −1.8 °C [88]. Because the flow in the vicinity of the southern border of the sea ice is much susceptive of the warming, one can regard the shrinkage as a high-sensitivity sensor for the global warming. Indeed, the moving average of the temperatures on Abashiri (see Rank #2 in Tables 5 and 6), a spot on the coastland of the Sea, versus the sea-ice force off the coast of the spot has demonstrated the negative correlation over the last thirty years [89].  Here it should be highlighted again that the coefficient of regression on the latter is 96% that on the former, demonstrating the influence due to the warming. Moreover, a comparison between the two plots shows that on the spot of 45° N the period of snow cover shortens for eight days in Figure 6b. Indeed, every year, ski resorts in Japan are hit by the snow shortage. Furthermore, the quality of their snow is degrading. To discuss this concern, we note that a pattern of atmospheric pressure in which the high-pressure area lies to the west and the low-pressure area to the east is characteristic of winter of the Japanese Islands. The physical mechanism of the pressure configuration can be explained as follows: First, because of the radiative cooling, temperatures in winter fall substantially on the inland areas of Siberia and China. Second, the enhanced density of the cooled air generates a descending air current, which is responsible for the emergence of the high pressure in the west of Hokkaido. Third, collision of two air masses coming from the North Pole and from the equator produces a large-scale stable vortex in the vicinity of the Kuriles and the southern area of the Kamchatka Peninsula; eventually the vortex is responsible for observing the low pressure in the east of Hokkaido. The warming as shown especially in the relatively high latitudes of the Northern Hemisphere (see Figure 2) is responsible for reducing the radiative cooling, i.e., reducing the strength of a seasonal wind from the Eurasian Continent to Hokkaido. In Table 2, one can confirm the reason, for instance, on Rank #13: Vladivostok (43°07′ N), #14: Irkutsk (52°16′ N), and #23: Chang-chun (43°54′ N). In response to the shrinkage of the Siberian cold mass moving southeast from the Maritime Territory (Primorsky Kray) of Russia to the Hidaka mountains, the world-renowned powdery snow on resorts in Hokkaido threatens to be in a critical condition.

Toward Detecting Synchronization
Of the combinations possible on the Northern Hemisphere, it has been found that the monthly averaged temperature patterns on spots in the northern district of Japan and those in North America (Table 7) as well as in Eastern Europe (Table 8) exhibit a synchronous behavior. First of all, it should be noted that Japan consists of the four major islands, i.e., Hokkaido, Honshu, Shikoku, and Kyushu. Of these, Hokkaido is separate from Honshu via the Blakiston's line on the midst of the Tsugaru Strait, across which the climate as well as the ecosystem change substantially. The line was found by a British zoologist and soldier, Thomas Wright Blakiston (1832-1891), who determined a zoogeographic border between the northern Siberian and the southern Manchurian subzone. While the line gives the northern limit of Asiatic black bears, Japanese monkeys, and Japanese squirrels, it gives the southern limit of brown bears, Asiatic pikas, and Eurasian striped squirrels. Specifically, with the Köppen's climatic classification, the former island (Hokkaido) belongs to Dfa whereas the latter three (Honshu, Shikoku, and Kyushu) belong to Cfa. Consequently, in both vegetation and soil, there are qualitative differences between the two regions across the strait. First, in Hokkaido there is no long spell of rainy weather in early summer, but in winter the surface is covered with snow. Second, the northernmost island is seldom visited by a typhoon because, as it propagates northward, the tropical cyclone generated in the vicinity of the Philippine Islands is converted gradually to an extratropical storm with its energy much more dispersed than the cyclone. With respect to features peculiar to the land, Sapporo, being located in Hokkaido, and the stations in North America and Europe share something in common in that they have a cool temperature climate. Here, we note that the Köppen's climatic classifications are Dfa for Boston, Chicago, and Bucuresti, with Cfb for Praha and Sofia and BSk for Denver.

Conclusions
For 116 observational stations in the world, sequential variations of the monthly distributions of meteorological data have been analyzed. For specific monthly data, precipitations as well as temperatures have been chosen, both of which are averaged over thirty years. Climate change can be detected through calculating the intersecting angle between two twelve-dimensional vectors being composed with monthly mean values. To enhance the precision and to avoid the production of a number of tie data, additional 21 components consisting of eleven first-order and ten second-order differences have been incorporated in each vector. With this extension, one eventually obtains the 33-dimensional vector. Subsequently the angle data for the entire stations have been analyzed statistically and been compared between the former  and the latter (1951-2010) period. With nonlinear as well as linear regressions, irrespective of the period and the hemisphere, the variation of the rotation angles has been found to show the exponential growth as a function of their latitudes. Furthermore, this trend has been shown to be stronger in the latter period, indicating that the snow/ice-albedo feedback gets stronger. In contrast to the World temperature, for the precipitation, no significant correlation has been found between the angle and the latitude. In order to examine the albedo feedback in more detail, a regional analysis has been carried out as well. Specifically, 75 stations in Japan, which range from 28°23′ N to 45°25′ N, have been chosen. Numerical results have shown that the warming is significant even for the relatively narrow range (19%) of the hemisphere. Finally, a synchronization of the monthly patterns of temperatures has been given between the northern district of Japan and both North America and Eastern Europe. To the author's knowledge, the methodology given in this paper is original and bears a comprehensive applicability also to the analysis of meteorological data other than the temperature and the precipitation. With this method, unknown synchronizations could be revealed. To conclude, the combination of vectorial and statistical analysis has allowed one to reveal the world-scale impendence distinctly.
In an essay published in 1957 [90], Ukichiro Nakaya , who is known worldwide for his achievement of artificial snow, issued a warning as a precursor of global warming. Along with earlier studies, including the IPCC's reports, and with the predictions by Earth simulators running all over the world, our results presented in this paper have demonstrated statistically that his presage was far from utterly unfounded. For global warming, we are now in a position to prepare to relieve the damage on the assumption that planet Earth will continue to become warmer, slowly but surely.
Supplementary Materials: The following are available online at www.mdpi.com/xxx/s1, Table S1: Observational stations on the Northern Hemisphere, Table S2: Observational stations on the Southern Hemisphere, Table S3: Top-thirty Northern Hemispheric stations in the rotation angle of the monthly distributions of precipitations from Period I  to II (1951II ( -1980, Table S4: Top-thirty Northern Hemispheric stations in the rotation angle of the monthly distributions of precipitations from Period II  to III (1981-2010), Table S5: The Ranking of Southern Hemispheric stations in the rotation angle of the monthly distributions of precipitations from Period I  to II (1951II ( -1980, Table S6: The Ranking of Southern Hemispheric stations in the rotation angle of the monthly distributions of precipitations from Period II  to III (1981-2010), Table S7: Observational stations in Japan. , Table S8: Top-twenty Japanese stations in the rotation angle of the monthly distributions of precipitations from Period I  to II (1951II ( -1980, Table S9: Top-twenty Japanese stations in the rotation angle of the monthly distributions of precipitations from Period II  to III (1981-2010).

Author Contributions:
The author confirms being the sole contributor of this work and approved it for publication.