Spectra of CO<Subscript>2</Subscript> and Water Vapour in the Marine Atmospheric Surface Layer

Spectra of CO₂ and water vapour fluctuations from measurements made in the marine atmospheric surface layer have been analyzed. A normalization of spectra based on Monin–Obukhov similarity theory, originally developed for wind speed and temperature, has been successfully extended also to CO₂ and humidity spectra. The normalized CO₂ spectra were observed to have somewhat larger contributions from low frequencies compared to humidity spectra during unstable stratification. However, overall, the CO₂ and humidity spectra showed good agreement as did the cospectra of vertical velocity with water vapour and CO₂ respectively. During stable stratification the spectra and cospectra displayed a well-defined spectral gap separating the mesoscale and small-scale turbulent fluctuations. Two-dimensional turbulence was suggested as a possible source for the mesoscale fluctuations, which in combination with wave activity in the vertical wind is likely to explain the increase in the cospectral energy for the corresponding frequency range. Prior to the analysis the turbulence time series of the density measurements were converted to time series of mixing ratios relative to dry air. Some differences were observed when the spectra based on the original density measurements were compared to the spectra based on the mixing ratio time series. It is thus recommended to always convert the density time series to mixing ratio before performing spectral analysis.     

Turbulence characteristics in the atmospheric surface layer during summer monsoon of 1997 over a semi-arid location in India

A land surface processes experiment (LASPEX) was conducted in the semi-arid region of Northwest India during January 1997–February 1998. Analysis of turbulent components of wind and air temperature collected in the surface layer (SL) at Anand (22°35′N, 72°55′E) during the Indian summer monsoon season from June to September 1997 is presented. Turbulent fluctuation of wind components and air temperature observed at Anand varied as a function of terrain features and stability of the surface layer. Under neutral conditions, the standard deviation of vertical velocity (σ _w ) and temperature (σ _T ) were normalized using respective surface layer scaling parameter u _* and T _* which fitted the expressions σ _w /u _* = 1.25 and σ _T /T _* ≈ 4. Micrometeorological spectrum of wind and temperature at 5 m above ground level (AGL) at Anand showed peaks at time scale of 1–3 min at the low-frequency end. The inertial sub-range characteristics (?2/3 slope) of the spectrum are exhibited mostly. However, in some occasions, slope of ?1 denoting brown noise was depicted by the wind and temperature spectrum, which indicated anisotropy in turbulence.     

Simplification and Validation of a Spectral-Tensor Model for Turbulence Including Atmospheric Stability

A spectral-tensor model of non-neutral, atmospheric-boundary-layer turbulence is evaluated using Eulerian statistics from single-point measurements of the wind speed and temperature at heights up to 100 m, assuming constant vertical gradients of mean wind speed and temperature. The model has been previously described in terms of the dissipation rate \(\epsilon \), the length scale of energy-containing eddies \(\mathcal {L}\), a turbulence anisotropy parameter \(\varGamma \), the Richardson number Ri, and the normalized rate of destruction of temperature variance \(\eta _\theta \equiv \epsilon _\theta /\epsilon \). Here, the latter two parameters are collapsed into a single atmospheric stability parameter z / L using Monin–Obukhov similarity theory, where z is the height above the Earth’s surface, and L is the Obukhov length corresponding to \(\{Ri,\eta _\theta \}\). Model outputs of the one-dimensional velocity spectra, as well as cospectra of the streamwise and/or vertical velocity components, and/or temperature, and cross-spectra for the spatial separation of all three velocity components and temperature, are compared with measurements. As a function of the four model parameters, spectra and cospectra are reproduced quite well, but horizontal temperature fluxes are slightly underestimated in stable conditions. In moderately unstable stratification, our model reproduces spectra only up to a scale \(\sim \) 1 km. The model also overestimates coherences for vertical separations, but is less severe in unstable than in stable cases.     

Numerical Simulations of the Microscale Heterogeneities of Turbulence Observed on a Complex Site

