An empirical formula to compute snow cover fraction in GCMs

There exists great uncertainty in parameterizing snow cover fraction in most general circulation models (GCMs) using various empirical formulae, which has great influence on the performance of GCMs. This work reviews the commonly used relationships between region-averaged snow depth (or snow water equivalent) and snow cover extent (or fraction) and suggests a new empirical formula to compute snow cover fraction, which only depends on the domain-averaged snow depth, for GCMs with different horizontal resolution. The new empirical formula is deduced based on the 10-yr (1978-1987) 0.5°× 0.5° weekly snow depth data of the scanning multichannel microwave radiometer (SMMR) driven from the Nimbus-7 Satellite. Its validation to estimate snow cover for various GCM resolutions was tested using the climatology of NOAA satellite-observed snow cover.     

Sensitivity of glaciation to initial snow cover,CO2, snow albedo,and oceanic roughness in the NCAR CCM

We present results from numerical experiments made with a GCM, the NCAR CCM1, that were designed to estimate the annual balance between snow-fall accumulation and ablation for geographically important land regions for a variety of conditions. We also attempt to assess the reliability of these results by investigating model sensitivity to changes in prescribed physical parameters. Experiments were run with an initial imposition of 1 m of (midwinter) snowcover over all northern hemisphere land points. Over Alaska, western Canada, Siberia, and the Tibetan Plateau the model tended to retain this snow cover through the summer and in some cases increase its depth as well. We define these regions as glaciation sensitive and note some correspondence between them and source regions for the Pleistocene ice sheets. An experiment with greatly reduced CO₂ (100 ppm) showed a tendency towards spontaneous glaciation, i.e., the model remained snow-covered throughout the summer over the same geographic regions noted above. With 200 ppm CO₂ (roughly equal to values at the last glacial maximum), snow cover over these regions did not quite survive the summer on a consistent basis. Combining 200 ppm CO₂ and 1 m of initial northern hemisphere snow cover yielded glaciation-sensitive conditions, agreeing remarkably well with locations undergoing glaciation during the Pleistocene. To assess the reliability of these results, we have determined minimal model uncertainty by varying two of the empirical coefficients in the model within physically plausible ranges. In one case surface roughness of all ocean gridpoints was reduced by an order of magnitude, leading to local 10% reductions in precipitation (snowfall), a change hard to distinguish from inherent model variability. In the other case, the fraction of a land grid square assumed to be occupied by snow cover for albedo purposes was varied from one-half to unity. Large changes occurred in the degree of summer melting, and in some cases the sign of the net balance changed as fractional snow cover was changed. We conclude that the model may be able to reveal regions sensitive to glaciation, but that it cannot yield a reliable quantitative computation of the magnitude of the net snow accumulation that can be implicitly or explicitly integrated through time.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     

An annual cycle of vegetation in a GCM. Part II: global impacts on climate and hydrology

The Met Office Hadley Centre Unified Model (HadAM3) with the tiled version of the Met Office Surface Exchange Scheme (MOSES2) land surface scheme is used to assess the impact of a comprehensive imposed vegetation annual cycle on global climate and hydrology. Two 25-year numerical experiments are completed: the first with structural vegetation characteristics (Leaf Area Index, LAI, canopy height, canopy water capacity, canopy heat capacity, albedo) held at annual mean values, the second with realistic seasonally varying vegetation characteristics. It is found that the seasonalities of latent heat flux and surface temperature are widely affected. The difference in latent heat flux between experiments is proportional to the difference in LAI. Summer growing season surface temperatures are between 1 and 4 K lower in the phenology experiment over a majority of grid points with a significant vegetation annual cycle. During winter, midlatitude surface temperatures are also cooler due to brighter surface albedo over low LAI surfaces whereas during the dry season in the tropics, characterized by dormant vegetation, surface temperatures are slightly warmer due to reduced transpiration. Precipitation is not as systematically affected as surface temperature by a vegetation annual cycle, but enhanced growing season precipitation rates are seen in regions where the latent heat flux (evaporation) difference is large. Differences between experiments in evapotranspiration, soil moisture storage, the timing of soil thaw, and canopy interception generate regional perturbations to surface and sub-surface runoff annual cycles in the model.     

