Sensitivity of air quality model prediction to parameterization of vertical eddy diffusivity

This paper investigates the effects of vertical eddy diffusivities derived from the 3 different planetary boundary layer (PBL) schemes on predictions of chemical components in the troposphere of East Asia. Three PBL schemes were incorporated into a regional air quality model (RAQM) to represent vertical mixing process and sensitivity simulations were conducted with the three schemes while other options are identical. At altitudes <2km, all schemes exhibit similar skill for predicting SO₂ and O₃, but more difference in the predicted NO_x concentration. The Gayno–Seaman scheme produces the smallest vertical eddy diffusivity (Kz) among all schemes, leading to higher SO₂ and NO_x concentrations near the surface than that from the other 2 schemes. However, the effect of vertical mixing on O₃ concentration is complex and varies spatially due to chemistry. The Gayno–Seaman scheme predicts lower O₃ concentrations than the other two schemes in the parts of northern China (generally VOC-limited) and higher ones in most parts of southern China (NO_x-limited). The Byun and Dennis scheme produces the largest mixing depth in the daytime, which bring more NO_x into upper levels, and the mixing depth predicted by the Gayno–Seaman scheme is the smallest, leading to higher NO_x and lower O₃ concentrations near the surface over intensive emission regions.     

The impact of chemical lateral boundary conditions on CMAQ predictions of tropospheric ozone over the continental United States

A sensitivity study is performed to examine the impact of lateral boundary conditions (LBCs) on the NOAA-EPA operational Air Quality Forecast Guidance over continental USA. We examined six LBCS: the fixed profile LBC, three global LBCs, and two ozonesonde LBCs for summer 2006. The simulated results from these six runs are compared to IONS ozonesonde and surface ozone measurements from August 1 to 5, 2006. The choice of LBCs can affect the ozone prediction throughout the domain, and mainly influence the predictions in upper altitude or near inflow boundaries, such as the US west coast and the northern border. Statistical results shows that the use of global model predictions for LBCs could improve the correlation coefficients of surface ozone prediction over the US west coast, but could also increase the ozone mean bias in most regions of the domain depending on global models. In this study, the use of the MOZART (Model for Ozone And Related chemical Tracers) prediction for CMAQ (Community Multiscale Air Quality) LBC shows a better surface ozone prediction than that with fixed LBC, especially over the US west coast. The LBCs derived from ozonesonde measurements yielded better O₃ correlations in the upper troposphere.     

Boundary Layer Evolution and its Influence on Ground-Level Ozone Concentrations

Accurate determination of the planetary boundary layer (PBL) height (i.e., mixing height (MH)) is critical to properly simulating pollutant levels with the grid-based photochemical models. In this study, the daytime mixing heights based on the parcel and bulk Richardson number methods are compared with those obtained directly from a numerical mesoscale meteorological model in an effort to evaluate the uncertainties in the estimation of the PBL evolution. Mixing heights are estimated from hourly outputs of meteorological variables of the Penn State/NCAR Mesoscale Model Version 3.3 (MM5V3) with two PBL schemes (Blackadar and Gayno-Seaman) during July 1999 over Philadelphia, PA. An analysis of the diurnal variation in the urban PBL and its influence on ground-level ozone (O₃) levels is presented in this paper. The results indicate that on average, the MHs determined from the bulk Richardson number were larger than those estimated from the parcel method. The MHs from the MM5V3 output were much smaller than those derived from the parcel and bulk Richardson number methods, especially for the Gayno-Seaman scheme that is based on turbulent kinetic energy. The MH and ground-level O₃ concentration have been found to be twice as much on episode days than on non-episode days. The average hourly MH growth rate and O₃ tendency (i.e. rate of change in O₃) were largest during the morning hours (0700 to 1000 eastern standard time (EST)), suggesting that vertical mixing contributes significantly to the accumulation of ground-level O₃ in urban areas in the morning hours.     

