The impact of climate change on the river rhine: A scenario study

This paper concerns the impact of human-induced global climate change on the River Rhine discharge. For this purpose a model for climate assessment, named ESCAPE, is coupled to a water balance model, named RHINEFLOW. From climate scenarios, changes in regional annual water availability and seasonal discharge in the River Rhine Basin are estimated. The climate scenarios are based on greenhouse gases emissions scenarios. An assessment is made for best guess seasonal discharge changes and for changes in frequencies of low and high discharges in the downstream reaches of the river. In addition, a quantitative estimation of the uncertainties associated with this guess is arrived at.The results show that the extent and range of uncertainty is large with respect to the best guess changes. The uncertainty range is 2–3 times larger for the Business-as-Usual than for the Accelerated Policies scenarios. This large range stems from the doubtful precipitation simulations from the present General Circulation Models. This scenario study showed the precipitation scenarios to be the key-elements within the present range of reliable climate change scenarios.For the River Rhine best guess changes for annual water availability are small according to both scenarios. The river changes from a present combined snow-melt-rain fed river to an almost entirely rain fed river. The difference between present-day large average discharge in winter and the small average discharge in autumn should increase for all scenarios. This trend is largest in the Alpine part of the basin. Here, winter discharges should increase even for scenarios forecasting annual precipitation decreases. Summer discharge should decrease. Best guess scenarios should lead to increased frequencies of both low and high flow events in the downstream (Dutch) part of the river. The results indicate changes could be larger than presently assumed in worst case scenarios used by the Dutch water management authorities.     

Stable boundary layers: Observations,models and variability part I: Modelling and measurements

Selected field measurements of evening stable boundary layers are presented in detail comparable with published Large Eddy Simulation results. Such models appear to match idealized theories more closely than do some boundary-layer observations. Any attempt to compare detailed observations with idealized models therefore highlights the variability of the real boundary layer.Here direct turbulence measurements across the stable boundary layer from a heterogeneous and an ideal site are contrasted. Recommendations are made for the information needed to distinguish heterogeneous and ideal cases.The companion paper (Part II) discusses further the issues of data, analysis in the presence of variability, and the effects of averaging over heterogeneous terrain.Part of UK Meteorological Office Atmospheric Process Research Division.     

Stable boundary layers: observations,models and variability part II: Data analysis and averaging effects

As argued in Part I (Derbyshire, 1995), variability is a key issue in stable boundary layers, and differences in variability between observations and idealized models may imply sytematic biases. Here we discuss how data analysis can be geared to allow for variability and thus consistency with models. Instrumental errors, smoothing methods and vertical discretization are considered. We then show how statistical averaging broadly improves the agreement of heterogeneous results in Part I with the Brost-Wyngaard closure. Recommendations are made for the information needed to analyze apparent differences between homogeneous and heterogeneous stable boundary layers.Part of UK Meteorological Office Atmospheric Processes Research Division     

Scavenging Efficiency of ‘Aerosol Carbon’ and Sulfate in Supercooled Clouds at Mt. Sonnblick (3106 m a.s.l., Austria)

Cloud water and interstitial aerosol samples collected at Mt. Sonnblick (SBO) were analyzed for sulfate and aerosol carbon to calculate in-cloud scavenging efficiencies. Scavenging efficiencies for sulfate (_SO) ranged from 0.52 to 0.99 with an average of 0.80. Aerosol carbon was scavenged less efficiently with an average value (_AC) of 0.45 and minimum and maximum values of 0.14 and 0.81, respectively. Both _SO and _AC showed a marked, but slightly different, dependence on the liquid water content (LWC) of the cloud. At low LWC, _SO increased with rising LWC until it reached a relatively constant value of 0.83 above an LWC of 0.3 g/m³. In the case of aerosol carbon, we obtained a more gradual increase of _AC up to an LWC of 0.5 g/m³. At higher LWCs, _{_} remained relatively constant at 0.60. As the differences between _SO and _A varied across the LWC range observed at SBO, we assume that part of the aerosol carbon was incorporated into the cloud droplets independently from sulfate. This hypothesis is supported by size classified aerosol measurements. The differences in the size distributions of sulfate and total carbon point to a partially external mixture. Thus, the different chemical nature and the differences in the size and mixing state of the aerosol particles are the most likely candidates for the differences in the scavenging behavior.     

