The Kyoto Protocol Emission Allocations: Windfall Surpluses for Russia and Ukraine

The emission targets adopted in the Kyoto Protocol¹ far exceed thelikely level of emissions from Russia and Ukraine. These countries could selltheir surplus if the Protocol is followedand industrialized countries establish an international emission tradingsystem. Critics have condemned the potentialsale and dubbed the surplus hot air because it does not represent anyreduction in emissions below the level thatwould have occurred anyway. Using the most recent, comprehensive regionalscenarios² for the emissions of carbon dioxide from the energysystem, we estimate that during the Protocol's2008–2012 budget period the surplus will range from 9 MtC (milliontons of carbon) to 900 MtC for Russia andfrom 3 MtC to 200 MtC for Ukraine. Even scenarios with high economic growthand carbon-intensive technologies donot exhaust the surplus before the budget period. In the central (middlecourse) scenario, the total carbon surplusexceeds 1000 MtC and is worth 22 to 170 billion U.S. dollars (4 to 34 billionU.S. dollars per year). This flow ofrevenues, which could exceed Russian earnings from natural gas exports($10 billion in 1997³), is comparable with the projectedtotal investmentsof the Russian energy system for 2008–2012. If directed towardslow-carbon infrastructure investments (e.g., gaspipelines), surplus transfers could reinforce and partially lock-indecarbonization of the world energy system.     

A Comparative Analysis of Transpiration and Bare Soil Evaporation

Transpiration E^v and bare soil evaporation E^b processes are comparatively analysed assuming homogeneous and inhomogeneous areal distributions of volumetric soil moisture content . For a homogeneous areal distribution of we use a deterministic model, while for inhomogeneous distributions a statistical-deterministic diagnostic surface energy balance model is applied. The areal variations of are simulated by Monte-Carlo runs assuming normal distributions of .The numerical experiments are performed for loam. In the experiments we used different parameterizations for vegetation and bare soil surface resistances and strong atmospheric forcing. According to the results theE^v()-^Eb() differences are great, especially in dry conditions. In spite of this, the available energy flux curves of vegetation A^v() and bare soil A^b() surfaces differ much less than the E^v() and ^Eb() curves. The results suggest that E^v is much more non-linearly related to environmental conditions than ^Eb. Both E^v and ^Eb depend on the distribution of , the wetness regime and the parameterization used. With the parameterizations, ^Eb showed greater variations than E^v. These results are valid when there are no advective effects or mesoscale circulation patterns and the stratification is unstable.     

Reactions of alkoxy radicals under atmospheric conditions: The relative importance of decomposition versus reaction with O2

The reactions of alkoxy radicals determine to a large extent the products formed during the atmospheric degradations of emitted organic compounds. Experimental data concerning the decompositions, 1,5-H shift isomerizations and reactions with O₂ of several classes of alkoxy radicals are inconsistent with literature estimations of their absolute or relative rate constants. An alternative, although empirical, method for assessing the relative importance under atmospheric conditions of the reactions of alkoxy radicals with O₂ versus decomposition was derived. This estimation method utilizes the differences in the heats of reaction, (H)=(H_{decomposition}–H_{O 2 reaction}), between these two reactions pathways. For (H)[22–0.5(H_{O 2 reaction})], alkoxy radical decomposition dominates over the reaction with O₂ at room temperature and atmospheric pressure of air, while for (H)[25-0.5(H_{O 2 reaction})], the O₂ reaction dominates over decomposition (where the units of H are in kcal mol^–1). The utility and shortcomings of this approach are discussed. It is concluded that further studies concerning the reactions of alkoxy radicals are needed.     

One-dimensional checkerboards and blending heights

We analyze the checkerboard problem of many alternating surfaces with different properties, on scales up to (say) 3,000 m. Power-law representations of the vertical profiles of mean wind speed and eddy diffusivity lead to solutions in terms of Kelvin and trigonometric functions.These solutions are used to determine blending heights (_*), where deviations from the mean of concentration, or of vertical flux density, fall to some small fraction, , of their value at the surface. Values of _*are important for regional and larger-scale meteorological models. In smaller scale micrometeorological studies, they may serve also as the top levels of surface boundary layers.An important result for both theoretical and experimental contexts is that deviations of flux persist with elevation much more strongly than those of concentration, so that, in general, _* should be based on flux rather than concentration. Representative values of _*, for = 0.05, are of order 5 and 30 m for surface pattern wavelengths of 10² and 10³ m, respectively. Values of _* are robust to changes in adopted power-law indices, and are independent of wind speed. Surface roughness has a mild but calculable effect.     

