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     

A Simple Mixing Length Formulation for the Eddy-Diffusivity Parameterization of Dry Convection

In this paper a simple mixing length formulation for the eddy-diffusivityparameterization of dry convection is suggested. The new formulation relates the mixinglength to the square root of the turbulent kinetic energy (e) and a time scale ( ):l = e. To close the parameterization the time scale is calculated as a functionof the boundary-layer height (h) and the convective velocity scale (w_*), h/w_*. Thesimpler approach of a constant time scale is also studied. The simulation of a case of dry atmosphericconvection with a one-dimensional boundary-layer model shows that the model with the new formulationreproduces quite well the main properties of the convective boundary layer. In particular,the entrainment is realistically represented by the new mixing length, which has the advantage of naturallydecreasing with the turbulent kinetic energy. Sensitivity studies to the surface flux and the lapserate, in the context of a simplified situation, show the robustness of the new formulation.     

Direct Measurement Of Shear Stress During Snow Saltation

Wind-tunnel experiments were carried out to measure the shear stressduring snow saltation. Shear stress acting on the snow surface, measured directly with a newly developed drag meter system, revealedthat the shear stress increased with the development ofsaltation. This result supports Owen's hypothesis that the saltationlayer acts as an increased roughness to the flow above the saltationlayer, leading to an increase in surface shear stress. To investigate the contribution of the grain borne shear stress_g and the fluid shear stress _f to the increment of the total shear stress _total, _g was calculated from the loss of horizontal momentum of saltating snowparticles. Since _g is the largest contribution to theincrement of _total, the collision of thesaltating particles is dominant for the shear stressmodification. The results qualitatively support the numericalsimulation reported by McEwan and Willetts.     

Cloud optical properties at Bergen (Norway) based on the analysis of long-term solar irradiance records

Summary Ground-based measurements of incoming solar irradiance and cloud observations during a 26 year period (1965–1990) at Bergen, Norway were used in conjunction with a comprehensive radiation model to infer the cloud optical depth under completely overcast conditions.Month-to-month and year-to-year (April through October) statistics of the cloud optical depth and observed cloud forms are presented. Some climate-related features, specifically, diurnal and seasonal variabilities in are examined. The effects of local cloudiness are pointed out and discussed. There appears to be a slight trend towards increasing cloud optical depth at noon during the warm period of the year. The possible uncertainties due to unknown size of cloud droplets are analyzed by model simulations. Possible directions for future research are suggested provided more meteorological and/or satellite information is available.With 9 Figures     

Second-order closure study of the covariance between chemically reactive species in the surface layer

A second-order modelling technique is used to investigate the influence of turbulence on chemical reactions. The covariance and variance equations for the NO-O₃-NO₂ system are developed as a function of the ratio of the timescale of turbulence ( _t ) and the timescale of chemistry (_Ch): the first Damköhler number ( _t /_Ch). Special attention is given to the calculation of the covariance between NO and O₃ normalized by the product of their means, the so-called intensity of segregation (I _S ). This parameter quantifies the state of mixing of two chemical species.The intensity of segregation is calculated as a function of the flux of NO and the first Damköhler number. The model results presented illustrate the importance of taking the effect of turbulence on chemical reactions into account for higher values of the NO flux, for values of the ratio O₃/NO larger than 12.5 and for values of the ratio _t /_CH larger than 0.1. For such cases, the effective reaction rates are slower than if the chemical species are assumed to be uniformly mixed.     

On the determination of the height of the Ekman boundary layer

The heighth of the Ekman turbulent boundary layer determined by the momentum flux profile is estimated with the aid of considerations of similarity and an analysis of the dynamic equations. Asymptotic formulae have been obtained showing that, with increasing instability,h increases as ¦¦^1/2 (where is the non-dimensional stratification parameter); with increasing stability, on the other hand,h decreases as ^–1/2. For comparison, a simple estimate of the boundary-layer heighth _u determined by the velocity profile is given. As is shown, in unstable stratification,h _u behaves asymptotically as ¦¦^–1, i.e., in a manner entirely different from that ofh .     

A dynamical analysis of the drag conditions at the sea surface

When applied to a sea surface, shortcomings are noted for the ordinary classification of drag conditions at rigid underlying surfaces according to the Reynolds roughness number Re _s . It is shown that in the case of mobile underlying surfaces, it would be more natural to use the dynamical classification of drag conditions according to the order of magnitude of the ratio ( = /) of the momentum flux toward the waves ( _w) to the viscous momentum flux through the surface ( _w). The relevant estimates of for the main stages of development of the wind waves indicate that the observed values of the drag coefficient of the sea surface correspond to the case of underdeveloped roughness.     

Modelling radiation quantities and photolysis frequencies in the troposphere

