South foehn in the Wipp Valley – Innsbruck region: Numerical simulations of the 24 October 1999 case (MAP-IOP 10)

Summary Numerical simulations of the 24 October 1999 south foehn (MAP-IOP 10) are performed with the Penn State/NCAR mesoscale model MM5 for the Wipp Valley and the adjacent parts of the Inn Valley. The model is run in a multiple-nest configuration, the area of interest being resolved at a mesh size of 800m in most experiments. The study serves to complement an earlier work in which typical flow features of the foehn in the Wipp Valley region were investigated by means of idealized simulations, assessing whether it is possible to reproduce the temporal evolution and the spatial structure of a particular foehn case. A further objective of the paper is to examine the dependence of the model performance on the horizontal resolution, giving some information which resolution will probably be needed for future high-resolution forecasts.An encouragingly large part of the observed flow features could be well reproduced in the simulations. Except for a small region to the east of Innsbruck, the foehn breakthrough is predicted correctly to within an hour. The spatial structure of the so-called pre-foehn, an enhanced westerly wind occurring at Innsbruck prior to the breakthrough of the foehn, also agrees very well with the observations. Moreover, the maximum extent of the foehn in the Inn Valley, the structure of the gravity wave field above the Wipp Valley and the upvalley progression of a shallow cold front in the evening are consistent with the observations. Except for a few places where the airmass boundary between the warm foehn air and the adjacent colder air is not captured correctly throughout the time, the simulated surface temperatures range within 2K of the observed values. Discrepancies between the model results and the observations are found in the vicinity of Innsbruck where a flow-splitting phenomenon induces a very complex flow pattern at low levels. Another source of problems is the surface potential temperature along the Wipp Valley. The observed potential-temperature increase between the Brenner Pass and Innsbruck, which appears to be related to turbulent vertical mixing of stably stratified air, is underestimated by the model. Reducing the horizontal resolution from 800m to 1.4km deteriorates the model performance in marginally resolved side valleys, but the results obtained for the Wipp Valley and the Inn Valley are still of high quality.     

A study on the number of snow days over Eurasia, Indian rainfall and seasonal circulations

Summary In this paper, the relationship between seasonal mean (June, July, August and September) monsoon circulation features and the midlatitude circulations in winter and spring seasons have been examined during contrasting years of more (less) number of snow days in winter/spring followed by deficient (excess) Indian Summer Monsoon Rainfall (ISMR) using NCEP/NCAR reanalyzed data for the period 1966–1994. The Historical Soviet Daily Snow Depth (HSDSD) version II data set has been used to calculate the number of days of snow over west and east Eurasia separately under three classes: class 1 for SD>5cm, class 2 for SD>10cm and class 3 for SD>50cm where SD stands for snow depth. Correlation coefficients are computed between the anomaly in the number of days of snow depth under the above three classes during winter/spring over west and east Eurasia and the subsequent ISMR. HSDSD data show that difference in the number of days of SD>10cm in two extreme years is most prominent in the west Eurasia in the months of January and April. Also the anomaly in the number of days of snow in January and April over west Eurasia has correlation coefficients of –0.69 and –0.56 with the following ISMR, respectively at 0.1% significance level when the SD is more than 10cm at all the stations. Results also show that low-level atmospheric temperature difference between two extreme years of snow days in winter is up to 10°C and the cooling persists up to spring season with a difference of 2°C. This cooling persistence may give rise to anomalous cyclonic circulations over the midlatitudes and tropics which may be responsible for weakening the monsoon circulation over India during the year of more snow days over west Eurasia.     

Dynamical and microphysical characteristics of Arctic clouds using integrated observations collected over SHEBA during the April 1998 FIRE.ACE flights of the Canadian Convair

