Climate change scenarios for surface temperature in Emilia-Romagna (Italy) obtained using statistical downscaling models

Summary Possible changes of mean climate and the frequency of extreme temperature events in Emilia-Romagna, over the period 2070–2100 compared to 1960–1990, are assessed. A statistical downscaling technique, applied to HadAM3P experiments (control, A2 and B2 scenarios) performed at the Hadley Centre, is used to achieve this objective. The method applied consists of a multivariate regression based on Canonical Correlation Analysis (CCA), using as possible predictors mean sea level pressure (MSLP), geopotential height at 500 hPa (Z500) and temperature at 850 hPa (T850), and as predictands the seasonal mean values of minimum and maximum surface temperature (T_min and T_max), 90^th percentile of maximum temperature (T_max90), 10^th percentile of minimum temperature (T_min10), number of frost days (T_nfd) and heat wave duration (HWD) at the station level. First, the statistical model is optimised and calibrated using NCEP/NCAR reanalysis to evaluate the large-scale predictors. The observational data at 32 stations uniformly distributed over Emilia-Romagna are used to compute the local predictands. The results of the optimisation procedure reveal that T850 is the best predictor in most cases, and in combination with MSLP, is an optimum predictor for winter T_max90 and autumn T_min10. Finally, MSLP is the best predictor for spring T_min while Z500 is the best predictor for spring T_max90 and heat wave duration index, except during autumn. The ability of HadAM3P to simulate the present day spatial and temporal variability of the chosen predictors is tested using the control experiments. Finally, the downscaling model is applied to all model output experiments to obtain simulated present day and A2 and B2 scenario results at the local scale. Results show that significant increases can be expected to occur under scenario conditions in both maximum and minimum temperature, associated with a decrease in the number of frost days and with an increase in the heat wave duration index. The magnitude of the change is more significant for the A2 scenario than for the B2 scenario.     

Variability of the Linke and Unsworth-Monteith turbidity parameters in Athens, Greece

Summary  This work studies the variability of the Linke (T _L ) and Unsworth-Monteith (T _U ) turbidity parameters in the urban atmosphere of Athens. Beam irradiance observations performed at the Actinometric Station of the National Observatory of Athens are used in the period 1975–1995. This study examines (i) the inter-annual variation of T _L and T _U , (ii) their mean seasonal variability, (iii) their monthly average variation, and (iv) their mean daily variation. Also, for various air-mass origins, the mean seasonal variations of T _L and T _U are given. The frequency of occurrence of the parameters in various ranges is also shown. An innovative interpretation of the impact of the Saharan dust on the radiative properties in the atmosphere of Athens is attempted. Received May 19, 2000 Revised October 17, 2000     

Monin–Obukhov Similarity Functions of the Structure Parameter of Temperature and Turbulent Kinetic Energy Dissipation Rate in the Stable Boundary Layer

The Monin–Obukhov similarity theory (MOST) functions f_ε and f_T, of the dissipation rate of turbulent kinetic energy (TKE). ε, and the structure parameter of temperature, C_T², were determined for the stable atmospheric surface layer using data gathered in the context of CASES-99. These data cover a relatively wide stability range, i.e. ζ=z/L of up to 10, where z is the height and L the Obukhov length. The best fits were given by f_ε = 0.8 + 2.5ζ and f_T= 4.7[ 1+1.6(ζ)^2/3], which differ somewhat from previously published functions. ε was obtained from spectra of the longitudinal wind velocity using a time series model (ARMA) method instead of the traditional Fourier transform. The neutral limit f_ε =0.8 implies that there is an imbalance between TKE production and dissipation in the simplified TKE budget equation. Similarly, we found a production-dissipation imbalance for the temperature fluctuation budget equation. Correcting for the production-dissipation imbalance, the ‘standard’ MOST functions for dimensionless wind speed and temperature gradients (φ_m and φ_m) were determined from f_ε and f_T and compared with the φ_m and φ_h formulations of Businger and others. We found good agreement with the Beljaars and Holtslag [J. Appl. Meteorol. 30, 327–341 (1991)] relations. Lastly, the flux and gradient Richardson numbers are discussed also in terms of f_ε and f_T.     