STAR (System for Transfer of Atmospheric Radiation) was developed to calculate accurately and efficiently the irradiance, the actinic flux, and the radiance in the troposphere. Additionally a very efficient calculation scheme to computer photolysis frequencies for 21 different gases was evolved. STAR includes representative data bases for atmospheric constituents, especially aerosol particles. With this model package a sensitivity study of the influence of different parameter on photolysis frequencies in particular of O₃ to Singlet D oxygen atoms, of NO₂, and of HCHO was performed. The results show the quantitative effects of the influence of the solar zenith angle, the ozone concentration and vertical profile, the aerosol particles, the surface albedo, the temperature, the pressure, the concentration of NO₂, and different types of clouds on the photolysis frequencies.Notation I _A(, ) actinic flux - I _H(, ) irradiance - L(, , , ) radiance - wavelength - azimuth angle - cosine of zenith angle - _s cosine of solar zenith angle - optical depth - _s scattering coefficient - _c extinction coefficient - _o single scattering albedo - p _mix mixed phase function - g _mix mixed asymmetry factor - J _gas photolysis frequency     

Wind stress over water waves: Field experiments on lake of Geneva

Summary A fixed platform (Fig.3), installed 100 m from the shoreline in 3 m water depth, was instrumented with velocity, temperature and wave-height sensors. 132 data (10 minutes averages) were analysed to calculate the wind stress; from these, 99 data were used to investigate the vertical distribution of the wind stress; all data are presented with Table 1.It was postulated that the total stress, _t being constant with height, is made up additively of two components, the wave-supporting stress, _w , and the turbulent stress, _c ; see Eq. 1. The vertical distribution of these two components is shown schematically in Fig. 1.The total stress, _t , evaluated outside the zone of wave influence, is given in the classical way with Fig. 4. The wave-supporting stress, _w (z), was evaluated from the data according to a relation proposed by Kitaigorodskii et al. (1984); it is given with Fig. 5. A height-dependency is clearly evident. The turbulent stress _c (z), was evaluated with data of the velocity gradient; it is given with Fig. 6. A height-dependency is not evident.The field data from the lake of Geneva give evidence that the additive relation of Eq. 1 seems to be justified.With 6 Figures     

A comparison of airborne eddy correlation and bulk aerodynamic methods for ocean-air turbulent fluxes during cold-air outbreaks

Four bulk schemes (LKB, FG, D and DB), with the flux-profile relationships of Liuet al. (1979), Francey and Garratt (1981), Dyer (1974), and Dyer and Bradley (1982), are derived from the viscous interfacial-sublayer model of Liuet al. These schemes, with stability-dependent transfer coefficients, are then tested against the eddy-correlation fluxes measured at the 50 m flight level above the western Atlantic Ocean during cold-air outbreaks. The bulk fluxes of momentum (), sensible heat (H), and latent heat (E) are found to increase with various von Kármán constants (k _M for k _H forH, andk _E forE). Except that the LKB scheme overestimates by 28% (46Wm^–2), on the average, the fluxes estimated by the four bulk schemes appear to be in fairly good agreement with those of the eddy correlation method (magnitudes of biases within 10% for , 17% forH, and 13% forE). The results suggest that the overall fluxes and surface-layer scaling parameters are best estimated by FG and thatk _H <k _E . On the average, the FG scheme underestimates by 10% (0.032N m^–2) andE by 4% (12Wm^–2), and overestimatesH by 0.3% (0.5W m^–2). The equivalent neutral transfer coefficients at 10 m height of the FG scheme compare well with some schemes of those tested by Blanc (1985).The relative importance of various von Kármán constants, dimensionless gradients and roughness lengths to the oceanic transfer coefficients is assessed. The dependence of transfer coefficients on wind speeds and roughness lengths is discussed. The transfer coefficients for andE agree excellently between LKB and FG. However, the ratio of the coefficient forH of LKB to that of FG, increasing with decreasing stability, is very sensitive to stability at low winds, but approaches the neutral value of 1.25 at high winds.     

Statistical properties of 49 years of rainfall rate events

Summary The paper presents a detailed investigation of magnitudes and properties pertaining to the large population of rainfall events recorded during 49 years at the Fabra Observatory in Barcelona. The study includes a statistical analysis of event durationT and rainfall quantityQ together with the statistical rainfall rate parameters: , ²(R) and maximum Rm within an event. The decorrelation time is also analysed. It is found thatQ, T, and can be well modelled by a lognormal distribution, but , ² and Rm are only so for a limited range of precipitation ratesR. The regression analysis between pairs of logarithms of the magnitudes investigated is generally good and a regression coefficient is often better than 0.9. Comparison with published work is also carried out. An attempt is made to discriminate between heavy and non-heavy rainfall rates, and the 50 mm/h threshold is used for the study. The twelve-monthly running average indicates that the rainfall amount has a small increasing trend over the fifty year period. However, this trend is reversed when considering heavy rains. Finally, the return period in years to exceeding a thresholdR within an event is also investigated and the distribution of the population of annual extremes is found to be Gumbel II.With 6 Figures     

