Investigating the Apparent Teleconnection between Cosmic-Ray Muon Flux over Jamaica and Large-Scale Climate Phenomena That Impact Caribbean

The apparent teleconnection between cosmic-ray muon flux over a base point in the Caribbean is discussed against the background of an extensive record of indices representing large-scale climatic phenomena, but limited cosmic-ray muon flux data. Many investigators have shown that large-scale climate phenomena influence sub-seasonal and seasonal climate variability, especially in the northern hemisphere and their impacts on the Caribbean are well documented. These climatic phenomena that impact the Caribbean include, but are not limited to, the El Nino Southern Oscillation, the Quasi-Biennial Oscillation, the North Atlantic Oscillation, and the Arctic Oscillation which is now being investigated. Although strong statistical correlation between variables over non-contiguous regions are not absolute as proof of teleconnections, the correlation strength can be used as an indication of its existence. The data gathered at the Mona Campus of the University of the West Indies, in Jamaica, using a simple QuarkNet 6000 muon detector over the period September 2011 to September 2013, showed an apparent significant relationship with these climatic indices. This suggests that cosmic-ray muon flux might be linked to the behavior of the climate phenomena and therefore can be used as a climate or meteorological index over the Caribbean.


Introduction: Muon Teleconnection
Outside of the well-known diurnal atmospheric cycles, other oscillation time scales vary from weeks to months, and years. The Madden-Julian Oscillation (MJO), El Nino and Southern Oscillation (ENSO) and Quasi-Biennial Oscillation (QBO) are just a few examples. It has been argued by researchers such as [1] and [2] that weather and climate phenomena (some of which are strongly correlated to each other), are also linked to cosmic-ray muon flux (CMF). The Caribbean based on its geographic, atmospheric and climatic features-some of which can be argued as being fixed or change very little over time-can be considered as being unique. For example, the highest mountains are found in the central Caribbean over Hispaniola and Jamaica [3] and sea surface temperatures above 28.5˚C referred to as the Caribbean Warm Pool constitute part of the second largest body of very warm water on earth [4]. Both of these contribute to cloud development, influence the movement of weather systems and the intensity of rainfall across the central Caribbean (especially). These and other features such as tropopause height and the geomagnetic cut off rigidities, do not have frequent significant change in values, but also have an impact on cosmic ray muon flux and make the region and sub-region unique.
Yet another significant feature in the Caribbean region is the mid-summer drought (MSD) which is a relative lull in rainfall activity during the period late June-early August of each year. In their effort to understand this phenomena [5] suggest a link between the MSD and sub-micron aerosols during this period. [6] also argued that a large aerosol flux of anthropogenic origin occurred during June, July and August. This they argue contains large amounts of submicron aerosol and suggest that this increased amount leads to a reduction in rainfall during the period. [7] explains that atmospheric aerosols affect precipitation by reducing the amount of radiation reaching the ground and therefore suppress precipitation in shallow and short-lived clouds. They also argued that this is due to increased heat added to the clouds in which they mix and leads to reduced cloud coverage as a result of evaporation. [8] argued a connection between galactic cosmic-ray (GCR) and ENSO trends, while [9] highlighted the emergence of Cosmoclimatology and a link between GCR, cloud condensation nuclei (CCN) and cloud formation based on satellite observed clouds. However, this has been very controversial as [10] in their review of this cloud-climate-cosmic ray flux hypothesis, found strong arguments on both sides, proving the need for more affirmative investigations. These arguments are of much interest to Caribbean researchers since CCN liberated by GCR are below normal sizes (<20 nm) and rainfall over the Caribbean over the last 20 to 30 years is generally on the decline.
Muons are sub-atomic particles produced when cosmic radiation enter the earth's atmosphere and interact with nuclei and other aerosols. The initial reaction generates Pions and Kaons which decay immediately to form muons, neutrinos, and electrons. Authors such as [11], [12], [13] and [14] have suggested Atmospheric and Climate Sciences that some teleconnections exist between cosmic-ray muons flux (CMF) and meteorological phenomena based on the physical relationship between atmospheric temperature and pressure. According to [15], over the past 30 years, the NAO has been trending towards a more positive phase, which indicates warmer and wetter conditions over eastern USA. [16] argued that the abnormally high solar activity between solar cycles 23 and 24 has caused the longest solar minima in over 80 years and an unusually small solar maxima of cycle 24. This has resulted in very weak solar activity, the highest flux of galactic cosmic rays and relatively small amounts of energetic particle events. [17] also argues the existence of a link between solar activity and cosmic-ray flux, by stating that a more active Sun has a stronger associated solar wind, which reduces the flux of cosmic rays reaching the Earth. [18] was among the first to show that the variation in cosmic ray muon flux (CMF) is conspicuously related to latitude, altitude and other factors such as the strength of geomagnetic cut off rigidity (GCOR). [19] have also shown that the diurnal anisotropy (unequal physical properties along different axes) of cosmic-ray has both latitudinal and longitudinal dependences and both arguments suggest that CMF can be specifically linked to parameters over a base point region.
Observed trends in the factors such as latitude, GCOR, tropopause height, temperature and pressure which have been shown to influence cosmic ray muon flux decrease polewards. However, air density and the general flux in cosmic ray muons increase polewards. [9] and [20] go further and also suggest that while clouds below 3 km are uncorrelated at mid and high latitudes with cosmic-rays, a strong correlation exists with clouds at low latitudes. This, therefore, gives impetus to this study as we seek a regional or Caribbean CMF perspective. Our hypothesis is therefore that extra-atmospheric factors can result in significant events such as ground level enhancements (GLE) and Forbush decrease (Fds). Apart from these phenomena the contribution to changes in cosmic ray muon flux would not be significant. The dominant factors that will modify the muon flux in the troposphere are therefore a) geospatial-that is the area of the base point region, the GCOR and the location of detector above sea level-and b) the geo-meteorological factors-that is the frequency and type of weather systems, the height of the tropopause and in our case the presence of the Caribbean warm pool (sea surface temperature anomalies). These are the factors that will impact atmospheric temperature and pressure and by extension the muon flux. These factors are also susceptible to being influenced by global and regional teleconnections. Furthermore, if the cosmic-ray muon flux can be linked to large scale circulation and by extension be used to explain the variability distinguished from small-scale localized effects, then as suggested by [21] it could also be used to account for 10% -20% of the total seasonal variance (depending on season considered).
Establishing linkages between atmospheric and climate variables within standardized base points with other elements from non-contiguous regions is described as teleconnection. According to [22], these base points can also be des-