SIMULTANEOUS DETERMINATION OF AEROSOL SIZE DISTRIBUTION AND REFRACTIVE INDEX AND SURFACE ALBEDO FROM RADIANCE－PART I: THEORY

Analyzing in detail the dependence of sky radiance on aerosol optical property and surface albedo, we present a new method for simultaneous determination of aerosol size distribution, its wavelength-dependent refractive index and surface albedo. The aerosol scattering phase function near 10o, its weighted phase function near 40o introduced in this paper, which can be inferred from sky radiance data, and the radiance near 90o are respectively used in retrieving the real part of refractive index, its imaginary part and surface albedo. Results in numerical experiments are satisfactory under different given conditions.     

Simultaneous determination of aerosol size distribution and refractive index and surface albedo from radiance- Part III: Parametrization and application

Based on the radiative transfer calculation results, three approximate expressions of the sky radiance in almucantar and its increment caused by surface albedo are presented. They are simple, but accurate enough. The dependence of the fitted aerosol scattering phase functions on refractive index is also studied, and its reasonable form is given. For Junge size distribution, the approximate equations of the phase functions with some special scattering angles are obtained. These approximate equations significantly simplify the retrieval algorithm of simultaneous determination of aerosol size distribution and its refractive index and surface albedo. This method can be realized with a microcomputer, and has been used to process and analyse the experimental data measured in Hefei of Anhui Province.     

Random cascade driven rainfall disaggregation for urban hydrology: An evaluation of six models and a new generator

Six variations of multiplicative random cascade models for generating fine-resolution (i.e., 5-minute interval) rainfall time series were evaluated for rainfall in Wroclaw, Poland. Of these variations, one included a new beta-normal generator for a microcanonical cascade. This newly proposed model successfully reproduces the statistical behavior of local 5-minute rainfalls, in terms of intermittency as well as variability. In contrast, both the canonical cascade models with either constant or time-scaled parameters and a microcanonical cascade model with a beta generator substantially underestimate 5-minute maximum rainfall intensities. The canonical models also fail to properly reproduce the intermittency of the rainfall process across a range of timescales. New observations are also made concerning the histograms of the breakdown coefficients (BDC). The tendency of the BDC histograms to have values exactly equal to 0.5 is identified and explained by the quality of pluviograph records. Moreover, the hierarchical evolution of BDC histograms from beta-like for long time steps to beta-normal histograms for short time steps is observed for the first time. The potential advantage is shown of synthetic high resolution rainfall time series generated by the revised microcanonical model for use in hydrology, especially hydrodynamic modelling of urban drainage networks.     

The past climate and sea-ice record from Iceland,Part 1: Data to A.D. 1780

A new reconstruction of the climate and sea-ice record for Iceland from medieval times to A.D. 1780 is presented, based on all available documentary sources. The importance of careful historical analysis to separate reliable from unreliable material is stressed, and these reconstructions are the first to have been produced using only reliable data. The major previous works on the subject (those of Thoroddsen, Koch, and Bergþórsson), which all include unreliable material, are discussed. Prior to A.D. 1600 the data are not considered to be full enough to permit a quantitative interpretation. For the period A.D. 1601 to 1780 decadal temperature and sea-ice indices are given.Although there is very little evidence for the first few centuries of settlement in Iceland (from c. 870 to c. 1170) the data suggest a fairly mild climatic period. Cold periods occurred around 1200, and at the end of the thirteenth century. The fourteenth century was very variable with a cold period in the 1350s to c. 1380. Between 1430 and c. 1560 there are very few contemporary sources and it is difficult to draw any conclusions on the climate during this time. The latter part of the sixteenth century was undoubtedly cold. From 1601 there are sufficient data to permit a decade by decade analysis. This shows a mild period between 1640 and 1670, and severe decades in the 1630s, 1690s, 1740s, and 1750s. Year to year and decade to decade variability is appreciable. The correlation between temperature and sea ice is not perfect but is still quite strong (similar to today). Because data have been gathered from different regions of Iceland it has been possible to demonstrate the spatial variability of Iceland's climate during the period 1601 to 1780. For example, during 1660 to 1700 there was a cooling in the north and west but warming in the south. The 1690s, the coldest decade of the Little Ice Age in Europe, was extremely cold in the west of Iceland, but less severe elsewhere.