We simulate the microscale heterogeneities of turbulent variables observed at a complex site for different wind directions. The atmospheric computational fluid dynamics (CFD) results are compared with an ensemble of 36 months of data collected at the experimental site SIRTA “Site Instrumental de Recherche par Télédétection Atmosphérique”, located near Paris (France) in a semi-urban environment. The experimental data show that the normalized turbulent kinetic energy (TKE) k/U ² (where k is TKE and U is the wind speed) at 10-m height, for two different locations, is highly dependent on wind direction and strongly influenced by trees. These measurements show a strong increase of the normalized TKE downstream of the forest canopies with a large variability within the 36-month period in part due to the variation of the tree foliage. The numerical simulations are carried out using the CFD code Code_Saturne with the standard k?ε closure, in neutral stratification. The buildings are taken into account explicitly in the mesh and the forested areas are modelled with two approaches: the classical roughness wall law and a drag porosity. A comparison has been performed between the calculated values and the median of measured values of the normalized TKE and the normalized friction velocity, for each wind sector of 10°. A very good agreement is obtained with the drag porosity model, whereas the classical roughness law leads to a strong underestimation downstream of the forested areas. However, this large improvement of the results using the drag porosity model can only be obtained with a refinement of the grid, especially in forested areas, and an accurate land-use map.     

Methods for Estimating Air–Sea Fluxes of CO<Subscript>2</Subscript> Using High-Frequency Measurements

The most direct method for flux estimation uses eddy covariance, which is also the most commonly used method for land-based measurements of surface fluxes. Moving platforms are frequently used to make measurements over the sea, in which case motion can disturb the measurements. An alternative method for flux estimation should be considered if the effects of platform motion cannot be properly corrected for. Three methods for estimating CO₂ fluxes are studied here: the eddy-covariance, the inertial-dissipation, and the cospectral-peak methods. High-frequency measurements made at the land-based Östergarnsholm marine station in the Baltic Sea and measurements made from a ship during the Galathea 3 expedition are used. The Kolmogorov constant for CO₂, used in the inertial-dissipation method, is estimated to be 0.68 and is determined using direct flux measurements made at the Östergarnsholm site. The cospectral-peak method, originally developed for neutral stratification, is modified to be applicable in all stratifications. With these modifications, the CO₂ fluxes estimated using the three methods agree well. Using data from the Östergarnsholm site, the mean absolute error between the eddy-covariance and inertial-dissipation methods is 0.25 μmol  m^?2 s^?1. The corresponding mean absolute error between the eddy-covariance and cospectral-peak methods is 0.26 μmol m^?2 s^?1, while between the inertial-dissipation and cospectral-peak methods it is 0.14 μmol m^?2 s^?1.     

The Influences of Thermodynamic Characteristics on Aerodynamic Roughness Length over Land Surface

It has previously been shown that aerodynamic roughness length changes significantly along with nearsurface atmospheric thermodynamic state; however, at present, this phenomenon remains poorly understood, and very little research concerning this topic has been conducted. In this paper, by using the data of different underlying surfaces provided by the Experimental Co-observation and Integral Research in Semi-arid and Arid Regions over North China, aerodynamic roughness length (z₀) values in stable, neutral, and unstable atmospheric stratifications are compared with one another, and the relationship between z₀ and atmospheric thermodynamic stability (ζ) is analyzed. It is found that z₀ shows great differences among the stable, neutral, and unstable atmospheric thermodynamic states, with the difference in z₀ values between the fully thermodynamic stable condition and the neutral condition reaching 60% of the mean z₀. Furthermore, for the wind speed range in which the wind data are less sensitive to z₀, the surface z₀ changes more significantly with ζ, and is highly correlated with both the Monin-Obukhov stability (ζ₀) and the overall Richardson number (R_ib), with both of their correlation coefficients greater than 0.71 and 0.47 in the stable and unstable atmospheric stratification, respectively. The empirical relation fitted with the experimental observations is quite consistent with the Zilitinkevich theoretical relation in the stable atmosphere, but the two are quite distinct and even show opposite variation tendencies in the unstable atmosphere. In application, however, verification of the empirical fitted relations by using the experimental data finds that the fitted relation is slightly more applicable than the Zilitinkevich theoretical relation in stable atmospheric stratification, but it is much more suitable than the Zilitinkevich relation in unstable atmospheric stratification.     