Measurements of solar radiation and estimation of optical depth in the high Andes of Peru

Summary During an expedition to the high Andes of Southern Peru in June–July 1977, measurements of direct solar radiation in four spectral bands (0.270–0.530–0.630–0.695–2.900 ) were conducted at six sites in elevations ranging from sea level to 5645 m. These measurements were evaluated in Langley plots to determine total optical depths () and irradiances at the top of the atmosphere. In addition, water vapor optical depths (_wv) were calculated from the mean radiosounding over Lima during the expedition, and Rayleigh (_ray) and ozone (_oz) optical depths were obtained from published tabulations. Subtracting _ray, _oz, and _wv from yielded estimates of aerosol optical depth _aer. The components _ray and _oz decrease from the shorter towards the longer wavelength bands and from the lower towards the higher elevation sites; _aer also decreases towards the higher elevations. Particularly pronounced is the decrease of _aer and from the lowlands of the Pacific coast to the highlands of the interior, reflecting the effect of a persistent lower-tropospheric inversion and the contrast from the marine boundary layer to the clear atmosphere of the high Andes.With 4 Figures     

Near-field dispersion from instantaneous sources in the surface layer

This paper considers the near-field dispersion of an ensemble of tracer particles released instantaneously from an elevated source into an adiabatic surface layer. By modelling the Lagrangian vertical velocity as a Markov process which obeys the Langevin equation, we show analytically that the mean vertical drift velocity w(t) is w()=bu _*(1–e ^–(1+)), where is time since release (nondimensionalized with the Lagrangian time scale at the source), b Batchelor's constant, and u _*, the friction velocity. Hence, the mean height and mean depth of the ensemble are calculated. Although the derivation is formally valid only when 1, the predictions for w, mean height and mean depth are consistent in the downstream limit ( 1) with surface-layer Lagrangian similarity theory and with the diffusion equation. By comparing the analytical predictions with numerical, randomflight solutions of the Langevin equation, the analytical predictions are shown to be good approximations at all times, both near-field and far-field.     

Microwave vegetation optical depth and inverse modelling of soil emissivity using Nimbus/SMMR satellite observations

Summary A radiative transfer model has been used to determine the large scale effective 6.6 GHz and 37 GHz optical depths of the vegetation cover. Knowledge of the vegetation optical depth is important for satellite-based large scale soil moisture monitoring using microwave radiometry. The study is based on actual observed large scale surface soil moisture data and observed dual polarization 6.6 and 37 GHz Nimbus/SMMR brightness temperatures over a 3-year period. The derived optical depths have been compared with microwave polarization differences and polarization ratios in both frequencies and with Normalized Difference Vegetation Index (NDVI) values from NOAA/AVHRR. A synergistic approach to derive surface soil emissivity from satellite observed brightness temperatures by inverse modelling is described. This approach improves the relationship between satellite derived surface emissivity and large scale top soil moisture fromR ²=0.45 (no correction for vegetation) toR ²=0.72 (after correction for vegetation). This study also confirms the relationship between the microwave-based MPDI and NDVI earlier described and explained in the literature.List of Symbols f frequency [Hz] - f _i(p) fractional absorption at polarizationp - h surface roughness - h h cos² - H horizontal polarization - n _i complex index of refraction - p polarization (H orV) - R _s microwave surface reflectivity - T _B(p) brightness temperature at polarizationp - T ^* normalized brightness temperature - T polarization difference (T _v-T _H) - T _s temperature of soil surface - T _c temperature of canopy - T _max daily maximum air temperature - T _min daily minimum air temperature - V vertical polarization - soil moisture distribution factor; also used for the constant to partition the influence of bound and free water components to the dielectric constant of the mixture - empirical complex constant related to soil texture - microwave transmissivity of vegetation (=e ^–) - ^* effective transmissivity of vegetation (assuming =0) - microwave emissivity - _s emissivity of smooth soil surface - _rs emissivity of rough soil surface - _vs emissivity of vegetated surface - soil moisture content (% vol.) - K dielectric constant [F·m^–1] - K _fw dielectric constant of free water [F·m^–1] - K _ss dielectric constant of soil solids [F·m^–1] - K _m dielectric constant of mixture [F·m^–1] - K _o permittivity of free space [8.854·10^–12 F·m^–1] - high frequency limit ofK _wf [F·m^–1] - wavelength [m] - incidence angle [degrees from nadir] - polarization ratio (T _H/T _V) - b soil bulk density [gr·cm^–3] - s soil particle density [gr·cm^–3] - R surface reflectivity in red portion of spectrum - NIR surface reflectivity in near infrared portion of spectrum - _eff effective conductivity of soil extract [mS·cm^–1] - vegetation optical depth - _6.6 vegetation optical depth at 6.6 GHz - ₃₇ vegetation optical depth at 37 GHz - ^* effective vegetation optical depth (assuming =0) - single scattering albedo of vegetation With 12 Figures     