Disagreements between gradient-diffusion and Lagrangian stochastic dispersion models,even for sources near the ground

It is well known that if turbulent mass convection is modelled as diffusion, errors result unless trajectories from the source (ath) to the point of observation (z _p ) comprise many statistically-independent segments (Taylor, 1921). We show that this is not guaranteed merely by the Lagrangian timescale () at the source being small (e.g., source at ground), but that a better criterion istmax[(h), (z _p )], wheret is a typical travel time toz _p .     

Arctic aerosol and Arctic climate: Results from an energy budget model

An energy budget model is used to study the effect on Arctic climate of optically active aerosol in the Arctic atmosphere. The dependence of the change in surface temperature on the vertical distribution of the aerosol and on the radiative properties of the aerosol-free atmosphere, the Arctic surface, and the aerosol, itself, are calculated. An extensive sensitivity analysis is performed to assess the degree to which the results of the model are dependent upon the assumptions underlying it.List of Symbols Used I ₀ Solar flux at the top of the Arctic Atmosphere (Arctic here means 70° N latitude to the pole) - a _S Surface albedo of the Arctic (a _S ^c is the value of surface albedo at which the sign of the surface temperature perturbation changes) - Reflection coefficient of the aerosol-free Arctic atmosphere - Absorption coefficient of the aerosol-free Arctic atmosphere - Transmission coefficient of the aerosol-free Arctic atmosphere - RI ₀ Total flux of sunlight reflected from the Arctic - A _A I ₀ Total flux of sunlight absorbed in the Arctic atmosphere - A _S I ₀ Total flux of sunlight absorbed at the Arctic surface - A _aer I ₀ Total flux of sunlight absorbed in the Arctic aerosol - Q _A Net atmospheric flow of energy, per unit of Arctic surface area, north across 70° N latitude - Q _S Net oceanic flow of energy, per unit of Arctic surface area, north across 70° N latitude - E Convective plus latent heat fluxes from surface to atmosphere - F _A Net flow of energy to the Arctic atmosphere - F _S Net flow of energy to the Arctic surface - T _A An effective temperature of the Arctic atmosphere - T _S Surface temperature of the Arctic - w Single-scattering albedo of the aerosol - t Optical depth of the aerosol - g Fraction of incident radiation scattered forward by the aerosol - Reflection coefficient of the aerosol - Absorption coefficient of the aerosol - Transmission coefficient of the aerosol - p,q Number of atmospheric layers and the inverse of the fraction of incident IR absorbed in each layer in the energy budget model - F,G,H Measures of the amount of IR-active atmosphere above the surface, the aerosol, and the clouds     

Similarity of atmospheric Reynolds shear stress and heat flux fluctuations over a rough surface

Horizontal u and vertical w velocity fluctuations have been measured together with temperature fluctuations in the atmospheric surface layer, at a small height above a wheat crop canopy. Marginal probability density functions are presented for both individual fluctuations u, w, and for the instantaneous Reynolds stress uw, and heat fluxes w and u. Probability density functions of the velocity fluctuations deviate less significantly from the Gaussian form than the probability density of temperature. There appears to be closer similarity between statistics of the instantaneous heat fluxes than between the momentum flux and either of the heat fluxes investigated. The mean momentum flux receives equal contributions from the events referred to as ejections and sweeps in laboratory boundary layers. Sweeps provide the largest contribution to the heat fluxes.     

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.     

Basis for integration of conventional observations of cloud into global nephanalyses

Water vapour is the largest and, radiatively, the most important trace gas in the Earth's atmosphere. Cloud amount and cloud optical depth feedbacks are, as yet, poorly understood and improvements in model parameterization schemes await an adequate observational data base. Satellite retrievals do not, and will not for some time, provide more than snapshot (a few months to a few years) records. Conventional surface-based observations of cloudiness could complement the global coverage offered by satellite retrievals if a sound relationship between the two observational measurements could be found. Observations underline the importance of the vertical dimension of clouds which affects the observational geometries from satellites and the surface. A new basis for the relationship between the (vertical) earthview of cloud amount and the (whole dome) skycover of cloud amount has been sought. Over four and a half thousand all-sky camera photographs, representing a considerable range of seasonal and climatological conditions, have been analyzed to give rise to a database from which predictive relationships for earthview, E, and skycover, S, have been established.Cubic functions are the most soundly based both physically and empirically. We find: S=0.647+2.192E–0.461E ²+0.037E ³ and E=-0.001+0.082S+0.269S ²–0.019S ³ for the prediction of skycover from earthview and earthview from skycover, respectively. If earthview is required from skycover observations then ES could be used with little additional error. Hence, conventional surface observations of skycover could be compared directly with satellite-derived earthview values. More importantly, these results do not support the widespread assumption that conventional (surface) observations of cloud amount always exceed the earthview could retrieval. Furthermore, climate model predictions of total cloud amount may also be interpreted via these relationships. GCM-predicted layer cloud amounts can be synthesized into modelled E values using the random overlap formula and hence it is possible to construct modelled S values which are directly comparable with conventionally observed cloud climatologies. The baseline observation of skycover by clouds therefore provides a valuable validation tool for both satellite programmes and climate models.     

