湍流通量参数化方案的非迭代方法研究

基于Högström (1996) 和Beljaars et al.(1991) 的研究工作, 沿用Louis et al.(1982) 和Launiainen (1995) 的思路, 本文采用多元回归分析方法, 研发了一种采用非迭代方法的湍流通量参数化方案。该方案直接用整体理查森数、 空气动力学粗糙度长度和热力学粗糙度长度对稳定度参数进行参数化, 从而避免了通过循环迭代计算Monin－Obukhov长度。该方案不仅有效地节省了CPU计算时间, 而且其计算结果与迭代方案 (BHH方案) 的计算结果非常接近。     

A mixed spectral finite-difference model for pollutant concentrations over a hill

Given incident logarithmic profiles of wind and pollutant concentration above a rough, absorbing surface, the three-dimensional distribution of pollutant concentration over a hill of gentle slope is determined from a linearized model. The model is applied in neutrally stratified flow, without chemistry, and is integrated using spectral methods in the horizontal and a finite-difference scheme in the vertical. This approach allows for flexibility in choosing a closure scheme and a variety of surface boundary conditions. This was not possible in the analytic approach of Padro (1987) who added pollutant concentration and flux to the MS3DJH/1 model of Walmsley et al. (1980). The present model requires as input the turbulent kinetic energy, E, dissipation, , and the perturbation vertical velocity, w, from the three-dimensional boundary-layer flow model of Beljaars et al. (1987), hereinafter referred to as MSFD, The latter model also supplies wind velocity perturbations at the upper boundary, as input to upper boundary conditions on the pollutant flux perturbations.The present study describes applications of the model to idealized terrain features: isolated two- and three-dimensional hills and ridges and an infinite series of ridges. (Application to real terrain, however, presents no difficulties.) Comparisons were made with different (though uniform) surface roughnesses. Tests were performed to examine the effect of upstream terrain features in the periodic domain and to illustrate the importance of the vertical resolution of the output for interpreting results from the sinusoidal terrain case.     

Shear stress results from a mixed spectral finite-difference model: Application to the Askervein Hill project data

Beljaars et al. (1987) developed a model for neutrally stratified boundary-layer flow over roughness changes and topography. It has been discovered that a constant parameter, , was missing in the algebraic-stress closure equations of their paper. This omission also occurred in the MSFD model code but only affects the Askervein Hill shear-stress results for the E-- turbulence closure in Beljaars et al. It also changes the stress results of Karpik (1988), but not his conclusions regarding the robustness of his improved numerical scheme. The present paper demonstrates the effect of the omission of the parameter, , and tests the sensitivity of the model to variations in its value. The new results are compared with the data and model results of Zeman and Jensen (1987).     

Do large-area-average temperature series have an urban warming bias?

Conclusions The arguments presented by Wood and his criticisms of the methods used by Jones et al. are largely fallacious and are generally based on misconceptions and unwarranted assumptions. This does not, of course, mean that the Jones et al. data are perfect. Jones et al. (1986a, b) have clearly stated that, at the regional level, there may be residual uncertainties in their gridded data set. Furthermore, they have noted that there is a residual uncertainty in the global-mean change since late last century of ±0.2 °C (Wigley et al., 1986), although the main reasons for this uncertainty do not relate to urban warming. Further studies of possible urban warming biases are certainly warranted, particularly at the regional scale. However, for averages over continental-scale areas upwards, it is unlikely that any significant urban warming bias remains.     

The effect of an oil slick on the albedo of the sea surface

Values of sea surface albedo estimated from a 3-day data set agree reasonably well with the findings of Payne (1972). However, when an oil slick moved over the observational site, the value of albedo suddenly jumped to 0.57. The effect of turbidity on the albedo over a lake has already been reported (Sadhuram et al., 1988).     

Numerical Simulation of Flow over Two-Dimensional Hills Using a Second-Order Turbulence Closure Model