AN E – ε – ℓ TURBULENCE CLOSURE SCHEME FOR PLANETARY BOUNDARY-LAYER MODELS: THE NEUTRALLY STRATIFIED CASE

The standard E – model generates aplanetary boundary layerthat appears to be much too deep. The cause of theproblem is traced to the equation for the dissipationrate () of turbulent kinetic energy (E), specifically theparameterization of dissipation production anddestruction. In the context of atmosphericboundary-layer modelling, we argue that a part of thedissipation production should be modelled as the inputto the spectral cascade from the energy-containingpart of the spectrum, with a characteristic length , while the equilibrium imbalancebetween local production and destruction ofdissipation is modelled as proportional toE2/E, as in the standard model. Wepropose an E – – turbulence closurescheme, in which both the mixing length, m, and are prescribed. The importance ofthe equation is diminished, though itstill determines the dissipation rate in the Eequation.     

Measurements of solar radiation and estimation of optical depth in the high Andes of Peru

Summary During an expedition to the high Andes of Southern Peru in June–July 1977, measurements of direct solar radiation in four spectral bands (0.270–0.530–0.630–0.695–2.900 ) were conducted at six sites in elevations ranging from sea level to 5645 m. These measurements were evaluated in Langley plots to determine total optical depths () and irradiances at the top of the atmosphere. In addition, water vapor optical depths (_wv) were calculated from the mean radiosounding over Lima during the expedition, and Rayleigh (_ray) and ozone (_oz) optical depths were obtained from published tabulations. Subtracting _ray, _oz, and _wv from yielded estimates of aerosol optical depth _aer. The components _ray and _oz decrease from the shorter towards the longer wavelength bands and from the lower towards the higher elevation sites; _aer also decreases towards the higher elevations. Particularly pronounced is the decrease of _aer and from the lowlands of the Pacific coast to the highlands of the interior, reflecting the effect of a persistent lower-tropospheric inversion and the contrast from the marine boundary layer to the clear atmosphere of the high Andes.With 4 Figures     

On the Response of the Greenland Ice Sheet to Greenhouse Climate Change

Numerical computations are performed with the three-dimensional polythermal ice-sheet model SICOPOLIS in order to investigate the possible impact of a greenhouse-gas-induced climate change on the Greenland ice sheet. The assumed increase of the mean annual air temperature above the ice covers a range from T = 1°C to 12°C, and several parameterizations for the snowfall and the surface melting are considered. The simulated shrinking of the ice sheet is a smooth function of the temperature rise, indications for the existence of critical thresholds of the climate input are not found. Within 1000 model years, the ice-volume decrease is limited to 10% of the present volume for T 3°C, whereas the most extreme scenario, T = 12°C, leads to an almost entire disintegration, which corresponds to a sea-level equivalent of 7 m. The different snowfall and melting parameterizations yield an uncertainty range of up to 20% of the present ice volume after 1000 model years.     

A Green Planet Versus a Desert World: Estimating the Maximum Effect of Vegetation on the Land Surface Climate

We quantify the maximum possible influence of vegetation on the global climate by conducting two extreme climate model simulations: in a first simulation (desert world), values representative of a desert are used for the land surface parameters for all non glaciated land regions. At the other extreme, a second simulation is performed (green planet) in which values are used which are most beneficial for the biosphere's productivity. Land surface evapotranspiration more than triples in the presence of the green planet, land precipitation doubles (as a second order effect) and near surface temperatures are lower by as much as 8 K in the seasonal mean resulting from the increase in latent heat flux. The differences can be understood in terms of more absorbed radiation at the surface and increased recycling of water. Most of the increase in net surface radiation originates from less thermal radiative loss and not from increases in solar radiation which would be expected from the albedo change. To illustrate the differences in climatic character and what it would imply for the vegetation type, we use the Köppen climate classification. Both cases lead to similar classifications in the extra tropics and South America indicating that the character of the climate is not substantially altered in these regions. Fundamental changes occur over Africa, South Asia and Australia, where large regions are classified as arid (grassland/desert) climate in the desert world simulation while classified as a forest climate in the green planet simulation as a result of the strong influence of maximum vegetation on the climate. This implies that these regions are especially sensitive to biosphere-atmosphere interaction.     