Stability and Mixing of a Vertical Axisymmetric Buoyant Jet in Shallow Water

The stability, mixing and effect of downstream control on axisymmetric turbulent buoyant jets discharging vertically into shallow stagnant water is studied using 3D Reynolds-averaged Navier–Stokes equations (RANS) combined with a buoyancy-extended k –ε model. The steady axisymmetric turbulent flow, temperature (or tracer concentration) and turbulence fields are computed using the finite volume method on a high resolution grid. The numerical predictions demonstrate two generic flow patterns for different turbulent heated jet discharges and environmental parameters (i) a stable buoyant discharge with the mixed fluid leaving the vertical jet region in a surface warm water layer; and (ii) an unstable buoyant discharge with flow recirculation and re-entrainment of heated water. A stratified counterflow region always appears in the far-field for both stable and unstable buoyant discharges. Provided that the domain radius L exceeds about 6H, the near field interaction and hence discharge stability is governed chiefly by the jet momentum length scale to depth ratio l_M/H, regardless of downstream control. The near field jet stability criterion is determined to be l_M/H = 3.5. A radial internal hydraulic jump always exists beyond the surface impingement region, with a 3- to 6-fold increase in dilution across the jump compared with vertical buoyant jet mixing. The predicted stability category, velocity and temperature/concentration fields are well-supported by experiments of all previous investigators.     

Numerical Investigation of the Effects of Boundary-Layer Evolution on the Predictions of Ozone and the Efficacy of Emission Control Options in the Northeastern United States

This paper examines the effects of two different planetary boundary-layer (PBL) parameterization schemes – Blackadar and Gayno–Seaman – on the predicted ozone (O₃) concentration fields using the MM5 (Version 3.3) meteorological model and the MODELS-3 photochemical model. The meteorological fields obtained from the two boundary-layer schemes have been used to drive the photochemical model to simulate O₃ concentrations in the northeastern United States for a three-day O₃ episodic period. In addition to large differences in the predicted O₃ levels at individual grid cells, the simulated daily maximum 1-h O₃ concentrations appear at different regions of the modeling domain in these simulations, due to the differences in the vertical exchange formulations in these two PBL schemes. Using process analysis, we compared the differences between the different simulations in terms of the relative importance of chemical and physical processes to O₃ formation and destruction over the diurnal cycle. Finally, examination of the photochemical model's response to reductions in emissions reveals that the choice of equally valid boundary-layer parameterizations can significantly influence the efficacy of emission control strategies.     

Transport and Diffusion of Ozone in the Nocturnal and Morning Planetary Boundary Layer of the Phoenix Valley

The evolution of ozone (O ₃) in the nocturnal and morning-transitional planetary boundary layer (PBL) of the Phoenix valley was measured as a part of the `Phoenix Sunrise Experiment 2001' of the U.S. Department of Energy conducted in June 2001. The goal of the field program was to study the transport, distribution and storage of ozone and its precursors in the urban boundary layer over a diurnal cycle. The ground level O <sub>3</sub> as well as mean meteorological variables and turbulence were measured over the entire period, and vertical profiling (using a tethered balloon) was made during the morning transition period. Approximately half of the observational days showed the usual diurnal cycle of high O <sub>3</sub> during the day and low O <sub>3</sub> at night, with nitrogen oxides (NO <sub>x</sub> = NO <sub>2</sub> + NO) showing an out of phase relationship with O <sub>3</sub>. The rest of the days were signified by an anomalous increase of O <sub>3</sub> in the late evening ( 2200 LST), concomitant with a sudden drop of temperature, an enhancement of wind speed and Reynolds stresses, a positive heat flux and a change of wind direction. NO <sub>x</sub> measurements indicated the simultaneous arrival of an `aged' air mass, which was corroborated by the wind predictions of a mesoscale numerical model. In all, the results indicate that the recirculation of O ₃ rich air masses is responsible for the said high-O ₃ events. Such air masses are produced during the transport of O ₃ precursors by the upslope flow toward mountainous suburbs during the day, and they return back to the city at night via downslope winds (i.e. mountain breeze). The corresponding flow patterns, and hence the high-O ₃ events, are determined by background meteorological conditions. The vertical profiling of O ₃ and flow variables during the morning transition points to a myriad of transport, mixing and chemical processes that determine the fate of tropospheric O ₃. How well such processes are incorporated and resolved in predictive O ₃ models should determine the accuracy of their predictions.     

