Numerical simulation of January 28 cold air outbreak during GALE Part I: The model and sensitivity tests of turbulence closures

A mesoscale planetary boundary layer (PBL) numerical model has been developed to study airflow over complex topography. Turbulence closures using the turbulent kinetic energy (TKE) and dissipation () equations are investigated in combination with the level 2.5 scheme of Mellor and Yamada (1982) to determine eddy diffusivities for momentum and heat. This modified E- closure is simpler than the level 3 one which requires more prognostic equations for moist turbulent transport.One-dimensional (1-D) model results show that the PBL mean flows under various stability conditions are not significantly sensitive to the modified Blackadar and Kolmogorov eddy mixing-length formulations used in this E- model, although the latter yields excessively large mixing lengths in the entrainment region of the upper PBL. Eddy mixing lengths in the Kolmogorov-type formulation can be better defined by introducing background dissipation. Using the same prognostic TKE equation, the 1-D model results are not significantly affected by different diagnostic formulations in the closures. The simulated results compare well with large-eddy simulations and those obtained using higher-order closure schemes including the level 3 one. The results are found to be insensitive to eddy Prandtl number, in contrast to the 2-D model results (see Part II).     

Urban effects of Chennai on sea breeze induced convection and precipitation

Doppler radar derived wind speed and direction profiles showed a well developed sea breeze circulation over the Chennai, India region on 28 June, 2003. Rainfall totals in excess of 100 mm resulted from convection along the sea breeze front. Inland propagation of the sea breeze front was observed in radar reflectivity imagery. High-resolution MM5 simulations were used to investigate the influence of Chennai urban land use on sea breeze initiated convection and precipitation. A comparison of observed and simulated 10m wind speed and direction over Chennai showed that the model was able to simulate the timing and strength of the sea breeze. Urban effects are shown to increase the near surface air temperature over Chennai by 3.0K during the early morning hours. The larger surface temperature gradient along the coast due to urban effects increased onshore flow by 4.0m s⁻¹. Model sensitivity study revealed that precipitation totals were enhanced by 25mm over a large region 150 km west of Chennai due to urban effects. Deficiency in model physics related to night-time forecasts are addressed.     

Observational and Numerical Study on the Influence of Large-Scale Flow Direction and Coastline Shape on Sea-Breeze Evolution

In this study radar, surface observations and numericalsimulations are used to examine the inland penetrationand intensity of the sea breeze during various large-scaleflow regimes along the curved coastline of the Carolinas,U.S.A. The results clearly indicate that the flow directionrelative to the curved coastline has a significant effecton the sea-breeze evolution.Overall, during northerly flow regimes alongthe curved North Carolina coast, observationsand numerical simulations show that the sea-breezefront has a tendency to remain close to the south-facingcoast. During these same flow regimes the frontmoves further inland relative to the east-facingcoast. The sea-breeze front during westerly flowcases progressed further inland relative to the southcoast and less so from the east-facing coastline.South-westerly flow allows the sea breeze to moveinland from both coastlines but the coastal shapeinfluence makes the inland penetration less fromthe easterly facing beaches. During periods of lightonshore flow (south-east), the sea breeze movesconsiderable distances inland but is not discernableuntil later in the afternoon. The simulations indicatedthat the sea-breeze intensity is greatest (least) when thelarge-scale flow direction has an offshore (onshore)component. Model results indicate the existence of astrong front well inland in the late afternoon duringlight onshore flow. Also noted was that the simulatedsea-breeze front develops earlier in the afternoon duringoffshore regimes and later in the day as the large-scaleflow becomes more onshore. It is concluded that thecoastline shape and coast-relative flow direction areimportant factors in determining how the sea-breezecirculation evolves spatially.     

Application of the E-ε closure model to simulations of mesoscale topographic effects

A mesoscale Planetary Boundary Layer (PBL) model with a simple turbulence closure scheme based on the turbulence kinetic energy (TKE) equation and the dissipation () equation is used to simulate atmospheric flow over mesoscale topography. Comparative studies with different parameterizations suggest that with a proper closure assumption for turbulence dissipation, the E-model can simulate the circulation induced by the mesoscale topography with results similar to those obtained using the E- model. On the other hand, the first-order closure using O'Brien's cubic interpolation for eddy diffusivities (K) generally produces much larger K profiles in the stable and the unstable regions, which is believed to be due to the overprediction of the height of the PBL. All models with the TKE equation yield quite similar ensemble mean fields, which are found to be little sensitive to the closure assumption for turbulence dissipation, though their predicted magnitudes of TKE and K may differ appreciably. A discussion on the diurnal evolution of the mesoscale topography-induced circulation and the spatial variations of the turbulence fluxes in the surface layer is also given based on the E- model results.     

