Parameterizations of cloud feedback in a radiative-convective model

The effect of cloud feedback on the response of a radiative-convective model to a change in cloud model parameters, atmospheric CO₂ concentration, and solar constant has been studied using two different parameterization schemes. The method for simulating the vertical distribution of both cloud cover and cloud optical thickness, which depends on the relative humidity and on the saturation mixing ratio of water vapor, respectively, is the same in both approaches, but the schemes differ with respect to modeling the water vapor profile. In scheme I atmospheric water vapor is coupled to surface parameters, while in scheme II an explicit balance equation for water vapor in the individual atmospheric layers is used. For both models the combined effect of feedbacks due to variations in lapse rate, cloud cover, and cloud optical thickness results in different relationships between changes in surface temperature, planetary temperature, and cloud cover. Specifically, for a CO₂ doubling and a 2% increase in solar constant, in both models the surface warming is reduced by cloud feedback, in contrast to no feedback, with the greater reduction in scheme I as compared to that of scheme II.     

Physical,radiative, and dynamical processes within a nighttime marine stratus cloud

Observations taken by aircraft and conventional platforms are used to investigate dynamical, physical, and radiative processes within a marine stratus cloud during the Canadian Atlantic Storms Program (CASP) II field project which took place over the east coast of Canada. Stratus which formed over the ocean on February 6, 1992 during the nighttime, is studied to analyze cloud top and base processes. The cloud was supercooled during the study period. Fluctuations and fluxes are calculated along constant flight altitude legs approximately 100 km long in space. The scales of structures larger than 5 km are removed from the analysis using a running average technique. Droplet spectra obtained by a forward scattering spectrometer probe (FSSP) were used in a 1-D radiative transfer model to calculate infrared (IR) fluxes and radiative heating rates. A heat conservation equation was used to estimate vertical air velocity (w _a ) within the cloud. The results showed that, because of a warmer ocean surface, significant moisture and heat were transferred from the ocean surface to the boundary layer. The cloud base was at about 400 m height and the top was at about 1.4 km.w _a at the cloud base was estimated about 5 cm s^–1. Strong IR cooling rate at the cloud top was calculated to be 75°C day^–1 for a 100 m thick layer. Negative skewness inw _a , suggesting narrow downdrafts, was likely due to radiative cooling at the cloud top. The entrainment velocity was found to be about 1.5 cm s^–1 at cloud top. Mean moisture and heat fluxes within the cloud were estimated to be comparable to those from the ocean surface. Vertical air velocity at the cloud top due to radiative cooling was found to be about –40 cm s^–1.     

Oscillations with planetary wave periods in variations in the ionospheric parameters over Eastern Siberia

The results of a periodogramanalysis of the variations in the ionospheric parameters, measured using the vertical radio sounding method at midlatitude Irkutsk observatory (Eastern Siberia), are presented. The 1984–1986 period of observations was used. It has been indicated that the statistically significant oscillations with periods typical of planetary waves are present in the variations in f ₀E_s, f _bE_s, h′E_s, f _min, f ₀F2, and h′F.

     

On the 18-day quasi-periodic oscillation in the ionosphere

The presence and persistence of an 18-day quasi-periodic oscillation in the ionospheric electron density variations were studied. The data of lower ionosphere (radio-wave absorption at equivalent frequency near 1 MHz), middle and upper ionosphere (critical frequencies f₀E and f₀F2) for the period 1970–1990 have been used in the analysis. Also, solar and geomagnetic activity data (the sunspot numbers Rz and solar radio flux F10.7 cm, and a_N index respectively) were used to compare the time variations of the ionospheric with the solar and geomagnetic activity data. Periodogram, complex demodulation, auto- and cross-correlation analysis have been used. It was found that 18-day quasi-periodic oscillation exists and persists in the temporal variations of the ionospheric parameters under study with high level of correlation and mean period of 18–19 days. The time variation of the amplitude of the 18-day quasi-periodic oscillation in the ionosphere seems to be modulated by the long-term solar cycle variations. Such oscillations exist in some solar and geomagnetic parameters and in the planetary wave activity of the middle atmosphere. The high similarities in the amplitude modulation, long-term amplitude variation, period range between the oscillation of investigated parameters and the global activity of oscillation suggests a possible solar influence on the 18-day quasi-periodic oscillation in the ionosphere.     