Double-average mean flow and local turbulence intensity profiles from PIV measurements for an open channel flow with rigid vegetation

This contribution presents particle image velocimetry measurements for an open channel stationary uniform and fully developed flow of water over a horizontal flat bed of uniform glass beads in presence of a staggered array of vertical cylindrical stems. The main objective was to explore and quantify the influence of the stems-to-flow relative submergence, h _v /h, over the mean flow and local turbulence intensities. A comparison with measurements for the non-vegetated flow over the same granular bed is presented. Results indicate a remarkable influence of h _v /h over the whole flow field. The time-average mean flow presents a strong spatial variation in the layer of the flow occupied by the stems. The local velocity fluctuations are strongly affected by the presence of the stems, with regions in between the stems where they reach peaks that are several times larger than those encountered in the flow in absence of vegetation. The turbulence intensity profiles are noticeably different when compared to those measured in the non-vegetated flow conditions. From previous works it was possible to derive an equation for the mean velocity, U _v , of the flow through the vegetated layer of height h _v . The prediction of this equation is in good agreement with the uniform value for the double-average longitudinal velocity profile in this layer. A final brief discussion about the possible impact of these vegetated-flow features on the sediment transport is presented.     

Experimental Study of the Criteria of Flow Laminarization in Two-Dimensional Dense Gas Plumes

Laminarization of flow in a two-dimensional dense gas plume was experimentally investigated in this study. The plume was created by releasing CO₂ through a ground-level line source into a simulated turbulent boundary layer over an aerodynamically rough surface in a meteorological wind tunnel. The bulk Richardson number (Ri_*), based on negative plume buoyancy, plume thickness, and friction velocity, was varied over a wide range so that the effects of stable stratification on plume laminarization could be observed. A variety of ambient wind speeds as well as three different sizes of roughness arrays were used so that possible effects of roughness Reynolds number (Re_*) on plume laminarization could also be identified. Both flow visualization methods and quantitative measurements of velocity and intermittency of turbulence were used to provide quantitative assessments of plume laminarization.Flow visualization provided an overall picture of how the plume was affected by the negative buoyancy. With increasing Ri_*, both the plume depth and the vertical mixing were significantly suppressed, while upstream propagation of the plume from the source was enhanced. The most important feature of the flow revealed by visualization was the laminarization of flow in the lower part of the plume, which appeared to be closely related to both Ri_* and Re_*.Measurements within the simulated dense gas plumes revealed the influence of the stable stratification on mean velocity and turbulence intensity profiles. Both the mean velocity and turbulence intensity were significantly reduced near the surface; and these reductions systematically depended on Ri_*. The roughness Reynolds number also had considerable influence on the mean flow and turbulence structure of the dense gas plumes.An intermittency analysis technique was developed and applied to the digitized instantaneous velocity signals. It not only confirmed the general flow picture within the dense plume indicated by the flow visualization, but also clearly demonstrated the changes of flow regime with variations in Ri_* and Re_*. Most importantly, based on this intermittency analysis, simple criteria for characterizing different flow regimes are formulated; these may be useful in predicting when plume laminarization might occur.     

Estimation of far-field horizontal and vertical turbulent diffusion coefficients from the concentration field of a wastewater plume near the Akashi Strait

