Characteristic of the Regional Rainy Season Onset over Vietnam: Tailoring to Agricultural Application

Owing to its unique position within multiple monsoon regimes, latitudinal extent, and complex topography, Vietnam is divided into seven agroclimatic zones, each with distinct rainy season characteristics. Variation in the dominant rainfall system across zones affects the rainfall climatology, the primary water resource for regional crops. This study explores the creation of an agronomic rainy season onset based on high-resolution rainfall data for each agroclimatic zone for applications in an agricultural context. Onset information has huge practical importance for both agriculture and the economy. The spatiotemporal characteristics of zonal onset date are analyzed using integrated approaches of spatial and interannual variability, temporal changes, and estimation of predictability using teleconnection with Niño 3.4 sea surface temperature anomalies (SSTA) for 1980 to 2010. Results suggest that northern and southern zones experience regional onset dates in May, while the central zones experience rainy season onset in late August. The regional variability of rainy season onset is lower in a single dominant monsoon regime (northern and southern zones) and higher in latitudinally extended zones on the border of monsoon regimes (central zones). The interannual variation in rainy season onset date is found to be approximately 2 weeks across all agroclimatic zones. The significant negative trend in rainy season onset date is found for Central Coast and South Central Coast zones, suggesting that the onset date shifted earlier for the entire period. In the decadal scale, the zonal mean onset date shifted later in the Northwest and earlier in the Central Highlands. Out of the seven climate zones, a significant positive correlation is only noticed in the Central Highlands and South zones between zonal mean onset date and Niño 3.4 SSTA for Dec-Jan-Feb, suggesting the potential of seasonal scale predictability of rainy season onset date with respect to preceding El Niño-Southern Oscillation (ENSO) events.


Introduction
In countries where rainy season rainfall is the main water resource to meet the requirements for agricultural production, reliable determination of the rainy season onset is crucial for agricultural planning as onset dates influence the time of land preparation, sowing and transplanting dates of major crops, mobilization of seed/crop, manpower, equipment, and it reduces the risk of planting and sowing too late or too early [1]. Any irregularity in rainfall during this period affects the prospect of yield production because of delayed transplanting and immature growth of crop owing to water scarcity, heavy infestation by weeds, and outbreaks of diseases and insect pests [2]. Although farmers use a range of non-scientific traditional strategies and criteria (e.g., the observation of the behavior of some birds or insects flowering of certain trees) to predict onset, these methods are unreliable, especially as global climate change alters long-term climatic patterns. Therefore, a scientific definition of the rainy season onset can aid yield increases among local farmers [3].
Defining the onset of the rainy season is an extremely complex problem and, thus, various definitions have been proposed. These definitions can be broadly classified into two groups: regional to large-scale methods where parameters measure large-scale atmospheric dynamics and local-scale methods using fixed precipitation thresholds selected based on the climatology of rainfall over the region of interest [4]. The local scale definition has direct application for agricultural production which depends on rainfall requirements of specific crops. Therefore, this local scale definition is often referred to as agronomic as it is based on crop-relevant local-scale daily rainfall features rather than regional-scale atmospheric circulation changes [5]. An agronomic definition can be estimated by finding the first rainy day which exceeds a particular rainfall threshold without a potential occurrence of crop-threating dry spells thereafter to account for false onsets [5].
With spurring growth in agricultural productivity in the 1990s, the sector became a key part of Vietnam's economy, now contributing 18.4% of Vietnam's Gross Domestic Product, and employing 54% of the working population [6]. In recent times, Vietnam became a leading global exporter of several important crops such as rice, coffee, cashew nuts, vegetables, and rubber [7]. Agricultural practices in Vietnam are crucially dependent on rainfall. Although most of the gross cultivated area is irrigated, irrigation itself is crucially dependent on monsoon rainfall as it largely depends on the use of surface water. Therefore, the date of onset of the rainy season is of particular importance to the agriculture sector in Vietnam. However, despite the fact that the investigation of rainy season onset dates has a huge practical importance for both agriculture and economy, targeted studies evaluating the timing of the rainy season in Vietnam have been very limited [8,9,10].
Studies investigate the climatology, trends and the predictability of rainy season onset dates in southern Vietnam [8] and the Central Highlands region [9] of Vietnam. Both studies also explored prediction of rainy season onset dates based on statistical modeling using the relationship between local rainy season onset dates and large-scale atmospheric variables (pressure and moist static energy gradients, outgoing long-wave radiation, wind fields and mean sea level pressures) over certain regions. Only one study [10] discussed the climatological pattern of onset dates for all regions of Vietnam but it is limited only to the summer monsoon season where the rainy season of the Central Highlands region is mainly influenced by local rain-producing weather systems [11].
Existing studies are limited either for a particular region of the country or for summer monsoon, while the climate of Vietnam substantially varies from north to south with seven climatic sub-regions which are strongly affected by multiple monsoon systems and local factors (more discussions can be found in the following sections). Thus, a more in-depth study was required to reveal an explicit understanding of spatiotemporal characteristics of the rainy season onset of Vietnam. For applications in an agricultural context, this study mainly focuses on the agronomic rainy season onset which has huge importance for agricultural activities, at different the agroclimatic zones of Vietnam. Nevertheless, this study is motivated to raise the following questions: How is the agronomic rainy season onset distributed over agroclimatic zones of Vietnam? How does rainy season onset vary over the spatiotemporal scale for each zone? Is there any predictability based on the teleconnection between the onset days for each zone and with sea surface temperature indices? Therefore, this paper aims to answer these issues by exploring the creation of an agronomic rainy season onset definition for all agroclimatic zones over Vietnam using recently developed high resolution Vietnam-gridded precipitation (VnGP) dataset [12] for the period 1980 to 2010. The novelty of this research lies in the fact that it is the first time the agronomic rainy season onset for each agroclimatic zone across Vietnam is analyzed by using integrated approaches of inter-annual variability, trend analysis and estimation of predictability teleconnection. The findings of this study offer a potential understanding of agronomic definition for the local stakeholders and decision makers for agriculture as well as for the Vietnamese National Hydro-Meteorological Service (VNHMS).
The remainder of this paper is organized as follows: After a short description of the agroclimatic setting of Vietnam in section 2, the description of the dataset is presented in section 3. Agronomic rainy season onset definition calculation is shown in section 4. Spatiotemporal characteristics of onset dates of each zone including interannual variability, temporal changes, together with estimation of teleconnections with El Niño-Southern Oscillation (ENSO) presented in section 5. The conclusions of this study follow in section 6.