基于低频振荡信号的中国南方冬半年持续性低温指数延伸期预报试验

利用1961-2009年36°N以南、108°E以东中国大陆191个站点逐日最低气温和NCEP/NCAR再分析日平均格点,研究与区域持续性低温事件有关的大气低频振荡信号,寻找可以在一定程度上表征不同类型区域持续性低温事件的指数,并尝试结合DERF2.0系统的预报产品进行持续性低温指数的延伸期预报试验。结果表明:（1）在研究范围内的区域持续性低温事件可以归纳为江北型、江南型和全区域型3类,其中江北型和江南型事件的环流背景差异体现在异常环流中心的纬度位置上,而全区域型事件属于增强型的江北型事件;（2）江北型和江南型区域平均最低气温时间序列的10-30 d低频分量的位相和强度变化与区域持续性低温事件的发生有显著关系,可以作为表征区域持续性低温事件指数和预报量;（3）100°-120°E范围内850 hPa温度场距平的经验正交函数分解前两个主模态具有显著的10-30 d变化周期,并且其空间结构分别与江北型和江南型事件的典型环流特征一致,前两个主模态时间系数能够作为持续性低温指数的预报因子;（4）检验结果表明,DERF2.0系统对上述预报因子有一定的预报能力。在延伸期预报时效内,利用统计学和动力学相结合的方法制作的持续性低温指数的预报效果好于模式直接预报的2 m气温,该预报方法有助于提升区域持续性低温事件的延伸预报能力。      

Humidity Effect and Its Influence on the Seasonal Distribution of Precipitation δ18O in Monsoon Regions

The humidity effect, namely the markedly positive correlation between the stable isotopic ratio in precipitation and the dew-point deficit △Td in the atmosphere, is put forward firstly and the relationships between the δ18O in precipitation and △Td are analyzed for the Urumqi and Kunming stations, which have completely different climatic characteristics. Although the seasonal variations in δ18O and △Td exhibit differences between the two stations, their humidity effect is notable. The correlation coefficient and its confidence level of the humidity effect are higher than those of the amount effect at Kunming, showing the marked influence of the humidity conditions in the atmosphere on stable isotopes in precipitation.Using a kinetic model for stable isotopic fractionation, and according to the seasonal distribution of meanmonthly temperature at 500 hPa at Kunming, the variations of the δ18O in condensate in cloud aresimulated. A very good agreement between the seasonal variations of the simulated mean δ18O and themean monthly temperature at 500 hPa is obtained, showing that the oxygen stable isotope in condensateof cloud experiences a temperature effect. Such a result is markedly different from the amount effect atthe ground. Based on the simulations of seasonal variations of δ18O in falling raindrops, it can be foundthat, in the dry season from November to April, the increasing trend with falling distance of δ18O in fallingraindrops corresponds remarkably to the great ATd, showing a strong evaporation enrichment function infalling raindrops; however, in the wet season from May to October, the δ18O in falling raindrops displaysan unapparent increase corresponding to the small ATd, except in May. By comparing the simulated meanδ18O at the ground with the actual monthly δ18O in precipitation, we see distinctly that the two monthlyδ18O variations agree very well. On average, the δ18O values are relatively lower because of the highlymoist air, heavy rainfall, small △Td and weak evaporation enrichment function of stable isotopes in thefalling raindrops, under the influence of vapor from the oceans; but they are relatively higher because of the dry air, light rainfall, great △Td and strong evaporation enrichment function in falling raindrops, under the control of the continental air mass. Therefore, the δ18O in precipitation at Kunming can be used to indicate the humidity situation in the atmosphere to a certain degree, and thus indicate the intensity of the precipitation and the strength of the monsoon indirectly. The humidity effect changes not only the magnitude of the stable isotopic ratio in precipitation but also its seasonal distribution due to its influence on the strength of the evaporation enrichment of stable isotopes in falling raindrops and the direction of the net mass transfer of stable isotopes between the atmosphere and the raindrops. Consequently, it is inferred that the humidity effect is probably one of the foremost causes generating the amount effect.     

热带太平洋表面水中CO2对El Niño事件响应的数值模拟

