Estimates of anthropogenic CO<Subscript>2</Subscript> uptake in a global ocean model

A global ocean general circulation model (L30T63) is employed to study the uptake and distribution of anthropogenic CO₂ in the ocean. A subgrid-scale mixing scheme called GM90 is used in the model. There are two main GM90 parameters including isopycnal diffusivity and skew (thickness) diffusivity. Sensitivities of the ocean circulation and the redistribution of dissolved anthropogenic CO₂ to these two parameters are examined. Two runs estimate the global oceanic anthropogenic CO₂ uptake to be 1.64 and 1.73 Pg C yr^-1 for the 1990s, and that the global ocean contained 86.8 and 92.7 Pg C of anthropogenic CO₂ at the end of 1994, respectively. Both the total inventory and uptake from our model are smaller than the data-based estimates. In this presentation, the vertical distributions of anthropogenic CO₂ at three meridional sections are discussed and compared with the available data-based estimates. The inventory in the individual basins is also calculated. Use of large isopycnal diffusivity can generally improve the simulated results, including the exchange flux, the vertical distribution patterns, inventory, storage, etc. In terms of comparison of the vertical distributions and column inventory, we find that the total inventory in the Pacific Ocean obtained from our model is in good agreement with the data-based estimate, but a large difference exists in the Atlantic Ocean, particularly in the South Atlantic. The main reasons are weak vertical mixing and that our model generates small exchange fluxes of anthropogenic CO₂ in the Southern Ocean. Improvement in the simulation of the vertical transport and sea ice in the Southern Ocean is important in future work.     

Long-term effects of anthropogenic CO<Subscript>2</Subscript> emissions simulated with a complex earth system model

A new complex earth system model consisting of an atmospheric general circulation model, an ocean general circulation model, a three-dimensional ice sheet model, a marine biogeochemistry model, and a dynamic vegetation model was used to study the long-term response to anthropogenic carbon emissions. The prescribed emissions follow estimates of past emissions for the period 1751–2000 and standard IPCC emission scenarios up to the year 2100. After 2100, an exponential decrease of the emissions was assumed. For each of the scenarios, a small ensemble of simulations was carried out. The North Atlantic overturning collapsed in the high emission scenario (A2) simulations. In the low emission scenario (B1), only a temporary weakening of the deep water formation in the North Atlantic is predicted. The moderate emission scenario (A1B) brings the system close to its bifurcation point, with three out of five runs leading to a collapsed North Atlantic overturning circulation. The atmospheric moisture transport predominantly contributes to the collapse of the deep water formation. In the simulations with collapsed deep water formation in the North Atlantic a substantial cooling over parts of the North Atlantic is simulated. Anthropogenic climate change substantially reduces the ability of land and ocean to sequester anthropogenic carbon. The simulated effect of a collapse of the deep water formation in the North Atlantic on the atmospheric CO₂ concentration turned out to be relatively small. The volume of the Greenland ice sheet is reduced, but its contribution to global mean sea level is almost counterbalanced by the growth of the Antarctic ice sheet due to enhanced snowfall. The modifications of the high latitude freshwater input due to the simulated changes in mass balance of the ice sheet are one order of magnitude smaller than the changes due to atmospheric moisture transport. After the year 3000, the global mean surface temperature is predicted to be almost constant due to the compensating effects of decreasing atmospheric CO₂ concentrations due to oceanic uptake and delayed response to increasing atmospheric CO₂ concentrations before.     

Consequences of Uncertainties in CO<Subscript>2</Subscript> Density for Estimating Net Ecosystem CO<Subscript>2</Subscript> Exchange by Open-path Eddy Covariance

Errors in the estimation of CO₂ surface exchange by open-path eddy covariance, introduced during the removal of density terms [Webb et al. Quart J Roy Meteorol Soc 106:85–100, (1980) - WPL], can happen both because of errors in energy fluxes [Liu et al. Boundary-Layer Meteorol 120:65–85, (2006)] but also because of inaccuracies in other terms included in the density corrections, most notably due to measurements of absolute CO₂ density (ρ _c ). Equations are derived to examine the propagation of all errors through the WPL algorithm. For an open-path eddy covariance system operating in the Sierra de Gádor in south-east Spain, examples are presented of the inability of an unattended, open-path infrared gas analyzer (IRGA) to reliably report ρ _c and the need for additional instrumentation to determine calibration corrections. A sensitivity analysis shows that relatively large and systematic errors in net ecosystem exchange (NEE) can result from uncertainties in ρ _c in a semi-arid climate with large sensible heat fluxes (H _s ) and (wet) mineral deposition. When ρ_c is underestimated by 5% due to lens contamination, this implies a 13% overestimation of monthly CO₂ uptake.     

