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 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.     

Outcome and value uncertainties in global-change policy

Choices among environmental policies can be informed by analysis of the potential physical, biological, and social outcomes of alternative choices, and analysis of social preferences among these outcomes. Frequently, however, the consequences of alternative policies cannot be accurately predicted because of substantial outcome uncertainties concerning physical, chemical, biological, and social processes linking policy choices to consequences. Similarly, assessments of social preferences among alternative outcomes are limited by value uncertainties arising from limitations of moral principles, the absence of economic markets for many environmental attributes, and other factors. Outcome and value uncertainties relevant to global-change policy are described and their magnitudes are examined for two cases: stratospheric-ozone depletion and global climate change. Analysis of information available in the mid 1980s, when international ozone regulations were adopted, suggests that contemporary uncertainties surrounding CFC emissions and the atmospheric response were so large that plausible ozone depletion, absent regulation, ranged from negligible to catastrophic, a range that exceeded the plausible effect of the regulations considered. Analysis of climate change suggests that, important as outcome uncertainties are, uncertainties about values may be even more important for policy choice.     

Transport-Dissipation Analytical Solutions to the E-∈Turbulence Model and their Role in Predictions of the Neutral ABL

E- turbulence model predictions of the neutralatmospheric boundary layer (NABL) are reinvestigated to determine thecause for turbulence overpredictions found in previous applications. Analytical solutions to the coupled E and equations for the case of steady balance between transport and dissipation terms, the dominant balance just below the NABL top, are derived. It is found that analytical turbulence profiles laminarizeat a finite height only for values of closure parameter ratio c ₂ /_e equal toor slightly greater than one, with laminarization as z for greater . The point = 2 is additionally foundthat where analytical turbulent length scale (l) profilesmade a transition from ones ofdecreasing ( < 2) to increasing ( > 2)values with height. Numerically predicted profiles near the NABL topare consistent with analytical findings. The height-increasingvalues of l predicted throughout the NABL with standard values ofclosure parameters thus appear a consequence of 2.5(> 2), implied by these values (c ₂ = 1.92, _{= 1.3}, _{e = 1}). Comparison of numericalpredictions with DNS data shows that turbulence overpredictions obtained with standard-valued parameters are rectifiedby resetting and _e to 1.1 and 1.6, respectively, giving, with c ₂ = 1.92, 1.3, and laminarization of the NABL's cappingtransport-dissipation region at a finite height.     

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.     

Tropical Glacier and Ice Core Evidence of Climate Change on Annual to Millennial Time Scales

This paper examines the potential of the stable isotopic ratios, ¹⁸O/¹⁶O ( ¹⁸O_ice)and ²H/¹H ( D_ice), preserved in mid to low latitude glaciers as a toolfor paleoclimate reconstruction. Ice cores are particularly valuable as they contain additional data, such as dust concentrations, aerosol chemistry, and accumulation rates, that can be combined with the isotopic information to assist with inferences about the regional climate conditions prevailing at the time of deposition. We use a collection of multi-proxy ice core histories to explore the ¹⁸O-climate relationship over the last 25,000 years that includes both Late Glacial Stage (LGS) and Holocene climate conditions. These results suggest that on centennial to millennial time scales atmospheric temperature is the principal control on the ¹⁸O_ice of the snowfall that sustains these high mountainice fields.Decadally averaged ¹⁸O_ice records from threeAndean and three Tibetan ice cores are composited to produce a low latitude ¹⁸O_ice history for the last millennium. Comparison ofthis ice core composite with the Northern Hemisphere proxy record (1000–2000A.D.) reconstructed by Mann et al. (1999) and measured temperatures(1856–2000) reported by Jones et al. (1999) suggests the ice cores have captured the decadal scale variability in the global temperature trends. These ice cores show a 20th century isotopic enrichment that suggests a large scale warming is underway at low latitudes. The rate of this isotopically inferred warming is amplified at higher elevations over the Tibetan Plateau while amplification in the Andes is latitude dependent with enrichment (warming) increasing equatorward. In concert with this apparent warming, in situobservations reveal that tropical glaciers are currently disappearing. A brief overview of the loss of these tropical data archives over the last 30 years is presented along with evaluation of recent changes in mean ¹⁸O_ice composition. The isotopic composition of precipitation should be viewed not only as a powerful proxy indicator of climate change, but also as an additional parameter to aid our understanding of the linkages between changes in the hydrologic cycle and global climate.     

Wind and temperature profile prediction for diabatic surface layers including strong inversion cases

