Developing and Demonstrating Climate Indicators for Monitoring the Changing Water Cycle

In this study, we developed a suite of spatially and temporally scalable Water Cycle Indicators (WCI) to examine the long-term changes inwater cycle variability and demonstrated their use over the contiguousUS (CONUS) during 1979–2013 using theMERRA reanalysis product. TheWCI indicators consist of six water balance variables monitoring the mean conditions and extreme aspects of the changing water cycle. The variables include precipitation (P), evaporation (E), runoff (R), terrestrial water storage (dS/dt), moisture convergence flux (C), and atmospheric moisture content (dW/dt). Means are determined as the daily total value, while extremes include wet and dry extremes, defined as the upper and lower 10th percentile of daily distribution. Trends are assessed for annual and seasonal indicators at several different spatial scales. Our results indicate that significant changes have occurred in most of the indicators, and these changes are geographically and seasonally dependent. There are more upward trends than downward trends in all eighteen annual indicators averaged over the CONUS. The spatial correlations between the annual trends in means and extremes are statistically significant across the country and are stronger for P, E, R, and C compared to dS/dt and dW/dt.


Introduction
The Intergovernmental Panel on Climate Change (IPCC) stated in its Fifth Assessment Report (AR5) that the global mean surface temperature has increased by about 0.85 ∘ C over the twentieth century and this increase is likely due to the enhancement of the greenhouse effect [1].There is growing evidence that anthropogenic warming would induce an intensification of the global hydrological cycle that produces increased atmospheric moisture content, faster evaporation, more precipitation and runoff, leading to exacerbation of hydrologic extremes of flooding and droughts, and redistribution of precipitation from the midlatitudes to high and low latitudes [2][3][4][5][6][7][8][9].Because of the central importance of water in the Earth system, any changes in mean, variability, and extremes of water cycle would have profound impact on our society and environment.It is essential to closely monitor water cycle changes and improve our knowledge of how the water cycle responds to climate change.
Tools for monitoring climate variability and change have been developed primarily through the use of indicators.
Indicators usually are measurements or calculations that represent the important features of the status and trend of climate system and are often used to guide both research and decisions [10].There are several commonly used climate change indictors across the US and the globe.Karl et al. [11] introduced a combined Climate Extremes Index (CEI) to quantify observed changes in multiple types of climate extremes within the contiguous US (CONUS).The CEI is composed of five components measuring the fraction area of the CONUS that experiences extremes outside the 90th/10th percentile in the maximum and minimum temperature, extreme 1-day precipitation, wet/dry days, and the Palmer Drought Severity Index (PDSI).A revised version of the CEI from Gleason et al. [12] is calculated operationally for different season and year by the US National Ocean and Atmospheric Administration (NOAA).The US Drought Monitor (USDM) is one of the widely used drought indices that incorporates multiple data sources and derivative products from local to national scales [13].The USDM utilizes a unique classification system from D0-D4 representing 30th, 20th, 10th, 5th, and 2nd percentile to measure the strength 2 Advances in Meteorology of drought.The Expert Team on Climate Change Detection and Indicators (ETCCDI) was established as a joint venture of the World Meteorological Organization Commission for Climatology (CCl) and the World Climate Research Program Climate Variability and Predictability (CLIVAR) program to develop climate change indicators focusing on extreme events [14].In total, 27 indicators are defined to provide a comprehensive overview of the changes in temperature and precipitation extremes.The ETCCD identifies most of the indicators by means of absolute, threshold, percentile, and duration values based on daily climate data.
Most studies, however, have addressed only the extreme aspect of water cycle with main focus on precipitation related climate change.To our knowledge, very few studies have been directed toward conducting a comprehensive analysis of changes in many aspects of water cycle in a warming climate.In this work, we propose a set of indicators to identify warming induced changes to several key water balance variables across the CONUS using a global reanalysis dataset.They are obtained at a grid point, a region, and the entire globe over annual, seasonal, and long-term time scales; that is, they are spatially and temporally scalable.This design flexibility is valuable for monitoring regional, continental, and global water cycle variability and hence it is important for decision-making associated with climate change.Six water balance variables are selected in this study, precipitation, evaporation, runoff, moisture convergence flux, atmospheric moisture content, and terrestrial water storage.To provide a general picture of water cycle changes, we detect trend in the indicators of these water variables and determine if existing trends are statistically significant.
The remaining paper is organized as follows.The methodology is briefly described in Section 2. The results are presented in Section 3, followed by the conclusions and discussions in Section 4.

