Recent changes of relative humidity: regional connections with land and ocean processes

We analyzed changes in surface relative humidity (RH) at the global scale from 1979 to 2014 using both observations and the ERA-Interim dataset. We compared the variability and trends in RH with those of land evapotranspiration and ocean evaporation in moisture source areas across a range of selected regions worldwide. The sources of moisture for each particular region were identified by integrating different observational data and model outputs into a Lagrangian approach. The aim was to account for the possible role of changes in air temperature over land, in comparison to sea surface temperature (SST), but also the role of land evapotranspiration and the ocean evaporation on RH variability. The results demonstrate that the patterns of the observed trends in RH at the global scale cannot be linked to a particular individual physical mechanism. Our results also stress that the different hypotheses that may explain the decrease in RH under a global warming scenario could act together to explain recent RH trends. Albeit with uncertainty in establishing a direct causality between RH trends and the different empirical moisture sources, we found that the observed decrease in RH in some regions can be linked to lower water supply from land evapotranspiration. In contrast, the empirical relationships also suggest that RH trends in other target regions are mainly explained by the dynamic and thermodynamic mechanisms related to the moisture supply from the oceanic source regions. Overall, while this work gives insights into the connections between RH trends and oceanic and continental processes at the global scale, further investigation is still desired to assess the contribution of both dynamic and thermodynamic factors to the evolution of RH over continental regions.


Introduction
Relative humidity (RH) is a key meteorological parameter that determines the aerodynamic component of the atmospheric evaporative demand (AED) (Wang and Dickinson, 2012;McVicar et al., 2012a).As such, changes in RH may significantly impact the evolution of the AED (Vicente-Serrano et al., 2014a), with particular implications for the intensity of the hydrological cycle (Sherwood, 2010), climate aridity (Sherwood and Fu, 2014), and severity of drought events (Rebetez et al., 2006;Marengo et al., 2008).
In a changing climate, temperature rise, as suggested by different climate scenarios, may impact the atmospheric humidity.Particularly, there is empirical evidence on the increase in the water vapor content at both the surface and upper tropospheric levels (Trenberth et al., 2005).In this context, numerous studies have supported the constant RH Published by Copernicus Publications on behalf of the European Geosciences Union.
scenario under global warming conditions (e.g., Dai, 2006;Lorenz and DeWeaver, 2007;Willett et al., 2008;McCarthy et al., 2009;Ferraro et al., 2015).In contrast, other studies supported the nonstationary behavior of RH, not only in continental areas located far from oceanic humidity (e.g., Pierce et al., 2013) but also in humid regions (e.g., Van Wijngaarden and Vincent, 2004).Assuming the stationary behavior of RH, the influence of RH on AED may be constrained, given that any possible change in AED would be mostly determined by changes in other aerodynamic variables (e.g., air temperature and wind speed) (McVicar et al., 2012a, b) or by changes in cloudiness and solar radiation (Roderick and Farquhar, 2002;Fan and Thomas, 2013).However, a range of studies have supported the nonstationary behavior of RH under global warming, giving insights into significant changes in RH over the past decades.A representative example is Simmons et al. (2010), who compared gridded observational and reanalysis RH data, suggesting a clear dominant negative trend in RH over the Northern Hemisphere since 2000.Also, based on a newly developed homogeneous gridded database that employed the most available stations from the telecommunication system of the WMO, Willett et al. (2014) found significant negative changes in RH, with strong spatial variability, at the global scale.This global pattern was also confirmed at the regional scale, but with different signs of change, including both negative (e.g., Vincent et al., 2007;Vicente-Serrano et al., 2014b, 2016;Zongxing et al., 2014) and positive trends (e.g., Shenbin et al., 2006;Jhajharia et al., 2009;Hosseinzadeh Talaee et al., 2012).
There are different hypotheses that explain the nonstationary evolution of RH under global warming conditions.One of these hypotheses is related to the slower warming of oceans in comparison to continental areas (Lambert and Chiang, 2007;Joshi et al., 2008).In particular, specific humidity of air advected from oceans to continents increases more slowly than saturation specific humidity increases over land (Rowell and Jones, 2006;Fasullo, 2010).This would decrease RH over continental areas, inducing an increase in AED and aridity conditions (Sherwood and Fu, 2014).Some studies employed global climate models (GCMs) to support this hypothesis under future warming conditions (e.g., Joshi et al., 2008;O'Gorman and Muller, 2010;Byrne and O'Gorman, 2013).
Another hypothesis to explain the nonstationary evolution of RH is associated with land-atmosphere feedback processes.Different studies indicated that atmospheric moisture and precipitation are strongly linked to moisture recycling in different regions of the world (e.g., Rodell et al., 2015).Thus, evapotranspiration may contribute largely to water vapor content and precipitation over land (Stohl and James, 2005;Bosilovich and Chern, 2006;Trenberth et al., 2007;Dirmeyer et al., 2009;van der Ent et al., 2010).Landatmospheric feedbacks may also have marked influence on atmospheric humidity (Seneviratne et al., 2006), given that soil drying can suppress evapotranspiration, reduce RH, and thus reinforce AED.All these processes would again reinforce soil drying (Seneviratne et al., 2002;Berg et al., 2016).
Indeed, it is very difficult to determine which hypothesis can provide an understanding of the observed RH trends at the global scale.Probably, the two hypotheses combined together can be responsible for the observed RH trends in some regions of the world (Rowell and Jones, 2006).In addition to the aforementioned hypotheses, some dynamic forces, which are associated with atmospheric circulation processes, can explain the nonstationary behavior of RH worldwide (e.g., Goessling and Reick, 2011).However, defining the relative importance of these physical processes in different world regions is quite challenging (Zhang et al., 2013;Laua and Kim, 2015).
The objective of this study is to compare the recent variability and trends in RH with changes in the two types of fluxes that affect RH: (i) vertical fluxes that were assessed using land evapotranspiration and precipitation and (ii) advection that was quantified using oceanic evaporation from moisture source areas.The novelty of this work stems from the notion that although different studies have already employed GCMs and different scenarios to explain the possible mechanisms behind RH changes under warming conditions, we introduce a new empirical approach that employs different observational datasets, reanalysis fields, and a Lagrangian-based approach, not only for identifying the continental and oceanic moisture areas for different target regions but also for exploring the relevance of the existing hypothesis to assess the magnitude, sign, and spatial patterns of RH trends in the past decades at the global scale.

