Feedbacks in the Land-Surface and Mixed-Layer Energy Budgets

A mixed-layer model of the surface energy budget and the planetary boundary layer (PBL) is developed, based on the prognostic equations for soil temperature, mixed layer potential temperature and specific humidity and the growth and abrupt collapse of the PBL. Detailed parameterizations of the longwave radiative fluxes are included. The feedbacks in the uncoupled (i.e. surface energy budget with non-responding PBL) and coupled land surface and atmospheric mixed-layer energy budgets are examined. A simplified, time continuous, version of the model, in which the specific humidity budget is the balance of evapotranspiration and dry-air entrainment, and the PBL height is given by the lifted condensation level, is shown to be in good agreement with the complete model. By forcing the simplified model with daily mean rather than periodic solar radiation, an equilibrium model state is achieved where the fluxes are in close agreement with the daily mean fluxes corresponding to the periodic forcing. The model also agrees favorably with measurements from the FIFE field experiment. Feedbacks are examined using the equilibrium model state. The uncoupled and coupled model sensitivities with respect to the minimal stomatal resistance and the atmospheric specific humidity not only differ in magnitude, but in sign as well. This results puts into question the extent to which uncoupled land-surface models that are forced with atmospheric variables may be used in sensitivity studies.     

Estimation of Surface-Layer Scaling Parameters in the Unstable Boundary Layer: Implications for Orographic Flow Speed-Up

A method is proposed for estimating the surface-layer depth \((z_s)\) and the friction velocity \((u_*)\) as a function of stability (here quantified by the Obukhov length, L) over the complete range of unstable flow regimes. This method extends that developed previously for stable conditions by Argaín et al. (Boundary-Layer Meteorol 130:15–28, 2009), but uses a qualitatively different approach. The method is specifically used to calculate the fractional speed-up \((\varDelta S)\) in flow over a ridge, although it is suitable for more general boundary-layer applications. The behaviour of \(z_s \left( L\right) \) and \(u_*\left( L\right) \) as a function of L is indirectly assessed via calculation of \(\varDelta S\left( L\right) \) using the linear model of Hunt et al. (Q J R Meteorol Soc 29:16–26, 1988) and its comparison with the field measurements reported in Coppin et al. (Boundary-Layer Meteorol 69:173–199, 1994) and with numerical simulations carried out using a non-linear numerical model, FLEX. The behaviour of \(\varDelta S\) estimated from the linear model is clearly improved when \(u_*\) is calculated using the method proposed here, confirming the importance of accounting for the dependences of \(z_s\left( L \right) \) and \(u_*\left( L \right) \) on L to better represent processes in the unstable boundary layer.     

Characteristics of carbonyl compounds in ambient air of Shanghai,China

The levels of carbonyl compounds in Shanghai ambient air were measured in five periods from January 2007 to October 2007 (covering winter, high-air-pollution days, spring, summer and autumn). A total of 114 samples were collected and eighteen carbonyls were identified. Formaldehyde, acetaldehyde and acetone were the most abundant carbonyls and their mean concentrations of 19.40 ± 12.00, 15.92 ± 12.07 and 11.86 ± 7.04 μg m⁻³ respectively, in the daytime for five sampling periods. Formaldehyde and acetaldehyde showed similar diurnal profiles with peak mixing ratios in the morning and early afternoon during the daytime. Their mean concentrations were highest in summer and lowest in winter. Acetone showed reversed seasonal variation. The high molecular weight (HMW, ≥C5) carbonyls also showed obvious diurnal variations with higher concentrations in the daytime in summer and autumn, while they were all not detected in winter. Formaldehyde and acetaldehyde played an important role in removing OH radicals in the atmosphere, but the contribution of acetone was below 1%. The carbonyls levels in high-air-pollution days were reported. More carbonyl species with higher concentrations were found in high-air-pollution days than in spring. These carbonyls were transported with other pollutants from north and northwest in March 27 to April 2, 2007 and then mixed with local sources. Comparing with Beijing and Guangzhou, the concentrations of formaldehyde and acetaldehyde in Shanghai were the highest, which indicated that the air pollution in Shanghai was even worse than expected.     

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.     

High-resolution climate simulation of the last glacial maximum

The climate of the last glacial maximum (LGM) is simulated with a high-resolution atmospheric general circulation model, the NCAR CCM3 at spectral truncation of T170, corresponding to a grid cell size of roughly 75 km. The purpose of the study is to assess whether there are significant benefits from the higher resolution simulation compared to the lower resolution simulation associated with the role of topography. The LGM simulations were forced with modified CLIMAP sea ice distribution and sea surface temperatures (SST) reduced by 1°C, ice sheet topography, reduced CO₂, and 21,000 BP orbital parameters. The high-resolution model captures modern climate reasonably well, in particular the distribution of heavy precipitation in the tropical Pacific. For the ice age case, surface temperature simulated by the high-resolution model agrees better with those of proxy estimates than does the low-resolution model. Despite the fact that tropical SSTs were only 2.1°C less than the control run, there are many lowland tropical land areas 4–6°C colder than present. Comparison of T170 model results with the best constrained proxy temperature estimates (noble gas concentrations in groundwater) now yield no significant differences between model and observations. There are also significant upland temperature changes in the best resolved tropical mountain belt (the Andes). We provisionally attribute this result in part as resulting from decreased lateral mixing between ocean and land in a model with more model grid cells. A longstanding model-data discrepancy therefore appears to be resolved without invoking any unusual model physics. The response of the Asian summer monsoon can also be more clearly linked to local geography in the high-resolution model than in the low-resolution model; this distinction should enable more confident validation of climate proxy data with the high-resolution model. Elsewhere, an inferred salinity increase in the subtropical North Atlantic may have significant implications for ocean circulation changes during the LGM. A large part of the Amazon and Congo Basins are simulated to be substantially drier in the ice age—consistent with many (but not all) paleo data. These results suggest that there are considerable benefits derived from high-resolution model regarding regional climate responses, and that observationalists can now compare their results with models that resolve geography at a resolution comparable to that which the proxy data represent.     

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.     

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.     