We study turbulent flow over two-dimensional hills. The Reynolds stresses are represented by a second-order closure model, where advection, diffusion, production and dissipation processes are all accounted for. We solve a full set of primitive non-hydrostatic dynamic equations for mean flow quantities using a finite-difference numerical method. The model predictions for the mean velocity and Reynolds stresses are compared with the measured data from a wind-tunnel experiment that simulates the atmospheric boundary layer. The agreement is good. The performance of the second-order closure model is also compared withthat of lower level turbulence models, including the eddy-viscositymodel and algebraic Reynolds stress models. It is concluded that thepresent closure is a considerable improvement over the other modelsin representing various physical effects in flow over hills. Thefeasibility of running a finite-difference numerical simulationincorporating a full second-order closure model on an IBM workstationis also demonstrated.     

Comment: On the need for validation of the Jones et al. temperature trends with respect to urban warming

The Jones et al. hemispheric and global temperature trends and the methodology used to detect and correct for urban heat island effects are examined in detail. The results of this review suggest that there is still the possibility of significant urban warming bias remaining in the hemispheric and global averages. The actual extent of any such bias is unknown but warrants full investigation, since undetected and/or uncorrected urban warming could account for some of the hemispheric and global warming reported by Jones et al. Several research activities are suggested for resolving this uncertainty, including the rigorous application of urban warming correction factors to the data set used by Jones et al. and the development of geographically representative rural temperature series.Office of Technology Assessment, United States Congress, Washington, D.C. 20510, U.S.A. The views expressed are those of the author and not necessarily those of the OTA, Technology Assessment Board, or U.S. Congress. The article is based on independent research of the author and not on work conducted for the OTA.     

Fluctuations of the horizontal wind under unstable conditions

The similarity relations for _u/u_* proposed by Panofskyet al. (1977) and Højstrup (1982) have been verified using eddy-correlation data collected during the EFEDA-experiment, conducted over the extensive plain of La Mancha (Spain), where vine plants form a primary crop. Also, the standard Monin-Obukhov relation is considered. It is found that the expressions by Panofskyet al. and Højstrup both yield almost identical results, and are better than the Monin-Obukhov expression. Also, _u measured with a cup anemometer obeys the similarity expressions well, provided that the effect of the variation of wind direction on propellor wind speed is accounted for. The relationship of Panofskyet al. works rather well even when the boundary-layer height scaleh is replaced by a fixed height,h _c. Best results were obtained forh _c=1800 m. This height scale is possibly associated with the horizontal variability in the surface sensible heat flux pattern.     

The effect of gas-phase chemistry on aqueous-phase sulfur dioxide oxidation rates

The rates and mechanisms of both gas and liquid phase reactions for the oxidation of sulfur dioxide play an important role in the production of atmospheric acids and aerosol particles. Rhodeet al. (1981) concluded that sulfate production rates were highly non-linear functions of sulfur dioxide emission rates. Their modelling study used an HO _x termination mechanism for the HO—SO₂ reaction in the gas-phase. Stockwell and Calvert (1983) determined that one of the products of the overall reaction of HO with sulfur dioxide was an HO₂ radical. The National Research Council (1983) using a version of the Rhodeet al. (1981) model modified to include HO₂ production from the HO—SO₂ reaction concluded that sulfate production becomes much more linear with respect to reductions in sulfur dioxide emissions. However, the cause of this increased linearity was not explained by the National Research Council report. It is demonstrated that the increased linearity is due to the coupling of gas-phases and aqueous phase chemistry. The gas-phase sulfur dioxide oxidation mechanism has a very significant effect on hydrogen perodide production rates.     

RETRIEVING OVER-OCEAN PRECIPITABLE WATER FROM THE SSM/I MEASUREMENTS WITH SEVERAL REGRESSION ALGORITHMS*