Choosing a Stabilization Target for CO<Subscript>2</Subscript>

The choice of stabilization target for CO₂ concentration depends on the following: what is considered to be dangerous anthropogenic interference with the climate system; the forcings that might arise from non-CO₂ gases; and the climate sensitivity. These three factors are specified here probabilistically, as probability density functions (pdfs), and combined to produce a pdf for the CO₂ concentration target. There is a probability of 17% that the stabilization target should be less than the present level, and the median target is 536 ppm. The effects of reducing the emissions of non-CO₂ gases and/or implementing adaptation strategies are considered probabilistically and shown to alter these figures significantly.     

Characterizing the annual-mean climatic effect of anthropogenic CO<Subscript>2</Subscript> and aerosol emissions in eight coupled atmosphere-ocean GCMs

This work examines the spatial patterns of the transient response of mean annual temperature and precipitation to CO₂ (or CO₂ plus aerosol or aerosol proxy) radiative forcing in eight coupled AOGCMs, generally for the period 1900–2099. Response patterns are characterized using empirical orthogonal functions (EOFs) and the quasi-EOFs of Harvey and Wigley (the first qEOF field, discussed here, is given by the correlation between local year-by-year temperature changes and the global mean temperature change). The first temperature EOF accounts for 80–95% of the space-time variation of the CO₂ run in all of the models, and is almost identical to qEOF₁ of the temperature response or to the temperature change pattern averaged over the last 30 years of the simulations. EOF₁ accounts for 80–95% of the space-time variation in the CO₂+aerosol runs in six of the eight models. The CO₂ response patterns of different models are highly correlated with one another (R ² generally >0.5), and are also highly correlated with the CO₂+aerosol response patterns (R ² 0.85 in all except one model). The difference between CO₂ and CO₂+aerosol runs can be represented by EOF₁ of the year-by-year differences, by qEOF₁ of the year-by-year differences, or by the difference in temperature averaged over the last 30 years of each run. In models where these representations are highly correlated with each other, they are also highly correlated with CO₂ EOF₁. In other cases, aerosol EOF₁ is modestly to highly correlated with control EOF₁ (i.e.: the year-by-year differences between CO₂ and CO₂+aerosol runs are dominated by internal variability), while aerosol qEOF₁ and the 30-year difference are highly correlated with each other. For all models, the decadal mean temperature change can be closely replicated by scaling the CO₂ EOF₁ pattern based on the global mean temperature changes (RMSE for the last decade is <6% of the RMS temperature change for CO₂ runs, <8% for CO₂+aerosol runs). The first EOF of the precipitation response to increasing CO₂ accounts for only 10–30% of the space-time variation, and is generally highly correlated (R ² up to 0.85) with control EOF₁. In all of the models, there is an increase in precipitation in the ITCZ and a decrease in bands at or near 30°S and 30°N. In many models there is an El Niño-like response, including a substantial decrease in precipitation over the Amazon. Global-mean precipitation increases in all models due to CO₂ forcing, but aerosols appear to have a disproportionally large effect in suppressing the increase compared to their effect in suppressing the warming. There is evidence in some models that the non-absorbing aerosols considered here reduce summer monsoon rainfall compared to the changes that would be expected based on the globally averaged effect of aerosols on precipitation. When regional precipitation changes over time are predicted by scaling a fixed precipitation-change pattern with the global mean temperature change, the global mean RMSE in the predicted change in decadal-mean precipitation is 25–35% of the global RMS precipitation changes by the end of the simulation.     

Quick measurement of CH<Subscript>4</Subscript>, CO<Subscript>2</Subscript> and N<Subscript>2</Subscript>O emissions from a short-plant ecosystem

