Adaptation of Pressure Based CFD Solvers for Mesoscale Atmospheric Problems

General purpose Computational Fluid Dynamics (CFD) solvers are frequently used in small-scale urban pollution dispersion simulations without a large extent of ver- tical flow. Vertical flow, however, plays an important role in the formation of local breezes, such as urban heat island induced breezes that have great significance in the ventilation of large cities. The effects of atmospheric stratification, anelasticity and Coriolis force must be taken into account in such simulations. We introduce a general method for adapting pressure based CFD solvers to atmospheric flow simulations in order to take advantage of their high flexibility in geometrical modelling and meshing. Compressibility and thermal stratification effects are taken into account by utilizing a novel system of transformations of the field variables and by adding consequential source terms to the model equations of incompressible flow. Phenomena involving mesoscale to microscale coupled effects can be analyzed without model nesting, applying only local grid refinement of an arbitrary level. Elements of the method are validated against an analytical solution, results of a reference calculation, and a laboratory scale urban heat island circulation experiment. The new approach can be applied with benefits to several areas of application. Inclusion of the moisture transport phenomena and the surface energy balance are important further steps towards the practical application of the method.     

Vertical and Horizontal Transport of Energy and Matter by Coherent Motions in a Tall Spruce Canopy

In the framework of the EGER (ExchanGE processes in mountainous Regions) project, the contribution of coherent structures to vertical and horizontal transports in a tall spruce canopy is investigated. The combination of measurements done in both the vertical and horizontal directions allows us to investigate coherent structures, their temporal scales, their role in flux transport, vertical coupling between the sub-canopy, canopy and air above the canopy, and horizontal coupling in the sub-canopy layer. The temporal scales of coherent structures detected with the horizontally distributed systems in the sub-canopy layer are larger than the temporal scales of coherent structures detected with the vertically distributed systems. The flux contribution of coherent structures to the momentum and sensible heat transport is found to be dominant in the canopy layer. Carbon dioxide and latent heat transport by coherent structures increase with height and reach a maximum at the canopy height. The flux contribution of the ejection decreases with increasing height and becomes dominant above the canopy level. The flux fraction transported during the sweep increases with height and becomes the dominant exchange process at the upper canopy level. The determined exchange regimes indicate consistent decoupling between the sub-canopy, canopy and air above the canopy during evening, nighttime and morning hours, whereas the coupled states and coupled by sweep states between layers are observed mostly during the daytime. Furthermore, the horizontal transport of sensible heat by coherent structures is investigated, and the heterogeneity of the contribution of coherent events to the flux transport is demonstrated. A scheme to determine the horizontal coupling by coherent structures in the sub-canopy layer is proposed, and it is shown that the sub-canopy layer is horizontally coupled mainly in the wind direction. The vertical coupling in most cases is observed together with streamwise horizontal coupling, whereas the cross-stream direction is decoupled.     

Interdecadal variations of cold air activities in Northeast China during springtime

Based on the daily mean temperature data of CN05.2 from 1961 to 2012, cold events (CEs) are first divided into two categories according to their duration: strong cold events (SCEs) and weak cold events (WCEs). Then, the characteristics of CEs, SCEs, and WCEs during springtime are investigated. The results indicate that in the pre-1990s epoch, ENSO and Arctic Oscillation events in the previous winter are closely related to SCEs in the following spring. The multidecadal variations of CEs, SCEs, and WCEs are obvious. The intensity trend for SCEs is significantly negative, but it seems less apparent for WCEs. Further analysis reveals that when both SCEs and WCEs occur, a typical East Asian trough in the 850- hPa wind field, whose northwesterly wind component invades Northeast China (NEC) and causes freezing days, can be found in every decade. For the SCEs, a cold vortex, with its center located over Okhotsk and northeasterly current affecting NEC, is found as an additional feature. For the WCEs, the cold vortex is located in Karafuto and its northwesterly airflow intrudes into NEC. As for the difference between SCEs and WCEs, the northwestern flow is weaker while the northeastern counterpart is stronger during the SCEs, in all decades. In the Takaya–Nakamura flux and divergence fields, for the SCEs, a divergence center exists over NEC; and over its downstream regions, a stronger divergence center appears, not like a wave train. However, the opposite is the case for the WCEs; moreover, the wave train appears clearly during the WCEs, which means that the wave energy can propagate and dissipate more easily during WCEs.     

Influence of coastal wind on surface ozone and nitrogen oxides in suburban Shanghai