The micrometeorological research program in Antarctica has provided extensive data on wind and temperature profile structure under strong to extreme inversion conditions (Dalrymple et al., 1966; Lettau et al., 1977). The basic similarity hypotheses and limiting conditions for prediction of diabatic surface layer profiles are summarized. The model by Businger et al. (1971) for dimensionless shear and temperature gradients is revised to conform with the new results for strong stability. A novel similarity hypothesis is introduced to complete the step from shear and gradient prediction to prediction of absolute wind speed, wind energy, and temperature on the basis of prescribed external factors of surface layer structure. The physics of interactions between predicted profile tilting and curving are discussed and used to explain several micrometeorological paradoxes, including that of the elevated minimum of air temperature observed occasionally near the active surface when the energy budget is of the nocturnal type.     

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.     

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     

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.     

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.     

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.     

Fluctuations and trends of annual precipitation in different climatic regions of Croatia

Summary Secular series of annual precipitation over Croatia have been studied at three stations representing the different climatic regions of Croatia: Osijek (continental precipitation climate), Zagreb-Gri (continental precipitation climate with maritime influence) and Crikvenica (maritime precipitation climate).The time series analysis has been deduced by using a quick test for stationarity according to Schönwiese and Malcher, moving average filters, the Mann-Kendall rank statistic, and a progressive test for trend according to Sneyers. From this analysis, the stationarity in time series of annual precipitation totals is stated over the entire interval at all locations except for a very short interval during the first decade of the twentieth century over the continental lowland (Osijek). A generally decreasing trend is present over the entire interval, but is statistically significant only in the continental lowland (Osijek).With 4 Figures     

Reactions of alkoxy radicals under atmospheric conditions: The relative importance of decomposition versus reaction with O2

The reactions of alkoxy radicals determine to a large extent the products formed during the atmospheric degradations of emitted organic compounds. Experimental data concerning the decompositions, 1,5-H shift isomerizations and reactions with O₂ of several classes of alkoxy radicals are inconsistent with literature estimations of their absolute or relative rate constants. An alternative, although empirical, method for assessing the relative importance under atmospheric conditions of the reactions of alkoxy radicals with O₂ versus decomposition was derived. This estimation method utilizes the differences in the heats of reaction, (H)=(H_{decomposition}–H_{O 2 reaction}), between these two reactions pathways. For (H)[22–0.5(H_{O 2 reaction})], alkoxy radical decomposition dominates over the reaction with O₂ at room temperature and atmospheric pressure of air, while for (H)[25-0.5(H_{O 2 reaction})], the O₂ reaction dominates over decomposition (where the units of H are in kcal mol^–1). The utility and shortcomings of this approach are discussed. It is concluded that further studies concerning the reactions of alkoxy radicals are needed.     

Australian public perception of the greenhouse issue

During 1987 and 1988 in Australia there have been two national meetings on the greenhouse effect and a campaign designed to increase public awareness. A study of the backgrounds, level of comprehension and attitudes of attendees at two state Greenhouse-88 meetings has been undertaken by means of a questionnaire survey and a set of personal interviews. Two crucial caveats pertain: (1) some of the questions reflect the prejudices of the author who is an atmospheric scientist and much of the interpretation is similarly tainted and (2) the respondents comprise a small, self-selected and probably highly motivated group. All the ensuing results should be viewed in the context of these caveats.Over 97% of the respondents believe that action should be taken now to alleviate or mitigate the effects of increased greenhouse gases. Despite the fact that the majority of the 321 respondents are professional people (73%) and that over 53% have tertiary level educational qualifications, there was a failure to grasp some fundamental issues. For example, only 120 (37%) correctly recognized that N₂ is not a greenhouse agent whilst also agreeing that CO₂, CH₄ and CFCs are greenhouse agents. On the other hand, the respondents generally demanded a relatively low level of confidence (50% to 70% certainty) about the greenhouse issue from scientists before action is taken. Sixty-four percent believe that life will be worse for them and/or their children in Australia in Greenhouse 2025 with the youngest age range being the second most pessimistic group about the future.Relatively little interest was shown in the possibility of obtaining more information on topics that interest climatic scientists such as the validity of measured temperature trends and inadequacies/errors in climate models but more information was desired on the social and economic implication and, interestingly, on the scientific background to the issues. Overall, teachers are perceived as trying to increase understanding; whereas politicians, multinational corporations, the media and some extreme environmentalists are perceived as often attempting to deceive intentionally. Scientists are seen as neither especially malevolent nor benign. A possible conclusion which might be drawn is that by attempting to simplify issues for public debate, scientists may significantly reduce, or even remove entirely, any chance of re-association of connected issues by members of the public. Perhaps more importantly, scientists need to recognize and learn to use the knowledge that the public may have the right views for the wrong reasons.     

