Climate change impact assessment on flow regime by incorporating spatial correlation and scenario uncertainty

Flooding risk is increasing in many parts of the world and may worsen under climate change conditions. The accuracy of predicting flooding risk relies on reasonable projection of meteorological data (especially rainfall) at the local scale. The current statistical downscaling approaches face the difficulty of projecting multi-site climate information for future conditions while conserving spatial information. This study presents a combined Long Ashton Research Station Weather Generator (LARS-WG) stochastic weather generator and multi-site rainfall simulator RainSim (CLWRS) approach to investigate flow regimes under future conditions in the Kootenay Watershed, Canada. To understand the uncertainty effect stemming from different scenarios, the climate output is fed into a hydrologic model. The results showed different variation trends of annual peak flows (in 2080–2099) based on different climate change scenarios and demonstrated that the hydrological impact would be driven by the interaction between snowmelt and peak flows. The proposed CLWRS approach is useful where there is a need for projection of potential climate change scenarios.

     

Arctic haze: Perturbation to the radiation field

Summary A model of the polluted arctic troposphere is constructed to estimate the magnitude and seasonal variation of the climate forcing function of arctic haze. Using a pill-box bathtub model for the Arctic and envisioning it to be filled with pollution from industrial sources in Eurasia, we estimate that maximum climate perturbation from arctic contamination occurs in the spring months. The major perturbation to the radiation budget is a lowering of the albedo (heating) of the earth-atmosphere system around the vernal equinox and is due to a trace amount (about 5% by mass) of black carbon associated with the removal-resistant submicron mode of aerosols. The black carbon over the reflecting polar ice/snow introduces a heating of about 1.5 degree per day into the haze layer.With 8 Figures     

Observation of gravity waves in a boreal forest

In this paper we report the results of the analysis of two 60-min wave events that occurred in a boreal aspen forest during the 1994 BOREAS (Boreal Ecosystems-Atmosphere Study) field experiment. High frequency wind and temperature data were provided by three 3-D sonic anemometer/thermometers and fourteen fine-wire thermocouples positioned within and above the forest. Wave phase speeds, estimated from information revealed by spectral analysis and linear plane wave equations, are 2.2 and 1.3 m s^-1 for the two events. The wavelengths are 130 m and 65 m respectively and are much larger than the vertical wave displacements. There is strong evidence from the present analysis and from the literature supporting our postulate that these waves are generated by shear instability. We propose that wind shear near the top of the stand is often large enough to reduce the gradient Richardson number below the critical value of 0.25 and thus is able to trigger the instability. When external conditions are favorable, the instability will grow into waves.     

Multiyear simulation of the African climate using a regional climate model (RegCM3) with the high resolution ERA-interim reanalysis

This study examines the ability of the latest version of the International Centre for Theoretical Physics (ICTP) regional climate model (RegCM3) to reproduce seasonal mean climatologies, annual cycle and interannual variability over the entire African continent and different climate subregions. The new European Center for Medium Range Weather Forecast (ECMWF) ERA-interim reanalysis is used to provide initial and lateral boundary conditions for the RegCM3 simulation. Seasonal mean values of zonal wind profile, temperature, precipitation and associated low level circulations are shown to be realistically simulated, although the regional model still shows some deficiencies. The West Africa monsoon flow is somewhat overestimated and the Africa Easterly Jet (AEJ) core intensity is underestimated. Despite these biases, there is a marked improvement in these simulated model variables compared to previous applications of this model over Africa. The mean annual cycle of precipitation, including single and multiple rainy seasons, is well captured over most African subregions, in some cases even improving the quality of the ERA-interim reanalysis. Similarly, the observed precipitation interannual variability is well reproduced by the regional model over most regions, mostly following, and sometimes improving, the quality of the ERA-interim reanalysis. It is assessed that the performance of this model over the entire African domain is of sufficient quality for application to the study of climate change and climate variability over the African continent.     

A Probabilistic Tsunami Hazard Study of the Auckland Region, Part II: Inundation Modelling and Hazard Assessment

Regional source tsunamis pose a potentially devastating hazard to communities and infrastructure on the New Zealand coast. But major events are very uncommon. This dichotomy of infrequent but potentially devastating hazards makes realistic assessment of the risk challenging. Here, we describe a method to determine a probabilistic assessment of the tsunami hazard by regional source tsunamis with an “Average Recurrence Interval” of 2,500-years. The method is applied to the east Auckland region of New Zealand. From an assessment of potential regional tsunamigenic events over 100,000 years, the inundation of the Auckland region from the worst 100 events is modelled using a hydrodynamic model and probabilistic inundation depths on a 2,500-year time scale were determined. Tidal effects on the potential inundation were included by coupling the predicted wave heights with the probability density function of tidal heights at the inundation site. Results show that the more exposed northern section of the east coast and outer islands in the Hauraki Gulf face the greatest hazard from regional tsunamis in the Auckland region. Incorporating tidal effects into predictions of inundation reduced the predicted hazard compared to modelling all the tsunamis arriving at high tide giving a more accurate hazard assessment on the specified time scale. This study presents the first probabilistic analysis of dynamic modelling of tsunami inundation for the New Zealand coast and as such provides the most comprehensive assessment of tsunami inundation of the Auckland region from regional source tsunamis available to date.     