Nitrate in the atmospheric boundary layer of the tropical South Pacific: Implications regarding sources and transport

Weekly bulk aerosol samples collected at Funafuti, Tuvalu (8°30S, 179°12E), American Samoa (14°15S, 170°35W), and Rarotonga (21°15S, 159°45W), from 1983 through most of 1987 have been analyzed for nitrate and other constituents. The mean nitrate concentration is about 0.11 g m^–3 at each of these stations: 0.107±0.011 g m^–3 at Funafuti; 0.116±0.008 at American Samoa; and 0.117±0.010 at Rarotonga. Previous measurements of mineral aerosol and trace metal concentrations at American Samoa are among the lowest ever recorded for the near-surface troposphere and indicate that this region is minimally affected by transport of soil material and pollutants from the continents. Consequently, the nitrate concentration of 0.11 g m^–3 can be regarded as the natural level for the remote marine boundary layer of the tropical South Pacific Ocean. In contrast, over the tropical North Pacific which is significantly impacted by the transport of material from Asia and North America, the mean nitrate concentrations are about three times higher, 0.29 and 0.36 g m^–3 at Midway and Oahu, respectively. The major sources of the nitrate over the tropical South Pacific are still very uncertain. A very significant correlation between the nitrate concentrations at American Samoa and the concentrations of ²¹⁰Pb suggests that transport from continental sources might be important. This continental source could be lightning, which occurs most frequently over the tropical continents. A near-zero correlation with ⁷Be indicates that the stratosphere and upper troposphere are probably not the major sources. A significant biogenic source would be consistent with the higher mean nitrate concentrations, 0.16 to 0.17 g m^–3, found over the equatorial Pacific at Fanning Island (3°55N, 159°20W) and Nauru (0°32S, 166°57E). The lack of correlation between nitrate and nss sulfate at American Samoa does not necessarily preclude an important role for marine biogenic sources.     

Vulnerability of Island Tropical Montane Cloud Forests to Climate Change, with Special Reference to East Maui, Hawaii

Island tropical montane cloud forests may be among the most sensitive of the world's ecosystems to global climate change. Measurements in and above a montane cloud forest on East Maui, Hawaii, document steep microclimatic gradients. Relatively small climate-driven shifts in patterns of atmospheric circulation are likely to trigger major local changes in rainfall, cloud cover, and humidity. Increased interannual variability in precipitation and hurricane incidence would provide additional stresses on island biota that are highly vulnerable to disturbance-related invasion of non-native species. Because of the exceptional sensitivity of these microclimates and forests to change, they may provide valuable listening posts for detecting the onset of human-induced global climate change.     

Adjusting to Policy Expectations in Climate Change Modeling

This paper surveys and interprets the attitudes of scientists to the use of flux adjustments in climate projections with coupled Atmosphere Ocean General Circulation Models. The survey is based largely on the responses of 19 climate modellers to several questions and a discussion document circulated in 1995. We interpret the responses in terms of the following factors: the implicit assumptions which scientists hold about how the environmental policy process deals with scientific uncertainty over human-related global warming; the different scientific styles that exist in climate research; and the influence of organisations, institutions, and policy upon research agendas. We find evidence that scientists' perceptions of the policy process do play a role in shaping their scientific practices. In particular, many of our respondents expressed a preference for keeping discussion of the issue of flux adjustments within the climate modeling community, apparently fearing that climate contrarians would exploit the issue in the public domain. While this may be true, we point to the risk that such an approach may backfire. We also identify assumptions and cultural commitments lying at a deeper level which play at least as important a role as perceptions of the policy process in shaping scientific practices. This leads us to identify two groups of scientists, pragmatists and purists, who have different implicit standards for model adequacy, and correspondingly are or are not willing to use flux adjustments.     