Summary The purpose of this study is to better understand the dynamical and microphysical processes within Arctic clouds, which occurred in April 1998 over the Surface Heat Budget of the Arctic Ocean (SHEBA) ship during the First ISCCP (International Satellite Cloud Climatology Project) Regional Experiment Arctic Cloud Experiment (FIRE.ACE). The observations from the four cases in the present study were collected by instruments mounted onboard the National Research Council (NRC) Convair, as well as, the National Oceanic and Atmospheric Administration (NOAA) Advanced Very High Resolution Radiometer (AVHRR) satellite, the SHEBA surface based NOAA Doppler radar (35GHz, Ka-Band), and the NOAA depolarization lidar (0.523µm) measurements. The aircraft observations were collected at 32Hz (3-m scale). The Meteorological Services of Canada (MSC) lidar (1.064µm) was operated onboard the Convair-580. The AVHRR observations, representing a 5-km horizontal resolution, were used to estimate particle size, phase, and optical thickness. Constant altitude flight legs were made at about 100m over the ocean surface. Vertical air velocity (w), reflectivity and Doppler velocity, and backscatter and depolarization ratio values were used to define the size of the important dynamical structures. Ice crystal number concentration (N_i), ice water content (IWC), droplet number concentration (N_d), liquid water content (LWC) and characteristic particle size and shape were summarized for each case. The effective radius (r_eff) values for liquid clouds obtained from in-situ and AVHRR observations were found comparable. The large variability in IWC can be due to undetected ice crystals at small size ranges. Mixed phased conditions in the AVHRR retrievals complicated the comparisons with in-situ data. N_i was found to be directly related to the history of the air-parcel dynamics e.g., w. The variability and differences in the parameters obtained from various platforms can be attributed to their instrumental capabilities, resolution, as well as the cloud development.     

Parameterizing turbulent diffusion through the joint probability density

The convective mass flux parameterization often used in meteorological modeling expresses the vertical flux of a transported scalar as proportional to the product of the difference in mean values of the scalar in updrafts and downdrafts and their characteristic velocity. The proportionality factor is a constant to be specified. We show that this proportionality factor also appears in the relaxed eddy accumulation technique of Businger and Oncley. That associates the surface-layer flux of a scalar with the product of the standard deviation of vertical velocity and the mean concentration difference between updrafts and downdrafts.We show that this constant (b) is determined uniquely by the joint probability density (jpd) of vertical velocity and the scalar. Using large-eddy simulation, we generate this jpd for a conservative scalar diffusing through a convective boundary layer. It has quite different forms in top-down and bottom-up diffusion geometries. The bottom-up jpd is fairly well represented by a jointly Gaussian form and implies b ~ 0.6, in good agreement with the surface-layer value reported by Businger and Oncley. The top-down jpd is strikingly non-Gaussian and gives b ~ 0.47. Updrafts carry the bulk of the scalar flux - 70% in the bottom-up case, 60% in the top-down case.The National Center for Atmospheric Research is sponsored by the National Science Foundation.     

Anomalous gradient winds in the subtropical jet stream and interpretations of forecast failures

Summary Dramatic examples of forecast failures in global models of moderate resolution (i.e., T106) have been shown to occur during periods of the negative phase of the Pacific/North American (PNA) pattern in the Northern Hemisphere winter. Specifically, in these periods forecast skills at 500hPa as measured by the standard anomaly correlation index dropped to rather low values by days 4 and 5 of the forecasts. This paper examines systematically some of the factors that may have contributed to the failure of these model forecasts.In particular, strong winds approaching intensities on the order of 100ms^–1 south of Japan at the 200hPa level were degraded by the initialization and data assimilation procedures of the models. These observed winds were found to be supergradient in nature and representative of the anomalous solution of the gradient wind equation. Procedures such as the multivariate optimum interpolation (with its geostrophic constraints) and the normal-mode initialization including several vertical modes apparently were factors that led to the degradation of these strong winds in the initial model states. In this paper, an analysis of these factors is presented, and it is shown that uninitialized analyses (with no constraints) based on a simple successive correction procedure can retain the strong winds evident in the observations. Forecasts thus performed appear to retain wave trains, a characteristic feature of negative PNA initial states, leading to a significant improvement in forecast skill.     

Air temperature retrieval from remote sensing data based on thermodynamics

Summary A new approach to retrieving air temperature from land surface temperature is presented. The new method is based on thermodynamics. Two important parameters, namely crop water stress index and aerodynamic resistance, were used to build a quantitative relationship between the land surface temperature and the ambient air temperature. The method was applied using MODIS satellite data for a location situated in the North China Plain. Comparing the measurement values at meteorological stations with air temperature, derived by the method for certain pixels, indicates that derived values can be obtained within an accuracy of 3°C for more than 80% of data processed. Sensitivity studies also suggest that inaccuracies associated with measurement error in the model variables are also within the 3°C range.     

