Search

The cloud mask is an essential product derived from satellite data. Whereas cloud analysis applications typically make use of information from cloudy pixels, many other applications require cloud-free conditions. For this reason many organizations have their own cloud masks tuned to serve their particular needs. Being a fundamental product, continuous quality monitoring and validation of these cloud masks are vital. This study evaluated the performance of the European Organisation for the Exploitation of Meteorological Satellites (EUMETSAT) Meteorological Products Extraction Facility cloud mask (MPEF), together with the Nowcasting Satellite Application Facility (SAFNWC) cloud masks provided by Meteo-France (SAFNWC/ MSG) and the Swedish Meteorological and Hydrological Institute (SAFNWC/PPS), in the high-latitude area of greater Helsinki in Finland. The first two used the Spinning Enhanced Visible and Infrared Imager (SEVIRI) instrument from the geostationary Meteosat-8 satellite, whereas the last used the Advanced Very High Resolution Radiometer (AVHRR) instrument on board the polar-orbiting NOAA satellite series. Ceilometer data from the Helsinki Testbed, an extensive observation network covering the greater Helsinki area in Finland, were used as reference data in the cloud mask comparison. A computational method, called bootstrapping, is introduced to account for the strong temporal and spatial correlation of the ceilometer observations. The method also allows the calculation of the confidence intervals (CI) for the results. This study comprised data from February and August 2006. In general, the SAFNWC/MSG algorithm performed better than MPEF. Differences were found especially in the early morning low cloud detection. The SAFNWC/PPS cloud mask performed very well in August, better than geostationary-based masks, but had problems in February when its performance was worse. The use of the CIs gave the results more depth, and their use should be encouraged. DOI: 10.1175/2010JAMC2442.1

A morphometric analysis of a southern European city and the derivation of relevant fluid dynamical parameters for use in urban flow and dispersion models are explained in this paper. Calculated parameters are compared with building statistics that have already been computed for parts of three northern European and two North American cities. The aim of this comparison is to identify similarities and differences between several building configurations and city types, such as building packing density, compact versus sprawling neighborhoods, regular versus irregular street orientation, etc. A novel aspect of this work is the derivation and use of digital elevation models (DEMs) for parts of a southern European city. Another novel aspect is the DEMs' construction methodology, which is low cost, low tech, and of simple implementation. Several building morphological parameters are calculated from the urban DEMs using image processing techniques. The correctness and robustness of these techniques have been verified through a series of sensitivity tests performed on both idealized building configurations, as well as on real case DEMs, which were derived using the methodology here. In addition, the planar and frontal area indices were calculated as a function of elevation. It is argued that those indices, estimated for neighborhoods of real cities, may be used instead of the detailed building geometry within urban canopy models as those indices together synthesize the geometric features of a city. The direct application of these results will facilitate the development of fast urban flow and dispersion models. DOI: 10.1175/2010JAMC2117.1

The Atmospheric Infrared Sounder (AIRS) and the Infrared Atmospheric Sounding Interferometer (IASI), together with the future Cross-track Infrared Sounder, will provide long-term hyperspectral measurements of the earth and its atmosphere at ~10 km spatial resolution. Quantifying the radiometric difference between AIRS and IASI is crucial for creating fundamental climate data records and establishing the space-based infrared calibration standard. Since AIRS and IASI have different local equator crossing times, a direct comparison of these two instruments over the tropical regions is not feasible. Using the Geostationary Operational Environmental Satellite (GOES) imagers as transfer radiometers, this study compares AIRS and IASI over warm scenes in the tropical regions for a time period of 16 months. The double differences between AIRS and IASI radiance biases relative to the GOES-11 and -12 imagers are used to quantify the radiance differences between AIRS and IASI within the GOES imager spectral channels. The results indicate that, at the 95% confidence level, the mean values of the IASI - AIRS brightness temperature differences for warm scenes are very small, that is, -0.0641 [+ or -] 0.0074 K, -0.0432 [+ or -] 0.0114 K, and -0.0095 [+ or -] 0.0151 K for the GOES-11 6.7-, 10.7-, and 12.0-[micro]m channels, respectively, and -0.0490 [+ or -] 0.0100 K, -0.0419 [+ or -] 0.0224 K, and -0.0884 [+ or -] 0.0160 K for the GOES-12 6.5-, 10.7-, and 13.3-[micro]m channels, respectively. The brightness temperature biases between AIRS and IASI within the GOES imager spectral range are less than 0.1 K although the AIRS measurements are slightly warmer than those of IASI. DOI: 10.1175/2009JAMC2218.1