Ozone content of the atmosphere. 1. Climatological characteristics of the total ozone content in the 100–110° E longitudinal zone

The results of a periodogramanalysis of the variations in the ionospheric parameters, measured using the vertical radio sounding method at midlatitude Irkutsk observatory (Eastern Siberia), are presented. The 1984–1986 period of observations was used. It has been indicated that the statistically significant oscillations with periods typical of planetary waves are present in the variations in f ₀E_s, f _bE_s, h′E_s, f _min, f ₀F2, and h′F.     

Partitioning resistance to overland flow on rough mobile beds

For overland flows transporting predominantly bed load over rough mobile beds without rainfall, resistance to flow f may be divided into four components: surface resistance f_s, form resistance f_f, wave resistance f_w, and bed‐mobility resistance f_m. In this study it is assumed that f = f_s + f_f + f_w + f_m, and an equation is developed for each component. The equations for f_s and f_f are borrowed from the literature, while those for f_w and f_m are developed from two series of flume experiments in which the beds are covered with various concentrations of large‐scale roughness elements. The first series consists of 65 experiments on fixed beds, while the second series contains 194 experiments on mobile beds. All experiments were performed on the same slope (S = 0·114) and with the same size of sediment (D = 0·00074 m). The equations for f_w and f_m are derived by a combination of dimensional analysis and regression analysis. The analyses reveal that the major controls of f_w and f_m are the Froude number F and the concentration of the roughness elements C_r. When the equations for f_w and f_m are summed, the C_r terms cancel out, leaving f_w+m = 0·63F^?2. An equation is developed that predicts total f, and the contributions of f_s, f_f, f_w and f_m to f are computed from the series 1 and 2 experiments. An analysis of the first series reveals that in clear‐water flows over fixed beds, f_w accounts for 52 per cent of f. A similar analysis of the second series indicates that in sediment‐laden flows over mobile beds f_w comprises 37 per cent and f_m 32 per cent of f, so that together f_w and f_m account for almost 70 per cent of f. Finally, regression analyses indicate that where F > 0·5, f_w and f_m each vary with F^?2 and f_w/f_m = 1·18. The equation developed here for predicting total f applies only to the range of hydraulic, sediment, and bed roughness conditions represented by the experimental data. With additional data from a broader range of conditions the same methodology as employed here could be used to develop a more general equation. Copyright © 2006 John Wiley & Sons, Ltd.     

WIND observations of coherent electrostatic waves in the solar wind

The time domain sampler (TDS) experiment on WIND measures electric and magnetic wave forms with a sampling rate which reaches 120 000 points per second. We analyse here observations made in the solar wind near the Lagrange point L1. In the range of frequencies above the proton plasma frequency f_pi and smaller than or of the order of the electron plasma frequency f_pe, TDS observed three kinds of electrostatic (e.s.) waves: coherent wave packets of Langmuir waves with frequencies f ≃ f_pe, coherent wave packets with frequencies in the ion acoustic range f_pi ≥ f ≥ f_pe, and more or less isolated non-sinusoidal spikes lasting less than 1 ms. We confirm that the observed frequency of the low frequency (LF) ion acoustic wave packets is dominated by the Doppler effect: the wavelengths are short, 10 to 50 electron Debye lengths λ_D. The electric field in the isolated electrostatic structures (IES) and in the LF wave packets is more or less aligned with the solar wind magnetic field. Across the IES, which have a spatial width of the order of ≃25_D, there is a small but finite electric potential drop, implying an average electric field generally directed away from the Sun. The IES wave forms, which have not been previously reported in the solar wind, are similar, although with a smaller amplitude, to the weak double layers observed in the auroral regions, and to the electrostatic solitary waves observed in other regions in the magnetosphere. We have also studied the solar wind conditions which favour the occurrence of the three kinds of waves: all these e.s. waves are observed more or less continuously in the whole solar wind (except in the densest regions where a parasite prevents the TDS observations). The type (wave packet or IES) of the observed LF waves is mainly determined by the proton temperature and by the direction of the magnetic field, which themselves depend on the latitude of WIND with respect to the heliospheric current sheet.     