Infrared observations and numerical modelling of grassland fires in the Northern Territory, Australia

Summary This paper reports on a small-scale pilot experiment held early in the dry season near Darwin, Australia, in which fine-scale observations of several prescribed fires were made using infrared digital video. Infrared imaging is used routinely to locate fires as infrared radiation suffers little attenuation as it propagates through the smoke that normally obscures visible imagery. However, until now, little use has been made of digital video imagery in analyzing the convective-scale structure of prescribed (or wild) fires. The advantage of digital video imagery is that the individual frames can be objectively analyzed to determine the convective motion in the plane viewed by the camera. The infrared imagery shows mostly rising plumes, much like convective clouds. The flow is highly convective, and the vertical transport of heat is confined to relatively narrow thermals. The updrafts range from a few ms^–1 to around 15ms^–1. A numerical model is used to simulate one of the prescribed fires at very high-resolution. For the most part, the model predictions compare well to the observations. The model produces plumes that are around 7m high, and spaced around 5m apart, which is similar to that observed. The model correctly predicts the mean rate of spread of the fire to be 1.3ms^–1. Perhaps the most serious limitations to using infrared observations of the type presented here are the difficulties in interpreting precisely the relationship between the observed infrared temperature field and the air temperature calculated by the model, and the exact connection between the infrared camera derived flow field and that calculated by the model.     

Application of a “Big-Tree” model to regional climate modeling: a sensitivity study

Summary ¶In order to better understand land-atmosphere interactions and increase the predictability of climate models, it is important to investigate the role of forest representation in climate modeling. Corresponding to the big-leaf model commonly employed in land surface schemes to represent the effects of a forest, a so called big-tree model, which uses multi-layer vegetation to represent the vertical canopy heterogeneity, was introduced and incorporated into the National Center for Atmospheric Research (NCAR) regional climate model RegCM2, to make the vegetation model more physically based. Using this augmented RegCM2 and station data for China during 1991 Meiyu season, we performed 10 experiments to investigate the effects of the application of the big-tree model on the summer monsoon climate.With the big-tree model incorporated into the regional climate model, some climate characteristics, e.g. the 3-month-mean surface temperature, circulation, and precipitation, are significantly and systematically changed over the model domain, and the change of the characteristics differs depending on the area. Due to the better representation of the shading effect in the big-tree model, the temperature of the lower layer atmosphere above the plant canopy is increased, which further influences the 850hPa temperature. In addition, there are significant decreases in the mean latent heat fluxes (within 20–30W/m²) in the three areas of the model domain.The application of the big-tree model influences not only the simulated climate of the forested area, but also that of the whole model domain, and its impact is greater on the lower atmosphere than on the upper atmosphere. The simulated rainfall and surface temperature deviate from the originally simulated result and are (or seem to be) closer to the observations, which implies that an appropriate representation of the big-tree model may improve the simulation of the summer monsoon climate.We also find that the simulated climate is sensitive to some big-tree parameter values and schemes, such as the shape, height, zero-plane displacement height and mixing-length scheme. The simulated local/grid differences may be very large although the simulated areal-average differences may be much lower. The area-average differences in the monthly-mean surface temperature and heat fluxes can amount to 0.5°C and 4W/m², respectively, which correspond to maximum local/grid differences of 3.0°C and 40W/m² respectively. It seems that the simulated climate is most sensitive to the parameter of the zero-plane displacement among the parameters studied.     

Changes in extreme daily mean temperatures in summer in eastern China during 1955–2000

Summary Statistical characteristics of extremely low and high daily mean temperatures in summer (June, July and August) in eastern China have been investigated. The extremely low temperatures are defined as those days with temperatures not exceeding the 10th percentile with respect to the reference period of 1961–90; similarly the extremely high temperatures are defined as those exceeding the 90th percentile. There are well-defined spatial structures in trends of the frequency of extremely low temperatures as well as of high temperature extremes. In the north region (i.e. northern and northeastern China) the linear trends of frequency of low and high temperature extremes are –1.09 and +1.23 days/10yr, respectively. For the southern portion of the study area, the trends are –1.32 and –2.32 days/10yr. Taking the study area as a whole, the linear trends are –0.76 days/10yr and +1.08 days/10yr, respectively. The changes of frequency of extreme temperatures are mainly related to the shift in the temperature means. There is a dominant anticyclonic pattern in the lower- to middle troposphere over East Asia in association with warmer conditions in the north region. For the south region there is a jump-like change in the summer mean temperature and the extreme temperature events in around 1976. The large-scale northwestern Pacific subtropical high plays an important role in the jump-like changes of the temperature extremes.     

