A climate change simulation starting from 1935

Due to restrictions in the available computing resources and a lack of suitable observational data, transient climate change experiments with global coupled ocean-atmosphere models have been started from an initial state at equilibrium with the present day forcing. The historical development of greenhouse gas forcing from the onset of industrialization until the present has therefore been neglected. Studies with simplified models have shown that this cold start error leads to a serious underestimation of the anthropogenic global warming. In the present study, a 150-year integration has been carried out with a global coupled ocean-atmosphere model starting from the greenhouse gas concentration observed in 1935, i.e., at an early time of industrialization. The model was forced with observed greenhouse gas concentrations up to 1985, and with the equivalent C0₂ concentrations stipulated in Scenario A (Business as Usual) of the Intergovernmental Panel on Climate Change from 1985 to 2085. The early starting date alleviates some of the cold start problems. The global mean near surface temperature change in 2085 is about 0.3 K (ca. 10%) higher in the early industrialization experiment than in an integration with the same model and identical Scenario A greenhouse gas forcing, but with a start date in 1985. Comparisons between the experiments with early and late start dates show considerable differences in the amplitude of the regional climate change patterns, particularly for sea level. The early industrialization experiment can be used to obtain a first estimate of the detection time for a greenhouse-gas-induced near-surface temperature signal. Detection time estimates are obtained using globally and zonally averaged data from the experiment and a long control run, as well as principal component time series describing the evolution of the dominant signal and noise modes. The latter approach yields the earliest detection time (in the decade 1990–2000) for the time-evolving near-surface temperature signal. For global-mean temperatures or for temperatures averaged between 45°N and 45°S, the signal detection times are in the decades 2015–2025 and 2005–2015, respectively. The reduction of the cold start error in the early industrialization experiment makes it possible to separate the near-surface temperature signal from the noise about one decade earlier than in the experiment starting in 1985. We stress that these detection times are only valid in the context of the coupled model's internally-generated natural variability, which possibly underestimates low frequency fluctuations and does not incorporate the variance associated with changes in external forcing factors, such as anthropogenic sulfate aerosols, solar variability or volcanic dust.     

Northern hemisphere atmospheric response to changes of atlantic ocean SST on decadal time scales: a GCM experiment

Analyses indicate that the Atlantic Ocean seasurface temperature (SST) was considerably colder at the beginning than in the middle of the century. In parallel, a systematic change in the North Atlantic sea-level pressure (SLP) pattern was observed. To find out whether the SST and SLP changes analyzed are consistent, which would indicate that the SST change was real and not an instrumental artifact, a response experiment with a low-resolution (T21) atmospheric GCM was performed. Two perpetual January simulations were conducted, which differ solely in the Atlantic Ocean (40° S-60° N) SST: the cold simulation utilizes the SSTs for the period 1904–1913; the warm simulation uses the SSTs for the period 1951–1960. Also, a control run with the model's standard SST somewhat between the cold and warm SST was made. For the response analysis, a rigorous statistical approach was taken. First, the null hypothesis of identical horizontal distributions was subjected to a multivariate significance test. Second, the level of recurrence was estimated. The multivariate statistical approaches are based on hierarchies of test models. We examined three different hierarchies: a scale-dependent hierarchy based on spherical harmonics (S), and two physically motivated ones, one based on the barotropic normal modes of the mean 300 hPa flow (B) and one based on the eigenmodes of the advection diffusion operator at 1000 hPa (A). The intercomparison of the cold and warm experiments indicates a signal in the geostrophic stream function that in the S-hierarchy is significantly nonzero and highly recurrent. In the A-hierarchy, the low level temperature field is identified as being significantly and recurrently affected by the altered SST distribution. The SLP signal is reasonably similar to the SLP change observed. Unexpectedly, the upper level stream-function signal does not appear to be significantly nonzero in the B-hierarchy. If, however, the pairs of experiments warm versus control and cold versus control are examined in the B-hierarchy, a highly significant and recurrent signal emerges. We conclude that the cold versus warm response is not a small disturbance that would allow the signal to be described by eigenmodes of the linear system. An analysis of the three-dimensional structure of the signal leads to the hypothesis that two different mechanisms are acting to modify the model's mean state. At low levels, local heating and advection are dominant, but at upper levels the extratropical signal is a remote responce to modifications of the tropical convection.This paper was presented at the International Conference on Modelling of Global Climate Change and Variability, held in Hamburg 11–15 September 1989 under the auspices of the Meteorological Institute of the University of Hamburg and the Max Planck Institute for Meteorology. Guest Editor for these papers is Dr. L. Dilmenil.AWI Publication no. 254     