The empirical relationships between M s, m b and M L for China and vicinity

Data from 753 earthquakes are used to determine a relationship between surface-wave magnitude (M _s) and bodywave magnitude (m _b), and from 541 earthquakes to determine a relationship between surface-wave magnitude (M _s) and local magnitude (M _L) for China and vicinity: M _s=0.9883 m _b-0.0420, M _s=0.9919 M _L-0.1773. The relationship of M _s versus m _b is obtained for 292 events occurred in the Chinese mainland in the time period from 1964 to 1996, 291 events occurred in Taiwan in the time period from 1964 to 1995 and 170 events occurred in the surrounding area. Standard deviation of the fitting is 0.445. Relationship of M _s versus M _L is obtained for 36 events occurred in the Chinese mainland, 293 events occurred in Taiwan, China and 212 events occurred in the surrounding area. The total amount is 541 events. Standard deviation of the fitting is 0.4673. The uncertainties of the converted M _s in different magnitude intervals can be estimated using complementary cumulative distribution function (CCDF). In the relationship of M _s versus m _b, taking ±0.25 as a range of uncertainties, in magnitude interval m _b 4.0–4.9, the probabilities for the converted M _s taken value less than (M _s-0.25) and more than (M _s+0.25) are 17% and 27% respectively. Similarly, we have probabilities for m _b 5.0–5.9 are 34% and 20% and that for m _b 6.0–6.9 are 11% and 47%. In the relationship of M _s versus M _L, if the range of uncertainties is still taken as ±0.25, the corresponding probabilities for magnitude interval M _L 4.0–4.9 are 22% and 38%, for M _L 5.0–5.9 are 20% and 15% and for magnitude interval M _L 6.0–6.9, are 15% and 29%, respectively. The relationships developed in this paper can be used for the conversion of one magnitude scale into another magnitude scales conveniently. The estimation of uncertainties described in this paper is more accurate and more objective than the usual estimation expressed by deviation. The estimations described in this paper indicate various dispersions in different magnitude intervals of original data. The estimations of uncertainties described by probabilities can be well connected with the total estimations of uncertainties in seismic hazard assessment.     

Modelling snowpack surface temperature in the Canadian Prairies using simplified heat flow models

Three practical schemes for computing the snow surface temperature T_s, i.e. the force–restore method (FRM), the surface conductance method (SCM), and the Kondo and Yamazaki method (KYM), were assessed with respect to T_s retrieved from cloud‐free, NOAA‐AVHRR satellite data for three land‐cover types of the Paddle River basin of central Alberta. In terms of R², the mean T_s, the t‐test and F‐test, the FRM generally simulated more accurate T_s than the SCM and KYM. The bias in simulated T_s is usually within several degrees Celsius of the NOAA‐AVHRR T_s for both the calibration and validation periods, but larger errors are encountered occasionally, especially when T_s is substantially above 0 °C. Results show that the simulated T_s of the FRM is more consistent than that of the SCM, which in turn was more consistent than that of the KYM. This is partly because the FRM considers two aspects of heat conduction into snow, a stationary‐mean diurnal (sinusoidal) temperature variation at the surface coupled to a near steady‐state ground heat flux, whereas the SCM assumes a near steady‐state, simple heat conduction, and other simplifying assumptions, and the KYM does not balance the snowpack heat fluxes by assuming the snowpack having a vertical temperature profile that is linear. Copyright © 2005 John Wiley & Sons, Ltd.     