We carried out a field study of the plume discharged by a near-shore wastewater outfall near the Akashi Strait, Japan. Using an Acoustic Doppler Current Profiler and a tow-body CTD, we measured the near-surface salinity and temperature fields in the region throughout an M₂ tidal cycle. We filtered the data in T–S space to remove water masses other than the wastewater, and then used the adiabatic mixing assumption to calculate the concentration of wastewater in the far field of this plume. Averaging the T–S fields of repeated surveys over a time period during which the tidal regime did not change substantially, allowed comparison of the time-averaged plume with the analytical solution for a plume diffusing in both the horizontal and vertical dimensions. The resulting vertical turbulent diffusion coefficients agreed well with those resulting from Thorpe scales determined via a vertically-profiling CTD, as well as with the canonical value for open channel flow of D _z = 0.067hu _*. The corresponding horizontal turbulent diffusion coefficients, however, were two orders of magnitude larger than those typically observed in straight channels, and an order of magnitude larger than those observed in meandering rivers. This is likely a result of enhanced horizontal mixing due to barotropic eddies generated by the interaction of strong tidal flow with headlands and levees, as well as due to the time-varying nature of tidal flow, and baroclinic spreading of the buoyant wastewater plume.     

Eddy diffusivity: a single dispersion analysis of high resolution drifters in a tidal shallow estuary

In an estuary, mixing and dispersion resulting from turbulence and small scale fluctuation has strong spatio-temporal variability which cannot be resolved in conventional hydrodynamic models while some models employs parameterizations large water bodies. This paper presents small scale diffusivity estimates from high resolution drifters sampled at 10 Hz for periods of about 4 h to resolve turbulence and shear diffusivity within a tidal shallow estuary (depth <3 m). Taylor’s diffusion theorem forms the basis of a first order estimate for the diffusivity scale. Diffusivity varied between 0.001 and 0.02 m²/s during the flood tide experiment. The diffusivity showed strong dependence (R² > 0.9) on the horizontal mean velocity within the channel. Enhanced diffusivity caused by shear dispersion resulting from the interaction of large scale flow with the boundary geometries was observed. Turbulence within the shallow channel showed some similarities with the boundary layer flow which include consistency with slope of 5/3 predicted by Kolmogorov’s similarity hypothesis within the inertial subrange. The diffusivities scale locally by 4/3 power law following Okubo’s scaling and the length scale scales as 3/2 power law of the time scale. The diffusivity scaling herein suggests that the modelling of small scale mixing within tidal shallow estuaries can be approached from classical turbulence scaling upon identifying pertinent parameters.     

Large-Eddy Simulations of the Atmospheric Boundary Layer Using a New Subgrid-Scale Model – II. Weakly and moderately stable cases

A series of simulations under weakly to moderately stable boundary layers (SBLs) have been performed using the proposed subgrid-scale (SGS) model implemented into the Terminal Area Simulation System (TASS). The proposed SGS model incorporates some aspects of the two-part eddy viscosity SGS model of Sullivan et al. (1994) and further refinements which include the dependence of SGS mixing length on stratification, two-part separation of the SGS eddy diffusivity of heat, and more realistic empirical forms of Monin–Obukhov similarity functions. The potential temperature profiles from simulations clearly show a three-layer structure: a stable surface layer of strong gradients, a middle layer of small gradients, and an inversion layer on the top. The wind speed profiles show the formation of low level jet (LLJ). However, the sub-layer structures under moderately SBLs differ from those under weakly SBLs. Both the momentum and heat fluxes decrease almost linearly in the lower part of the SBL. The near surface values of the normalized turbulent kinetic energy (TKE/u _* ²) in all simulations are about 4 which is much less than the typical value of 5.5 under the neutral condition. The decay of turbulence first occurs in the area with large values of Richardson number (R _i<0.2). Generally, instantaneous values of the TKE and R _i at the various grid points are negatively correlated, but there is not a unique relationship between the two parameters.     