Combining improved injector, gas line and valve-driving models, a gas chromatograph (GC) equipped with Hydrogen Flame Ionization Detector (FID) and Electron Capture Detector (ECD), can measure CH4, CO2, and N2O simultaneously in an air sample in four minutes. Test results show that the system has high sensitivity, resolution, and precision; the linear response range of the system meets the requirement of flux measurements in situ. The system is suitable for monitoring fluxes of the main greenhouse gases in a short-plant field since it is easy to use, efficacious, and constant and reliable in collecting data.     

CO<Subscript>2</Subscript> and non-CO<Subscript>2</Subscript> radiative forcings in climate projections for twenty-first century mitigation scenarios

Climate is simulated for reference and mitigation emissions scenarios from Integrated Assessment Models using the Bern2.5CC carbon cycle–climate model. Mitigation options encompass all major radiative forcing agents. Temperature change is attributed to forcings using an impulse–response substitute of Bern2.5CC. The contribution of CO₂ to global warming increases over the century in all scenarios. Non-CO₂ mitigation measures add to the abatement of global warming. The share of mitigation carried by CO₂, however, increases when radiative forcing targets are lowered, and increases after 2000 in all mitigation scenarios. Thus, non-CO₂ mitigation is limited and net CO₂ emissions must eventually subside. Mitigation rapidly reduces the sulfate aerosol loading and associated cooling, partly masking Greenhouse Gas mitigation over the coming decades. A profound effect of mitigation on CO₂ concentration, radiative forcing, temperatures and the rate of climate change emerges in the second half of the century.     

A nocturnal atmospheric loss of CH<Subscript>2</Subscript>I<Subscript>2</Subscript> in the remote marine boundary layer

Ocean emissions of inorganic and organic iodine compounds drive the biogeochemical cycle of iodine and produce reactive ozone-destroying iodine radicals that influence the oxidizing capacity of the atmosphere. Di-iodomethane (CH₂I₂) and chloro-iodomethane (CH₂ICl) are the two most important organic iodine precursors in the marine boundary layer. Ship-borne measurements made during the TORERO (Tropical Ocean tRoposphere Exchange of Reactive halogens and Oxygenated VOC) field campaign in the east tropical Pacific Ocean in January/February 2012 revealed strong diurnal cycles of CH₂I₂ and CH₂ICl in air and of CH₂I₂ in seawater. Both compounds are known to undergo rapid photolysis during the day, but models assume no night-time atmospheric losses. Surprisingly, the diurnal cycle of CH₂I₂ was lower in amplitude than that of CH₂ICl, despite its faster photolysis rate. We speculate that night-time loss of CH₂I₂ occurs due to reaction with NO₃ radicals. Indirect results from a laboratory study under ambient atmospheric boundary layer conditions indicate a k _CH2I2+NO3 of ≤4 × 10^?13 cm³ molecule^?1 s^?1; a previous kinetic study carried out at ≤100 Torr found k _CH2I2+NO3 of 4 × 10^?13 cm³ molecule^?1 s^?1. Using the 1-dimensional atmospheric THAMO model driven by sea-air fluxes calculated from the seawater and air measurements (averaging 1.8 +/? 0.8 nmol m^?2 d^?1 for CH₂I₂ and 3.7 +/? 0.8 nmol m^?2 d^?1 for CH₂ICl), we show that the model overestimates night-time CH₂I₂ by >60 % but reaches good agreement with the measurements when the CH₂I₂ + NO₃ reaction is included at 2–4 × 10^?13 cm³ molecule^?1 s^?1. We conclude that the reaction has a significant effect on CH₂I₂ and helps reconcile observed and modeled concentrations. We recommend further direct measurements of this reaction under atmospheric conditions, including of product branching ratios.     

CO<Subscript>2</Subscript> dynamics on the shelf of the East Siberian Sea

Expedition data obtained in the coastal-shelf zone of the East Siberian Sea in September 2003, 2004, and 2008 are generalized. Studies of carbonate system in water and CO₂ fluxes between ocean and atmosphere in this region confirmed that it was reasonable to divide the water area studied into two biogeochemical provinces and that the ecosystem of its coastal part is mainly of a heterotrophic nature. In different years, the extent of water supersaturation in carbon dioxide in the East Siberian Sea and the area of the CO₂ release significantly changed. Geographic localization of the atmosphere action centers over the Arctic and their intensity were main determining factors; that told both on the formation of a basic character of the atmospheric and hydrological processes and on the dynamics of the CO₂ exchange between water and air.     