Surge in sulphur and halogen degassing from Ambrym volcano,Vanuatu

Volcanoes provide important contributions to atmospheric budgets of SO₂ and reactive halogens, which play significant roles in atmospheric oxidative capacity and radiation. However, the global source strengths of volcanic emissions remain poorly constrained. These uncertainties are highlighted here by the first measurements of gas emission rates from Ambrym volcano, Vanuatu. Our initial airborne ultraviolet spectroscopic measurements made in January 2005 indicate fluxes of 18–270 kg s^-1 of SO₂, and 62–110 g s^-1 of BrO, into the atmosphere, placing Ambrym amongst the largest known contemporary point sources of both these species on Earth. We also estimate high Cl and F fluxes of ~8–14 and ~27–50 kg s^-1, respectively, for this period. Further observations using both airborne and spaceborne remote sensing reveal a fluctuating SO₂ output between 2004 and 2008, with a surge in the first half of 2005, and underline the substantial contribution that a single passively degassing volcano can make to the atmospheric budget of sulfur and halogens.     

Pan evaporation modelling and changing attribution analysis on the Tibetan Plateau (1970–2012)

Pan evaporation (E_p) is an important indicator of water and energy and the decline of E_p has been reported in many regions over the last decades. The climate and E_p are dependent on each other. In this study, the temporal trends of E_p and main E_p drivers, namely mean air temperature (T_a), wind speed (u), global solar radiation (R_s), net long‐wave radiation(R_nl) and vapour pressure deficit (D) from 1970 to 2012, were calculated on the basis of 26 meteorological stations on the Tibetan Plateau. The arithmetic average of E_p from 26 stations decreased with the rate of ?11.91 mm a^?2; the trends of R_s, R_nl, T_a, u and D were ?1.434 w m^?2 decade^?1, 0.2511 w m^?2 decade^?1, 0.3590°C decade^?1, ?0.2376 m s^?1 decade^?1 and 9.523 Pa decade^?1, respectively. The diffuse irradiance is an essential parameter to model E_p and quantify the contribution of climatic factors to changing E_p. 60 724 observations of R_s and diffuse solar irradiance (R_d) from seven of the 26 stations were used to develop the correlation between the diffuse fraction (R_d/R_s), and the clearness index (R_s/R_o). On the basis of the estimation of the diffuse component of R_s and climatic data, we modified the PenPan model to estimate Chinese micro‐pan evaporation (E_p) and assess the attribution of E_p dynamics using partial derivatives. The results showed that there was a good agreement between the observed and calculated daily E_p values. The observed decrease in E_p was mostly due to declining wind speed (?13.7 mm a^?2) with some contributions from decreasing solar irradiance (?3.1 mm a^?2); and the increase of temperature had a large positive effect (4.55 mm a^?2) in total whilst the increase of R_nl had insignificant effect (0.35 mm a^?2) on E_p rates. The change of E_p is the net result of all the climatic variables. Copyright © 2014 John Wiley & Sons, Ltd.     

A bootstrap calculation of confidence regions for proportions of sediment contributed by different source areas in a ‘fingerprinting’ model