频散效应对近地面层风速廓线的影响

应用含湍流频散效应的近地面层的运动方程求解了不同层结下的风速廓线，着重分析了湍流频散效应对近地面层平稳运动的影响。分析指出：湍流的频散效应对经典的幂律廓线一对数修正，该修正在不稳定层结时比稳定层时明显；利用相似理论也得到了该常数。     

Stably stratified flow in a marine atmospheric surface layer

Data from the marine atmospheric surface layer have been analysed. The data set consists of about two weeks with tower measurements up to 31 m of mean profiles of wind, temperature, and humidity, together with 20 Hz turbulence data. Mean wind, temperature, and humidity profiles up to 2000 m are also available from pibal trackings and radio soundings. Wave height was measured at 2 Hz, using an inverted echo-sounder.It was found from pibal wind profiles that low level jets were present during 2/3 of the measurements, having their maxima in the height interval 40 to 300 m. Here only data from the remaining 1/3 of the measurements, without low level jets, have been analysed.Non-dimensional wind and temperature gradients agree with results over homogeneous land surfaces as regards stability dependence during stable conditions that prevailed during this experiment. Linear regression gave _m = 1 + 6.8z/L and _m = 1 + 8.3z/L. No significant sea wave influence was found. The same was vrue for me dimensionless standard deviations of the three wind components, except for the vertical component. The expected wind speed dependence was found for the neutral drag coefficient, givingC _dN = 0.109U + 0.33 at 10 m, and a dependence on the wave parameter,C/u _*, was confirmed. Note, however, that the data set was restricted to low and moderate wind speeds and that stratification was mainly stable.Power spectra, non-dimensionalized according to suface-layer theories, do not follow the expected stability dependence. It was shown that this may be a consequence of the presence of gravity waves in the stable marine boundary layer. Indicators of gravity waves were found in most runs. The TKE budget agrees with findings over homogeneous land areas. The pressure transport term was found to be a source of energy also for near neutral conditions.     

Turbulence spectra of CO2, water vapor,temperature and wind velocity fluctuations over a crop surface

A prototype rapid-response CO₂ sensor was used in conjunction with a Lyman-alpha hygrometer, fine-wire thermocouples and a three-dimensional drag anemometer to measure CO₂, humidity, temperature and wind velocity fluctuations. Measurements were made over a soybean crop grown on relatively flat terrain near Mead, Nebraska.Temperature, humidity and CO₂ spectra measured under near neutral conditions were most similar in shape while longitudinal velocity (U) spectra appeared to be somewhat broader. Peaks occurred around f = 0.02 to 0.06 (where f is the non-dimensional frequency) in all spectra except for vertical velocity which had a peak near f = 0.5. As thermal stability changed from neutral to stable, spectra exhibited lower peaks, were narrower in shape and were shifted towards higher frequencies. Opposite behavior was observed with increasing instability.Cospectra for momentum, heat, water vapor and CO₂ measured during neutral conditions had similar peak frequencies (near f = 0.15). Among the cospectra studied, CO₂ and water vapor cospectra had the greatest similarity in shape.Published as Paper No. 7481, Journal Series, Nebraska Agricultural Experiment Station. The work reported here was conducted under Nebraska Agricultural Experiment Station Project 27-003 and Regional Research Project 11-33.Post-doctoral Research Associate and Professor, respectively, Center for Agricultural Meteorology and Climatology, Institute of Agriculture and Natural Resources, University of Nebraska, Lincoln, NE 68583-0728.     

Turbulent transport of carbon dioxide over a paddy field