The validation of various schemes for parameterizing evaporation from bare soil for use in meteorological models: A numerical study using in situ data

Parameterization of evaporation from a non-plant-covered surface is very important in the hierarchy strategy of modelling land surface processes. One of the representations frequently used in its computation is the resistance formulation. The performance of the evaporation schemes using the , , and their combination resistance approaches to parameterize evaporation from bare soil surfaces is discussed. For that purpose, the nine schemes, based on a different dependence of and on volumetric soil moisture content and its saturated value, are used.The tests of performances of the considered schemes are based on time integrations by the land surface module (BARESOIL) using observed data. The 23 data sets at a bare surface experimental site in Rimski anevi, Yugoslavia on chernozem soil, were used for the resistance algorithm evaluation. The quality of the schemes was compared with the observed values of the latent heat flux using several statistical parameters.     

Similarity of atmospheric Reynolds shear stress and heat flux fluctuations over a rough surface

Horizontal u and vertical w velocity fluctuations have been measured together with temperature fluctuations in the atmospheric surface layer, at a small height above a wheat crop canopy. Marginal probability density functions are presented for both individual fluctuations u, w, and for the instantaneous Reynolds stress uw, and heat fluxes w and u. Probability density functions of the velocity fluctuations deviate less significantly from the Gaussian form than the probability density of temperature. There appears to be closer similarity between statistics of the instantaneous heat fluxes than between the momentum flux and either of the heat fluxes investigated. The mean momentum flux receives equal contributions from the events referred to as ejections and sweeps in laboratory boundary layers. Sweeps provide the largest contribution to the heat fluxes.     

Shear convection

The paper discusses convection in the presence of wind shear, a condition analysed previously by Zilitinkevich (1971). This region (between the forced and fully developed convective layers) was also considered by Betchov and Yaglom (1971). In the present paper the author endeavors to develop a consistent analysis from the basic hypothesis of a very weak interaction between the vertical convective motions and mechanical turbulence, employing a new similarity model of the turbulent regime. Additional experimental data are introduced. Unlike the notation used in the references quoted above, this regime is termed shear convection rather than free convection. The latter is traditionally regarded as synonymous with the terms pure or fully developed convection.     

Carbon dioxide emissions in a methane economy

Increasing reliance on natural gas (methane) to meet global energy demands holds implications for atmospheric CO₂ concentrations. Analysis of these implications is presented, based on a logistic substitution model viewing energy technologies like biological species invading an econiche and substituting in case of superiority for existing species. This model suggests gas will become the dominant energy source and remain so for 50 years, peaking near 70 percent of world supply. Two scenarios of energy demand are explored, one holding per capita consumption at current levels, the second raising the global average in the year 2100 to the current U.S. level. In the first (efficiency) scenario concentrations peak about 450 ppm, while in the second (long wave) they near 600 ppm. Although projected CO₂ concentrations in a methane economy are low in relation to other scenarios, the projections confirm that global climate warming is likely to be a major planetary concern throughout the twenty-first century. A second finding is that data on past growth of world per capita energy consumption group neatly into two pulses consistent with longwave theories in economics.     

Orbital forcing of Arctic climate: mechanisms of climate response and implications for continental glaciation