Structure of the marine atmospheric boundary layer during two cold air outbreaks of varying intensities: GALE 86

Aircraft, surface, upper air and satellite measurements have been used to observe the evolution and growth of the convective Marine Atmospheric Boundary Layer (MABL) offshore of North Carolina in close proximity to the Gulf Stream, during the intense cold air outbreak of 28 January 1986 and the moderate event of 12 February 1986, as part of the Genesis of Atlantic Lows Experiment (GALE). Air mass modification processes, driven primarily by the ocean-atmosphere exchanges of surface turbulent sensible and latent heat fluxes, caused the overlying air mass to warm and moisten as it advected over the warmer waters of the eastern United States continental shelf. Maximum observed near-surface total heat fluxes were 1045 and 811 W·m^–2 over the core of the Gulf Stream, for 28 January and 12 February 1986, respectively. The observed changes in the overlying air mass occurred almost instantaneously as the ambient flow traversed different underlying SST conditions.The turbulent structure showed a buoyancy-dominated MABL below approximately 0.8z/h. However, shear was also observed to be an important production term above 0.8z/h and below 0.1z/h for the 28 January 1986 event. Dissipation of turbulent kinetic energy was the dominant destruction term in the budgets, but vertical transport of energy was a strong contributor below 0.5z/h, above which this term became a source of turbulent energy. Additionally, the normalized standard deviations of the horizontal velocity components showed a near-equal contribution to the turbulence, while the vertical velocity components displayed the characteristic mid-layer maximum profile observed for a convective, well-mixed boundary layer.     

A study of mean boundary-layer structures over the Arabian Sea and the Bay of Bengal during active and break monsoon periods

Observations from research ships which took part in the Indo-Soviet Monsoon Experiment of 1977 (MONSOON 77) and the International Monsoon Experiments (MONEX 79) over the central Arabian Sea and the north central Bay of Bengal were analyzed to study the mean wind and temperature structure of the monsoon boundary layer during active and break conditions. Mean profiles of wind speed and direction along with virtual potential temperature obtained by averaging data from several research ships during 1977 and 1979 indicate that onset conditions were associated with substantial increases in wind speed over the Arabian Sea and a shift to strong southwest flow. Monsoon onset was also characterized by near-neutral to slightly unstable temperature profiles in the lowest kilometer. Break conditions in 1977 in which the monsoon trough moved northward and substantial (5 mb) pressure rises were noted over the Arabian Sea show wind speeds typically decreasing from approximately 18 m s^–1 during active conditions to roughly 8 m s ^–1. Temperature profiles during break conditions are similar to those observed in pre-monsoon conditions in that the boundary layer is observed to be generally much more stable up to 900 mb. Above 900 mb, profiles of virtual potential temperature show little variation.Analysis of latent and sensible heat fluxes during June 1977 calculated by the bulk aerodynamic method indicates values of latent heat flux during active conditions to be roughly two to three times larger than those during break conditions. Sensible heat flux shows an increase from approximately 20 to 80 W m ^–1 during the onset of the monsoon. Surface fluxes of water vapor indicate the importance of water vapor transport over the ship observation region in the central Arabian Sea during active conditions. Onset of the monsoon over the Arabian Sea is accompanied by an increase in the surface moisture flux by a factor of about two. Time histories of precipitable water show decreases of approximately 15% from active to break periods.     

Surface layer turbulence spectra and dissipation rates during low winds in tropics

Spectral characteristics of surface layer turbulence in an urban atmosphere are investigated. The observations used for this purpose represent low wind conditions in the tropics. The normalized power spectral shapes exhibit the usual characteristics in the inertial subrange and obey Monin-Obukhov scaling. However, the low-frequency behaviours do not conform to the previous observed relations. For horizontal components, large energy is contained in the low frequencies in contrast to the vertical component where roll-off to zero frequency is faster.The turbulent kinetic energy dissipation rate estimated from the spectra using Kolmogorov's inertial subrange law is found to be isotropic unlike the velocity variances. The expressions for the dimensionless dissipation rate do not seem to work well in low winds in an urban atmosphere. For the data considered, the dissipation rate exhibits a power law relationship with the mean windspeed and the friction velocity.