Turbulent fluctuations in CO₂ concentrations over a paddy field are measured by a fastresponse device with an open sensing path. This IR device coupled with a sonic anemometer constitutes an eddy correlation instrument to measure CO₂ fluxes. Three experiments were conducted in the surface layer over paddy 90 cm high. The stability (z – d)/L ranged from -0.14 to 0.20, where L denotes the Monin-Obukhov length.CO₂ power spectra show the range of applicability of the -2/3 power law to be between f = 0.2 and f = 2, where f is the frequency normalized by wind speed and height. The cospectral estimate between CO₂ and vertical component of wind speed ranging from f = 0.005 to f = 2 shows a peak at about f = 0.15 under near-neutral stratification.Hourly means of CO₂ flux measured by the eddy correlation method increase with intensity of net radiation. The maximum value of downward flux of CO₂ rises to 0.6 mg cm^-2 hr^-1 over the paddy field at the stage of ear emergence.Some turbulence statistics relating to the CO₂ transport are evaluated: the correlation coefficient between CO₂ and vertical velocity is about -0.3, and that between CO₂ and humidity attains -0.7 ~ -0.8 under unstable stratification; nondimensional gradients _c for CO₂ and _m for wind speed are 0.89 and 0.99, respectively.     

Low-frequency and high-frequency changes in temperature and effective humidity during the Holocene in south-central Sweden: implications for atmospheric and oceanic forcings of climate

An integrated use of independent palaeoclimatological proxy techniques that reflect different components of the climate system provides a potential key for functional analysis of past climate changes. Here we report a 10,000 year quantitative record of annual mean temperature (T _ann), based on pollen-climate transfer functions and pollen-stratigraphical data from Lake Flarken, south-central Sweden. The pollen-based temperature reconstruction is compared with a reconstruction of effective humidity, as reflected by a δ¹⁸O record obtained on stratigraphy of lacustrine carbonates from Lake Igelsjön, c. 10 km from Lake Flarken, which gives evidence of pronounced changes in effective humidity. The relatively low T _ann, and high effective humidity as reflected by a low evaporation/inflow ratio suggest a maritime early Holocene climate (10,000–8,300 cal year BP), seemingly incompatible with the highly seasonal solar insolation configuration. We argue that the maritime climate was due to the stronger-than-present zonal flow, enhanced by the high early Holocene sea-surface temperatures in the North Atlantic. The maritime climate mode was disrupted by the abrupt cold event at 8,200 cal year BP, followed at 8,000 cal year BP by a stable Holocene Thermal Maximum. The latter was characterized by T _ann values about 2.5°C higher than at present and markedly dry conditions, indicative of stable summer-time anti-cyclonic circulation, possibly corresponding with modern blocking anticyclonic conditions. The last 4,300 year period is characterized by an increasingly cold, moist, and unstable climate. The results demonstrate the value of combining two independent palaeoclimatic proxies in enhancing the reliability, generality, and interpretability of the palaeoclimatic results. Further methodological refinements especially in resolving past seasonal climatic contrasts are needed to better understand the role of different forcing factors in driving millennial-scale climate dynamics.     

Comparison of a simple logarithmic and equivalent neutral wind approaches for converting buoy-measured wind speed to the standard height: special emphasis to North Indian Ocean

The difference between the transferred wind speed to 10-m height based on the equivalent neutral wind approach (U _n) and the logarithmic approach (U _log) is studied using in situ observations from the Indian, Pacific, and Atlantic Oceans, with special emphasis given to the North Indian Ocean. The study included U _n ? U _log variations with pressure, relative humidity, wind speed, air temperature, and sea surface temperature (SST). U _n ? U _log variation with respect to air temperature (T _a) reveals that U _n ? U _log is out of phase with air temperature. Further analysis found that U _n ? U _log is in phase with SST (T _s) ? T _a and varies between ?1.0 and 1.0 m/s over the North Indian Ocean, while for the rest of the Oceans, it is between ?0.3 and 0.8 m/s. This higher magnitude of U _n ? U _log over the North Indian Ocean is due to the higher range of T _s ? T _a (?4 to 6 °C) in the North Indian Ocean. Associated physical processes suggested that the roughness length and friction velocity dependence on the air–sea temperature difference contributes to the U _n ? U _log difference. The study is further extended to evaluate the behavior of U _n ? U _log under cyclonic conditions (winds between 15 and 30 m/s), and it was found that the magnitude of U_n ? U _log varies 0.5–1.5 m/s under the cyclonic wind conditions. The increasing difference with the wind speed is due to the increase in the momentum transfer coefficient with wind speed, which modifies the friction velocity significantly, resulting in U _n higher than U _log. Thus, under higher wind conditions, U _n ? U _log can contribute up to half the retrieval error (5 % of the wind speed magnitude) to the satellite validation exercise.     

