Water vapor transport for spring persistent rains over southeastern China based on five reanalysis datasets

In this study, atmospheric water vapor transport was analyzed to study the changes in spring persistent rainfall (SPR) over southeastern China from 1980 to 2012. The performances of five sets of reanalysis data in reproducing the climatology, the long-term trend and interannual variability of the SPR were synthetically evaluated. To understand the mechanisms dominating SPR variation, the major components of moisture budget, including vertical moisture advection, horizontal moisture advection and evaporation, were examined based on the five reanalysis datasets. The results show that all five reanalysis datasets reproduce the climatology of the SPR reasonably well. Strong westerly wind flow over the southern Tibetan Plateau and southwest wind flow over the western North Pacific are the two main channels remotely supplying water vapor for the SPR. Locally, moisture budget diagnosis shows that the SPR is primarily contributed by the evaporation and vertical moisture advection terms. The SPR shows a decreasing tendency [−0.38 mm day−1 (10 year)−1] during 1980–2012 along with strong interannual variation, which is reasonably captured by all five reanalysis datasets. Vertical moisture advection dominates the SPR variation, and the decreasing trend of the SPR is mainly due to the weakening of ascending motion during this period. In an El Niño decaying spring, the anomalous lower-tropospheric anticyclone over the western North Pacific intensifies vertical upward motion and leads to more precipitation. Both the Japanese 55-year reanalysis (JRA55) and the European Centre for Medium-Range Weather Forecast interim reanalysis (ERAIM) show higher skill in reproducing the climatology and changes of the SPR.