文中用一个带生物泵的三维全球海洋碳循环模式模拟了热带太平洋表面水中CO2总量(TCO2)在ElNi     

SHEBA flux–profile relationships in the stable atmospheric boundary layer

Measurements of atmospheric turbulence made during the Surface Heat Budget of the Arctic Ocean Experiment (SHEBA) are used to examine the profile stability functions of momentum, φ _m , and sensible heat, φ _h , in the stably stratified boundary layer over the Arctic pack ice. Turbulent fluxes and mean meteorological data that cover different surface conditions and a wide range of stability conditions were continuously measured and reported hourly at five levels on a 20-m main tower for 11 months. The comprehensive dataset collected during SHEBA allows studying φ _m and φ _h in detail and includes ample data for the very stable case. New parameterizations for φ _m (ζ) and φ _h (ζ) in stable conditions are proposed to describe the SHEBA data; these cover the entire range of the stability parameter ζ = z/L from neutral to very stable conditions, where L is the Obukhov length and z is the measurement height. In the limit of very strong stability, φ _m follows a ζ ^1/3 dependence, whereas φ _h initially increases with increasing ζ, reaches a maximum at ζ ≈ 10, and then tends to level off with increasing ζ. The effects of self-correlation, which occur in plots of φ _m and φ _h versus ζ, are reduced by using an independent bin-averaging method instead of conventional averaging.     

Comparative tests and analyses on monsoon and regional precipitation

Summary Due to the strong variation of the Asian monsoon, many countries in Asia often suffer from serious natural disasters. Droughts and floods appear in East China frequently related to the large anomalies of the two branches of East Asian monsoon. Based on rainfall data recorded by 336 Chinese stations in the 1980s, two distinctly opposite rainfall types over East China in summer (JJA) are discovered. Correspondly, 850 hPa anomalous wind fields in the Eastern Hemisphere are also possessed by two types of converse patterns in spring (MAM) and summer (JJA). The onset time and intensity of the Somali jet and the two branches of Southeast and Indian monsoon are quite different. Furthermore, in the 500 hPa geopotential height anomaly fields in spring and summer, the variations of the previous general atmospheric circulation (in spring) are closely correlated to the two kinds of conversely distributed rainfall in summer. These two types of rainfall are also related to two types of conversely distributed sea surface temperature anomaly (SSTA) in the equatorial Pacific and Indian Oceans. To investigate model capabilities refelecting the above observed features, eight numerical experiments are carried out using the IAP 2-L AGCM, with observed monthly mean global SSTs as external forcing and observed atmospheric data on February 15 as initial conditions. The simulated distributions of rainfall anomalies over East China in summer are in good accordance with observations. With conversely distributed SSTAs in the equatorial Pacific, the simulated 850 hPa anomalous wind fields and the 500 hPa geopotential height anomalies are also conversely distributed, and are closely related to the two types of simulated rainfall anomalies over East China. The cross equatorial wind varies in strength, space and time. The simulated distributions of anomalous 500 hPa geopotential height in spring and the anomalous wind at 850 hPa in spring and summer are quite similar to observations.With 10 Figures     

Skill as a function of time scale in ensembles of seasonal hindcasts

Forecast skill as a function of time lead and time averaging is examined in two 6-member ensembles of seasonal hindcasts. One ensemble is produced with the second generation general circulation model of the Canadian Centre for Climate Modelling and Analysis (GCM2) and the other with a reduced resolution version of the numerical weather prediction model of the Canadian Meteorological Centre (SEF). The integrations are initiated from the NCEP/NCAR reanalyzed data. Monthly sea surface temperature anomalies observed prior to the forecast period are maintained throughout the forecast season. A statistical forecast improvement technique, based on the singular value decomposition of forecast and reanalyzed fields, is discussed and evaluated. A simple analogue of the hindcast integrations is used to examine the behavior of two common skill scores, the correlation skill score and the explained variance skill score. The maximal skill score and the corresponding optimal forecast in this analogue are identified. The total skill of the optimal forecast is a sum of two terms, one associated with the initial conditions and the other with the lower boundary forcing. The two sources of skill operate on different time scales, with initial conditions being more important in the first one-two weeks and the atmospheric response to the boundary forcing becoming more dominant for longer time leads and time averages. This suggests that these sources of skill should be considered separately in forecast optimization. The statistical technique is moderately successful in improving the skill of monthly to seasonal forecasts of 500 hPa height (Z ₅₀₀) and 700 hPa temperature (T ₇₀₀) in the Northern Hemisphere and in the North Pacific/North America sector. The improvement is better when the forecasts for the first week and for the rest of the season are optimized separately. The SEF model produces better Z ₅₀₀ and T ₇₀₀ forecasts than GCM2 in the first one-two weeks whereas GCM2 performs slightly better at longer time leads. The skill of zero time lead forecast decays rapidly with averaging interval for time averages up to about 30–45 days and stabilizes, or even rises, for longer time averages. Excluding the first week from seasonal forecasts results in substantial degradation of predictive skill. Received: 1 November 1999 / Accepted: 24 May 2000     