A collocated SSM/I and radiosonde measurement data set provided by the NASDA(Japan) was used to retrieve the total precipitable water(PW) over oceans.The retrieval results obtained with several regression algorithms were compared against the radiosonde measurements.It is shown that:(a) the routinely operational algorithm of Alishouse et al.(1990) yields significant underestimation in high PW regime and overestimation in low PW regime;(b) a cubic correction by Colton and Poe(1994) is not sufficient and globally improves slightly the retrieval results;and(c) the regression algorithm with the form of brightness temperature(T_b) function In(280-T_b) gives a little largely scattered retrievals in whole PW range but without considerable over-and underestimates in low and high PW regimes.To improve the estimation of the oceanic precipitable water from the SSM/I measurements,a composite algorithm with different forms of T_b function in low.medium and high PW regimes is proposed and tested.     

Pan coefficient (K p) estimation under uncertainty on fetch

The FAO Penman–Monteith (F-PM) method is a frequently applied approach for calculating the daily reference evapotranspiration (ET₀). This method requires long records of meteorological data, which makes it quite hard to employ in locations with no or limited available data. Evaporation pans are widely used to estimate the reference ET₀, but this method requires reliable estimates of the pan coefficient (K _p). The objectives of this study were to determine the proper values of monthly and annual K _p, as well as the best method among those available for the estimation of K _p values in the study area. Measured weather data from 1992 to 2006 were obtained from 18 stations in the North and Northwest of Iran. Daily ET₀ calculated using methods by Bernardo et al. and Pereira et al. were compared with those calculated by the F-PM method. The employed methods at all stations, except those located in the north of the study area with high relative humidity, overestimated the ET₀ compared to the F-PM method. The constant parameters of these methods were optimized by a trial and error scheme to minimize the root mean square error. The results indicated that modified K _p coefficients from Bernardo et al.’s method ranged between 0.41 and 0.87 and the optimal coefficient of Pereira et al.’s method ranged between 0.49 and 0.95. Modified monthly K _p from Bernardo et al.’s method ranged between 0.3 and 1.07 and those from Pereira et al.’s method ranged between 0.4 and 1.18. Modified K _p of the methods by Bernardo et al. and Pereira et al. showed the higher estimation accuracy of daily ET₀ values. In general, the performance of the modified K _p of Bernardo et al.’s method was higher than Pereira et al.’s method for all stations. Thus, in the study region and under the same climatic conditions [in areas with only pan evaporation (E _p) records], the use of climatic monthly modified K _p to calculate ET₀ based on class A E _p is recommended.     

Incorporation of Non-Linear Effective Cross-Section Parameterizations into a Fast Photolysis Computation Code (Fast-J)

We have investigated several methods of incorporating laboratorymeasurements of effective cross-sections (quantum yield timesabsorption cross-section) that are non-linear functions oftemperature, pressure, and wavelength, into the Fast-J photolysis codeof Wild et al. (2000). The best method was found to beparameterization of the solar weighted effective cross-section usingthe form of the effective cross-section, with numerical fitting of theparameterization constants. Using this technique on the acetone dataof Gierczak et al. (1998) produces aparameterization which fits all points to within 1%. This methodshould also benefit other photolysis codes with large wavelength bins.     

OH-initiated degradation of several hydrocarbons in the atmosphere simulation chamber SAPHIR

Degradation of isoprene, m-xylene, n-octane, propene, and methacrolein by hydroxyl radicals has been studied in the simulation chamber SAPHIR under burden of trace gases as they are typical for the moderately polluted planetary boundary layer. Measured time series of the hydrocarbon mixing ratios and the OH concentrations were used to determine the rate constants. The hydrocarbons were measured with gas chromatography and proton transfer reaction mass spectrometry. OH was measured with the Jülich DOAS (differential optical absorption spectroscopy) instrument. In all cases except methacrolein good agreement was found with the reference rate constants taken from the Master Chemical Mechanism (MCM3.1). The data for methacrolein are consistent with the results of Karl et al. (J. Atmos. Chem 55, 2006, doi:) who reported a 12% smaller value. The degradation of hydrocarbons provides an independent method to analyse precision and accuracy of the OH measurements. A precision of better than 4% over a period of nearly 4 months was found. The accuracy is within the limitations given by the light absorption cross section of OH. Both results are consistent with earlier results by Hausmann et al. (J. Geophys. Res. 102:16011–16022, 1997).     