When ‘fingerprinting’ is used to identify what proportions P_s (s = 1, …, g) of suspended sediment come from g different source areas, measures of the uncertainties in estimates of the P_s are also required. These uncertainties are influenced by two kinds of correlation whose effects are rarely recognized in the literature. These are (i) correlation between the estimated P_s because they must add to 1 and (ii) correlation between the geochemical tracers measured in sediment samples. This paper uses bootstrap procedures to identify joint confidence regions for the estimated proportions [responding to correlation of type (i)] and to explore alternatives to the ‘standard’ least‐squares criterion used to estimate the proportions when tracer measurements are correlated [correlation of type (ii)]. Using a limited dataset with three sediment source areas for illustration (g = 3), results were obtained from 5000 bootstrap samples, using two criteria (standard and generalized least squares, GLS) with two inequality constraints: (a) 0 ≤ P_s ≤ 1, where P_s is the fraction of suspended sediment contributed by the s‐th source area (s = 1, 2, 3) and (b) 0 < P_s < 1, which, the paper argues, better represents reality. Approximate 95% confidence regions for the P_s, given by the two criteria and two inequality constraints, were compared. Using inequality constraint (a), the confidence region given by the GLS criterion was slightly smaller than that given by the standard; using constraint (b), the two confidence regions' boundaries were almost identical, suggesting that the effects of correlations between tracers were not large for the dataset used. For both criteria, the scatter amongst estimated proportions P_s obtained by bootstrapping was large, raising issues concerning the efficiency of sampling and the allocation of sampling effort, both in source areas and in transported suspended sediment. The results suggest that apparently small differences in the constraints applied to the proportions P_s can give quite different numerical results. Copyright © 2014 John Wiley & Sons, Ltd.     

Atmospheric carbon dioxide and the long-term control of the Earth’s climate

A CO₂-weathering model has been used to explore the possible evolution of the Earth’s climate as the Sun steadily brightened throughout geologic time. The results of the model calculations can be described in terms of three, qualitatively different, “Megaclimates”. Mega-climate 1 resulted from a period of rapid outgassing in the early Archean, with high, but declining, temperatures caused by the small weathering rates on a largely water-covered planet. Mega-climate 2 began about 3 Gyear ago as major continental land masses developed, increasing the weathering rate in the early Proterozoic and thereby depleting the atmospheric CO₂ concentration. This process produced the first Precambrian glaciations about 2.3 Gyear ago. During Mega-climate 2, evolutionary biological processes increased the surface weatherability in incremental steps and plate tectonics modulated the CO₂ outgassing rate with an estimated period of 150 Myear (approximately one-half the period for the formation and breakup of super continents). Throughout Mega-climate 2 the surface temperature was controlled by variations in the atmospheric CO₂ level allowing transitions between glacial and non-glacial conditions. The results of the model for Mega-climate 2 are in agreement with the occurrence (and absence) of glaciations in the geologic record. Extending the model to the future suggests that CO₂ control of the Earth’s temperature will no longer be able to compensate for a solar flux that continues to increase. The present level of atmospheric CO₂ is so small that further reduction in CO₂ cannot prevent the Earth from experiencing Mega-climate 3 with steadily increasing surface temperatures caused by the continued brightening of the Sun. During Mega-climate 3, the main danger to the biosphere would come not from an increasing temperature but from a decreasing (rather than an increasing) CO₂ level which could, in time, fall below 0.5 PAL, causing serious damage to the biosphere. Fortunately, the rates of change due to solar brightening are slow enough that Mega-climate 3 appears to pose no threat to the biosphere for the next 0.5-2 Gyear.     

Earthquake damage detection in the Imperial County Services Building III: Analysis of wave travel times via impulse response functions

The majority of structural health monitoring methods are based on detecting changes in the modal properties, which are global characteristics of the structure, and are not sensitive to local damage. Wave travel times between selected sections of a structure, on the other hand, are local characteristics, and are potentially more sensitive to local damage. In this paper, a structural health monitoring method based on changes in wave travel times is explored using strong motion data from the Imperial Valley Earthquake of 1979 recorded in the former Imperial County Services (ICS) Building, severely damaged by this earthquake. Wave travel times are measured from impulse response functions computed from the recorded horizontal seismic response in three time windows—before, during, and after the largest amplitude response, as determined from previous studies of this building, based on analysis of novelties in the recorded response. The results suggest initial spatial distribution of stiffness consistent with the design characteristics, and reduction of stiffness following the major damage consistent with the spatial distribution of the observed damage. The travel times were also used to estimate the fundamental fixed-base frequency of the structure f₁ (assuming the building deformed as a shear beam), and its changes during this earthquake. These estimates are consistent with previous estimates of the soil–structure system frequency, f_sys, during the earthquakes (f₁<f_sys as expected from soil–structure interaction studies), and with other estimates of frequency (f₁ from ETABS models, and f_sys from ambient vibration tests, and “instantaneous” f₁ from high-frequency pulse propagation).     