Land use and local climate: A case study near Santa Cruz, Bolivia

Summary The Tierras Bajas regions of eastern Santa Cruz, Bolivia have undergone among the most rapid rates of concentrated deforestation during the 1980s and 1990s. We investigate the sensitivity of local climate to these land cover changes as observed from Landsat images acquired between 1975 and 1999. The Simple Biosphere model (SiB2) is used to assess the effects of both morphological and physiological changes in vegetation and the implications for fluxes of water, energy and carbon between the vegetation and the atmosphere during the rainy season.Conversion from tropical forest to cropland implicates morphological changes in vegetation as the primary drivers for a daily maximum warming of about 2°C and a slight nighttime cooling, suggesting that clearing of tropical forests for agricultural use may increase the diurnal temperature range, mainly by increasing the maximum temperature. On the other hand, the conversion of wooded grassland to cropland resulted in a similar daily warming and drying but exclusively due to vegetation physiological activity.The area-averaged monthly mean response for each conversion type resulted in a warming of about 0.6°C for the conversion of broadleaf evergreen and 1.2°C for conversion of wooded grassland. These temperature differences represent an augmentation in the local heat source associated with a reduction in evapotranspiration due to land cover conversion and do not reflect variations forced by changes in atmospheric circulation.When averaged over the entire domain, the effect of landscape conversion results in a reduction of the latent heat flux and an increase in sensible heat flux, producing a large-scale apparent heat source of 0.5°C during January. This warming is in line with an increasing trend observed in monthly mean temperature in Santa Cruz, Bolivia during the same period.     

Amazonian forest dieback under climate-carbon cycle projections for the 21st century

Summary The first GCM climate change projections to include dynamic vegetation and an interactive carbon cycle produced a very significant amplification of global warming over the 21st century. Under the IS92a business as usual emissions scenario CO₂ concentrations reached about 980ppmv by 2100, which is about 280ppmv higher than when these feedbacks were ignored. The major contribution to the increased CO₂ arose from reductions in soil carbon because global warming is assumed to accelerate respiration. However, there was also a lesser contribution from an alarming loss of the Amazonian rainforest. This paper describes the phenomenon of Amazonian forest dieback under elevated CO₂ in the Hadley Centre climate-carbon cycle model.     

The UHF wind profiler at Vienna airport – data quality control and comparisons to rawinsonde data

Summary A 1290MHz wind profiler (Radian Lap-3000), at present one of three operational wind profilers in Austria, is operated at Vienna airport. In spite of quality assurance procedures as consensus averaging included in the data evaluation process from profiler raw data, some spurious peaks of wind speed and unrealistic changes of the wind vector in time or height occur in the wind measurements. This is especially true for sampling intervals of only 5 minutes which are used to resolve the temporal evolution of summer thunderstorms and frontal passages. Averaging periods of only a few minutes are rather the lower limit apt for wind profiler observations and result in a low data availability of 28%, whereas about 55% of data (relative to the maximum height range according to the parameter setting) are available for 10 to 30 minutes profiles.Approaches to a posteriori quality control using checks for automatic error detection are proposed and tested on a one and a half year data-set: Flagging data when the three-dimensional wind divergence exceeds a predefined limit (0.5s^–1) is in most cases successful in combination with thresholds for wind speed (2 times the median of the daily data-set) or wind shear (0.2s^–1).The wind profiler data is compared to wind profiles from the next radiosonde station where soundings are launched 4 times a day at Hohe Warte, approx. 20km northwest, at the hill-side of the Viennese Woods. Deviations of about 1m s^–1 in wind speed are found between the observations of the two systems. Differences between the wind profiles within the boundary layer can be explained by local differences in the wind regime observed at the airport and the radiosounding – blocking effects of the Viennese Woods during south-easterly flow. Comparing the profiler data to radiosoundings on a monthly basis gives a tool to monitor the profiler performance.     

Remote sensing of the dynamic stability of the atmospheric boundary layer