Observational RH dataset
We employed the monthly RH HadISDH dataset, available through http://www.metoffice.gov.uk/hadobs/hadisdh(last access: 25 June 2018; Willet et al., 2014).We selected only those series with no more than 20 % of missing values over the period 1979-2014.In order to fill these gaps, we created a standardized regional series for each station using the most correlated series with each target series.In order to avoid biases, mostly originated from differences in the distribution parameters (mean and variance) between the candidate and the objective data series, a bias correction was applied to the candidate data.The data of the candidate series were rescaled to match the statistical distribution of the observed series to be filled, based on the overlapping period among them.Overall, a final dataset of 3462 complete stations spanning different regions worldwide and covering the period 1979-2014 was employed in this work.

Reanalysis RH dataset
Daily data of dew point (T d ), air temperature (T ), and surface pressure (P mst ) at a spatial interval of 0.5 • were obtained from the ERA-Interim, covering the period 1979-2014 (http://www.ecmwf.int/en/research/climate-reanalysis/era-interim, last access: 25 June 2018; Dee et al., 2011).To calculate RH, we followed the formulation used by Willett et al. (2014) for the HadISDH RH dataset.Based on the selected variables, we calculated the daily RH following Buck (1981): where e is the actual vapor pressure in hectopascals and e s is the saturated vapor pressure in hectopascals.As a function of the wet bulb air temperature (T w ) in degrees Celsius, e is estimated following two different equations with respect to water and ice.If T w is above 0 • C, e is calculated as where T d is the dew point temperature in degrees Celsius. where where P mst is the pressure at the height level.
T w is obtained according to Jensen et al. (1990): where and T is the 2 m air temperature in degrees Celsius e s is obtained by substituting T for T d .

Land precipitation and land air temperature
We employed the gridded land precipitation and surface air temperature data (TS v.3.23),provided by the Climate Research Unit (CRU; UK), at a 0.5 • spatial interval for the period 1979-2014 (Harris et al., 2014).

Ocean evaporation and continental evapotranspiration data
To quantify the temporal variability and trends in land evapotranspiration and oceanic evaporation, we employed two different datasets.First, the oceanic evaporation was quantified using the Objectively Analyzed air-sea Fluxes (OAFlux) product (Yu et al., 2008) from 1979 to 2014, which was used to analyze recent variability and changes in evaporation from global oceans (Yu, 2007).To account for land evapotranspiration, we employed the Global Land Evaporation Amsterdam Model (GLEAM) (version 3.0a) (http://www.gleam.eu/,last access: 25 June 2018; Miralles et al., 2011) from 1980 to 2014.This dataset has been widely validated using in situ measurements of surface soil moisture and evaporation across the globe (Martens et al., 2017).

Relative humidity trends
We assessed the seasonal (boreal cold season: October-March; boreal warm season: April-September) and annual trends in RH for 1979-2014 using two different global datasets (HadISDH and ERA-Interim).To quantify the magnitude of change in RH, we used a linear regression analysis between the series of time (independent variable) and RH series (dependent variable).To assess the statistical significance of the detectable changes, we applied the nonparametric Mann-Kendall statistic (Zhang et al., 2001) to account for any possible influence of serial autocorrelation on the robustness of the defined trends (Hamed and Rao, 1998;Yue and Wang, 2004).The statistical significance of the trend was tested at the 95 % confidence level (p < 0.05).
Following the trend analysis results, we selected those regions that showed a high agreement between HadISDH and ERA-Interim datasets in terms of the sign and magnitude of RH changes.Nevertheless, we also extended our selection to some other regions, with low station density in the HadISDH dataset.This decision was simply motivated by the consistent changes found over these regions, as suggested by the ERA-Interim dataset.For all the defined regions, we identified the oceanic and continental moisture sources by means of the FLEXPART Lagrangian model.

Identification of continental and oceanic moisture sources
We used the FLEXPART v9.0 particle dispersion model fed with the ERA-Interim reanalysis data.According to this model, the atmosphere is divided homogeneously into threedimensional finite elements (hereafter "particles"); each represents a fraction of the total atmospheric mass (Stohl and James, 2004).These particles may be advected backward or forward in time using three-dimensional wind taken from the ERA-Interim data every time step, with superimposed stochastic turbulent and convective motions.The rates of increase (e) and decrease (p) in moisture (e-p) along the trajectory of each particle were calculated via changes in the specific humidity (q) with time (e-p = mdq / dt), where m is the mass of the particle.A description of this methodology is detailed in Stohl and James (2004).
The FLEXPART dataset used in this study was provided by a global experiment in which the entire global atmosphere was divided into approximately 2.0 million particles.The tracks were computed using the ERA-Interim reanalysis data at 6 h intervals, at a 1 • horizontal resolution, and at a vertical resolution of 60 levels from 0.1 to 1000 hPa.For each particular target region, all the particles were tracked backward in time, and its position and specific humidity (q) were recorded every 6 h.All areas where the particles gained humidity (E-P > 0) along their trajectories towards the target region can be considered to be "sources of moisture".In contrast, all areas with lost humidity (E-P < 0) are considered to be "sinks".
We followed the methodology of Miralles et al. (2016), in which an optimal lifetime of vapor in the atmosphere was calculated to reproduce the sources of moisture and we defined the climatological spatial extent of each source region corresponding to a particular target region by applying a 95th percentile criterion computed for the annual and seasonal (boreal summer and winter) positive (E-P ) field (Vázquez et al., 2016).Then, for each year of the period, we estimated the total moisture supply from each source region.Also from FLEXPART simulations, we obtained the fractions of moisture from the continental and oceanic sources annually and for each cold and warm season.The purpose was to compare with the results obtained on the role of the land evapotranspiration and ocean evaporation in RH variability and trends.

Relationship between RH and the selected land-oceanic climate variables
Based on defining the spatial extent of each moisture source region, we calculated annual, warm, and cold season regional series for ocean evaporation and land evapotranspiration using the OAFlux and GLEAM datasets, respectively.The regional series of ocean evaporation and land evapotranspiration were created using a weighted average based on the seasonal and annual fields of (E-P ) > 0. This approach allows the creation of a time series that better represents the interannual variability in ocean evaporation and land evapotranspiration in the source(s) of moisture for each defined region.Following the same approach, we also calculated the regional series of SST corresponding to each oceanic moisture source region.Likewise, we calculated the regional series of land precipitation and air temperature for each target region using the CRU TS v.3.23 dataset, and the ratio between air temperature in the target region and SST in the source region.
For each target region, we related the regional series of seasonal and annual RH to the corresponding regional time series of all aforementioned climatic variables.However, to limit the possible influence of the trends presented in the data itself on the computed correlations, we de-trended the series of the climate variables prior to calculating the correlation.We also assessed changes in the regional series of the different variables; their statistical significance was tested by means of the modified Mann-Kendall test at the 95 % level.For each target region, we summarized the results of the magnitude of change in RH as well as other investigated variables at the seasonal and annual scales.However, to facilitate the comparison among the different variables and the target regions worldwide, we transformed the amount of change of each variable to percentages.
Finally, we also computed the association between RH and land evapotranspiration at the annual and seasonal scales using the available gridded evapotranspiration series.While a pixel-to-pixel comparison does not produce a reliable assessment of the possible contribution of land evapotranspiration to RH changes, given that the source of moisture can apparently be far from the target region, we still believe that this association can give insights into the possible relationship between land evapotranspiration and RH changes.

Trends in relative humidity
There are regional differences where the precipitation dominates (negative values) over the evaporation (positive values), from the ERA-Interim dataset (Supplement Fig. S1). Figure 1 summarizes the magnitude of change in RH for the boreal cold and warm seasons and at the annual scale for the period between 1979 and 2014.For HadISDH, it is noted that the available RH stations are unevenly distributed over the globe, with higher density in the midlatitudes of the Northern Hemisphere.Nevertheless, the available stations show coherent and homogeneous spatial patterns of RH changes (Fig. S2).The ERA-Interim dataset showed magnitudes of change close to those suggested by HadISDH.In addition, the ERA-Interim also provides information on RH changes in regions with low density of RH observations (e.g., eastern Amazonia, eastern Sahel, and Iran), suggesting a dominant RH decrease across these regions.For the boreal warm season, a clear tendency towards a reduction in RH was observed in vast regions of the world, including (mostly the Iberian Peninsula, France, Italy, Turkey, and Morocco) eastern Europe and the western part of Russia.This reduction was also noted in South America, with a general homogeneous pattern over Peru and Bolivia and a strong decrease over central Argentina.Conversely, the positive evolution of RH observed during the cold season across Canada and Scandinavia was reinforced during the boreal warm season.ERA-Interim also revealed a strong RH decrease over the whole Amazonian region and the western Sahel.A wide range of these regions exhibited statistically significant trends from 1979 to 2014 (Fig. S3).
With these complex spatial patterns of RH changes, there is a global dominant negative trend (Fig. 2).This pattern was observed using both the HadISDH and the ERA-Interim datasets, although there is marked spatial bias in data availability of the HadISDH.The relationships among the magnitudes of change in RH, as suggested by the HadISDH dataset versus the ERA-Interim dataset, show a relatively high correlation (mostly above 0.55).Given this high consistency between the HadISDH and the ERA-Interim datasets in terms of both the magnitude and sign of change in RH and also in the interannual variations (Figs.S4-S6), we decided to re-strict our subsequent analysis to the ERA-Interim dataset, recalling its denser global coverage compared to the HadISDH.
As RH is mostly dependent on changes in specific humidity (q), there is a dominant high correlation between the interannual variability in RH and q (Fig. S7).In accordance, the magnitude of observed change in these two variables showed a strong agreement for 1979-2014.Figure 3 summarizes the magnitude of change in specific humidity (q) as well as changes in specific humidity necessary to maintain constant RH as recorded in 1979.Specific humidity showed the strongest decrease in southwestern North America, the Amazonian region, southern South America, and the Sahel regions: a spatial pattern that is similar to the RH pattern.Given the evolution of air temperature for 1979-2014, these regions exhibited a deficit of water vapor on the order of −2 g kg −1 in order to maintain constant RH.

Spatial patterns of the dependency between RH and climate variables
We selected a range of regions (N = 14) worldwide (Fig. 4).For these selected regions, we assessed the connection between RH and some relevant climatic variables for the pe- riod 1979-2014.In addition, we defined the oceanic and continental sources of moisture corresponding to these regions using the FLEXPART model.We assessed the optimal lifetime for each region: during 4 days for regions 1-5 and 7-11, during 5 days for regions 6 and 12-13, and during 7 days for region 5. Figures 7-9 show some examples of the dependency between RH and different climate variables at the annual scale.Results for all regions at the seasonal and annual scales are presented in the Supplement.

Western Sahel
Figure 5a, b illustrate RH trends in the western Sahel using the HadISDH and ERA-Interim datasets.We also showed the distribution of the average annual moisture sources (E-P in millimeters) over this region for 1979-2014.As illustrated, the atmospheric moisture mostly comes from the western Sahel region itself, in addition to some oceanic sources located in the central eastern Atlantic Ocean.At the seasonal scale, there are some differences in the location and the intensity of the moisture sources, with more oceanic contribution during the boreal warm season.However, in both cases, the continental moisture seems to be the key source of humidity in the region (Figs.S23 and S37). Figure 5c shows different scatter plots summarizing the relationships between the de-trended annual series of RH and those of relevant climate variables (e.g., precipitation, air temperature, and SST).The interannual variability in RH in the region is correlated to changes in the total annual precipitation and the total annual land evapotranspiration in the continental source region.Specifically, the correlation between the de-trended annual RH and precipitation and land evapotranspiration is generally above 0.8 (p < 0.05).In contrast, RH shows negative correlations with air temperature and SST ratio over the oceanic source.While the correlation is statistically insignificant (p > 0.05), it suggests that higher differences between air temperature and SST reinforce lower annual RH.At the seasonal scale, we found similar patterns (Figs.S23 and S37).Nevertheless, in the warm season, a significant negative correlation with air temperature and SST ratio was observed.This pattern concurs with the significant increase in specific humidity (q) for 1979-2014; this is probably related to the high increase in land evapotranspiration (19.5 %, p < 0.05).

La Plata region
Figure 6 summarizes the corresponding results, but for the La Plata region (South America).Results indicate a general decrease in RH at the annual and seasonal scales using both the HadISDH observational data and the ERA-Interim dataset.
As depicted, the main humidity sources are located in the same region, combined with some other neighboring conti- nental areas over South America.A similar finding was also observed on the seasonal scale (Figs.S27 and S41).Similar to the western Sahel region, we found a significant association between the interannual variations in RH and precipitation and the land evapotranspiration in the continental source region.Similarly, we did not find any significant correlation between RH changes and the interannual variability in the oceanic evaporation in the oceanic source region or the ratio between air temperature in the continental target region and SST in the oceanic source region.Again, we found a negative correlation between RH and the air temperature / SST ratio, though it is statistically insignificant at the annual scale (p > 0.05).In the La Plata region, we noted a strong decrease in RH (−6.21 % decade −1 ) for 1979-2014, which agrees well with the strong decrease in absolute humidity.

Southwestern North America
Results for southwestern North America are also illustrated in Fig. 7.In accordance with both previously studied examples (western Sahel and La Plata), this region also exhibited a strong and positive relationship between the interannual vari-ability in RH and precipitation and land evapotranspiration.This pattern was also recorded for the boreal warm and cold seasons (Figs.S30 and S43).In this region, we found a strong negative trend in RH for 1979-2014, which concurs with the significant decrease in absolute humidity.We noted a significant increase in air temperature and air temperature-to-SST ratio, while a negative and statistically significant decrease in land evapotranspiration in the continental sources of moisture was observed.

Other regions
Other regions of the world (see Supplement) also showed strong dependency between the interannual variability in RH and that of land  tive correlations with RH in particular regions, including the western Sahel, La Plata, the west coast of the USA, centraleastern Europe, India, central North America, and the Amazonian region.This finding suggests that higher differences between air temperature in the target area and SST in the oceanic moisture region would favor decreased RH.

Land and ocean contribution to RH trends
Establishing a direct influence of land evapotranspiration on RH is a challenging task, also including any attempt to directly compare these influences with the possible contribution from oceanic evaporation and moisture transport.This is primarily because, apart from very humid regions, the increase in land evapotranspiration could be driven by increased precipitation, which is accompanied by anomalous RH conditions.This dependency explains the correlation found between precipitation and land evapotranspiration in some regions worldwide (Figs.S47 to S49).In cold and humid regions, land evapotranspiration is also related to the interannual variability in the AED (Figs.S50 to S52).Correspondingly, the magnitude of the oceanic evaporation may be insufficient to explain RH anomalies in the target region.
Taken together, the transport of moisture to any target region is a fundamental process.Hence, we assessed the contribution of land and ocean to precipitation, as represented by (E-P ).Overall, results reveal important differences among the analyzed regions, with statistically insignificant correlations found between the interannual variations in RH in some regions and ocean and land contribution to precipita-tion (Figs. 8 and 9).A similar pattern was observed at the seasonal scale, albeit with greater contribution during the cold season, especially in the regions where precipitation is mostly driven by western flows (e.g., western North America, western Europe, and Scandinavia) (Figs.S53-S56).Conversely, the land contribution to precipitation is rather complex, with strong spatial differences.At the annual scale, a positive and significant contribution to precipitation (E-P ) is found in regions 3 (central-eastern Europe), 6 (India), 7 (China), 9 (La Plata), and 11 (central US).This suggests that there are no generalized patterns in terms of the contribution of ocean and land to the interannual variability in RH, making it difficult to attribute RH trends to a unique driver.
Figure 10 illustrates the evolution of the land contribution (E-P ) to precipitation in the different target regions.We noted positive and significant changes in regions 2 (Scandinavia), 3 (central-eastern Europe), 5 (western Sahel), and 14 (eastern Sahel).A contradictory behavior is observed for region 9 (La Plata), with a statistically significant downward trend in land contribution (E-P ) to precipitation.At the seasonal scale, results suggest considerable differences (Figs.S57 and S58), with no clear positive or negative trends.
Changes in RH were more associated with those of land evapotranspiration across the selected regions (Fig. 11).In contrast, changes in annual RH did not correlate significantly with the observed changes in precipitation, air temperature-SST, and oceanic evaporation.The observed patterns were similar for both the warm and the cold seasons (Figs.S59 and  S60).These positive and significant correlations do not im-  Earth Syst.Dynam., 9, 915-937, 2018 www.earth-syst-dynam.net/9/915/2018/ply causation between land evapotranspiration and RH variations over space and time.Nevertheless, these findings suggest a role of land evapotranspiration in explaining the observed RH trends.Specifically, for many regions and at different temporal scales (i.e., seasonal and annual), changes in land contribution to precipitation show statistically significant positive correlation with changes in evapotranspiration and precipitation (Fig. S61).Again, this correlation does not imply a true causal relationship between RH variability and evapotranspiration, given the strong coupling between many of these controlling variables (e.g., precipitation, RH, and land evapotranspiration).The agreement between changes in the land contribution to precipitation and changes in land evapotranspiration, but also the fact that in many regions land contribution (E-P ) to precipitation is important (contributions close to 50 % or even higher; Table S1 in the Supplement), suggests the possible role of land evapotranspiration in explaining the spatial patterns of RH changes.Nevertheless, there is some uncertainty in this attribution since the good spatial agreement between RH, precipitation, and land evapotranspiration makes it difficult to accurately define the most dominant variable(s) that control the temporal variability in RH (Figs.S62-S64).
The complex spatial trends for different variables add another source of uncertainty to proper attribution of RH changes.Figure 12 illustrates the spatial distribution of the magnitude of change in annual and seasonal land evapotranspiration at the global scale from 1979 to 2014.As depicted, the spatial patterns of land evapotranspiration changes resemble those of RH in some regions (refer to Fig. 1), including -for example -the Canadian region, southwestern North America, and the western and eastern Sahel.Nevertheless, other regions showed a divergent pattern between both variables (e.g., the Guinea Gulf in Nigeria and Cameroon), where we noted a strong increase in land evapotranspiration, as opposed to RH changes.The Eurasian continent showed the main divergences between both variables.For example, in the western Europe, we noted a dominant RH decrease, which was not observed for land evapotranspiration.A similar pattern was observed over eastern China, with a dominant RH negative trend and a positive land evapotranspiration.Overall, the lack of significant spatial association between the magnitude of trends in RH and the magnitude of trends in evapotranspiration can be seen in the context of the strong spatial diversity of trends in these two variables at both annual and seasonal scales (Fig. S65).This complexity is similar to that found for trends in precipitation (Fig. S66).Thus, results suggest that while the variability in precipitation, RH, and land evapotranspiration show strong interannual associations, their observed trends are completely decoupled over space.This high spatial variability in trends at the global scale confirms that direct attribution of observed RH changes to land contributions is a challenge and quite a complex task.
In relation to the ocean evaporation, it is also quite difficult to establish this connection.In particular, it is not feasible to identify the moisture sources and the ocean contribution of the precipitation for each 0.5 • land pixel at the global scale.However, we believe that the analysis of the evolution of SST and oceanic evaporation for 1979-2014 and the evolution of the oceanic contribution to precipitation can give indications on some relevant patterns.Figure 13 illustrates the spatial distribution of the magnitude of change of annual and seasonal SST and oceanic evaporation.Figure S67 shows the spatial distribution of trend significance.As depicted, complex spatial patterns and high variability in the trends were observed, particularly for oceanic evaporation.Furthermore, the spatial distribution of the magnitude of change in annual and seasonal oceanic evaporation was not related to the SST changes (Fig. S68).Thus, although some regions showed positive changes in the oceanic evaporation, the amount of increase was much lower than that found for SST, which suggests that only SST changes do not drive evaporation changes (Fig. S69, Table S2).

Relative humidity trends
We assessed the temporal variability and trends in RH at the global scale using a dense observational network of meteorological stations (HadISDH) and reanalysis data (ERA-Interim).Results revealed high agreement of the interannual variability in RH using both datasets for 1979-2014.Recent studies have suggested a dominant decrease in observed RH during the last decade (e.g., Simmons et al., 2010;Willet et al., 2014).Our study suggests dominant negative trends of RH using the HadISDH dataset.This decrease is mostly linked to the temporal evolution of RH during the boreal warm season.In accordance with the HadISDH dataset, the ERA-Interim revealed dominant negative RH trends, albeit with a lower percentage of the total land surface compared to the HadISDH dataset.These differences cannot be attributed to the selected datasets, given that both mostly agree on the magnitude and sign of changes in RH.
Observed changes in RH were closely related to the magnitude and the spatial patterns of specific humidity changes.Results demonstrate a general deficit of specific humidity to maintain constant RH in large areas of the world, including central and southern North America, the Amazon and La Plata basins in South America, and the eastern Sahel.In other regions, RH increased in accordance with higher specific humidity.Studies suggest that changes in air temperature could partly cancel the effects of the atmospheric humidity to explain RH changes (e.g., McCarthy and Tuomi, 2004;Wright et al., 2010;Sherwood, 2010).Nevertheless, although air temperature trends showed spatial differences at the global scale over the past 4 decades (IPCC, 2013), our results show that air temperature is not the main driver of the observed changes of RH globally.There is a close resemblance between relative and specific humidity trends at the global scale, suggesting that specific humidity is the main driver in the magnitude and spatial pattern of RH trends.

Contribution of continental areas to changes in RH
Overall, there is an agreement among the interannual variability in precipitation, land evapotranspiration in the continental moisture source, and the interannual variability in RH in different regions.Nevertheless, considering gridded datasets at the global scale, we found that this good agreement is restricted only to the arid and semiarid regions.In humid regions, soil moisture is not constrained so land evapotranspiration variability is mostly driven by changes in the AED (Stephenson, 1990).This makes it difficult to unravel the possible direct contribution of land evapotranspiration to the variability in RH using statistical approaches and empirical information, particularly with the strong coupling among these variables.Land evapotranspiration is closely related to precipitation variability in arid and semiarid regions; increased land evapotranspiration tends to be caused primarily by increased precipitation, which is accompanied by corresponding RH anomalies.Also, RH may affect land evapotranspiration, both in arid and humid regions, given its impor-Earth Syst.Dynam., 9, 2018 www.earth-syst-dynam.net/9/915/2018/tant contribution to the aerodynamic component of the AED (Wang et al., 2012;Vicente-Serrano et al., 2014a).Nevertheless, although the interannual variability in these three variables is coupled in some regions, its long-term trends may strongly differ as a consequence of changes in precipitation, increasing influence of the AED, and also changes in land and atmospheric contribution to RH and precipitation.The 14 analyzed regions, in which FLEXPART was applied, show a relevant continental contribution (E-P ) to precipitation The average contribution is generally below 40 % of the total precipitation in some regions (e.g., western Europe, Scandinavia, or western North America), but exceeds 50 % in specific regions (e.g., Sahel and eastern China).
Therefore, it is reasonable to consider that changes in the contribution of continental areas to precipitation may affect land evapotranspiration processes and ultimately affect RH variability.Thus, our results suggest an influence of landatmosphere water feedbacks and recycling processes on RH trends.This is simply because more available soil humidity under favorable atmospheric and land conditions would result in more evapotranspiration and accordingly higher air moisture (Eltahir and Bras, 1996;Dominguez et al., 2006;Kunstmann and Jung, 2007).Recalling that the ocean surface contributes about 84 % of the water evaporated over the Earth (Oki, 2005), oceanic evaporation is highly important for continental precipitation (Gimeno et al., 2010).However, the continental humidity sources can be more im-portant than oceanic sources in some regions (e.g., the Sahel) (Wei et al., 2016a).In this context, our results concur with previous works.For example, numerous model-based studies have supported an influence of land evaporation processes on air humidity, and precipitation over land surface (e.g., Bosilovich and Chern, 2006;Dirmeyer et al., 2009;Goessling and Reick, 2011) recycling is strongly important in some regions of the world, such as China and central Asia, the western part of Africa, and central South America (Pfahl et al., 2014;van der Ent et al., 2010, and in general in semiarid and desert areas worldwide (Miralles et al., 2016).All these studies assessed the role of continental evapotranspiration on average precipitation conditions, with few studies focusing on the possible impacts of changes in soil moisture and evapotranspiration on RH.In Europe, Rowell and Jones ( 2006) concluded for future climate scenarios that reduced evaporation in summer will drop RH and hence continental rainfall.
Although our study was limited to specific regions across the world, results indicate that humidity in the analyzed regions substantially originates over continents.This finding concurs with some regional studies that defined sources of moisture (e.g., Nieto et al., 2014;Gimeno et al., 2010;Drumond et al., 2014;Ciric et al., 2016).Overall, the spatial differences in the possible attribution of the observed changes in RH to changes in land evapotranspiration are important.Nevertheless, in some regions our results suggest a robust contribution of land processes to the interannual variability www.earth-syst-dynam.net/9/915/2018/Earth Syst.Dynam., 9, 915-937, 2018 in RH.A representative example is La Plata, where a strong decrease (−6.6 % from 1979 to 2014) in RH is suggested by both observations and ERA-Interim.This region did not exhibit a significant trend in precipitation, but conversely there is a significant decrease in absolute humidity and land evapotranspiration.In the La Plata region, the oceanic evaporation in the source region has shown a significant increase (6.33 %) since 1979.However, this increase seems to be insufficient to favor an increase in RH values.Herein, although the average oceanic contribution (E-P ) to precipitation is slightly higher (54 %), compared to the continental contribution (46 %), the interannual variability in RH is positively correlated with the interannual variability in land E-P rather than oceanic E-P .Moreover, this region exhibited a significant decrease in land contribution to precipitation between 1979 and 2014.
In the western US, the large decrease in RH (−6.2 %) corresponded to a large decrease in the absolute humidity (−0.58 g kg −1 ).However, it is difficult to attribute this pattern to changes in land evapotranspiration, given the low (37 %) continental contribution to precipitation.Wei et al. (2016b) showed that the transport of atmospheric moisture from the Pacific is the main contributor to the interannual variability in precipitation in the region.In this sense, we found a strong relationship between the interannual variability in RH and the oceanic E-P , albeit with insignificant trends in the oceanic and continental contribution (E-P ) to precipitation in the region.Here, in the absence of significant changes in oceanic evaporation and contribution to precipitation, it is reasonable to consider that the decrease in absolute humidity is linked to the atmospheric circulation that controls moisture transport in the region (Wei et al., 2016b), but the decrease in RH could also be favored by the trend toward higher differences between air temperature in the land target region and the SST in the source region.
Therefore, the relationships among RH, land evapotranspiration variability, and changes in the contribution of continental areas can be extremely complex.Similarly, the relationships between RH and land evapotranspiration are rather complex and cannot be easily interpreted.
Nonetheless, albeit with the strong uncertainty existing at the global scale, it is reasonable to consider that changes in soil moisture budget, combined with land evapotranspiration, could somehow impact the observed RH trends.In the same context, there is strong evidence that low levels of soil moisture and land evapotranspiration are usually accompanied by a reinforcement of low RH, particularly during drought episodes.Under these circumstances, the suppression of the latent heat flows from the soil to the atmosphere would enhance soil and vegetation warming and sensible heat, inducing air temperature rise and the reinforcement of heat waves (Hirschi et al., 2011).Seneviratne et al. (2002) showed that vegetation control on transpiration might contribute significantly to enhancement of summer drying, particularly when soil water is limited.Other studies confirmed this finding, employing both observational data (e.g., Hirschi et al., 2011) and model outputs (e.g., Seneviratne et al., 2006;Fischer et al., 2007).

Contribution of oceans to changes in RH
We indicated strong differences among the 14 analyzed regions in terms of the contribution of Oceanic water bodies to RH variability.In some regions (e.g., western North America, India, Scandinavia, and Amazonia), the interannual variability in RH is closely related to oceanic E-P , indicating that changes in oceanic evaporation, combined with the processes of atmospheric moisture transport to target regions, play a main role in explaining changes in RH.Moisture advection is the main driver of precipitation variability (Trenberth, 1999;Wei et al., 2016a).Nevertheless, there is some uncertainty in recent trends in moisture advections from oceanic areas (Zahn and Allan, 2013).As such, it is difficult to determine -in general terms -whether the strong complexity of RH changes identified at the global scale is driven by changes in moisture advections from oceanic areas, given that -in the context of changing climate -SST, oceanic evaporation, and oceanic contribution to precipitation in the target regions can jointly account for the possible influence of oceans on RH variability in some regions.
Different modeled climate studies suggested strong differences between land and ocean RH trends, as a consequence of the different warming rates between oceanic and continental areas (e.g., Joshi et al., 2008;Dessler and Sherwood, 2009;O'Gorman and Muller, 2010).As the warming rates are generally slower over oceans, the specific humidity of air advected from oceans to continents would increase more slowly than the saturation-specific humidity over land, causing a reduction in RH (Rowell and Jones, 2006).Due to this effect, RH will not remain constant in areas located very far from humidity sources, as warmer air temperatures under limited moisture humidity would reduce RH (Pierce et al., 2013).This physical process could explain observed RH trends in some regions (e.g., Amazonia, where around 68 % of the average atmospheric moisture originates over oceanic sources).In this region RH is correlated with the oceanic E-P , although there is no change in the oceanic contribution to precipitation from 1979 to 2014.Moreover, RH strongly decreased at the annual scale (−7.7 %), as a consequence of a decrease in the absolute humidity (−0.56 g kg −1 ), accompanied by an increase in air temperature (1.09 • C).While SST in the oceanic source region slightly increased by 0.33 • C, other variables (i.e., oceanic evaporation, precipitation, and land evapotranspiration) did not exhibit significant changes.Herein, we believe that the mechanisms proposed by Sherwood and Fu (2014) to explain decreased RH are applicable to our study.They attributed RH decrease to subsaturated oceanic moisture supply, which compensates for air temperature increase.In this study, this mechanism is also capable of explaining RH variability, given that the difference in warming rates between SST in the oceanic moisture source Earth Syst.Dynam., 9, 915-937, 2018 www.earth-syst-dynam.net/9/915/2018/region and air temperature of Amazonia is increasing.This feature is also supported by the significant correlation found between the interannual variability in RH and the ratio between air temperature over land and SST.A similar pattern is found for the eastern Sahel, a region in which continental recycling is particularly important (Wei et al., 2016a).This region witnessed a strong decrease in RH, but it was not compensated for by increased precipitation.Although there is no significant correlation between RH and the ratio between air temperature and SST, the latter variable shows a significant increase.This could reinforce the drying effect of the suggested trend toward a lower moisture supply from the oceanic source, especially with the significant negative trend in oceanic contribution to precipitation in this region.
In other regions, there are no empirical relationships that confirm the impact of oceanic moisture supply processes on RH changes.A possible explanation of these contrasting findings could be related to the low differences in the warming rates between the oceanic sources and the majority of the continental target areas.We found that -in most cases -these differences were not strong enough to generate a clear effect at the global scale, particularly with the available number of observations.
Although oceanic evaporation is decisive on moisture supply to continental regions (Gimeno et al., 2010), several processes, which are not considered in this study, may strongly affect RH and precipitation in continental areas.A global warming signal does not necessarily imply above-normal oceanic evaporation.Here, we indicated that oceanic evaporation trends for 1979-2014 showed strong spatial variability, with dominant positive trends.Nevertheless, large areas also exhibited insignificant trends and even negative evaporation trends.While SST increase is mainly associated with radiative processes, evaporation processes are mainly controlled by a wide range of meteorological variables that impact the aerodynamic and radiative components of the AED rather than SST alone (McVicar et al., 2012b).Thus, changes in solar radiation and wind speed can also influence the evaporation evolution (Yu, 2007;Kanemaru and Masunaga, 2013).
Due to the unlimited water availability over oceans, air vapor pressure deficit would be expected to be driven by the Clausius-Clapeyron relation.Nevertheless, the observed evolution indicates that the global mean precipitation and evaporation is not scaled by the Clausius-Clapeyron relation (Held and Soden, 2006).As such, given the slow oceanic evaporation trends in large regions of the world, RH trends in some of the analyzed regions could be seen in the context of a lower water supply to maintain constant RH, particularly with air temperature rise.
Herein, we have not considered possible changes in other variables that could explain the weak relationship among RH, oceanic evaporation, and oceanic moisture contribution, for example, the transport mechanisms between the source and target regions, which are key variables in some regions like the western North America (Wei et al., 2016b).More-over, moisture source regions are not stationary, as the intensity of humidity can vary greatly from one year to another (Gimeno et al., 2013).This aspect could be another source of uncertainty in the explanatory factors of current RH trends.Furthermore, other different factors that control atmospheric humidity and RH have not been addressed in this study.Sherwood (1996) indicated that RH distributions are strongly controlled by dynamical fields rather than local air temperatures.This suggests that atmospheric circulation processes could largely affect the temporal variability and trends in RH.A range of studies indicate noticeable RH changes in response to low-frequency atmospheric oscillations, such as the Atlantic Multidecadal Oscillation (AMO) and El Niño-Southern Oscillation (e.g., McCarthy and Toumi, 2004;Zhang et al., 2013), regional circulation (Wei et al., 2016a, b), and changes in the Hadley cell (HC) (Hu and Fu, 2007).Wright et al. (2010) employed a global climate model under double CO 2 concentrations to show that tropical and subtropical RH is largely dependent on a poleward expansion of the Hadley cell: a deepening of the height of convective detrainment, a poleward shift of the extratropical jets, and an increase in the height of the tropopause.Also, Laua and Kim (2015) assessed changes in the HC under CO 2 warming from the Coupled Model Intercomparison Project Phase 5 (CMIP5) model projections.They suggest that strengthening of the HC induces atmospheric moisture divergence and reduces tropospheric RH in the tropics and subtropics; this spatial pattern resembles those areas with negative RH trends in the Northern Hemisphere as well as the Southern Hemisphere, as described in this study.

Conclusions
This study analyzed relative humidity (RH) trends at the global scale using observations and ERA-40 data.It extended further to link RH trends with a range of variables, which can give indications on the possible oceanic and continental contribution to RH trends.Our results suggest significant RH trends over many regions worldwide and a generally dominant negative trend at the global scale.
The strong diversity in the observed RH trends highlights the complex and divergent role of different physical processes, including both dynamic and thermodynamic components.In general the supply of specific humidity is a main source of the observed RH trends since there is a high agreement between RH and specific humidity trends.This finding suggests that the evolution of specific humidity in vast areas of the world has not provided the necessary humidity to maintain constant RH according to the observed warming trends.
This study also analyzed the possible contribution of continental and oceanic moisture supply to explaining the magnitude and spatial patterns of RH trends.For this purpose, 14 regions were defined and the contribution of continental www.earth-syst-dynam.net/9/915/2018/Earth Syst.Dynam., 9, 915-937, 2018 and oceanic sources to RH was assessed using a Lagrangian scheme.Results indicate that no single physical mechanism can be responsible for the observed trends in RH at the global scale.Our findings stress two hypotheses, (i) land water supply and (ii) insufficient oceanic moisture supply to maintain constant continental RH, could act together to explain recent RH trends.However, although it is quite difficult to establish a direct causality between RH and different underlying processes using different empirical sources, the observed decrease in RH in some regions (e.g., La Plata) can be linked to the lower water supply from land evapotranspiration.In other regions, the empirical relationships suggest dynamic and thermodynamic mechanisms related to moisture supply from oceanic source regions (e.g., Amazonia and western North America).Taken together, these physical mechanisms could coexist in some regions given the strong relationship found among precipitation, RH, and land evapotranspiration.
Overall, this study confirms the strong complexity of determining a general physical process that may explain the complex spatial patterns of RH trends, particularly at the global scale.As such, further research is still needed to unravel the complex physical factors driving the dominant RH negative trends over large continental regions.The availability of long-term historical and reanalyses data and the advancement of modeling approaches are potential assets in any future research to explore whether the land and oceanic processes drive the observed RH trends.This is important, given that understanding current RH is relevant in hydroclimatic research due to its impacts on atmospheric evaporative demand, crop development and yield, forest fire risk, and bioclimatic comfort in addition to other hydrological processes.
Author contributions.SMVS, RN, and LG designed the research.SMVS and RN performed research.SMVS drafted the paper.CAM, AD, AEK, FDC, MTB, and MPG analyzed data, and all authors contributed to writing the paper.
Competing interests.The authors declare that they have no conflict of interest.Special issue statement.This article is part of the special issue "The 8th EGU Leonardo Conference: From evaporation to precipitation: the atmospheric moisture transport".It is a result of the 8th EGU Leonardo Conference, Ourense, Spain, 25-27 October 2016.

Figure 1 .
Figure 1.Spatial distribution of the magnitude of change of RH (% per decade) over the period 1979-2014 from the HadISDH (a) and ERA-Interim datasets (b).Results are provided for the boreal cold (October-March) and warm (April-September) seasons and annually.

Figure 2 .
Figure 2. Relative frequencies (%) of the RH magnitude of change in the HadISDH and ERA-Interim datasets.Color bar plots represent the percentage of stations (from HadISDH) and world regions (from ERA-Interim) with statistically significant positive trends at p < 0.05 (blue), statistically insignificant positive trends (cyan), statistically insignificant negative trends (orange), and statistically significant negative trends (red).

Figure 3 .
Figure 3. Spatial distribution of the seasonal and annual magnitudes of change in specific humidity (g kg 1 ) (a) and the deficit/surplus of specific humidity to maintain constant RH at the levels of 1979 (b) according to land air temperature evolution (from the CRU TS v.3.23 dataset) for 1979-2014.

Figure 4 .
Figure 4. Spatial distribution of the selected 14 world regions, based on the high consistency in RH trends between the HadISDH and the ERA-Interim datasets.These regions were selected for the identification of the oceanic and land humidity sources by means of the FLEXPART scheme.

Figure 5 .
Figure 5. (a) Annual RH humidity trends in the western Sahel (region 6); (b) average (E-P ) > 0 at the annual scale to identify the main humidity sources in the region (mm yr −1 ).(c) Relationships between the de-trended annual RH and the de-trended annual variables for 1979-2014.(d) Annual evolution of the different variables corresponding to the western Sahel region.The magnitude of change and their corresponding statistical significance are indicated for each variable.

Figure 11 .
Figure 11.Relationships between the average annual magnitudes of change in RH identified for each one of the 14 analyzed regions and the annual magnitudes of change in precipitation, the ratio between air temperature and SST, oceanic evaporation, and land evapotranspiration.

Figure 12 .
Figure 12.Spatial distribution of the magnitude of change in the annual and seasonal land evapotranspiration (1979-2014) and their corresponding statistical significance.

Figure 13 .
Figure 13.Annual and seasonal magnitude of change of SST and OAFlux oceanic evaporation for 1979-2014.
evapotranspiration in the land moisture sources.Some examples include western Europe, centraleastern Europe, southeastern Europe, Turkey, India, and the eastern Sahel.In contrast, other regions showed a weak cor- relation between the temporal variability in RH and land evapotranspiration in the moisture source region.A representative example is China, which witnessed a strong decrease in RH for 1979-2014.Nevertheless, the air temperature / SST ratio in the oceanic moisture sources also exhibited negawww.earth-syst-dynam.net/9/915/2018/Earth Syst.Dynam., 9, 915-937, 2018