Free convection scaling over the ocean derived from buoy measurements during gate

The structure of atmospheric turbulence in the surface layer over the open ocean is examined under conditions of local free convection. The raw data consist of profile and fluctuation measurements of wind and temperature as obtained from a meteorological buoy. For near neutral conditions and for waves running approximately along the wind direction, wave-induced wind fluctuations can be described by a simplified linear theory based on Miles (1957). In this case, the spectrum of wind velocity is given as the sum of two parts; for the turbulent part, the parameterization as obtained by Kaimal et al. (1972) applies, while the wave-induced part is parameterized using a simplification of Miles' linear theory. For cases of local free convection, the measurements of the vertical component of the wind velocity are well described by similarity theory; as expected, _w /(-uw)^1/2 is proportional to (- z/L)^1/3. In order to scale the longitudinal wind velocity component, it seems to be reasonable to extend the list of relevant parameters by the height of the mixed layer z _i. We obtain _u /(- uw)^1/2 (z/z _i)^1/3(- z/L)^1/3 with only a poor correlation coefficient of r = 0.6. Overall, the results of local free convection scaling obtained from direct measurements show good agreement with those obtained from profile measurements. A comparison between direct and indirect determination of turbulent fluxes of momentum shows an unexplained difference of about 20%. This discrepancy is mainly due to a gap in the uw-cospectrum at the swell frequency.     

Mass transfer to bean leaves

A vapour of radio-lead (²¹²Pb) has been used to measure the Sherwood number, Sh, of model leaves at various angles of incidence,, to the airstream in a wind tunnel. The results for=0 are compared with Pohlhausen's formula and the results for 0, with Powell's experiments. The local values of Sh on the upwind and downwind sides of discs have been obtained. For leaves in the canopy, Sh was found to be about 25% greater than would be predicted by applying Pohlhausen's equation without correction for orientation.     

In-water and remote measurements of ocean color

Spectral measurements of downwelling irradiance, E _d(), above the surface, and of upwelling irradiance just below the surface, E _u(), allow computation of spectral values of the diffuse reflectance R() = E _u()/E _d(); this yields full information about the true color and brightness of the ocean. Typical results are presented and interpreted for waters very different in turbidity and phytoplankton content. Conversely, the possibility of infering the water content from R() data at selected wavelengths is examined in terms of the respective number of equations and unknowns. The necessary use of assumptions and of empirical laws is emphasized.The magnitude of the useful signal emerging from the water, and the magnitude of the additional signals due to specular reflexion at the interface and to atmospheric scattering, are compared on the basis of spectroradiometric measurements performed within and above the sea, from different altitudes. These unwanted signals are dominant, causing a drastic change in the spectral composition of the light received by a remote sensor. The evaluation of the atmospheric effect must be very precise in order to recover the marine signal with a sufficient accuracy for a meaningful application of any kind of algorithm.This work was supported in part by the Centre National d'Exploitation des Océans (under contract CNEXO 77/1695) and in part by the Centre National de la Recherche Scientifique (RCP 247 & ERA 278).     