Teleconnection Base Point and Domains
In the broadest sense, Caribbean countries are those which border the Caribbean Sea. The sub-divisions include the Greater Antilles, which include Jamaica in the northwestern section, and the Lesser Antilles which include Barbados located outside the chain of islands in the east. The area designated as "base point" is one established by the Meteorological Service, Jamaica (MSJ) as the "severe weather grid" (hereinafter the SWG). This area is bounded by coordinates 14˚N/ 71˚W, 14˚N/83˚W, 21˚N/83˚W and 21˚N/71˚W with Jamaica almost at its centre [24]. Several important studies of the region have described different geographical areas as the Caribbean. For example [25] Figure 1 is a schema of the region with SWG.
This study therefore, had 3 main components. The first was to establish a base region (the SWG) as described earlier (Figure 1). This specifies a designated reference area for the data gathered and provides a Caribbean context. The   . Of all these indices the link between the arctic oscillation AO and Caribbean weather or climate is least or never been explored. The indices representing these climatic phenomena were sourced from various research data sites associated with the US National Oceanic and Atmospheric Administration (NOAA). The period of the study coincided with the mini-maximum of solar cycle 24 and represents a random sample. To assist with linking the various regions involved Figure 2 shows a map of the globe showing the source region of the various indices. The following sections give background on the climate indices used to establish teleconnection, this is followed by discussion regarding data, method, results, and conclusion based on the statistical correlation.