Scattering and anelastic attenuation of seismic energy in northern Greece

The relative contribution of scattering (Q _s ^–1 ) and intrinsic (Q _i ^–1 ) attenuation to the totalS-wave attenuation for the frequencies of 1.5, 3.0, 6.0 and 12.0 Hz has been studied by applying the radiative energy transfer theory, Data of local earthquakes which occurred in northern Greece and were recorded by the permanent telementered network of the Geophysical Laboratory of the University of Thessaloniki have been used. The results show that in this area the scattering attenuation is dominant over all frequencies while intrinsic attenuation is significantly lower. The estimatedQ _s ^–1 andQ _i ^–1 values have frequency dependences off ^–0.72 andf ^–0.45, respectively. The frequency dependence ofQ _s ^–1 is the same as that of the codaQ _c ^–1 , obtained by applying the single scattering model, which probably implies that the frequency dependence of the coda wave attenuation is attributed to the frequency dependence of the scattering attenuation.Q _c ^–1 values are very close to scattering attenuation for short lapse times, (10–20 sec), and intermediate between scattering and intrinsic attenuation for the longer lapse times, (50–100 sec). This difference is explained as the result of the depth-dependent attenuation properties and the multiple scattering effects.     

CO2 absorption of sandy soil induced by rainfall pulses in a desert ecosystem

Soil CO₂ flux is strongly influenced by precipitation in many ecosystem types, yet knowledge of the effects of precipitation on soil CO₂ flux in semi‐arid desert ecosystems remains insufficient, particularly for sandy soils. To address this, we investigated the response of sandy soil CO₂ flux to rainfall pulses in a desert ecosystem in northern China during August–September 2011. Significant changes (P < 0.05) were found in diel patterns of soil CO₂ flux induced by small (2.1 mm), moderate (12.4 mm) and large (19.7 mm) precipitation events. Further analysis indicated that rainfall pulses modified the response of soil CO₂ flux to soil temperature, including hysteresis between soil CO₂ flux and soil temperature, with F_s higher when T_s was increasing than when T_s was decreasing, and the linear relationship between them. Moreover, our results showed that rainfall could result in absorption of atmospheric CO₂ by soil, possibly owing to mass flow of CO₂ induced by a gradient of gas pressure between atmosphere and soil. After each precipitation event, soil CO₂ flux recovered exponentially to pre‐rainfall levels with time, with the recovery times exhibiting a positive correlation with precipitation amount. On the basis of the amounts of precipitation that occurred at our site during the measurement period (August–September), the accumulated rain‐induced carbon absorption evaluated for rainy days was 1.068 g C m^?2; this corresponds approximately to 0.5–2.1% of the net primary production of a typical desert ecosystem. Thus, our results suggest that rainfall pulses can strongly influence carbon fluxes in desert ecosystems. Copyright © 2014 John Wiley & Sons, Ltd.     

Estimation of land surface evapotranspiration over complex terrain based on multi‐spectral remote sensing data