Effect of Data Assimilation Parameters on The Optimized Surface CO<Subscript>2</Subscript> Flux in Asia

In this study, CarbonTracker, an inverse modeling system based on the ensemble Kalman filter, was used to evaluate the effects of data assimilation parameters (assimilation window length and ensemble size) on the estimation of surface CO₂ fluxes in Asia. Several experiments with different parameters were conducted, and the results were verified using CO₂ concentration observations. The assimilation window lengths tested were 3, 5, 7, and 10 weeks, and the ensemble sizes were 100, 150, and 300. Therefore, a total of 12 experiments using combinations of these parameters were conducted. The experimental period was from January 2006 to December 2009. Differences between the optimized surface CO₂ fluxes of the experiments were largest in the Eurasian Boreal (EB) area, followed by Eurasian Temperate (ET) and Tropical Asia (TA), and were larger in boreal summer than in boreal winter. The effect of ensemble size on the optimized biosphere flux is larger than the effect of the assimilation window length in Asia, but the importance of them varies in specific regions in Asia. The optimized biosphere flux was more sensitive to the assimilation window length in EB, whereas it was sensitive to the ensemble size as well as the assimilation window length in ET. The larger the ensemble size and the shorter the assimilation window length, the larger the uncertainty (i.e., spread of ensemble) of optimized surface CO₂ fluxes. The 10-week assimilation window and 300 ensemble size were the optimal configuration for CarbonTracker in the Asian region based on several verifications using CO₂ concentration measurements.     

Interpreting CO<Subscript>2</Subscript> Fluxes Over a Suburban Lawn: The Influence of Traffic Emissions

Turf-grass lawns are ubiquitous in the United States. However direct measurements of land–atmosphere fluxes using the eddy-covariance method above lawn ecosystems are challenging due to the typically small dimensions of lawns and the heterogeneity of land use in an urbanised landscape. Given their typically small patch sizes, there is the potential that CO₂ fluxes measured above turf-grass lawns may be influenced by nearby CO₂ sources such as passing traffic. In this study, we report on two years of eddy-covariance flux measurements above a 1.5 ha turf-grass lawn in which we assess the contribution of nearby traffic emissions to the measured CO₂ flux. We use winter data when the vegetation was dormant to develop an empirical estimate of the traffic effect on the measured CO₂ fluxes, based on a parametrised version of a three-dimensional Lagrangian footprint model and continuous traffic count data. The CO₂ budget of the ecosystem was adjusted by 135gCm⁻² in 2007 and by 134gCm⁻² in 2008 to determine the natural flux, even though the road crossed the footprint only at its far edge. We show that bottom-up flux estimates based on CO₂ emission factors of the passing vehicles, combined with the crosswind-integrated footprint at the distance of the road, agreed very well with the empirical estimate of the traffic contribution that we derived from the eddy-covariance measurements. The approach we developed may be useful for other sites where investigators plan to make eddy-covariance measurements on small patches within heterogeneous landscapes where there are significant contrasts in flux rates. However, we caution that the modelling approach is empirical and will need to be adapted individually to each site.     

Influence of dynamic vegetation on climate change arising from increasing CO<Subscript>2</Subscript>

A dynamic global vegetation model (DGVM) is coupled to an atmospheric general circulation model (AGCM) to investigate the influence of vegetation dynamics on climate change under conditions of global warming. The model results are largely in agreement with observations and the results of previous studies in terms of the present climate, present potential vegetation, present net primary productivity (NPP), and pre-industrial carbon budgets. The equilibrium state of climate properties are compared among pre-industrial, doubled, and quadrupled atmospheric CO₂ values using DGVM–AGCM and current AGCM with fixed vegetation to evaluate the influence of dynamic vegetation change. We also separated the contributions of temperature, precipitation and CO₂ fertilization on vegetation change. The results reveal an amplification of global warming climate sensitivity by 10% due to the inclusion of dynamic vegetation. The total effects of elevated CO₂ and climate change also lead to an increase in NPP and vegetation coverage globally. The reduction of albedo associated with this greening results in enhanced global warming. Our separation analysis indicates that temperature alters vegetation at high latitudes such as Siberia or Alaska, where there is a switch from tundra to forest. On the other hand, CO₂ fertilization provides the largest contribution to greening in arid/semi-arid region. Precipitation change did not cause any drastic vegetation shift.     

