Real Data Assimilation Using the Local Ensemble Transform Kalman Filter (LETKF) System for a Global Non-hydrostatic NWP model on the Cubed-sphere

An ensemble data assimilation system using the 4-dimensional Local Ensemble Transform Kalman Filter is implemented to a global non-hydrostatic Numerical Weather Prediction model on the cubed-sphere. The ensemble data assimilation system is coupled to the Korea Institute of Atmospheric Prediction Systems Package for Observation Processing, for real observation data from diverse resources, including satellites. For computational efficiency in a parallel computing environment, we employ some advanced software engineering techniques in the handling of a large number of files. The ensemble data assimilation system is tested in a semi-operational mode, and its performance is verified using the Integrated Forecast System analysis from the European Centre for Medium-Range Weather Forecasts. It is found that the system can be stabilized effectively by additive inflation to account for sampling errors, especially when radiance satellite data are additionally used.     

A dynamic link between the basin-scale and zonal modes in the Tropical Indian Ocean

Summary The interannual variability of sea surface temperature (SST) anomalies in the tropical Indian Ocean is dominated mainly by a basin-scale mode (BM) and partly by an east–west contrast mode (zonal mode, ZM). The BM reflects the basin-scale warming or cooling and is highly correlated with El Nino with 3- to 6-month lags, while the ZM is marginally correlated with El Nino with 9-month lags.During an El Nino, large-scale anomalous subsidence over the maritime continent occurs as a result of an eastward shift in the rising branch of the Walker circulation suppresses convection over the eastern Indian Ocean, allowing more solar radiation over the eastern Indian Ocean. At the same time, the anomalous southeasterly wind over the equatorial Indian Ocean forces the thermocline over the western Indian Ocean to deepen, especially in the southern part. As a result, SST over the whole basin increases. As El Nino decays, the subsidence over the maritime continent ceases and so does the anomalous southeasterly wind. However, the thermocline perturbation does not quickly shoal back to normal because of inertia and it disperses as Rossby waves. These Rossby waves are reflected back as an equatorial Kelvin wave, causing deepening of the thermocline in the eastern Indian Ocean, and preventing SSTs from cooling in that region. Moreover, the weaker wind speed of the monsoon circulation results in less latent heat loss, and thus warms the eastern Indian Ocean. These two processes therefore help to maintain warm SSTs over the eastern Indian Ocean until fall. During the fall, the warm SST over the eastern Indian Ocean and the cold SST over the western Indian Ocean are enhanced by air–sea interaction and the ZM returns. The ZM dissipates through the seasonal reversal of the monsoon atmospheric circulation and the boundary-reflected Kelvin wave. In the same manner, a basin-scale cooling in the tropical Indian Ocean can induce the ZM warming in the west and cooling in the east.     

Arctic Climate Changes Based on Historical Simulations (1900‒2013) with the CAMS-CSM

The Chinese Academy of Meteorological Sciences Climate System Model (CAMS-CSM) is a newly developed global climate model that will participate in the Coupled Model Intercomparison Project phase 6. Based on historical simulations (1900?2013), we evaluate the model performance in simulating the observed characteristics of the Arctic climate system, which includes air temperature, precipitation, the Arctic Oscillation (AO), ocean temperature/salinity, the Atlantic meridional overturning circulation (AMOC), snow cover, and sea ice. The model?data comparisons indicate that the CAMS-CSM reproduces spatial patterns of climatological mean air temperature over the Arctic (60°?90°N) and a rapid warming trend from 1979 to 2013. However, the warming trend is overestimated south of the Arctic Circle, implying a subdued Arctic amplification. The distribution of climatological precipitation in the Arctic is broadly captured in the model, whereas it shows limited skills in depicting the overall increasing trend. The AO can be reproduced by the CAMS-CSM in terms of reasonable patterns and variability. Regarding the ocean simulation, the model underestimates the AMOC and zonally averaged ocean temperatures and salinity above a depth of 500 m, and it fails to reproduce the observed increasing trend in the upper ocean heat content in the Arctic. The large-scale distribution of the snow cover extent (SCE) in the Northern Hemisphere and the overall decreasing trend in the spring SCE are captured by the CAMS-CSM, while the biased magnitudes exist. Due to the underestimation of the AMOC and the poor quantification of air–sea interaction, the CAMS-CSM overestimates regional sea ice and underestimates the observed decreasing trend in Arctic sea–ice area in September. Overall, the CAMS-CSM reproduces a climatological distribution of the Arctic climate system and general trends from 1979 to 2013 compared with the observations, but it shows limited skills in modeling local trends and interannual variability.     