Methodology
2.1.Data.Reanalyses have the advantage of generating complete global coverage over long time periods, thus providing a potential source for studying long-term changes in water cycle especially for regions where the station data are spare or unavailable.Reanalysis data have recently been used in assessing the trends in climate extreme indicators in conjunction with station and model data [15][16][17].We are thus motivated to utilize the reanalysis product to monitor the changes in the different aspects of the water cycle.
The water cycle related variables are obtained from the Modern Era Retrospective-Analysis for Research and Applications (MERRA [18]), a state-of-the-art reanalysis that provides the multiyear estimates of atmospheric conditions and land surface fields.MERRA reanalysis is generated by the NASA Global Modeling and Assimilation Office (GMAO) using the Goddard Earth Observing System version 5 (GEOS-5) atmospheric GCM through a threedimensional variational (3DVAR) system.Besides in situ observations, MERRA also assimilates satellite observations, including Special Sensor Microwave Imager (SSM/I) radiances, TIROS Operational Vertical Sounder (TOVS) radiances, Atmospheric Infrared Sounder (AIRS) radiances, and scatterometer-based wind retrievals.MERRA uses Catchment Land Surface Model to improve prediction of the land surface processes through explicit modeling of subgrid-scale soil moisture variability and its effect on runoff and evaporation.The soil and vegetation parameters are obtained from the NASA GEOS-5 GCM.The MERRA data are available for the entire globe at 1/2 ∘ ×2/3 ∘ spatial resolution with 72 hybrid vertical levels from 1979 to present.
MERRA has the unique Incremental Analysis Update (IAU) procedure, which could reduce the precipitation spindown problem and improve the stratospheric circulation.But the water vapor increment in MERRA is found to be affected by the introduction of Advanced Microwave Sounding Unit (AMSU-A) series in late 1998, which directly influences precipitation [19].

Definition of Water Cycle
Indicators.The goal of this paper is to develop a suite of climate indicators to quantify the changes of water cycle in the context of warming climate.IPCC AR5 documented the observed twentieth century changes in the mean states of the different aspects of water cycle globally and regionally [1].Their impact on society is felt through changes in variability and intensity.In addition to the mean, frequency, intensity, and duration of extreme events have also experienced significant changes as manifested by their adverse effect on infrastructure and community.Therefore, the Water Cycle Indicators (WCI) are developed in the study primarily for monitoring the mean conditions and extreme aspects of the changing water cycle.Specifically, the WCI indictors contain six water cycle components involved in the terrestrial and atmospheric water balance equations: where  is terrestrial water storage term,  is precipitation,  is evaporation,  is runoff,  is atmospheric moisture content, and  is the vertically integrated moisture convergence flux.The expressions of  and  are given by where  is constant,   is surface pressure,  is specific humidity, and V is wind speed.For each of six water balance variables in (1), we define three types of indicators: mean, dry and wet extremes.The mean indicator is calculated as a normalized anomaly: where   denotes the lower or upper 10th percentile value at time  and  and   indicate the mean and standard deviation for corresponding dry or wet extremes.The complete set of 18 WCI indicators is listed in Table 1.The defined indicators are aimed for the assessment of changes in many aspects of the water cycle including means and extremes.They are also helpful for a better understanding of climate change and gaining insight into the impact of changes on human activities and the environment.Additionally, these indicators could act as the reference data to evaluate future climate projections and hence to inform potential consequences of climate change to the public and policy makers.It is worthy to mention that this suite of indicators is scalable in space and time.They can be used to assess the broad scale changes in the water cycle.We mostly focus on our discussion to assessment of annual and seasonal aspects of regional change over the CONUS, while we must point out that the WCI indicators are also capable of global and continental analysis.The mean and extremes indicators are calculated from daily MERRA variables on annual and seasonal basis.For threshold indicators, the values of upper or lower 10th percentile are computed from daily time series over each calendar year or individual season during 1979-2013.

