Stable Boundary-Layer Scaling Regimes: The Sheba Data

Turbulent and mean meteorological data collected at five levels on a 20-m tower over the Arctic pack ice during the Surface Heat Budget of the Arctic Ocean experiment (SHEBA) are analyzed to examine different regimes of the stable boundary layer (SBL). Eleven months of measurements during SHEBA cover a wide range of stability conditions, from the weakly unstable regime to very stable stratification. Scaling arguments and our analysis show that the SBL can be classified into four major regimes: (i) surface-layer scaling regime (weakly stable case), (ii) transition regime, (iii) turbulent Ekman layer, and (iv) intermittently turbulent Ekman layer (supercritical stable regime). These four regimes may be considered as the basic states of the traditional SBL. Sometimes these regimes, especially the last two, can be markedly perturbed by gravity waves, detached elevated turbulence (‘upside down SBL’), and inertial oscillations. Traditional Monin–Obukhov similarity theory works well in the weakly stable regime. In the transition regime, Businger–Dyer formulations work if scaling variables are re-defined in terms of local fluxes, although stability function estimates expressed in these terms include more scatter compared to the surface-layer scaling. As stability increases, the near-surface turbulence is affected by the turning effects of the Coriolis force (the turbulent Ekman layer). In this regime, the surface layer, where the turbulence is continuous, may be very shallow (< 5 m). Turbulent transfer near the critical Richardson number is characterized by small but still significant heat flux and negligible stress. The supercritical stable regime, where the Richardson number exceeds a critical value, is associated with collapsed turbulence and the strong influence of the earth’s rotation even near the surface. In the limit of very strong stability, the stress is no longer a primary scaling parameter.     

Did Humankind Prevent a Holocene Glaciation?

Recently, W.F. Ruddiman (2003, Climatic Change, Vol. 61, pp. 261–293) suggested that the anthropocene, the geological epoch of significant anthropospheric interference with the natural Earth system, has started much earlier than previously thought (P. I. Crutzen and E. F. Stoermer, 2000, IGBP Newsletter, Vol. 429, pp. 623–628). Ruddiman proposed that due to human land use, atmospheric concentrations of CO₂ and CH₄ began to deviate from their natural declining trends some 8000 and 5000 years ago, respectively. Furthermore, Ruddiman concluded that greenhouse gas concentrations grew anomalously thereby preventing natural large-scale glaciation of northern North America that should have occurred some 4000–5000 years ago without human interference. Here we would like to comment on (a) natural changes in atmospheric CO₂ concentration during the Holocene and (b) on the possibility of a Holocene glacial inception. We substantiate our comments by modelling results which suggest that the last three interglacials are not a proper analogue for Holocene climate variations. In particular, we show that our model does not yield a glacial inception during the last several thousand years even if a declining trend in atmospheric CO₂ was assumed.     

The earth system model of intermediate complexity CLIMBER-3α. Part I: description and performance for present-day conditions

We herein present the CLIMBER-3α Earth System Model of Intermediate Complexity (EMIC), which has evolved from the CLIMBER-2 EMIC. The main difference with respect to CLIMBER-2 is its oceanic component, which has been replaced by a state-of-the-art ocean model, which includes an ocean general circulation model (GCM), a biogeochemistry module, and a state-of-the-art sea-ice model. Thus, CLIMBER-3α includes modules describing the atmosphere, land-surface scheme, terrestrial vegetation, ocean, sea ice, and ocean biogeochemistry. Owing to its relatively simple atmospheric component, it is approximately two orders of magnitude faster than coupled GCMs, allowing the performance of a much larger number of integrations and sensitivity studies as well as longer ones. At the same time its oceanic component confers on it a larger degree of realism compared to those EMICs which include simpler oceanic components. The coupling does not include heat or freshwater flux corrections. The comparison against the climatologies shows that CLIMBER-3α satisfactorily describes the large-scale characteristics of the atmosphere, ocean and sea ice on seasonal timescales. As a result of the tracer advection scheme employed, the ocean component satisfactorily simulates the large-scale oceanic circulation with very little numerical and explicit vertical diffusion. The model is thus suited for the study of the large-scale climate and large-scale ocean dynamics. We herein describe its performance for present-day boundary conditions. In a companion paper (Part II), the sensitivity of the model to variations in the external forcing, as well as the role of certain model parameterisations and internal parameters, will be analysed.     

Time-scale and state dependence of the carbon-cycle feedback to climate