Two approaches to producing a hail climatology for Finland are compared. The first approach is based on 70 yr of hail reports from different sources (newspapers, storm spotters, and other volunteers). The second is derived primarily from radar data. It is shown that a selection of newspaper articles of hail damage covering a period of 70 yr provides a good overview of the typical monthly and diurnal distribution of hail occurrence over the country. Radar data covering five summers (2001-05) provide another data source, but with different potential sources of errors. The two distinct methods compared in this paper give roughly the same results in describing the hail climatology of Finland, which gives additional confidence in each of the methods. On the basis of both methods, most hailstones are observed in the afternoon, 1400-1600 local time. The hail 'season' extends from May to early September with maximum occurrences in June, July, and August. This means that hail is most frequently observed when the convective energy available for storm growth is at its diurnal or seasonal peak. The length of the hail season is the same according to both radar and newspaper data. The main difference emerges in relation to July and August events: 37% of news about hail events is published in newspapers in late July but only 8% in early August, whereas for radar data the numbers are more evenly distributed, 33% and 18%, respectively. This can be partially explained by sociological factors--July is the main holiday month in Finland, when outdoor activities in more remote areas are more popular. DOI: 10.1175/2009JAMC2116.1

This paper details a number of existing formulations used in Gaussian models in a clear and usable way, and provides a comparison within a single framework--the Gaussian plume and puff models of the air quality modeling system Polyphemus. The emphasis is made on the comparison between 1) the parameterizations to compute the standard deviations and 2) the plume rise schemes. The Gaussian formulas are first described and theoretically compared. Their evaluation is then ensured by comparison with the observations as well as with several well-known Gaussian and computational fluid dynamics model performances. The model results compare well to the other Gaussian models for two of the three parameterizations for standard deviations, Briggs's and similarity theory, while Doury's shows a tendency to underestimate the concentrations because of a large horizontal spread. The results with the Kincaid experiment point out the sensitivity to the plume rise scheme and the importance of an accurate modeling of the plume interactions with the inversion layer. Using three parameterizations for the standard deviations and the same number of plume rise schemes, the authors were able to highlight a large variability in the model outputs. DOI: 10.1175/2009JAMC2160.1

In this study, robust parametric regression methods are applied to temperature and precipitation time series in Switzerland and the trend results are compared with trends from classical least squares (LS) regression and nonparametric approaches. It is found that in individual time series statistically outlying observations are present that influence the LS trend estimate severely. In some cases, these outlying observations lead to an over-/underestimation of the trends or even to a trend masking. In comparison with the classical LS method and standard nonparametric techniques, the use of robust methods yields more reliable trend estimations and outlier detection.

