Effect of land management on ecosystem carbon fluxes at a subalpine grassland site in the Swiss Alps

Summary The influence of agricultural management on the CO₂ budget of a typical subalpine grassland was investigated at the Swiss CARBOMONT site at Rigi-Seebodenalp (1025m a.s.l.) in Central Switzerland. Eddy covariance flux measurements obtained during the first growing season from the mid of spring until the first snow fall (17 Mai to 25 September 2002) are reported. With respect to the 10-year average 1992–2001, we found that this growing season had started 10 days earlier than normal, but was close to average temperature with above-normal precipitation (100–255% depending on month). Using a footprint model we found that a simple approach using wind direction sectors was adequate to classify our CO₂ fluxes as being controlled by either meadow or pasture. Two significantly different light response curves could be determined: one for periods with external interventions (grass cutting, cattle grazing) and the other for periods without external interventions. Other than this, meadow and pasture were similar, with a net carbon gain of –128±17g Cm^–2 on the undisturbed meadow, and a net carbon loss of 79±17g Cm^–2 on the managed meadow, and 270±24g Cm^–2 on the pasture during 131 days of the growing season, respectively. The grass cut in June reduced the gross CO₂ uptake of the meadow by 50±2% until regrowth of the vegetation. Cattle grazing reduced gross uptake over the whole vegetation period (37±2%), but left respiration at a similar level as observed in the meadow.     

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.     

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.     

SODAR-based estimation of TKE and momentum flux profiles in the atmospheric boundary layer: Test of a parameterization model

Summary A simple parameterization for the estimation of turbulent kinetic energy (TKE) and momentum flux profiles under near-neutral stratification based on sodar measurements of the vertical velocity variance has been tested using data from the LINEX-2000 experiment. Measurements included operation of a phased-array Doppler sodar DSDPA.90 and of a sonic anemometer USA-1 mounted at a meteorological tower at a height of 90m. Good agreement has been found between the TKE and momentum flux values derived from the sonic and sodar data (with correlation coefficients r>0.90 and a slope of the regression lines of about 1.01.1) suggesting the possible use of sodar measurements of _w ² to derive turbulence parameter profiles above the tower range.     

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.     

Recent warming in a small region with semi-oceanic climate, 1949–1998: what is the ground truth?

Summary Trends of monthly air temperature extremes were investigated in five meteorological stations of the Grand-Duchy of Luxembourg during the period 1949–1998. The application of an innovative homogenization method based on the concept of relative homogeneity to climatic time series allows identifying multiple break points, as well as correcting data series in an objective and robust statistical way. The rise of maximum temperature (Tmax) has occurred at a rate of 1.5 times that of the minimum temperature (Tmin) in winter (+1.4°C versus +0.9°C) and summer (+1.4°C versus +0.8°C). No trend in temperature extremes was found in autumn, while spring was affected by a small warming (+0.3°C) of Tmin and no change in Tmax resulting in a decrease of the diurnal temperature range (DTR) (–0.3°C). In spring, a strong positive linear relationship between Tmin warming and local terrain slope could be found. Comparison to new-gridded large-scale climatologies indicates generally close agreement to temperature trends during the 1949–1998 period, while a lower local warming was observed in summer during the post-1975 period following the changing-point year of atmospheric circulation over North-western Europe. This study shows that the question of data homogeneity is not trivial and should receive careful attention before quantifying historical temperature trends and identifying their spatial patterns at regional scale.     

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.     

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.     

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.     

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.     

Tropical cyclone environments over the northeastern Pacific, including mid-level dry intrusion cases

