On the development of summer thunderstorms in the city of Sao Paulo: Mean meteorological characteristics and pollution effect

This study investigates how the summer thunderstorms developed over the city of Sao Paulo and if the pollution might affect its development or characteristics during the austral summer (December-January-February-March, DJFM months). A total of 605 days from December 1999 to March 2004 was separated as 241 thunderstorms days (TDs) and 364 non-thunderstorm days (NTDs). The analyses are performed by using hourly measurements of air temperature (T), web-bulb temperature (Tw), surface atmospheric pressure (P), wind velocity and direction, rainfall and thunder and lightning observations collected at the Meteorological Station of the University of Sao Paulo in conjunction with aerosol measurements obtained by AERONET (Aerosol Robotic Network), and the NCEP-DOE (National Centers for Environmental Prediction Department of Energy) reanalysis and radiosondes. The wind diurnal cycle shows that for TDs the morning flow is from the northwest rotating to the southeast after 16: 00 local time (LT) and it remains from the east until the night. For the NTDs, the wind is well characterized by the sea-breeze circulation that in the morning has the wind blowing from the northeast and in the afternoon from the southeast. The TDs show that the air temperature diurnal cycle presents higher amplitude and the maximum temperature of the day is 3.2 degrees C higher than in NTDs. Another important factor found is the difference between moisture that is higher during TDs. In terms of precipitation, the TDs represent 40% of total of days analyzed and those days are responsible for more than 60% of the total rain accumulation during the summer, for instance 50% of the TDs had more than 15.5mm day(-1) while the NTDs had 4 mm day(-1). Moreover, the rainfall distribution shows that TDs have higher rainfall rate intensities and an afternoon precipitation maximum; while in the NTDs there isn`t a defined precipitation diurnal cycle. The wind and temperature fields from NCEP reanalysis concur with the local weather station and radiosonde observations. The NCEP composites show that TDs are controlled by synoptic circulation characterized by a pre-frontal situation, with a baroclinic zone situated at southern part of Sao Paulo. In terms of pollution, this study employed the AERONET data to obtain the main aerosol characteristics in the atmospheric column for both TDs and NTDs. The particle size distribution and particle volume size distribution have similar concentrations for both TDs and NTDs and present a similar fine and coarse mode mean radius. In respect to the atmospheric loading, the aerosol optical depth (AOD) at different frequencies presented closed mean values for both TDs and NTDs that were statistically significant at 95% level. The spectral dependency of those values in conjunction with the Angstrom parameter reveal the higher concentration of the fine mode particles that are more likely to be hygroscopic and from urban areas. In summary, no significant aerosol effect could be found on the development of summer thunderstorms, suggesting the strong synoptic control by the baroclinic forcing for deep convective development. (C) 2010 Published by Elsevier B. V.

The mean evolution and variability of the Asian summer monsoon: Comparison of ECMWF and NCEP-NCAR reanalyses

The behavior of the Asian summer monsoon is documented and compared using the European Centre for Medium-Range Weather Forecasts (ECMWF) Reanalysis (ERA) and the National Centers for Environmental Prediction-National Center for Atmospheric Research (NCEP-NCAR) Reanalysis. In terms of seasonal mean climatologies the results suggest that, in several respects, the ERA is superior to the NCEP-NCAR Reanalysis. The overall better simulation of the precipitation and hence the diabatic heating field over the monsoon domain in ERA means that the analyzed circulation is probably nearer reality. In terms of interannual variability, inconsistencies in the definition of weak and strong monsoon years based on typical monsoon indices such as All-India Rainfall (AIR) anomalies and the large-scale wind shear based dynamical monsoon index (DMI) still exist. Two dominant modes of interannual variability have been identified that together explain nearly 50% of the variance. Individually, they have many features in common with the composite flow patterns associated with weak and strong monsoons, when defined in terms of regional AIR anomalies and the large-scale DMI. The reanalyses also show a common dominant mode of intraseasonal variability that describes the latitudinal displacement of the tropical convergence zone from its oceanic-to-continental regime and essentially captures the low-frequency active/break cycles of the monsoon. The relationship between interannual and intraseasonal variability has been investigated by considering the probability density function (PDF) of the principal component of the dominant intraseasonal mode. Based on the DMI, there is an indication that in years with a weaker monsoon circulation, the PDF is skewed toward negative values (i,e., break conditions). Similarly, the PDFs for El Nino and La Nina years suggest that El Nino predisposes the system to more break spells, although the sample size may limit the statistical significance of the results.

A comparison of recent reanalysis datasets using objective feature tracking: Storm tracks and tropical easterly waves

Data from four recent reanalysis projects [ECMWF, NCEP-NCAR, NCEP - Department of Energy ( DOE), NASA] have been diagnosed at the scale of synoptic weather systems using an objective feature tracking method. The tracking statistics indicate that, overall, the reanalyses correspond very well in the Northern Hemisphere (NH) lower troposphere, although differences for the spatial distribution of mean intensities show that the ECMWF reanalysis is systematically stronger in the main storm track regions but weaker around major orographic features. A direct comparison of the track ensembles indicates a number of systems with a broad range of intensities that compare well among the reanalyses. In addition, a number of small-scale weak systems are found that have no correspondence among the reanalyses or that only correspond upon relaxing the matching criteria, indicating possible differences in location and/or temporal coherence. These are distributed throughout the storm tracks, particularly in the regions known for small-scale activity, such as secondary development regions and the Mediterranean. For the Southern Hemisphere (SH), agreement is found to be generally less consistent in the lower troposphere with significant differences in both track density and mean intensity. The systems that correspond between the various reanalyses are considerably reduced and those that do not match span a broad range of storm intensities. Relaxing the matching criteria indicates that there is a larger degree of uncertainty in both the location of systems and their intensities compared with the NH. At upper-tropospheric levels, significant differences in the level of activity occur between the ECMWF reanalysis and the other reanalyses in both the NH and SH winters. This occurs due to a lack of coherence in the apparent propagation of the systems in ERA15 and appears most acute above 500 hPa. This is probably due to the use of optimal interpolation data assimilation in ERA15. Also shown are results based on using the same techniques to diagnose the tropical easterly wave activity. Results indicate that the wave activity is sensitive not only to the resolution and assimilation methods used but also to the model formulation.

Analysis of moisture variability in the European Centre for Medium-Range Weather Forecasts 15-year reanalysis over the tropical oceans

We compare European Centre for Medium-Range Weather Forecasts 15-year reanalysis (ERA-15) moisture over the tropical oceans with satellite observations and the U.S. National Centers for Environmental Prediction (NCEP) National Center for Atmospheric Research 40-year reanalysis. When systematic differences in moisture between the observational and reanalysis data sets are removed, the NCEP data show excellent agreement with the observations while the ERA-15 variability exhibits remarkable differences. By forcing agreement between ERA-15 column water vapor and the observations, where available, by scaling the entire moisture column accordingly, the height-dependent moisture variability remains unchanged for all but the 550–850 hPa layer, where the moisture variability reduces significantly. Thus the excess variation of column moisture in ERA-15 appears to originate in this layer. The moisture variability provided by ERA-15 is not deemed of sufficient quality for use in the validation of climate models.

Inter-comparison of stratospheric mean-meridional circulation and eddy mixing among six reanalysis data sets

The stratospheric mean-meridional circulation (MMC) and eddy mixing are compared among six meteorological reanalysis data sets: NCEP-NCAR, NCEP-CFSR, ERA-40, ERA-Interim, JRA-25, and JRA-55 for the period 1979–2012. The reanalysis data sets produced using advanced systems (i.e., NCEP-CFSR, ERA-Interim, and JRA-55) generally reveal a weaker MMC in the Northern Hemisphere (NH) compared with those produced using older systems (i.e., NCEP/NCAR, ERA-40, and JRA-25). The mean mixing strength differs largely among the data products. In the NH lower stratosphere, the contribution of planetary-scale mixing is larger in the new data sets than in the old data sets, whereas that of small-scale mixing is weaker in the new data sets. Conventional data assimilation techniques introduce analysis increments without maintaining physical balance, which may have caused an overly strong MMC and spurious small-scale eddies in the old data sets. At the NH mid-latitudes, only ERA-Interim reveals a weakening MMC trend in the deep branch of the Brewer–Dobson circulation (BDC). The relative importance of the eddy mixing compared with the mean-meridional transport in the subtropical lower stratosphere shows increasing trends in ERA-Interim and JRA-55; this together with the weakened MMC in the deep branch may imply an increasing age-of-air (AoA) in the NH middle stratosphere in ERA-Interim. Overall, discrepancies between the different variables and trends therein as derived from the different reanalyses are still relatively large, suggesting that more investments in these products are needed in order to obtain a consolidated picture of observed changes in the BDC and the mechanisms that drive them.

Spurious shifts in the pattern of diurnal variation of sea level pressure of reanalysis datasets

The main purpose of this work is to report the presence of spurious discontinuities in the pattern of diurnal variation of sea level pressure of the three reanalysis datasets from: the National Centers for Environmental Prediction (NCEP) and National Center for Atmospheric Science (R1), the NCEP and Department of Energy (R2), and the European Centre for Medium Range Weather Forecasting (ERA-40). Such discontinuities can be connected to the major changes in the global observing system that have occurred throughout reanalyses years. In the R1, the richest period in discontinuities is 1956-1958, coinciding with the start of modern radiosonde observation network. Rapid increase in the density of surface-based observations from 1967 also had an important impact on both R1 and ERA-40, with larger impact on R1. The reanalyses show discontinuities in the 1970s related to the assimilation of radiances measured by the Vertical Temperature Profile Radiometer and TIROS-N Operational Vertical Sounders onboard satellites. In the ERA-40, which additionally assimilated Special Sensor Microwave/Imager data, there are discontinuities in 1987-1989. The R1 also presents further discontinuities, in 1988-1993 likely connected to replacement/introduction of NOAA-series satellites with different biases, and to the volcanic eruption of Mount Pinatubo in June 1991, which is known to have severely affected measurements of infrared radiances for several years. The discontinuities in 1996-1998 might be partially connected to change in the type of radiosonde, from VIZ-B to VIZ-B2. The R2, which covers only satellite era (1979-on), shows discontinuities mainly in 1992, 1996-1997, and 2001. The discontinuities in 1992 and 2001 might have been caused by change in the satellite measurements and those in 1996-1997 by some changes in land-based observations network. © 2012 Springer-Verlag.

BRAZILIAN OFFSHORE WAVE CLIMATE BASED ON NWW3 REANALYSIS

This paper provides a description of the wave climate off the Brazilian coast based on an eleven-year time series (Jan/1997-Dec/2007) obtained from the NWW3 operational model hindcast reanalysis. Information about wave climate in Brazilian waters is very scarce and mainly based on occasional short-term observations, the present analysis being the first covering such temporal and spatial scales. To define the wave climate, six sectors were defined and analyzed along the Brazilian shelf-break: South (W1), Southeast (W2), Central (W3), East (W4), Northeast (W5) and North (W6). W1, W2 and W3 wave regimes are determined by the South Atlantic High (SAH) and the passage of synoptic cold fronts; W4, W5 and W6 are controlled by the Intertropical Convergence Zone (ITCZ) and its meridional oscillation. The most energetic waves are from the S, generated by the strong winds associated to the passage of cold fronts, which mainly affect the southern region. Wave power presents a decrease in energy levels from south to north, with its annual variation showing that the winter months are the most energetic in W1 to W4, while in W5 and W6 the most energetic conditions occur during the austral summer. The information presented here provides boundary conditions for studies related to coastal processes, fundamental for a better understanding of the Brazilian coastal zone.

Variability and trends of downward surface global solar radiation over the Iberian Peninsula based on ERA-40 reanalysis

A climate study of the incidence of downward surface global solar radiation (SSRD) in the Iberian Peninsula (IP) based primarily on ERA-40 reanalysis is presented. NCEP/NCAR reanalysis and ground-based records from several Portuguese and Spanish stations have been also considered. The results showthat reanalysis can capture a similar inter-annual variability as compared to ground-based observations, especially on a monthly basis, even though annual ERA-40 (NCEP/NCAR) values tend to underestimate (overestimate) the observations with a mean relative difference of around 20Wm–2 (40Wm–2). On the other hand, ground-based measurements in Portuguese stations during the period 1964–1989 show a tendency to decrease until the mid-1970s followed by an increase up to the end of the study period, in line with the dimming/brightening phenomenon reported in the literature. Nevertheless, there are different temporal behaviours as a greater increase since the 1970s is observed in the south and less industrialized regions. Similarly, the ERA-40 reanalysis shows a noticeable decrease until the early 1970s followed by a slight increase up to the end of the 1990s, suggesting a dimming/brightening transition around the early 1970s, earlier in the south and centre and later in the north of the IP. Although there are slight differences in the magnitude of the trends as well as the turning year of the dimming/brightening periods, the decadal changes of ERA-40 fairly agree with the ground-based observations in Portugal and Spain, in contrast to most of the literature for other regions of the world, and is used in the climatology of the SSRD in the study area. NCEP/NCAR reanalysis does not capture the decadal variations of SSRD in the IP. The results show that part of the decadal variability of the global radiation in the IP is related to changes in cloud cover (represented in ERA-40).

A reanalysis of 409 European-Ancestry and African American schizophrenia pedigrees reveals significant linkage to 8p23.3 with evidence of locus heterogeneity

The detection and replication of schizophrenia risk loci can require substantial sample sizes, which has prompted various collaborative efforts for combining multiple samples. However, pooled samples may comprise sub-samples with substantial population genetic differences, including allele frequency differences. We investigated the impact of population differences via linkage reanalysis of Molecular Genetics of Schizophrenia 1 (MGS1) affected sibling-pair data, comprising two samples of distinct ancestral origin: European (EA: 263 pedigrees) and African-American (AA: 146 pedigrees). To exploit the linkage information contained within these distinct continental samples, we performed separate analyses of the individual samples, allowing for within-sample locus heterogeneity, and the pooled sample, allowing for both within-sample and between-sample heterogeneity. Significance levels, corrected for the multiple tests, were determined empirically. For all suggestive peaks, stronger linkage evidence was obtained in either the EA or AA sample than the combined sample, regardless of how heterogeneity was modeled for the latter. Notably, we report genomewide significant linkage of schizophrenia to 8p23.3 and evidence for a second, independent susceptibility locus, reaching suggestive linkage, 29 cM away on 8p21.3. We also detected suggestive linkage on chromosomes 5p13.3 and 7q36.2. Many regions showed pronounced differences in the extent of linkage between the EA and AA samples. This reanalysis highlights the potential impact of population differences upon linkage evidence in pooled data and demonstrates a useful approach for the analysis of samples drawn from distinct continental groups.

Coupled ocean-atmosphere inter-decadal modes in the tropics

Homogencous upper air data for 50 years (1949-1998) from the National Centers for Environmental Prediction/National Center for Atmospheric Research (NCEP/NCAR) reanalysis project, sea surface temperatures and sea level pressure are used to bring out the three dimensional structure of two dominant decadal/multi-decadal variations in the tropics. The global three dimensional modes represent generalized forms of inter-decadal modes studied earlier only with surface data. In the vertical, both modes show approximate first baroclinic structures over the tropics. The Walker circulation associated with the multidecadal mode has a wavenumber two structure in the zonal direction. It is shown that the magnitude of major ascending and descending motions associated with the multi-decadal Hadley and Walker circulations, are comparable to those associated with the dominant inter-annual mode. Implications of these large scale global circulations associated with the low frequency oscillations in modulating regional climate on a inter-annual time scale are discussed.

Correction for Bias in Downscaling GCM Simulations for Hydrologic Impact Assessment

In a statistical downscaling model, it is important to remove the bias of General Circulations Model (GCM) outputs resulting from various assumptions about the geophysical processes. One conventional method for correcting such bias is standardisation, which is used prior to statistical downscaling to reduce systematic bias in the mean and variances of GCM predictors relative to the observations or National Centre for Environmental Prediction/ National Centre for Atmospheric Research (NCEP/NCAR) reanalysis data. A major drawback of standardisation is that it may reduce the bias in the mean and variance of the predictor variable but it is much harder to accommodate the bias in large-scale patterns of atmospheric circulation in GCMs (e.g. shifts in the dominant storm track relative to observed data) or unrealistic inter-variable relationships. While predicting hydrologic scenarios, such uncorrected bias should be taken care of; otherwise it will propagate in the computations for subsequent years. A statistical method based on equi-probability transformation is applied in this study after downscaling, to remove the bias from the predicted hydrologic variable relative to the observed hydrologic variable for a baseline period. The model is applied in prediction of monsoon stream flow of Mahanadi River in India, from GCM generated large scale climatological data.

Development and Evaluation of Statistical Downscaling Models for Monthly Precipitation

Several statistical downscaling models have been developed in the past couple of decades to assess the hydrologic impacts of climate change by projecting the station-scale hydrological variables from large-scale atmospheric variables simulated by general circulation models (GCMs). This paper presents and compares different statistical downscaling models that use multiple linear regression (MLR), positive coefficient regression (PCR), stepwise regression (SR), and support vector machine (SVM) techniques for estimating monthly rainfall amounts in the state of Florida. Mean sea level pressure, air temperature, geopotential height, specific humidity, U wind, and V wind are used as the explanatory variables/predictors in the downscaling models. Data for these variables are obtained from the National Centers for Environmental Prediction-National Center for Atmospheric Research (NCEP-NCAR) reanalysis dataset and the Canadian Centre for Climate Modelling and Analysis (CCCma) Coupled Global Climate Model, version 3 (CGCM3) GCM simulations. The principal component analysis (PCA) and fuzzy c-means clustering method (FCM) are used as part of downscaling model to reduce the dimensionality of the dataset and identify the clusters in the data, respectively. Evaluation of the performances of the models using different error and statistical measures indicates that the SVM-based model performed better than all the other models in reproducing most monthly rainfall statistics at 18 sites. Output from the third-generation CGCM3 GCM for the A1B scenario was used for future projections. For the projection period 2001-10, MLR was used to relate variables at the GCM and NCEP grid scales. Use of MLR in linking the predictor variables at the GCM and NCEP grid scales yielded better reproduction of monthly rainfall statistics at most of the stations (12 out of 18) compared to those by spatial interpolation technique used in earlier studies.

Equivalent-barotropic definition of tropospheric mean temperature

An equivalent-barotropic (EB) description of the tropospheric temperature field is derived from the geostrophic empirical mode (GEM) in the form of a scalar function Gamma(p, phi), where p is pressure and phi is 300-850-mb thickness. Baroclinic parameter phi plays the role of latitude at each longitudinal section. Compared with traditional Eulerian-mean methods, GEM defines a mean field in baroclinic streamfunction space with a time scale much longer than synoptic variability. It prompts an EB concept that is only based on a baroclinic field. Monthly GEM fields are diagnosed from NCEP-NCAR reanalysis data and account for more than 90% of the tropospheric thermal variance. The circumglobal composite of GEM fields exhibits seasonal, zonal, and hemispheric asymmetries, with larger rms errors occurring in winter and in the Northern Hemisphere (NH). Zonally asymmetric features and planetary deviation from EB are seen in the NH winter GEM. Reconstruction of synoptic sections and correlation analysis reveal that the tropospheric temperature field is EB at the leading order and has a 1-day phase lag behind barotropic variations in extratropical regions.

Role of snow depth in spring of Tibetan Plateau in onset of South China Sea summer monsoon

The role of snow depth of Tibetan Plateau in the onset of South China Sea summer monsoon and the influence of ENSO on snow depth of Tibetan Plateau are investigated with use of data from ECMWF reanalysis and NCEP/NCAR reanalysis. The results are as follows: (1) The snow depth data from ECMWF reanalysis are tested and reliable, and can be used to study the influence of snow depth of Tibetan Plateau on the onset of South China Sea summer monsoon; (2) Anomaly of snow depth of Tibetan Plateau causes anomaly in air temperature and its contrast between the Indian Ocean and the continent resulting in easterly wind anomaly over 500 hPa and hence as well as in the atmospheric circulation in the lower layer. For the year of negative anomaly of snow depth a westerly wind anomaly with a cyclone pair takes place, while for positive anomaly of snow depth an easterly anomaly occurs with an anticyclone pair; (3) While positive anomaly of SST occurs in the eastern Pacific Ocean, positive anomaly of air pressure also takes place over the eastern Indian Ocean and the South China Sea, causing stronger meridional pressure gradient between the ocean and continent and then westerly wind anomaly. At the same time, the atmospheric pressure increases in the northern Tibetan Plateau, northerly wind gets stronger, and subtropical front strengthens. All of these are favorable for snowfall over Tibetan Plateau.