The Quasi-Biennial Oscillation (QBO)
The QBO is a downward propagating easterly and westerly wind regime that occur within the tropical stratosphere and troposphere. This climate phenomenon have both tropical and extra-tropical effects with a wide range of impacts extending as far as the mesosphere [27]. The effects on the Caribbean region is likely from a combination of the tropospheric and stratospheric QBO influencing the number and intensity of tropical cyclones over the Caribbean. For example, the seasonal number of tropical storms and hurricanes over the Atlantic area, are above normal when the stratospheric QBO is in its westerly phase when compared with its easterly phase [28]. Figure 3 indicates the trend in the stratospheric Quasi-Biennial Oscillation Index (QBOI) during the research period. The QBO index (QBOI) used in this study was derived from the USA Na- Both the easterlies and westerlies transitions show a strong preference to occur during April to June. It is noteworthy that the QBO influences the breakdown of stratospheric polar vortices during the winter and that the polar vortex in the stratosphere, affects surface weather patterns [27]. [29] goes further to argue that when the QBO is easterly (QBOE) or negative phase in the lower stratosphere, stronger Madden Julian Oscillation (MJO) activity is likely during the northern hemisphere winter. This relationship they argue has become more evident since the 1980s when a cooling trend in the equatorial lower stratosphere and a warming trend in the equatorial upper troposphere became more evident. This is likely to play an important role in the production of cosmic ray muons and the flux to the surface since cosmic ray muon production starts in the lower stratosphere.

El Nino Southern Oscillation (ENSO)
The changes in weather to extreme levels over the globe, from droughts to floods, storms to light stagnant winds and from blizzards to heat waves, have  [25] and [26]. The trend in the ONI and SOI over the research period are shown in Figure 4 and Figure 5 below.  Darwin. This is typical of the El Nino episodes over the eastern tropical Pacific. [30] found that ENSO affected hurricane activities by shifting the location of large scale convection. This is due to the upper tropospheric winds over the Caribbean and tropical Atlantic becoming more westerly. When the lower level tropospheric wind speeds are relatively small, this results in enhanced vertical wind shear and reduced hurricane activity during El Nino hurricane season over the Caribbean. The situation however, is reversed during La Nina episodes.

North Atlantic Oscillation (NAO)
The North Atlantic Oscillation (NAO) is a large scale seesawing movement of atmospheric mass between the North Atlantic subtropical high pressure system and the Icelandic low pressure system [31]. [32] also showed that long term upward trend in NAO accounted for much of the surface air temperature warming poleward of 20˚N. The changes in circulation associated with changes in NAO index [32] explained, was based on a difference in sea level pressure between winters with an index value greater than ±1 mb. This scenario was similar to the Southern Oscillation (SO) since both involve pressure changes and are linked to sea surface warming. The NAO is however, recognized as a regional phenomenon that influences weather over the Atlantic, with its depth limited to the troposphere. It has been commonly stated that the NAO index, when in a positive phase, results in mild and wet winter conditions over the eastern USA while colder, snowier weather conditions accompany the negative phase. [31] concluded that the influence of the NAO on the rainfall variability pattern over Barbados was dependent on the phase and intensity of ENSO events. While [26] argued that NAO was negatively associated with Caribbean rainfall when the preceding winter was considered and positively associated with rainfall over the southeastern Caribbean (especially during the late summer) when monthly rainfall was considered. Figure 6 shows the trend in the NAO during the study period.
This investigation therefore, is in part, a continuation of the search for teleconnections. In this case, it is between climate phenomena and an extra-atmospheric one (cosmic rays) and its secondary components (cosmic ray muons) that is influenced by atmospheric conditions. For example, the mechanisms that excite and sustain the North Atlantic Oscillation (NAO)-one of the notable and earliest discovered patterns-are still being investigated." As [22] stated most eloquently regarding the NAO, "although the structure of the prominent teleconnection patterns has been known for some time, the reasons for their origin are not yet well understood.

Arctic Oscillation (AO)
The Arctic Oscillation (AO) is a seesaw in the atmospheric pressure between the Arctic basin and the surrounding zonal ring [33]. [33] also argued that the AO is similar to the NAO with a slightly larger horizontal scale and a higher degree of zonal symmetry. These authors further argued that the NAO is a regional manifestation of the AO, since the AO is hemispheric. The strength of this argument is that the NAO and AO share similar characteristics due to the geographic area covered and the atmospheric depth occupied by both. [34] suggest however, that the AO extends into the stratosphere and that this phenomenon is linked to the sudden warming of the stratosphere (aka sudden stratospheric warming or SSW) and the QBO. Both the QBO and the NAO have been shown to influence Caribbean weather due to stratospheric-tropospheric coupling. There is however, a paucity of work seeking to determine the impact or linkages between the AO and Caribbean weather. This study and others serve to fill that knowledge gap. Figure 7 shows the trend in the AOI over the research period and the web source where the data was acquired.
A positive AOI indicates stronger winds associated with the polar vortex which results in cold air closer to the Polar region while a negative AOI is an indication of weaker westerlies associated with the Polar vortex resulting in cold arctic air moving further south and storms tracts remaining further south. According to [35]   AO occur, while a 0% chance of extreme freeze over is likely when extreme values of positive AO occur. We also suggest that during the episodic freeze overs over Florida that the frontal systems that trigger these events are likely to also affect the central Caribbean and Jamaica.