Performance of a multi-RCM ensemble for South Eastern South America

The ability of four regional climate models to reproduce the present-day South American climate is examined with emphasis on La Plata Basin. Models were integrated for the period 1991–2000 with initial and lateral boundary conditions from ERA-40 Reanalysis. The ensemble sea level pressure, maximum and minimum temperatures and precipitation are evaluated in terms of seasonal means and extreme indices based on a percentile approach. Dispersion among the individual models and uncertainties when comparing the ensemble mean with different climatologies are also discussed. The ensemble mean is warmer than the observations in South Eastern South America (SESA), especially for minimum winter temperatures with errors increasing in magnitude towards the tails of the distributions. The ensemble mean reproduces the broad spatial pattern of precipitation, but overestimates the convective precipitation in the tropics and the orographic precipitation along the Andes and over the Brazilian Highlands, and underestimates the precipitation near the monsoon core region. The models overestimate the number of wet days and underestimate the daily intensity of rainfall for both seasons suggesting a premature triggering of convection. The skill of models to simulate the intensity of convective precipitation in summer in SESA and the variability associated with heavy precipitation events (the upper quartile daily precipitation) is far from satisfactory. Owing to the sparseness of the observing network, ensemble and observations uncertainties in seasonal means are comparable for some regions and seasons.     

Numerical integration of a linear barotropic model using three methods of treating meteorological and gravitational modes separately

Summary The analytical solution of a linear barotropic model is derived, including details of the quasi-geostrophic initialization procedure. The prognostic equations are integrated using three different methods of treating the meteorological and gravitational modes separately. These are a semi-Eulerian, semi-implicit (EI) technique, a semi-Lagrangian, semi-implicit (LI) procedure, and a split-explicit (SE) method. The stability criteria and phase speeds are derived for each of the three techniques.The following theoretical conclusions are derived. Of course, in actual numerical integrations particularly those using more complex models, the results are not so unequivocal.The stability of the EI procedure is governed by the CFL criterion for the meteorological mode. Gravity waves have no effect on the timestep but move more slowly than the analytical waves. The LI method is unconditionally stable with respect to both meteorological and gravitational modes. There is thus no timestep restriction. However, the gravity waves have the same reduced phase speed as in the EI technique. The SE procedure has CFL timestep criteria for both the meteorological and gravitational calculations. However, its gravity wave phase speeds are relatively accurate. Moreover, it is the only one of the three methods that handles the nearly-compensating pressure gradient and Coriolis forces together. From the point of computational efficiency, the LI technique is probably the best.     

A New Method of Significance Testing for Correlation-Coefficient Fields and Its Application

Correlation-coefficient fields are widely used in short-term climate prediction research. The most frequently used significance test method for the correlation-coefficient field was proposed by Livezey, in which the number of significantcorrelation lattice(station) points on the correlation coherence map is used as the statistic. However, the method is based on two assumptions:(1) the spatial distribution of the lattice(station) points is uniform;and(2) there is no correlation between the physical quantities in the correlation-coefficient field. However, in reality, the above two assumptions are not valid.Therefore, we designed a more reasonable method for significance testing of the correlation-coefficient field. Specifically, a new statistic, the significant-correlation area, is introduced to eliminate the inhomogeneity of the grid(station)-point distribution, and an empirical Monte Carlo method is employed to eliminate the spatial correlation of the matrix.Subsequently, the new significance test was used for simultaneous correlation-coefficient fields between intensities of the atmospheric activity center in the Northern Hemisphere and temperature/precipitation in China. The results show that the new method is more reasonable than the Livezey method.     

A study of two mass-consistent models: Problems and possible solutions

Two existing mass-consistent models, COMPLEX and NOABL, were tested in three regions in the UK: Devon, the Northern Pennines and Shetland. In order to solve problems arising with both models, a number of modifications were made, leading to the development of a new model, MC-3. First, the entire wind field was scaled down (or up) by a certain percentage to make the wind predictions at the predictor stations as close to the observed values as possible. Second, modifications to the non-divergent process used in COMPLEX were made. Site roughness-length and anemometer height corrections were also made. The new model worked well in terms of mean wind speed predictions, and gave a significant improvement over the predictions of COMPLEX and NOABL.