Atmospheric effects in the remote sensing of phytoplankton pigments

We investigate the accuracy with which relevant atmospheric parameters must be estimated to derive phytoplankton pigment concentrations (chlorophyll a plus phaeophytin a ) of a given accuracy from measurements of the ocean's apparent spectral radiance at satellite altitudes. The analysis is limited to an instrument having the characteristics of the Coastal Zone Color Scanner scheduled to orbit the Earth on NIMBUS-G. A phytoplankton pigment algorithm is developed which relates the pigment concentration (C) to the three ratios of upwelling radiance just beneath the sea surface which can be formed from the wavelengths () 440, 520 and 550 nm. The pigment algorithm explains from 94 to 98% of the variance in log₁₀ C over three orders of magnitude in pigment concentration. This is combined with solutions to the radiative transfer equation to simulate the ocean's apparent spectral radiance at satellite altitudes as a function of C and the optical properties of the aerosol, the optical depth of which is assumed to be proportioned to ^-n . A specific atmospheric correction algorithm, based on the assumption that the ocean is totally absorbing at 670 nm, is then applied to the simulated spectral radiance, from which the pigment concentration is derived. Comparison between the true and derived values of C show that: (1) n is considerably more important than the actual aerosol optical thickness; (2) for C 0299-1 0.2 g l^-1 acceptable concentrations can be determined as long as n is not overestimated; (3) as C increases, the accuracy with which n must be estimated, for a given relative accuracy in C, also increases; and (4) for C greater than about 0.5 g 1^-1, the radiance at 440 nm becomes essentially useless in determining C. The computations also suggest that if separate pigment algorithms are used for C 1gl^-1 and C 1 gl^-1, accuracies considerably better than ±± in log C can be obtained for C 1 g l^-1 with only a coarse estimate of n, while for C 10 gl^-1, this accuracy can be achieved only with very good estimates of n.Contribution No. 387 from the NOAA/ERL Pacific Environmental Laboratory.On leave from Department of Physics, University of Miami, Coral Gables, Florida.     

Spatial Variability of Turbulent Fluxes in the Roughness Sublayer of an Even-Aged Pine Forest

The spatial variability of turbulent flow statistics in the roughness sublayer (RSL) of a uniform even-aged 14 m (= h) tall loblolly pine forest was investigated experimentally. Using seven existing walkup towers at this stand, high frequency velocity, temperature, water vapour and carbon dioxide concentrations were measured at 15.5 m above the ground surface from October 6 to 10 in 1997. These seven towers were separated by at least 100m from each other. The objective of this study was to examine whether single tower turbulence statistics measurements represent the flow properties of RSL turbulence above a uniform even-aged managed loblolly pine forest as a best-case scenario for natural forested ecosystems. From the intensive space-time series measurements, it was demonstrated that standard deviations of longitudinal and vertical velocities (_u, _w) and temperature (_T) are more planar homogeneous than their vertical flux of momentum (u_* ²) and sensible heat (H) counterparts. Also, the measured H is more horizontally homogeneous when compared to fluxes of other scalar entities such as CO₂ and water vapour. While the spatial variability in fluxes was significant (>15 %), this unique data set confirmed that single tower measurements represent the canonical structure of single-point RSL turbulence statistics, especially flux-variance relationships. Implications to extending the moving-equilibrium hypothesis for RSL flows are discussed. The spatial variability in all RSL flow variables was not constant in time and varied strongly with spatially averaged friction velocity u_*, especially when u_* was small. It is shown that flow properties derived from two-point temporal statistics such as correlation functions are more sensitive to local variability in leaf area density when compared to single point flow statistics. Specifically, that the local relationship between the reciprocal of the vertical velocity integral time scale (I_w) and the arrival frequency of organized structures (/h) predicted from a mixing-layer theory exhibited dependence on the local leaf area index. The broader implications of these findings to the measurement and modelling of RSL flows are also discussed.     

An observational study of heat fluxes and their relationships with net radiation