The nature and significance of the feedback of changes in terrestrial vegetation on global atmospheric and climatic change

The potential feedback on global atmospheric and climate change of climate-driven changes in terrestrial vegetation is examined by systematically relating the surface exchanges of energy, mass and momentum to two dimensions of vegetation, structure and taxonomy, such that the significance of climate driven changes in these characteristics can be assessed. A detailed quantitative understanding of this feedback is an important prerequisite to realistic and dynamic representations of the Earth's surface within general circulation and biological models (GCMs and GBMs). Without realistic representations of terrestrial vegetation within these models, any forecasts of future climates by these models must be suspect.Several general conclusions are drawn. The first is that the indirect feedbacks, those associated with the clouds and aerosols of the planetary boundary layer, appear to be very powerful but as yet their behaviour and connections with the underlying surface are both poorly understood and captured within GCMs.The physical structure of vegetation, the disposition of biomass in 3-D, is the characteristic that most strongly influences the exchange of momentum (via aerodynamic roughness) and solar radiation (via albedo). Vegetation structure and species composition determine the most important of the mass exchanges, evapotranspiration. Of all of the surface exchanges, the parameterization of evapotranspiration (E ) and the simulation of the water balance over time is the most critical.Lastly, the problems of scaling and spatial heterogeneity, the sub-grid variability of the modellers, looms as a difficult, but not insoluble, problem. It remains a critical problem however, and the detailed parameterization of the various big leaf models stands in absurd contrast to the simplistic generalization of the spatial heterogeneity of terrestrial landscapes.Plant ecologists can contribute to the task of improving the representation of vegetated landscapes within GCMs. There is need to simply and unify the way in which vegetation can be grouped at landscape scales. A classification that is based on function rather than phylogeny is required. The definition of Vegetation Functional Types (VFTs) would expedite research on both the impact of, and feedback on, climate change.     

Numerical integration errors in calculated tropospheric photodissociation rate coefficients

Tropospheric photodissociation rate coefficients (J values) were calculated for NO₂, O₃, HNO₂, CH₂O, and CH₃CHO using high spectral resolution (0.1 mm wavelength increments), and compared to the J values obtained with numerically degraded resolution (=1, 2, 4, 6, 8, and 10 nm, and several commonly used nonuniform grids). Depending on the molecule, substantial errors can be introduced by the larger increments. Thus for =10 nm, errors are less than 1% for NO₂, less than 2% for HNO₂, +6.5% to -16% for CH₂O, -6.9% to +24% for CH₃CHO, and -24% to +110% for O₃. The errors for CH₂O arise from the fine structure of its absorption spectrum, and are prevalently negative (underestimate of J). The errors for O₃, and to a lesser extent for CH₃CHO, arise mainly from under-resolving the overlap of the molecular action spectrum and the tropospheric actinic flux in the wavelength region of stratospheric ozone attenuation. The sign of those errors depends on whether the actinic flux is averaged onto the grid before or after the radiative transfer calculation. In all cases studied, grids with 2 nm produced errors no larger than 5%.     

The greening of the McGill Paleoclimate Model. Part I: Improved land surface scheme with vegetation dynamics

The formulation of a new land surface scheme (LSS) with vegetation dynamics for coupling to the McGill Paleoclimate Model (MPM) is presented. This LSS has the following notable improvements over the old version: (1) parameterization of deciduous and evergreen trees by using the models climatology and the output of the dynamic global vegetation model, VECODE (Brovkin et al. in Ecological Modelling 101:251–261 (1997), Global Biogeochemical Cycles 16(4):1139, (2002)); (2) parameterization of tree leaf budburst and leaf drop by using the models climatology; (3) parameterization of the seasonal cycle of the grass leaf area index; (4) parameterization of the seasonal cycle of tree leaf area index by using the time-dependent growth of the leaves; (5) calculation of land surface albedo by using vegetation-related parameters, snow depth and the models climatology. The results show considerable improvement of the models simulation of the present-day climate as compared with that simulated in the original physically-based MPM. In particular, the strong seasonality of terrestrial vegetation and the associated land surface albedo variations are in good agreement with several satellite observations of these quantities. The application of this new version of the MPM (the green MPM) to Holocene millennial-scale climate changes is described in a companion paper, Part II.

---

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.