Towards the detection and attribution of an anthropogenic effect on climate

It has been hypothesized recently that regional-scale cooling caused by anthropogenic sulfate aerosols may be partially obscuring a warming signal associated with changes in greenhouse gas concentrations. Here we use results from model experiments in which sulfate and carbon dioxide have been varied individually and in combination in order to test this hypothesis. We use centered [R (t)] and uncentered [C (t)] pattern similarity statistics to compare observed time-evolving surface temperature change patterns with the model-predicted equilibrium signal patterns. We show that in most cases, the C (t) statistic reduces to a measure of observed global-mean temperature changes, and is of limited use in attributing observed climate changes to a specific causal mechanism. We therefore focus on R (t), which is a more useful statistic for discriminating between forcing mechanisms with different pattern signatures but similar rates of global mean change. Our results indicate that over the last 50 years, the summer (JJA) and fall (SON) observed patterns of near-surface temperature change show increasing similarity to the model-simulated response to combined sulfate aerosol/CO₂ forcing. At least some of this increasing spatial congruence occurs in areas where the real world has cooled. To assess the significance of the most recent trends in R (t) and C (t), we use data from multi-century control integrations performed with two different coupled atmosphere-ocean models, which provide information on the statistical behavior of 'unforced' trends in the pattern correlation statistics. For the combined sulfate aerosol/CO₂ experiment, the 50-year R (t) trends for the JJA and SON signals are highly significant. Results are robust in that they do not depend on the choice of control run used to estimate natural variability noise properties. The R (t) trends for the CO₂-only signal are not significant in any season. C (t) trends for signals from both the CO₂-only and combined forcing experiments are highly significant in all seasons and for all trend lengths (except for trends over the last 10 years), indicating large global-mean changes relative to the two natural variability estimates used here. The caveats regarding the signals and natural variability noise which form the basis of this study are numerous. Nevertheless, we have provided first evidence that both the largest-scale (global-mean) and smaller-scale (spatial anomalies about the global mean) components of a combined CO₂/anthropogenic sulfate aerosol signal are identifiable in the observed near-surface air temperature data. If the coupled-model noise estimates used here are realistic, we can be highly confident that the anthropogenic signal that we have identified is distinctly different from internally generated natural variability noise. The fact that we have been able to detect the detailed spatial signature in response to combined CO₂ and sulfate aerosol forcing, but not in response to CO₂ forcing alone, suggests that some of the regional-scale background noise (against which we were trying to detect a CO₂-only signal) is in fact part of the signal of a sulfate aerosol effect on climate. The large effect of sulfate aerosols found in this study demonstrates the importance of their inclusion in experiments designed to simulate past and future climate change. Received: 10 November 1994 / Accepted: 19 July 1995     