Diurnal gut pigment rhythm and metabolic rate of<Emphasis Type="Italic"> Calanus euxinus</Emphasis> in the Black Sea

The vertical distribution, diel gut pigment content and oxygen consumption of Calanus euxinus were studied in April and September 1995 in the Black Sea. Gut pigment content of C. euxinus females was associated with diel vertical migration of the individuals, and it varied with depth and time. Highest gut pigment content was observed during the nighttime, when females were in the chlorophyll a (chl a) rich surface waters, but significant feeding also occurred in the deep layer. Gut pigment content throughout the water column varied from 0.8 to 22.0 ng pigment female^–1 in April and from 0.2 to 21 ng pigment female^–1 in September 1995. From the diel vertical migration pattern, it was estimated that female C. euxinus spend 7.5 h day^–1 in April and 10.5 h day^–1 in September in the chl a rich surface waters. Daily consumption by female C. euxinus in chl a rich surface waters was estimated by taking into account the feeding duration and gut pigment concentrations. Daily carbon rations of female C. euxinus, derived from herbivorous feeding in the euphotic zone, ranged from 6% to 11% of their body carbon weight in April and from 15% to 35% in September. Oxygen consumption rates of female and copepodite stage V (CV) C. euxinus were measured at different temperatures and at different oxygen concentrations. Oxygen consumption rates at oxygen-saturated concentration ranged from an average of 0.67 g O₂ mg^–1 dry weight (DW) h^–1 at 5°C to 2.1 g O₂ mg^–1 DW h^–1 at 23°C for females, and ranged from 0.48 g O₂ mg^–1 DW h^–1 at 5°C to 1.5 g O₂ mg^–1 DW h^–1 at 23°C for CVs. The rate of oxygen consumption at 16°C varied from 0.62 g O₂ mg^–1 DW h^–1 at 0.65 mg O₂ l^–1 to 1.57 g O₂ mg^–1 DW h^–1 at 4.35 mg O₂ l^–1 for CVs, and from 0.74 g O₂ mg^–1 DW h^–1 at 0.57 mg O₂ l^–1 to 2.24 g O₂ mg^–1 DW h^–1 at 4.37 mg O₂ l^–1 for females. From the oxygen consumption rates, daily requirements for the routine metabolism of females were estimated, and our results indicate that the herbivorous daily ration was sufficient to meet the routine metabolic requirements of female C. euxinus in April and September in the Black Sea.Communicated by O. Kinne, Oldendorf/Luhe     

A Computational Fluid Dynamics Approach for Urban Area Transport and Dispersion Modeling

A simulation tool has been developed to model the wind fields, turbulence fields, and the dispersion of Chemical, Biological, Radiological and Nuclear (CBRN) substances in urban areas on the building to city blocks scale. A Computational Fluid Dynamics (CFD) approach has been taken that naturally accounts for critical flow and dispersion processes in urban areas, such as channeling, lofting, vertical mixing and turbulence, by solving the steady-state, Reynolds-Averaged Navier–Stokes (RANS) equations. Rapid generation of high quality cityscape volume meshes is attained by a unique voxel-based model generator that directly interfaces with common Geographic Information Systems (GIS) file formats. The flow and turbulence fields are obtained by solving the steady-state RANS equations using a collocated, pressure-based approach formulated for unstructured and polyhedral mesh elements. Turbulence modeling is based upon the Renormalization Group variant of the k–ε model (k–ε RNG). Neutrally buoyant simulations are made by prescribing velocity boundary condition profiles found by a power–law relationship, while turbulence quantities boundary conditions are defined by a prescribed mixing length in conjunction with the assumption of turbulence equilibrium. Dispersion fields are computed by solving an unsteady transport equation of a dilute gas, formulated in a Eulerian framework, using the velocity and turbulence fields found from the steady-state RANS solution. In this paper the model is explained and detailed comparisons of predicted to experimentally obtained velocity, turbulence and dispersion fields are made to neutrally stable wind tunnel and hydraulic flume experiments.     