A mechanism of Alpine lee cyclogenesis as revealed by a quasigeostrophic variational filter

Summary Mechanisms associated with Alpine lee cyclogenesis during the early phase of their generation are investigated using a variational quasigeostrophic filter technique. It was possible to extract the quasigeostrophic signal from the available analyzed real data set.The results presented here are for the 11–12 March 1982, an example of so-called orographically induced lee cyclogenesis. Non-quasigeostrophic fields, calculated as a difference between observations and the quasigeostrophic fields, show significant magnitudes indicating the possible importance of non-quasigeostrophic processes. A dipole structure in the residual geopotential field was observed, similar to the results of numerical model experiments. Also, a strong upper-level non-quasigeostrophic divergence was found in the Alpine region 24 hours prior to lee cyclogenesis, lasting for 6–12 hours. On the other hand, quasigeostrophic results indicate only a local effect of mountain slopes, suggesting possibly a dominant role of the low-level blocking. A hypothetical scenario of Alpine lee cyclogenesis is proposed, based on results obtained here.With 14 Figures     

Arctic hazes in summer over Greenland and the North American Arctic: II. Nature and concentrations of accumulation-mode and giant particles

Airborne measurements made during August 1985 over Greenland and its environs show that both accumulation-mode (0.1 m D2.0 m) and giant (D2 m) particles were present in relatively high concentrations in arctic haze layers and that the accumulation-mode particles dominated light scattering. Particles with diameters (D) between 1 and 4 m consisted predominately of mixed materials, small and dense inclusions, and probably organic compounds containing sulfur. Many of the particles from 0.1 to 1 m in diameter were also of mixed composition, with sulfuric acid, ammonium sulfate and organics probably the dominant constituents.     

Time-dependent greenhouse warming computations with a coupled ocean-atmosphere model

Climate changes during the next 100 years caused by anthropogenic emissions of greenhouse gases have been simulated for the Intergovernmental Panel on Climate Change Scenarios A (business as usual) and D (accelerated policies) using a coupled ocean-atmosphere general circulation model. In the global average, the near-surface temperature rises by 2.6 K in Scenario A and by 0.6 K in Scenario D. The global patterns of climate change for both IPCC scenarios and for a third step-function 2 x CO₂ experiment were found to be very similar. The warming delay over the oceans is larger than found in simulations with atmospheric general circulation models coupled to mixed-layer models, leading to a more pronounced land-sea contrast and a weaker warming (and in some regions even an initial cooling) in the Southern Ocean. During the first forty years, the global warming and sea level rise due to the thermal expansion of the ocean are significantly slower than estimated previously from box-diffusion-upwelling models, but the major part of this delay can be attributed to the previous warming history prior to the start of present coupled ocean-atmosphere model integration (cold start).     

Trace gas exchange at the air/water interface: Measurements of mass accommodation coefficients

A liquid jet of 90 m diameter and variable length has been utilized to determine absorption rates and, hence, mass accommodation coefficients , of atmospheric trace gases. The compounds investigated are HCl (0.01), HNO₃ (0.01), N₂O₅ (0.005), peroxyacetyl nitrate (>0.001), and HONO (0.005). It is concluded that the absorption of these trace gases by liquid atmospheric water is not significantly retarded by interfacial mass transport. The strengths and limitations of the liquid jet technique for measuring mass accommodation coefficients are explored.     

A study of the sensitivity of bare soil evaporation schemes to soil surface wetness,using the coupled soil moisture and surface temperature prediction model,BARESOIL