Monthly and daily products of the Global Precipitation Climatology Project (GPCP) are evaluated through a comparison with Finnish Meteorological Institute (FMI) gauge observations for the period January 1995-December 2007 to assess the quality of the GPCP estimates at high latitudes. At the monthly scale both the final GPCP combination satellite--gauge (SG) product is evaluated, along with the satellite-only multisatellite (MS) product. The GPCP daily product is scaled to sum to the monthly product, so it implicitly contains monthly-scale gauge influence, although it contains no daily gauge information. As expected, the monthly SG product agrees well with the FMI observations because of the inclusion of limited gauge information. Over the entire analysis period the SG estimates are biased low by 6% when the same wind-loss adjustment is applied to the FMI gauges as is used in the SG analysis. The interannual anomaly correlation is about 0.9. The satellite-only MS product has a lesser, but still reasonably good, interannual correlation (~0.6) while retaining a similar bias due to the use of a climatological bias adjustment. These results indicate the value of using even a few gauges in the analysis and provide an estimate of the correlation error to be expected in the SG analysis over ocean and remote land areas where gauges are absent. The daily GPCP precipitation estimates compare reasonably well at the 1[degrees] latitude x 2[degrees] longitude scale with the FMI gauge observations in the summer with a correlation of 0.55, but less so in the winter with a correlation of 0.45. Correlations increase somewhat when larger areas and multiday periods are analyzed. The day-to-day occurrence of precipitation is captured fairly well by the GPCP estimates, but the corresponding precipitation event amounts tend to show wide variability. The results of this study indicate that the GPCP monthly and daily fields are useful for meteorological and hydrological studies but that there is significant room for improvement of satellite retrievals and analysis techniques in this region. It is hoped that the research here provides a framework for future high-latitude assessment efforts such as those that will be necessary for the upcoming satellite-based Global Precipitation Measurement (GPM) mission.

Cloud-type classification based on multispectral satellite imagery data has been widely researched and demonstrated to be useful for distinguishing a variety of classes using a wide range of methods. The research described here is a comparison of the classifier output from two very different algorithms applied to Geostationary Operational Environmental Satellite (GOES) data over the course of one year. The first algorithm employs spectral channel thresholding and additional physically based tests. The second algorithm was developed through a supervised learning method with characteristic features of expertly labeled image samples used as training data for a 1-nearest-neighbor classification. The latter's ability to identify classes is also based in physics, but those relationships are embedded implicitly within the algorithm. A pixel-to-pixel comparison analysis was done for hourly daytime scenes within a region in the northeastern Pacific Ocean. Considerable agreement was found in this analysis, with many of the mismatches or disagreements providing insight to the strengths and limitations of each classifier. Depending upon user needs, a rule-based or other postprocessing system that combines the output from the two algorithms could provide the most reliable cloud-type classification.

A number of Earth remote sensing satellites are currently carrying passive microwave radiometers. A variety of different retrieval algorithms are used to estimate surface rain rates over the ocean from the microwave radiances observed by the radiometers. This study compares several different satellite algorithms with each other and with independent data from rain gauges on ocean buoys. The rain gauge data are from buoys operated by the NOAA Pacific Marine Environmental Laboratory. Potential errors and biases in the gauge data are evaluated. Satellite data are from the Tropical Rainfall Measuring Mission Microwave Imager and from the Special Sensor Microwave Imager instruments on the operational Defense Meteorological Satellite Program F13, F14, and F15 satellites. These data have been processed into rain-rate estimates by the NASA Precipitation Measurement Mission and by Remote Sensing Systems, Inc. Biases between the different datasets are estimated by computing differences between long-term time averages. Most of the satellite datasets agree with each other, and with the gauge data, to within 10% or less. The biases tend to be proportional to the mean rain rate, but the geographical patterns of bias vary depending on the choice of data source and algorithm. Some datasets, however, show biases as large as about 25%, so care should be taken when using these data for climatological studies.

The characterization of high clouds as performed from selected spaceborne observations is assessed in this article by employing a number of worldwide ground-based lidar multiyear datasets as reference. Among the latter, the ground lidar observations conducted at Lannion, Bretagne (48.7[degrees]N, 3.5[degrees]W), and Palaiseau, near Paris [the Site Instrumental de Recherche par Te1edetection Atmospherique (SIRTA) observatory: 48.7[degrees]N, 2.2[degrees]E], both in France, are discussed in detail. High-cloud altitude statistics at these two sites were found to be similar. Optical thicknesses disagree, and possible reasons were analyzed. Despite the variety of instruments, observation strategies, and methods of analysis employed by different lidar groups, high-cloud optical thicknesses from the Geoscience Laser Altimeter System (GLAS) on board the Ice, Cloud and land Elevation Satellite (ICESat) were found to be consistent on the latitude band 40[degrees]-60[degrees]N. Respective high-cloud altitudes agree within 1 km with respect to those from ground lidars at Lannion and Palaiseau; such a finding remains to be verified under other synoptic regimes. Mean altitudes of high clouds from Lannion and Palaiseau ground lidars were compared with altitudes of thin cirrus from the Television and Infrared Observation Satellite (TIROS) Operational Vertical Sounder (TOVS) Path-B 8-yr climatology for a common range of optical thicknesses (0.1-1.4). Over both sites, the annual altitude distribution of thin high clouds from TOVS Path-B is asymmetric, with a peak around 8-9.5 km, whereas the distribution of high clouds retrieved from ground lidars seems symmetric with a peak around 9.5-11.5 km. Additional efforts in standardizing ground lidar observation and processing methods, and in merging high-cloud statistics from complementary measuring platforms, are recommended.