Predictability of SST forced climate signals in two atmospheric general circulation models

 The predictability of atmospheric responses to global sea surface temperature (SST) anomalies is evaluated using ensemble simulations of two general circulation models (GCMs): the GENESIS version 1.5 (GEN) and the ECMWF cycle 36 (ECM). The integrations incorporate observed SST variations but start from different initial land and atmospheric states. Five GEN 1980–1992 and six ECM 1980–1988 realizations are compared with observations to distinguish predictable SST forced climate signals from internal variability. To facilitate the study, correlation analysis and significance evaluation techniques are developed on the basis of time series permutations. It is found that the annual mean global area with realistic signals is variable dependent and ranges from 3 to 20% in GEN and 6 to 28% in ECM. More than 95% of these signal areas occur between 35 °S–35 °N. Due to the existence of model biases, robust responses, which are independent of initial condition, are identified over broader areas. Both GCMs demonstrate that the sensitivity to initial conditions decreases and the predictability of SST forced responses increases, in order, from 850 hPa zonal wind, outgoing longwave radiation, 200 hPa zonal wind, sea-level pressure to 500 hPa height. The predictable signals are concentrated in the tropical and subtropical Pacific Ocean and are identified with typical El Ni?o/ Southern Oscillation phenomena that occur in response to SST and diabatic heating anomalies over the equatorial central Pacific. ECM is less sensitive to initial conditions and better predicts SST forced climate changes. This results from (1) a more realistic basic climatology, especially of the upper-level wind circulation, that produces more realistic interactions between the mean flow, stationary waves and tropical forcing; (2) a more vigorous hydrologic cycle that amplifies the tropical forcing signals, which can exceed internal variability and be more efficiently transported from the forcing region. Differences between the models and observations are identified. For GEN during El Ni?o, the convection does not carry energy to a sufficiently high altitude, while the spread of the tropospheric warming along the equator is slower and the anomaly magnitude smaller than observed. This impacts model ability to simulate realistic responses over Eurasia and the Indian Ocean. Similar biases exist in the ECM responses. In addition, the relationships between upper and lower tropospheric wind responses to SST forcing are not well reproduced by either model. The identification of these model biases leads to the conclusion that improvements in convective heat and momentum transport parametrizations and basic climate simulations could substantially increase predictive skill. Received: 25 April 1996 / Accepted: 9 December 1996     

A Comparative Study of the Atmospheric Layers below First Lifting Condensation Level for Instantaneous Pre-Monsoon Thunderstorm Occurrence at Agartala (23o30’N, 91o15’E) and Ranchi (23o14’N, 85o14’E) of India

An attempt has been made to investigate the role of vertical wind shear, corrective instability and the thermodynamic parameter (θes - θe) below the first lifting condensation level (FLCL) in the occurrence of instanta-neous premonsoon thunderstorm over Agartala (AGT) and Ranchi (RNC) at 12 GMT Radiosonde data of 1988 have been utilized here. The study has however been confined to 1000 hPa-500 hPa range at most Here the convectively unstable layers with positive vertical wind shear upto 500 hPa have been termed as ‘Fa?vourable Layers’ (FL) and the level at which an initially stable layer turns out to be convectively unstable for the first time has been termed as ‘Transition Level’ (TL). It is observed that the changes in vertical wind shear are positive at TL at the time of occurrence of thunderstorm (TS) and the corresponding change is negative on fair-weather situa?tion Moreover, the 90% confidence interval for (θes - θe) reveals that for AGT the upper layer thermodynamic characteristic is important at the time of occurrence of TS whereas for RNC, the value of (θes - θe) at the surface is much more effective     

The mid-latitude horizontal and vertical structure of the zonally asymmetric intraseasonal and interannual ozone variability in boreal winters

 A singular value decomposition (SVD) is used to calculate SVD-selected fields of ozone and geopotential, which exhibit maximum covariance, from the observed zonally asymmetric total ozone field and that of the three-dimensional geopotential field thus leaving almost purely dynamical induced variations in the remaining ozone field. This procedure was applied to Total Ozone Mapping Spectrometer data (TOMS) and to geopotential values from the National Centers for Environmental Prediction in the boreal mid-latitudes in the winter months of 1979–1992. Intraseasonal variability (December–February) and trend-eliminated interannual winter mean variability of total ozone and geopotential are analyzed. The first four modes of SVD analysis explain more than 70% of the covariance for the intraseasonal variability and more than 80% of that for the interannual variability. The vertical structure of geopotential regression maps reveals a clear wave-1 pattern for modes one and two and a wave-2 pattern for modes three and four. These patterns show differently but generally westward tilted phases and are more complex at heights below 70 hPa. Further a linear transport model of a conservative tracer was applied to each individual geopotential mode found by the SVD analysis in connection with an observed height and latitude dependent zonal mean ozone distribution. The model results of total ozone reproduce the spatial patterns of the SVD-selected total ozone field quite well whereas their magnitudes are variously underestimated. This method allows us to assess the vertical distribution of the contribution of single modes to the total ozone variability. Maximum contributions are found between 150 and 70 hPa. Smaller amplitude maxima are found around 10 hPa, which result from contributions of horizontal advection of ozone alone. These results reflect an expected dynamical link between the variability of the zonally asymmetric parts of geopotential and ozone. Received: 7 November 1997 / Accepted: 10 June 1998     