Surface ozone, NO, NO₂, and NO _x were measured at a coastal site (Shihua) and a nearby inland site (Zhujing) in suburban Shanghai for the whole year of 2009. More days with ozone pollution in a longer time range were observed at the coastal site than the inland site. The diurnal variations of NO _x concentrations were obviously higher at Zhujing station, while those of ozone concentrations were higher at Shihua station, indicating their different air pollution conditions. Coastal wind has significant influence on the levels and characteristics of the air pollutants. The ozone concentrations during maritime winds (MW) were much higher than those during continental winds (CW) at each of the site, while the NO and NO₂ concentrations were both opposite. The ozone concentrations at Shihua station were much higher than those at Zhujing station, while the NO and NO₂ concentrations were both opposite. The ozone concentrations at both of the two sites showed a distinct “weekend effects” and “weekdays effects” patterns during CW and MW, respectively. Correlation analysis of the pollutants showed that, the compounds during MW were more age than those during CW, and the compounds at Shihua were more age than those at Zhujing. The air pollutions at both of the two sites are mainly associated with the pollutants emitted in this region instead of long range transport.     

Vertical distributions of aerosol optical properties during haze and floating dust weather in Shanghai

A comparative study on the vertical distributions of aerosol optical properties during haze and floating dust weather in Shanghai was conducted based on the data obtained from a micro pulse lidar.There was a distinct difference in layer thickness and extinction coefficient under the two types of weather conditions.Aerosols were concentrated below 1 km and the aerosol extinction coefficients ranged from 0.25 to 1.50km^-1 on haze days.In contrast,aerosols with smaller extinction coefficients(0.20 0.35 km^-1) accumulated mainly from the surface to 2 km on floating dust days.The seasonal variations of extinction and aerosol optical depth（AOD） for both haze and floating dust cases were similar greatest in winter,smaller in spring,and smallest in autumn.More than 85%of the aerosols appeared in the atmosphere below 1 km during severe haze and floating dust weather.The diurnal variation of the extinction coefficient of haze exhibited a bimodal shape with two peaks in the morning or at noon,and at nightfall,respectively.The aerosol extinction coefficient gradually increased throughout the day during floating dust weather.Case studies showed that haze aerosols were generated from the surface and then lifted up,but floating dust aerosols were transported vertically from higher altitude to the surface.The AOD during floating dust weather was higher than that during haze.The boundary layer was more stable during haze than during floating dust weather.     

Are Urban-Canopy Velocity Profiles Exponential?

Using analyses of data from extant direct numerical simulations and large-eddy simulations of boundary-layer and channel flows over and within urban-type canopies, sectional drag forces, Reynolds and dispersive shear stresses are examined for a range of roughness densities. Using the spatially-averaged mean velocity profiles these quantities allow deduction of the canopy mixing length and sectional drag coefficient. It is shown that the common assumptions about the behaviour of these quantities, needed to produce an analytical model for the canopy velocity profile, are usually invalid, in contrast to what is found in typical vegetative (e.g. forest) canopies. The consequence is that an exponential shape of the spatially-averaged mean velocity profile within the canopy cannot normally be expected, as indeed the data demonstrate. Nonetheless, recent canopy models that allow prediction of the roughness length appropriate for the inertial layer’s logarithmic profile above the canopy do not seem to depend crucially on their (invalid) assumption of an exponential profile within the canopy.     

Climatic regime shift and decadal anomalous events in China

Climatic time series from historical documents and instrumental records from China showed temporal and regional patterns in the last two to three centuries, including two multidecadal oscillations at quasi-20-year and quasi-70-year timescales revealed by signal analysis from wavelet transform. Climatic anomalous events on the decadal timescale were identified based on the two oscillations when their positive (or negative) phases coincide with each other to amplify amplitude. The coldest event occurred in the decade of 1965–1975 in eastern China, while the periods of 1920–1930, 1940–1950, and 1988–2000 appeared to be warmer in most parts of China. For the precipitation series in northern China, the dry anomalous event was found in the late 1920s, while the wet anomalous event occurred in the 1950s. A severe drought in 1927–1929 in northern China coincided with the anomalous warm and dry decade, caused large-scale famine in nine provinces over northern China. Climatic anomalous events with a warm-dry or cold-wet association in the physical climate system would potentially cause severe negative impacts on natural ecosystem in the key vulnerable region over northern China. The spatial pattern of summer rainfall anomalies in the eastern China monsoon region showed an opposite variations in phase between the Yellow River Valley (North China) and the mid-low Yangtze River Valley as well as accompanied the shift of the northernmost monsoon boundary. Climatic regime shifts for different time points in the last 200 years were identified. In North China, transitions from dry to wet periods occurred around 1800, 1875, and 1940 while the transitions from wet to dry periods appeared around 1840, 1910, and the late 1970s. The reversal transition in these time points can also be found in the lower Yangtze River. Climatic regime shifts in China were linked to the interaction of mid- and low latitude atmospheric circulations (the westerly flow and the monsoon flow) when they cross the Tibetan Plateau in East Asia.     

Scaling of the Vertical Spreading of a Plume of a Passive Tracer in a Rough-Wall Neutral Boundary Layer

