Multiscale analysis of vegetation surface uxes: from seconds to years
|
|
- Christopher Richard
- 6 years ago
- Views:
Transcription
1 Advances in Water Resources ) 1119± Multiscale analysis of vegetation surface uxes: from seconds to years Gabriel Katul a, *, Chun-Ta Lai a, Karina Schafer a, Brani Vidakovic b, John Albertson c, David Ellsworth a, Ram Oren a a School of the Environment, Duke University, Box 90328, Durham, NC , USA b Institute of Statistics and Decision Sciences, Duke University, Durham, NC , USA c Department of Environmental Sciences, University of Virginia, Charlottesville, VA 22903, USA Received 5 June 2000; received in revised form 7 December 2000; accepted 19 March 2001 Abstract The variability in land surface heat H), water vapor LE), and CO 2 or net ecosystem exchange, NEE) uxes was investigated at scales ranging from fractions of seconds to years using eddy-covariance ux measurements above a pine forest. Because these uxes signi cantly vary at all these time scales and because large gaps in the record are unavoidable in such experiments, standard Fourier expansion methods for computing the spectral and cospectral statistical properties were not possible. Instead, orthonormal wavelet transformations OWT are proposed and used. The OWT are ideal at resolving process variability with respect to both scale and time and are able to isolate and remove the e ects of missing data or gaps) from spectral and cospectral calculations. Using the OWT spectra, we demonstrated unique aspects in three appropriate ranges of time scales: turbulent time scales fractions of seconds to minutes), meteorological time scales hour to weeks), and seasonal to interannual time scales corresponding to climate and vegetation dynamics. We have shown that: 1) existing turbulence theories describe the short time scales well, 2) coupled physiological and transport models e.g. CANVEG) reproduce the wavelet spectral characteristics of all three land surface uxes for meteorological time scales, and 3) seasonal dynamics in vegetation physiology and structure inject strong correlations between land surface uxes and forcing variables at monthly to seasonal time scales. The broad implications of this study center on the possibility of developing low-dimensional models of land surface water, energy, and carbon exchange. If the bulk of the ux variability is dominated by a narrow band or bands of modes, and these modes ``resonate'' with key state and forcing variables, then low-dimensional models may relate these forcing and state variables to NEE and LE. Ó 2001 Published by Elsevier Science Ltd. 1. Introduction Temporal variability in land surface sensible H) and latent LE) heat uxes and net ecosystem exchange NEE) occur over a wide range of scales ranging from fractions of seconds to decades. These uxes are interwoven at the most fundamental level, the leaf stomata. Describing the dynamics of H, LE, and NEE requires resolving an extensive high-dimensional spectrum of variability [2]. To date, no one model is able to reproduce the variability in these uxes across such spectrum of scales. Over seconds, biosphere±atmosphere exchange of water vapor q,heat T a, and carbon dioxide C) are governed by complex turbulent eddy motion, rich in spectral properties. At hours to monthly scales, much of the variability is induced by interactions between * Corresponding author. Tel.: ; fax: address: gaby@duke.edu G. Katul). weather patterns, bulk turbulent ow characteristics, eco-physiological and biochemical characteristics of the canopy, and the hydrologic state of the ecosystem. The latter variable is strongly subject to biological feedbacks by means of stomatal closure and xylem cavitation [40]. At monthly to annual time scales, synoptic weather patterns superimposed on seasonal variations are responsible for changes in ecosystem properties and, consequently, dominate the variability in these land surface uxes. This multiscale nature of land surface uxes is now motivating the development of a new generation of canopy±atmosphere models accompanied by long-term measurements [11]. The objective of this study is to explore structural aspects and controls on the temporal spectra of H, LE, and NEE from fractions of seconds to years. A case study is made using measurements of H, LE, and NEE collected over three years at 10 Hz above a 17 year old loblolly pine stand at Duke Forest, near Durham, North Carolina. The Duke Forest site is part of a long-term NEE monitoring initiative /01/$ - see front matter Ó 2001 Published by Elsevier Science Ltd. PII: S )00029-X
2 1120 G. Katul et al. / Advances in Water Resources ) 1119±1132 aimed at representing the ecological behaviour of managed pine plantation in the Southeastern US. About 70% of the total forested area in the Southeastern US 35: Ha is pines, and 41% of this area is pine plantations younger than 25 years [33]. The site index ˆ 16 m at 25 years) of the Duke Forest pine plantation also represents the mode of managed pine plantations in the Southeast. 2. Methods of analysis We propose multiscale methods that can analyze the spectral and cospectral properties of H, LE and NEE. Traditional spectral and cospectral analysis in the Fourier domain is primarily complicated by random gaps in the ux time series measurements and the convolution between these uxes and the suite of hydrologic and environmental variables at multiple scales. Wavelet transformations can localize the gaps in the ux record and eliminate them from the spectral and cospectral calculations. Additionally, the multiscale interactions between the uxes and forcing variables are naturally resolved in the standard multiresolution formulation of OWT [13,14,30,31,42,43,47]. Towards this end, all spectral and cospectral properties will be analyzed using OWT, as described in the following section Fast wavelet transform FWT) The Haar wavelet basis is selected for: i) its temporal di erencing characteristics, ii) its locality in the time domain, iii) its short support which eliminates any edge e ects in the transformed series, and iv) its wide use in atmospheric turbulence research e.g. see review in [18]), and v) its computational e ciency. The Haar wavelet coe cients d and coarse grained signal or scaling function coe cients) S m can be calculated using the usual multiresolution decomposition on a measurement vector f x j originally stored in S 0 for m ˆ 0 to M 1, where S 0 ˆ 0 and M ˆ log 2 N is the total number of scales. Throughout, n is time averaging any variable n over intervals ranging from 20 to 30 min. The decomposition produces N 1 wavelet coe cients de- ning the Haar OWT of f x j [26,27] Wavelet spectra and cospectra The N 1 discrete Haar wavelet coe cients satisfy the energy conservation analogous to Parseval's identity) X N 1 jˆ0 f j 2 ˆ XM mˆ1 2 M m X 1 iˆ0 d m iš 2 ; 1 where m is the scale index and i the time or space) index. The total energy T E contained in a scale R m is computed from the sum of the squared wavelet coef- cients in Eq. 1) recall S 0 ˆ0) at a scale index m using T E R m ˆ 1 N 2 M m X 1 iˆ0 d m iš 2 ; where R m ˆ 2 m dy with dy ˆ fs 1 ) is the time scale, and f s is the sampling frequency. Hence, the power spectral density function P R m is computed by dividing T E in Eq. 2) by the change in frequency DR m ˆ 2 m dy 1 ln 2 so that P R ˆhh d m iš 2 iidy m ; 3 ln 2 where hh ii in Eq. 3) is averaging overall values of the position index iš at scale index m. Notice that the wavelet power spectrum at R m is directly proportional to the average of the squared wavelet coe cients at that scale. Because the power at R m is determined by averaging many squared wavelet coe cients as evidenced in Eq. 3), the wavelet power spectrum is generally smoother than its Fourier counterpart and does not require windowing and tapering. Such averaging in the wavelet domain is statistically stable because of the whitening properties of wavelets later discussed in greater details). The Haar wavelet cospectrum of two variables e.g. w 0 and c 0, vertical velocity and CO 2 uctuations, respectively) can be calculated from their respective wavelet coe cient vectors d w m and d c m at scale R m using P wc R ˆhhd m w m where w 0 c 0 ˆ 1 X M N mˆ1 iš d m c ln 2 2 M m X 1 iˆ0 2 iš iidy ; 4 d m w išd m c iš : 5 While w 0 and c 0 are turbulent uctuations of vertical velocity and CO 2 concentration, with w 0 ˆ c 0 ˆ 0 see [18,19,21]), the above formulations in Eqs. 4) and 5) are valid for all three scalars Robustness to frequency shifts and gaps in time series To illustrate the robustness of OWT to time±frequency shifts and to the presence of gaps in the time series, we consider two synthetic experiments. The rst experiment considers a frequency change within the time series and contrast its detection by Fourier spectra and wavelet half-plane imaging. The second experiment considers a fractional Brownian motion fbm) time series, randomly infected by large gaps typical in longterm ux monitoring), and contrasts the wavelet spectra
3 G. Katul et al. / Advances in Water Resources ) 1119± for the original and gap-perturbed synthetic signals. For the rst experiment, consider the two signals y 1 t and y 2 t shown in Fig. 1. These two signals are constructed as follows: y 1 t ˆ0:5 sin 4t 0:5sin 32t and y 2 t ˆ sin 4t ; t 6 p; sin 32t ; t > p: Notice in Fig. 1 that the Fourier power spectra for these two signals are nearly identical suggesting that Fourier spectra alone are not well suited to discern these frequency shifts. The squared wavelet amplitudes are shown in the wavelet time±frequency plane in Fig. 1. In contrast to the Fourier power spectrum, the energy is resolved with respect to both scale or frequency) and temporal extent in the wavelet half-plane. This is an important quality for studying geophysical phenomenon that possess temporal variability in the scale-wise distribution of energy for example, studies of interdecadal drought variability). An additional challenge in analyzing long-term ux measurements is the frequent occurrence of random gaps in the measurement record. In fact, many long-term ux monitoring sites within the AmeriFlux and EuroFlux programs have a data recovery rate below 70% and gap events whose duration exceeds 10% of the record size [6]. The occurrence of large gaps precludes the use of traditional Fourier type power spectra. Wavelets, and OWT in particular, o er distinct advantages here. The wavelet coe cients are localized in space, and orthogonal in their relationship to each other, thus limiting the e ects of a record gap to the coe cients local to the gap. In computing the global spectra the a ected coe cients are simply omitted from the averaging operator. To illustrate, consider a class of non-stationary stochastic processes, known as fractional Brownian motion fbm) whose statistical properties in the Fourier and wavelet domains are well known see Appendix A for a brief review). Fig. 2 a) shows an fbm time series with a Hurst exponent H e ˆ 1=3. At random, three unevenly spaced, large gaps are injected into this time series Fig. 2 b)). The gaps collectively represent 30% of the sampling duration for this experiment a representative value for the Duke Forest AmeriFlux site). The wavelet image of these two time series is also shown in Fig. 2 c) and d). Notice that the gaps are well localized and can be readily detected in the wavelet half-plane. The wavelet power spectra of these two fbm series are compared in Fig. 3 with excellent agreement between the original and gapinfected time series. In the gap-infected time series, the hh ii in the power spectrum de nition is applied only to the non-gap-infected coe cients at scale index m. Fig. 1. a) Synthetic time series with and without time±frequency shifts; b) standard Fourier power spectra for the two series; c) squared wavelet coe cients in the time±frequency half-plane. Fig. 2. The e ects of gaps on wavelet spectral properties: a) a standard fbm 1/3) time series; b) a standard fbm 1/3) time series with 30% gaps in the record; c) and d) the wavelet coe cients in the time± frequency half-plane.
4 1122 G. Katul et al. / Advances in Water Resources ) 1119±1132 Fig. 3. Comparison between wavelet spectra for an fbm 1/3) signal with gaps and no gaps. For reference, the standard Fourier power spectrum for the original fbm 1/3) time series and the theoretical power-law ˆ )5/3) associated with a Hurst exponent H e ˆ 1=3 are also shown Whitening property of wavelet transformations An additional desirable feature of wavelets in dealing with time series measurements is their decorrelation property [8,43,47]. Informally speaking, when correlated measurements in the time domain are transformed to wavelet domain, the wavelet coe cients exhibit little or no autocorrelation. Because of this property, wavelets are often called approximate Karhunen±Loeve transformations. In terms of a random process X t, the Karhunen±Loeve transformation is P i Z i/ i t, where Z i are uncorrelated random variables and / i t are functions from an orthonormal system, see [43]. To illustrate the decorrelation property of wavelets, we consider two standard test functions: blocks and doppler, for a sample size N ˆ The two time series are shown in Fig. 4 top panels). These two test-signals represent frequency shifts doppler) and rapid transients blocks) in time. The lower panels show the autocorrelation function for both time series dotted) as well as the autocorrelation function for their Haar wavelet coe cients at the nest scale or level). Notice that the autocorrelation function of the Haar wavelet coe cients decays very rapidly with lag for both signals when compared to the autocorrelation function in the time domain. This rapid decorrelation amongst coe cients is attributed to the so-called whitening property of wavelet transformations. This property is desirable when computing Haar spectral and cospectral functions in Eqs. 3) and 4) because uncorrelated wavelet coe cients provide more stable averages or statistics) when contrasted to correlated and long-range dependent coe cients. Having discussed the methods of analysis, the experimental setup for measuring H, LE, and NEE as well as the physiological variables and forcing functions is brie y described. Fig. 4. Two standard test-signals of Donoho±Johnstone Doppler and Blocks). The autocorrelation q s as a function of lag s is shown for the time series dotted) and their wavelet coe cients solid) at the nest level. 3. Experiment 3.1. Study site Hydrometeorological and eco-physiological measurements were conducted at the Blackwood Division of the Duke Forest near Durham, North Carolina N; W, elevation ˆ 163 m). These measurements, initiated in August of 1997 and continuing to the present date, are part of the AmeriFlux) long-term ux monitoring initiative. Here, we focus on a long-term period from mid-august 1997 to August The site is a loblolly pine Pinus taeda L.) forest, planted in It extends 300±600 m in the east±west direction and 1000 m in the north±south direction. The mean canopy height h c was :5 m) in The topographic variations are small terrain slope changes < 5%) so that the in uence of topography on the turbulent transport can be neglected [12] Eddy-covariance measurements The turbulent uxes of CO 2, water vapor, and heat) between the land and the atmosphere were measured by an eddy-covariance system comprised of a CO 2 =H 2 O infrared gas analyzer LICOR-6262, LI-COR, LINCOLN, NE, USA), a triaxial sonic anemometer CSAT3, CAMPBELL SCIENTIFIC, LOGAN, UT, USA), and a krypton hygrometer KH 2 O, Campbell Scienti c). The anemometer and hygrometer were positioned 15 m above the ground surface and anchored on a horizontal bar extending 1.5 m away from the walkup tower. The hygrometer was used to assess and support corrections of) the tube attenuation e ects in the infrared gas analyzer. The open path hygrometer also provided a reference against which the lag time of the closed path
5 G. Katul et al. / Advances in Water Resources ) 1119± Fig. 5. Eddy-covariance measured time series of uxes H, LE, and NEE as well as the spatially averaged soil moisture content h in the top 30 cm and net radiation R n from August 1997 to August The horizontal line in the moisture content time series is the critical at which further reduction in h induces bulk canopy stomatal closure after [36]). infrared gas analyzer could be measured, as discussed in [16,17]. Analog signals from these instruments were sampled at 10 Hz using a Campbell Scienti c 21X data logger. All the digitized signals were transferred to a portable computer via an optically isolated RS232 interface for future processing. Raw 10 Hz measurements were processed to de ne 30 min average turbulent uxes using the procedures described in [16,17,23]. The approximate data recovery rate of the eddy-covariance measurements, shown in Fig. 5, was 65% with downtime due to storms, hurricanes, pump failure, non-optimal wind directions, loss of electrical power, on-site calibration of gas analyzer, and periodic maintenance and factory calibration of instruments. The annual NEE and LE for this stand vary from 580 to 670 g C m 2 yr 1 and from 520 to 580 mm yr 1 using a 12 month moving average) Other meteorological variables In addition to the ux measurements, a T a =RH probe HMP35C, CAMPBELL SCIENTIFIC) was positioned at 15.5 m to measure the mean air temperature T and relative humidity. A Fritchen-type net radiometer and a quantum sensor Q7 and LI-190SA, respectively, LI- COR) were installed to measure net radiation R n and PAR, respectively. All meteorological variables were sampled at 1 s and averaged every 30 min using a 21X Campbell Scienti c datalogger Hydrologic measurements Precipitation was measured by a tipping bucket gage TEXAS ELECTRONICS, MODEL 525, DALLAS, TX) above the canopy every 1 s and integrated every 30 min. The gage has a 476 cm 2 circular opening catchment area. Soil moisture content h was measured at 24 locations by an array of Campbell Scienti c CS 615 soil moisture content vertical rods, each 30 cm in length. The mean of all 24 rods is also shown in Fig. 5. The calibration of these rods is described in [35]. We note that when soil moisture drops below 0.2, cavitation in the xylem-root system is likely to occur thereby reducing the stomatal conductance [29,35]. Such low soil moisture content regimes in the root zone coincide with some of the extended summer droughts recently experienced in the Southeast. In fact, based on our three year record at the site, the time fraction in which the root-zone experienced such low h or drought) is about 16% as evidenced from Fig CO 2 =H 2 O pro les within the canopy A multi-port system was installed to measure the mean CO 2 C and H 2 O q concentration inside the canopy at 10 levels 0.1, 0.5, 1.5, 3.5, 5.5, 7.5, 9.5, 11.5, 13.5 and 15.5 m). Each level was sampled for 1 min 45 s sampling and 15 s purging) at the beginning, the middle, and the end of a 30 min sampling duration. The ow rate within the tube internal diameter ˆ 1=6 00 ) was 0:9 l min 1 to dampen all turbulent uctuations. The modeled concentration pro les within the canopy volume were compared to these measurements every 30 min to assess the accuracy of the modeled source distribution and dispersion calculations. These measurements were also used to investigate the storage ux terms within the canopy volume Leaf-level physiological parameters The leaf-level physiological parameters were determined from gas exchange measurements by a portable infra-red gas analyzer system for CO 2 and H 2 O CIRAS-1, PP-SYSTEMS, HERTS., U.K.) operated in open ow mode with a 5.5 cm long leaf chamber and an integrated CO 2 gas supply system. The chamber was modi ed with an attached Peltier cooling system to maintain chamber temperature near ambient atmospheric temperature. The data were collected for upper canopy foliage at heights 95% of the total tree height, accessed with a system of towers and mobile, vertically telescoping lifts [5]. These measurements were collected over a broad range of environmental and hydrologic conditions spanning a period of two and a half years May 1997±October 1999) and were used to de ne values of all the canopy physiological parameters included in a one-dimensional vegetation±atmosphere coupled photosynthesis±turbulence model hereafter referred to as CANVEG).
6 1124 G. Katul et al. / Advances in Water Resources ) 1119± Plant area density The vertical distribution of the plant area density PAI) was measured by gap fraction techniques following the theory presented in [34]. A pair of optical sensors with hemispherical lenses LAI-2000, LI-COR) was used for canopy light transmittance measurements from which gap fraction and plant area densities were calculated. The measurements were made at 1 m intervals from the top of the canopy to 1 m above the ground to produce the vertical pro le in plant area index. PAI measurements were made within two weeks of each of the measurement periods at peak PAI, from the same tower used for measuring the ow statistics, CO 2 concentration pro les and turbulent uxes. 4. Results and discussion The H, LE, NEE, h, and R n time series, shown in Fig. 5, clearly exhibit rich variability at multiple scales. When discussing the spectral properties of land surface uxes it is helpful to consider three fundamental time scales: fractions of seconds to minutes turbulent time scales), hours to tens of days meteorological time scales), and months to years seasonal time scales). At the very ne time scales, much of the ux variability is induced by turbulent motion whose production is partly governed by the existing mean meteorological conditions. At the hour to days time scale, the interaction between the physiological and biophysical characteristics of the canopy and its micro-climate become the dominant modes of ux variability. Over these ne- and mid-range time scales, the vegetation and its physiological characteristics are assumed, in a rst-order analysis, to be static. At longer time scales, seasonal patterns and vegetation dynamics phenology) share importance with interannual climate variability. The interaction between available energy, hydrologic status, and the vegetation function control the critical modes of temporal variability in heat, water vapor and CO 2 uxes. Below, we present a framework for representing the processes at their respective time scale Turbulent time scales At time scales shorter than an hour, the processes governing the land surface uxes are complex eddy motion. Spectrally, the eddy motion is commonly decomposed into three regimes: production, inertial cascade, and dissipation. The production time scales are the scales in which turbulence extracts energy from the meteorological forcing, with in uences by canopy morphology and mean meteorological conditions, particularly the mean longitudinal velocity U. Scale-wise energy distribution in this range can be well described by h c, the friction velocity u ˆ u 0 w 0 1=2, and the characteristic scalar concentration uctuation scale c ˆ w 0 u 0 =u [38]. Near the canopy±atmosphere interface i.e. z=h c ˆ 1), the friction velocity a variable related to the mechanical production of turbulence) is strongly correlated with U a variable related to the mean meteorological forcing) with U=u 3:3 [19,38] for near-neutral conditions. In the vicinity of z=h c ˆ 1 for rough canopies, the ow is predominantly nearneutral e.g., >70% of the time for this study). Recent advances in hydrodynamic stability theory suggested that much of the organized motion near the canopy± atmosphere interface, through which the turbulence is produced, is attributed to Kelvin±Helmholtz type instability which scales with h c and u. In fact, such theory explained well the spectral peaks in scalar uxes e.g. [19]) via a mixing layer analogy. The latter analogy challenged the common paradigm that the structure of turbulence near the canopy±atmosphere interface resembles a rough-wall boundary layer. The role of this organized eddy motion on mass and heat transfer cannot be overemphasized, with more than 90% of the net transport occurring in 10% of the time fraction [21,23]. In the mid-range of turbulent eddy scales we encounter the onset of an inertial cascade, described well by Kolmogorov [9,15,25] type scaling or K41). This scaling is evident in many scalar spectra and cospectra, with P wc proportional to Rm 7=3 as in Fig. 6, with the exception of w 0 T 0, for which the power-laws reported near the canopy±atmosphere interface vary from Rm 5=3 to Rm 7=3. The departure from Rm 7=3 has been attributed to the active role of T 0 in buoyant production or destruction of turbulent kinetic energy, active at these ne scales. At yet smaller scales, the action of molecular viscosity becomes dominant and the turbulent uctuations are dissipated away as heat. These small scales are in fact ner than the nest scales 10 Hz) measured in this study Meteorological time scales In this range of time scales, the variations in land surface uxes are due to the strong and often non-linear interaction between the physiological, radiative, biochemical, and drag properties of the canopy and its microclimate. To represent such interactions, we developed a onedimensional multilevel canopy-vegetation model that resolves the key processes in such interaction. The model, now known as CANVEG after [1], is described in [28] but is brie y reviewed below. In its primitive form, the model considers three coupled one-dimensional scalar conservation equations at a level within the canopy volume of thickness dz, in the absence of advective transport:
7 G. Katul et al. / Advances in Water Resources ) 1119± (a) (b) (c) Fig. 6. Measured Haar wavelet ux cospectra for w 0 T 0 a), w 0 q 0 b), w 0 c 0 c), as a function of frequency Hz) for 41 half-hour runs. For comparisons, the w 0, T 0, q 0, and c 0 time series are normalized to zero-mean and unit variance. The 7=3 power-law, as derived from K41's locally homogeneous and isotropic turbulence, is also shown. ot ot ˆ of T oz S T ; oq ot ˆ of q oz S q; 6 oc ot ˆ of c oz S c; where, F T ˆ w 0 T 0, F q ˆ w 0 q 0,andF c ˆ w 0 c 0 are the turbulent uxes at a level z and time t; S T ; S q, and S c are the heat, water vapor, and CO 2 sources S ˆ source, S ˆ sink) within a canopy layer of thickness dz. At z=h c > 1, q C p F T ˆ H, L v F q ˆ LE, and F c ˆ NEE, where L v is the latent heat of vaporization, q is the mean air density, and C p is the speci c heat capacity of dry air under constant pressure. Eq. 6) contains nine unknowns. Using Lagrangian uid mechanics [41], it is possible to establish a relationship between the sources and mean concentrations via Z Z T z; t ˆ P z; t jz 0 ; t 0 S T z 0 ; t 0 dz 0 dt 0 ; Z Z q z; t ˆ P z; t jz 0 ; t 0 S q z 0 ; t 0 dz 0 dt 0 ; 7 Z Z C z; t ˆ P z; t jz 0 ; t 0 S c z 0 ; t 0 dz 0 dt 0 ; where P z; t jz 0 ; t 0 is the transition probability density function that can be estimated from the velocity statistics within the canopy volume via a dispersion matrix as discussed in [37]. Higher-order Eulerian closure models for estimating these velocity statistics within the canopy from the drag and foliage distribution are described elsewhere e.g. [20,24]) and are used to compute P z; t jz 0 ; t 0 or the dispersion matrix. The boundary conditions used in modeling the ow statistics inside the canopy needed to compute P z; t jz 0 ; t 0 are presented in Table 1. In short, Lagrangian uid mechanics principles provide three additional equations to the set of three equations in 6). The Lagrangian frame of reference or Eq. 7) does not uniquely de ne such relationships and an equivalent formulation in the Eulerian frame of reference is also possible e.g. [22]). The later formulation is based on higher-order turbulence closure principles for scalar transport. To mathematically ``close'' this problem de ned by Eqs. 6) and 7), three additional equations are needed and are derived by assuming the transfer of scalars from the leaf to the atmosphere is governed by Fickian diffusion through the stomatal cavities and within a thin boundary layer at the leaf surface. This approximation leads to: S T z; t ˆa z T i z; t T z; t ; r b z; t S q z; t ˆa z q i z; t q z; t r b z; t r s z; t ; 8 S c z; t ˆa z C i z; t C z; t r b z; t r s z; t ; where r b and r s are the boundary layer and stomatal resistances, respectively, and T i ; q i ; C i are the stomatalcavity temperature, water vapor and carbon dioxide concentrations, respectively. Hence, while the Fickian di usion approximation in Eq. 8) provided the necessary equations to close the primitive equations, it introduces ve additional unknowns: r b ; r s ; T i ; q i ; C i.in the CANVEG approach in [28], these unknowns are determined from the leaf-level energy balance at each layer within the canopy) and radiation attenuation model of Norman [3], the photosynthesis model of Farquhar [7], the boundary layer resistance r b model of [39], and the stomatal conductance 1=r s model of Collatz [4], given by: 1 ˆ m RH A b; 9 r s C s where, m and b are species-speci c parameters, A is the carbon assimilation rate which depends on S c and leaf area density a z related to LAI via LAI ˆ R h a z dz), 0 and RH z; t and C s are the mean leaf relative humidity and CO 2 concentration at time t and level z. The coupling between the radiative transfer, energy balance, and
8 1126 G. Katul et al. / Advances in Water Resources ) 1119±1132 Table 1 Key eco-physiological, radiative, and drag properties of the pine forest used in the CANVEG model Variable Value Comments Leaf area index August, :52 m 2 m 2 Measured using LAI-2000 July, :31 m 2 m 2 Measured using LAI-2000 Canopy height h c 14 m Measured in 1999 Characteristic length for calculating r b m Measured Clumping factor 0.8 Estimated in [28] Maximum Rubisco capacity V cmax 64 lmol m 2 leaf s 1 Estimated from porometry Ball±Berry slope parameter m 5.9 Estimated from porometry Ball±Berry intercept parameter b Estimated from porometry Maximum quantum e ciency 0:08 mol mol 1 Assumed as in [28] Leaf absorptivity for PAR 0.83 Assumed as in [28] Drag coe cient 0.2 Estimated from [20] At z ˆ 1:15h c r u =u 2.4 See [28] r v =u 1.9 See [28] r w =u 1.25 See [28] d U=dz u =k z d See [28] At z ˆ 0 dr u =dz 0 See [28] dr v =dz 0 See [28] dr w =dz 0 See [28] U 0 See [28] The boundary conditions for longitudinal, lateral, and vertical velocity standard deviations r u, r v, r w, respectively) used in the second-order closure model are also shown. The mean wind speed U boundary conditions are also presented, with d representing the zero-plane displacement, iteratively calculated by the closure model, and k ˆ 0:4 is von Karman's constant. the Farquhar photosynthesis model requires all three scalars to be simultaneously considered in Eqs. 6)± 8) and at each canopy layer. These equations are solved for F T z; t, F q z; t, F c z; t, T z; t, q z; t, C z; t, S T z; t, S q z; t, S c z; t, r s z; t, r b z; t, T i, q i, and C i for speci ed mean wind speed U, mean air relative humidity, mean air temperature, mean CO 2 concentration, and PAR in the free space above the canopy every 30 min. The parameters for the Farquhar photosynthesis model, Eq. 9), and the radiation/energy balance are shown in Table 1. The details of implementing these radiative, energy, and eco-physiological schemes in Eqs. 6)± 9) are discussed in [28,29]. While several multilayer radiatiative-photosynthesis models have been proposed e.g. [32,44,45]), the present CANVEG model has the added advantage of fully integrating the turbulent transport dynamics and the microclimate to the sources and sinks via Eqs. 6) and 7). Hence, canopy morphology via a z ) a ects the radiation attenuation, the ow statistics within the canopy, the energy balance, and the scalar source strength at each level in a consistent manner. To illustrate the CANVEG model performance, a comparison between modeled and measured mean CO 2 concentration inside the canopy volume is shown in Fig. 7. The over-all features of CO 2 buildup during night-time and the well-mixed features during day-time are well reproduced by the CANVEG model. Using a similar approach but with prescribed within-canopy velocity statistics, other studies e.g. [10]) also report good agreement between measured and modeled CO 2 concentration pro les for a boreal forest. We restate that many of the coupled photosynthesis±radiative transfer schemes e.g. [32,44,45]) cannot predict the space±time evolution of the mean CO 2 concentration, temperature, and water vapor concentration, but rather require such concentration measurements as input. A comparison between eddy-covariance measured and modeled values of H, LE, and NEE, on a 30 min time step, are shown in Fig. 8 for the same periods as in Fig. 7 with good agreement between measured and modeled uxes. To assess the spectral recovery of measured H, LE, and NEE by CANVEG, the measured and modeled spectra in Fig. 9 are also compared, with reasonable agreement between these spectra for a wide range of time scales 1 h to 21 days). We restate that in the above calculations, the vegetation was assumed static i.e. a z and all the physiological parameters of the Collatz [4] model, drag properties, and radiative characteristics including foliage clumping are not permitted to vary in time). If the temporal evolution of these canopy physical and physiological properties are known or speci ed, then the CANVEG model can be used to reproduce the full spectrum of H, LE, and NEE from hours to years Seasonal to annual time scales While the CANVEG approach can be integrated to seasonal time scales, the temporal dynamics of the
9 G. Katul et al. / Advances in Water Resources ) 1119± Fig. 7. Comparison between measured and CANVEG modeled mean CO 2 concentration for two periods: August 1999 a) and July 2000 b). physiological parameters required by such an approach are rarely measured at such temporal resolution [5]. For example, one of the key physiological parameters needed in the Collatz [4] model adopted in our CAN- VEG is the maximum Rubisco capacity per unit leaf area V cmax. For pines, V cmax depends on leaf nitrogen concentration, which varies at monthly to annual time scales. Furthermore, the leaf area density pro le as well as LAI) also varies at seasonal time scales thereby altering the radiation attenuation, clumping, and the fraction of sunlit foliage throughout the canopy. Additionally, a z dynamics commonly result in phenological and physiological changes with over-wintering foliage having di erent m and V cmax when compared to newly grown foliage. In short, the land surface ux variability at monthly and longer time scales is attributed to changes in forcing and meteorological variables e.g. R n ), hydrologic state e.g. h), eco-physiological properties e.g. V cmax ; m), and ecological properties e.g. a z all convolved at multiple time scales. While such
10 1128 G. Katul et al. / Advances in Water Resources ) 1119±1132 Fig. 8. Comparison in the time domain between eddy-covariance measured and CANVEG modeled H, LE, and NEE time series for the two periods: August 1999 a) and July 2000 b). physiological changes are likely to pose major challenges to modeling water and carbon uxes at longer time scales, a recent study on six European coniferous forests suggests that long-term water and carbon uxes can be modeled without using detailed physiological and site speci c information [46]. In [46], arti cial neural networks were used to select a minimal set of input variables to model LE and NEE without resorting to detailed physiological information. The appeal to a neural network type analysis is to empirically derive highly non-linear relationships between forcing variables and uxes. Here, rather than producing such empirical non-linear relationships between forcing variables and uxes, we explore the characteristic time scales which exhibit the largest variability in uxes and forcing variables as well as the scales which exhibit strongest wavelet spectral correlations. For this reason, we investigate both measured wavelet spectral and cospectral characteristics of the land surface uxes with respect to R n energy forcing) and soil moisture h hydrologic state) in Figs. 10 and 11. In Fig. 10, it is clear that at hourly time scales, the three scalar uxes exhibit larger variability than R n and h particularly NEE). More important from Fig. 10 is that the seasonal variability for all the variables is at least one order of magnitude more energetic or important) than the diurnal or daily variability. That is, when annual carbon or water budgets are required, it is the seasonal variability that is critical to resolve. We restate that this pine forest is evergreen and retains foliage and water ux as well as carbon uptake) yearround. In all the cospectral calculations, the interactions amongst the variables are strongest for seasonal 200 d) time scales. For all three land surface uxes, the seasonal time scales are at least one order of magnitude more im- Fig. 9. Wavelet spectral comparison between eddy-covariance measured and CANVEG modeled H, LE, and NEE for the time series shown in Fig. 8. Fig. 10. Measured Haar wavelet ux spectra for H, LE, NEE, R n, and h time series in Fig. 5. For spectral intercomparisons, all the time series are normalized to zero mean and unit variance.
11 G. Katul et al. / Advances in Water Resources ) 1119± Fig. 11. Measured Haar wavelet ux cospectra for the time series in Fig. 5. For comparisons, the cospectra are normalized by covariances between signals in the Fig. 5 caption. portant than the diurnal or daily time scales for the correlations considered. Interestingly, time scales at which the variability in h ``resonates'' most with variability in LE are quite localized to seasonal time scales despite several drought events i.e. h < 0:2 that induced-stomatal closure at shorter time scales e.g., Fig. 5). The correlation at the annual scale is a hydrologic interaction induced mainly by the combined drought and high LE in the summer and the high soil moisture and low LE in the winter despite an even precipitation distribution at monthly time scales). Note that the sample size of the soil moisture record is different from the LE record e.g. Fig. 5). Such a discrepancy in sampling sizes prevent the use of Fourier-type cospectral calculations, and further lends support to the use of OWT for such long-term ux analysis. To further explore the relationship between R n and LE, we consider their correlation at 30 min and 30 days averaging intervals in Fig. 12. While the relationship between LE and R n is poor on 30 min time scale, the opposite is true at monthly or longer) time scales. This increased correlation between R n and LE, may appear obvious at rst glance, as summer months have higher maximum daily net radiation, longer day-lengths, and higher LAI when compared to winter months. Less apparent in this strong correlation between R n and LE is that the temporal dynamics of LAI and R n are not synchronous as shown in Fig. 12. Note that the rapid change in LAI occurs over one to two months with little change thereafter. On the other hand, monthly R n exhibits variability amongst all months. In other words, while monthly LE is strongly related to monthly R n, monthly latent heat ux per unit leaf area can be large or small for a given R n as evidenced from Fig. 12. The time scale of foliage dynamics may be further ampli ed for deciduous forests in highly seasonal latitudes. Another Fig. 12. Relationship between R n and LE at multiple time scales: a) the relationship for the 30 min time scale; b) variations of LE per unit leaf area at monthly time scales with R n at the canopy-top; c) monthly time scales with the solid line representing equilibrium evaporation; d) monthly variation of LAI with R n. interesting point is that for R n > 120 W m 2 the ole=or n is well described by equilibrium evaporation D= D c, where D is the slope of the saturation vapor pressure±temperature curve, and c is the psychrometric constant. However, for monthly R n < 80 W m 2 i.e. winter months), we note that ole=or n < D= D c coinciding with low foliage and temperature regulation on stomata. For pines, stomatal closure occurs when the air temperature drops below 10 C [5]. Over the entire year, we found that LE < D= D c R n suggesting that the decrease in LE ux by stomata exceeds the net increase in LE because of nite atmospheric drying power. Finally, we note that the increased correlation between LE and R n with increasing time scales is not restricted to water vapor exchange. We also found similar patterns for NEE and H as evidenced in Fig. 13. Fig. 13. Variations of sensible heat H a) and net ecosystem exchange NEE b) on R n for monthly time scales.
12 1130 G. Katul et al. / Advances in Water Resources ) 1119± Conclusions This study, the rst to explore the measured wavelet spectra and cospectra of land surface uxes from fractions of seconds to three years, demonstrated the following: Orthonormal wavelet transformation provides a robust framework for analyzing the spectral and cospectral properties of long-term ux records that exhibit 1) frequency shifts in time, and 2) multiple gaps or missing data. The wavelet spectra of the three land surface uxes carbon, water, and heat) demonstrated that the largest variability occurs at seasonal time scales. In fact, the energy or activity at seasonal time scales is at least one order of magnitude larger than variability at diurnal 12 h) time scales. As expected, the measured soil moisture spectrum exhibits large activities at seasonal time scales with little activities at sub-daily scales. Three broad categories of time scale models must be considered when describing temporal dynamics of land surface uxes: 1) turbulent time scales seconds to an hour), 2) meteorological time scales hours to days), in which physiological properties, radiative transfer, and diurnal boundary layer dynamics interact, and 3) seasonal time scales weeks to year). The wavelet cospectral comparison suggests that for short turbulent time scales, the Kolmogorov theory K41 describes well observed spectral power-laws. For longer time scales characteristic of turbulent production 100±1000 s), the mixing layer analogy of Raupach [38] reproduces well the observed cospectral peaks in the turbulent ux of water and carbon. The CANVEG model reproduced well the H, LE, and NEE spectra over a broad intermediate range of time scales hours to days), as long as the vegetation characteristics can be assumed static. We showed that on seasonal time scales, the interaction between radiative forcing and land surface uxes is substantially enhanced when compared to the 30 min time scales consistent with a recent study on six European conifer sites. This enhancement is, in part, attributed to the intricate balance between the ux per unit leaf area, leaf area dynamics, and the seasonal dynamics of R n. The broad implications of this study center on the possibility of developing low-dimensional models of land surface water, energy, and carbon exchange. If the bulk of the ux variability is dominated by a narrow band or bands of modes, and these modes ``resonate'' with key state and forcing variables, then low-dimensional models may relate these forcing and state variables to NEE and LE. In fact, current use of remotely sensed data to infer the land surface uxes in combination with low-dimensional models bene t from such an analysis in two ways: 1) it identi es the time scales at which the remotely sensed observations are likely to ``resonate'' with the uxes and hence the time scales at which remotely sensed observations must be resolved, and 2) it identi es the level of predictability that these models can o er given the observed degree of spectral correlation between these observations and the scalar uxes. Acknowledgements The authors would like to thank Yavor Parashkevov and Cheng-I Hsieh for their assistance in the data collection and Judd Edburn for his overall help at the Duke Forest. This research was supported by the Department of Energy DOE) Na- tional Institute for Global Environmental Change NIGEC) through the NIGEC Southeast Regional Center at the University of Alabama in Tuscaloosa DE-FC03-90ER61010) and the FACE-FACTSproject, the National Science Foundation NSF) Grants DMS , EAR , and BIR at Duke University. The MATLAB programs can be obtained from the authors on request. The Duke Forest data sets can be accessed from /facts-1/facts.html. Appendix A. Properties of fractional Brownian motion The fbm is a Gaussian, zero mean, continuous, nonstationary process, indexed by the parameter H e Hurst exponent, 0 < H e < 1) such that B H 0 ˆ0; and B H t r B H t N 0; r 2 H j r j2he ; where r 2 H ˆ C 1 2H e cos ph e ph e and C x ˆR 1 0 t x 1 e t dt is the C function. From its de nition, the sample paths of fbm satisfy the scaling equation B H at ḓ a H B H t, and hjb H t r B H t j n iˆc fbm jrj nhe ; where the C fbm is r 2 H 2 n=2 C n 1=2 C 1=2 and where ḓ is equal in distribution [8]. The Fourier transformation of the scaling equation is proportional to 1=jxj n He 1=2 with an average power spectrum proportional to 1=jxj 2He 1.
13 G. Katul et al. / Advances in Water Resources ) 1119± References [1] Baldocchi D, Meyers T. On using micrometeorological and biogeochemical theory to evaluate carbon dioxide, water vapor and trace gas uxes over vegetation: a perspective. Agric For Meteorol 1998;90:1±25. [2] Baldocchi D, Falge E, Wilson K. A spectral analysis of biosphere± atmosphere trace gas ux densities and meteorological variables across hour to year time scales. Agric For Meteorol 2001;107: 1±27. [3] Campbell G, Norman J. An introduction to environmental biophysics. New York: Springer; [4] Collatz G, Ball J, Grivet C, Berry JA. Physiological and environmental regulation of stomatal conductance, photosynthesis and transpiration: a model that includes a laminar boundary layer. Agric For Meteorol 1991;54:107±36. [5] Ellsworth DS. Seasonal CO 2 assimilation and stomatal limitations in a Pinus taeda canopy with varying climate. Tree Physiol 2000;20:435±45. [6] Falge E, Baldocchi D, Olson R, Anthoni P, Aubinet M, Bernhofer Ch, et al. Gap lling strategies for defensible annual sums of net ecosystem exchange. Agric For Meteorol 2001;107:43±69. [7] Farquhar GD, Von Caemmerer S, Berry JA. A biochemical model of photosynthetic CO 2 assimilation in leaves of C3 species. Planta 1980;149:78±90. [8] Flandrin P. Time-scale analyses and self-similar stochastic processes. In: Byrnes, et al., editors. Wavelets and their applications. NATO ASI Ser, Vol p. 121±42. [9] Frisch U. Turbulence. Cambridge: Cambridge University Press; [10] Gu L, Shugart H, Fuente J, Black TA, Shewchuk S. Micrometeorological, biophysical exchanges and NEE decomposition in a two-story boreal forest ± development and test of an integrated model. Agric For Meteorol 1999;94:123±48. [11] Hutjes RWA, Kabat P, Running SW, Shuttleworth WJ, Field C, Bass B, et al. Biospheric aspects of the hydrologic cycle. J Hydrol 1998;213:1±21. [12] Kaimal JC, Finnigan JJ. In: Atmospheric boundary layer ows: their structure and measurements. Oxford: Oxford Press; p [13] Katul GG, Parlange MB, Chu CR. Intermittency, local isotropy, and non-gaussian statistics in atmospheric surface layer turbulence. Phys Fluids 1994;6:2480±92. [14] Katul GG, Parlange MB. Analysis of land surface heat uxes using the orthonormal wavelet approach. Water Resour Res 1995;31:2743±9. [15] Katul GG, Hsieh CI, Sigmon J. Energy-inertial scale interaction for temperature and velocity in the unstable surface layer. Bound Layer Meteorol 1997;82:49±80. [16] Katul GG, Hsieh CI, Kuhn G, Ellsworth D, Nie D. The turbulent eddy motion at the forest±atmosphere interface. J Geophys Res 1997;102:13409±21. [17] Katul GG, Oren R, Ellsworth D, Hsieh CI, Phillips N, Lewin K. A Lagrangian dispersion model for predicting CO 2 sources, sinks, and uxes in a uniform loblolly pine Pinus taeda L.) stand. J Geophys Res 1997;102:9309±21. [18] Katul GG, Vidakovic B. Identi cation of low-dimensional energy containing/ ux transporting eddy motion in the atmospheric surface layer using wavelet thresholding methods. J Atmos Sci 1998;54:91±103. [19] Katul GG, Geron CD, Hsieh CI, Vidakovic B, Guenther AB. Active turbulence and scalar transport near the land±atmosphere interface. J Appl Meteorol 1998;37:1533±46. [20] Katul GG, Albertson JD. An investigation of higher order closure models for a forested canopy. Bound Layer Meteorol 1998;89:47± 74. [21] Katul GG, Albertson JD. Low dimensional turbulent transport mechanics near the forest atmosphere interface. In: Muller, Vidakovic, editors. Bayesian inference in wavelet based models. Lect Notes Stat. Berlin: Springer; [22] Katul GG, Albertson JD. Modeling CO 2 sources, sinks, and uxes within a forest canopy. J Geophys Res 1999;104:6081±91. [23] Katul GG, Hsieh CI, Bowling D, et al. Spatial variability of turbulent uxes in the roughness sublayer of an even-aged pine forest. Bound Layer Meteorol 1999;93:1±28. [24] Katul GG, Chang W-H. Principal length scales in second-order closure models for canopy turbulence. J Appl Meteorol 1999;38:1631±48. [25] Kolmogorov AN. The local structure of turbulence in incompressible viscous uid for very large Reynolds numbers. Dokl Akad Nauk SSSR 1941;30:301±5. [26] Kumar P, Foufoula-Georgiou E. Wavelet analysis in geophysics: An introduction. In: Foufoula-Georgiou E, Kumar P, editors. Wavelets in geophysics. New York: Academic Press; p. 1±43. [27] Kumar P, Foufoula-Georgiou E. A multicomponent decomposition of spatial rainfall elds: 2. self-similarity in uctuations. Water Resour Res 1993;29:2533±44. [28] Lai CT, Katul GG, Ellsworth DS, Oren R. Modeling vegetation± atmosphere CO 2 exchange by a coupled Eulerian±Lagrangian approach. Bound Layer Meteorol 2000;95:91±122. [29] Lai CT, Katul GG, Oren R, Ellsworth D, Schafer K. Modeling CO 2 and water vapor turbulent ux distribution within a forest canopy. J Geophys Res 2000;105:26333±51. [30] Mallat S. A wavelet tour of signal processing. New York: Academic Press; p [31] Mallat S. A theory for multiresolution signal decomposition: The wavelet representation. IEEE Trans Patt Anal Mach Intell 1989;11:674±93. [32] Medlyn BE, Badeck FW, DePury DGG, Barton CVM, Broadmeadow M, Ceulemans R, et al. E ects of elevated CO 2 on photosynthesis in European forest species: a meta-analysis of model parameters. Plant Cell Environ 1999;22:1475±95. [33] Mickler RA. Southern pine forests of eastern North America. In: Mickler RA, Fox SA, editors. The productivity and sustainability of southern forest ecosystems in a changing environment. Ecol Studies, Vol New York: Springer; p. 18±40. [34] Norman JM, Welles J. Radiative transfer in an array of canopies. Agron J 1983;75:481±8. [35] Oren R, Phillips N, Katul GG, Ewers BE, Pataki DE. Scaling xylem sap ux and soil water balance and calculating variance: a method for partitioning water ux in forests. Ann Sci For 1998;55:191±216. [36] Oren R, Ewers B, Todd P, Phillips N, Katul GG. Water balance delineates the soil layer in which moisture a ects canopy conductance. Ecol Appl 1998;8:990±1002. [37] Raupach MR. Canopy transport processes. In: Ste en WL, Denmead OT, editors. Flow and transport in the natural environment: advances and applications. Berlin: Springer; p. 95±127. [38] Raupach MR, Finnigan JJ, Brunet Y. Coherent eddies and turbulence in vegetation canopies: the mixing-layer analogy. Bound Layer Meteorol 1996;78:351±82. [39] Schuepp PH. Tansley review no. 59: leaf boundary layers. New Phytol 1993;125:477±507. [40] Sperry JS, Adler FR, Campbell GS, Comstock JP. Limitation of plant water use by rhizosphere and xylem conductance: results from a model. Plant, Cell, Environ 1998;21:347±59. [41] Thomson DJ. Criteria for the selection of stochastic models of particle trajectories in turbulent ows. J Fluid Mech 1987;180:529±56.
Supplement of Upside-down fluxes Down Under: CO 2 net sink in winter and net source in summer in a temperate evergreen broadleaf forest
Supplement of Biogeosciences, 15, 3703 3716, 2018 https://doi.org/10.5194/bg-15-3703-2018-supplement Author(s) 2018. This work is distributed under the Creative Commons Attribution 4.0 License. Supplement
More informationA NOTE ON THE CONTRIBUTION OF DISPERSIVE FLUXES TO MOMENTUM TRANSFER WITHIN CANOPIES. Research Note
A NOTE ON THE CONTRIBUTION OF DISPERSIVE FLUXES TO MOMENTUM TRANSFER WITHIN CANOPIES Research Note D. POGGI Dipartimento di Idraulica, Trasporti ed Infrastrutture Civili, Politecnico di Torino, Torino,
More informationAn objective method for determining principal time scales of coherent eddy structures using orthonormal wavelets
PII: S 0 3 0 9-1708(98)00046-3 Advances in Water Resources Vol. 22, No. 6, pp 561±566, 1999 Ó 1999 Elsevier Science Ltd Printed in Great Britain. All rights reserved 0309-1708/99/$ ± see front matter An
More informationThe structure of canopy turbulence and its implication to scalar dispersion
The structure of canopy turbulence and its implication to scalar dispersion Gabriel Katul 1,2,3 & Davide Poggi 3 1 Nicholas School of the Environment and Earth Sciences, Duke University, USA 2 Department
More informationSub-canopy. measurements in. Turbulenssista ja turbulenttisista pystyvoista mäntymetsän n latvuston alapuolella
Sub-canopy measurements in Hyytiälä,, SMEAR II station Samuli Launiainen s Master thesis Turbulenssista ja turbulenttisista pystyvoista mäntymetsän n latvuston alapuolella TKE-yht yhtälö latvuston sisäll
More informationRoughness Sub Layers John Finnigan, Roger Shaw, Ned Patton, Ian Harman
Roughness Sub Layers John Finnigan, Roger Shaw, Ned Patton, Ian Harman 1. Characteristics of the Roughness Sub layer With well understood caveats, the time averaged statistics of flow in the atmospheric
More informationP2.1 DIRECT OBSERVATION OF THE EVAPORATION OF INTERCEPTED WATER OVER AN OLD-GROWTH FOREST IN THE EASTERN AMAZON REGION
P2.1 DIRECT OBSERVATION OF THE EVAPORATION OF INTERCEPTED WATER OVER AN OLD-GROWTH FOREST IN THE EASTERN AMAZON REGION Matthew J. Czikowsky (1)*, David R. Fitzjarrald (1), Osvaldo L. L. Moraes (2), Ricardo
More informationDependence of evaporation on meteorological variables at di erent time-scales and intercomparison of estimation methods
Hydrological Processes Hydrol. Process. 12, 429±442 (1998) Dependence of evaporation on meteorological variables at di erent time-scales and intercomparison of estimation methods C.-Y. Xu 1 and V.P. Singh
More informationFlux Tower Data Quality Analysis in the North American Monsoon Region
Flux Tower Data Quality Analysis in the North American Monsoon Region 1. Motivation The area of focus in this study is mainly Arizona, due to data richness and availability. Monsoon rains in Arizona usually
More informationContents. 1. Evaporation
Contents 1 Evaporation 1 1a Evaporation from Wet Surfaces................... 1 1b Evaporation from Wet Surfaces in the absence of Advection... 4 1c Bowen Ratio Method........................ 4 1d Potential
More informationSpatial scales resolved. Four hierarchical scales: Leaf. Plant. Canopy. Landscape
Spatial scales resolved Four hierarchical scales: Leaf Plant Canopy Landscape Canopy embedded in a turbulent air stream - extracting momentum and releasing/uptaking scalars at the solid-fluid interface.
More information6 VORTICITY DYNAMICS 41
6 VORTICITY DYNAMICS 41 6 VORTICITY DYNAMICS As mentioned in the introduction, turbulence is rotational and characterized by large uctuations in vorticity. In this section we would like to identify some
More information3. Carbon Dioxide (CO 2 )
3. Carbon Dioxide (CO 2 ) Basic information on CO 2 with regard to environmental issues Carbon dioxide (CO 2 ) is a significant greenhouse gas that has strong absorption bands in the infrared region and
More informationActive Turbulence and Scalar Transport near the Forest Atmosphere Interface
VOLUME 37 JOURNAL OF APPLIED METEOROLOGY DECEMBER 1998 Active Turbulence and Scalar Transport near the Forest Atmosphere Interface GABRIEL G. KATUL School of the Environment and Center for Hydrologic Science,
More informationMeasuring Carbon Using Micrometeorological Methods
Measuring Carbon Using Micrometeorological Methods Theory Theory Sink = Vertical Transport + Change in Horizontal Transport + Change in Storage + Diffusion Storage within walls Diffusion across walls Transport
More informationPeter Molnar 1 and Paolo Burlando Institute of Environmental Engineering, ETH Zurich, Switzerland
Hydrology Days 6 Seasonal and regional variability in scaling properties and correlation structure of high resolution precipitation data in a highly heterogeneous mountain environment (Switzerland) Peter
More informationGapfilling of EC fluxes
Gapfilling of EC fluxes Pasi Kolari Department of Forest Sciences / Department of Physics University of Helsinki EddyUH training course Helsinki 23.1.2013 Contents Basic concepts of gapfilling Example
More informationEvapotranspiration. Here, liquid water on surfaces or in the very thin surface layer of the soil that evaporates directly to the atmosphere
Evapotranspiration Evaporation (E): In general, the change of state from liquid to gas Here, liquid water on surfaces or in the very thin surface layer of the soil that evaporates directly to the atmosphere
More informationA Note on the Estimation of Eddy Diffusivity and Dissipation Length in Low Winds over a Tropical Urban Terrain
Pure appl. geophys. 160 (2003) 395 404 0033 4553/03/020395 10 Ó Birkhäuser Verlag, Basel, 2003 Pure and Applied Geophysics A Note on the Estimation of Eddy Diffusivity and Dissipation Length in Low Winds
More informationEnvironmental Fluid Dynamics
Environmental Fluid Dynamics ME EN 7710 Spring 2015 Instructor: E.R. Pardyjak University of Utah Department of Mechanical Engineering Definitions Environmental Fluid Mechanics principles that govern transport,
More informationThis is the first of several lectures on flux measurements. We will start with the simplest and earliest method, flux gradient or K theory techniques
This is the first of several lectures on flux measurements. We will start with the simplest and earliest method, flux gradient or K theory techniques 1 Fluxes, or technically flux densities, are the number
More informationSPATIAL VARIABILITY OF TURBULENT FLUXES IN THE ROUGHNESS SUBLAYER OF AN EVEN-AGED PINE FOREST
SPATIAL VARIABILITY OF TURBULENT FLUXES IN THE ROUGHNESS SUBLAYER OF AN EVEN-AGED PINE FOREST GABRIEL KATUL 1,, CHENG-I HSIEH 1, DAVID BOWLING 2, KENNETH CLARK 3, NARASINHA SHURPALI 4, ANDREW TURNIPSEED
More informationAnalysis of meteorological measurements made over three rainy seasons in Sinazongwe District, Zambia.
Analysis of meteorological measurements made over three rainy seasons in Sinazongwe District, Zambia. 1 Hiromitsu Kanno, 2 Hiroyuki Shimono, 3 Takeshi Sakurai, and 4 Taro Yamauchi 1 National Agricultural
More informationAssimilation of satellite fapar data within the ORCHIDEE biosphere model and its impacts on land surface carbon and energy fluxes
Laboratoire des Sciences du Climat et de l'environnement Assimilation of satellite fapar data within the ORCHIDEE biosphere model and its impacts on land surface carbon and energy fluxes CAMELIA project
More informationEvapotranspiration. Andy Black. CCRN Processes Workshop, Hamilton, ON, Sept Importance of evapotranspiration (E)
Evapotranspiration Andy Black CCRN Processes Workshop, Hamilton, ON, 12-13 Sept 213 Importance of evapotranspiration (E) This process is important in CCRN goals because 1. Major component of both terrestrial
More informationDiscussion Reply to comment by Rannik on A simple method for estimating frequency response corrections for eddy covariance systems. W.J.
Agricultural and Forest Meteorology 07 (200) 247 25 Discussion Reply to comment by Rannik on A simple method for estimating frequency response corrections for eddy covariance systems W.J. Massman USDA/Forest
More informationComparing independent estimates of carbon dioxide exchange over 5 years at a deciduous forest in the southeastern United States
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 106, NO. D24, PAGES 34,167 34,178, DECEMBER 27, 2001 Comparing independent estimates of carbon dioxide exchange over 5 years at a deciduous forest in the southeastern
More informationSeasonal and interannual relations between precipitation, soil moisture and vegetation in the North American monsoon region
Seasonal and interannual relations between precipitation, soil moisture and vegetation in the North American monsoon region Luis A. Mendez-Barroso 1, Enrique R. Vivoni 1, Christopher J. Watts 2 and Julio
More informationEffects of sub-grid variability of precipitation and canopy water storage on climate model simulations of water cycle in Europe
Adv. Geosci., 17, 49 53, 2008 Author(s) 2008. This work is distributed under the Creative Commons Attribution 3.0 License. Advances in Geosciences Effects of sub-grid variability of precipitation and canopy
More informationConvective Fluxes: Sensible and Latent Heat Convective Fluxes Convective fluxes require Vertical gradient of temperature / water AND Turbulence ( mixing ) Vertical gradient, but no turbulence: only very
More informationTemperature and light as ecological factors for plants
PLB/EVE 117 Plant Ecology Fall 2005 1 Temperature and light as ecological factors for plants I. Temperature as an environmental factor A. The influence of temperature as an environmental factor is pervasive
More informationAtmospheric Moisture, Precipitation, and Weather Systems
Atmospheric Moisture, Precipitation, and Weather Systems 6 Chapter Overview The atmosphere is a complex system, sometimes described as chaotic in nature. In this chapter we examine one of the principal
More informationWind and turbulence experience strong gradients in vegetation. How do we deal with this? We have to predict wind and turbulence profiles through the
1 2 Wind and turbulence experience strong gradients in vegetation. How do we deal with this? We have to predict wind and turbulence profiles through the canopy. 3 Next we discuss turbulence in the canopy.
More informationApproaches in modelling tritium uptake by crops
Approaches in modelling tritium uptake by crops EMRAS II Approaches for Assessing Emergency Situations Working Group 7 Tritium Accidents Vienna 25-29 January 2010 D. Galeriu, A Melintescu History Different
More informationModeling nighttime ecosystem respiration from measured CO 2 concentration and air temperature profiles using inverse methods
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 111,, doi:1.129/25jd5976, 26 Modeling nighttime ecosystem respiration from measured CO 2 concentration and air temperature profiles using inverse methods Jehn-Yih
More informationSimplified expressions for adjusting higher-order turbulent statistics obtained from open path gas analyzers
Boundary-Layer Meteorol (27) 22:25 26 DOI.7/s546-6-95- ORIGINAL PAPER Simplified expressions for adjusting higher-order turbulent statistics obtained from open path gas analyzers M. Detto G. G. Katul Received:
More informationSkin temperature perturbations induced by surface layer turbulence above a grass surface
WATER RESOURCES RESEARCH, VOL. 34, NO. 5, PAGES 1265 1274, MAY 1998 Skin temperature perturbations induced by surface layer turbulence above a grass surface Gabriel G. Katul, 1,2 John Schieldge, 3 Cheng-I
More informationSoil moisture variability and scale-dependency of nonlinear parameterizations in coupled land±atmosphere models
Advances in Water Resources 24 2001) 1143±1157 www.elsevier.com/locate/advwatres Soil moisture variability and scale-dependency of nonlinear parameterizations in coupled land±atmosphere models Deborah
More informationDiagnosing the Climatology and Interannual Variability of North American Summer Climate with the Regional Atmospheric Modeling System (RAMS)
Diagnosing the Climatology and Interannual Variability of North American Summer Climate with the Regional Atmospheric Modeling System (RAMS) Christopher L. Castro and Roger A. Pielke, Sr. Department of
More information5B.1 DEVELOPING A REFERENCE CROP EVAPOTRANSPIRATION CLIMATOLOGY FOR THE SOUTHEASTERN UNITED STATES USING THE FAO PENMAN-MONTEITH ESTIMATION TECHNIQUE
DEVELOPING A REFERENCE CROP EVAPOTRANSPIRATION CLIMATOLOGY FOR THE SOUTHEASTERN UNITED STATES USING THE FAO PENMAN-MONTEITH ESTIMATION TECHNIQUE Heather A. Dinon*, Ryan P. Boyles, and Gail G. Wilkerson
More informationEcosystem-Climate Interactions
Ecosystem-Climate Interactions Dennis Baldocchi UC Berkeley 2/1/2013 Topics Climate and Vegetation Correspondence Holdredge Classification Plant Functional Types Plant-Climate Interactions Canopy Microclimate
More informationComparative Plant Ecophysiology
Comparative Plant Ecophysiology 2. Plant traits and climate factors that form bases for eco- physiological comparison 3. Life form comparisons of: Stomatal conductance Photosynthesis Xylem Anatomy Leaf
More informationFlux Tower Data Quality Analysis. Dea Doklestic
Flux Tower Data Quality Analysis Dea Doklestic Motivation North American Monsoon (NAM) Seasonal large scale reversal of atmospheric circulation Occurs during the summer months due to a large temperature
More informationRadiation transfer in vegetation canopies Part I plants architecture
Radiation Transfer in Environmental Science with emphasis on aquatic and vegetation canopy medias Radiation transfer in vegetation canopies Part I plants architecture Autumn 2008 Prof. Emmanuel Boss, Dr.
More informationImproving canopy processes in the Community Land Model using Fluxnet data: Assessing nitrogen limitation and canopy radiation
Improving canopy processes in the Community Land Model using Fluxnet data: Assessing nitrogen limitation and canopy radiation Gordon Bonan, Keith Oleson, and Rosie Fisher National Center for Atmospheric
More informationEstimating global and local scaling exponents in turbulent flows using discrete wavelet transformations
PHYSICS OF FLUIDS VOLUME 13, NUMBER 1 JANUARY 2001 Estimating global and local scaling exponents in turbulent flows using discrete wavelet transformations Gabriel Katul a) School of the Environment, Duke
More informationThe role of soil moisture in influencing climate and terrestrial ecosystem processes
1of 18 The role of soil moisture in influencing climate and terrestrial ecosystem processes Vivek Arora Canadian Centre for Climate Modelling and Analysis Meteorological Service of Canada Outline 2of 18
More informationTemporal change of some statistical characteristics of wind speed over the Great Hungarian Plain
Theor. Appl. Climatol. 69, 69±79 (2001) 1 Department of Meteorology, University of Debrecen, Hungary 2 Department of Climatology and Landscape Ecology, University of Szeged, Hungary Temporal change of
More informationTurbulence is a ubiquitous phenomenon in environmental fluid mechanics that dramatically affects flow structure and mixing.
Turbulence is a ubiquitous phenomenon in environmental fluid mechanics that dramatically affects flow structure and mixing. Thus, it is very important to form both a conceptual understanding and a quantitative
More informationGEOG415 Mid-term Exam 110 minute February 27, 2003
GEOG415 Mid-term Exam 110 minute February 27, 2003 1 Name: ID: 1. The graph shows the relationship between air temperature and saturation vapor pressure. (a) Estimate the relative humidity of an air parcel
More informationCarbon Input to Ecosystems
Objectives Carbon Input Leaves Photosynthetic pathways Canopies (i.e., ecosystems) Controls over carbon input Leaves Canopies (i.e., ecosystems) Terminology Photosynthesis vs. net photosynthesis vs. gross
More informationCarbon Assimilation and Its Variation among Plant Communities
Carbon Assimilation and Its Variation among Plant Communities Introduction By, Susan Boersma, Andrew Wiersma Institution: Calvin College Faculty Advisor: David Dornbos Currently, global warming remains
More informationAppendix A. Stemflow
Appendix A Stemflow The amount of stemflow S f is directly related to P. Table A.1 shows the regression results for the four experimental sites at which stemflow was measured. To investigate a possible
More informationSoil Water Atmosphere Plant (SWAP) Model: I. INTRODUCTION AND THEORETICAL BACKGROUND
Soil Water Atmosphere Plant (SWAP) Model: I. INTRODUCTION AND THEORETICAL BACKGROUND Reinder A.Feddes Jos van Dam Joop Kroes Angel Utset, Main processes Rain fall / irrigation Transpiration Soil evaporation
More informationCHAPTER 2 DATA AND METHODS. Errors using inadequate data are much less than those using no data at all. Charles Babbage, circa 1850
CHAPTER 2 DATA AND METHODS Errors using inadequate data are much less than those using no data at all. Charles Babbage, circa 185 2.1 Datasets 2.1.1 OLR The primary data used in this study are the outgoing
More informationAgricultural and Forest Meteorology 2915 (2000) 1 27
Abstract Agricultural and Forest Meteorology 2915 (2000) 1 27 A spectral analysis of biosphere atmosphere trace gas flux densities and meteorological variables across hour to multi-year time scales Dennis
More informationThe Spatial Variability of Energy and Carbon Dioxide Fluxes at the Floor of a Deciduous Forest
University of Nebraska - Lincoln DigitalCommons@University of Nebraska - Lincoln U.S. Bureau of Land Management Papers U.S. Department of the Interior 2001 The Spatial Variability of Energy and Carbon
More informationImproved Fields of Satellite-Derived Ocean Surface Turbulent Fluxes of Energy and Moisture
Improved Fields of Satellite-Derived Ocean Surface Turbulent Fluxes of Energy and Moisture First year report on NASA grant NNX09AJ49G PI: Mark A. Bourassa Co-Is: Carol Anne Clayson, Shawn Smith, and Gary
More informationBoundary layer equilibrium [2005] over tropical oceans
Boundary layer equilibrium [2005] over tropical oceans Alan K. Betts [akbetts@aol.com] Based on: Betts, A.K., 1997: Trade Cumulus: Observations and Modeling. Chapter 4 (pp 99-126) in The Physics and Parameterization
More informationLength Scale Analysis of Surface Energy Fluxes Derived from Remote Sensing
1212 JOURNAL OF HYDROMETEOROLOGY VOLUME 4 Length Scale Analysis of Surface Energy Fluxes Derived from Remote Sensing NATHANIEL A. BRUNSELL* Department of Plants, Soils, and Biometeorology, Utah State University,
More informationModeling CO 2 sources, sinks, and fluxes within a forest canopy
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 104, NO. D6, PAGES 6081 6091, MARCH 27, 1999 Modeling CO 2 sources, sinks, and fluxes within a forest canopy Gabriel G. Katul School of the Environment, Duke University,
More informationAn experimental investigation of the thermal performance of an asymmetrical at plate heat pipe
International Journal of Heat and Mass Transfer 43 (2000) 2657±2668 www.elsevier.com/locate/ijhmt An experimental investigation of the thermal performance of an asymmetrical at plate heat pipe Y. Wang,
More informationLecture 7: The Monash Simple Climate
Climate of the Ocean Lecture 7: The Monash Simple Climate Model Dr. Claudia Frauen Leibniz Institute for Baltic Sea Research Warnemünde (IOW) claudia.frauen@io-warnemuende.de Outline: Motivation The GREB
More informationTerrestrial land surfacesa pot pourri
CALTECH JPL Center for Climate Sciences March 26, 2018 Terrestrial land surfacesa pot pourri Graham Farquhar Australian National University What do we want from our models? Timescale is a key issue What
More informationApril 1990 T. Watanabe and J. Kondo 227. The Influence of Canopy Structure and Density. upon the Mixing Length within and above Vegetation
April 1990 T. Watanabe and J. Kondo 227 The Influence of Canopy Structure and Density upon the Mixing Length within and above Vegetation By Tsutomu Watanabe and Junsei Kondo Geophysical Institute, Tohoku
More informationUsing mathematical inverse theory to estimate respiratory and photosynthetic fluxes in a heterogeneous conifer canopy
Using mathematical inverse theory to estimate respiratory and photosynthetic fluxes in a heterogeneous conifer canopy John M. Zobitz with David R. Bowling, Frederick R. Adler, James P. Keener, Jerome Ogée
More informationThe in uence of a forest canopy on top-down and bottom-up diffusion in the planetary boundary layer
Q. J. R. Meteorol. Soc. (2003), 129, pp. 1415 1434 doi: 10.1256/qj.01.175 The in uence of a forest canopy on top-down and bottom-up diffusion in the planetary boundary layer By E. G. PATTON 1, P. P. SULLIVAN
More informationEVALUATING LAND SURFACE FLUX OF METHANE AND NITROUS OXIDE IN AN AGRICULTURAL LANDSCAPE WITH TALL TOWER MEASUREMENTS AND A TRAJECTORY MODEL
8.3 EVALUATING LAND SURFACE FLUX OF METHANE AND NITROUS OXIDE IN AN AGRICULTURAL LANDSCAPE WITH TALL TOWER MEASUREMENTS AND A TRAJECTORY MODEL Xin Zhang*, Xuhui Lee Yale University, New Haven, CT, USA
More informationEvapotranspiration: Theory and Applications
Evapotranspiration: Theory and Applications Lu Zhang ( 张橹 ) CSIRO Land and Water Evaporation: part of our everyday life Evapotranspiration Global Land: P = 800 mm Q = 315 mm E = 485 mm Evapotranspiration
More information282 Journal of the Meteorological Society of Japan Vol. 60, No. 1. A Theory and Method of Long-Range Numerical
282 Journal of the Meteorological Society of Japan Vol. 60, No. 1 A Theory and Method of Long-Range Numerical Weather Forecasts By Chao Jih-Ping, Guo Yu-Fu and Xin Ru-Nan Institute of Atmospheric Physics,
More informationThe Ocean-Atmosphere System II: Oceanic Heat Budget
The Ocean-Atmosphere System II: Oceanic Heat Budget C. Chen General Physical Oceanography MAR 555 School for Marine Sciences and Technology Umass-Dartmouth MAR 555 Lecture 2: The Oceanic Heat Budget Q
More informationLecture 3A: Interception
3-1 GEOG415 Lecture 3A: Interception What is interception? Canopy interception (C) Litter interception (L) Interception ( I = C + L ) Precipitation (P) Throughfall (T) Stemflow (S) Net precipitation (R)
More informationThe Heat Budget for Mt. Hope Bay
The School for Marine Science and Technology The Heat Budget for Mt. Hope Bay Y. Fan and W. Brown SMAST, UMassD SMAST Technical Report No. SMAST-03-0801 The School for Marine Science and Technology University
More informationRegional dry-season climate changes due to three decades of Amazonian deforestation
In the format provided by the authors and unedited. SUPPLEMENTARY INFORMATION DOI:./NCLIMATE Regional dry-season climate changes due to three decades of Amazonian deforestation Jaya problemkhanna by using
More informationLecture 3. Turbulent fluxes and TKE budgets (Garratt, Ch 2)
Lecture 3. Turbulent fluxes and TKE budgets (Garratt, Ch 2) The ABL, though turbulent, is not homogeneous, and a critical role of turbulence is transport and mixing of air properties, especially in the
More informationMeteorology. Chapter 15 Worksheet 1
Chapter 15 Worksheet 1 Meteorology Name: Circle the letter that corresponds to the correct answer 1) The Tropic of Cancer and the Arctic Circle are examples of locations determined by: a) measuring systems.
More information7. Basics of Turbulent Flow Figure 1.
1 7. Basics of Turbulent Flow Whether a flow is laminar or turbulent depends of the relative importance of fluid friction (viscosity) and flow inertia. The ratio of inertial to viscous forces is the Reynolds
More informationOCN 401. Photosynthesis
OCN 401 Photosynthesis Photosynthesis Process by which carbon is reduced from CO 2 to organic carbon Provides all energy for the biosphere (except for chemosynthesis at hydrothermal vents) Affects composition
More informationONE OF THE MOST WIDELY USED APPROACHES to modeling
Analysis of Scaling-Up Resistances from Leaf to Canopy Using Numerical Simulations Adriana C. Furon, Jon S. Warland,* and Claudia Wagner-Riddle ABSTRACT A multi-layer model, combining Lagrangian dispersion
More informationBenefits of NT over CT. Water conservation in the NT benefits from reduced ET and runoff, and increased infiltration.
Benefits of NT over CT Water conservation in the NT benefits from reduced ET and runoff, and increased infiltration. Weed control. Increased water and root penetration Uniform stands. Typically 4 to 8
More informationMay 3, :41 AOGS - AS 9in x 6in b951-v16-ch13 LAND SURFACE ENERGY BUDGET OVER THE TIBETAN PLATEAU BASED ON SATELLITE REMOTE SENSING DATA
Advances in Geosciences Vol. 16: Atmospheric Science (2008) Eds. Jai Ho Oh et al. c World Scientific Publishing Company LAND SURFACE ENERGY BUDGET OVER THE TIBETAN PLATEAU BASED ON SATELLITE REMOTE SENSING
More informationOctober 1991 J. Wang and Y. Mitsuta 587 NOTES AND CORRESPONDENCE. Turbulence Structure and Transfer Characteristics
October 1991 J. Wang and Y. Mitsuta 587 NOTES AND CORRESPONDENCE Turbulence Structure and Transfer Characteristics in the Surface Layer of the HEIFE Gobi Area By Jiemin Wang Lanzhou Institute of Plateau
More information1 BASIC CONCEPTS AND MODELS
1 BASIC CONCEPTS AND ODELS 1.1 INTRODUCTION This Volume III in the series of textbooks is focused on applications of environmental isotopes in surface water hydrology. The term environmental means that
More informationAQRP Monthly Technical Report
AQRP Monthly Technical Report PROJECT TITLE PROJECT PARTICIPANTS Improving Modeled Biogenic Isoprene Emissions under Drought Conditions and Evaluating Their Impact on Ozone Formation Qi Ying, Gunnar W.
More informationA new lidar for water vapor and temperature measurements in the Atmospheric Boundary Layer
A new lidar for water vapor and temperature measurements in the Atmospheric Boundary Layer M. Froidevaux 1, I. Serikov 2, S. Burgos 3, P. Ristori 1, V. Simeonov 1, H. Van den Bergh 1, and M.B. Parlange
More informationSimulating Carbon and Water Balances in the Southern Boreal Forest. Omer Yetemen, Alan Barr, Andrew Ireson, Andy Black, Joe Melton
Simulating Carbon and Water Balances in the Southern Boreal Forest Omer Yetemen, Alan Barr, Andrew Ireson, Andy Black, Joe Melton Research Questions: How will climate change (changes in temperature and
More informationNote the diverse scales of eddy motion and self-similar appearance at different lengthscales of the turbulence in this water jet. Only eddies of size
L Note the diverse scales of eddy motion and self-similar appearance at different lengthscales of the turbulence in this water jet. Only eddies of size 0.01L or smaller are subject to substantial viscous
More informationTURBULENT STATISTICS OF NEUTRALLY STRATIFIED FLOW WITHIN AND ABOVE A SPARSE FOREST FROM LARGE-EDDY SIMULATION AND FIELD OBSERVATIONS
TURBULENT STATISTICS OF NEUTRALLY STRATIFIED FLOW WITHIN AND ABOVE A SPARSE FOREST FROM LARGE-EDDY SIMULATION AND FIELD OBSERVATIONS HONG-BING SU 1,, ROGER H. SHAW 1,KYAWTHAPAWU 1, CHIN-HOH MOENG 2 and
More informationChristopher L. Castro Department of Atmospheric Sciences University of Arizona
Spatiotemporal Variability and Covariability of Temperature, Precipitation, Soil Moisture, and Vegetation in North America for Regional Climate Model Applications Christopher L. Castro Department of Atmospheric
More informationA FIRST INVESTIGATION OF TEMPORAL ALBEDO DEVELOPMENT OVER A MAIZE PLOT
1 A FIRST INVESTIGATION OF TEMPORAL ALBEDO DEVELOPMENT OVER A MAIZE PLOT Robert Beyer May 1, 2007 INTRODUCTION Albedo, also known as shortwave reflectivity, is defined as the ratio of incoming radiation
More informationONE DIMENSIONAL CLIMATE MODEL
JORGE A. RAMÍREZ Associate Professor Water Resources, Hydrologic and Environmental Sciences Civil Wngineering Department Fort Collins, CO 80523-1372 Phone: (970 491-7621 FAX: (970 491-7727 e-mail: Jorge.Ramirez@ColoState.edu
More informationThe Planetary Boundary Layer and Uncertainty in Lower Boundary Conditions
The Planetary Boundary Layer and Uncertainty in Lower Boundary Conditions Joshua Hacker National Center for Atmospheric Research hacker@ucar.edu Topics The closure problem and physical parameterizations
More informationFor the operational forecaster one important precondition for the diagnosis and prediction of
Initiation of Deep Moist Convection at WV-Boundaries Vienna, Austria For the operational forecaster one important precondition for the diagnosis and prediction of convective activity is the availability
More information% FOREST LEAF AREA. Figure I. Structure of the forest in proximity of the Proctor Maple Research Center -~--~ ~
NTRODUCTON There is a critical need to develop methods to address issues of forest canopy productivity and the role of environmental conditions in regulating forest productivity. Recent observations of
More informationAssimilation of satellite derived soil moisture for weather forecasting
Assimilation of satellite derived soil moisture for weather forecasting www.cawcr.gov.au Imtiaz Dharssi and Peter Steinle February 2011 SMOS/SMAP workshop, Monash University Summary In preparation of the
More informationThe North Atlantic Oscillation: Climatic Significance and Environmental Impact
1 The North Atlantic Oscillation: Climatic Significance and Environmental Impact James W. Hurrell National Center for Atmospheric Research Climate and Global Dynamics Division, Climate Analysis Section
More informationObservation and Modeling of Net Ecosystem Carbon Exchange Over Canopy
Chapter 1 Observation and Modeling of Net Ecosystem Carbon Exchange Over Canopy Tomo omi Kumagai* Institute for Space-Earth Environmental, Nagoya University, Chikusa-ku, Nagoya 464-861, Japan Summary...269
More informationData Analysis and Modeling with Stable Isotope Ratios. Chun-Ta Lai San Diego State University June 2008
Data Analysis and Modeling with Stable Isotope Ratios Chun-Ta Lai San Diego State University June 2008 Leaf water is 18 O-enriched via transpiration δ 18 O vapor : -12 H 2 16 O H 2 18 O δ 18 O leaf : +8
More informationTwo requirements for evaporation from natural surfaces. The atmospheric boundary layer. Vertical structure of the ABL. Bloc 1. Hydrologie quantitative
3.1 Introduction Figure 3.1 M1 SDE MEC558 Hydrologie continentale et ressources en eau Continental Hydrology and Water Resources Two requirements for evaporation from natural surfaces Bloc 1. Hydrologie
More informationPrecipitation processes in the Middle East
Precipitation processes in the Middle East J. Evans a, R. Smith a and R.Oglesby b a Dept. Geology & Geophysics, Yale University, Connecticut, USA. b Global Hydrology and Climate Center, NASA, Alabama,
More information