Global and regional variability in a coupled AOGCM

 The variability of near surface temperature on global and regional spatial scales and interannual time scales from a 1000 year control integration of the Hadley Centre coupled model (HADCM2-CTL) are compared with the observational record of surface temperature. The model succeeds in reproducing the observed patterns of natural variability, with high variability over the northern continents and low variability over much of the tropics. The model global mean variability has similar strength to observed global mean variability on time scales less than 20 years. The warming seen in the historical record is outside the range of natural variability as simulated in HADCM2-CTL. The model has El-Ni?o/Southern Oscillation (ENSO)-like behaviour with a central Pacific, peak to peak, strength of approximately 3 K. Changes in near surface temperature in the central Pacific are strongly correlated with changes in near surface temperature over most of the tropics, large regions of the extra-tropics and changes in tropical ocean upper 250 m heat content. Tropospheric temperature changes and tropical surface pressure changes are also strongly correlated with changes in the central Pacific surface temperature. Oceanic regions show significant departures from an AR1 or first order Markov behaviour in the Northwest Atlantic, Northwest Pacific and Arctic oceans. The Northwest Atlantic region has large amounts of variability over periods greater than 50 years. This variability is associated with a jump in the strength of North Atlantic meridional stream function. The spectra of the Western European and Continental US land regions are not significantly different from an AR1 process. The flow through the Drake Passage has an interannual standard deviation of approximately 2.5 Sv with significant departures from an AR1 process at time scales greater than 40 years. Winter northern hemispheric 500 hPa geopotential height shows some evidence of multiple regimes but no year to year persistence of these regimes. Received: 31 January 1996/Accepted: 22 July 1996     

1998年夏季HUBEX/GAME期间热量和水汽收支(英)

By using the high-resolution GAME reanalysis data, the heat and moisture budgets during the period of HUBEX/GAME in the summer of 1998 are calculated for exploring the thermodynamic features of Meiyu over the Changjiang-Huaihe (CH) valley. During the CH Meiyu period, an intensive vertically-integrated heat source and moisture sink are predominant over the heavy rainfall area of the CH valley, accompanied by strong upward motion at 500 hPa. The heat and moisture budgets show that the main diabatic heating component is condensation latent heat released by rainfall. As residual terms, the evaporation and sensible heating are relatively small. Based on the vertical distribution of the heat source and moisture sink, the nature of the rainfall is mixed, in which the convective rainfall is dominant with a considerable percentage of continuous stratiform rainfall. There are similar time evolutions of the main physical parameters(〈Q1〉,〈Q2〉,and vertical motion ω at 500 hPa).The time variations of〈Q1〉and〈Q2〉are in phase with those of -ω500, and have their main peaks within the CH Meiyu period. This shows the influence of the heat source on the dynamic structure of the atmosphere. The wavelet analyses of those time series display similar multiple timescale characteristics. During the CH Meiyu period, both the synoptic scale(～6 days) and mesoscale (～2 days and ～12 hours) increase obviously and cause heavy rainfall as well as the appearances of the maxima of the main physical parameters. Among them, the mesoscale systems are the main factors.     

Climate change impacts on dryland cropping systems in the Central Great Plains,USA