The Agroclimatic Setting of Vietnam
Vietnam, the easternmost country on the Indochina Peninsula in Southeast Asia, is bordered by China to the north, Laos and Cambodia to the west, the Gulf of Thailand to the south and the Gulf of Tonkin and the South China Sea to the east. From the Tropic of Cancer (approximately 23.26N) to the deep tropics (approximately 8.58N), Vietnam extends 1650-km north to south and has an extensive coastline of 3,260 km stretching from the Gulf of Tonkin in the north to the Gulf of Thailand in the south. This long, thin, S-shaped country has a very complex topography with the mountainous North, except for the Red River valley and the coastal plain; the central highlands and coast; and the low-lying South, home to the Mekong River Delta.
The total annual rainfall is distributed unevenly throughout the year, where approximately 70 percent of the rainfall occurs during the main rainy season from May to September/October [6]. The northern region (+20N) rainy season is from May to September with the highest monthly rainfall in August. The central region (11N-20N) rainy season has an extended bimodal rainy season from May to November with rainfall peaks in July and October. The rainy season shifts shorter and later in the northern part of this elongated region. The central region has varied and unusual monsoon characteristics because of the proximity of the central highlands, eastern coast, and the effects of the adjacent mountains in Cambodia and Laos to the west. The southern region (8N-11N) has a similarly extended rainy season from May to November with near constant high magnitude rainfall from June to October. Therefore, climate regimes in Vietnam are strongly affected by both monsoon systems and local factors [13]. Due to complex topography, an elongated latitudinal extent, and the influence of multiple monsoon systems, like South Asian summer monsoon, East Asian winter monsoon as well as North-Western Pacific monsoons, Vietnam's climate varies significantly from north to south with several sub-climate regions [13]. Vietnam is divided into seven terrestrial sub-regions with distinct climate characteristics based on the duration of the rainy season, the timing of peak rainfall, and the assessment of rainfall station data by VNHMS [10]. These zones are referred to as Northwest (denoted as Z1), Northeast (denoted as Z2), North Central (denoted as Z3), Central Coast (denoted as Z4), South Central Coast (denoted as Z5), Central Highlands (denoted as Z6), and South (denoted as Z7) (Figure 1.a).
Almost 35% of the total national land area, including arable land, permanent crop land, and permanent meadow land, is used for agricultural production. Nationally, rice, coffee, tea, and pepper are the predominant crops, but the choice of the crops is specialized by agroecological zone [7]. The Northwest and Northeast zones are mountainous areas with often deficient transportation facilities, poor market access, and limited irrigation systems. The agricultural production in the northern part of Vietnam is mostly for industrial crops (tea and rubber) and food subsistence purposes considering that these areas are mountainous with limited irrigation systems. Cash crops (mostly coffee) are mainly produced in the Central Highlands and the South zones. Rice production mostly is concentrated in the two delta regions (Red River Delta in the North Central zone and Mekong River Delta in the South zone). Agricultural production accounts for up to 95% of total water withdrawals in Vietnam while 49% of the total agricultural land is serviced by an irrigation scheme which is primarily designed for rice in the two deltas [7]. Irrigation is crucially dependent on monsoon rainfall because it largely depends on the use of surface water. Any aberration in rainfall during the onset time affects the prospects of a good yield as farmers make many agricultural decisions based on the local monsoon onset, including sowing dates and fertilizer timing. Therefore, it is very important to investigate the rainy season onset characteristic focusing on each sub-climate region.