Estimate of Trend.
To determine whether the trends in WCI indicators are significant, the nonparametric Mann-Kendall (MK) [20,21] test is adapted and applied.The MK test is widely used in assessing the significance of monotonic trends in hydrometeorological time series [22][23][24].This technique is a rank-based procedure, robust to the influence of outliers.It does not make any assumptions on the distribution of the data but requires that individual data values are uncorrelated.This condition is met by taking into account serial correlation in the residuals when conducting the statistical significance test.Under the null hypothesis  0 that sample data,   ( = 1, 2, . . ., ), are independent and identically distributed random variables, the MK statistic  is defined: where sgn assigns 1, −0, or −1 depending on whether   −   is positive, 0, or negative.The standardized statistic  is approximated by the standard normal distribution with a mean of 0 and variance of 1: where the variance   is given by where  is the number of tied groups and   is the number of data points in th tied group.If || >  1−/2 , the null hypothesis  0 is rejected at significance level .A positive (negative) value of  indicates an upward (downward) monotone trend for the data   .The magnitude of the trend is calculated as the slope  of the Kendal-Theil robust line [25]: where  is time and the slope  is expressed as for  = 2, 3, . . ., .The intercept  is solved by

Probability Density Function (PDF)
. In order to determine if there was a temporal change in the distributions of the above indicators over the study period, we compare the probability distributions of given indicator over two independent 17-year periods, 1979-1995 and 1996-2012.For each WCI, PDF was obtained by binning annual (or seasonal) values and calculating the corresponding probabilities.The sum of area below PDF curve across all bins is normalized to 1.We examine if the distribution for the period 1979-1995 differs from the distribution for 1996-2012.The nonparametric Kolmogorov-Smirnoff test is used to test whether the probability distributions of one specific indicator from two time periods are significantly different or not.

Results
One of the appealing features of the WCI indicators is the ability to analyze long-term changes in the water cycle at various temporal and spatial scales.In this section, we present spatial pattern of changes in WCI indicators and temporal variation of spatially averaged indicators to assess the annual and seasonal characteristics of water cycle changes over the CONUS.2. The formula and red line represent the Kendall-Theil robust line in each panel.Nationally averaged, mean  has slightly decreased by ∼0.008/year (bot not significant) since 1979.Although several studies found a nationwide increase of annual total precipitation [26][27][28], they are based on different station data spanning a variety of time periods.