Progress in understanding how terrestrial ice volume is linked to Earths orbital configuration has been impeded by the cost of simulating climate system processes relevant to glaciation over orbital time scales (10³–10⁵ years). A compromise is usually made to represent the climate system by models that are averaged over one or more spatial dimensions or by three-dimensional models that are limited to simulating particular snapshots in time. We take advantage of the short equilibration time (10 years) of a climate model consisting of a three-dimensional atmosphere coupled to a simple slab ocean to derive the equilibrium climate response to accelerated variations in Earths orbital configuration over the past 165,000 years. Prominent decreases in ice melt and increases in snowfall are simulated during three time intervals near 26, 73, and 117 thousand years ago (ka) when aphelion was in late spring and obliquity was low. There were also significant decreases in ice melt and increases in snowfall near 97 and 142 ka when eccentricity was relatively large, aphelion was in late spring, and obliquity was high or near its long term mean. These glaciation-friendly time intervals correspond to prominent and secondary phases of terrestrial ice growth seen within the marine ¹⁸O record. Both dynamical and thermal effects contribute to the increases in snowfall during these periods, through increases in storm activity and the fraction of precipitation falling as snow. The majority of the mid- to high latitude response to orbital forcing is organized by the properties of sea ice, through its influence on radiative feedbacks that nearly double the size of the orbital forcing as well as its influence on the seasonal evolution of the latitudinal temperature gradient.     

An observational study of heat fluxes and their relationships with net radiation

Experimental evidence indicates that the diurnal behaviour of the fluxes of heat into the ground and into the atmosphere versus the net flux of radiation can be modelled by closed curves, the hourly values folowing one another in either a clockwise or counter clockfashion. A general formulation to express the different heat fluxes as a function of net radiation is proposed. This relationship between the different heat fluxes and can be expressed as a sum of three terms: the first indicates a direct proportionality to , the second gives the deviation from linearity and depends on /t, and the third gives the value of the flux when = 0. The formulae are then expressed as a function of time and the ratios between the heat fluxes and are evaluated. A comparison with the approximations generally used shows that the latter may be considered as particular cases of the more general equations proposed here.     

Ozone production due to emissions from vegetation burning

Ozone has been observed in elevated concentrations by satellites over areas previously believed to be background. There is meteorological evidence, that these ozone plumes found over the Atlantic Ocean originate from vegetation fires on the African continent.In a previous study (DECAFE-88), we have investigated ozone and assumed precursor compounds over African tropical forest regions. Our measurements revealed large photosmog layers at altitudes from 1.5 to 4 km. Both chemical and meteorological evidence point to savanna fires up to several thousand km upwind as sources.Here we describe ozone mixing ratios observed over western Africa and compare ozone production ratios from different field measurement campaigns related to vegetation burning. We find that air masses containing photosmog ingredients require several days to develop their oxidation potential, similar to what is known from air polluted by emissions from fossil fuel burning. Finally, we estimate the global ozone production due to vegetation fires and conclude that this source is comparable in strength to the stratospheric input.     

Public expectation as an element of human perception of climate change

Human expectations regarding weather and climate sometimes lead to perceptions of climate change which are not supported by observational evidence. This paper analyses two very characteristic complaints about current climate in Switzerland, i.e., the lack of snow in winter and the lack of sunshine in summer, through a statistical investigation of climatological data. As one major problem of public perception of climate in mid-latitude regions is linked to the strong variability of the climatic parameters, the paper suggests means of presenting climatic data which include a measure of this variability. Such presentations would help overcome the common confusion between the terms weather and climate, and stress the fact that short-term extreme events are not necessarily indicative of a long-term shift in climate.     

Occurrence of roll circulations in a shallow boundary layer

Meteorological measurements taken at the Näsudden wind turbine site during slightly unstable conditions have been analyzed. The height of the convective boundary layer (CBL) was rather low, varying between 60 and 300 m. Turbulence statistics near the ground followed Monin-Obukhov similarity, whereas the remaining part of the boundary layer can be regarded as a near neutral upper layer. In 55% of the runs, horizontal roll vortices were found. Those were the most unstable runs, with -z _i/L > 5. Spectra and co-spectra are used to identify the structures. Three roll indicators were identified: (i) a low frequency peak in the spectrum of the lateral component at low level; (ii) a corresponding increase in the vertical component at mid-CBL; (iii) a positive covariance {ovvw} together with positive wind shear in the lateral direction (V/z) in the CBL. By applying these indicators, it is possible to show that horizontal roll circulations are likely to be a common phenomenon over the Baltic during late summer and early winter.