Peter S. Bakwin

A global terrestrial monitoring network integrating tower fluxes, flask sampling, ecosystem modeling and EOS satellite data

Abstract Accurate monitoring of global scale changes in the terrestrial biosphere has become acutely important as the scope of human impacts on biological systems and atmospheric chemistry grows. For example, the Kyoto Protocol of 1997 signals some of the dramatic socioeconomic and political decisions that may lie ahead concerning CO2 emissions and global carbon cycle impacts. These decisions will rely heavily on accurate measures of global biospheric changes Schimel 1998 , IGBP TCWG 1998 . An array of national and international programs have inaugurated global satellite observations, critical field measurements of carbon and water fluxes, and global model development for the purposes of beginning to monitor the biosphere. The detection by these programs of interannual variability of ecosystem fluxes and of longer term trends will permit early indication of fundamental biospheric changes which might otherwise go undetected until major biome conversion begins. This article describes a blueprint for more comprehensive coordination of the various flux measurement and modeling activities into a global terrestrial monitoring network that will have direct relevance to the political decision making of global change.

Phase and amplitude of ecosystem carbon release and uptake potentials as derived from FLUXNET measurements

As length and timing of the growing season are major factors explaining differences in carbon exchange of ecosystems, we analyzed seasonal patterns of net ecosystem carbon exchange (FNEE) using eddy covariance data of the FLUXNET data base (http://www-eosdis.ornl.gov/FLUXNET). The study included boreal and temperate, deciduous and coniferous forests, Mediterranean evergreen systems, rainforest, native and managed temperate grasslands, tundra, and C3 and C4 crops. Generalization of seasonal patterns are useful for identifying functional vegetation types for global dynamic vegetation models, as well as for global inversion studies, and can help improve phenological modules in SVAT or biogeochemical models. The results of this study have important validation potential for global carbon cycle modeling. The phasing of respiratory and assimilatory capacity differed within forest types: for temperate coniferous forests seasonal uptake and release capacities are in phase, for temperate deciduous and boreal coniferous forests, release was delayed compared to uptake. According to seasonal pattern of maximum nighttime release (evaluated over 15-day periods, Fmax) the study sites can be grouped in four classes: (1) boreal and high altitude conifers and grasslands; (2) temperate deciduous and temperate conifers; (3) tundra and crops; (4) evergreen Mediterranean and tropical forests. Similar results are found for maximum daytime uptake (Fmin) and the integral net carbon flux, but temperate deciduous forests fall into class 1. For forests, seasonal amplitudes of Fmax and Fmin increased in the order tropical C3-crops>temperate deciduous forests>temperate conifers>boreal conifers>tundra ecosystems. Due to data restrictions, our analysis centered mainly on Northern Hemisphere temperate and boreal forest ecosystems. Grasslands, crops, Mediterranean ecosystems, and rainforests are under-represented, as are savanna systems, wooded grassland, shrubland, or year-round measurements in tundra systems. For regional or global estimates of carbon sequestration potentials, future investigations of eddy covariance should expand in these systems.

Long-Term Carbon Dioxide Fluxes from a Very Tall Tower in a Northern Forest: Flux Measurement Methodology

Methodology for determining fluxes of CO 2 and H2O vapor with the eddy-covariance method using data from instruments on a 447-m tower in the forest of northern Wisconsin is addressed. The primary goal of this study is the validation of the methods used to determine the net ecosystem exchange of CO 2. Two-day least squares fits coupled with 30-day running averages limit calibration error of infrared gas analyzers for CO2 and H2O signals to 2%‐3%. Sonic anemometers are aligned with local streamlines by fitting a sine function to tilt and wind direction averages, and fitting a third-order polynomial to the residual. Lag times are determined by selecting the peak in lagged covariance with an error of 1.5%‐2% for CO2 and 1% for H2O vapor. Theory and a spectral fit method allow determination of the underestimation in CO2 flux ( ,5% daytime, ,12% nighttime) and H2O vapor flux ( ,21%), which is due to spectral degradation induced by long air-sampling tubes. Scale analysis finds 0.5-h flux averaging periods are sufficient to measure all flux scales at 30-m height, but 1 h is necessary at higher levels, and random errors in the flux measurements due to limited sampling of atmospheric turbulence are fairly large ( 15%‐20% for CO2 and 20%‐40% for H2O vapor at lower levels for a 1-h period).

Summertime photochemistry of the troposphere at high northern latitudes

The budgets of O3, NOx (NO+NO2), reactive nitrogen (NOy), and acetic acid in the 0–6 km column over western Alaska in summer are examined by photochemical modeling of aircraft and ground-based measurements from the Arctic Boundary Layer Expedition (ABLE 3A). It is found that concentrations of O3 in the region are regulated mainly by input from the stratosphere, and losses of comparable magnitude from photochemistry and deposition. The concentrations of NOx (10–50 ppt) are sufficiently high to slow down O3 photochemical loss appreciably relative to a NOx-free atmosphere; if no NOx were present, the lifetime of O3 in the 0–6 km column would decrease from 46 to 26 days because of faster photochemical loss. The small amounts of NOx present in the Arctic troposphere have thus a major impact on the regional O3 budget. Decomposition of peroxyacetyl nitrate (PAN) can account for most of the NOx below 4-km altitude, but for only 20% at 6-km altitude. Decomposition of other organic nitrates might supply the missing source of NOx. The lifetime of NOy, in the ABLE 3A flight region is estimated at 29 days, implying that organic nitrate precursors of NOx could be supplied from distant sources including fossil fuel combustion at northern mid-latitudes. Biomass fire plumes sampled during ABLE 3A were only marginally enriched in O3; this observation is attributed in part to low NOx emissions in the fires, and in part to rapid conversion of NOx to PAN promoted by low atmospheric temperatures. It appears that fires make little contribution to the regional O3 budget. Only 30% of the acetic acid concentrations measured during ABLE 3A can be accounted for by reactions of CH3CO3 with HO2 and CH3O2. There remains a major unidentified source of acetic acid in the atmosphere.

Climatologies of NOx and NOy: A comparison of data and models

Abstract Climatologies of tropospheric NOx (NO + NO2) and NOy (total reactive nitrogen: NOx + N03 + 2 × N2O5 + HNO2 + HNO3 + HNO4 + ClONO2 + PAN (peroxyacetylnitrate) + other organic ni trates) have been compiled from data previously published and, in most cases, publicly archived. Emphasis has been on non-urban measurements, including rural and remote ground sites, as well as aircraft data. Although the distribution of data is sparse, a compilation in this manner can begin to provide an understanding of the spatial and temporal distributions of these reactive nitrogen species. The cleanest measurements in the boundary layer are in Alaska, northern Canada and the eastern Pacific, with median NO mixing ratios below 10 pptv, NOx below 50 pptv, and NOy below 300 pptv. The highest NO values (greater than 1 ppbv) were found in eastern North America and Europe, with correspondingly high NOy (∼ 5 ppbv). A significantly narrower range of concentrations is seen in the free troposphere, particularly at 3–6 km, with NO typically about 10 pptv in the boreal summer. NO increases with altitude to ∼ 100 pptv at 9–12 km, whereas NOy does not show a trend with altitude, but varies between 100 and 1000 pptv. Decreasing mixing ratios eastward of the Asian and North American continents are seen in all three species at all altitudes. Model-generated climatologies of NOx and NOy from six chemical transport models are also presented and are compared with observations in the boundary layer and the middle troposphere for summer and winter. These comparisons test our understanding of the chemical and transport processes responsible for these species distributions. Although the model results show differences between them, and disagreement with observations, none are systematically different for all seasons and altitudes. Some of the differences between the observations and model results may likely be attributed to the specific meteorological conditions at the time that measurements were made differing from the model meteorology, which is either climatological flow from GCMs or actual meteorology for an arbitrary year. Differences in emission inventories, and convection and washout schemes in the models will also affect the calculated NOα and NOy distributions.

Influence of advection on measurements of the net ecosystem‐atmosphere exchange of CO2 from a very tall tower

In most studies of the net ecosystem-atmosphere exchange of CO2 (NEE) using tower-based eddy covariance (EC) systems it has been assumed that advection is negligible. In this study we use a scalar conservation budget method to estimate the contribution of advection to NEE measurements from a very tall tower in northern Wisconsin. We examine data for June-August 1997. Measured NEE0, calculated as the sum of the EC flux plus the rate of change of storage below the EC measurement level, is expected to be constant with measurement height, and we take the differences between levels as a measure of advection. We find that the average difference in total advection DFCadtot between 30 and 122 m is as large as 6 mmol m 22 s 21 during the morning transition from stable to convective conditions and the average difference DFCadtot between 122 and 396 m is as large as 4 mmol m 22 s 21 during daytime. For the month of July, advection between 30 and 122 m is 27% of the diurnally integrated NEE0 at 122 m, and advection between 122 and 396 m accounts for 5% of the NEE0 observed at 396 m. The observed differences of advection often have significant correlation with the vertical integral of wind speed within the same layer. This indicates that the horizontal advection contribution to NEE could be significant. Direct observations of the vertical gradient in CO2 show that DFCadtot cannot be explained by vertical advection alone. It is hypothesized that differing flux footprints and pooling of CO2 in the heterogeneous landscape causes the advection contribution. The magnitudes of the total advection component FCadtot of NEE at the 30 m level are roughly estimated by a linear extrapolation. A peak in F Cadtot at 30 m of ;3 mmol m 22 s 21 during the morning transition is predicted for all three months. The July integrated F Cadtot is estimated to be 10% of the diurnally integrated NEE0 at 30 m.

Determination of the isotopic(13C/12C) discrimination by terrestrial biology from a global network of observations

We analyze data from the National Oceanic and Atmospheric Administration / Climate Monitoring and Diagnostics Laboratory global air sampling network in order to extract the signatures of isotopic (13C/12C) discrimination by the terrestrial biota and of fossil fuel combustion for the regions surrounding the sampling sites. We utilize measurements of carbon monoxide (CO) to give an estimate of the contribution of fossil fuel combustion to the short-term variability of carbon dioxide (CO2). In general, variations of CO2 are more strongly dominated by biological exchange, so the isotopic signature of fossil fuel combustion, while consistent with inventory estimates, is not well constrained by the observations. Conversely, results for isotope discrimination by the terrestrial biosphere are not strongly dependent on assumptions about fossil fuel combustion. Our analysis appears valid primarily for stations fairly near continental source/sink regions, particularly for midlatitude regions of the northern hemisphere. For these stations we derive a mean discrimination of −16.8 per mil (‰), with site-to-site variability of 0.8‰ (1 standard deviation) and with little or no consistent latitudinal gradient.

What is the concentration footprint of a tall tower

Studies that have attempted to estimate sources and sinks of trace gases such as CO2 with inverse calculations unanimously identify the lack of continental stations as a prime obstacle. Continental stations have traditionally been avoided because of the difficulty of interpretation due to large time-variability of trace substance mixing ratios. Large variability is caused by the proximity to the strongly variable sources in space and time and the complicated airflow within the lowermost 100–200 m of the planetary boundary layer. To address the need for continental stations and to overcome the problems associated with them, the National Oceanic and Atmospheric Administration Climate Monitoring and Diagnostics Laboratory started in 1992 to measure CO2 and other trace gases on tall television transmission towers [Bakwin et al., 1995]. An essential question in connection with these tower measurements is the area around the tower from which fluxes substantially contribute to the observed short-term variability of trace gas mixing ratios. We present here a simple data and back trajectory-based method to estimate the fraction of the observed short-term variability explained by a localized flux around a tall television transmission tower in Wisconsin in dependence of its location relative to the tower (the concentration “footprint”). We find that the timescale over which the imprint of surface fluxes on air parcels before its arrival at the tower are still discernible in the mixing ratio variations observed at the tower is of the order of 1.5 days. Based on this timescale and the characteristics of air parcel trajectories, we infer a spatial extent of the footprint of the order of 106 km2, or roughly a tenth of the area of the United States. This is encouraging evidence that tall tower measurements may be useful in global inversions and may also have implications for monitoring fluxes of anthropogenic trace substances on regional scales.

Regional budgets for nitrogen oxides from continental sources: Variations of rates for oxidation and deposition with season and distance from source regions

Measurements of nitrogen deposition and concentrations of NO, NO2, NOy (total oxidized N), and O3 have been made at Harvard Forest in central Massachusetts since 1990 to define the atmospheric budget for reactive N near a major source region. Total (wet plus dry) reactive N deposition for the period 1990–1996 averaged 47 mmol m−2 yr−1 (126 μmol m−2 d−1, 6.4 kg N ha−1 yr−1), with 34% contributed by dry deposition. Atmospheric input adds about 12% to the N made available annually by mineralization in the forest soil. The corresponding deposition rate at a distant site, Schefferville, Quebec, was 20 mmol m−2 d−1 during summer 1990. Both heterogeneous and homogeneous reactions efficiently convert NOx to HNO3 in the boundary layer. HNO3 is subsequently removed rapidly by either dry deposition or precipitation. The characteristic (e-folding) time for NOx oxidation ranges from 0.30 days in summer, when OH radical is abundant, to ∼1.5 days in the winter, when heterogeneous reactions are dominant and O3 concentrations are lowest. The characteristic time for removal of NOx oxidation products (defined as NOy minus NOx) from the boundary layer by wet and dry deposition is ∼1 day, except in winter when it decreases to 0.6 day. Biogenic hydrocarbons contribute to N deposition through formation of organic nitrates but are also precursors of reservoir species, such as peroxyacetylnitrate, that may be exported from the region. A simple model assuming pseudo first-order rates for oxidation of NOx, followed by deposition, predicts that 45% of NOx in the northeastern U.S. boundary layer is removed in 1 day during summer and 27% is removed in winter. It takes 3.5 and 5 days for 95% removal in summer and winter, respectively.

Long-Term Observations of the Dynamics of the Continental Planetary Boundary Layer

Time series of mixed layer depth, zi, and stable boundary layer height from March through October of 1998 are derived from a 915-MHz boundary layer profiling radar and CO 2 mixing ratio measured from a 447-m tower in northern Wisconsin. Mixed layer depths from the profiler are in good agreement with radiosonde measurements. Maximum zi occurs in May, coincident with the maximum daytime surface sensible heat flux. Incoming radiation is higher in June and July, but a greater proportion is converted to latent heat by photosynthesizing vegetation. An empirical relationship between zi and the square root of the cumulative surface virtual potential temperature flux is obtained ( r 2 5 0.98) allowing estimates of zi from measurements of virtual potential temperature flux under certain conditions. In fair-weather conditions the residual mixed layer top was observed by the profiler on several nights each month. The synoptic mean vertical velocity (subsidence rate) is estimated from the temporal evolution of the residual mixed layer height during the night. The influence of subsidence on the evolution of the mixed, stable, and residual layers is discussed. The CO 2 jump across the inversion at night is also estimated from the tower measurements.