Spectral nudging sensitivity experiments in a regional climate model

In this study, the scale selective bias correction (SSBC) method described by Kanamitsu et al. (2010) is further examined by considering the full wind nudging and the vertically weighted damping coefficient. The 2001 June?CJuly?CAugust RSM simulation over a relatively large domain covering much of the Asian continent, the northern part of Australia, and the Indian and western Pacific oceans was the main focus. The full wind nudging shows wind fields closer to the driving global analysis. However, it leads to significantly distorted fields (e.g., temperature and geopotential height) aloft, accompanying excessive precipitation over the western Pacific. The gradual reduction of vorticity nudging from the model top to the ground surface improves rainfall patterns without a discernible distortion of large-scale fields. Further evaluation of a 10-year-summer simulation over East Asia confirmed that this revised SSBC method improves the monsoonal rainfall against the method of Kanamitsu et al. It is therefore concluded that vorticity nudging alleviates largescale errors by maintaining the near geostrophic balance between mass and winds. The reduction of this nudging factor in the lower troposphere allows the ageostrophic component of wind to develop as in nature, which leads to the improvement of precipitation.     

Carbon tetrachloride lifetimes and emissions determined from daily global measurements during 1978–1985

Seven years of daily gas chromatographic measurements of CCl₄ at the five globally distributed ALE/GAGE surface sites are reported. It is determined that CCl₄ has been accumulating in the atmosphere at a rate of 1.3±0.1%/yr over the period 1978–1985 and that the releases of CCl₄ into the atmosphere have remained fairly constant, with the smallest releases in 1981–1982. Using an inversion scheme based on a nine box model of the atmosphere, we infer a CCl₄ lifetime of approximately 40 yr, an inventory on 1 July 1978 of (2.08±0.07)×10⁹ kg and an average rate of release over the period 1978–1985 of (0.9±0.9)×10⁷ kg/yr. These results produce excellent agreement with a release scenario derived from global production estimates for CCl₄ and the major CCl₄ byproduct, the chlorofluorocarbons. However, to obtain this consistency, it is necessary that our current ALE/GAGE absolute calibration standard be reduced approximately 25% thus bringing it into agreement with measurements by Yokohata et al. (1985) and Hanst et al. (1975).     

Internal boundary-layer height formulae — A comparison with atmospheric data

The height of the internal boundary layer (IBL) downwind of a step change in surface roughness is computed using formulae of Elliott (1958), Jackson (1976) and Panofsky and Dutton (1984). The results are compared with neutral-stratification atmospheric data extracted from the set of wind-tunnel and atmospheric data summarized by Jackson (1976) as well as neutral-stratification data presented by Peterson et al. (1979) and new data measured at Cherrywood, Ontario. It is found that the Panofsky-Dutton formulation gives the least root-mean-square (RMS) absolute errors for atmospheric applications.     

A new method for estimation of near surface specific humidity over global oceans

Summary A new methodology for deriving daily averages of near surface specific humidity (Q _a ) is developed. Remotely sensed parameters, total water vapor (W), boundary-layer water vapor (W _b ), and sea-surface temperature (SST) are used to derive Q _a . Genetic algorithm (GA) is used to find the empirical function relating the input (W, W _b , and SST) and output (Q _a ) parameters. The input data consist of 2 years (1999–2000) of daily W, W _b , from SSM/I (Special Sensor Microwave Imager), and SST data from AVHRR (Advanced Very High Resolution Radiometer). COADS (Comprehensive Ocean Atmospheric Data Set) observations of Q _a are used to develop and evaluate the new methodology. The performance of the algorithm is measured with COADS observations, which are not used in the development phase. The global mean rms error for daily averages Q _a is 1.5 ± 0.40 g/kg. Slightly higher discrepancies between Q _a derived from the new method and COADS observations are found over the Northwestern Pacific, North Atlantic oceans and Arabian Sea. This method improves upon the humidity retrieval of Liu (1986), Schulz et al (1993), and Chou et al (1997).     

