1. Introduction
New Caledonia is a group of islands composing one of the largest and most populated territories in the tropical South Pacific. Located between 18°–23°S and 163°–169°E, about 1500 km east of Australia, it is made up of a long and narrow main island (the “Grande Terre” or Mainland) and other smaller islands. The main island is roughly 400 km long, 50 km wide, and features a central mountain range averaging 800 m in altitude, with the two highest peaks reaching around 1600 m. On a climatological average, New Caledonia is swept by the South Pacific southeasterly trade winds that are almost parallel to the 45° orientation of the Grande Terre (Lefèvre et al. 2010), as shown in Fig. 1. Although winter temperatures frequently drop below 10°C, it experiences a tropical climate with an average temperature of 18°C or higher every month of the year. The annual cycle of rainfall can be decomposed into a warm, wet season from November to April, peaking in the austral summer, and a dry, cold season from May to October. Rainfall is a paramount matter in New Caledonia. At the end of the dry season, agricultural activities and drinking water resources are particularly vulnerable to drought hazards on the West Coast plains. Subsequent rainy spells often come as a relief in the fight against wildfires considered as the main cause of degradation of primary forests and their unique ecosystem (Barbero et al. 2011). Reducing New Caledonia’s vulnerability to heavy precipitating events is also a major concern as some inhabited areas are exposed to landslides and flash flooding (Governement de la Nouvelle-Calédonie 2016).
Long-term mean 10-m wind velocity (vectors) and intensity (shadings) around New Caledonia.
Citation: Monthly Weather Review 148, 1; 10.1175/MWR-D-19-0177.1
Heavy rainfall generally occurs on short time scales ranging from a few hours to a few days. Within these time scales, the diurnal cycle is expected to play an important role as it influences the timing of rainfall initiation or intensification. In the tropics, the diurnal cycle of convection and corresponding wind patterns control the timing of diurnal precipitation peaks (Mori et al. 2004). Over the open ocean, a precipitation maximum is generally reported at night or in the early morning (Yang and Smith 2006; Kikuchi and Wang 2008), resulting from convective instability created by nocturnal longwave radiative cooling in the mid- and upper troposphere. Over islands, studies in the Maritime Continent and the western Pacific have shown that provided the land surface is large enough, tropical islands are rainier than the open ocean and exhibit a diurnal cycle with a precipitation peak in the afternoon due to local forcings (Saito et al. 2001; Yang and Smith 2006; Kikuchi and Wang 2008; Vincent and Lane 2016). The main forcing is the sea breeze (Kikuchi and Wang 2008; Cronin et al. 2015). Owing to the smaller heat capacity of land surfaces compared to the ocean, the magnitude of the warming and cooling of the island with the day–night cycle is much larger than that of the ocean. A sea breeze develops in the afternoon when the land is warmer than the ocean and reverses into a land breeze at night when the land cools. Using an idealized model experiment of a flat, circular, low heat-capacity island, Cronin et al. (2015) have proposed a mechanism involving this land/sea-breeze cycle to explain the specific features of rainfall on tropical islands. The convergence at the sea-breeze front is responsible for the formation of low-level convective clouds around noon, before these shallow clouds develop into deeper precipitating convective cells in the afternoon that eventually merge together and converge at the center of the island, triggering heavy rainfall in the late afternoon or early evening. Evidence of this process over real islands in the Maritime Continent is also shown in Qian (2008). Owing to the heavy afternoon rainfall, the island experiences heavier rainfall than the surrounding ocean. However, the size of the island and its orography also influence the rainfall pattern: according to Sobel et al. (2011), the larger the island, the more rainfall and for large islands, mountains increase precipitation. Indeed, in the case of mountainous islands, orographic effects can lead to further enhancement of the afternoon peak in the diurnal cycle. First, a mountain range tends to regularly experience afternoon cumulus formation because of the elevated surface heating, and is therefore a favorable place of its own for a strong diurnal cycle. Second, its interaction with the sea breeze is likely to intensify cumulus formation in the afternoon through orographic lifting (Yang and Chen 2008; Vincent and Lane 2016). As for New Caledonia, Lefèvre et al. (2010) showed that the main island exhibits a strong diurnal cycle in the wind that is related to a diurnal cycle in surface heating. They also highlighted an asymmetric sea-breeze pattern between the windward and leeward sides. Consequently, precipitation on the island is also expected to exhibit an asymmetric pattern in average rainfall, that has already been documented (Météo-France 2008), and a diurnal cycle that will be the purpose of this study.
At longer time scales, the climate in New Caledonia is influenced by large-scale meteorological patterns located both in the tropics and the midlatitudes of the Southern Hemisphere. Midlatitude highs and lows, and especially the Kermadec high, control the orientation of the prevailing trade winds, while tropical influences impact New Caledonia through the seasonal variations in the location of the South Pacific convergence zone (SPCZ, Griffiths et al. 2003). Being located to the south of the mean position of the SPCZ, New Caledonia’s weather is impacted by the frequent southward shifts of the SPCZ in the wet season, which bring increased precipitation (Lefort 2005; Moron et al. 2016). Global modes of variability then modulate these climatological processes. At the interannual time scale, El Niño–Southern Oscillation (ENSO) influences the orientation of the SPCZ (Vincent et al. 2011; Salinger et al. 2014) and the probability of occurrence of tropical cyclones, that is increased in La Niña phase and decreased in El Niño phase (Leroy and Wheeler 2008). At the subseasonal time scale, the Madden–Julian oscillation (MJO; Madden and Julian 1971) impacts wind patterns (Lefort 2005; Lefèvre et al. 2010), while it is also responsible for shifts in the SPCZ (Matthews 2012; Moron et al. 2016) and changes in the probability of observing a tropical cyclone (Leroy and Wheeler 2008). As a result, Leroy (2006) highlights an impact of the MJO on New Caledonia rainfall, and Moron et al. (2016) also show a modulation of rainfall by the subseasonal variability, although they need to use a tailor-made OLR-based intraseasonal index centered on New Caledonia to demonstrate it.
Moreover, the MJO not only impacts the daily rainfall but also the diurnal cycle around islands, as documented by several studies based on observations (Chen and Houze 1997; Sui et al. 1997; Tian et al. 2006; Suzuki 2009; Hidayat and Kizu 2009; Rauniyar and Walsh 2011; Oh et al. 2012; Kanamori et al. 2013; Peatman et al. 2014) and model simulations (Birch et al. 2016; Vincent and Lane 2017). These studies analyze the direct relationship between the MJO phase and the diurnal cycle of precipitation or convection with composite or case study analysis. Most of them focus on the Maritime Continent area. Their conclusions vary depending on mesoscale phenomena such as orographic effects and sea breeze, specific to the location under study. But they all highlight the differentiated impacts of MJO variability on the diurnal cycle between land and ocean. In particular, through the analysis of TRMM data, Peatman et al. (2014) come to the conclusion that, in the Maritime Continent, 80% of the modulation of the precipitation by the MJO over land result from MJO-induced changes of the amplitude of the diurnal cycle.
Finally, synoptic variability also influences daily rainfall and diurnal cycle over islands. Wind regimes enable to classify the synoptic atmospheric circulation around the area of interest into a few specific and recurrent patterns. They can be related to the large-scale, intraseasonal MJO variability by considering how the MJO modulates their probability of occurrence. They can also be related to the island mesoscale effects that are responsible for variations in the diurnal cycle. Wind regimes were used to study the diurnal cycle, for instance by Qian et al. (2013) over Borneo, by Moron et al. (2015) over the whole Maritime Continent, and more recently by Hopuare et al. (2019) over the Pacific island of Tahiti. In New Caledonia, Lefort (2005) and Leroy (2006) have shown that atmospheric conditions could be classified into a small number of regimes. Lefèvre et al. (2010) have also highlighted the link between wind regimes and the phases of the MJO. With a six-state classification performed on daily New Caledonia rainfall, Moron et al. (2016) have even managed to explain the relationships between the local precipitation and the large-scale atmospheric dynamics at various time scales. However, the link between the subseasonal variability and the diurnal cycle of precipitation in New Caledonia has not been investigated so far and the present study intends to fill this gap.
The following sections document the influence of the MJO and wind regimes on the diurnal cycle of rainfall over the main island of New Caledonia during the wet season, and investigate the relationship between wind regimes and MJO. Section 2 introduces the data and methodology used in this article. Section 3 is devoted to the change in diurnal cycle with synoptic and large-scale circulation. It first presents the influence of the MJO on rainfall and its diurnal cycle in New Caledonia, and that of wind regimes on the spatial variations of the diurnal cycle over the island, independently from the MJO’s influence. Section 4 investigates the link between the MJO and the wind regimes and its impacts on the diurnal cycle. Finally, the main conclusions of the study are presented in section 5.
2. Data and methodology
a. Data
As this study is carried out for the wet season, only the months from November to April will be considered in all datasets.
1) Large-scale data
Since the diurnal cycle of precipitation involves interactions between the large-scale atmospheric dynamics and the island’s specific topography, surface wind is presumably the best field to determine wind regimes and investigate the dynamical signal of the MJO as seen by the island. For this purpose, the daily averaged 10-m winds from the ERA-Interim reanalysis (Dee et al. 2011) are used at 1° resolution over a small domain encompassing New Caledonia (25°–17°S, 161°–170°E). The period considered is 1 January 1979–31 December 2017 (39 years). ERA-Interim 10-m winds are also considered at a lower resolution of 2° on a large tropical domain (30°S–0°, 50°–240°E) in order to represent the circulation anomalies associated to the MJO, along with OLR anomalies computed from the daily averaged OLR from NOAA’s polar-orbiting satellites (Liebmann and Smith 1996).
2) New Caledonia rainfall data
A comprehensive study of the diurnal cycle of precipitation requires the use of high-frequency rainfall data, ideally with an hourly or shorter time step. This work relies on the Météo-France hourly rain gauge data for 28 stations scattered all over the main island of New Caledonia. These stations were chosen among the available stations according to three criteria: they provided hourly rain rates, they had a sufficiently long record over the year and they were deemed reliable by the Météo-France climatological services in New Caledonia. These stations are listed in Table A1 and shown on the map in Fig. 2. We deliberately choose to focus this study on the main island owing to the smaller sizes and the very different topographic profiles of the other islands, which make them likely to exhibit a diurnal cycle that would be governed by different mechanisms. The period of available data is station-dependent and some stations present periods of missing values. To include as much data as possible, the period considered is 1 January 1991–31 December 2017 (27 years), 1991 being the year when the first hourly records started.
Average precipitation at the 28 stations in wet season (NDJFMA). Gray lines indicate elevation levels every 200 m from 200 to 1000 m, retrieved from Shuttle Radar Topography Mission (SRTM) data.
Citation: Monthly Weather Review 148, 1; 10.1175/MWR-D-19-0177.1
3) MJO index
Wheeler and Hendon (2004) have defined a set of two indices to represent the propagation of the MJO and classify every day with significant MJO activity into one of the eight phases representing the location of the convective envelope. These two Real-time Multivariate MJO indices (RMM1 and RMM2) are computed operationally at the Australian Bureau of Meteorology (BoM) from NOAA OLR data (as a proxy for cloud amount) and NCEP2 reanalysis zonal winds at lower (850 hPa) and upper (200 hPa) levels of the atmosphere. In this study, the indices’ values were retrieved directly from the BoM website for the period from 1 January 1979 to 31 December 2017. The (RMM1, RMM2) phase space is separated into eight quadrants and the location of a point indicates the MJO phase. In this phase space, a “canonical” MJO event will be represented with a counterclockwise trajectory, the phase for each day going from 1 to 8. Note that when the amplitude of the MJO [i.e., the norm of the (RMM1, RMM2) vector], is less than 1, the MJO is considered to be weak (no MJO event) and the classification of the days into phases is no longer relevant.
b. Wind regime classification
Wind regime classification aims at sorting out the regional wind maps around New Caledonia (i.e., in the selected 25°–17°S, 161°–170°E domain) into a few specific and distinct circulation patterns. For this purpose, we use the same method as Lefèvre et al. (2010), Moron et al. (2015) and Hopuare et al. (2019), which is based on a k-means clustering algorithm. Zonal and meridional components of the daily mean 10-m wind data at the 90 grid points in the domain were used for each day of the wet season. Standardized anomalies were computed at each grid point by subtracting the long-term mean and dividing by the corresponding standard deviation. Principal component analysis was performed to filter out the smallest scales of spatial variability and the first two principal components, accounting for 69% of the total variance, were retained. The 7069 days of the 39 wet seasons were partitioned in this two-dimensional space using the k-means clustering algorithm with Euclidean distance. The optimal number of clusters was chosen using a classifiability index as described by Michelangeli et al. (1995). The k-means clustering algorithm takes the user-provided k number of clusters to choose k initial centroids (the seeds) at random among the points to classify. The clusters are created by gathering each point to its closest cluster, and a new centroid is computed for each cluster. This step is repeated until the algorithm converges (i.e., no more changes in the clusters with the next iteration). The classifiability index is computed for each number of clusters k. For a given k, N different partitions are computed. The similarity between two partitions is taken as the largest of the k correlations (in the initial wind field space) between the pairs of corresponding centroids. The classifiability index is the average of the similarities in N(N − 1)/2 possible pairs of partitions. In this study, the classifiability index was calculated for a number of clusters k from 2 to 10 and N = 50. The closer the classifiability index is to 1, the less the partition depends on the random initial seeds and the more robust it is.
c. Quantifying the difference between two diurnal cycles
The first term on the right-hand side of Eq. (3) represents the contribution of the difference in daily means, and the second term is the contribution of the difference in diurnal variability.
d. Spatial classification of rain gauge stations
A classification of the 28 rain gauge stations was performed in order to distinguish several geographical zones which exhibit significantly different diurnal cycles of rainfall. Similar to any clustering problem, such a classification requires the definition of an accurate metric of the difference between the items being classified (i.e., rain gauge stations). Our main purpose in this study is to illustrate, for each geographical zone, the influence of the wind regime on the diurnal cycle. For each station, we computed the smoothed composite diurnal cycle for each wind regime by averaging hourly rainfall over the days belonging to that regime and filtering out all harmonics with periods below 8 h. Given that four wind regimes (see section 3) were finally retained, each station is characterized by 96 parameters corresponding to the 24 hourly rain rate values of its smoothed composite diurnal cycle in each of the four regimes. The metric used to quantify the difference between two stations was the square-rooted sum of the squared differences in this 96-parameter space. This metric is relevant for this case as it takes into account both the average rainfall at the station and the diurnal variation: if two stations receive similar rainfall on average, they will be close to each other but they will be even closer if, for instance, they also have similar diurnal peak. The clustering algorithm used for this classification is hierarchical clustering. At each step, the two closest clusters according to Ward’s minimum variance method were merged together. The optimal number of clusters could be determined thanks to the dendrogram to ensure that the clusters were well separated from each other.
e. Statistical testing of the influence of MJO on the frequency of wind regimes
Because they are generally assumed to follow a normal distribution with mean 0 and standard deviation 1, the Pearson’s residuals with values less than −2 or greater than +2 are considered to indicate the phases for which the deviation from the expected distribution of the wind regimes under the null hypothesis is significant. These phases are the largest contributors to the rejection of the chi-square null hypothesis. Given their regime-phase symmetry, Pearson’s residuals similarly indicate if a regime is significantly more or less frequent within a phase. The over or underrepresentation at the 95% level of a given regime within a MJO phase was also tested with bootstrapping. This method leads to conclusions that are identical to the Pearson’s residuals approach.
f. Impact of the interaction between the MJO and wind regimes on the diurnal cycle of precipitation
To understand how the interaction between the MJO and wind regimes impacts the diurnal cycle, a method is proposed based on two simple but opposite models.
This corresponds to considering that the characteristics of the regimes change between phases of the MJO, but the frequency of these regimes do not change. Equation (13) also shows that the deviations of the diurnal cycle riϕ can be taken from the long-term mean
We can expect the observed diurnal cycles to follow a pattern somewhat in between the two extreme simple models 1 and 2, and we can actually determine which model is closest to the observations.
This decomposition into 2 × Nw + 1 terms, which are scalar products of anomalous diurnal cycles against the anomalous diurnal cycle of the MJO phase, represents the contributions of all phenomena responsible for the change in the diurnal cycle within the MJO phase. The first term on the right-hand side of Eq. (16) represents the difference to the long-term mean diurnal cycle accounted for by the change in frequency of each of the Nw regimes, following Model 1. The second term on the right-hand side of Eq. (16) terms represents the difference to the long-term mean diurnal cycle accounted for by the change in the diurnal cycle of the regime within the MJO phase relative to the average diurnal cycle of the regime, according to Model 2. Here s is the second-order residual term.
3. Influence of large-scale dynamics on the diurnal cycle of rainfall
a. Climatology of rainfall in New Caledonia
The average daily rainfall received in the wet season by the 28 selected rain gauge stations is represented in Fig. 2. The rainiest stations are all along the east coast and at specific locations inland that correspond to mountain peaks or passes. Stations along the west coast are comparatively drier. This east–west contrast is usually explained by two factors. First, with the presence of the central mountain range and the east-southeasterly climatological direction of the trade winds (illustrated in Fig. 1), orographic lifting on the windward side combined with foehn effect on the lee side drive a wetter climate along the east coast and a drier climate along the west coast (Météo-France 2008). Second, rainfall in New Caledonia can also be triggered by the advection of warm and moist tropical air masses by northerly and easterly winds that therefore preferentially hit the east coast (Moron et al. 2016).
The diurnal cycle of observed rainfall is represented in Fig. 3 for both the warm, wet (November–April) and the cool, dry (May–October) seasons, averaged over all stations shown in Fig. 2. While there is a strong diurnal cycle in the warm season, with a marked afternoon peak at 1600 local time and a secondary peak at night, rainfall in the dry season shows very little diurnal variation. As detailed in section 1, sea breeze and orographic effects are likely to explain the afternoon maximum in warm convection, and therefore in precipitation, over land (Cronin et al. 2015). Nighttime convection over the open ocean and advection of precipitating clouds by large-scale winds can presumably account for the secondary, nocturnal maximum. During the dry season, the sea surface temperature is lower and the troposphere is much drier than in the warm season: these inhibit deep convection, resulting in lower precipitation rates. The diurnal solar forcing is also smaller in that season, so it creates weaker land/sea surface temperature gradient and sea breeze, which is rarely strong enough to break through the convective inhibition. The diurnal cycle is therefore considered to be one major player to explain differences in rainfall rates between wet and dry season.
Mean diurnal cycle of rainfall (mm h−1) averaged over the 28 stations in the wet season (blue) and the dry season (orange).
Citation: Monthly Weather Review 148, 1; 10.1175/MWR-D-19-0177.1
b. Impact of the MJO on the diurnal cycle of rainfall in New Caledonia
The influence of the MJO on the mean precipitation in New Caledonia (averaged over the 28 stations) is displayed in Fig. 4a, along with the 95% confidence interval determined from a bootstrap procedure (for each phase, 10 000 draws of the days in the phase are performed with replacement). Daily mean rainfall variation oscillates with the MJO, following the propagation of the MJO convective envelope. Phases 4, 5, and 6 are the rainiest phases, and in particular, they are significantly rainier at the 95% confidence level than phases 1 and 2 which are the driest phases. This is in agreement with the MJO-induced OLR anomalies over New Caledonia as shown in Fig. 5: the largest negative OLR anomalies over New Caledonia, associated with intensified convection, are found in phases 4, 5, and 6, while the largest positive OLR anomalies, associated with suppressed convection, correspond to phases 1 and 2. Figure 5 also shows that New Caledonia is located at the southern edge of the strong, equatorial MJO signal and it experiences convection anomalies in advance compared to the equatorial band. Positive or negative MJO-induced convection anomalies extend southeastward from the heart of the equatorial convective patch. For instance, New Caledonia experiences intensified convection when the equatorial convective envelope is located over the Maritime Continent (phase 5), while it experiences reduced convection when the equatorial convective envelope is roughly north of New Caledonia (phase 7).
(a) Daily mean rainfall (mm day−1) averaged on the 28 stations for each MJO phase. (b) Amplitude (mm h−1) of the spatially averaged composite diurnal cycle for each MJO phase. Error bars indicate the 95% confidence interval determined from a bootstrap procedure.
Citation: Monthly Weather Review 148, 1; 10.1175/MWR-D-19-0177.1
NOAA OLR anomalies [OLRa as defined in Wheeler and Hendon (2004)] (shadings, in W m−2) and 10-m wind anomalies (vectors, m s−1) in ERA-Interim for each MJO phase. The black rectangle indicates New Caledonia’s region shown in Fig. 1.
Citation: Monthly Weather Review 148, 1; 10.1175/MWR-D-19-0177.1
The influence of the MJO on the diurnal cycle is summarized in Fig. 4b which shows the evolution of the diurnal amplitude with the MJO phase, along with its 95% bootstrap confidence interval. The diurnal amplitude is defined as the difference between the maximum and the minimum of the composite diurnal cycle of each MJO phase. The impact of the MJO on the diurnal amplitude is less clear than that on average precipitation considering the error bars are larger, but it also exhibits noteworthy features. The diurnal amplitude is largest for phase 4–6, which means that the diurnal cycle tends to peak approximately one phase in advance compared to the mean rainfall (phase 4 versus phase 5). This is reminiscent of Peatman et al. (2014)’s results, which documented a peak in diurnal amplitude with a one-phase lead before the arrival of the MJO envelope over the Maritime Continent. As for the time of the daily rainfall maximum, no significant modulation by the phase of the MJO has been observed (not shown). Therefore, the change in the timing of the afternoon peak in relation with the MJO will not be further investigated in this study.
The squared distance between the long-term mean diurnal cycle
Squared distance between the composite diurnal cycle for each MJO phase rϕ (weak MJO days not included) and the mean diurnal cycle
Citation: Monthly Weather Review 148, 1; 10.1175/MWR-D-19-0177.1
Following Eq. (3), the modifications of the diurnal cycle in Fig. 6 are also split into a contribution from the change in daily mean (in dark gray) and a contribution from the change in the hourly anomalies relative to the daily mean (in light gray). The modification of the diurnal cycle by the MJO appears to be primarily driven by the modification of the daily mean, while the modification of the diurnal variability only plays a secondary role. This is consistent with the results in Fig. 4 in which the impact of the MJO is clearer on the daily mean than on the diurnal amplitude. However, this decomposition neglects the influence of the diurnal cycle on the mean and the contribution from the diurnal anomalies should be considered a lower bound of the contribution from the diurnal cycle. The contribution of the diurnal variability is noteworthy for the five MJO phases when the deviation from the long-term mean diurnal cycle is particularly significant (e.g., in phase 4).
c. Identification of wind regimes
The surface wind anomalies associated with the MJO propagation (Fig. 5) suggest an influence of the MJO on the circulation around New Caledonia. This is an incentive to investigate how the effect of the MJO on the diurnal cycle of rainfall over New Caledonia is mediated through the synoptic wind regimes, as in Leroy (2006), Lefèvre et al. (2010) and Moron et al. (2015). The wind regimes were determined using ERA-Interim 10-m winds as explained in section 2b. The computation of the classifiability index for 2 to 10 clusters is illustrated in Fig. 7. Given the excellent classifiability for four clusters which is also a reasonable number of wind regimes, four wind regimes were retained. This number is consistent with the previous findings of Leroy (2006) and Lefèvre et al. (2010), although the classifications were carried out on different datasets.
Classifiability index for each number of clusters, determined from the two first PCs of ERA-Interim 10-m winds in the New Caledonia domain (25°–17°S, 161°–170°E).
Citation: Monthly Weather Review 148, 1; 10.1175/MWR-D-19-0177.1
The composite wind maps for each regime are represented in Fig. 8. These regimes can be described as follows, in order of increasing average wind speed:
Northerly regime: the flow crossing the island is slowed down on the windward side, while the wind flowing to the west remains slightly stronger. It is the least frequent regime (14.3% of wet-season days) and has the weakest wind speed.
Southerly regime: the flow comes south-southeasterly and is stronger south of the island than north of it. It has a frequency of 19.3% of wet-season days.
Easterly regime: the flow is straight from the east and slows down slightly in the island’s wake. It is one of the two most common regime, with a frequency of 33.2%.
Trade wind regime: the flow is very regular, almost along the main axis of the island. This flow is very similar to, although faster than, the seasonal mean. It is the other most common regime (33.2% of wet-season days) and it has the strongest wind speed.
Composite maps of ERA-Interim 10-m winds (arrows) and wind speed (shadings) for the four wet season regimes.
Citation: Monthly Weather Review 148, 1; 10.1175/MWR-D-19-0177.1
Since wind regimes are supposed to represent recurrent patterns that may last for a few days, the probabilities of transition between regimes are also important additional information about their behavior. Table 1 displays these probabilities. The only transition that is significantly more probable than a random draw at the 95% confidence level is the transition from the northerly to the southerly regime. This corresponds to the eastward propagation of a cyclonic perturbation over the island. All other probabilities that are significantly more probable than a random draw are the probabilities of persistence of each regime. All other transitions are less likely than a random draw.
(top) Probabilities of transition between wind regimes. Row: regime in day d. Column: regime in day d + 1. The probabilities are expressed relative to the regime in day d (rows add up to 1). Bold: transition more probable at the 95% level than in 10 000 random sequences. Italics: transition less probable at the 95% level than in 10 000 random sequences. (bottom) Number and average duration for each regime spell. All statistics are computed for the 1979–2017 NDJFMA period.
d. Influence of the wind regimes on rainfall and its diurnal cycle
1) Spatially averaged influence
Figure 9a displays the diurnal cycle for each regime, averaged over all 28 stations. In terms of average rainfall, the northerly regime is by far the wettest, the easterly and southerly regimes exhibit a similar, intermediate behavior, and the trade wind regime is the driest. The northerly regime can be interpreted as a southward shift of the SPCZ identified in previous studies (Lefort 2005; Leroy 2006). It corresponds to the arrival of warm, moist air from the equatorial band resulting in precipitation over the island. The difference in rainfall averages between wind regimes generally coincides with a difference in amplitude of the diurnal cycle. The afternoon peak in the northerly regime is larger in amplitude than that of the easterly and southerly regime, while the trade wind regime exhibits an even flatter cycle. Another interesting feature for the northerly regime is the secondary maximum during the night. This secondary peak is presumably due to the advection of precipitating cloud from the equatorial oceanic regions. Similar to the MJO phases (Fig. 6), the distance between the average diurnal cycle for each regime and the long-term mean diurnal cycle is illustrated in Fig. 9b. The two regimes for which there is a significant difference are the northerly and the trade regimes, that is to say the more extreme (wettest and driest) ones, in agreement with Fig. 9a.
(a) Spatially averaged diurnal cycle of rainfall for each wind regime (mm h−1). (b) Squared distance between the composite diurnal cycle for each wind regime ri and the mean diurnal cycle
Citation: Monthly Weather Review 148, 1; 10.1175/MWR-D-19-0177.1
2) Spatial classification of rain gauge stations
The presence of a mountain range and the asymmetry in the island’s mean rainfall pattern (see Fig. 2) suggests that the influence of the wind regimes on the diurnal cycle might not be the same for all stations. This matter is investigated by dividing the island into distinct geographical zones following the procedure of spatial clustering detailed in section 2d. The resulting dendrogram is represented in Fig. 10a. It shows that a classification into 3 geographical clusters is relevant in order to obtain distinct, well-separated clusters. The location of the stations in the 3 clusters is represented in Fig. 10b. The clustering algorithm separates the stations into geographically consistent zones: the mountain range (“Range”), the east coast (“East”) and the west coast (“West”). The mountain stations of the Range can be identified by their altitude since they are higher than any other station in East or West. The smoothed, average diurnal cycles for each cluster are represented in Fig. 11. The mountain range and the east coast are both wetter than the west coast, as could be expected from Fig. 1. In terms of diurnal variability, the west coast and the mountain range exhibit a similar behavior with an afternoon peak as in the all-station average diurnal cycle (Fig. 3). On the other hand, the east coast features a flatter diurnal cycle with a minimum in the morning around 1000–1100 and slightly more rain from the midafternoon until the end of the night.
(a) Dendrogram of the hierarchical clustering of the 28 rain gauge stations. (b) Location of the 3 clusters retained with elevation levels indicated by gray lines.
Citation: Monthly Weather Review 148, 1; 10.1175/MWR-D-19-0177.1
Diurnal cycles of rainfall (mm h−1) for each geographical zone.
Citation: Monthly Weather Review 148, 1; 10.1175/MWR-D-19-0177.1
Two hypotheses can be put forward to explain the flatter diurnal cycle in the East. First, the east coast is on the windward side most of the time so it presumably undergoes precipitation that is mostly related to the large-scale wind dynamics (cloud formation due to orographic lifting, cloud advection from the open sea). This kind of precipitation is therefore independent from the local diurnal variability. Second, the local diurnal variability of winds and solar radiation is itself very weak on the east coast. This aspect has already been thoroughly investigated with numerical modeling by Lefèvre et al. (2010). Because of the higher mean cloud cover on the east coast, less solar radiation reaches the surface. The sea breeze then tends to be far weaker than that on the west coast due to this difference in insolation, but also to a difference in sea surface temperatures. Indeed, the sea surface temperature is also warmer off the east coast than off the west coast because New Caledonia is located at the edge of the western Pacific warm pool. This creates a smaller land/sea surface temperature contrast and a weaker sea breeze on the east coast. As a result, rainfall on the east coast shows less diurnal variations, while the west coast experiences a pronounced diurnal cycle of rainfall owing to its strong sea breeze, even if its average precipitation is lower than in the other two parts of the island. As for the mountain range, the diurnal cycle of thermal convection and the interaction with the sea breeze are likely to explain most of its strong afternoon peak.
3) Geographical influence of the wind regimes on the diurnal cycle
Figure 12 shows the influence of the wind regimes for each geographical cluster of stations. The influence on rainfall in clusters Range and West (Figs. 12a,c) is very similar to the influence that is observed in Fig. 9 for the diurnal cycle averaged over all 28 stations. In other words, the rainiest regime is the northerly regime, followed by the southerly and easterly regimes, and finally by the driest trade-wind regime. The diurnal amplitude is also similarly related to the average rainfall: the rainier the regime, the larger the diurnal amplitude. However, the influence of the wind regimes on the east coast (Fig. 12b) shows an altogether different behavior. While the easterly and trade wind regimes are associated with a relatively flat diurnal cycle with a minimum in late morning, the northerly and southerly regimes differ significantly from the average east coast diurnal cycle (Fig. 11). The southerly regime, that is quite dry on the east coast, exhibits a diurnal variability with a marked afternoon peak. The northerly regime is also characterized by a similar afternoon peak and by another secondary peak at the end of the night that has already been identified in Fig. 9.
Spatially averaged diurnal cycle of rainfall for each wind regime (mm h−1): (a) mountain range, (b) east coast, and (c) west coast.
Citation: Monthly Weather Review 148, 1; 10.1175/MWR-D-19-0177.1
We can propose a few tentative explanations for the specific behavior of the diurnal cycle on the east coast in the southerly and northerly regimes, although the mechanisms that explain this sensitivity would deserve a thorough investigation that is beyond the scope of this article. In section 3d(2), two hypotheses have been put forward to explain the flatter average diurnal cycle on the east coast: the continuous rainfall due to large-scale dynamics, and the weaker insolation and sea-breeze diurnal variability. These hypotheses are strongly linked to the fact the east coast is the average windward coast. In the southerly regime, however, the east coast becomes the leeward coast. Sheltered by the mountain range, a sea-breeze circulation can develop and the absence of clouds favors a diurnal cycle of surface heating and afternoon convection. As a result, the diurnal cycle of rainfall on the east coast in the southerly regime exhibits an afternoon peak similar to the diurnal cycle on the west coast in the easterly regime. During the northerly regime, very weak winds and high humidity generally enhance the diurnal cycle of convection all over the island and the east coast is also affected, hence the afternoon peak it also exhibits in this case. The secondary nocturnal peak, on the other hand, can be attributed to the advection of precipitating clouds from oceanic regions further north.
4. MJO impact on the diurnal cycle via wind regimes
The findings in section 3 illustrate separately the impact of the MJO and wind regimes on the diurnal cycle. These two points can be linked, first by considering the influence of the MJO on the wind regimes, and next by considering whether the MJO influence on the diurnal cycle of rainfall is mediated by the wind regimes.
a. Frequency variability of wind regimes with the MJO
We first investigate whether the MJO influences the frequency of occurrence of the wind regimes, in line with Model 1 introduced in section 2f. Following the methodology described in section 2e, we perform a statistical analysis of these frequencies of occurrence in the eight MJO phases. Table 2 shows the p values of the chi-square test applied separately for each regime. The p value is always less than 5%, which enables to reject the null hypothesis that the regime frequencies do not depend on the MJO phase. Table 3 indicates the percentage of the regimes per MJO phase. Percentages corresponding to a significantly more frequent regime within the MJO phase (according to the values of Pearson’s residuals) are in bold, those corresponding to a significantly less frequent regime are in italics.
The p value of the χ2 test assessing the distribution of each wind regime among the MJO phases (weak MJO included).
Percentage of wind regimes by MJO phase. Bold: the regime is significantly more frequent in the phase (Pearson’s residual greater than +2). Italics: the regime is significantly less frequent in the MJO phase (Pearson’s residual less than −2).
A striking fact from Table 3 is that the easterly, northerly and southerly regimes seem to succeed each other in the MJO cycle in terms of overrepresentation. Indeed, the easterly regime occur significantly more frequently in phases 2, 3, and 4, when the MJO convective envelope is still in the Indian Ocean or over the Maritime Continent. On the contrary, it is significantly less frequent in phases 6, 7, and 8 when the convective envelope is over the Central Pacific or weak and further east. The northerly regime occurs significantly less frequently in the phases 1, 2, and 3 when convection is suppressed over New Caledonia, and it occurs significantly more frequently in the wettest phases 5 and 6. As for the southerly regime, it is also significantly less frequent in the dry phases 2 and 3 and significantly more frequent in phase 6. Trade winds are frequent in all phases, except in the phases 5 and 6 of intensified convection over New Caledonia, when they are significantly less frequent. These findings are in agreement with the wind anomalies associated to the MJO phases as represented in Fig. 5. In phases 2 and 3, New Caledonia experiences easterly MJO low-level wind anomalies. These MJO anomalies rotate counterclockwise around New Caledonia with the succession of phases: they have a northeasterly orientation in phase 4, a northerly orientation in phase 5, a westerly orientation in phase 6 and a southerly orientation in phase 7. In phases 8 and 1, which are phases for which the MJO is less active over New Caledonia, the anomalies are southeasterly, with the same orientation as the climatological trade winds. Considering these observations and the transition probabilities in Table 1, it can be deduced that there is a sequence of wind regimes (easterly → northerly → southerly) that is due to the propagation of the MJO wind anomalies. This sequence will not necessarily be observed exactly for every MJO event, because fast intraseasonal and synoptic variability allow all regimes to occur in all phases, but it corresponds to a canonical evolution that the atmosphere tends to follow during such events. It is typical of the MJO dynamical signature (Rui and Wang 1990; Zhang 2005), with easterly winds east of the convective signal and two off-equatorial cyclonic gyres propagating along with the equatorial convective center. The southern gyre is responsible for the wind anomalies around New Caledonia, with a northerly component when New Caledonia is in the poleward branch of the gyre and a southerly component when it is in the equatorward branch.
b. How does the interaction between MJO and wind regimes impact the diurnal cycle?
As seen in Figs. 4 and 6, the MJO modifies the diurnal cycle of precipitation over New Caledonia during its most active and suppressed phases. As for wind regimes, Figs. 8 and 11 also demonstrate their impact on the diurnal cycle. How do these two large-scale influences interact? And how does this interaction help understand the impact of the MJO on the diurnal cycle? The decomposition of the anomalous diurnal cycles proposed in section 2f with two simple models, for the particular case where the number of wind regimes is Nw = 4, is used to answer these questions.
The contributions of the 9 terms in the squared distance to the long-term mean diurnal cycle, detailed in Eq. (16), are represented in Fig. 13. According to Fig. 6, the deviations from the long-term mean diurnal cycle are weak for phases 3, 7, and 8. Consequently, Fig. 13 only represents the decomposition for the five other MJO phases (1, 2, 4, 5, and 6). The contributions are shown as fractions of
Decomposition of the squared distances between the MJO composite cycles rϕ and the long-term mean diurnal cycle
Citation: Monthly Weather Review 148, 1; 10.1175/MWR-D-19-0177.1
The contributions of frequency changes
Phases 1 and 2 (Fig. 13, first and second panels) exhibit a similar behavior: the dominant term in the modification of the diurnal cycle is the change in diurnal cycle for the northerly regime
5. Summary and discussion
This study demonstrates that the Madden–Julian oscillation is associated with rainfall anomalies in New Caledonia that can be noted mostly in the daily average but also, to a lesser extent, in the diurnal variability. These two components are in fact closely entangled together when it comes to modulation by the MJO because the variations in the diurnal amplitude result in variations in the daily mean. MJO phases 1 and 2 are responsible for a weaker diurnal cycle while it is enhanced in phases 4, 5, and 6. The link between the large-scale intraseasonal variability and the local-scale diurnal cycle of rainfall is established through the use of wind regimes as an intermediate tool. Four regimes were identified to classify the surface wind circulation in a small domain encompassing New Caledonia. They are also associated with differences in the rainfall average and in the diurnal cycle, depending on where the air masses come from and how they interact with the island’s topography. Owing to this complex topography, the island can also be objectively separated into three distinct geographical zones that exhibit differences in rainfall amounts and diurnal cycles: the mountain range, the east coast and the west coast. The diurnal cycle in these zones might be subject to quantitative but also qualitative changes depending on the regime, the most notable ones being for the east coast during the northerly and southerly regimes.
From a statistical point of view, the MJO appears to modulate the probability of occurrence of the wind regimes. Its propagation brings about a counterclockwise rotation of surface wind anomalies around New Caledonia that is suggestive of an easterly–northerly–southerly sequence. This sequence is also discernible in the frequencies of occurrence of the wind regimes within the MJO phases, some regimes being significantly more frequent in some phases while the others are significantly less frequent. However, this undeniable relationship does not necessarily mean that the influence of the MJO on the diurnal cycle of precipitation is entirely mediated by the change in frequency of wind regimes. Indeed, Moron et al. (2016, see Fig. 12) also highlight the limitations of the frequentist interpretation for the influence of subseasonal variability when it comes to rainfall itself, except for the wettest days.
This is why this study proposes a framework to assess the different contributions to the sensitivity of the diurnal cycle to the MJO. This framework consists in a decomposition into 1) terms associated with the modification of the frequency of occurrence of the wind regimes by the MJO, and 2) terms due to changes of the characteristics of the wind regimes within the phases of the MJO. It shows that the modification of the regimes’ diurnal cycle of precipitation due to the MJO (Model 2) eventually proves to account for the largest part of the MJO diurnal cycle anomalies, and that the modification of regime frequencies (Model 1) plays a lesser role. This raises the question of how the MJO brings changes in the wind regimes’ diurnal cycles. We can provide several hypotheses which would require further examination. First, the MJO-induced large-scale convective conditions (enhanced, neutral or suppressed) might have different effects on the different regimes but are not directly related to the wind patterns. Second, the regime circulation anomalies in a given MJO phase are presumably involved through modifications of the wind orientation and/or speed.
The application of this framework for New Caledonia therefore highlights complex interactions that go far beyond the simple model assuming that the MJO mostly modulates the frequencies of wind regimes, contrary to what was previously expected considering past work on other regions (Ichikawa and Yasunari 2008; Hopuare et al. 2019). The complexity of these interactions even suggests that other tools besides wind regimes should be considered in order to fully decompose and understand the MJO signal that impacts New Caledonia at a local scale.
Acknowledgments
We acknowledge support from the Pacific Fund of the Ministère des Affaires Étrangères (France) through the research program PLUVAR. All data analysis in this study was performed using the open source language R. The clustering algorithms were carried out with the functions “kmeans” and “hclust” from the R base package “stats.” MJO RMM values were obtained from the Australian Bureau of Meteorology: http://www.bom.gov.au/climate/mjo/graphics/rmm.74toRealtime.txt. We would also like to thank the two anonymous reviewers whose comments helped improve this study.
APPENDIX
Rain Gauge Station Information
Table A1 lists the main characteristics (location, altitude, and recording period) of the rain gauge stations where precipitation data used in this study were measured.
List of Météo-France rain gauge stations and coordinates.
REFERENCES
Barbero, R., V. Moron, M. Mangeas, M. Despinoy, and C. Hély, 2011: Relationships between MODIS and ATSR fires and atmospheric variability in New Caledonia (SW Pacific). J. Geophys. Res., 116, D21110, https://doi.org/10.1029/2011JD015915.
Birch, C. E., S. Webster, S. C. Peatman, D. J. Parker, A. J. Matthews, Y. Li, and M. E. E. Hassim, 2016: Scale interactions between the MJO and the western Maritime Continent. J. Climate, 29, 2471–2492, https://doi.org/10.1175/JCLI-D-15-0557.1.
Chen, S. S., and R. A. Houze, 1997: Diurnal variation and life-cycle of deep convective systems over the tropical Pacific warm pool. Quart. J. Roy. Meteor. Soc., 123, 357–388, https://doi.org/10.1002/qj.49712353806.
Cronin, T. W., K. A. Emanuel, and P. Molnar, 2015: Island precipitation enhancement and the diurnal cycle in radiative-convective equilibrium: Island rainfall enhancement in radiative-convective equilibrium. Quart. J. Roy. Meteor. Soc., 141, 1017–1034, https://doi.org/10.1002/qj.2443.
Government de la Nouvelle-Calédonie, G., 2016: Dossier sur les risques majeurs de la Nouvelle-Calédonie. Direction de la Sécurité Civile et de la Gestion des Risques de la Nouvelle-Calédonie, accessed 20 November 2019, https://securite-civile.nc/documents/dossier-sur-les-risques-majeurs-de-la-nouvelle-caledonie.
Dee, D. P., and Coauthors, 2011: The ERA-Interim reanalysis: Configuration and performance of the data assimilation system. Quart. J. Roy. Meteor. Soc., 137, 553–597, https://doi.org/10.1002/qj.828.
Griffiths, G. M., M. J. Salinger, and I. Leleu, 2003: Trends in extreme daily rainfall across the South Pacific and relationship to the South Pacific Convergence Zone. Int. J. Climatol., 23, 847–869, https://doi.org/10.1002/joc.923.
Hidayat, R., and S. Kizu, 2009: Influence of the Madden-Julian Oscillation on Indonesian rainfall variability in austral summer. Int. J. Climatol., 30, 1816–1825, https://doi.org/10.1002/joc.2005.
Hopuare, M., M. Guglielmino, and P. Ortega, 2019: Interactions between intraseasonal and diurnal variability of precipitation in the South Central Pacific: The case of a small high island, Tahiti, French Polynesia. Int. J. Climatol., 39, 670–686, https://doi.org/10.1002/joc.5834.
Ichikawa, H., and T. Yasunari, 2008: Intraseasonal variability in diurnal rainfall over New Guinea and the surrounding oceans during austral summer. J. Climate, 21, 2852–2868, https://doi.org/10.1175/2007JCLI1784.1.
Kanamori, H., T. Yasunari, and K. Kuraji, 2013: Modulation of the diurnal cycle of rainfall associated with the MJO observed by a dense hourly rain gauge network at Sarawak, Borneo. J. Climate, 26, 4858–4875, https://doi.org/10.1175/JCLI-D-12-00158.1.
Kikuchi, K., and B. Wang, 2008: Diurnal precipitation regimes in the global Tropics. J. Climate, 21, 2680–2696, https://doi.org/10.1175/2007JCLI2051.1.
Lefèvre, J., P. Marchesiello, N. C. Jourdain, C. Menkes, and A. Leroy, 2010: Weather regimes and orographic circulation around New Caledonia. Mar. Pollut. Bull., 61, 413–431, https://doi.org/10.1016/j.marpolbul.2010.06.012.
Lefort, T., 2005: Starting up medium-range forecasting for New Caledonia in the South-West Pacific Ocean—A not so boring tropical climate. ECMWF Newsletter, No. 102, ECMWF, Reading, United Kingdom, 2–7.
Leroy, A., 2006: Utilisation des Prévisions Saisonnières en Nouvelle-Calédonie. Note de la DP 6, Météo-France, 171 pp.
Leroy, A., and M. C. Wheeler, 2008: Statistical prediction of weekly tropical cyclone activity in the Southern Hemisphere. Mon. Wea. Rev., 136, 3637–3654, https://doi.org/10.1175/2008MWR2426.1.
Liebmann, B., and C. A. Smith, 1996: Description of a complete (interpolated) outgoing longwave radiation dataset. Bull. Amer. Meteor. Soc., 77, 1275–1277.
Madden, R. A., and P. R. Julian, 1971: Detection of a 40-50 day oscillation in the zonal wind in the tropical Pacific. J. Atmos. Sci., 28, 702–708, https://doi.org/10.1175/1520-0469(1971)028<0702:DOADOI>2.0.CO;2.
Matthews, A. J., 2012: A multiscale framework for the origin and variability of the South Pacific Convergence Zone. Quart. J. Roy. Meteor. Soc., 138, 1165–1178, https://doi.org/10.1002/qj.1870.
Météo-France, 2008: Atlas Climatique de la Nouvelle-Calédonie. Accessed 10 December 2018, http://www.meteo.nc/en-savoir-plus/comprendre-la-meteo/publications/atlas-climatique-de-la-nouvelle-caledonie-sommaire.
Michelangeli, P.-A., R. Vautard, and B. Legras, 1995: Weather regimes: Recurrence and quasi stationarity. J. Atmos. Sci., 52, 1237–1256, https://doi.org/10.1175/1520-0469(1995)052<1237:WRRAQS>2.0.CO;2.
Mori, S., and Coauthors, 2004: Diurnal land–sea rainfall peak migration over Sumatera Island, Indonesian Maritime Continent, observed by TRMM satellite and intensive rawinsonde soundings. Mon. Wea. Rev., 132, 2021–2039, https://doi.org/10.1175/1520-0493(2004)132<2021:DLRPMO>2.0.CO;2.
Moron, V., A. W. Robertson, J.-H. Qian, and M. Ghil, 2015: Weather types across the Maritime Continent: From the diurnal cycle to interannual variations. Front. Environ. Sci., 2, 65, https://doi.org/10.3389/fenvs.2014.00065.
Moron, V., R. Barbero, and A. W. Robertson, 2016: Subseasonal-to-interannual variability of rainfall over New Caledonia (SW Pacific). Climate Dyn., 46, 2449–2468, https://doi.org/10.1007/s00382-015-2712-0.
Oh, J.-H., K.-Y. Kim, and G.-H. Lim, 2012: Impact of MJO on the diurnal cycle of rainfall over the western Maritime Continent in the austral summer. Climate Dyn., 38, 1167–1180, https://doi.org/10.1007/s00382-011-1237-4.
Peatman, S. C., A. J. Matthews, and D. P. Stevens, 2014: Propagation of the Madden-Julian Oscillation through the Maritime Continent and scale interaction with the diurnal cycle of precipitation: MJO propagation and scale interaction with the diurnal cycle. Quart. J. Roy. Meteor. Soc., 140, 814–825, https://doi.org/10.1002/qj.2161.
Qian, J.-H., 2008: Why precipitation is mostly concentrated over islands in the Maritime Continent. J. Atmos. Sci., 65, 1428–1441, https://doi.org/10.1175/2007JAS2422.1.
Qian, J.-H., A. W. Robertson, and V. Moron, 2013: Diurnal cycle in different weather regimes and rainfall variability over Borneo associated with ENSO. J. Climate, 26, 1772–1790, https://doi.org/10.1175/JCLI-D-12-00178.1.
Rauniyar, S. P., and K. J. E. Walsh, 2011: Scale interaction of the diurnal cycle of rainfall over the Maritime Continent and Australia: Influence of the MJO. J. Climate, 24, 325–348, https://doi.org/10.1175/2010JCLI3673.1.
Rui, H., and B. Wang, 1990: Development characteristics and dynamic structure of tropical intraseasonal convection anomalies. J. Atmos. Sci., 47, 357–379, https://doi.org/10.1175/1520-0469(1990)047<0357:DCADSO>2.0.CO;2.
Saito, K., T. Keenan, G. Holland, and K. Puri, 2001: Numerical simulation of the diurnal evolution of tropical island convection over the Maritime Continent. Mon. Wea. Rev., 129, 378–400, https://doi.org/10.1175/1520-0493(2001)129<0378:NSOTDE>2.0.CO;2.
Salinger, M. J., S. McGree, F. Beucher, S. B. Power, and F. Delage, 2014: A new index for variations in the position of the South Pacific convergence zone 1910/11-2011/2012. Climate Dyn., 43, 881–892, https://doi.org/10.1007/s00382-013-2035-y.
Sobel, A. H., C. D. Burleyson, and S. E. Yuter, 2011: Rain on small tropical islands. J. Geophys. Res., 116, D08102, https://doi.org/10.1029/2010JD014695.
Sui, C.-H., X. Li, K.-M. Lau, and D. Adamec, 1997: Multiscale air–sea interactions during TOGA COARE. Mon. Wea. Rev., 125, 448–462, https://doi.org/10.1175/1520-0493(1997)125<0448:MASIDT>2.0.CO;2.
Suzuki, T., 2009: Diurnal cycle of deep convection in super clusters embedded in the Madden-Julian Oscillation. J. Geophys. Res., 114, D22102, https://doi.org/10.1029/2008JD011303.
Tian, B., D. E. Waliser, and E. J. Fetzer, 2006: Modulation of the diurnal cycle of tropical deep convective clouds by the MJO. Geophys. Res. Lett., 33, L20704, https://doi.org/10.1029/2006GL027752.
Vincent, C. L., and T. P. Lane, 2016: Evolution of the diurnal precipitation cycle with the passage of a Madden–Julian Oscillation event through the Maritime Continent. Mon. Wea. Rev., 144, 1983–2005, https://doi.org/10.1175/MWR-D-15-0326.1.
Vincent, C. L., and T. P. Lane, 2017: A 10-year austral summer climatology of observed and modeled intraseasonal, mesoscale, and diurnal variations over the Maritime Continent. J. Climate, 30, 3807–3828, https://doi.org/10.1175/JCLI-D-16-0688.1.
Vincent, E. M., M. Lengaigne, C. E. Menkes, N. C. Jourdain, P. Marchesiello, and G. Madec, 2011: Interannual variability of the South Pacific Convergence Zone and implications for tropical cyclone genesis. Climate Dyn., 36, 1881–1896, https://doi.org/10.1007/s00382-009-0716-3.
Wang, S., and A. H. Sobel, 2017: Factors controlling rain on small tropical islands: Diurnal cycle, large-scale wind speed, and topography. J. Atmos. Sci., 74, 3515–3532, https://doi.org/10.1175/JAS-D-16-0344.1.
Wheeler, M. C., and H. H. Hendon, 2004: An all-season real-time multivariate MJO index: Development of an index for monitoring and prediction. Mon. Wea. Rev., 132, 1917–1932, https://doi.org/10.1175/1520-0493(2004)132<1917:AARMMI>2.0.CO;2.
Wilks, D. S., 2006: Statistical Methods in the Atmospheric Sciences. 2nd ed. International Geophysics Series, Vol. 100, Academic Press, 648 pp.
Yang, Y., and Y.-L. Chen, 2008: Effects of terrain heights and sizes on island-scale circulations and rainfall for the Island of Hawaii during HaRP. Mon. Wea. Rev., 136, 120–146, https://doi.org/10.1175/2007MWR1984.1.
Yang, S., and E. A. Smith, 2006: Mechanisms for diurnal variability of global tropical rainfall observed from TRMM. J. Climate, 19, 5190–5226, https://doi.org/10.1175/JCLI3883.1.
Zhang, C., 2005: Madden–Julian Oscillation. Rev. Geophys., 43, RG2003, https://doi.org/10.1029/2004RG000158.