The Critical Richardson Number and Limits of Applicability of Local Similarity Theory in the Stable Boundary Layer

Measurements of atmospheric turbulence made over the Arctic pack ice during the Surface Heat Budget of the Arctic Ocean experiment (SHEBA) are used to determine the limits of applicability of Monin–Obukhov similarity theory (in the local scaling formulation) in the stable atmospheric boundary layer. Based on the spectral analysis of wind velocity and air temperature fluctuations, it is shown that, when both the gradient Richardson number, Ri, and the flux Richardson number, Rf, exceed a ‘critical value’ of about 0.20–0.25, the inertial subrange associated with the Richardson–Kolmogorov cascade dies out and vertical turbulent fluxes become small. Some small-scale turbulence survives even in this supercritical regime, but this is non-Kolmogorov turbulence, and it decays rapidly with further increasing stability. Similarity theory is based on the turbulent fluxes in the high-frequency part of the spectra that are associated with energy-containing/flux-carrying eddies. Spectral densities in this high-frequency band diminish as the Richardson–Kolmogorov energy cascade weakens; therefore, the applicability of local Monin–Obukhov similarity theory in stable conditions is limited by the inequalities Ri < Ri _cr and Rf < Rf _cr. However, it is found that Rf _cr  =  0.20–0.25 is a primary threshold for applicability. Applying this prerequisite shows that the data follow classical Monin–Obukhov local z-less predictions after the irrelevant cases (turbulence without the Richardson–Kolmogorov cascade) have been filtered out.     

Air–Sea Interaction Features in the Baltic Sea and at a Pacific Trade-Wind Site: An Inter-comparison Study

A systematic comparison of wind profiles and momentum exchange at a trade wind site outside Oahu, Hawaii and corresponding data from the Baltic Sea is presented. The trade wind data are to a very high degree swell dominated, whereas the Baltic Sea data include a more varied assortment of wave conditions, ranging from a pure growing sea to swell. In the trade wind region swell waves travel predominantly in the wind direction, while in the Baltic, significant cross-wind swells are also present. Showing the drag coefficient as a function of the 10-m wind speed demonstrates striking differences for unstable conditions with swell for the wind-speed range 2 m s^?1 < U ₁₀ < 7 m s^?1, where the trade-wind site drag values are significantly larger than the corresponding Baltic Sea values. In striking contrast to this disagreement, other features studied are surprisingly similar between the two sites. Thus, exactly as found previously in Baltic Sea studies during unstable conditions and swell, the wind profile in light winds (3 m s^?1) shows a wind maximum at around 7–8 m above the water, with close to constant wind speed above. Also, for slightly higher wind speeds (4 m s^?1 < U ₁₀ < 7 m s^?1), the similarity between wind profiles is striking, with a strong wind-speed increase below a height of about 7–8 m followed by a layer of virtually constant wind speed above. A consequence of these wind-profile features is that Monin–Obukhov similarity is no longer valid. At the trade-wind site this was observed to be the case even for wind speeds as high as 10 m s^?1. The turbulence kinetic energy budget was evaluated for four cases of 8–16 30- min periods at the trade-wind site, giving results that agree very well with corresponding figures from the Baltic Sea.     

Evaluation of the antarctic surface wind climate from ERA reanalyses and RACMO2/ANT simulations based on automatic weather stations