Dataset
The gridded data product named Vietnam-gridded precipitation data (VnGP) is used for all the rainy season onset calculations. The dataset was built from Vietnamese National Hydro-Meteorological Service's 481 rain gauge stations. The Spheremap interpolation method, a modified Shepard's interpolation for spherical coordinates, was used to construct the gridded products from rain gauge stations. These gridded products are generated and archived by the Department of Meteorology and Climate Change, VNU University of Science, Vietnam at daily time scale with Preprints (www.preprints.org) | NOT PEER-REVIEWED | Posted: 8 December 2020 doi:10.20944/preprints202012.0202.v1 spatial resolutions of 0.1° and 0.25°, covering the period from 1980 to 2010. In the present study, the 0.25° version for the entire 31 years (1980 to 2010) is used. The detailed procedure of the generation of VnGP has been discussed in detail by [12] and data are available at http://danida.vnu.edu.vn/cpis/en/content/gridded-precipitation-data-of-vietnam.html To investigate the influence of El Niño-Southern Oscillation (ENSO) over zonal rainy season onset, the sea surface temperature anomalies (SSTA) is obtained from the second version of the optimum interpolation (OI) analysis (OIv2) produced by the National Oceanic and Atmospheric Administration (NOAA) [14]. The OIv2 has a spatial resolution of 1° and a daily timescale, with monthly means from December 1981 to January 2020, ongoing. In this study, we use the monthly mean SST anomaly from 1981 to 2010 averaged over the Niño 3.4 region of the Pacific Ocean (120ºW-170ºW, 5ºS-5ºN). The Niño 3.4 SSTA data is obtained from the Data Library of International Research Institute (https://iridl.ldeo.columbia.edu/SOURCES/.Indices/.nino/.NCEP_OIv2/.NINO34/index.html#info).

Defining Agronomic Rainy Season Onset based on Local-Scale Climatic Conditions
In this work, the local rainy season onset for each agroclimatic zone is calculated using the agronomic definition proposed by [5]. As described in Section 1, a local definition like the one used here is of great relevance and can have numerous applications in the development of operational products in agriculture. This agronomic definition is advantageous in regions where dry and rainy seasons are well defined [4]. As all the Vietnamese agroclimatic zones have well-defined dry and rainy seasons, this agronomic definition is advantageous in Vietnam. This agronomic definition, originally developed using historical precipitation from India, allows for the calculation of monsoon timing for any location using five basic parameters: the amount of rainfall during the first and last wet spell, its duration, the duration and intensity of the post-onset dry spells, and the length of the period in which these dry spells are searched. In the present study, the original parameters were adjusted to satisfy agronomic definition parameters based on the climatic features of Vietnam. All parameters for the initial wet spell are based on discussions with local experts at VNU University of Science and Vietnamese National Hydro-Meteorological Service (VNHMS). In short, this method defines onset as the first wet day of the first 5-day period with an average daily rainfall equal or larger than the climatological 5-day wet spell for the season and without a 7-day dry spell (no wet days) during the following 20 days to account for false onsets.
As the focus of this paper is the rainy season, which is not limited by monsoon for the central highlands, a wet day is considered any day with at least 1 mm/day to account for sufficient rainfall to penetrate the soil and reach planted seeds. A requirement for five days of 5mm/day rainfall is translated to 25 mm of rainfall in five days with at least four wet days. The reason behind choosing a higher rainfall threshold to define a wet day rather than the commonly used wet day definition (>0.1mm/day) is a low rainfall threshold makes agronomic definition less suitable for agricultural decisions [5].