Summary The performance of evaporation schemes with and approach and their combination within resistance representation of evaporation from bare soil surface is discussed. For this purpose nine schemes, based on different functions of or , on the ratio of the volumetric soil moisture content and its saturated value are used.The quality of the chosen schemes has been evaluated using the results of time integration by the coupled soil moisture and surface temperature prediction model, BARESOIL, using in situ data. A sensitivity analysis was made using two sets of data derived from the volumetric soil moisture content of the top soil layer. One with values below the wilting point (0.17 m³m^–3) and the second with values above 0.20m³m^–3. Data sets were obtained at the experimental site Rimski anevi, Yugoslavia, from the bare surface of a chernozem soil.With 4 Figures     

A study of a fair weather boundary layer in TOGA-COARE: Parameterization of surface fluxes in large scale and regional models for light wind conditions

The study focuses on a way to parameterize the effect of subgrid scale convective motions on surface fluxes in large scale and regional models for the case of light surface winds. As previously proposed, these subgrid effects are assumed to scale with the convection intensity through the relationship: where is the mean velocity of the wind, U₀ the velocity of the mean wind, w_* the free convection velocity, and an empirical coefficient to be determined. Both observations and numerical simulation are presently used to determine the free convection coefficient .Large eddy simulation of a fair weather convective boundary layer case observed during TOGA-COARE is performed. Comparisons between observations and the simulation of surface properties and vertical profiles in the planetary boundary layer are presented. The simulated vertical turbulent fluxes of heat, moisture and buoyancy range well within estimates from aircraft measurements.The most important result is that the true free convection coefficient , directly estimated from simulation, leads to a value of 0.65, smaller than the ones estimated from temporal and spatial variances. Using observations and simulation, estimates of from temporal and spatial variances are obtained with similar values 0.8. From both theoretical derivations and numerical computations, it is shown that estimates of the true from variances are possible but only after applying a correction factor equal to 0.8. If this correction is not used, is overestimated by about 25%. The time and space sampling problem is also addressed in using numerical simulations.     

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.     

Equatorial disturbances,monsoons, and 1982–1983 ENSO

Summary Interannual modes are described in terms of three-month running mean anomaly winds (u,v), outgoing longwave radiation (OLR), and sea surface temperature (T _* ). Normal atmospheric monsoon circulations are defined by long-term average winds (u _n,v _n) computed every month from January to December. Daily winds are grouped into three frequency bands, i.e., 30–60 day filtered winds (u _L,v _L); 7–20 day filtered winds (u _M,v _M); and 2–6 day filtered winds (u _S,v _S). Three-month running mean anomaly kinetic energy (signified asK _L , K _M , andK _S , respectively) is then introduced as a measure of interannual variation of equatorial disturbance activity. Interestingly, all of theseK _L , K _M , andK _S perturbations propagate slowly eastward with same phase speed (0.3 ms^–1) as ENSO modes. Associated with this eastward propagation is a positive (negative) correlation between interannual disturbance activity (K _L , K _M , K _S ) and interannualu (OLR) modes. Namely, (K _L , K _M , K _S ) becomes more pronounced than usual nearly simultaneously with the arrival of westerlyu and negativeOLR (above normal convection) perturbutions. In these disturbed areas with (K _L , K _M , K _S >0), upper ocean mixing tends to increase, resulting in decreased sea surface temperature, i.e.T _* 0. Thus, groups (not individual) of equatorial disturbances appear to play an important role in determiningT _* variations on interannual time scales. HighestT _* occurs about 3 months prior to the lowestOLR (convection) due primarily to radiational effects. This favors the eastward propagation of ENSO modes. The interannualT _* variations are also controlled by the prevailing monsoonal zonal windsu _n, as well as the zonal advection of sea surface temperature on interannual time scales. Over the central Pacific, all of the above mentioned physical processes contribute to the intensification of eastward propagating ENSO modes. Over the Indian Ocean, on the other hand, some of the physical processes become insignificant, or even compensated for by other processes. This results in less pronounced ENSO modes over the Indian Ocean.With 10 FiguresContribution No. 89-6, Department of Meteorology, University of Hawaii, Honolulu, Hawaii.     

Bioaccumulation of Cobalt in Parmelia sulcata