Tilt corrections over complex terrain and their implication for CO<Subscript>2</Subscript> transport

Observational data from sonic anemometers are commonly rotated from sonic to streamline coordinates, a procedure that is called tilt correction. Tilt corrections are often used to post-process air velocity data collected from sonic anemometers to allow objective interpretation of air flow data relative to the Earth. Since streamline coordinates depend on dynamical characteristics of the flow, the tilt correction depends not only on temporal and spatial variations of the flow, but also on local circulations. We found that ensemble- averaged slope flows are approximately parallel to the terrain slope close to the ground within the canopy layer, but not above, due to the influence of the diurnal variation of local vertical circulations. As a result, the diurnal variation of the observed vertical velocity in streamline coordinates at 21.5 m above the ground over 11-m tall forest canopies can be opposite to that calculated from the continuity equation. To estimate CO₂ transport over sloping terrain, a workable reference coordinate system is needed such that multiple sonic anemometers have a common reference relative to the Earth. Streamline coordinate systems can be the choice of the common reference coordinate system only if flow, at least ensemble-averaged flow, is parallel to terrain slopes. The choice of the reference coordinate system and its implication in investigation of CO₂ transport are discussed. The National Center for Atmospheric Research is sponsored by the National Science Foundation.     

Seasonal and diurnal dynamics of CO<Subscript>2</Subscript> emission from oligotrophic peat soil surface

The results of research of diurnal and seasonal dynamics of CO₂ emission from the oligotrophic swamp surface in the southern taiga subzone of Western Siberia in 2005–2007 are under consideration. During the summertime, the intensity of CO₂ emission increases from spring to the midsummer and then decreases by the fall. A mean CO₂ emission value was 118 mg CO₂/(m² hour). The analysis of diurnal dynamics of CO₂ emission showed that the maximum CO₂ flux is observed at 16:00, while the minimum, at 07:00. Mean amplitude of diurnal variations of the CO₂ emission is 74 mg CO₂/(m² hour). The relations established between air temperature and CO₂ flux allowed calculating carbon dioxide emission for the periods between measurements. It was found that in the summertime, the period between 10:00 and 13:00 was optimal for measuring CO₂ emission with a chamber method.     

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.     

Simulation of Non-Homogeneous CO<Subscript>2</Subscript> and Its Impact on Regional Temperature in East Asia

Carbon dioxide (CO₂) is an important greenhouse gas that influences regional climate through disturbing the earth’s energy balance. The CO₂ concentrations are usually prescribed homogenously in most climate models and the spatiotemporal variations of CO₂ are neglected. To address this issue, a regional climate model (RegCM4) is modified to investigate the non-homogeneous distribution of CO₂ and its effects on regional longwave radiation flux and temperature in East Asia. One-year simulation is performed with prescribed surface CO₂ fluxes that include fossil fuel emission, biomass burning, air–sea exchange, and terrestrial biosphere flux. Two numerical experiments (one using constant prescribed CO₂ concentrations in the radiation scheme and the other using the simulated CO₂ concentrations that are spatially non-homogeneous) are conducted to assess the impact of non-homogeneous CO₂ on the regional longwave radiation flux and temperature. Comparison of CO₂ concentrations from the model with the observations from the GLOBALVIEW-CO₂ network suggests that the model can well capture the spatiotemporal patterns of CO₂ concentrations. Generally, high CO₂ mixing ratios appear in the heavily industrialized eastern China in cold seasons, which probably relates to intensive human activities. The accommodation of non-homogeneous CO₂ concentrations in the radiative transfer scheme leads to an annual mean change of–0.12 W m^–2 in total sky surface upward longwave flux in East Asia. The experiment with non-homogeneous CO₂ tends to yield a warmer lower troposphere. Surface temperature exhibits a maximum difference in summertime, ranging from–4.18 K to 3.88 K, when compared to its homogeneous counterpart. Our results indicate that the spatial and temporal distributions of CO₂ have a considerable impact on regional longwave radiation flux and temperature, and should be taken into account in future climate modeling.     