Search Window for Onset Definition
Finding the first wet day of the period which meets the subsequent onset date conditions (defined in Section 4) is necessary for the calculation of agronomic onset date in a given year; thus, determining the date of the beginning of the search the window in which to look for a day which satisfies the onset definition is a crucial step. In this study, the climatological mean monthly rainfall pattern is used to determine the beginning of the rainy season for each agroclimatic zone. As mentioned in Section 2, the predominant and intersecting monsoon systems over the different regions of Vietnam influence the onset of the rainy season across the agroclimatic zones. The East Asian Monsoon is the primary monsoonal influence on the Northern zones of Vietnam, while the South Asian monsoon dominates the southern zones. The central coastal zones and the central highlands are located between these two systems and, thus, are dominated by local rainfall patterns. The variation in the dominant rainfall system across zones affects the monthly rainfall climatology and contributes to the regional primary crop; for example, coffee in the central highlands is suited to more dispersed rainfall throughout the year with a later onset of the rainy season. The monthly mean rainfall for each agroclimatic zone along with country average is shown in Figure 1. It is seen that the main rainy season for northern and southern zones is during May to September, while the peak in rainfall is in July to August. For central zones it shows a decrease in rain relative to the rest of the country in the summer months and a later rainy season from August to December with a peak in rainfall in October. Further, the figure shows that the mean monthly rainfall for the whole country of Vietnam has an extended rainy season of moderate rainfall which slowly increases from May to October, and a period of low rainfall from December to April. For the 1980-2010 time period, the majority of agroclimatic zones, excluding the northwestern mountains and the central coast, experience rainy season onset after mid-April and, thus, the search period begins around this time for these zones (Z2, Z3, Z6, Z7). The Northwest (Z1) has a consistently earlier rainy season onset, while the Central Coast and South Central Coast zones (Z4, Z5) are found to have a significantly delayed rainy season onset, and thus the search window shifts forward and backward, respectively, for these zones. Additionally, Summer monsoon rainfall in Vietnam has a wet-dry phase shift of 2-10 days [10] and, thus, the dry spell period for the false start parameter is seven consecutive days within the twenty-day period following the first wet day. This requirement is to avoid false-start onsets. The search window for the rainy season onset date for each zone is 60-days to capture interannual rainy season onset date variation.
The onset date coverage is determined by the percentage of years that satisfy the agronomic onset definition criteria within the search period for each grid point for the entire 31-year study period ( Figure 2). Years that do not satisfy the agronomic definition and, thus, are not included, may have had an anomalously early or late onset or unusual rainfall patterns. High onset date coverage suggests suitability of the agronomic definition to each of the agroclimatic zones. The initial date for the search window was determined using this method of testing possible dates and assessing suitability based on coverage. The earliest date which had high coverage was deemed suitable to capture the onset of the rainy season while representing its variability of the 31-year time period. The Northwest zone (Z1) is found to have the earliest zonal rainy season onset date and the beginning day of the search window for the rainy season is set to April 15 th (Figure 2a) to fully capture the interannual variability of the rainy season onset. Other beginning days of the search window were tested and found that the April 25 th date (Figure 2b) captures the largest spread of onset date for the majority of zones (Z2, Z3, Z6, Z7). Lastly, due to latitudinal extent and unique topography, the Central Coast and South Central Coast zones (Z4, Z5) are found to have a much more variable rainy season and a later onset date. The beginning day of the search window for these zones is determined to be August 5 th (Figure 2c).