A new view of dispersion in well-mixed estuaries

A two-dimensional simulation of Delaware estuary hydrodynamics has been constructed. This simulation has been achieved through a rational estimate of the character of natural turbulence. Non-homogeneous velocities, on the cross-section, are employed in two-dimensional, laterally homogeneous species mass balances. In turn, concentration profiles are interpreted in the form of classical, one-dimensional dispersion coefficients. Variation of dispersion as a function of both freshwater inflow and longitudinal distance was generated. Variation of dispersion in time within a tidal cycle was found to be insignificant while no significant variation from one tidal cycle to the next has been detected.The modeling process involves the solution of tractable equations by implicit numerical methods and is capable of being excited by a wide range of input conditions.A study of the sensitivity of dispersion due to vertical mass diffusion revealed that longitudinal mixing characteristics are inversely proportional to vertical eddy diffusivity and analysis of the numerical results showed the dispersion coefficient is essentially insensitive to variation of longitudinal mass diffusivity. This leads to the conclusion that turbulent diffusivity of mass in the longitudinal direction may be taken as constant for most purposes in the study of a two-dimensional species mass balance model.A field program was carried out near the Delaware Memorial Bridge to collect velocity profiles. Substantial portions of the scheme have been verified (i.e. one- and two-dimensional tidal dynamic models) through the use of these data.     

The three-dimensionality of obstructed shear flows

Obstructed shear flows (i.e. those over permeable media) are common in the environment. An archetypal example, flow over a submerged vegetation canopy, is investigated here. Like any flow through complex geometry, canopy flows are characterised by strong spatial gradients. The focus of this experimental study is the three-dimensionality of aquatic canopy flow, in particular that of the coherent interfacial vortices that govern mixing into and out of the canopy. It is shown here that the vortices have a finite lateral scale that is comparable to their vertical scale; both are of the order of the drag length scale of the canopy, (C _D a)⁻¹, where a is the frontal area density and C _D is a bulk drag coefficient. The finite lateral extent of the vortices generates strong lateral hydrodynamic gradients, both instantaneously and in the long-term. The instantaneous gradients, which can contribute greatly to the dispersion of dissolved and particulate species, are far more pronounced. Finally, the potential for canopies to generate differential roughness secondary circulation is examined. In the consideration of vertical scalar transport, this circulation can be of the same order as turbulent diffusion.     

A statistical model to forecast the profile of the index structure constant $${C_N^2}$$

A statistical study was performed over 145 profiles of meteorological balloons, equipped with microthermal sensors, from the ground to the midstratosphere. This study put into evidence the lognormal distribution of the fluctuations of the structure constant of temperature, the fluctuations of the buoyancy force and the vertical shear of the horizontal wind speed. We show that these quantities, computed over a large scale (100 m), are correlated up to the midstratosphere. A model is adjusted to estimate the optical turbulence strength from the macroscopic meteorological parameters. The model performances, to estimate and forecast the strength and the altitude of the optical turbulent layers, are quantified and compared with other already defined models. This model was shown to have the best performances and put into evidence a new relation to describe the turbulence on a large scale.     

Turbulence closure models and their application in RAMS

Turbulence closures are fundamental for modelling the atmospheric diffusion in numerical codes and the resulting eddy diffusivities are key parameters in describing the transport and dispersion in the boundary layer. In this work, four turbulence closure schemes have been applied for reproducing a neutral flow over schematic complex terrain using the meteorological model RAMS. Two of the closures, a one-equation (E-l) and a two-equations (E-) model, have been implemented in RAMS in alternative to the ones originally available. In these cases, an analytical method based on the similarity theory for the atmospheric surface layer and boundary layer is adopted to calculate the empirical constants of the turbulence closures. Some examples of numerical studies performed to simulate the flow and turbulence over a 3-D hill in wind-tunnel experiment in neutral stratification are presented and discussed. An intercomparison of simulations related to different closures is considered by analysing the main features of the flow over the hill and by comparing calculated vertical profiles of turbulent kinetic energy with measured ones.     

