Parameterization of Albedo over Heterogeneous Surfaces in Coupled Land-Atmosphere Schemes for Environmental Modeling. Part I: Theoretical Background

In grid-based environmental models, the underlying surface consists of patches of solid and liquid parts and different plant communities, creating a very heterogeneous picture in the grid cell. In these cases, numerical modelers usually use a simple arithmetic average to determine the grid-cell albedo, a key variable in the parameterization of the land-surface radiative transfer over the grid cell. The object of this paper is to consider the assumptions for aggregating the albedo over a very heterogeneous surface where various surfaces occur at different heights, and, then propose a method for deriving a general expression for it. The suggested expression for the albedo is compared with the conventional approach, for the two-patches grid-cell with a simple geometrical distribution and different heights of its components. A numerical test is performed to compare the two approaches by numerical simulation of the evolution of the surface temperature over the particular grid-cell. Specifically, a one-dimensional land-surface model was applied to an isolated rocky grid-cell with a hole in the center; the model was forced with meteorological observations taken on July 17, 1999 in Philadelphia, PA.     

An Approach for the Aggregation of Aerodynamic Surface Parameters in Calculating the Turbulent Fluxes over Heterogeneous Surfaces in Atmospheric Models

Experimental evidence indicates that there is a significant departure of the wind profile above the underlying surface consisting of patches of solid and liquid parts, and plant communities with different morphological from that predicted by the logarithmic relationship, which gives the values larger than those observed. This situation can seriously affect the transfer of momentum, heat and water vapor from the surface fluxes into the atmosphere.The object of this paper is to generalize the calculation of the exchange of momentum between the atmosphere and a very heterogeneous surface, find a general equation for the wind speed profile in a roughness sublayer under neutral conditions, and, then, derive aggregated roughness length and displacement height over the grid cell. The suggested expression for the wind profile is compared with some earlier approaches, using a common parameterization of aerodynamic parameters over the grid cell, and the observations obtained at an experimental site in Philadelphia, PA.     

Occupational noise in printing companies

The extent of noise in five printing companies in Novi Sad, Serbia, was determined using TES-1358A Sound Analyzer with RS-232 Interface. The data on equivalent A-level (dBA), as well as, maximum and minimum sound pressure levels were collected. It was found that folders and offset printing units are the predominant noise sources, with the average L (eq) levels of 87.66 and 82.7 dBA, respectively. Forty percent of the machines produced noise levels above the limiting threshold level of 85 dBA, allowed by law. The noise in all printing companies was dominated by higher frequency noise, and the maximum level mostly appeared at 4,000 Hz. For offset printing machines and folders, the means of L (eq) levels exceeded the permissible levels given by NR-80 curve at higher frequencies. There are no published studies of occupational noise and hearing impairment of workers exposed to hazardous noise in printing industry in Serbia. More extensive studies are needed to determine the exact impact of noise on the workers. Technical and organizational measures in order to control noise and prevent noise exposure, and general hearing conservation program to protect workers, should be introduced in printing industry.     

Chemical transport models

Background, aim, and scope  Improving the parameterization of processes in the atmospheric boundary layer (ABL) and surface layer, in air quality and chemical transport models. To do so, an asymmetrical, convective, non-local scheme, with varying upward mixing rates is combined with the non-local, turbulent, kinetic energy scheme for vertical diffusion (COM). For designing it, a function depending on the dimensionless height to the power four in the ABL is suggested, which is empirically derived. Also, we suggested a new method for calculating the in-canopy resistance for dry deposition over a vegetated surface. Materials and methods  The upward mixing rate forming the surface layer is parameterized using the sensible heat flux and the friction and convective velocities. Upward mixing rates varying with height are scaled with an amount of turbulent kinetic energy in layer, while the downward mixing rates are derived from mass conservation. The vertical eddy diffusivity is parameterized using the mean turbulent velocity scale that is obtained by the vertical integration within the ABL. In-canopy resistance is calculated by integration of inverse turbulent transfer coefficient inside the canopy from the effective ground roughness length to the canopy source height and, further, from its the canopy height. Results  This combination of schemes provides a less rapid mass transport out of surface layer into other layers, during convective and non-convective periods, than other local and non-local schemes parameterizing mixing processes in the ABL. The suggested method for calculating the in-canopy resistance for calculating the dry deposition over a vegetated surface differs remarkably from the commonly used one, particularly over forest vegetation. Discussion  In this paper, we studied the performance of a non-local, turbulent, kinetic energy scheme for vertical diffusion combined with a non-local, convective mixing scheme with varying upward mixing in the atmospheric boundary layer (COM) and its impact on the concentration of pollutants calculated with chemical and air-quality models. In addition, this scheme was also compared with a commonly used, local, eddy-diffusivity scheme. Simulated concentrations of NO₂ by the COM scheme and new parameterization of the in-canopy resistance are closer to the observations when compared to those obtained from using the local eddy-diffusivity scheme. Conclusions  Concentrations calculated with the COM scheme and new parameterization of in-canopy resistance, are in general higher and closer to the observations than those obtained by the local, eddy-diffusivity scheme (on the order of 15–22%). Recommendations and perspectives  To examine the performance of the scheme, simulated and measured concentrations of a pollutant (NO₂) were compared for the years 1999 and 2002. The comparison was made for the entire domain used in simulations performed by the chemical European Monitoring and Evaluation Program Unified model (version UNI-ACID, rv2.0) where schemes were incorporated.