Summary
Depending on the phase and severity of these indices researchers have shown that the Caribbean can be impacted by large-scale climate phenomena. What is also clear is that a connection has been made between large-scale climate phenomena centred over distant regions and the Caribbean but there has been no attempt that the author recognized after review the available literature to link these with cosmic ray muon flux over the Caribbean. Importantly the production of cosmic ray muon is initiated in the lower stratosphere and the flux of muons is regulated by the depth of the atmosphere whose expansion and contraction is regulated by the temperature in the lower stratosphere and upper troposphere. Several randomly chosen climate phenomena were used along with a randomly chosen period to determine the possible teleconnection. The method, results and conclusions follow.

Data
The data captured was for the period September 2011 to October 2012. The bin width used in this study was the 5-minute bin width and consequently, the data for this entire period was referred to as Muon5. This dataset subsequently di-  Table 1 below indicates the large-scale phenomena and the name of the indices used in this study.
In considering the teleconnection between base points, the geomagnetic cut off rigidity (GCOR) values are also important. This is a quantitative measure of the shielding provided by the Earth's magnetic field [37]. Both [37] and [38] highlighted the importance of this parameter and argued that the GCOR is an ordinate that describes the charge particle access at any location within the magnetosphere. [37] further stated that cut off rigidity helps to predict the ener- Another large-scale climate phenomenon of interest was the sudden stratospheric warming (SSW), which was not considered in this study. However, the SSW does not have an index but have been found to be positively correlated to cosmic-ray muon flux in the lower stratosphere and inversely correlated to Jamaican dry season rainfall in separate studies.

Instrumentation
In the absence of a Caribbean-wide cosmic-ray muon detection network, a single detector was used in this study. The QuarkNet 6000 cosmic-ray muon detector is similar to the one used by [41] and further described in [42]. The various pieces of the detector are shown in Figure 8.
In Figure 8, the actual apparatus used during the investigation is presented.
This apparatus involves the electronic devices and structural material to enable a stable setup of the cosmic ray muon detector as described by [41] and [43]. This system provided a cost-effective method for obtaining high-quality data based on coincidence theory when compared to a more sophisticated and expensive  The miniature computer and laptop were used to remotely access and view the data collected by the detector. This detector is similar to the type used in [41], however, in [41] a method where the orientation of the detector could be varied to determine the direction of muons impacting the detector was also employed. In current study, the detector orientation was fixed. There were no modifications to the primary equipment setup, however, structural material was used only to house and secure the detector parts. Section 3.3 explains the method used to arrive at our conclusion.

Method
As indicated in the previous sections there are many atmospheric and oceanic oscillations that impact meteorological and climatic systems. The advent of the field of muon Diagnostics facilitated the exploration of a link between cosmic-ray muon flux, atmospheric and extra atmospheric phenomena. This study is therefore motivated by: 1) the increasing view in the research community that cosmic-rays play an important role in our changing atmospheric and climate systems, 2) the paucity of information about the behavior of this ubiquitous sub-atomic particle, its intensity, impact and relationship with meteorological and climatic phenomena over the Caribbean, and, 3) the need to pursue additional explanations to those used for Caribbean climate and weather variability. The cosmic-ray muon flux data captured were first divided into several bin widths, ranging from 5-minute to 24-hours. This unfiltered period is being viewed as having as much as 70 percent of the subatomic particles measured at the surface being muons [40]. The CMF represented a period when 100 percent of the data collected were muons. The strength of the resulting correlation was then used to indicate the link or possible teleconnection between the cosmic ray muon flux (CMF+ and CMF) and the large scale climate indices. In current study, we acknowledge that other atmospheric and extra-atmospheric factors are likely to have played a role in the enhancement or suppression of the influence of the large-scale climatic phenomena used and their subsequent the relationship between the monthly Muon5 (Mu5M), [44]. This wider examination was not the focus of this study but will be pursued later. Since rainfall is a fundamental atmospheric variable, the mean monthly rainfall over the island was therefore used to evaluate the strength of the relationship between the climatic indices and cosmic-ray muon flux for the base point.