---

     

An Updated Coupled Model for Land-Atmosphere Interaction. Part I： Simulations of Physical Processes

A new two-way land-atmosphere interaction model （R42_AVIM） is fulfilled by coupling the spectral atmospheric model （SAMIL_R42L9） developed at the State Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics, Institute of Atmospheric Physics, Chinese Academy of Sciences （LASG/IAP/CAS） with the land surface model, Atmosphere-Vegetation-Interaction-Model （AVIM）. In this coupled model, physical and biological components of AVIM are both included. Climate base state and land surface physical fluxes simulated by R42_AVIM are analyzed and compared with the results of R42_SSIB [which is coupled by SAMIL_R42L9 and Simplified Simple Biosphere （SSIB） models]. The results show the performance of the new model is closer to the observations. It can basically guarantee that the land surface energy budget is balanced, and can simulate June-July-August （JJA） and December-January- February （DJF） land surface air temperature, sensible heat flux, latent heat flux, precipitation, sea level pressure and other variables reasonably well. Compared with R42_SSIB, there are obvious improvements in the JJA simulations of surface air temperature and surface fluxes. Thus, this land-atmosphere coupled model will offer a good experiment platform for land-atmosphere interaction research.     

An Updated Coupled Model for Land-Atmosphere Interaction. Part II： Simulations of Biological Processes

In Part I, the authors succeeded in coupling the spectral atmospheric model （SAMIL_R42L9） developed at the State Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics, Institute of Atmospheric Physics, Chinese Academy of Sciences （LASG/IAP/CAS） with the land surface model, Atmosphere-Vegetation-Interaction-Model （AVIM） and analyzed the climate basic state and land surface physical fluxes simulated by R42_AVIM. In this Part Ⅱ, we further evaluate the simulated results of the biological processes, including leaf area index （LAI）, biomass and net primary productivity （NPP） etc. Results indicate that R42_AVIM can simulate the global distribution of LAI and has good consistency with the monthly mean LAI provided by Max Planck Institute for Meteorology. The simulated biomass corresponds reasonably to the vegetation classifications. In addition, the simulated annual mean NPP has a consistent distribution with the data provided by IGBP and MODIS, and compares well with the work in literature. This land-atmosphere coupled model will offer a new experiment tool for the research on the two-way interaction between climate and biosphere, and the global terrestrial ecosystem carbon cycle.     

Trend analysis of climatic variables and their relation to snow cover and water availability in the Central Himalayas: a case study of Langtang Basin,Nepal

Theoretical and Applied Climatology - This paper investigates the long-term monthly, seasonal, and annual trends of climatic variables, snow cover extent, and discharge in Langtang Basin, Central...     

MODIS/Terra observed snow cover over the Tibet Plateau: distribution, variation and possible connection with the East Asian Summer Monsoon (EASM)

Based on the snow cover fraction (SCF) data acquired from the Moderate Resolution Imaging Spectroradiometer (MODIS) on the NASA Terra spacecraft from 2000–2006, statistical analyses are performed to explore the spatial and temporal distribution and variation of the snow cover over the Tibetan Plateau (TP). It is found that the snow persistence over the TP varies in different elevation ranges generally becomes longer with increases in the terrain elevation. In addition, the spatial distribution of the snow cover not only depends on the elevation but also varies with terrain features, such as aspect, slope, and curvature in the local areas. With 7-year observational data, seasonal and interannual variability of snow cover has been detected. There are slight decreasing trends in SFCs from 2000–2006. With MODIS satellite snow-cover fraction data and the National Centers for Environmental Predictions and U.S. Department of Energy NCEP/DOE reanalysis II dataset, the relationship between snow cover anomalies over the TP and the East Asian Summer Monsoon (EASM) is examined. Results indicate that the onset of the EASM is closely associated with snow cover anomalies in the spring. Specifically, a positive (negative) snow cover anomaly is followed by a later (earlier) onset of the EASM.     