Experimental evidence indicates that the diurnal behaviour of the fluxes of heat into the ground and into the atmosphere versus the net flux of radiation can be modelled by closed curves, the hourly values folowing one another in either a clockwise or counter clockfashion. A general formulation to express the different heat fluxes as a function of net radiation is proposed. This relationship between the different heat fluxes and can be expressed as a sum of three terms: the first indicates a direct proportionality to , the second gives the deviation from linearity and depends on /t, and the third gives the value of the flux when = 0. The formulae are then expressed as a function of time and the ratios between the heat fluxes and are evaluated. A comparison with the approximations generally used shows that the latter may be considered as particular cases of the more general equations proposed here.     

Major wet interval in white mountains medieval warm period evidenced inδ 13C of bristlecone pine tree rings

A long ¹³C chronology was developed from bristlecone pine (Pinus longaeva) at the Methuselah Walk site in the White Mountains of California. The chronology represents cellulose from five-year ring groups pooled from multiple radii of multiple trees. The most dramatic isotopic event in the chronology appears from A.D. 1080–1129, when ¹³C values are depressed to levels ~ 2 below the mean for the period A.D. 925–1654. This isotopic excursion appears to represent a real event and is not an artifact of sampling circumstances; in fact, a similar excursion occurs in a previously-reported, independent ¹³C chronology from bristlecone pine. By carbon isotope fractionation models, the shift to low ¹³C values is consistent with abundant soil moisture, permitting leaf stomata to remain open, and allowing ready access of CO₂ from which carbon fixation may discriminate more effectively against¹³C in favor of¹²C. According to this model, the¹³C-depleted 50-yr isotopic excursion represents the wettest period in the White Mountains in the past 1000 yr, during which isotope-reconstructed July Palmer Drought Severity Indices averaged ~ +2.2.     

An improved calibration for tethered balloon wind measurements

A previously published technique for using tethered spherical balloons as anemometers for measuring light low-level winds has been further developed. Earlier data on the relationship between the aerodynamic drag coefficient and the Reynolds number of spherical rubber balloons were combined with a large number of new data and re-analysed; and the errors in the relationship were estimated. The results allowed a more accurate calculation of wind speed from the deflection of a tethered balloon from the vertical. When combined with a new technique for calculating the effects of the tether, this enabled light to moderate low-level winds at fixed heights up to 600 m or more to be measured with simple, cheap, and readily mobile equipment; and a slight modification of the technique allowed measurement of winds in and above fog. Wind speeds measured by the ballon technique showed reasonably good agreement with measurements by an anemometer carried beneath the balloon.Glossary of Symbols a, b, c Coefficients in the relationship between lnC _d and lnR - A Quantity under square root in solution for lnV whena0 - C _d Wind drag coefficient for balloon - C _dc Value ofC _d given by calibration curve of Table I - D Dynamic wind pressure force on balloon - F Buoyant free lift of balloon with load - Re Reynold's number of balloon (sphere) - R = Re/10⁵ - r Radius of sphere - T Tension in tether - V Wind speed - ₈₃() =(lnC _dc -lnC _d ) when 83° , or 0 for other - Error in lnC _d - Elevation of tether where attached to balloon - Elevation of balloon from ground tether point - Molecular viscosity of air - Ratio of circumference to diameter of circle - Density of air     

The prediction of dust erosion by wind: An interactive model