Trends in Annual Indicators
Dry extreme of  is an indicator of extremely dry condition and exhibits a tendency toward drier conditions over southwest and northeast regions and conversely wetter conditions over the remaining country.However, these changes are small and mostly are not statistically significant.This general pattern is also reflected in the national time series (Figure 2) showing a nonsignificant decrease from 1979 to 2013.For  wet extreme, a potential flood indicator, the trend patterns are remarkably similar to the changes of mean .It can be seen that 19.9% of the nation has experienced significant increase of intense  concentrated over the eastern and southern coastal regions as opposed to 24.2% with significant decrease predominantly occurred over northwest and central regions between 1979 and 2013.Similar to mean , the nationally averaged  wet extreme shows a decline over the past several decades with steeper slope than mean and dry extreme indicators.Earlier observational studies of trends   in extreme indicators have revealed significant changes in extreme precipitation across the country in the regional [26][27][28] and global [14,29,30] analysis.However, their definitions and thresholds of extremes are quite different from this study and thus make it difficult to directly compare our results with previous findings.
To further determine the relation between the mean and the extremes, we also calculated the spatial correlations between the annual trends in mean and extremes for each indicator over the CONUS (Table 3) using Pearson productmoment correlation.It is noted that the correlation coefficients for mean  with respect to dry and wet extremes are 0.39 and 0.93, respectively, which is consistent with their trends changes in Figures 1 and 2. The high correlation between the mean and wet extreme of precipitation is also in agreement with the findings of [26], who showed that most of the increase in annual precipitation is associated with more intense and frequent precipitation over the CONUS during the last sixty years.
The  indicators exhibit similar patterns of changes to the trends of mean and wet extreme of  over the 35-year period.A vast area of the Midwest also exhibits significant increase in dry extreme, leading to a larger percentage of the nation showing positive trend in dry extreme (35.1%) than mean and wet extreme (22.6% and 21%).Furthermore, each averaged  indicator has stronger decreasing trend than the corresponding  indicator with mean  changing at the fastest rate (Figure 2).For  indicators, the correlation is higher between mean and wet extreme than dry extreme (Table 3).Such relationship is in agreement with  indicators but the correlation between mean and dry extreme is bigger in  than  due to more positive changes occurring in dry extreme of .
Spatial distribution and correlation of trends in  indicators agree well with the  and  indicators (Figure 1).However, trends of  have smaller magnitude and lower level of statistical significance.This difference is also seen in Table 2, which shows smaller percentage of the nation identified as having significant trends in comparison with  and  indicators.In addition, the nationally averaged mean and wet extreme of  show increasing trend in contrast to the downward trend in both  and  indicators.This is probably because more precipitation turns into runoff as  shows smaller decreasing rate than .The significant increase of mean and dry extreme of  has been observed over parts of central and southwest regions.These patterns are apparently opposite to the most of the indicators in above-mentioned three water balance variables.This is somewhat expected according to water balance equations (1).The wet extreme of  shows a widespread and significant increase over much of the country, leading to more significant positive changes than significant negative change compared to the other two  indicators (Table 2).Mean  is highly correlated with each of extreme indices but dry extreme has a larger coefficient than wet extreme (0.83∼0.64).Such relationship is also noted in the national series (Figure 2), which shows decreasing trend in mean and dry extreme indicators and increasing trend in wet extreme.
The trends in two water storage variables, / and /, differ considerably from the other indicators.Most notably the mean trends appear spatially more heterogeneous and statistically not significant, resulting in fewer areas of significant changes (Table 2) and smaller national trends (Figure 2).In addition, dry extreme displays opposite changes with wet extreme for both / and /, as evident from their spatial patterns, significance percentage, and national series.For these two water storage indicators, the mean and extremes have weaker relations compared to the other four water variables as evidenced by smaller and even nonsignificant correlation coefficients, for instance, mean and dry extreme of /.
Given the trends patterns identified by each of WCI indicators, we are interested in looking at the general picture of changes in water cycle variability.Figure 3 shows the sum of signs of annual trends from all 18 indicators across the county.Regions with positive sums are located along a large area over eastern coast, parts of north, northwest coast, and west with the maximum value exceeding 12.In contrast, there are widespread and coherent negative values appearing over much of central, southwest, and large portion of northwest, where the maximum sum lies above 12.This contrast pattern is pretty much consistent with annual trends in most of the indicators and further reinforces that water cycle displays varied geographical features associated with climate change.
To further illustrate the spatial scalability of the WCI indicators, we show in Figure 4 the time variations of indicators for water balance variables at two different spatial domains, a small region over the east coast (RG1 hereafter) and one grid in the central region (RG2 hereafter), denoted by red box and red dot in Figure 3, respectively.The region RG1 locates over the eastern region where positive changes of many WCI indicators have been observed, while RG2 locates in the swath of negative changes.Not surprisingly, there is significant increase in most of the indicators over the RG1 with several mean and dry extreme indicators showing decline over the 35-year period.There is much more significant decreasing trend than significant increasing trend in the time series of indicators over RG2.

Trends in Seasonal Indicators.
Besides spatial scalability, the WCI indicators are also scalable on different temporal regimes.In this section, we examine the trends in seasonal WCI indicators associated with climate change over the CONUS.Seasons are defined as follows: winter (December-January-February), spring (March-April-May), summer (June-July-August), and fall (September-October-November).The spatial patterns of trends (Figures 5 and 6) and time series (Figure 7) are presented for all indicators during summer and fall over the period of 1979-2013.Seasonal trends in all indicators are generally consistent with annual results but differing in magnitude and significance.We found that the spatial trends in  and  indicators during fall tend to agree well with the corresponding annual patterns, suggesting that the changes observed in annual indicators mainly are attributed by fall and, to a less extent, summer.For , , and / indicators, annual trends are primarily affected more by summer changes, while spring (not shown) contributes to the annual change of /.The trends in dry extremes of  and  are mostly nonsignificant in all seasons with less than 9% of the nation showing significant changes in each season (Table 2).In general, the most significant changes occur in summer with the least significant changes appearing in fall.In addition, more areas exhibit positive changes for , , and / indicators during winter (not shown) and for  indicators during fall.
Generally, the national series show a mixed of upward and downward trends across the seasons although these trends are mostly nonsignificant (Figure 7).Overall, there are more positive changes than negative changes in winter and fall where at least two upward trends are observed in each of the , , , and / indicators.Notably, fall has the most positive trends in the mean and dry extreme indicators, while the most increased wet extreme events have occurred in winter.Large number of decreasing trends are found in spring and summer with four to five negative trends for each of mean and extremes indicators.
Throughout all seasons, trends in wet extremes of , , , and / are highly correlated with the means than the corresponding dry extreme (Table 3), which is opposite to  indicators.The strongest relation occurs in  indicators across all seasons with coefficients constantly above 0.7.The correlation coefficients for / are considerably small and nonsignificant except in fall, indicating the weakest relation between the mean conditions and the extremes.
Like annual analysis, we also combine the signs of trends of all 18 indicators for all seasons as shown in Figure 8. Overall, seasonal sums exhibit similar patterns to those observed in annual results.It is evident that summer generally appears to bear the most similarity to annual pattern with discernible differences in winter where parts of central and southwest show larger positive values.The results are in general consistent with the pattern of trends from individual index and further highlight the seasonal aspects of water cycle changes.