Atmospheric blockings in Western Siberia. Part 1. Detection features,objective criteria,and their comparison

Based on the NCEP/NCAR and ERA-Interim reanalysis archives, we investigated different methods to detect atmospheric blocking events in Western Siberia. Two criteria were studied that are based on calculating meridional gradients of the 500 hPa height and potential temperature at the dynamic tropopause. A situation is considered blocking, when it features a gradient inversion of the investigated characteristics. Additionally, we performed a synoptic analysis of individual blocking events.     

An ocean wind-wave climatology for the Southern Brazilian Shelf. Part II: Variability in space and time

Previously validated model results were used to characterize the wave climate over the Southern Brazilian Shelf (SBS). The low mean significant wave height over the western South Atlantic shelves was shown together with examples of cyclone-induced extreme wave fields and other typical wave conditions. The mean offshore spectra showed a bimodal shape with a predominance of S/SSW and ENE/E waves with distinctive interannual rising periods in wave energy density. Along-shelf wave energy gradients were seen near the coast with higher energy located off capes and coastal projections and energy minima between them. A considerable drop in wave energy suggests the 40 m depth as the mean wave base and consequently the lower limit of the SBS shoreface. The upper shoreface mean wave energy density varied abruptly along the shelf in response to differences in bottom declivities. The large and shallow shoreface was responsible for an intense refraction of the waves and hence very small angles of attack. Additionally, it was shown the sheltering effect caused by capes and coastal projections and a remarkable north/south energy asymmetry between them, caused by a windowing on the wave propagation to the shore. Altogether, it was possible to state that bottom friction plays a major role in wave differentiation along the SBS shoreface, thus suggesting that shelf morphology might indeed be more important to generate wave variability than the offshore wave variation itself.     

An ocean wind-wave climatology for the Southern Brazilian Shelf. Part I: Problem presentation and model validation

Geomorphological features and sparse visual wave measurements, presented in previous works, point out to the possible existence of alongshore wave energy gradients over the Southern Brazilian Shelf (SBS). After describing the two numerical models utilized and the model grids and settings, the present study evaluated the proposed validation of the basin-scale results against orbital altimetry and the regional-scale results against the two available wave-buoy data, discussing the extent to which the model reproduced local reality. The basin-scale model results of significant wave height could be considered as in good agreement with observations, presenting remarkable similarities with observed altimetry. Regarding the regional-scale modeling, both significant wave height and peak wave periods were considered as in good agreement with observations. The peak wave directions, however, were classified as in poor agreement, once the ESE waves were erroneously reproduced as E. Increased spatial resolution certainly played a role, but it was the lowering of bottom friction that represented the major improvement in the coastal grid simulations.     

Automated mechanism generation. Part 1: mechanism development and rate constant estimation for VOC chemistry in the atmosphere

A framework for automated mechanism generation for modeling atmospheric chemistry at the mechanistic level was developed. In part 1, categorization of reactions into reaction families and determination of rate coefficients using a hierarchical approach that uses experimental data and kinetic correlations are described. The main focus was to develop kinetic correlations for estimating rate coefficients that are not available experimentally, and the main correlation used was the Evans–Polanyi relationship that relates the activation energy to the heat of reaction. A hierarchical scheme for calculating heats of reaction and other thermodynamic properties was developed. The rate constants calculated using the proposed correlations are in most cases within an order of magnitude of available experimental values, and 82% are within a factor of five.     

High-Latitude Springtime Photochemistry. Part I: NOx, PAN and Ozone Relationships