Possible climatic impacts of land cover transformations,with particular emphasis on tropical deforestation

The climatic impact of albedo changes associated with land-surface alterations has been examined. The total surface global albedo change resulting from major land-cover transformations (i.e. deforestation, desertification, irrigation, dam-building, urbanization) has been recalculated, modifying the estimates of Sagan et al., (1979). Tropical deforestation (11.1 million ha yr^-1, or 0.6% yr^-1, Lanly, 1982) ranks as a major cause of albedo change, although uncertainties in the areal extent of desertification could conceivably render this latter process of similar significance. The maximum total global albedo change over the last 30 yr for the various processes lies between 0.000 33 and 0.000 64, corresponding to a global temperature decrease of between 0.06 K and 0.09 K (scaled from the 1-D radiative convective model of Hansen et al., 1981), which falls well below the interannual and longer period variability.An upper bound to the impact of tropical deforestation was obtained by concentrating all vegetation change into a single region. The magnitude of this modification is equivalent to 35–50 yr of global deforestation at the current rate, but centered on the Brazilian Amazon. The climatic consequences of such tropical deforestation were simulated, using the GISS GCM (Hansen et al., 1983). In the simulation, a total area of 4.94 × 10⁶ km² of tropical moist forest was removed and replaced by a grass/crop cover. Although surface albedo increased from 0.11 to 0.19, the effect upon surface temperature was negligible. However, other climate parameters were altered. Rainfall decreased by 0.5–0.7 mm day^-1 and both evapotranspiration and total cloud cover were reduced. The absence of a temperature decrease in spite of the increased surface albedo arises because the reduction in evapotranspiration has offset the effects of radiative cooling. The decrease in cloud cover also counteracts the increase in surface albedo. These locally significant changes had no major impact on regional (Hadley or Walker cells) or the global circulation patterns.We conclude that the albedo changes induced by current levels of tropical deforestation appear to have a negligibly small effect on the global climate.     

An aircraft study of turbulence dissipation rate and temperature structure function in the unstable marine atmospheric boundary layer

The dissipation rate of turbulent kinetic energy, , and the temperature structure function parameter, C _T ², have been measured over water from the near surface (Z = 3 m) to the top of the boundary layer. The near surface values of and C _T ² were used to calculate the velocity and temperature Monin-Obukhov scaling parameters u _* and T _*. The data collected during unstable lapse rates were used to evaluate the feasibility of extrapolating the values of and C _T ² as a function of height with empirical scaling formulae. The dissipation rate scaling formula of Wyngaard et al. (l971 a) gave a good fit to an average of the data for Z < 0.8 Z _i. In the surface layer the scaling formula of Wyngaard et al. (1971b) disagreed with the C _T ² values by as much as 50%. This disagreement is due to an unexpected reduction in the measured values of C _T ² forZ < 30 m. At this point it is not clear if the discrepancy is a unique property of the marine boundary layer or if it is simply some unknown instrumental or analytical problem. The mixed layer scaling results were similar to the overland results of Kaimal et al. (1976).     

Evidence for zonally-trapped propagating waves in the eastern atlantic from satellite sea surface temperature observations

It has been hypothesized (Moore et al., 1978; O'Brien et al., 1978), that equatorial upwelling and subsequent coastal upwelling on the eastern boundary of the Atlantic Ocean are the result of eastward propagating equatorially trapped Kelvin waves in the Atlantic. Concurrent satellite and ship sea surface temperature observations taken during the GATE experiment permit validation of the satellite data as well as relating sea surface temperature (SST) variability to the local current dynamics. A method based on cross-correlations and cross-spectra of the SST field at various locations is utilized to test the Kelvin wave hypothesis. Significant periodic variation of time lags in the SST variability in the eastern Atlantic is observed by the spectral techniques. Satellite data for the 1974 summer show periodic variability which fits either eastward or westward propagating waves with 1 m s^-1 phase speed, i.e., SST supports the quasi-continuous presence of Kelvin or Yanai waves. We find no evidence for a seasonally solitary eastward propagating signal in the eastern Atlantic from SST.