Decadal Changes in the PDFs.
To further investigate the long-term changes in water cycle variability, we present the PDFs for annual WCI indicators for the periods 1979-1995 and 1996-2012 in Figure 9.For each index, time series from across the country are used to generate PDFs for each of the two periods.We found that the distributions of all the indicators from different time periods are significantly different.The figure indicates a small shift toward drying conditions associated with a shift toward decreased extremes in the  indicators between 1996 and 2012 when compared to 1979-1995.There are also visible shifts in  corresponding to decreased mean and wet extreme along with increased dry extreme.The distributions of the  indicators show a marked reduction in the mean and extremes during the most recent decade.For /, the PDFs exhibit little change in the mean accompanied with slight increase in dry extreme and decrease in wet extreme.Mean and dry extreme of  and / have significantly decreased over the recent 17-year period but the corresponding wet extremes have increased.
For all the indicators, we also calculated the PDFs in all the seasons over two different periods.During each season, the PDFs for the period 1979-1995 are statistically significant compared to those for the 1996-2012 period for every index.The changes in general are more pronounced in summer (Figure 10) with least changes generally occurring in fall (Figure 11).The  PDFs for summer show a similar pattern of change to annual values although of larger magnitude.Means of  in spring (not shown) and fall also exhibit a tendency toward drying conditions that is consistent with the annual change, whereas winter mean  shows a shift toward wetter

Conclusions and Discussions
This study presented a suite of spatially and temporally scalable climate indictors to examine the long-term changes in water cycle variability across the CONUS over the period 1979-2013 using MERRA.The indicators consist of six water balance variables monitoring the mean conditions and extreme aspects of the following: precipitation, evaporation, runoff, vertically integrated moisture convergence flux, and terrestrial and atmospheric water storage changes.The means  are defined as daily total value, while extremes are determined as upper (wet) and lower (dry) tails of daily distribution based on the 10th percentile value.
Trends and statistical significance are estimated for annual and seasonal indicators across different spatial scales.The annual analysis found statistically significant increase along northeast, southeast, west north central, and northwest coastal areas of the CONUS but significant decrease over the remaining country in most of , , and  indicators during the 35-year period.The significant increase in  indicators was detected over much of central and southwest with a large portion of the eastern coast regions showing more positive changes for wet extreme.Throughout the country, these four indicators generally have more significant negative changes than significant positive changes, which is also reflected in their averaged time series showing decreasing trends except in  and .The mean / and / indicators do not show significant changes although the patterns are mixed, consistent with their smaller national trends.The opposite changes have been observed in extremes for both / and / as evident from the corresponding spatial pattern and temporal variability.The results also reveal that the national trends in means are greater than trends in dry extremes except in  and / and smaller than trends in wet extremes except in  although most of these trends are not statistically significant.
In general, there are statistically significant correlations between the trends in means and extremes across the country.The spatial correlations between trends in extremes and trends in means are stronger for , , and  than for / and /.Trends in , , and  means are highly correlated with the corresponding wet extremes than the corresponding dry extremes, while mean  has higher correlation with dry extreme than wet extreme.The strong relationship between the trends in means and extremes implies that the same driver is responsible for the long-term changes of means and extremes.
A partition of entire dataset between 1979 and 2013 into two equal time periods suggests that the probability distributions of all eighteen indicators are significantly different between 1979-1995 and 1996-2012 over the nation.The  indicators exhibit a tendency toward drying conditions in mean associated with shift toward less intense extremes in the most decade.The distributions for  indicators show shift towards decreased mean and wet extreme along with increased dry extreme.There is also visible reduction in the mean and extremes of  during 1996-2012 compared to 1979-1995.For /, , and /, apparent shifts are observed in the distributions of mean and extremes.
Seasonal trends show similar patterns to those observed in annual results but with different magnitude and statistical significance.For instance, annual trends of  and  mostly   except for / which has smaller and even nonsignificant correlation coefficients.Summer has the most pronounced changes in the probability distributions while the least changes appear in fall.