UV attenuation in the cloudy atmosphere

Ultraviolet (UV) energy absorption plays a very important role in the Earth–atmosphere system. Based on observational data for Beijing, we suggest that some atmospheric constituents utilize or transfer UV energy in chemical and photochemical (C&P) reactions, in addition to those which absorb UV energy directly. These constituents are primarily volatile organic compounds (VOCs) emitted from both vegetative and anthropogenic sources. The total UV energy loss in the cloudy atmosphere for Beijing in 1990 was 78.9 Wm⁻². This attenuation was caused by ozone (48.3 Wm⁻²), other compounds in the atmosphere (26.6 Wm⁻²) and a scattering factor (4.0 Wm⁻²). Our results for a cloudy atmosphere in the Beijing area show that the absorption due to these other compounds occurs largely through the mediation of water vapor. This fraction of energy loss has not been fully accounted for in previous models. Observations and previous models results suggest that 1) a cloudy atmosphere absorbs 25∼30 Wm⁻² more solar shortwave radiation than models predict; and 2) aerosols can significantly decrease the downward mean UV-visible radiation and the absorbed solar radiation at the surface by up to 28 and 23 Wm⁻², respectively. Thus, quantitative study of UV and visible absorption by atmospheric constituents involved in homogeneous and heterogeneous C&P reactions is important for atmospheric models.     

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.     

Interdecadal change in the connection between Hadley circulation and winter temperature in East Asia

Based on NCEP/NCAR reanalysis data, the interdecadal variability of Hadley circulation （HC） and its association with East Asian temperature in winter are investigated. Results indicate that the Northern Hemisphere winter HC underwent apparent change in the 1970s, with transition occurring around 1976/77. Along with interdecadal variability of HC, its linkage to surface air temperature （SAT） in East Asia also varied decadally, from weak relations to strong relations. Such a change may be related to the interaction between HC and the atmospheric circulation system over the Philippines, which is associated with the East Asian winter monsoon （EAWM）. Before the 1970s, the connection between HC and the anticyclonic circulation around the Philippines was insignificant, but after the late 1970s their linkage entered a strong regime. The intensification of this connection may therefore be responsible for the strong relations between HC and East Asian winter temperatures after the late 1970s.     

Implementation of a surface runoff model with Horton and Dunne mechanisms into the regional climate model RegCM_NCC

A surface runoff parameterization scheme that dynamically represents both Horton and Dunne runoff generation mechanisms within a model grid cell together with a consideration of the subgrid-scaie soil heterogeneity, is implemented into the National Climate Center regional climate model （RegCM_NCC）. The effects of the modified surface runoff scheme on RegCMANCC performance are tested with an abnormal heavy rainfall process which occurred in summer 1998. Simulated results show that the model with the original surface runoff scheme （noted as CTL） basically captures the spatial pattern of precipitation, circulation and land surface variables, but generally overestimates rainfall compared to observations. The model with the new surface runoff scheme （noted as NRM） reasonably reproduces the distribution pattern of various variables and effectively diminishes the excessive precipitation in the CTL. The processes involved in the improvement of NRM-simulated rainfall may be as follows： with the new surface runoff scheme, simulated surface runoff is larger, soil moisture and evaporation （latent heat flux） are decreased, the available water into the atmosphere is decreased; correspondingly, the atmosphere is drier and rainfall is decreased through various processes. Therefore, the implementation of the new runoff scheme into the RegCMANCC has a significant effect on results at not only the land surface, but also the overlying atmosphere.     

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.