In situ respiration of benthic communities in Castle Harbor,Bermuda

In situ measurements of community respiration were made at two stations in Castle Harbor, Bermuda, during April and May, 1971. Total community respiration was 20.67 and 19.11 ml O₂ m^-2 h^-1 at Stations 1 and 2, respectively, in April. In May, respiration increased with water temperature to 26.99 and 24.56 ml O₂ m^-2 h^-1. Significant differences (P<0.05) existed between stations and sampling periods. Bacterial respiration was estimated from sediment treatment with streptomycin-SO₄. Values ranged from 7.71 to 8.72 ml O₂ m^-2 h^-1 in April and May, respectively. No significant difference existed between sampling periods or stations. No detectable chemical O₂ demand of the sediment, determined by a formalin treatment, was found. Total community respiration was further compartmentalized into macrofaunal, meiofaunal, and microfaunal-microfloral components. Both the macrofaunal and meiofaunal respiration was negligible compared to that estimated for the microfaunal-microfloral component.Contribution No. 2708 from the Woods Hole Oceanographic Institution, and Contribution No. 552 from the Bermuda Biology Station. This study was supported in part by National Science Foundation Grants GZ 1508 and GB 16161.     

Experimental analysis of Froude number effect on air entrainment in the hydraulic jump

A hydraulic jump is characterized by strong energy dissipation and mixing, large-scale turbulence, air entrainment, waves, and spray. Despite recent pertinent studies, the interaction between air bubbles diffusion and momentum transfer is not completely understood. The objective of this paper is to present experimental results from new measurements performed in a rectangular horizontal flume with partially developed inflow conditions. The vertical distributions of the void fraction and the air bubbles count rate were recorded for inflow Froude number Fr ₁ in the range from 5.2 to 14.3. Rapid detrainment process was observed near the jump toe, whereas the structure of the air diffusion layer was clearly observed over longer distances. These new data were compared with previous data generally collected at lower Froude numbers. The comparison demonstrated that, at a fixed distance from the jump toe, the maximum void fraction C _max increases with the increasing Fr ₁. The vertical locations of the maximum void fraction and bubble count rate were consistent with previous studies. Finally, an empirical correlation between the upper boundary of the air diffusion layer and the distance from the impingement point was derived.     

A Simple Model to Determine In Situ Mixing Height Growth from Surface Measurements

A generic In Situ Mixing Height Growth (IMG) model, capable of predicting the real-time growth of the mixed layer and its diurnal evolution from routinely observed simple surface meteorological is developed. The algorithm for the determination of temporally growing daytime mixing height includes both convective and mechanical turbulence contributions. It accounts for neutral as well as height varying potential temperature gradients above the mixed layer. For thermally stable and mechanically dominated unstable night time Atmospheric Boundary Layer (ABL) the module uses similarity formulae based on the wind velocity [1]), the Monin—Obukhov length [2], and the Coriolis parameter. In the convective case simple slab model is integrated, based on initial lapse rate and the surface heat flux. The lapse rate is evaluated on the basis of vertical atmospheric stability, surface type and surface temperature. This differentiates the IMG model from other existing mixing height models that need initial measured lapse rate for calculation. IMG model is site specific as it calculates the radiative incoming heat flux depending on the solar declination estimates based on-site latitude and longitude. The IMG model is applied to calculate mixing height for India by using surface data (viz. wind speed, surface temperature, surface type) from 152 surface meteorological stations. Results have been evaluated with radiosonde mixing height data procured from 18 upper air stations. Sensitivity analysis of the model with respect to various parameters is performed. The model is formulated after reviewing presently available radiosonde mixing height data in India and can satisfactorily provide an alternative means of estimating mixing height for air pollution dispersion models.