Agricultural systems models are essential tools to assess potential climate change (CC) impacts on crop production and help guide policy decisions. In this study, impacts of projected CC on dryland crop rotations of wheat-fallow (WF), wheat-corn-fallow (WCF), and wheat-corn-millet (WCM) in the U.S. Central Great Plains (Akron, Colorado) were simulated using the CERES V4.0 crop modules in RZWQM2. The CC scenarios for CO₂, temperature and precipitation were based on a synthesis of Intergovernmental Panel on Climate Change (IPCC 2007) projections for Colorado. The CC for years 2025, 2050, 2075, and 2100 (CC projection years) were super-imposed on measured baseline climate data for 15–17 years collected during the long-term WF and WCF (1992–2008), and WCM (1994–2008) experiments at the location to provide inter-annual variability. For all the CC projection years, a decline in simulated wheat yield and an increase in actual transpiration were observed, but compared to the baseline these changes were not significant (p > 0.05) in all cases but one. However, corn and proso millet yields in all rotations and projection years declined significantly (p < 0.05), which resulted in decreased transpiration. Overall, the projected negative effects of rising temperatures on crop production dominated over any positive impacts of atmospheric CO₂ increases in these dryland cropping systems. Simulated adaptation via changes in planting dates did not mitigate the yield losses of the crops significantly. However, the no-tillage maintained higher wheat yields than the conventional tillage in the WF rotation to year 2075. Possible effects of historical CO₂ increases during the past century (from 300 to 380 ppm) on crop yields were also simulated using 96 years of measured climate data (1912–2008) at the location. On average the CO₂ increase enhanced wheat yields by about 30%, and millet yields by about 17%, with no significant changes in corn yields.     

Decadal time scale variability recorded in the Quelccaya summit ice core δ18O isotopic ratio series and its relation with the sea surface temperature

 The spectral characteristics of the δ¹⁸O isotopic ratio time series of the Quelccaya ice cap summit core are investigated with the multi taper method (MTM), the singular spectrum analysis (SSA) and the wavelet transform (WT) techniques for the 500 y long 1485–1984 period. The most significant (at the 99.8% level) cycle according to the MTM F-test has a period centered at 14.4 y while the largest variance explaining oscillation according to the SSA technique has a period centered at 12.9 y. The stability over time of these periods is investigated by performing evolutive MTM and SSA on the 500 y long δ¹⁸O series with a 100 y wide moving window. It is shown that the cycles with largest amplitude and that the oscillations with largest extracting variance have corresponding periods aggregated around 13.5 y that are very stable over the period between 1485 and 1984. The WT of the same isotopic time series reveals the existence of a main oscillation around 12 y which are also very stable in time. The relation between the isotopic data at Quelccaya and the annual sea surface temperature (SST) field anomalies is then evaluated for the overlapping 1919–1984 period. Significant global correlation and significant coherency at 12.1 y are found between the isotopic series and the annual global sea surface temperature (GSST) series. Moreover, the correlation between the low (over 8 y) frequency component of the isotopic time series and the annual SST field point out significant values in the tropical North Atlantic. This region is characterized by a main SST variability at 12.8 y. The Quelccaya δ¹⁸O isotopic ratio series may therefore be considered as a good recorder of the tropical North Atlantic SSTs. This may be explained by the following mechanism: the water vapor amount evaporated by the tropical North Atlantic is function of the SST. So is the water vapor δ¹⁸O isotopic ratio. This water vapor is advected during the rainy season by northeast winds and precipitates at the Quelccaya summit with its tropical North Atlantic isotopic signature. It is also suggested from this described stability of the decadal time scale variability observed in the Quelccaya isotopic series, that the decadal time scale GSST variability was also stable during the last five centuries. Received: 12 February 1997 / Accepted: 9 September 1997     

ENSO 与中国东部夏季降水的关联

计算1 月减6 月El Niño 3.4 指数与6—8 月平均200、850 hPa 风场的相关矢量,分析中等或强ElNiño/La Niña 事件后的夏季(6—8 月)中国东部降水异常分布、西太平洋副热带高压异常特征。结果表明,对ENSO 的响应,无论高、底层大气环流还是西太平洋副热带高压,1970 年代中期气候突变后变为更敏感。主要表现在:对衰减的El Niño 的响应,夏季南亚高压偏东,西太平洋副热带高压偏强、偏西、偏南,印度季风、南海季风减弱,黄河下游以南副热带季风增强。黄河中下游及以南形成异常环流辐合带,由El Niño 导致的降水正异常最有可能出现在这一西南-东北的带状区域。对衰减的La Niña 响应大致相反。      