Multi-Source Emission Determination Using an Inverse-Dispersion Technique

Inverse-dispersion calculations can be used to infer atmospheric emission rates through a combination of downwind gas concentrations and dispersion model predictions. With multiple concentration sensors downwind of a compound source (whose component positions are known) it is possible to calculate the component emissions. With this in mind, a field experiment was conducted to examine the feasibility of such multi-source inferences, using four synthetic area sources and eight concentration sensors arranged in different configurations. Multi-source problems tend to be mathematically ill-conditioned, as expressed by the condition number κ. In our most successful configuration (average κ = 4.2) the total emissions from all sources were deduced to within 10% on average, while component emissions were deduced to within 50%. In our least successful configuration (average κ = 91) the total emissions were calculated to within only 50%, and component calculations were highly inaccurate. Our study indicates that the most accurate multi-source inferences will occur if each sensor is influenced by only a single source. A “progressive” layout is the next best: one sensor is positioned to “see” only one source, the next sensor is placed to see the first source and another, a third sensor is placed to see the previous two plus a third, and so on. When it is not possible to isolate any sources κ is large and the accuracy of a multi-source inference is doubtful.     

Winter precipitation variability over Korean Peninsula associated with ENSO

In this study, winter precipitation variability associated with the El Niño-Southern Oscillation (ENSO) over the Korean Peninsula was investigated using a 5-pentad running mean data because significant correlation pattern cannot be revealed using seasonal-mean data. It was found a considerably significant positive correlation between Niño3 sea-surface temperature and precipitation during early winter (from Mid-November to early-December), when the correlation coefficient is close to 0.8 in early-December; the correlation is distinctively weakened during late winter. It is demonstrated that such sudden intraseasonal change in relation to ENSO is associated with the presence of anticyclonic flow over the Kuroshio extension region (Kuroshio anticyclone). In early winter, there is strong southerly wind over the Korean Peninsula, which is induced by the Philippine Sea anticyclone and Kuroshio anticyclone. However, in January, although the Philippine Sea anticyclone develops further, the Kuroshio anticyclone suddenly disappears; as a result, the impact of ENSO is considerably weakened over the Korean Peninsula. These results indicate that the Kuroshio anticyclone during El Niño peak phase plays a critical role by strongly affecting Northeast Asia climate, including the Korean Peninsula. In addition, it is also found that there are distinctive interdecadal changes of the relationship between ENSO and precipitation over the Korean Peninsula. In particular, the strong correlation in early winter is clearer in the recent 30 years than that in the previous period of 1950–1979.     

How effective is the new generation of GPM satellite precipitation in characterizing the rainfall variability over Malaysia?

We investigated the potential of the new generation of satellite precipitation product from the Global Precipitation Mission (GPM) to characterize the rainfall in Malaysia. Most satellite precipitation products have limited ability to precisely characterize the high dynamic rainfall variation that occurred at both time and scale in this humid tropical region due to the coarse grid size to meet the physical condition of the smaller land size, sub-continent and islands. Prior to the status quo, an improved satellite precipitation was required to accurately measure the rainfall and its distribution. Subsequently, the newly released of GPM precipitation product at half-hourly and 0.1° resolution served an opportunity to anticipate the aforementioned conflict. Nevertheless, related evidence was not found and therefore, this study made an initiative to fill the gap. A total of 843 rain gauges over east (Borneo) and west Malaysia (Peninsular) were used to evaluate the rainfall the GPM rainfall data. The assessment covered all critical rainy seasons which associated with Asian Monsoon including northeast (Nov. - Feb.), southwest (May - Aug.) and their subsequent inter-monsoon period (Mar. - Apr. & Sep. - Oct.). The ability of GPM to provide quantitative rainfall estimates and qualitative spatial rainfall patterns were analysed. Our results showed that the GPM had good capacity to depict the spatial rainfall patterns in less heterogeneous rainfall patterns (Spearman’s correlation, 0.591 to 0.891) compared to the clustered one (r = 0.368 to 0.721). Rainfall intensity and spatial heterogeneity that is largely driven by seasonal monsoon has significant influence on GPM ability to resolve local rainfall patterns. In quantitative rainfall estimation, large errors can be primarily associated with the rainfall intensity increment. 77% of the error variation can be explained through rainfall intensity particularly the high intensity (> 35 mm d^-1). A strong relationship between GPM rainfall and error was found from heavy (~35 mm d^-1) to violent rain (160 mm d^-1). The output of this study provides reference regarding the performance of GPM data for respective hydrology studies in this region.