Advances in Meteorology
It is worthy to mention that the changes in the water cycle could pose a significant challenge to societies and ecosystems.For instance, as extreme precipitation increases, both dry and wet conditions can intensify, leading to prolonged drought or flash flooding.Drought is one of the most costly natural hazards affecting agriculture crops, vegetation cover, and water supplies.On the other hand, heavy rainfall increases runoff and can destroy infrastructure, contaminate water bodies, and threaten public health.

Advances in Meteorology
The indicators proposed in this study are valuable in monitoring regional and global changes in water cycle variability.However, the results from reanalysis data should be interpreted with caution.Reanalyses are essentially model output constrained by available observations.Variables such as precipitation are usually not directly assimilated to the forecast model.Furthermore, reanalyses are subject to inhomogeneities due to changes in the satellite observing system.These issues may affect the estimated trends in mean and extreme indicators presented here.
Our results cannot be directly compared with the previous trend analysis [26][27][28] as they use station data with different time periods based on different methodology.On the other hand, the station data only represent the value of a point location, while the reanalysis data are representative of the area average of a grid.In addition, in situ observations are often affected by the changes in observation time, instrumentation, and station location.Therefore, there is need to develop the effective approaches for comparing the trends derived from station and reanalysis data.Such comparison is useful for improving the representation of relevant processes in reanalysis products.

Figure 1 : 3 Figure 2 :
Figure 1: Trends (per year) in the annual WCI indicators of mean, dry extreme, and wet extreme of , , , /, , and / for 1979-2013 in MERRA.Hatching indicates trends that are statistically significant at the 5% level.

Figure 3 :
Figure 3: The sum of signs of annual trend changes of all 18 WCI indicators over the CONUS.

Figure 5 :
Figure 5: Same as Figure 1 but for summer.

Figure 6 :
Figure 6: Same as Figure 1 but for fall.

Figure 7 :
Figure 7: Same as Figure 2 but for winter, spring, summer, and fall, indicated by purple, green, red, and blue.The Kendall-Theil robust line is denoted by the same light color of the corresponding season.

Figure 8 :
Figure 8: Same as Figure 3 but for winter, spring, summer, and fall.

Figure 10 :
Figure 10: Same as Figure 9 but for summer.

Figure 11 :
Figure 11: Same as Figure 9 but for fall.

Table 1 :
A suite of 18 indicators defined for mean, dry and wet extremes of , , , /, , and /.

Table 2 :
Percentage of significant upward (indicated by ↑) and significant downward (indicated by ↓) trends over the CONUS for indicators of , , , /, , and /.

Table 2 ,
. Figure1shows the spatial maps of trend in mean and extremes indicators for six water balance components over the CONUS for the period 1979-2013 derived from MERRA.Regions where trend is statistically significant at the 5% level are denoted by hatching.20.2% of the CONUS shows significant upward trend in mean  while 29.6% of the nation shows significant downward trend.The temporal evolution of each indicator averaged over the nation during 1979-2013 period is shown in Figure [27]annual mean  shows a significant tendency toward wetter conditions over a large part of east coast and west north central areas, whereas significant drying trends are located over central and west coast during the past 35 years.Previous observational studies by Kunkel et al.[26]reported local increase of annual  in the central and southwest regions over 1931-1990 period, while Groisman et al.[27]found the increase in the central region with decrease over the southwestern area during 1900-2002.Since these studies used different time periods for trend analysis, their results cannot be directly compared with our results.According to

Table 3 :
Spatial correlation between the trend in mean and extremes indicators of , , , /, , and / over the CONUS.The coefficients are highlighted with asterisk mark if they are statistically nonsignificant at 5% level. 6 resemble the changes in fall, while annual patterns of , , and / are more consistent with summer changes.The most significant changes occur in summer whereas the least significant changes tend to occur in fall.Overall, there are more positive national trends than negative national trends in winter and fall where each of , , , and / indicators exhibits at least two upward trends.Seasonally, trends in means are highly correlated with trends in extremes