Six locations across mainland Portugal were selected for exposing Parmelia sulcata, for a one-year period (8 months for one site), with simultaneous measurement of total (dry + wet) deposition (one-month periods). The exposed lichens and the total (dry + wet) deposition were analysed for cobalt contents by INAA (instrumental neutron activation analysis) and ICP-MS (inductively coupled plasma mass spectroscopy), respectively. The designated wet deposition was evaluated through the collected water volume; the designated dry deposition was assessed after the (dried) residual mass of the wet deposition. An excellent agreement between Co contents in exposed lichens and the cumulative (1) Co contents in the dry deposition, (2) dry deposition, and (3) wet deposition has been found for the locations with alternate drought and precipitation months, high dry deposition, and high Co contents in the latter. Continuous rainfall was found to hinder the Co accumulation in the lichen due to its release from the lichen and/or lower Co contents in the dry deposition. At three locations, P. sulcata Co contents, after subtraction of the background (before exposure), equalled or exceeded the Co contents in the cumulative dry deposition at the end of the exposure time. The optimal exposure period for this species likely depends on the exposure conditions.     

Hybrid modeling in meteorological applications Part I: Concepts and approaches

Summary Hybrid modeling entails the combination of a numerical weather prediction model and a symbolic model. The latter uses symbolic objects, their characterizing attributes and sets of behavioral constrains which prescribe changes in the states of these objects as functions of time, space, and other imposed quantitative or heuristic conditions. Integration of these two modeling components for an on-line, real-time, operational system is feasible only if both the numerical and the symbolic model can be executed in a distributed mode, i.e. at a user's location rather than in a central weather service office. This condition entails the design of a numerical model that can run on relatively inexpensive desktop workstations or high-end PCs. Given such a capability, the output from the numerical model can be used to satisfy a number of behavioral constraints of objects (such as thunderstorm, blizzard, etc.), defined in the symbolic model. These constraints can be embedded invisibly as functions of time and pixel location on the computer screen, to be called upon as soon as the respective object is activated, e.g. by placing an icon on the screen.To make such a hybrid weather prediction model responsive to details in topography, it will have to be able to interface with a geographic information system (GIS) database. Since such databases can be very voluminous, management procedures for indexing and rapid information retrieval have to be instituted. The approach discussed here involves restructuring of given GIS data into B+ trees.The hybrid prediction model whose design is described in this paper, executes very quickly on a PC (e.g. a 33 megahertz Intel 80486 chip based machine). It allows assimilation of locally generated observational data to improve forecast quality, and can respond to queries of a highly specialized nature in support of tactical decisions within the time frame between nowcasting (3 hours) and 24 hours.With 6 Figures     

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.     

ALPEX-Simulation

Summary In a project ALPEX-Simulation, sponsored by the Österreichischer Fond zur Förderung der wissenschaftlichen Forschung (FWF), all eight cases of ALPEX-SOP cyclones were numerically simulated with a fine mesh isentropic model of the atmosphere. These numerical simulations in six-hourly intervals allow a deeper insight into the synoptics and dynamics of the cyclogeneses in the Western Mediterranean, especially into the genesis of the two basic types of cyclones: the so-called Überströmungs-type and Vorderseiten-type. In the first phase of cyclogenesis of the Überströmungs-type, the blocking and flow splitting of the cold air due to the Alps and the canalization between the Alps and the Massif Central are important. Cold air flows cyclonically around the western part of the Alps, creating a vorticity maximum at the south western edge of the Alpine, bow and leads also to an enhanced PV. In connection with warm air in the Mediterranean, a strong baroclinic zone is generated. The interaction between the arriving PV maximum in the upper troposphere and the enhanced PV at the bottom leads to cyclogenesis in the Western Mediterranean. In the case of the Vorderseiten-type warm air advection dominates with the exception of a shallow layer of cold air in the inner Po-Valley, which is shielded by the Alpine ridge. A well-pronounced PV maximum builds up and couples with the PV maximum arriving at upper levels, even before the cold air, coming from the north-west, has surrounded the Alps. The cold air only intensifies the development by raising the baroclinity. Therefore, the Vorderseiten-cyclogenesis is an orographically modified cyclogenesis, in the course of which the cyclonic development is triggered by the Alps, whereas the Überströmungs-cyclogenesis is an orographically induced cyclogenesis i.e. a true lee cyclogenesis.With 14 FiguresDied in a tragic traffic accident on June 6, 1993.     

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.     