Springtime measurements of NO_x, ozone, PAN,J(NO₂), and other compounds were made near Ny-Ålesund,Svalbard (78°54N, 11°53E), in 1994 and Poker Flat,Alaska (65°08N, 147°29W), in 1995. At Svalbard medianmixing ratios for PAN and NO_x of 237 and 23.7 pptv,respectively, were observed. The median mixing ratios at Poker Flat for PANand NO_x were 79.5 and 85.9 pptv, respectively. These data areused to estimate thermal PAN decomposition using several differentapproaches. At Svalbard PAN decomposition was very small, while at PokerFlat up to 30 pptv/h PAN decomposed. At both sites the NO_x/PANratio increased with temperature between –10 and 20°C implyingthat PAN decomposition is an important NO_x source. In-situozone production was calculated from the measured NO, NO₂,O₃, J(NO₂), and temperature data, using thesteady state assumption Median ozone production was 605 pptv/h at PokerFlat, and one order of magnitude smaller at Svalbard during the daytime.Only at Poker Flat could a direct influence on the diurnal ozone cycle beobserved from in-situ production. These results imply that PAN decompositionis a major source of NO_x in the high latitude troposphere, andthat this contributes to the observed spring maximum in surface ozone.     

An intercomparison of model-simulated in extreme rainfall and temperature events during the last half of the twentieth century. Part 1: mean values and variability

In this study we examine the performance of eight of the IPCC AR4 global coupled climate models used in the WCRP CMIP3 Multimodel Dataset, as well as their ensemble mean, in simulating annual indices of extreme temperature and precipitation climate events in South America. In this first part we focus on comparing observed and modeled mean values and interannual variability. Two extreme temperature indices based on minimum temperature (warm nights and frost days) and three indices of extreme precipitation (R95t, R10 and consecutive dry days), obtained both from meteorological stations during 1961–2000 and model outputs, were compared. The number of warm nights are better represented by models than the FD. The interannual variability pattern is also in good agreement with the observed values. For precipitation, the index that is best represented by the models is the R95t, which relates the extreme precipitation to local climate. The maximum of dryness observed over the central Argentinian Andes or the extensive dry season of the Amazon region could not be represented by any model.     

An annual cycle of vegetation in a GCM. Part I: implementation and impact on evaporation

The sensitivity of evaporation to a prescribed vegetation annual cycle is examined globally in the Met Office Hadley Centre Unified Model (HadAM3) which incorporates the Met Office Surface Exchange Scheme (MOSES2) as the land surface scheme. A vegetation annual cycle for each plant functional type in each grid box is derived based on satellite estimates of Leaf Area Index (LAI) obtained from the nine-year International Satellite Land Surface Climatology Project II dataset. The prescribed model vegetation seasonality consists of annual cycles of a number of structural vegetation characteristics including LAI as well as canopy height, surface roughness, canopy water capacity, and canopy heat capacity, which themselves are based on empirical relationships with LAI. An annual cycle of surface albedo, which in the model is a function of soil albedo, surface soil moisture, and LAI, is also modelled and agrees reasonably with observed estimates of the surface albedo annual cycle. Two 25-year numerical experiments are completed and compared: the first with vegetation characteristics held at annual mean values, the second with prescribed realistic seasonally varying vegetation. Initial analysis uncovered an unrealistically weak relationship between evaporation and vegetation state that is primarily due to the insensitivity of evapotranspiration to LAI. This weak relationship is strengthened by the adjustment of two MOSES2 parameters that together improve the models LAI-surface conductance relationship by comparison with observed and theoretical estimates. The extinction coefficient for photosynthetically active radiation, k _par , is adjusted downwards from 0.5 to 0.3, thereby enhancing the LAI-canopy conductance relationship. A canopy shading extinction coefficient, k _sh , that controls what fraction of the soil surface beneath a canopy is directly exposed to the overlying atmosphere is increased from 0.5 to 1.0, which effectively reduces soil evaporation under a dense canopy. When the experiments are repeated with the adjusted parameters, the relationship between evaporation and vegetation state is strengthened and is more spatially consistent. At nearly all locations, the annual cycle of evaporation is enhanced in the seasonally varying vegetation experiment. Evaporation is stronger during the peak of the growing season and, in the tropics, reduced transpiration during the dry season when LAI is small leads to significantly lower total evaporation.