Validation of Column-Averaged Dry-Air Mole Fraction of CO<Subscript>2</Subscript> Retrieved from OCO-2 Using Ground-Based FTS Measurements

In order to correctly use the column-averaged atmospheric CO₂ dry-air mole fraction (XCO₂) data in the CO₂ flux studies, XCO₂ measurements retrieved from the Orbiting Carbon Observatory-2 (OCO-2) in 2015 were compared with those obtained from the global ground-based high-resolution Fourier Transform Spectrometer (FTS) participating in the Total Carbon Column Observing Network (TCCON). The XCO₂ retrieved from three observing modes adopted by OCO-2, i.e., nadir, target, and glint, were separately validated by the FTS measurements at up to eight TCCON stations located in different areas. These comparisons show that OCO-2 glint mode yields the best qualitative estimations of CO₂ concentration among the three operational approaches. The overall results regarding the glint mode show no obvious systematic biases. These facts may indicate that the glint concept is appropriate for not only oceans but also land regions. Negative systematic biases in nadir and target modes have been found at most TCCON sites. The standard deviations of XCO₂ retrieved from target and nadir modes within the observation period are similar, and larger than those from glint mode. We also used the FTS site in Beijing, China, to assess the OCO-2 XCO₂ in 2016. This site is located in a typical urban area, which has been absent in previous studies. Overall, OCO-2 XCO₂ agrees well with that from FTS at this site. Such a study will benefit the validation of the newly launched TanSat products in China.     

Testing the Clausius–Clapeyron constraint on changes in extreme precipitation under CO<Subscript>2</Subscript> warming

Increases in extreme precipitation greater than in the mean under increased greenhouse gases have been reported in many climate models both on global and regional scales. It has been proposed in a previous study that whereas global-mean precipitation change is primarily constrained by the global energy budget, the heaviest events can be expected when effectively all the moisture in a volume of air is precipitated out, suggesting the intensity of these events increases with availability of moisture, and significantly faster than the global mean. Thus under conditions of constant relative humidity one might expect the Clausius–Clapeyron relation to give a constraint on changes in the uppermost quantiles of precipitation distributions. This study examines if the phenomenon manifests on regional and seasonal scales also. Zonal analysis of daily precipitation in the HadCM3 model under a transient CO₂ forcing scenario shows increased extreme precipitation in the tropics accompanied by increased drying at lower percentiles. At mid- to high-latitudes there is increased precipitation over all percentiles. The greatest agreement with Clausius–Clapeyron predicted change occurs at mid-latitudes. This pattern is consistent with other climate model projections, and suggests that regions in which the nature of the ambient flows change little give the greatest agreement with Clausius–Clapeyron prediction. This is borne out by repeating the analyses at gridbox level and over season. Furthermore, it is found that Clausius–Clapeyron predicted change in extreme precipitation is a better predictor than directly using the change in mean precipitation, particularly between 60°N and 60°S. This could explain why extreme precipitation changes may be more detectable then mean changes.     

Practical applicability of high frequency correction theories to CO<Subscript>2</Subscript> flux measured by a closed-path system

The theoretical correction of CO₂ fluxes for high frequency attenuation in closed-path systems was re-summarized and its applicability examined using both measurements obtained at an Asiaflux forest site and empirical transfer functions used in previous studies. For our measurement system, the theoretical transfer function was applicable to high frequency correction, even when condensation occurred in the sampling line. Further, in respect to some measurement systems described in previous studies, it was found that the theoretical function was potentially applicable along with the empirical functions used. Meanwhile, in some systems significant errors could not be resolved by re-estimation of the theories. In these systems, because of undefined buffering effects, the actual response lag time decided by the maximum covariance method or by measurement of the system response time using tracer gas was significantly different from the lag time calculated from the tube dimensions and the measured flow rate. If the average flow rate calculated by the actual lag time was used to determine the theoretical function, the theoretical function became closer to, and sometimes agreed with, the empirical function. Any remaining deviation from each function might be associated with pressure fluctuations, but this problem was unable to be examined here. The results suggested that an empirical formulation for each site is considered applicable rather than a theoretical approach, although the theories are being developed to practical application.