Summary This study uses a 1°×1° lat/long dataset, extracted from ECMWF re-analyses for the 15-year period 1979–1993 (ERA-15), to diagnose the synoptic-scale kinematic, thermodynamic and moisture environments in the vicinity of named tropical cyclones (TCs) in the eastern North Pacific. Based on the NCDC best track dataset, TCs are partitioned into one of three categories: weak (W), strong (S) or intensifying (I). In total, 63TCs are examined: 8Ws and 20Is at point A (maximum intensification) and 11Ws, 13Ss and 11Is at point B (maximum frequency). Composite maps are compiled for all five groups, and six individual case studies are examined, four for extreme TC cases and two for cases involving dry air intrusions.For the most part, peak values and patterns of composited ERA-15 variables display circulation, thermodynamic and moisture characteristics that are compatible with the strength represented by a groups classification. Intercomparison between Ws and Is at points A and B yielded larger conditional instability of low-level air parcels and upper-level outflow within the region of maximum intensification (point A).The intrusions of dry versus moist mid-level air are addressed for each storm with the assistance of 72-hour backward trajectories. Trajectory density maps indicate two preferred paths of air parcels that reach the environment of W storms at point A on the 700 and 500hPa levels. The first one crossed Central America in the region of the Isthmus of Tehuantepec and the second one south of the Central American mountains. Several storms revealed that these trajectories were associated with dry air intrusions into the larger storm area, and this might be one reason for their weak status at point A. One documented example is Kevin (1985). By the time it reached point B, the dry air was replaced by air that was moist and Kevin intensified, although it remained a W system. In contrast, Narda (1989) received a dry air intrusion from Central Mexico at 500hPa as a weak storm at point B and did not intensify. Despite possible analyses problems, the documentation in this study of mid-level dry air intrusions into eastern Pacific TCs from the Mexican-Central American region suggests a hitherto unexploited forecast potential. Received January 15, 2002; revised November 28, 2002; accepted December 19, 2002 Published online: May 8, 2003     

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.     

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.     

Outgoing long-wave radiation over the Tropical Pacific and Atlantic Ocean and Indian summer monsoon rainfall

Summary In this paper, the interannual variability of satellite derived outgoing longwave radiation (OLR) is examined in relation to the Indian summer monsoon rainfall (June to September total rainfall; ISMR). Monthly grid point OLR field over the domain i.e. the tropical Pacific and Atlantic region (30°N to 30°S, 110°E to 10°W) and the ISMR for the period 1974–2001 are used for the study. A strong and significant north–south dipole structure in the correlation pattern is found between the ISMR and the OLR field over the domain during January. This dipole is located over the west Pacific region with highly significant negative (positive) correlations over the South China Sea and surrounding region (around north-east Australia). The dipole weakens and moves northwestward during February and disappears in March. During the month of May, the OLR over the central Atlantic Ocean shows a significant positive relationship with the ISMR. These relationships are found to be consistent and robust during the period of analysis and can be used in the prediction of the ISMR.A multiple regression equation is developed, using the above results, for prediction of the ISMR and the empirical relationships are verified using an independent data set. The results are encouraging for the prediction of the ISMR. The composite annual cycle of the OLR, over the west Pacific regions during extreme ISMR is found to be useful in the prediction of extreme summer monsoon rainfall conditions over the Indian subcontinent.     

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.     

Statistical analysis and inter-comparison of the solar UVB, UVA and global radiation for Beer Sheva and Neve Zohar (Dead Sea), Israel

Summary A statistical analysis and inter-comparison of the solar UVB, UVA and global radiation for Beer Sheva and Neve Zohar (Dead Sea) are presented utilizing data measured from January 1995 through December 2002. Beer Sheva is located approximately 65km to the west of and 700m above Neve Zohar. The monthly average hourly and daily values for all radiation types at both sites are reported. The standard errors of the monthly average daily values have been calculated in order to ascertain whether the average daily radiation intensities are representative, i.e. if the magnitude of the standard error is less than the inherent measurement uncertainty of the instruments, and, thereby, justify an inter-comparison between the two sites. The relative magnitude of the global, UVB and UVA radiation intensity at the two sites is attributed to the enhanced scattering of the incident solar radiation at the Dead Sea location due to the longer optical path length it must traverse to arrive at the Dead Sea, the lowest terrestrial point on earth. The degree of attenuation of solar radiation due to the scattering phenomena is inversely proportional to the wavelength raised to some power and, consequently, it is greatest for UVB and negligible for global radiation.     

Sensitivity of freezing drizzle formation in stably stratified clouds to ice processes