Introduction
Spring persistent rainfall (SPR) is a unique synoptic and climatic phenomenon in East Asia, spanning the 13th to 27th pentads of the solar year (Tian and Yasunari 1998). The distribution and intensity of the SPR rain belt are subject to the topography of the "Nanling" (24-26.5°N, 110-116°E) and "Wuyi" (25-29°N, 116-119°E) Mountains Wu 2007, 2008;Zhang et al. 2013). Over southeastern China, boreal spring rainfall, featured by high frequency (Huang et al. 2015b), accounts for approximately 35% of the annual accumulation (shown in Fig. 1b). As a climatic phenomenon, the SPR is the major rainy period before the Meiyu/Baiu in early summer over the middle and lower reaches of the Yangtze River. Spring is the important period in the transition of atmospheric circulation from winter monsoon state to summer, significantly influencing the onset of the summer monsoon Pan and Jiang 2014;Wu et al. 2015). Since southeastern China is one of the most important economic agricultural regions, Abstract In this study, atmospheric water vapor transport was analyzed to study the changes in spring persistent rainfall (SPR) over southeastern China from 1980 to 2012. The performances of five sets of reanalysis data in reproducing the climatology, the long-term trend and interannual variability of the SPR were synthetically evaluated. To understand the mechanisms dominating SPR variation, the major components of moisture budget, including vertical moisture advection, horizontal moisture advection and evaporation, were examined based on the five reanalysis datasets. The results show that all five reanalysis datasets reproduce the climatology of the SPR reasonably well. Strong westerly wind flow over the southern Tibetan Plateau and southwest wind flow over the western North Pacific are the two main channels remotely supplying water vapor for the SPR. Locally, moisture budget diagnosis shows that the SPR is primarily contributed by the evaporation and vertical moisture advection terms. The SPR shows a decreasing tendency [−0.38 mm day −1 (10 year) −1 ] during 1980-2012 along with strong interannual variation, which is reasonably captured by all five reanalysis datasets. Vertical moisture advection dominates the SPR variation, it is necessary to understand variations of the SPR and responsible mechanisms from the climate respective.
Much effort has been devoted to studying the SPR. Tian and Yasunari (1998) first proposed the concept of "SPR". They found that the east-west thermal contrast between the Indochinese Peninsula and the western North Pacific, which is due to the time-lag of seasonal warming in the spring, is the main mechanism of SPR formation. The mechanical and thermal forcing of the Tibetan Plateau (TP) also plays primary roles in the formation of the SPR (Wan and Wu 2007;Wan et al. 2009). For interannual variability, previous studies have suggested a strong linkage between the intensity and frequency of the SPR over southeastern China and the East Asian subtropical jet (i.e., EASJ; Wen et al. 2007;Deng et al. 2014). The oceanic thermal condition over the equatorial Pacific and the South China Sea can also influence the variability of the SPR (Yang and Lau 2004;Qiang and Yang 2013;Chen et al. 2014;Wu et al. 2015).  suggested that there is a strong relationship between the Multivariate ENSO Index in the preceding winter and the SPR. In an El Niño decaying spring, an anticyclonic circulation anomaly still exists in the middle-low-level of the troposphere over the western North Pacific. These circulation anomalies favor strengthened southwestern warm and moist flow and enhanced convergence of water vapor fluxes over southeastern China, leading to a significant increase in precipitation . In recent decades, southeastern China has experienced a significant decrease in precipitation in the boreal spring and has frequently suffered from persistent drought, and the droughts are consistent with anomalous overturning circulation and local anomalous descent motion (Xin et al. 2006;Sun and Yang 2012;Feng et al. 2014;Zhu et al. 2014;Li et al. 2016), but the long-term behavior of SPR precipitation in different reanalysis datasets remains unestablished.
Various reanalysis datasets have been extensively used in climatological studies. The reanalysis datasets are created via an unchanging data assimilation scheme and the most advanced operational numerical models. Under the process of data assimilation by observations, the atmospheric circulation in the reanalysis datasets is usually considered the optimal estimate. Comparing different datasets can provide robust knowledge on actual climate variability and climate change. For example, Bosilovich et al. (2008) used five reanalysis datasets to evaluate global precipitation in different regions. Lin et al. (2014) evaluated the changes in global monsoon precipitation by using multiple reanalysis datasets. Ma and Zhou (2016) compared seven reanalysis datasets and found that the tropical Pacific Walker Fig. 1 Climatology of spring persistent rains (SPR) over southeastern China based on GPCP data and Chinese rain gauge observation dataset: a spatial distributions of SPR based on GPCP data (unit: mm day −1 ); b the same as a, but derived from Chinese rain gauge observation dataset; c percentage of climatological SPR in the annual total based on GPCP (unit: %); d the same as c, but for Chinese rain gauge observation dataset. The rectangle indicates the major rainfall center circulation has strengthened and shifted westward in the recent decades. Previous studies often considered one or two reanalysis datasets as the optimal observation. At present, the uncertainty among different reanalysis datasets in describing the SPR change remains unknown. In this study, we use five reanalysis datasets to estimate uncertainties. We aim to answer the following questions: (1) how well can five reanalysis datasets capture climatology, the long-term trend and interannual variability of the SPR during the past 30 years? (2) Which components of the moisture budget dominate climatology and changes in the SPR? We will show that vertical moisture advection dominates the SPR variation, both JRA55 and ERAIM show relatively higher skill in reproducing the climatology and changes of the SPR and are thus recommended for climate variability and change studies. The remainder of this paper is organized as follows. The observational datasets, five reanalysis datasets, and analysis methods are described in Sect. 2. In Sect. 3, we evaluate the performance of the five reanalysis datasets in reproducing the climatology and changes of SPR, and investigate the dominant component of the moisture budget for the aspect of climatology, the long-term trend and interannual variability of the SPR. Finally, a summary is given in Sect. 4.