We have devised a partial differential equation for the prediction of dust concentration in a thin layer near the ground. In this equation, erosion (detachment), transport, deposition and source are parameterised in terms of known quantities. The interaction between a wind prediction model in the boundary layer and this equation affects the evolution of the dust concentration at the top of the surface layer. Numerical integrations are carried out for various values of source strength, ambient wind and particle size. Comparison with available data shows that the results appear very reasonable and that the model should be subjected to further development and testing.Notation (x, y, z, t) space co-ordinates and time (cm,t) - u, v components of horizontal wind speed (cm s^–1) - u _g, v_g components of the geostrophic wind (cm s^–1) - V=(u²+v²)^1/2 (cm s^–1) - (û v)= 1/(h – k) _k ^h(u, v)dz(cm s^–1) - V _* friction velocity (cm s^–1) - z ₀ roughness length (cm) - k ₁ von Karman constant =0.4 - V _d deposition velocity (cm s^–1) - V _g gravitational settling velocity (cm s^–1) - h height of inversion (cm) - k height of surface layer (cm) - potential temperature (°K) - _gr potential temperature at ground (°K) - _K potential temperature at top of surface layer (°K) - P pressure (mb) - P ₀ sfc pressure (mb) - C _p/C_v - (t)= /z lapse rate of potential temperature (°K cm^–1) - A(z) variation of wind with height in transition layer - B(z) variation of wind with height in transition layer - Cd drag coefficient - C _HO transfer coefficient for sensible heat - C dust concentration (g m^–3) - C _K dust concentration at top of surface layer (g m^–3) - D(z) variation with height of dust concentration - u, v, w turbulent fluctuations of the three velocity components (cm s^–1) - A ₁ constant coefficient of proportionality for heat flux =0.2 - Ri Richardson number - g gravitational acceleration =980 cm s^–2 - Re Reynolds number = - D _s thickness of laminar sub-layer (cm) - v molecular kinematic viscosity of air - coefficient of proportionality in source term - dummy variable - t time step (sec) - n time index in numerical equations On sabbatical leave at University of Aberdeen, Department of Engineering, September 1989–February 1990.     

Final results of the CONDORS convective diffusion experiment

TheConvectiveDiffusionObserved byRemoteSensors (CONDORS) field experiment conducted at the Boulder Atmospheric Observatory used innovative techniques to obtain three-dimensional mappings of plume concentration fields, /Q, of oil fog detected by lidar and chaff detected by Doppler radar. It included extensive meteorological measurements and, in 1983, tracer gases measured at a single sampling arc. Final results from ten hours of elevated and surface release data are summarized here. Many intercomparisons were made. Oil fog /Q measured 40m above the arc are mostly in good agreement withSF ₆ values, except in a few instances with large spacial inhomogeneities over short distances. After a correction scheme was applied to compensate for the effect of its settling speed, chaff dy/Q agreed well with those of oil except in two cases of oil fog hot spots. Mass or frequency distribution vs. azimuth or elevation angle comparisons were made for chaff, oil, and wind, with mostly good agreements. Spacial standard deviations, _y and _z, of chaff and oil agree overall and are consistent at short range with velocity standard deviations _vand _w 0.6w^* (the convective scale velocity), as measured atz>100m. Surface release _y is enhanced up to 60% at smallx, consistent with the Prairie Grass measurements and with larger _v and reduced wind speed measured near the surface. Decreased _y at small dimensionless average times is also noted. Finally, convectively scaled dy, C _y, were plotted versus dimensionlessx andz for oil, chaff, and corrected chaff for each 30–60 min period. Aggregated CONDORSC _y fields compare well with laboratory tank and LES numerical simulations; surface-released oil fog compares expecially well with the tank experiments. However, large deviations from the norm occurred in individual averaging periods; these deviations correlated strongly with anomalies in measured distributions.On assignment to the US Environmental Protection Agency, Atmospheric Research and Exposure Assessment Laboratory, RTP, NC.     

On the scavenging of gaseous nitrogen compounds by large and small rain drops: II. Wind tunnel and theoretical studies of the simultaneous uptake of NH3, SO2 and CO2 by water drops

An experimental investigation of the simultaneous absorption of NH₃ and SO₂ from the ambient atmosphere by freely falling water drops has been carried out in the Mainz vertical wind tunnel. The experimental results were found to be in good agreement with the results derived from computations with the Kronig-Brink convective diffusion model and also with a model which assumes a drop to be well mixed at all times. Encouraged by this agreement, these computation schemes for the uptake of gas by single drops where incorporated in a pollution washout model with realistic SO₂, NH₃ and CO₂ gas profiles. This model allows an entire raindrop size distribution to fall through a gas layer. The results of this plume-model show that the SO₂ uptake is strongly dependent on the NH₃ concentration in the atmosphere and on the rainrate. We also find that the small drops contribute more towards the washout of these gases. In the case of simultaneous presence of NH₃ and SO₂, desorption of these gases is negligible.     