Connections between Climate Phenomena
The climate Normals for the period Jan 1961-Dec 1990 are considered by the World Meteorological Organization (WMO) as reference normals for climate change analysis [45]. An examination of the correlation statistics before and after Dec 1990 confirmed the strength of the correlations among the climate indices used. Table 3

Teleconnections with Cosmic-Ray Muon Flux
Significant information can be retrieved from underground and underwater muon investigations. For example, a direct correlation has been shown to exist between underground muon intensities and muon production in the stratosphere due to the behaviour of Pions and Kaons to temperature changes, [46]. [47] has shown that the rate of high energy cosmic-ray muons are strongly correlated to upper air temperature during short term atmospheric events know as sudden stratospheric warmings. This leads to increase muon rate when Pion and Kaon decay, rather than interact as it is the decay process that lead to muons being formed. [48] also argues that there is a strong connection between sudden atmosphere. These changes affect, for example, the atmospheric temperature (with increases of up to 50˚C), winds (by changing wind direction and speed) and electron density from the earth surface to the thermosphere. They further argue that there is a connection between space weather and SSWs and that SSWs are a significant source of variability in the earth's thermosphere and ionosphere (especially the high ionospheric conductivity regions of the equatorial low latitudes).
Since the muons are the most abundant charged particle at the earth's surface composing up to 70 percent of cosmic-rays [43], [46] and [49], a strong argument can be made for teleconnection between the various atmospheric phenomena, not only during the stratospheric sudden warming periods but for all seasons. This is so since cosmic ray muon flux is impacted by the changes in temperature and pressure during the different seasons. By traversing the atmosphere at the nanoscale it is a prime candidate to be used for baro-thermal assessment. [22] also argues that although teleconnection patterns often evolve on a sub-monthly timescales, their spatial patterns are typically characterized in monthly and seasonal data. The results are presented in Table 4 and Table 5 below for the respective periods.
In Table 4, the values generated, revealed significant correlation coefficient values (CCV) between all the indices. All the CCVs for the CMF+ were positive except for the ONI which was negatively correlated with all other large-scale climatic indices including the unfiltered muon flux. A similar correlation assessment  Table 5.
The results in Table 5  compared with those during the study period, only very small differences were observed. Table 6 below explains the comparison.
The coefficients in Table 6 shows a similarity between the indices in CMF+