Land surface evapotranspiration (ET) plays an important role in energy and water balances. ET can significantly affect the runoff yield of a basin and the available water resources in mountainous areas. The existing models to estimate ET are typically applicable to plains, and excessive data are required to calculate the surface fluxes accurately. This study established a simple and practical model capable of depicting the surface fluxes, while using relatively less parameters. Considering the complex terrain, solar radiation was corrected by importing a series of topographic factors. The water deficit index, a measure of land surface wetness, was calculated by applying the f_c (vegetation fractional cover)‐T_rad (land surface temperature) framework in the two‐source trapezoid model for evapotranspiration model to mountainous areas after corrections of temperature based on altitude variations. The model was successfully applied to the Kaidu River Basin, a basin with few gauges located in the east Tien Shan Mountains of China. Based on the time scale extensions, ET was analyzed at different time scales from 2000 to 2013. The results demonstrated that the corrected solar radiation and water deficit index were reasonably distributed in space and that this model is applicable to ungauged catchments, such as the Kaidu River Basin.     

Changes in Earth’s Energy Flows and Clouds in 228-Year Simulation with a High-Resolution AGCM

We have examined long-term changes in Earth’s energy flows at top of the atmosphere (TOA) and at Earth’s surface (land and ocean) by using 228-year simulation of a high-resolution global atmosphere model, MRI-AGCM3.2. It is found that the net downward short wave (SW) radiation (absorbed solar radiation, ASR) at TOA significantly increases during twenty-first century in agreement with a previous study. However, in the present study, the reason for the change is an increase in clear sky SW absorption by increased water vapor in the atmosphere, while it is a decrease in cloud amount in the previous study. It is also found that the long wave (LW) cloud radiative forcing for atmosphere is positive and increasing during twenty-first century in agreement with a previous study. The reason for the change in the present study is an increase in absorption by water vapor of the downward LW radiation emitted from clouds, while it is reductions of cloud amount in the middle troposphere in the previous study.     

On the water vapour absorption in the 8–13 μm spectral region for different atmospheric conditions

Experimental measurements of the absorption coefficient of atmospheric water vapour, at wavelengths between 8 and 13 m, are examined on the basis of atmospheric models describing the meteorological conditions of the observed atmospheres to obtain estimates of the foreign-broadening absorption coefficient for homogeneous paths. The results show that the variable contribution given by unresolved lines predominates on the continuum term due to wing effects of remote lines, even for rather high spectral resolutions.Associated with estimates of the self-broadening absorption coefficient, as proposed byRoberts et al. (1976), these data are applied to a wide set of atmospheric models, corresponding to various latitudes and seasons, indicating that the two components of the atmospheric absorption coefficient are closely related to the surface temperature. Because of the variable weights given by foreign- and self-broadening terms, the atmospheric absorption coefficient turns out to increase with the surface temperature with different rates at various wavelengths.     

Towards quantifying the negative feedback regulation of peatland evaporation to drought

The frequency and intensity of drought is projected to increase within the boreal region under future climatic conditions. Peatlands are widely considered to regulate water loss under drought conditions, increasing surface resistance (r_s) and reducing evaporative losses. This maintains peat moisture content, increasing the resilience of these globally important carbon stores. However, the magnitude and form of this important negative feedback response remains uncertain. To address this, we monitored the response of r_s to drought within four peat cores under controlled meteorological conditions. When the water‐table was dropped to a depth of 0.30 m and the humidity reduced to ≤40%, a step shift in r_s from ~50 s m^‐1 up to 1000 s m^‐1 was observed within burned and unburned peat, which virtually shuts down evaporation, limiting water loss. We show that measured near‐surface tension cannot be used to directly calculate this transition in peat surface resistance. However, empirical relationships that account for strong vertical variations in tension through the near‐surface and/or disequilibrium between pore air and near‐surface pore water pressure provide the potential to incorporate this negative feedback response into peatland ecohydrological models. Further observations are necessary to examine this response under dynamic atmospheric conditions. We suggest that the link between surface temperature and evaporation provides potential to further examine this feedback in either burned peatlands or peatlands with a low vascular vegetation cover. Copyright © 2013 John Wiley & Sons, Ltd.