Data description
The observational monthly precipitation data used in this study is the Global Precipitation Climatology Project (GPCP) version 2.2 dataset (Huffman et al. 1997;Adler et al. 2003). GPCP v2.2 is derived from a mix of satellite estimates of precipitation over ocean and land, and then adjusted by rain gauge measurements from land and atoll stations. It is available monthly and covers 1980-2012. The spatial resolution of GPCP v2.2 is 2.5° × 2.5°. We also use rain gauge observation data over China to make the results based on GPCP convincible. Rain-gauge data of daily precipitation from 1980 to 2012 were obtained from the China Meteorological Administration (CMA; http://data.cma.cn/en). The CMA dataset contains 756 stations and is updated through recent years. To facilitate analysis, the rain gauge data were interpolated onto a 1.5° × 1.5° grid.
The five reanalysis datasets used in this study include monthly precipitation, evaporation, surface pressure, vertical velocity, zonal and meridional components of wind and specific humidity: These reanalysis datasets differ in the atmospheric forecast model, data assimilation method, model resolution, and model physics. Their characteristics are summarized in Table 1. Each reanalysis dataset has a different time span but all cover our analysis period of 1980-2012. Results of each reanalysis are derived from their original horizontal resolution. The ensemble mean of five reanalysis datasets is calculated after each reanalysis dataset is interpolated into a global 2.5° × 2.5° grid using conservative interpolation (using the local area-conservative  1979-2012 1958-present 1979-present 1900-2010 1979-present  In addition, we use the ENSO index derived from the CPC in National Oceanic and Atmospheric Administration (NOAA; information available online at http://www. cpc.noaa.gov/products/analysis_monitoring/ensostuff/ensoyears.shtml). Warm and cold periods are determined based on a threshold of ±0.5 °C for the 3-month running mean of Extended Reconstructed Sea Surface Temperature version 4 (ERSST v4; Huang et al. 2015a) anomalies in the Niño3.4 region (5°N-5°S, 120°-170°W), and based on the 1981-2010 period.

Analysis methods
To understand the mechanisms dominating SPR variation over southeastern China, we perform a moisture budget analysis following previous studies (Held and Soden 2006;Seager et al. 2010;Chou and Lan 2012;Huang et al. 2013;Feng and Zhou 2012;Lin et al. 2014;Ma and Zhou 2015). The moisture budget equation is (Eq. 1), where P is precipitation, E is evaporation, q is specific humidity and ⇀ V is wind vector. is a residual term including transient eddies (Trenberth and Guillemot 1995;Zhou and Yu 2005) and contributions from surface processes due to topography (Seager et al. 2010). All values in Eq. 1 are a 3 month mean of March-April-May (MAM). "⟨⟩" denotes a vertical mass integration through the whole troposphere (Eq. 2), where g is gravitational acceleration, p s is surface pressure and p T is the pressure of the tropopause, taken as 100 hPa. − t ⟨q⟩ in Eq. 1 is the time derivative of specific humidity q, which can generally be ignored as its seasonal mean value is much smaller than other terms.−∇ ⋅ (q ⇀ V) is the convergence of integrated moisture flux. According to the mass conservation equation, −∇ ⋅ (q ⇀ V) can be divided into two terms (Chou and Lan 2012;Chou et al. 2013): a vertical moisture advection term −⟨ q p ⟩ (simply written as − ⟨wdq⟩) and a horizontal moisture advection term −⟨ ⇀ V h ⋅ ∇ h q⟩ (simply written as − ⟨vdq⟩). Therefore, Eq. 1 can be reformulated as (Eq. 3), where is the vertical pressure velocity and ⇀ V h is the horizontal wind vector. If is assumed to be zero at the surface and at the tropopause (i.e., ⟨ p q⟩ = 0), then from the mass conservation equation we can obtain , so the vertical moisture advection is equal to the horizontal flow convergence of moisture term (Seager et al. 2010). Vertical moisture advection is mainly contributed by the convergence in the lower troposphere that promotes the vertical transport of moisture and then leads to rain formation.
To understand mechanisms that induce SPR anomalies, according to Eq. 3, the precipitation anomalies can be decomposed into (Eq. 4), Here, ()′ represents the departure from climatology (see definition below). The residual term ′ involves transient and nonlinear processes.
Following previous studies (Chou and Lan 2012;Chou et al. 2013;Huang et al. 2013;Ma and Zhou 2015), the changes of vertical moisture advection − ⟨ p q ⟩ can be further divided as (Eq. 5), where (̄) denotes the climatology from 1980 to 2012. The first term on the right of Eq. 5 represents changes in q with unchanged, conventionally called the thermodynamic component of − ⟨ p q ⟩ � , contributed by the changes in water vapor (simply written as − ⟨wdq � ⟩); the second term involves changes in with q constant, associated with changes in pressure velocity, which is mainly induced by atmospheric circulation changes (simply written as− ⟨w � dq⟩), and often called the dynamic component of − ⟨ p q ⟩ � ; the third term is a nonlinear term that is induced by changes in q and , and is much smaller than the thermodynamic and dynamic components. As in a previous study (Chou et al. 2009(Chou et al. , 2013Ma and Zhou 2015), the contribution of horizontal moisture advection to precipitation changes is relatively small, so it is not divided into three terms as the vertical one. Therefore, Eq. 4 can be reformulated as Eq. 6,

Climatology of the SPR
In this section, the mean state of the SPR and atmospheric water vapor transport from five reanalysis datasets from to 2012 are evaluated. The GPCP precipitation shows that during MAM, the major rainfall center is in southeastern China, where precipitation intensity could be as large as 7 mm day −1 (Fig. 1a). Furthermore, the precipitation in MAM accounts for approximately 30-35% of the annual accumulation (Fig. 1b), which means that the SPR contributes considerably to the total annual precipitation. Both observational rain gauge data and GPCP reveal coincident results with small disagreements, indicates that the results based on GPCP is reliable. We use GPCP to carry on our further study, because GPCP contains precipitation over ocean and the resolution of GPCP is more comparable with reanalysis datasets. The box (20-32°N, 108-122°E) in Fig. 1 enclosing the major rainfall center is used to perform the moisture budget analysis. The climatological mean of the SPR over southeastern China for the five reanalysis datasets and their ensemble mean are shown in Fig. 2. Table 2 shows the area-averaged precipitation intensity (mm day −1 ) and pattern correlation coefficients (PCCs) over the SPR major rainfall center for the five reanalysis datasets. The five reanalysis datasets reproduce the spatial patterns of GPCP reasonably well. PCCs are 0.66, 0.61, 0.62, 0.49 and 0.52 for JRA25, JRA55, ERAIM, ERA20C and NCEP2, respectively. All The mean state of the SPR and bias over southeastern China a JRA25, b ERAIM, c NCEP2, d JRA55, e ERA20C, f the multi reanalysis ensemble mean results, g JRA25 bias, h ERAIM bias, i NCEP2 bias, j JRA55 bias, k ERA20C bias, l the multi reanalysis ensemble mean bias. The bias is calculated as the difference with GPCP. The rectangle indicates the major rainfall center Table 2 The area-averaged precipitation intensity (mm day −1 ) and pattern correlation coefficients (PCCs) over the SPR major rainfall center from five reanalysis datasets and the ensemble mean for climatology and the El Niño decaying spring (i.e., 1983, 1987, 1992 and 1998 five reanalysis datasets show a westward shift of the rainfall center. The western boundary (defined as the longitude of the 5 mm day −1 contour line) in the reanalysis datasets extends to almost 105°E, approximately 5° westward compared to observation (~110°E). In addition, ERAIM, ERA20C and NCEP2 have an obvious northward shift in reproducing the spring precipitation. The JRA55 is superior in reproducing rainfall intensity averaged over the major rainfall center among the five reanalysis datasets. The ensemble mean shows the best performance in reproducing rainfall spatial distribution (PCC = 0.68), higher than any single reanalysis dataset. The climatological water vapor transport that contributes to the SPR is shown in Fig. 3. All five reanalysis datasets indicate that there are two main channels remotely supplying water vapor for the SPR over southeastern China. One is the water vapor transported by strong westerly wind flow over the southern Tibetan Plateau (TP); the other water vapor source is from southwestern wind flow over the western North Pacific (WNP). These two moisture sources merge approximately at 110°E and converge over southeastern China, providing adequate water vapor for the SPR. The spatial distribution of integrated water vapor convergence is similar to that of precipitation.
To understand the mechanisms dominating local SPR precipitation, the area-averaged precipitation, as well as other components of the moisture budget in MAM over the major rainfall center is shown in Fig. 4. All five reanalysis datasets show that locally the SPR is primarily contributed by evaporation and vertical moisture advection component. The pattern of evaporation is uniform with an intensity of approximately 3 mm day −1 (Fig. 4c). Four reanalysis datasets (JRA25, JRA55, ERAIM and ERA20C) show similar results in reproducing evaporation, except for NCEP2. The evaporation in NCEP2 is approximately 0.8 mm day −1 stronger than that of the other four reanalysis datasets (Fig. 4a). The pattern of vertical moisture advection is uneven, anchored by the "Nanling" (24-26.5°N, 110-116°E) and "Wuyi" (25-29°N, 116-119°E) Mountains (Fig. 4d), and similar to the precipitation distribution, indicating that the "Nanling-Wuyi" Mountains in southeastern China can enhance the convergence of water vapor flux and then favor SPR formation. The vertical moisture advection in ERA20C is approximately 0.9 mm day −1 weaker than that of the other four reanalysis datasets (Fig. 4a). The horizontal moisture advection and residual components contribute little to the total SPR precipitation.

The long-term trend of the SPR
The long-term trend of the SPR precipitation shows a significant decrease (Fig. 5a, b). In observations, southeastern China has experienced a significant drying trend over the past 33 years . The land and coastline in southeastern China have suffered the most obvious decrease, with most of this region exhibiting a significant negative trend of more than 0.6 mm day −1 (10 year) −1 . The patterns of the trend derived from the ensemble mean of the five reanalysis datasets are consistent with observation. However, more regions in southeastern China experienced a decreasing trend in the reanalysis datasets compared to the observation. The time series of the SPR precipitation anomalies averaged over the major rainfall center are shown in Fig. 5e. All five reanalysis datasets are consistent with observation. The observed SPR precipitation shows a negative tendency of −0.38 mm day −1 (10 year) −1 at a 5% significance level from 1980 to 2012. The trend is −0.79, −0.46, −0.41 and −0.65 mm day −1 (10 year) −1 in the JRA25, JRA55, ERAIM and NCEP2 reanalysis datasets, respectively, all exceeding the 1% significance level (Table 3). The trend in ERA20C shows a weaker decreasing trend [−0.26 mm day −1 (10 year) −1 ], which is still statistically significant at the 10% level. Both JRA55 and ERAIM are more similar to observation than the other reanalysis datasets (Table 3). In additional to the long-term trend, observation also shows a robust year by year variability. JRA55 shows the highest correlation coefficient (0.96) with observation, and the ensemble mean of five reanalysis datasets also shows a relatively high correlation coefficient (0.92). NCEP2 shows the lowest skill among the datasets, but the correlation coefficient (0.78) is still statistically significant at the 1% level (Table 3).
To better understand the long-term trend of SPR precipitation, we performed a moisture budget analysis to examine the trend of each term in the moisture budget. The linear trend of SPR precipitation and other moisture budget components are shown in Fig. 6. The evaporation in all five reanalysis datasets shows a weak increasing trend from 1980 to 2012. The decreasing trend of SPR precipitation is dominated by change in the vertical moisture advection, which shows a significant decreasing trend in all five reanalysis datasets. The PCC between these two terms for the ensemble mean is 0.83. We decompose the change in vertical moisture advection into three terms (Eq. 5), and find that the dynamic term dominates the change in vertical moisture advection. The spatial pattern of the linear trend for the Fig. 4 Climatological SPR precipitation and other moisture budget components in GPCP data and reanalysis datasets (unit: mm day −1 ): a moisture budget components averaged over the major rainfall center, b multi reanalysis ensemble mean precipitation, c multi reanalysis ensemble mean evaporation, d multi reanalysis ensemble mean vertical moisture advection (a) (b) (c) (d)  The time series of vertical moisture advection anomalies and thermodynamic and dynamic components are shown in Fig. 7. The vertical moisture advection anomalies term in the ensemble mean of the five reanalysis datasets shows a negative trend at a rate of −0.50 mm day −1 (10 year) −1 from 1980 to 2012, exceeding the 1% significance level (Fig. 7a). We further examine the long-term trends of the thermodynamic and dynamic components. The dynamic term shows a negative trend of −0.48 mm day −1 (10 year) −1 (Fig. 7b), but the thermodynamic term shows no significant long-term trend (Fig. 7c). Hence, the negative long-term trend of SPR from 1980 to 2012 is due to a weakening of circulation.
The weakening of dynamic component is also evidenced by the vertical pressure velocity of 500 hPa. All five reanalysis datasets show a weakening trend from 1980 to 2012 at the 1% significance level (Fig. 7d). The weakening of ascending motion is consistent with the decreasing trend of SPR precipitation. Then we further examine the large-scale atmospheric circulation associated with the weakening of ascending motion (Fig. 7e). The anomalous descending motions are over southeastern China and surrounding regions, the long-term linear trend of ω500 is around 0.03 Pa s −1 (33 year) −1 over these regions, corresponding with the low-level northwesterly anomalies over eastern China and South China Sea. There are studies (Li et al. 2016) suggested that the intensified latent heating released by the convection over South China Sea and the Philippine Sea could build and reinforce the heat-induced atmospheric circulation, suppress the ascending motion over southeastern China, thus leads to the decreasing trend of the SPR.

Interannual variability of the SPR
Regarding interannual variability, standard deviations of the SPR over southeastern China are shown in Fig. 5c, d. In observation, the spatial distribution of the standard deviation is along the coastline of southeastern China and the center is located at the Pearl River Estuary (approximately 22.5°N, 113.7°E). The maximum value in the center of standard deviations is over 1.8 mm day −1 in observation (Fig. 5c). The pattern of the ensemble mean is similar to observation, but with a weaker intensity (Fig. 5d).
ENSO is an important driving factor for the interannual variability of the SPR over southeastern China Feng and Li 2011;Li et al. 2013). The influences of El Niño and La Niña are asymmetric, and only the positive or warm events have significant influences on the SPR (Feng and Li 2011). To investigate the influence of El Niño on the SPR, we select 4 El Niño years (i.e., 1983, 1987, 1992 and 1998; only in these 4 years the Niño3.4 index was higher than +1.0 °C in preceding winter and could still reach up to +1.0 °C during the decaying spring) based on the Niño3.4 index (Trenberth 1997) to make a composite analysis. Per observation, the spatial distributions of composite precipitation anomalies (Fig. 8b) are consistent with the standard deviation (Fig. 5c, PCC is up to 0.87). The PCCs are 0.67, 0.59, 0.70, 0.53 and 0.19 for JRA25, JRA55, ERAIM, ERA20C and NCEP2, respectively ( Table 2). The ensemble mean of the five reanalysis datasets captures the precipitation anomalies with a PCC value of 0.64 (Fig. 8c).
To understand which processes mainly contribute to positive precipitation anomalies, we examine each component of the moisture budget anomalies in an El Niño decaying spring (Fig. 8a). The PCC between the anomalies of vertical moisture advection and precipitation is 0.77, so the positive precipitation anomalies in the El Niño decaying spring are mainly caused by the positive vertical moisture advection anomalies (Fig. 8d). The changes in horizontal moisture advection are relatively small and show negative anomalies. The contribution of evaporation is negligible. We further examine the thermodynamic and dynamic components of vertical moisture advection. Both have positive contributions, but the dynamic component dominates (Fig. 8a). Thus, the positive precipitation anomalies in the El Niño decaying spring are mainly caused by the enhancement of upward motion.
What is the large-scale circulation pattern that controls the enhanced upward motion? The vertically integrated water vapor transport flux and its convergence in the El Niño decaying spring are shown in Fig. 9. The anomalous lower-tropospheric anticyclone over the western North Pacific (i.e., WNPAC) establishes and matures during El Niño years, resulting from a Rossby-wave response to suppressed convective heating, coupled with in situ ocean surface cooling and the subsidence forced remotely by the central Pacific warming (Wang et al. 2000;Wang and Zhang 2002). The convergence of water vapor transport flux in JRA55 has an evident center around the Pearl River Estuary, consistent with the spatial distribution of positive precipitation anomalies (Fig. 9b). Although the five reanalysis datasets exhibit a spread in simulating the convergence of water vapor transport flux, they are highly consistent in reproducing the WNPAC in the El Niño decaying spring (Fig. 9). All five reanalysis datasets feature a strengthened southwestern water vapor transport at the western flank of the WNPAC, which brings more moisture from the tropical Pacific and South China Sea and converges over southeastern China, and then leads to more precipitation in these regions.

Summary
The SPR over southeastern China is a unique climatic phenomenon before onset of the East Asian summer monsoon.
In this study, we analyze the climatology, the long-term trend and interannual variability of the SPR derived from five reanalysis datasets compared to GPCP data. The moisture budget is diagnosed to better understand the formation and changes of the SPR. The major results are summarized in Fig. 10. The main findings are summarized as follows: 1. Climatology (see Fig. 10a  tion. JRA55 has a rainfall intensity (5.22 mm day −1 ) closest to observation (5.25 mm day −1 ) among the five reanalysis datasets [4.94 (4.54-5.32) mm day −1 ]. All five reanalysis datasets exhibit a strong westerly wind flow over the southern Tibetan Plateau and southwest wind flow over the western North Pacific that are the two main channels remotely supplying water vapor for the SPR. Locally, the moisture budget diagnosis shows that the SPR is primarily contributed by evaporation [2.98 (2.69-3.67) mm day −1 ] and vertical moisture advection [2.54 (1.83-2.79) mm day −1 ]. Area-averaged horizontal moisture advection and residual component contribute little to the total SPR precipitation. NCEP2 overestimates evaporation (3.67 mm day −1 ) and ERA20C underestimates vertical moisture advec-tion (1.83 mm day −1 ) compared with other reanalysis datasets. 2. The long-term trend (see Fig. 10b

Discussion
How to understand the difference of precipitations derived from five reanalysis data? The five reanalysis datasets differ in atmospheric forecast model, model resolution, convective parameterization schemes, data assimilation algorithm and satellite data processing (Table 1), which could cause different results derived from different reanalysis (Bosilovich et al. 2008;Lin et al. 2014). For example, the data assimilation methodology of ERAIM and JRA55 is fourdimensional variational analysis (4D-Var), which is the state-of-the-art data assimilation method. What's more, the model resolution of ERAIM (T255L60) and JRA55 (TL319L60) are relatively higher than other three reanalysis datasets (T106L40 for JRA25, T159L91 for ERA20C and T62L28 for NCEP2). In our results, both ERAIM and JRA55 have relatively higher skill in reproducing the climatology and interannual variability of the SPR, indicating that the model resolution and data assimilation algorithm could influence the performance of reanalysis. In addition, the ERAIM and JRA55 also assimilate the radiance products in satellite data. The ERA20C uses the same data assimilation method as ERAIM and uses a relatively higher vertical resolution, but it just assimilated observational surface pressures, mean sea level pressures and surface marine winds. That's the reason why the performance of ERAIM is better than ERA20C in reproducing the climatology and interannual variability of the SPR. What are the mechanisms for the decreasing trend of the SPR? Although this is beyond the scope of this study, we note that several mechanisms have been suggested. The dynamics of spring rainfall reduction in recent decades over southeastern China is associated with anomalous atmospheric circulation induced by the La Niña-like pattern with warming in the tropical Indian Ocean and western Pacific region (Qiu et al. 2009;Zhu et al. 2014). In addition, the reduction of late spring precipitation concurs with the upper-level cooling, which is associated with an anomalous meridional cell over southeastern China and is found to strongly link to the North Atlantic Oscillation (Xin et al. 2006(Xin et al. , 2008. The aerosol forcing Moisture budget components from 1980 to 2012 over the major rainfall center for the five reanalysis datasets: a Moisture budget components for climatology, the bars indicate precipitation, evaporation, vertical moisture advection, horizontal moisture advection and the residual term in the ensemble mean of five reanalysis datasets (unit: mm day −1 ), different symbols represent results of each reanalysis dataset, the red dash line indicates climatological SPR precipitation in GPCP. b The same as a, but the bars indicate the longterm trend of precipitation, evaporation, vertical moisture advection, the dynamic term, the thermodynamic term and the horizontal moisture advection term [unit: mm day −1 (10 year) −1 ], the red dash line indicates the long-term trend of the SPR in GPCP. c The same as b, but for composite results in the El Niño decaying springs is also a potential mechanism. The anthropogenic aerosols could make a shortened duration and reduced rainfall amount to the SPR (Deng et al. 2014). Modeling study of the effect of anthropogenic aerosols shows that the forcing of anthropogenic aerosols could induce the descending motion over southern China and lead the decreasing trend of the SPR (Kim et al. 2007;Hu and Liu 2013). Although the mechanisms dominating the weakening trend of SPR deserves further study, the current analysis presents an observational metric for examining the suggested mechanisms and may serve as a basis for model simulation of spring drought over southeastern China.
Why do only the positive or warm events of ENSO have significant influences on spring rainfall? There exists an asymmetry of the western North Pacific (WNP) low-level atmospheric circulation anomalies between the El Niño and La Niña during the cold seasons from winter to spring. Due to the longitudinal shifting of anomalous heating associated with El Niño and La Niña and the amplitude asymmetry of the associated sea surface temperature anomalies (SSTAs) in the WNP, the anomalous WNP cyclone center during La Niña shifts westward relative to an anomalous WNP anticyclone (WNPAC) center during El Niño (Wu et al. 2010). Following the asymmetric response of WNP low level atmospheric circulations, the enhanced precipitation over the southeastern China during El Niño years are more significant (Feng and Li 2011;Guo et al. 2017). Hence in our study, we focus on the influence of El Niño on the SPR.
Finally, there is an increasing trend of available reanalysis datasets. How to reliably use these reanalysis data sets has been of great concern to the climate research community. Based on our evaluation, among the five sets of reanalysis data we evaluated, both the JRA55 and ERAIM show higher skills in reproducing the climatology and changes of the SPR in comparison to the observations, and are thus recommenced for climate variability and climate change researches with focus on East Asia.