Summary
Oceanic and atmospheric oscillations are reflected in various parameters such as wind, temperature pressure and density over different regions globally. The wind direction as well as the speed is of significant importance especially in the upper troposphere and lower stratosphere. [50] suggest that the stratospheric-  [29] suggest that climate change promotes coupling between large-scale enhanced deep convection and that the emergence of this relationship coincides with the recent cooling trend in the lower equatorial stratosphere and the warming trend in the upper equatorial troposphere. [47] and [51] however argued that an increase in temperature in the stratosphere results in a decrease in air density, which reduces the chances of interaction of Pions and Kaons, which then results in a greater portion decaying to form muons. This, therefore can result in a higher muon flux at the earth's surface. In this work, we also argue that the observed downward trend in naturally occurring unfiltered cosmic-ray muon flux (CMF+) and cosmic ray muon flux (CMF) isolated from the background noise also implies general cooling in the upper troposphere and or lower stratosphere. The observed trend in muon flux over the base point within the Caribbean also coincides with lower solar activity over the past 30 plus years which coincides with a 30 year cooling trend in the lower stratosphere and upper troposphere, a 30 positive trend in the NAO and a general reduction in rainfall over the Caribbean over the past 20 -30 years. This link to the lower stratosphere is very important since this is the area of initial muon production. [29] also argued that a link exists between the NAO, the QBO and large regions of enhanced deep convection in the tropics known as the madden Julian oscillations (MJO). Interannual variation of MJO activity, due to the QBO they argued, is much greater than the variation due to the tropical sea surface temperatures, which is the usual source for interannual variation of tropical convective activity. This tropical connectivity also plays an important role as control mechanisms for the flux of muons to the surface [2].
The changes in phase of climate indices (representing a base point) from positive, neutral or negative, indicate the level of influence on weather, that is, their influence on the parameters such as temperature, pressure, wind and rainfall. More storms are therefore likely over the Caribbean during the positive phase of NAO and less during the negative phase. More storms are also likely during the La Nina phase over the Caribbean and less storms during El Nino phase. More storms are likely when the QBO is in its westerly (QBOW) or positive phase than when the QBO is easterly (QBOE) or its negative phase. The convective activity and presence of storms affect the depth of the atmosphere which also plays an important role in the flux of muons over any region. An assessment of the study period using 4 separate 6-months time spans was carried out in order to capture the most dominant behaviour of each climatic index. Table 8 outlines the dominant impact of the climatic index over the 4 separate 6-month impact periods.
The assessment revealed various changes in phase over the period, however the cosmic ray muon flux being the only one remaining unchanged during period 1 and period 2. This consistent tendency regardless of the change in other  climate indices in the muon flux was referred to as cumulative teleconnection effect or CTE. It should be noted that since no climatic variable acts independently, their impact are therefore based on one or more phenomena influencing the weather condition over any base point region, which then results in either an enhancement or suppression of weather and or the other indices. There is therefore evidence that a teleconnection may exist between the indices from the various source regions, the cosmic-ray muon flux and weather conditions within the severe weather grid (SWG). Section 6 outlines the conclusions drawn from this study.

Conclusions
The strength of teleconnection was based on the statistical correlation between indices representing climatic phenomena for incontiguous regions and a base point was investigated and discussed. The relationship between these large-scale climatic indices and two (2) periods of cosmic-ray muon flux, which was gathered consecutively from a single detector, was used for the investigation. The first period was unfiltered and composed of up to 70% cosmic-ray muon flux (CMF+) while the second consecutive period (the test period) was filtered using coincidence theory and thus contained only the flux of cosmic ray muons (CMF).  and, Arctic and North Atlantic oscillations (AO-NAO). Over the study period the ONI-SOI was uncorrelated during the CMF period however the AOI-NAOI remained strong. Since the connection between the NAO and the Caribbean is well documented, this trend suggests the need to explore a teleconnection between cosmic ray muon flux, weather over the Caribbean and the Arctic oscillation.
3) The correlation coefficient values over the CMF+ period (Sep 2011-Oct 2012) were all positive except for the ONI which was negative while those over the CMF period (Nov 2012-Sep 2013) were all negative except for ONI which was positive. Of those climate indices that were positive during the CMF+ period the SOI was robust but not significant (0.43). This was also the case during the CMF period where the SOI was not significantly correlated (−0.05). We concluded that this reduced correlation coefficient values may be partly due to the reduced intensity in cosmic-ray muon flux generally along with impact of climate phenomena. This therefore meant that the teleconnection between SOI and CMF+ was robust but not significant and that during the CMF period there was no teleconnection. This was likely because the source region was the farthest from the Caribbean. 4) The change in correlation coefficient values, generally from positive over the CMF+ period to negative values for the CMF period, was attributed mainly to increased sensitivity due to increased muon flux component which resulted in a phase shift. This sensitivity could also be partly due to increase instrument sensitivity as a change was made to detect muons only.
Our focus on a base point in the Caribbean region, points not only to the uniqueness of the study, but also to the need for Caribbean researchers to determine alternate mechanisms which may be impacting the region's weather and our changing climate system. Teleconnections between weather and climatic indices are therefore strongly suggested, based on the impact on constituent elements such as temperature, wind, atmospheric density and rainfall but the time is now ripe for enthusiastic assessment of impacts of an extra-atmospheric origin. A strong statistical link was also observed between the cosmic-ray muon flux observed from the data gathered at the UWI Mona with most of the large-scale indices over the study period (short term). [48] also argues that there is a strong connection between sudden stratospheric warmings (SSW) and extensive changes throughout the earth's atmosphere. Since an index is not readily available for the SSW but it has been linked to cosmic ray muon flux (CMF),