Summary The relative strength of the stabilizing effect of buoyancy and the destabilizing effect of velocity shear in a stratified shear flow, such as a stable atmospheric boundary layer, is measured by the gradient Richardson number, Ri_g. The boundary layer static stability, as described by the buoyancy frequency, N, can be calculated from the virtual potential temperature gradient derived from RASS temperature profiles. The mean wind profiles from a sodar can be used to calculate the mean vertical velocity shear. In combination these profilers are potentially a powerful tool for the remotely sensing the dynamic stability of the boundary layer. However, experience shows that the combinations of two experimentally derived quantities, like N and shear, may give highly variable results. On the other hand, a simple sensitivity analysis shows that reasonable estimates of Ri_g are achievable over a range of conditions in the stable nocturnal boundary layer. To test this conclusion, high spatial and temporal resolution temperature and velocity soundings were obtained above 50m in the stable nocturnal boundary layer using a 920MHz continuous wave Radio Acoustic Sounding System (RASS) and 1.875kHz and 5.00kHz Doppler sodars. Examples of the evolution of Ri_g are presented from 24 hours of observations of the boundary layer in Canberra, on the tablelands in south- eastern Australia. Most of the boundary layer had Ri_g between 0.1 and 1. Thus, it was marginally dynamically stable, even with the gradient Richardson number calculated from finite differences over a vertical interval of 68m. A comparison of the results from the two sodars showed that the velocity shear increased significantly when the vertical differencing interval was decreased from 68m to 20m.     

Rain height in relation to 0°C isotherm height for satellite communication over the Indian Subcontinent

Summary ¶The 0°C isotherm height, a parameter needed for the estimation of attenuation of microwave and millimetre wave for earth-space communication, has been estimated for different stations spread over India. The variations of 0°C isotherm height for different seasons over these stations are presented. Attenuations of radio wave due to rain at frequencies 10GHz and above have also been estimated for few stations using the 0°C isotherm height so derived. The results are useful for radio systems designers.     

Environmental aerosols and their effect on the Earth’s local fair-weather electric field

Summary The Earths local fair-weather electric field is significantly affected by small ions present in the atmosphere. These ions are typically smaller than 0.001µm and occur in concentrations from 500 to 600cm^–3 in air. Attachment to larger aerosol particles may severely decrease the mobility of these atmospheric ions resulting in an increased local electric field. The number concentration of environmental aerosol particles in the size range 0.1 to 5.0µm was measured with two automatic laser scattering particle counters. The Earths electric field was monitored with an electric fieldmeter. Measurements were made in clean air and in an environment highly polluted by wood smoke. The electric field was found to be positively correlated to the aerosol number concentration. During one 24-hour period of measurement, the electric field increased from 180 to about 280Vm^–1 as the number concentration of aerosols larger than 0.1µm increased from about 2000 to 9000cm^–3. The number concentrations of aerosols larger than 0.1 and 0.3µm were both found to be positively correlated with the Earths electric field with correlation coefficients of 70% and 61%, respectively.Present address: School of Physical Sciences, Queensland University of Technology, Brisbane 4001, Australia.     

The recent characteristics of Asian dust and haze events in Seoul, Korea

Summary The study on the characteristics of aerosol in Seoul during springtime from 1998 to 2003 is performed by the size-resolved number concentrations of aerosol. Asian dust events occur in spring most frequently, but it has been often observed in wintertime since 1999. Since 2000, the number of Asian dust days has been increasing, and the intensity has been more severe until 2002. However, there were only 3 dust days in Seoul during the spring of 2003, since the synoptic cyclone was relatively not intense enough to rise and transport dust to Korean peninsula, and the air stream was usually tiled to north of Korean peninsula. In addition, the precipitation was relatively plentiful and the air temperature was cold enough not to keep dry soil condition.Haze is the suspended particles in the air, reducing visibility by scattering light, and it is often a mixture of aerosols and photochemical smog. Dry particles with diameters of the order of 0.1µm, are small enough to scatter short wavelengths of light. Haze occurs well in winter and spring, and severe haze is observed in the afternoon. The occurrence frequency of haze has been decreasing since 2000 except in May of 2003.During Asian dust events from 1998 to 2003, the number concentration of aerosol with diameters from 0.3µm to 0.5µm decreases notably, but that larger than 1µm increases rapidly. On the other hand, for the haze events the number concentration from 0.3µm to 0.5µm increases notably, but that larger than 1µm decreases.     