Effects of water vapour on the structure parameter of the refractive index for near-infrared radiation

The refractive index of air (n) mainly depends on temperature andwater vapour content. For near-infrared radiation, temperature isthe main determining factor. To determinethe structure parameter of temperature ( ) fromthe structure parameter of the refractive index( ), the influence ofwater vapour content on n needs to be taken into account as acorrection. Three levels of approximation are presented. The firstinvolves the standard deviations of T and q (_T and_q) as well as the correlation coefficient betweenT and q (R_Tq). The second approximation involvesR_Tq and the Bowen ratio (), and the last uses only the Bowen ratio.The latter is the classical Bowen ratio correction.Evaluation of the validity of the assumptions used in the derivationreveals that a large error may be introduced (for situations withR_Tq < 0.8) by replacing _T/_qwith the Bowen ratio.Consequently, estimating from using the Bowen ratiocan contain errors of 5–40% for || < 1. Adding R_Tqas extra information reduces the errors only marginally. For || > 1,the correction is small, and all three approximations give errors of lessthan 1% in . When is used to compute the sensible heat flux,the influence of the quality of the correction forwater vapour fluctuations on the measured energy balance is small: forsmall ||, the correction is large, but the absolute value ofthe heat flux is small, whereas for large ||, the correctionis insignificant.     

Propagation of the sea-surface temperature anomalies in the Tropical Atlantic

Summary The evolving modes of the sea-surface temperature (SST) in the Tropical Atlantic on the short interannual (IA) timescale were obtained by performing the extended empirical orthogonal function (EEOF) analyses on this variable separately for the 106-year (1871–1976) and 20-year (1881–1900; 1901–1920; 1921–1940; 1941–1960) periods. The equatorial and inter-hemispheric patterns manifest in the first EEOF mode of each analysis as part of the short IA evolution of the SST anomalies in the Tropical Atlantic. Another outstanding feature of the first EEOF mode of each analysis concerns the propagations of the SST anomalies in the meridional direction within the 20°N–20°S band and in the zonal direction in the sector 40°W–20°W. For all analyses, the SST anomalies propagate northward from the equator to 15°N and southward from 20°N to 15°N, with the same sign anomalies merging approximately at 15°N. On the other hand, the SST anomalies propagate westward in the sector 40°W–20°W with a propagation rate close to that of the phase speed of the fastest baroclinic Rossby wave in the ocean. So, the observed propagations of the SST anomalies in the 20°N–20°S band might result from the combined effect of the surface oceanic currents in this band and the baroclinic Rossby waves in the ocean.     

On Trajectory Curvature as a Selection Criterion for Valid Lagrangian Stochastic Dispersion Models

Recently Wilson and Flesch (Boundary-Layer Meteorology, 84, 411-426, 1997) suggested that the average increment d z to the orientation = arctan(w/u) of the Lagrangian velocity-fluctuation vector can be used to distinguish the better Lagrangian stochastic models within the well-mixed class. Here it is demonstrated that the specification of d z constitutes neither a sufficient or universally applicable criterion to distinguish the better Lagrangian stochastic models within the well-mixed class. The hypothesis made by Wilson and Flesch that Lagrangian stochastic models with /P_E irrotational are zero-spin models, having d z=0, is proven     

On some new statistical characteristics of the wind “Koshava”

Summary Koshava is a gusty wind of moderate to strong intensity, blowing from the south-eastern direction, over the area of the Republic of Serbia. It is caused by the interaction between the synoptic circulation and the orography of the Carpathian and the Balkan mountains. The Koshava wind can damage buildings, factories and industrial plants or city infrastructure. Therefore it is important to estimate its gust and the gustiness factor on the basis of the measured data.This paper discusses a statistical analysis of wind data in the maximum influence area of the Koshava wind in the periods of maximum duration of Koshava. The focus of the paper is the examination of urban and suburban effects on Koshava wind and the correlation between the instantaneous maximum wind speed and the hourly mean wind speed. The best fitting with various empirical distributions is proposed.With 10 Figures