Summary This study examines the impact of ice formation and growth processes on freezing drizzle formation in stably stratified clouds. In particular we investigate the reason why freezing drizzle is rarely observed in clouds with top temperatures less than –15°C. We also investigate the sensitivity of freezing drizzle formation to the Hallett Mossop secondary ice process (Hallet and Mossop, 1974). The evaluation is performed by simulating cloud formation over a two-dimensional idealized mountain using a detailed microphysical scheme. The height and width of the two-dimensional mountain were designed to produce an updraft pattern with extent and magnitude similar to documented freezing drizzle cases. The simulations show that: (i) drizzle formation is very sensitive to the ice crystal concentration, with a significant reduction in the area over which drizzle forms and the maximum drizzle water content as the cloud top temperature decreases below –10°C, and (ii) secondary ice crystal formation has a significant effect on drizzle formation at cloud top temperatures below –10°C.The above two factors are likely the main cause for the lack of freezing drizzle at cloud top temperatures less than –15°C. We also found that neglecting the depletion of ice forming nuclei resulted in considerable overestimation of the ice crystal concentration and suppression of drizzle, even for the –10°C case.     

The convective boundary layer over pasture and forest in Amazonia

Summary The coupling between different types of surface (tropical forest or grass) and the Convective Boundary Layer (CBL) has been investigated using observational (rawinsoundings) data collected over Rondônia in southwest Amazonia. The data reported here support the notion that deforestation may modify the dynamics of the boundary layer, in particular during the dry season. In this period the sensible heat fluxes are very high over pasture, creating a CBL around 550m deeper compared to that over the forest. The measurements showed the height of the fully developed CBL for pasture to be 1650m, compared to around 1100m for forest. During the wet season the height of the CBL is lower than during the dry season and has the same height (around 1000m) for forest and pasture sites. The CBL over pasture is hotter and drier than over forest during the dry season, but during the wet season the air temperatures and humidities are similar. Comparing the CBL growth during the dry and wet season, there is evidence that the CBL properties over the forest are not dependent on the surface characteristics, but over the pasture they are.     

Feature-based predictability of 500 hpa height in the Australia-New Zealand region

Summary Feature-based predictability stratifies forecast model errors on the basis of individual weather systems. We examine only one level and chose very simple categories: high, cut-off low, trough and block. European Centre for Medium-range Weather Forecasts data are used. We emphasize systems found near Australia and New Zealand during winter 1987. Calculations for the preceding summer and fall and for other midlatitude regions of the southern hemisphere yield similar results.The approach herein is fully automated and simple to implement. Features are identified in the verification field. Then an error calculation is made on moving grids that are each centered upon and contain one system.The total error is the root mean square difference (RMS) between forecast and verification. The structural error is the RMS difference when the forecast and verification small grids are independently centered upon the corresponding feature in each field. The difference between the structural and total errors, called the locational error, is typically a quarter of the total. Even when normalized by presistence, highs are forecast better than troughs; cut-off lows and troughs have similar errors. The distance between the forecast and observed positions is typically 3° longitude west and 0.5° latitude north of where features should be at 72 hours. The model has a systematic bias of too small amplitude of vorticity. No relation is found between skill and jet stream splitting.With 4 Figures     

Objective climatology of 500-hPa cyclones in central and east mediterranean region during warm-dry period of the year

Summary ¶An Objective 500-hPa cyclone detection and analysis is performed during the warm-dry period (16 April to 15 October) of the year for the Central and Eastern Mediterranean Region (30°N–50°N, 5°E–35°E). The 40-year NCEP/NCAR reanalyzes gridded data of geopotential height and temperature employed in the study, enable a climatic approach with a spatial (2.5°×2.5°) and a 6-hour temporal (00, 06, 12, 18 UTC) resolution. The occurrence frequency for the entire region shows a high interannual variability without significant trends. The frequency maximization over land in the middle of the warm period indicates a possible relationship with land-air temperature difference. Three primary activity centers (Turkey, Black Sea and Genoa) are objectively determined in the frequency domain as local maxima and their intensity is assessed in terms of average geopotential height and temperature values. The Turkey center, though the most frequent, consists of shallow lows in contrast to the Black Sea center, which is second in frequency, but on average, contains the deepest lows amongst the three. The Genoa center, which comes third in frequency, shows the highest variability in the intensity of the lows. Frequency distributions in the area, during the three sub periods (Early, Middle and Late Warm Period) and during the four synoptic hours, provide insight in to the origins of these systems at the three locations and reveal secondary centers, such as the centers in the Adriatic (at 18 UTC), Northern Greece (during Middle Warm Period) and the South Aegean Sea (during Late Warm Period).Received February 20, 2002; revised November 4, 2002; accepted January 16, 2003 Published online July 30, 2003