Hailstorms in Northern Greece: synoptic patterns and thermodynamic environment

Summary ¶Two cyclonic vortices close to each other, a binary cyclone or binary system, tend to rotate cyclonically relative to one another and to merge, i.e. the Fujiwhara effect. The point vortex model that represents barotropic binary cyclones predicts their rotation features as follows. The rotation rate is proportional linearly to the sum of the cyclones intensities and inversely to the square of their separation distance while the more intense cyclone rotates slower. Our earlier observational analysis of 1423 mid-latitude binary cyclones (Ziv and Alpert, 1995) showed a reasonable fit to theory, except for the absence of a correlation between individual speeds and intensities within the binary systems, and a reversal of the inverse rotation-separation relationship at the range of 1400–1800km.This study is the first attempt to describe the mid-latitude binary systems using potential vorticity concepts (PV thinking), which implies that a binary interaction takes place between the 3-D flow patterns induced by upper-PV or surface-thermal anomalies rather than by the surface cyclones alone. It is argued that the upper-anomalies dominate the rotation process, and hence the rotational speeds of the interacting surface cyclones are more closely correlated with the relative intensities of their corresponding upper-level anomalies rather than with their own intensities, as reflected in weather charts. Data analysis indicates that mid-latitude binary cyclones are normally associated with at least one upper-PV anomaly. This explains the absence of a correlation between the rotation speed and the intensity of the surface cyclones there.A unique type of a mid-latitude binary system is identified, in which one cyclone coincides with an upper major PV-anomaly and the other moves along the periphery of the former. Such a binary system is entitled here the Contact Binary System (CBS), in contrast with remote interacting systems implied by the point vortex theory.Analytical considerations yield an increase in the rotation rate with separation for CBSs of separation smaller than 1000–1500km, in contrast to the normal decrease with R ². The contribution of CBSs is suggested here to explain the abnormal increase in rotation rate at 1400–1900km range.     

Reduction of forecast error for global numerical weather prediction by The Florida State University (FSU) Superensemble

Summary The skill of the FSU Superensemble technique as applied to global numerical weather prediction is evaluated extensively in this paper. The global mass and motion fields for year 2000 and precipitation over the domain 55S to 55N for year 2001, as predicted by the Superensemble, the ensemble member models, and the mean of the ensemble members, are evaluated by standard statistical measures of skill to determine the performance of the Superensemble in relation to the other models. The member models are global forecast models from 5 of the worlds operational forecast centers in addition to the FSU global spectral model. For precipitation 5 additional versions of the FSU global model are utilized in the ensemble, as defined by different initial conditions provided by various physical initialization algorithms. Statistical parameters calculated for the mass and motion fields include root mean square (RMS) error, systematic error (or bias), and anomaly correlation. These are applied to the mean sea level pressure, 500hPa heights, and the wind fields at 850hPa and 200hPa. Statistical parameters that were calculated for precipitation include RMS error, correlation, equitable threat score (ETS), and a special definition of bias appropriate for the precipitation field. For the mass and motion fields the performance of the Superensemble was considered for the annual global case, as well as for each hemisphere (north and south) and for each of the four seasons. For precipitation only the annual case was considered over the domain cited above.For the mass and motion fields the RMS calculations showed the Superensemble to be superior (to have the smallest total forecast error) in all comparisons to the ensemble member models, and to be superior to the ensemble mean in the vast majority of comparisons. Performance in comparison to the other models was generally better in the Southern Hemisphere than in the Northern Hemisphere, and better in the transition seasons of fall and spring than in the extreme seasons of winter and summer. The Superensemble had the best success with mean sea level pressure, followed in order by 500hPa geopotential heights, 850hPa winds, and 200hPa winds.In the calculations of 500hPa geopotential height anomaly correlation the Superensemble had higher scores in all comparisons to the ensemble member models, as well as higher scores in the majority of comparisons to the ensemble mean. As with the RMS error results, the Superensemble performed better in the Southern Hemisphere than in the Northern Hemisphere, and better in fall than in summer, in comparison to the other models. The superior anomaly correlation scores of the Superensemble attest to the ability of the model to forecast daily perturbations from the climatological means, perturbations that are associated with transient synoptic scale features, given the horizontal resolution in the forecast models.In terms of systematic error reduction the Superensemble produces its most impressive results. Annual global mean sea-level pressure systematic errors for day 5 forecasts are generally in the range of ±1hPa (compared to errors as high as 8hPa in other models), and day 2 forecasts of 500hPa geopotential height produced systematic errors generally in the range of ±10 meters (compared to errors as high as 60 meters in other models). The Superensemble was able to reduce systematic errors in forecasts of a variety of important features in the global mass and motion fields: surface equatorial trough, wave amplitude in geopotential heights at 500hPa, trade winds and Somali Jet at 850hPa, mid-latitude westerlies, subtropical jet, and Tropical Easterly Jet (TEJ) at 200hPa.In terms of forecasting precipitation the Superensemble outperforms all ensemble member models and the ensemble mean in terms of RMS error, correlation coefficient, equitable threat score, and bias. The superior correlation scores indicate that the Superensemble is more reliable than the other models in predicting perturbations in the area distribution of precipitation, perturbations that are essentially associated with migrant synoptic scale disturbances, considering the horizontal resolution of the forecast models.The Superensemble is a valuable tool for significantly improving upon the global model forecasts of the worlds operational forecast centers. These forecasts are used daily as important guidance in making weather forecasts in all regions of the world. This paper will demonstrate that the Superensemble improves upon the ensemble member model forecasts: (1) in a statistical sense considering broad areas of the globe, (2) in a synoptic climatology sense through focus on the improved forecasts of climatological features seen in the global mass and motion fields, (3) in a synoptic sense through use of anomaly correlation and correlation coefficient where improvement is demonstrated in the forecasts of perturbations from mean fields which are essentially associated with transient synoptic scale disturbances.     