A continental scale evaluation of Antarctic surface winds is presented from global ERA-40 and ERA-Interim reanalyses and RACMO2/ANT regional climate model at 55 and 27 km horizontal resolution, based on a comparison with observational data from 115 automatic weather stations (AWS). The Antarctic surface wind climate can be classified based on the Weibull shape factor k _w . Very high values (k _w  > 3) are found in the interior plateaus, typical of very uniform katabatic-dominated winds with high directional constancy. In the coast and all over the Antarctic Peninsula the shape factors are similar to the ones found in mid-latitudes (k _w  < 3) typical of synoptically dominated wind climates. The Weibull shape parameter is systematically overpredicted by ERA reanalyses. This is partly corrected by RACMO2/ANT simulations which introduce more wind speed variability in complex terrain areas. A significant improvement is observed in the performance of ERA-Interim over ERA-40, with an overall decrease of 14 % in normalized mean absolute error. In escarpment and coastal areas, where the terrain gets rugged and katabatic winds are further intensified in confluence zones, ERA-Interim bias can be as high as 10 m s^?1. These large deviations are partly corrected by the regional climate model. Given that RACMO2/ANT is an independent simulation of the near-surface wind speed climate, as it is not driven by observations, it compares very well to the ERA-Interim and AWS-115 datasets.     

Observational evidence for the Monin-Obukhov similarity under all stability conditions

Data collected in the surface layer in a northern suburban area of Nanjing from 15 November to 29 December 2007 were analyzed to examine the Monin-Obukhov similarity for describing the turbulent fluctu- ations of 3D winds under all stability conditions and to obtain the turbulence characteristics under different weather conditions. The results show that the dimensionless standard deviations of turbulent velocity com- ponents (σ u /u* , σ v /u* , σ w /u * ) and dimensionless turbulent kinetic energy (TKE) can be well described by "1/3" power law relationships under stable, neutral, and unstable conditions, with σ u /u * > σ v /u * > σ w /u* . Land use and land cover changes mainly impact dimensionless standard deviations of horizontal component fluctuations, but they have very little on those of the vertical component. The dimensionless standard devi- ations of wind components and dimensionless TKE are remarkably affected by different weather conditions; the deviations of horizontal wind component and dimensionless TKE present fog day > clear sky > overcast > cloudy; the trend of the vertical wind component is the reverse. The surface drag coefficient at a Nan- jing suburban measurement site during the observation period was obviously higher than at other reported plains and plateau areas, and was approximately one order larger in magnitude than the reported plains areas. Dimensionless standard deviation of temperature declined with increasing |z /L| with an approximate "-1/3" slope in unstable stratification and "-2/3" slope in stable stratification.     

Trends in middle- and upper-level tropospheric humidity from NCEP reanalysis data

The National Centers for Environmental Prediction (NCEP) reanalysis data on tropospheric humidity are examined for the period 1973 to 2007. It is accepted that radiosonde-derived humidity data must be treated with great caution, particularly at altitudes above the 500 hPa pressure level. With that caveat, the face-value 35-year trend in zonal-average annual-average specific humidity q is significantly negative at all altitudes above 850 hPa (roughly the top of the convective boundary layer) in the tropics and southern midlatitudes and at altitudes above 600 hPa in the northern midlatitudes. It is significantly positive below 850 hPa in all three zones, as might be expected in a mixed layer with rising temperatures over a moist surface. The results are qualitatively consistent with trends in NCEP atmospheric temperatures (which must also be treated with great caution) that show an increase in the stability of the convective boundary layer as the global temperature has risen over the period. The upper-level negative trends in q are inconsistent with climate-model calculations and are largely (but not completely) inconsistent with satellite data. Water vapor feedback in climate models is positive mainly because of their roughly constant relative humidity (i.e., increasing q) in the mid-to-upper troposphere as the planet warms. Negative trends in q as found in the NCEP data would imply that long-term water vapor feedback is negative—that it would reduce rather than amplify the response of the climate system to external forcing such as that from increasing atmospheric CO₂. In this context, it is important to establish what (if any) aspects of the observed trends survive detailed examination of the impact of past changes of radiosonde instrumentation and protocol within the various international networks.     

Estimation of atmospheric surface layer parameters and numerical simulation using MM5 at Mangalore, West Coast of India