Snow cover has a strong effect on the surface and lower atmosphere in NWP models. Because the progress of in situ observations has stalled, satellite-based snow analyses are becoming increasingly important. Currently, there exist several products that operationally map global or continental snow cover. In this study, satellite-based snow cover analyses from NOAA, NASA, and the European Organisation for the Exploitation of Meteorological Satellites (EUMETSAT), and NWP snow analyses from the High-Resolution Limited-Area Model (HIRLAM) and ECMWF, were compared using data from January to June 2006. Because no analyses were independent and since available in situ measurements were already used in the NWP analyses, no independent ground truth was available and only the consistency between analyses could be compared. Snow analyses from NOAA, NASA, and ECMWF were similar, but the analysis from NASA was greatly hampered by clouds. HIRLAM and EUMETSAT deviated most from other analyses. Even though the analysis schemes of HIRLAM and ECMWF were quite similar, the resulting snow analyses were quite dissimilar, because ECMWF used the satellite information of snow cover in the form of NOAA analyses, while HIRLAM used none. The differences are especially prominent in areas around the snow edge where few in situ observations are available. This suggests that NWP snow analyses based only on in situ measurements would greatly benefit from inclusion of satellite-based snow cover information.

A procedure to accurately resample spaceborne and ground-based radar data is described and then is applied to the measurements taken from the Tropical Rainfall Measuring Mission (TRMM) precipitation radar (PR) and the ground-based Weather Surveillance Radar-1988 Doppler (WSR-88D or WSR) for the validation of the PR measurements and estimates. Through comparisons with the well-calibrated, non-attenuated WSR at Melbourne, Florida, for the period 1998-2007, the calibration of the PR aboard the TRMM satellite is checked using measurements near the storm top. Analysis of the results indicates that the PR, after taking into account differences in radar reflectivity factors between the PR and WSR, has a small positive bias of 0.8 dB relative to the WSR, implying a soundness of the PR calibration in view of the uncertainties involved in the comparisons. Comparisons between the PR and WSR reflectivities are also made near the surface for evaluation of the attenuation-correction procedures used in the PR algorithms. It is found that the PR attenuation is accurately corrected in stratiform rain but is underestimated in convective rain, particularly in heavy rain. Tests of the PR estimates of rainfall rate are conducted through comparisons in the overlap area between the TRMM overpass and WSR scan. Analyses of the data are made both on a conditional basis, in which the instantaneous rain rates are compared only at those pixels at which both the PR and WSR detect rain, and an unconditional basis, in which the area-averaged rain rates are estimated independently for the PR and WSR. Results of the conditional rain comparisons show that the PR-derived rain is about 9% greater and 19% less than the WSR estimates for stratiform and convective storms, respectively. Overall, the PR tends to underestimate the conditional mean rain rate by 8% for all rain categories, a finding that conforms to the results of the area-averaged rain (unconditional) comparisons.