Breaking global temperature records after Mt. Pinatubo

The probability of setting global temperature records is reconsidered in light of cooling due to the Mt. Pinatubo eruption. The cooling resets temperature by moving it away from the top of its 100-year range. Depending on the statistical model for temperature, this reset can lead to a much lower probability for a record in the next few years. The exercise illustrates how record setting depends on the underlying model, the current record value, and - if there is serial correlation- the current temperature.     

A Comparative Analysis of Transpiration and Bare Soil Evaporation

Transpiration E^v and bare soil evaporation E^b processes are comparatively analysed assuming homogeneous and inhomogeneous areal distributions of volumetric soil moisture content . For a homogeneous areal distribution of we use a deterministic model, while for inhomogeneous distributions a statistical-deterministic diagnostic surface energy balance model is applied. The areal variations of are simulated by Monte-Carlo runs assuming normal distributions of .The numerical experiments are performed for loam. In the experiments we used different parameterizations for vegetation and bare soil surface resistances and strong atmospheric forcing. According to the results theE^v()-^Eb() differences are great, especially in dry conditions. In spite of this, the available energy flux curves of vegetation A^v() and bare soil A^b() surfaces differ much less than the E^v() and ^Eb() curves. The results suggest that E^v is much more non-linearly related to environmental conditions than ^Eb. Both E^v and ^Eb depend on the distribution of , the wetness regime and the parameterization used. With the parameterizations, ^Eb showed greater variations than E^v. These results are valid when there are no advective effects or mesoscale circulation patterns and the stratification is unstable.     

Periodicity of annual precipitation in different climate regions of Croatia

Summary The periodicity of a 100-year series of annual precipitation over Croatia has been studied by means of power spectrum analysis at 3 stations representing the different climatic regions of Croatia. The annual precipitation variance spectra in the continental lowland (Osijek) and at the north East Adriatic coast (Crikvenica) can be fitted by Markov white noise continuum, but in the transitional region between the Dinaric Alps and the Pannonian lowland (Zagreb-Gri) a non-white noise continuum is necessary. Quasi-periodic oscillations appear in two spectra ranges: short (2.2 and 4.7 years) and medium (25.0 and 33.3 years). These results are compared with those of other authors for other parts of the Europe.With 2 Figures     

Application of a “Big-Tree” model to regional climate modeling: a sensitivity study

Summary ¶In order to better understand land-atmosphere interactions and increase the predictability of climate models, it is important to investigate the role of forest representation in climate modeling. Corresponding to the big-leaf model commonly employed in land surface schemes to represent the effects of a forest, a so called big-tree model, which uses multi-layer vegetation to represent the vertical canopy heterogeneity, was introduced and incorporated into the National Center for Atmospheric Research (NCAR) regional climate model RegCM2, to make the vegetation model more physically based. Using this augmented RegCM2 and station data for China during 1991 Meiyu season, we performed 10 experiments to investigate the effects of the application of the big-tree model on the summer monsoon climate.With the big-tree model incorporated into the regional climate model, some climate characteristics, e.g. the 3-month-mean surface temperature, circulation, and precipitation, are significantly and systematically changed over the model domain, and the change of the characteristics differs depending on the area. Due to the better representation of the shading effect in the big-tree model, the temperature of the lower layer atmosphere above the plant canopy is increased, which further influences the 850hPa temperature. In addition, there are significant decreases in the mean latent heat fluxes (within 20–30W/m²) in the three areas of the model domain.The application of the big-tree model influences not only the simulated climate of the forested area, but also that of the whole model domain, and its impact is greater on the lower atmosphere than on the upper atmosphere. The simulated rainfall and surface temperature deviate from the originally simulated result and are (or seem to be) closer to the observations, which implies that an appropriate representation of the big-tree model may improve the simulation of the summer monsoon climate.We also find that the simulated climate is sensitive to some big-tree parameter values and schemes, such as the shape, height, zero-plane displacement height and mixing-length scheme. The simulated local/grid differences may be very large although the simulated areal-average differences may be much lower. The area-average differences in the monthly-mean surface temperature and heat fluxes can amount to 0.5°C and 4W/m², respectively, which correspond to maximum local/grid differences of 3.0°C and 40W/m² respectively. It seems that the simulated climate is most sensitive to the parameter of the zero-plane displacement among the parameters studied.