Spatiotemporal Characteristics of Regional Onset Date
Following the above-mentioned procedure, for each agroclimatic zone, the onset dates are estimated for each grid point independently for each year during the time period . The zonal mean of the onset dates for each agroclimatic zone are estimated by taking the mean of the onset dates across the time period (31 years) for each grid point within that agroclimatic zone, and, then, averaging the mean onset date across all grid points. It is found that the zonal mean of rainy season onset for the Northeast (Z2), North Central (Z3), and South (Z7) zones is in mid to late May (May 20 th , May 21 th and May 14 th respectively). However, for the Northwest (Z1) and Central Highlands (Z6), the mean onset days are in early May (May 9 th and May 11 th , respectively). These onset dates are closely related to the summer monsoon system (as described in sec.5.1). Although these results are in agreement with earlier studies [8;10], they are not exactly the same as a different definition is used to determine the onset dates. Lastly, the Central Coast and South Central Coast zones (Z4, Z5) have a noticeably late rainy season onset in late August (August 24 th and August 29 th , respectively). This is likely due to their unique topography, between the Cambodian mountain range to the west and the sea to the east, and position on the border of interacting monsoon systems. The zonal spatial variability of onset dates, or the variation of mean onset date between grid points across each agroclimatic zone, is calculated to assess how coherent the onset dates are across each agroclimatic zone. The spatial variability is calculated by taking the standard deviation of all mean onset dates (mean of the onset dates across the time period for each grid point) across all grid points in a given agroclimatic zone. It is expected that the mean onset dates occur at a similar time scale across the agroclimatic zone (i.e., low zonal spatial variability of onset dates); however, it is found that the mean onset dates vary from 4 to 8 days within each zone. Moreover, spatial variability is lower in more compact zones with a single dominant monsoon regime (Z1, Z2, Z3, Z6, Z7) at around 4-6 days and higher for latitudinally extended zones on the border of monsoon regimes (Z4, Z5) at 7-8 days.
The zonal interannual variability of onset date, or temporal variability, is also estimated to illustrate the year-to-year variation of onset dates in a given agroclimatic zone. This is calculated by taking the standard deviation of the onset dates across the time period for each grid point in a given zone, and, finally, averaging the standard deviation of onset date across all grid points. Interannual temporal variability of rainy season onset was found to be approximately two weeks across all agroclimatic zones. Specifically, the majority of zones (Z1, Z2, Z3, Z6, Z7) have an interannual variability between 11 to 13 days, while the Central Coast and South Central Coast zones (Z4, Z5) have a slightly higher interannual variability of more than 14 days. This is notable because an interannual variation of approximately two weeks presents a challenge to the prediction of rainy season at a seasonal timescale. Table 1 summarizes the mean rainy season onset day, interannual temporal variability, and spatial variability for each climatic zone.