Estimates of snow accumulation and volume in the Swiss Alps under changing climatic conditions

Summary ¶Snow is a key feature of mountain environments in terms of the controls it exerts on hydrology, vegetation, and in terms of its economic significance (e.g. for the ski industry). Its quantification in a changing climate is thus important for various environmental and economic impact assessments. Based on observational analysis, surface energy balance modeling, and the latest data from high-resolution regional climate models, this paper investigates the possible changes in snow volume and seasonality in the Swiss Alps. An average warming of 4°C as projected for the period 2071–2100 with respect to current climate suggests that snow volume in the Alps may respond by reductions of at least 90% at altitudes close to 1000m, by 50% at 2000m, and 35% at 3000m. In addition, the duration of snow cover is sharply reduced in the warmer climate, with a termination of the season 50–60 days earlier at high elevations above 2000–2500m and 110–130 days earlier at medium elevation sites close to the 1000m altitude. The shortening of the snow season concerns more the end (spring) rather than the beginning (autumn), so that it should be expected that snow melt will intervene much earlier in the season than under current conditions. The results of this study are of relevance to the estimations of the impacts that the projected warming may have on the amount and timing of water in hydrological basins, on the start of the vegetation season, and on the financial status of many mountain resorts.     

Analysis of the Saharan dust regional transport

Summary The Balkan Peninsula is situated in the impact zone of Saharan dust storms. The case of Saharan dust transport to Belgrade in the period of 14–17 April 1994 is analyzed using the Eta model for synoptic and meso scale processes. Air back trajectories are calculated at six model levels from 434 up to 5129m with horizontal grid resolutions of 1°×1° and 10×10. Following cyclonic circulation the dust was picked up from North Africa, and transported over Mediterranean. Simultaneously, according to the analysis of the three lowest trajectories, transport of trace metals from Macedonia and southern Serbia by the Koshava wind might be dominant in the observed episode. Turbulent flow enhanced the coagulation process of initially clean dust particles with particles containing Pb and Cd. The coagulation and scavenging processes below and in clouds increased deposition rates of Pb and Cd in Belgrade in the course of wet removal, and consequently trough resuspension processes. Dry deposition samples contained characteristic particles up to 30µm in diameter with Fe content of 11 to 15 atomic% and significant ratio Si/Fe of 3 to 5, determined for selected single particles by the SEM/EDX method. Following dry and wet deposition of Cd and Pb, a residual effect of dustfall is noticed throughout the vegetational period.