Atmospheric surface layer meteorological observations obtained from 20-m-high meteorological tower at Mangalore, situated along the west coast of India are used to estimate the surface layer scaling parameters of roughness length (z _o) and drag coefficient (C _D), surface layer fluxes of sensible heat and momentum. These parameters are computed using the simple flux–profile relationships under the framework of Monin–Obukhov (M–O) similarity theory. The estimated values of z _o are higher (1.35–1.54 m) than the values reported in the literature (>0.4–0.9 m) probably due to the undulating topography surrounding the location. The magnitude of C _D is high for low wind speed (<1.5 m s^?1) and found to be in the range 0.005–0.03. The variations of sensible heat fluxes (SHF) and momentum fluxes are also discussed. Relatively high fluxes of heat and momentum are observed during typical days on 26–27 February 2004 and 10–11 April 2004 due to the daytime unstable atmospheric conditions. Stable or near neutral conditions prevail after 1700 h IST with negative SHF. A mesoscale model PSU/NCAR MM5 is run using a high-resolution (1 km) grid over the study region to examine the influence of complex topography on the surface layer parameters and the simulated fluxes are compared with estimated values. Spatial variations of the frictional velocity (u _*), C _D, surface fluxes, planetary boundary layer (PBL) height and surface winds are noticed according to the topographic variations in the simulation.     

Turbulent exchange above a pine forest,I: Fluxes and gradients

Measurements of gradients of wind, temperature and humidity and of the corresponding turbulent fluxes have been carried out over a sparse pine forest at Jädra»s in Sweden. In order to ascertain that correct gradient estimates were obtained, two independent measuring systems were employed: one system with sensors at 10 fixed levels on a 51 m tower and another with reversing sensors for temperature and humidity, covering the height interval 23 to 32 m. Turbulent fluxes were measured at three levels simultaneously. Data from three field campaigns: in June 1985, June 1987 and September 1987 have been analyzed. The momentum flux is found on the average to be virtually constant from tree top level, at 20 to 50 m. The average fluxes of sensible and latent heat are not so well behaved. The ratio of the non-dimensional gradients of wind and temperature to their corresponding values under ideal conditions (low vegetation) are both found to be small immediately above the canopy (about 0.3 for temperature and 0.4 for wind). With increasing height, the ratios increase, but the values vary substantially with wind direction. The ratios are not found to vary systematically with stability (unstable stratification only studied). The ratio of the non-dimensional humidity gradient to the corresponding non-dimensional potential temperature gradient (equivalent to k _h /k _w ) is found to be unity for (z – d)/L _v less than about –0.1 and about 1.4 for near neutral stratification, but the scatter of the data is very large.     

Multiple temporal scale variability during the twentieth century in global carbon dynamics simulated by a coupled climate–terrestrial carbon cycle model

A coupled climate–carbon cycle model composed of a process-based terrestrial carbon cycle model, Sim-CYCLE, and the CCSR/NIES/FRCGC atmospheric general circulation model was developed. We examined the multiple temporal scale functions of terrestrial ecosystem carbon dynamics induced by human activities and natural processes and evaluated their contribution to fluctuations in the global carbon budget during the twentieth century. Global annual net primary production (NPP) and heterotrophic respiration (HR) increased gradually by 6.7 and 4.7%, respectively, from the 1900s to the 1990s. The difference between NPP and HR was the net carbon uptake by natural ecosystems, which was 0.6 Pg C year^?1 in the 1980s, whereas the carbon emission induced by human land-use changes was 0.5 Pg C year^?1, largely offsetting the natural terrestrial carbon sequestration. Our results indicate that monthly to interannual variation in atmospheric CO₂ growth rate anomalies show 2- and 6-month time lags behind anomalies in temperature and the NiNO₃ index, respectively. The simulated anomaly amplitude in monthly net carbon flux from terrestrial ecosystems to the atmosphere was much larger than in the prescribed air-to-sea carbon flux. Fluctuations in the global atmospheric CO₂ time series were dominated by the activity of terrestrial vegetation. These results suggest that terrestrial ecosystems have acted as a net neutral reservoir for atmospheric CO₂ concentrations during the twentieth century on an interdecadal timescale, but as the dominant driver for atmospheric CO₂ fluctuations on a monthly to interannual timescale.