Quantitative reconstruction of the last interglacial vegetation and climate based on the pollen record from Lake Baikal, Russia

Changes in mean temperature of the coldest (T _c) and warmest month (T _w), annual precipitation (P _ann) and moisture index (α) were reconstructed from a continuous pollen record from Lake Baikal, Russia. The pollen sequence CON01-603-2 (53°57′N, 108°54′E) was recovered from a 386 m water depth in the Continent Ridge and dated to ca. 130–114.8 ky BP. This time interval covers the complete last interglacial (LI), corresponding to MIS 5e. Results of pollen analysis and pollen-based quantitative biome reconstruction show pronounced changes in the regional vegetation throughout the record. Shrubby tundra covered the area at the beginning of MIS 5e (ca. 130–128 ky), consistent with the end of the Middle Pleistocene glaciation. The late glacial climate was characterised by low winter and summer temperatures (T _c ~ −38 to −35°C and T _w~11–13°C) and low annual precipitation (P _ann~300 mm). However, the wide spread of tundra vegetation suggests rather moist environments associated with low temperatures and evaporation (reconstructed α~1). Tundra was replaced by boreal conifer forest (taiga) by ca. 128 ky BP, suggesting a transition to the interglacial. Taiga-dominant phase lasted until ca. 117.4 ky BP, e.g. about 10 ky. The most favourable climate conditions occurred during the first half of the LI. P _ann reached 500 mm soon after 128 ky BP. However, temperature changed more gradually. Maximum values of T _c ~ −20°C and T _w~16–17°C are reconstructed from about 126 ky BP. Conditions became gradually colder after ca. 121 ky BP. T _c dropped to ~ −27°C and T _w to ~15°C by 119.5 ky BP. The reconstructed increase in continentality was accompanied by a decrease in P _ann to ~400–420 mm. However, the climate was still humid enough (α~0.9) to support growth of boreal evergreen conifers. A sharp turn towards a dry climate is reconstructed after ca. 118 ky BP, causing retreat of forest and spread of cool grass-shrub communities. Cool steppe dominated the vegetation in the area between ca. 117.5 ky and 114.8 ky BP, suggesting the end of the interglacial and transition to the last glacial (MIS 5d). Shift to the new glaciation was characterised by cool and very dry conditions with T _c ~ −28 to −30°C, T _w~14–15°C, P _ann~250 mm and α~0.5.     

2009 年我国汛期降水形势总结与三个常用模式预报效果检验

利用实况24小时降水、形势场资料及T213、T639、Japan模式降水、形势场的预报资料,对2009年汛期(5-9月,下同)中国降水时空分布进行分析,并对T213、T639、Japan三个常用模式对2009年汛期的天气形势、降水及其影响系统的预报做主客观检验,以期得出2009年汛期降水分布特点及三个模式的降水预报效果对比.结果表明:(1) 2009年汛期华南地区降水量为全国之最,长江中下游和西南东部地区其次,东北和华北地区再次.(2)从TS评分看,Japan模式的小雨～大雨量级评分较高,T639模式暴雨～大暴雨量级评分较高；T213模式对华北地区暴雨、大暴雨量级降水预报评分高于Japan和T639模式.(3)从降水预报偏差看,T213模式对华北预报明显偏强,T639模式对华北预报强度较为适中,两模式对其他区域中等以下强度降水预报偏强,对强降水预报偏弱；T639对中等以下强度降水预报偏强程度明显小于T213,而对强降水除华南和东北区域外,预报偏弱程度明显大于T213；Japan模式预报偏差随降水量级增大而减小,对大雨以上各量级预报均明显偏弱,且偏弱程度明显大于T213、 T639.(4)由代表性形势场预报检验结果可知,除T213对500hPa高度场、850hPa温度场预报效果好于其他两模式外,各模式预报效果相差不大.(5)三个模式对500 hPa副高总体预报偏东、偏北、偏强,但Japan预报效果明显好于T213、T639.(6) T639模式对台风和低涡的预报相对较好,T213较差.