Temporal Changes of Regional Onset Date
To study the temporal changes in the onset dates for each agroclimatic zone across the time period, the Mann-Kendall (MK) test [15], a rank-based non-parametric test, is performed on the regionally averaged onset dates across all grid points (average onset dates over all grid points for each year) at the 90% statistical significance level. The MK test is widely used for trend analysis because of its robustness for non-normally distributed data and low sensitivity to an abrupt change. However, the Mann-Kendall test gives the direction (a positive or negative value of MK's test statistic indicates the rainy season onset date is trending later or earlier respectively) but not the magnitude of the significant trends. Therefore, in addition, the Sen's slope estimator [16], a non-parametric test to identify the slope of the trend, is calculated to estimate the magnitude of change. The MK's test statistic and magnitude of Sen's slope for the entire study period  and each decade (1980s, 1990s, 2000s, respectively) for each agroclimatic zone is shown in Table 2. It is found that for the entire 31 years, a significant negative trend in rainy season onset date is only observed for the Central Coast and South Central Coast zones (Z4, Z5) out of all seven zones. The negative trend for both the zones represent that the onset dates shifted earlier across the entire period. The Sen's slope for the South Central Coast (Z5) has a greater magnitude than that of the Central Coast (Z4), 0.44 compared to 0.34, 8 of 12 which shows that the annual zonal mean rainy season onset date in the South Central Coast (Z5) shifted earlier more rapidly than in the Central Coast (Z4). Mostly there is no significant trend found for individual decades for each zone except for the Northwest and Central Highlands (Z1, Z6) in the 1990s. For the second decade (1991)(1992)(1993)(1994)(1995)(1996)(1997)(1998)(1999)(2000), a significant positive and negative trend with Sen's slope 1.647 and -2.011 are observed for the Northwest and Central Highlands (Z1, Z6), respectively. This means that the zonal mean onset date shifted later in the northwest and earlier in the Central Highlands in this decade. For a visual representation, the yearly zonal mean onset dates with a linear trend line for the entire period (198-2010) and each decade is presented for each agroclimatic zone in Figure 3.

Seasonal Predictability of Regional Onset Date
Prediction of the rainy season onset would be crucial for agricultural and economic practices in each agroclimatic zone for Vietnam. The influence of the Pacific Ocean climate drivers, namely El Niño-Southern Oscillation (ENSO), on the rainfall in Vietnam has been studied extensively [17; 18]. In this section, we explore the teleconnections between the tropical Pacific SST (as represented by the Niño 3.4 SSTA) anomalies and the rainy season onset for each agroclimatic zone to investigate potential sources of seasonal predictability. For the Northern (Z1, Z2 and Z3), and Southern zones (Z6, Z7) where the mean onset dates are in May, Pearson's correlation coefficient (r), a measure of linear relationship, is calculated between zonal onset date and Niño 3.4 SSTA for the seasonal mean of December to February (DJF) period. As the onset date for the Central Coast and South Central Coast zones (Z4, Z5) are in late August, Niño 3.4 SSTA for the seasonal mean of June to August (JJA) period is used to calculate the same correlation coefficient. The student's T test is performed to determine statistical significance of the correlation coefficient at the 95% level. The scatter plot between zonal onset dates and Niño 3.4 SSTA along with the Pearson's correlation coefficient is shown in Figure 4. Out of the seven agroclimatic zones, a significant positive correlation is only noticed in the Central Highlands and the Southern zones (Z6, Z7) between zonal mean onset date and Niño 3.4 SSTA for DJF. This suggests that the El Niño (La Niña) events result in a delayed (early) Preprints (www.preprints.org) | NOT PEER-REVIEWED | Posted: 8 December 2020 doi:10.20944/preprints202012.0202.v1 rainy season onset date in these two agroclimatic zones. This finding suggests the potential of seasonal prediction of rainy season onset date in these zones with respect to preceding ENSO events. As both the zones are significantly important for agricultural production in Vietnam (the South zone is the location of the Mekong River Delta) in terms of coffee and rice cultivation, seasonal prediction of rainy season onset would be very useful to many Vietnamese farmers in these agroclimatic zones.

Concluding Remarks
The primary goal of this study is to define a rainy season onset tailored to an agricultural application over the seven primary agroclimatic zones of Vietnam and assess its spatiotemporal characteristics by using integrated approaches of spatial and interannual variability, temporal changes, and estimation of predictability using teleconnection with El Niño-Southern Oscillation (ENSO) for 1980 to 2010. The novelty of this research is that a tailored agronomic definition based on high-resolution rainfall data is used for the first time over Vietnam at a regional climate scale. The results of this study would be highly relevant to local farmers and decision makers in Vietnam where the productivity and profitability of the agricultural sector is highly dependent on rainfall and the timing of the rainy season.
In summary, the major findings of the present study are:  Rainy season onset date varies across the agroclimatic zone. For the Northern and Southern zones, the zonal mean of rainy season onset dates is primarily in May which associated with summer monsoon system. Due to their unique topography and position on the border of interacting monsoon systems, Central Coast and South Central coast have rainy season onset in late August.  The zonal spatial variability of onset dates, define as the standard deviation zonal mean onset dates across all grid points in a given agroclimatic zone, is vary from 4 to 8 days within each agroclimatic zone. However, the variability is lower in a single dominant monsoon regime (northern and southern zones) and higher in latitudinally extended zones on the border of monsoon regimes (Central Coast and South Central Coast).  The interannual variation in rainy season onset date is approximately 2 weeks (11 to 14 days) across all agroclimatic zones of Vietnam. This is result suggest that the prediction of rainy season onset is highly challenging work due to such high temporal variation.  The trend analysis for the zonal onset dates of each agroclimatic zone has been analyzed by using the Mann-Kendall (MK) Test to gather with Sen's Slope Estimator for the entire 31-year period and each decade. The rainy season onset date is significantly trended earlier in the Central Coast and South Central Coast zones. For 1990s (second decade) the zonal mean onset date shifted later in the northwest and earlier in the central highlands at 90% significance level.  The teleconnections between Niño 3.4 SSTA for Dec-Jan-Feb season and onset date in the Central Highlands and South zones is found (positive) to be statistically significant at 95% level. This finding suggests the potential of seasonal prediction of rainy season onset date in these zones with respect to preceding ENSO events. Considering the fact that the seasonal prediction of rainy season onset at each climate zones would be very beneficial to the agriculture sector; the future scope of the present work includes the development of an operational forecast system using state-of-art general circulation models outputs. Funding: This research work and the APC was funded by ACToday (https://iri.columbia.edu/actoday/), the first of Columbia University's World Projects (https://worldprojects.columbia.edu/). E.B. acknowledges the support from the Climate and Society MA program at Columbia University.