The influence of surface roughness on the dispersion of a passive scalar in a rough wall turbulent boundary layer has been studied using wind-tunnel experiments. The surface roughness was varied using different sizes of roughness elements, and different spacings between the elements. Vertical profiles of average concentration were measured at different distances downwind of the source, and the vertical spread of the plume was computed by fitting a double Gaussian profile to the data. An estimate of the integral length scale is derived from the turbulence characteristics of the boundary layer and is then used to scale the measured values of plume spread. This scaling reduces the variability in the data, confirming the validity of the model for the Lagrangian integral time scale, but does not remove it entirely. The scaled plume spreading shows significant differences from predictions of theoretical models both in the near and in the far field. In the region immediately downwind of the source this is due to the influence of the wake of the injector for which we have developed a simple model. In the far field we explain that the differences are mainly due to the absence of large-scale motions. Finally, further downwind of the source the scaled values of plume spread fall into two distinct groups. It is suggested that the difference between the two groups may be related to the lack of dynamical similarity between the boundary-layer flows for varying surface roughness or to biased estimates of the plume spread.     

The influence of tropical Pacific forcing on the Arctic Oscillation

The potential role of tropical Pacific forcing in driving the seasonal variability of the Arctic Oscillation (AO) is explored using both observational data and a simple general circulation model (SGCM). A lead–lag regression technique is first applied to the monthly averaged sea surface temperature (SST) and the AO index. The AO maximum is found to be related to a negative SST anomaly over the tropical Pacific three months earlier. A singular value decomposition (SVD) analysis is then performed on the tropical Pacific SST and the sea level pressure (SLP) over the Northern Hemisphere. An AO-like atmospheric pattern and its associated SST appear as the second pair of SVD modes. Ensemble integrations are carried out with the SGCM to test the atmospheric response to different tropical Pacific forcings. The atmospheric response to the linear fit of the model’s empirical forcing associated with the SST variability in the second SVD modes strongly projects onto the AO. Idealized thermal forcings are then designed based on the regression of the seasonally averaged tropical Pacific precipitation against the AO index. Results indicate that forcing anomalies over the western tropical Pacific are more effective in generating an AO-like response while those over the eastern tropical Pacific tend to produce a Pacific-North American (PNA)-like response. The physical mechanisms responsible for the energy transport from the tropical Pacific to the extratropical North Atlantic are investigated using wave activity flux and vorticity forcing formalisms. The energy from the western tropical Pacific forcing tends to propagate zonally to the North Atlantic because of the jet stream waveguide effect while the transport of the energy from the eastern tropical Pacific forcing mostly concentrates over the PNA area. The linearized SGCM results show that nonlinear processes are involved in the generation of the forced AO-like pattern.     

Atmospheric chemistry of hydrogen fluoride

Although a large volume of monitoring and computer simulation data exist for global coverage of HF, study of HF in the troposphere is still limited to industry whose primary interest is the safety and risk assessment of HF release because it is a toxic gas. There is very limited information on atmospheric chemistry, emission sources, and the behavior of HF in the environment. We provide a comprehensive review on the atmospheric chemistry of HF, modeling the reactions and transport of HF in the atmosphere, the removal processes in the vertical layer immediately adjacent to the surface (up to approximately 500 m) and recommend research needed to improve our understanding of atmospheric chemistry of HF in the troposphere. The atmospheric chemistry, emissions, and surface boundary layer transport of hydrogen fluoride (HF) are summarized. Although HF is known to be chemically reactive and highly soluble, both factors affect transport and removal in the atmosphere, the chemistry can be ignored when the HF concentration is at a sufficiently low level (e.g., 10 ppmv). At a low concentration, the capability for HF to react in the atmosphere is diminished and therefore the species can be mathematically treated as inert during the transport. At a sufficiently high concentration of HF (e.g., kg/s release rate and thousands of ppm), however, HF can go through a series of rigorous chemical reactions including polymerization, depolymerization, and reaction with water to form molecular complex. As such, the HF species cannot be considered as inert because the reactions could intimately influence the plume’s thermodynamic properties affecting the changes in plume temperature and density. The atmospheric residence time of HF was found to be less than four (4) days, and deposition (i.e., atmosphere to surface transport) is the dominant mechanism that controls the removal of HF and its oligomers from the atmosphere. The literature data on HF dry deposition velocity was relatively high compared to many commonly found atmospheric species such as ozone, sulfur dioxide, nitrogen oxides, etc. The global average of wet deposition velocity of HF was found to be zero based on one literature source. Uptake of HF by rain drops is limited by the acidity of the rain drops, and atmospheric particulate matter contributes negligibly to HF uptake. Finally, given that the reactivity of HF at a high release rate and elevated mole concentration cannot be ignored, it is important to incorporate the reaction chemistry in the near-field dispersion close to the proximity of the release source, and to incorporate the deposition mechanism in the far-field dispersion away from the release source. In other words, a hybrid computational scheme may be needed to address transport and atmospheric chemistry of HF in a range of applications. The model uncertainty will be limited by the precision of boundary layer parameterization and ability to accurately model the atmospheric turbulence.