Glaciers and Climate in the Andes between the Equator and 30° S: What is Recorded under Extreme Environmental Conditions?

Sublimation and melt disturb the environmental information obtained from ice core records in the Andes. In two case studies we demonstrate to what extent these post-depositional processes may remove major parts of the accumulated snow cover. Dark ash layers from the Tungurahua eruption changed the albedo of surface snow on Chimborazo glacier (6268 m, 1°30 S,78°36 W, Ecuador) between two ice core drilling campaigns and forced substantial melt. Re-distribution and washout of the chemical constituents shifted the concentration profiles obtained in December 1999 as compared to an equivalent core drilled in December 2000. The stable isotope records showed that approximately the water equivalent (weq) of an annual layer had melted, and that the percolating melt water penetrated within the firn layer to a depth of at least 16.5 m without refreezing. In the second example, from a site on the dry axis between the tropical and extra-tropical precipitation belts, significant loss of accumulated snow layers occurred by sublimation. A surface experiment at Cerro Tapado glacier (5536 m, 30°08 S,69°55 W, Chile) revealed that losses of 2 mm weq (5 mm snow) per day occurred during the dry period following the 1997/98 El Niño. This loss generally included the entire surface layer enriched in stable isotopes, and thus caused minimal disturbance of the isotopic signature (and hence climatic information) of the net accumulation, yet chemical constituents again experienced considerable changes in concentration. From annual layer counting and direct dating it is obvious that the major part of the accumulated ice on both glaciers is younger than 100 years; however, isotopic and chemical variations at least in the basal ice from Cerro Tapado clearly point to climate conditions different from the recent centuries. This evidence is supported by mass balance considerations derived from a glacier-climate model. The possibility of a third type of disturbance aside from sublimation and melting – in this case a significant hiatus in the environmental chronology – also deserves consideration for other icecore records from this region. Potential disruptions or discontinuities need to be carefully evaluated given the profound changes in climatic and glaciological conditions since the Last Glacial Maximum throughout Holocene times.     

A radiation and energy balance study of mature forest and clear-cut sites

Components of the radiation and energy balances were measured over a clear-cut area and a mature, mixed forest during the summer of 1981 at the Petawawa National Forestry Institute, Chalk River, Ontario. The work concentrated on the clear-cut site which supported a canopy layer composed primarily of bracken fern and logging remnants.Forty days of radiation data were collected at the clear-cut site. After the first four weeks of measurements (the green season), most of the ferns quickly died, and their foliage changed appearance from a green to brownish colour (the brown season). The daily mean reflection coefficient of solar radiation determined over the green season was 0.20 and decreased to 0.13 for the brown season. The corresponding value for the forest was 0.13, based on a limited amount of data. The clear-cut site received 11% and 21% less net radiation than the forest on a 24-hr and daylight-hours basis, respectively, as a consequence of the higher reflection coefficient and larger daytime longwave radiation emission.A reversing temperature difference measurement system (RTDMS), incorporating ten-junction thermopiles was employed at each site in order to determine Bowen ratios () via differential psychrometry. Both systems performed well, especially the RTDMS over the forest which was capable of resolving very small differences of temperature, typically less than 0.2 °C over a height of 3 m. The mean hourly Bowen ratio, calculated from values from 0800 to 1600 hr, varied from 0.2 to 1.0 for the forest and from 0.4 to 0.8 for the clear-cut site in the green season.A significant canopy heat storage component of the energy balance, Q _S , was found at the clear-cut site. In the early morning, a portion of the available energy was used to heat the biomass materials and air within the canopy layer. The stored heat within the canopy was released later in the day, increasing the available energy total.The daily mean value of the Priestley-Taylor coefficient (Priestley and Taylor, 1972) for the green season at the clear-cut site was 1.14, and individual values tended to increase during wet surface conditions and decrease when the surface dried. The daylight mean value during dry canopy conditions at the forest was 1.05, and much higher values occurred when the canopy was wet. The enhancement of for the wet forest was a result of the evaporation of intercepted rain (which is not limited by stomatal resistance) and the concomitant transfer of sensible heat to the forest.