Studies in the past several years have documented the climatology of warm-season precipitation-episode statistics (propagation speed, span, and duration) over the United States using a national composited radar dataset. These climatological studies have recently been extended to other continents, including Asia, Africa, and Australia. However, continental regions outside the United States have insufficient radar coverage, and the newer studies have had to rely on geostationary satellite data at infrared (IR) frequencies as a proxy for rainfall. It is well known that the use of IR brightness temperatures to infer rainfall is subject to large errors. In this study, the statistics of warm-season precipitation episodes derived from radar and satellite IR measurements over the United States are compared and biases introduced by the satellite data are evaluated. It is found that the satellite span and duration statistics are highly dependent upon the brightness temperature threshold used but with the appropriate choices of thresholds can be brought into good agreement with those based upon radar data. The propagation-speed statistics of satellite events are on average -4 m [s.sup.-1] faster than radar events and are relatively insensitive to the brightness temperature threshold. A simple correction procedure based upon the difference between the steering winds for the precipitation core and the winds at the level of maximum anvil outflow is developed.

The diurnal cycle of outgoing longwave radiation (OLR) computed by a climate model provides a powerful test of the numerical description of various physical processes. Diurnal cycles of OLR computed by version 3 of the Hadley Centre Atmospheric Model (HadAM3) are compared with those observed by the Earth Radiation Budget Satellite (ERBS) for the boreal summer season (June-August). The ERBS observations cover the domain from 55[degrees]S to 55[degrees]N. To compare the observed and modeled diurnal cycles, the principal component (PC) analysis method is used over this domain. The analysis is performed separately for the land and ocean regions. For land over this domain, the diurnal cycle computed by the model has a root-mean-square (RMS) of 11.4 W [m.sup.-2], as compared with 13.3 W [m.sup.-2] for ERBS. PC-1 for ERBS observations and for the model are similar, but the ERBS result has a peak near 1230 LST and decreases very slightly during night, whereas the peak of the model result is an hour later and at night the OLR decreases by 7 W [m.sup.-2] between 2000 and 0600 LST. Some of the difference between the ERBS and model results is due to the computation of convection too early in the afternoon by the model. PC-2 describes effects of morning/afternoon cloudiness on OLR, depending on the sign. Over ocean in the ERBS domain, the model RMS of the OLR diurnal cycle is 2.8 W [m.sup.-2], as compared with 5.9 W [m.sup.-2] for ERBS. Also, for the model, PC-1 accounts for 66% of the variance, while for ERBS, PC-1 accounts for only 16% of the variance. Thus, over ocean, the ERBS results show a greater variety of OLR diurnal cycles than the model does.

This study employs principal component analysis to correct tropical precipitation estimates in the University of Wisconsin--Milwaukee (UWM)/Comprehensive Ocean-Atmosphere Dataset (COADS) data. The idea was to use a matrix made up of the other variables in the set, to reduce the dimensionality of the matrix by considering a small number of principal components, and then to regress precipitation to these principal components. The results indicate that, although some information on precipitation could be restored by this method, overall the resulting precipitation estimates are not reliable. This result is traced to the intrinsic complexity of precipitation and possibly to a newly discovered bias in the UWM/COADS data.

Byline: Per-Erik Johansson The performance of two meteorological pre-processors, which calculate the surface fluxes of heat and momentum from routine meteorological observations, was studied during two winter months at high latitudes, characterised by low solar altitudes and snow-covered surfaces, i.e. high albedo values. This, together with low temperatures and low wind speeds, causes stable stratification (often with inversion) to be predominant. The pre-processors studied are those in use by the Swedish and Finnish Meteorological Institutes. Comparisons were made between the pre-processors and between each one of the pre-processors and direct turbulence measurements of the fluxes. In addition, the net radiation predicted by the two schemes was compared with measurements. The results show that the calculated momentum fluxes agree quite well with the measurements, but that the pre-processors need to be further developed to give reasonable estimates of the surface heat flux in the studied conditions. The friction velocities calculated by the two pre-processors agree quite well with each other, but there is a large discrepancy between the two schemes for sensible heat flux. One reason for the difference between measured and calculated heat flux values is the parameterisation of net radiation in the pre-processors. The net radiation estimated by the pre-processors was shown to differ from the measured values. This is not surprising, bearing in mind that the pre-processors have been calibrated against measurements at mid-latitudes and with grass-covered surfaces.