Climate and atmospheric CO₂ concentration are intimately coupled in the Earth system: CO₂ influences climate through the greenhouse effect, but climate also affects CO₂ through its impact on the amount of carbon stored on land and in the ocean. The change in atmospheric CO₂ as a response to a change in temperature ( $\varDelta CO_{2}/\varDelta T$ ) is a useful measure to quantify the feedback between the carbon cycle and climate. Using an ensemble of experiments with an Earth system model of intermediate complexity we show a pronounced time-scale dependence of $\varDelta CO_{2}/\varDelta T$ . A maximum is found on centennial scales with $\varDelta CO_{2}/\varDelta T$ values for the model ensemble in the range 5–12 ppm °C^?1, while lower values are found on shorter and longer time scales. These results are consistent with estimates derived from past observations. Up to centennial scales, the land carbon response to climate dominates the CO₂ signal in the atmosphere, while on longer time scales the ocean becomes important and eventually dominates on multi-millennial scales. In addition to the time-scale dependence, modeled $\varDelta CO_{2}/\varDelta T$ show a distinct dependence on the initial state of the system. In particular, on centennial time-scales, high $\varDelta CO_{2}/\varDelta T$ values are correlated with high initial land carbon content. A similar relation holds also for the CMIP5 models, although for $\varDelta CO_{2}/\varDelta T$ computed from a very different experimental setup. The emergence of common patterns like this could prove to usefully constrain the climate–carbon cycle feedback.     

The Critical Richardson Number and Limits of Applicability of Local Similarity Theory in the Stable Boundary Layer

Measurements of atmospheric turbulence made over the Arctic pack ice during the Surface Heat Budget of the Arctic Ocean experiment (SHEBA) are used to determine the limits of applicability of Monin–Obukhov similarity theory (in the local scaling formulation) in the stable atmospheric boundary layer. Based on the spectral analysis of wind velocity and air temperature fluctuations, it is shown that, when both the gradient Richardson number, Ri, and the flux Richardson number, Rf, exceed a ‘critical value’ of about 0.20–0.25, the inertial subrange associated with the Richardson–Kolmogorov cascade dies out and vertical turbulent fluxes become small. Some small-scale turbulence survives even in this supercritical regime, but this is non-Kolmogorov turbulence, and it decays rapidly with further increasing stability. Similarity theory is based on the turbulent fluxes in the high-frequency part of the spectra that are associated with energy-containing/flux-carrying eddies. Spectral densities in this high-frequency band diminish as the Richardson–Kolmogorov energy cascade weakens; therefore, the applicability of local Monin–Obukhov similarity theory in stable conditions is limited by the inequalities Ri < Ri _cr and Rf < Rf _cr. However, it is found that Rf _cr  =  0.20–0.25 is a primary threshold for applicability. Applying this prerequisite shows that the data follow classical Monin–Obukhov local z-less predictions after the irrelevant cases (turbulence without the Richardson–Kolmogorov cascade) have been filtered out.     

Systematic Errors of High–Precision Photometric Catalogues

A very accurate imitation of Hipparcos and Tycho Hp, B _T, andV _T magnitudes was made using W, B, V, R magnitudes from the Tien Shan photometric catalogue. The calculated magnitudes were compared to the observed ones. It is shown that there are systematic differences between calculated and observed magnitudes. The systematic errors are supposed to be bound up with the sky scanning procedure on the Hipparcos satellite. Polynomials in powers of coordinates have been proposed to take into account the systematic errors. 6558 stars have been found to be appropriate high-precision photometric standards. This revised version was published online in July 2006 with corrections to the Cover Date.     

Intergranular diamonds derived from partial melting of crustal rocks at ultrahigh-pressure metamorphic conditions

By reporting for the first time intergranular diamond in quartz–feldspar (Qtz–Kfs) aggregates, the processes of metamorphic diamond formation have to be reconsidered. Based on their Kfs/Qtz ratio, the texture of these aggregates are proposed to result from ‘granitic’ melt with a calculated composition that corresponds well with that of experimental data for the pelitic system. Taking into account experiments on CO₂ solubility in silicate melt under ultrahigh‐pressure conditions, a granitic melt is further suggested to act as a crystallization medium as well as a transport medium for producing metamorphic diamond.     

Southern East Siberia Pliocene–Quaternary faults: Database, analysis and inference

This paper presents the first release of an Informational System（IS）devoted to the systematic collection of all available data relating to Pliocene-Quaternary faults in southern East Siberia,their critical analysis and their seismotectonic parameterization.The final goal of this project is to form a new base for improving the assessment of seismic hazard and other natural processes associated with crustal deformation.The presented IS has been exploited to create a relational database of active and conditionally active faults in southern East Siberia（between 100°-114° E and 50°-57° N）whose central sector is characterized by the highly seismic Baikal rift zone.The information within the database for each fault segment is organized as distinct but intercorrelated sections（tables,texts and pictures,etc.）and can be easily visualized as HTML pages in offline browsing.The preliminary version of the database distributed free on disk already highlights the general fault pattern showing that the Holocene and historical activity is quite uniform and dominated by NE-SW and nearly E-W trending faults;the former with a prevailing dip-slip normal kinematics,while the latter structures are left-lateral strike-slip and oblique-slip（with different proportion of left-lateral and normal fault slip components）.These faults are mainly concentrated along the borders of the rift basins and are the main sources of moderate-to-strong（M≥5.5）earthquakes on the southern sectors of East Siberia in recent times.As a whole,based on analyzing the diverse fault kinematics and their variable spatial distribution with respect to the overall pattern of the tectonic structures formed and/or activated during the late Pliocene-Quaternary,we conclude they were generated under a regional stress field mainly characterized by a relatively uniform NW-SE tension,but strongly influenced by the irregular hard boundary of the old Siberian craton.The obtained inferences are in an agreement with the existing models of the development of