Monin–Obukhov Similarity Functions of the Structure Parameter of Temperature and Turbulent Kinetic Energy Dissipation Rate in the Stable Boundary Layer

The Monin–Obukhov similarity theory (MOST) functions f_ε and f_T, of the dissipation rate of turbulent kinetic energy (TKE). ε, and the structure parameter of temperature, C_T², were determined for the stable atmospheric surface layer using data gathered in the context of CASES-99. These data cover a relatively wide stability range, i.e. ζ=z/L of up to 10, where z is the height and L the Obukhov length. The best fits were given by f_ε = 0.8 + 2.5ζ and f_T= 4.7[ 1+1.6(ζ)^2/3], which differ somewhat from previously published functions. ε was obtained from spectra of the longitudinal wind velocity using a time series model (ARMA) method instead of the traditional Fourier transform. The neutral limit f_ε =0.8 implies that there is an imbalance between TKE production and dissipation in the simplified TKE budget equation. Similarly, we found a production-dissipation imbalance for the temperature fluctuation budget equation. Correcting for the production-dissipation imbalance, the ‘standard’ MOST functions for dimensionless wind speed and temperature gradients (φ_m and φ_m) were determined from f_ε and f_T and compared with the φ_m and φ_h formulations of Businger and others. We found good agreement with the Beljaars and Holtslag [J. Appl. Meteorol. 30, 327–341 (1991)] relations. Lastly, the flux and gradient Richardson numbers are discussed also in terms of f_ε and f_T.     

Tropical cyclones and polar lows: Velocity,size, and energy scales,and relation to the 26°C cyclone origin criteria

The goal of this paper is to quantitatively formulate some necessary conditions for the development of intense atmospheric vortices. Specifically, these criteria are discussed for tropical cyclones (TC) and polar lows (PL) by using bulk formulas for fluxes of momentum, sensible heating, and latent heating between the ocean and the atmosphere. The velocity scale is used in two forms: (1) as expressed through the buoyancy flux b and the Coriolis parameter l_c for rotating fluids convection, and (2) as expressed with the cube of velocity times the drag coefficient through the formula for total kinetic energy dissipation in the atmospheric boundary layer. In the quasistationary case the dissipation equals the generation of the energy. In both cases the velocity scale can be expressed through temperature and humidity differences between the ocean and the atmosphere in terms of the reduced gravity, and both forms produce quite comparable velocity scales. Using parameters b and l_c, we can form scales of the area and, by adding the mass of a unit air column, a scale of the total kinetic energy as well. These scales nicely explain the much smaller size of a PL, as compared to a TC, and the total kinetic energy of a TC is of the order 10¹⁸-10¹⁹ J. It will be shown that wind of 33 m s^-1 is produced when the total enthalpy fluxes between the ocean and the atmosphere are about 700 W m^-2 for a TC and 1700 W m^-2 for a PL, in association with the much larger role of the latent heat in the first case and the stricter geostrophic constraints and larger static stability in the second case. This replaces the mystical role of 26^oC as a criterion for TC origin. The buoyancy flux, a product of the reduced gravity and the wind speed, together with the atmospheric static stability, determines the rate of the penetrating convection. It is known from the observations that the formation time for a PL reaching an altitude of 5--6 km can be only a few hours, and a day, or even half a day, for a TC reaching 15--18 km. These two facts allow us to construct curves on the plane of T^s and ΔT=T_s-T_a to determine possibilities for forming an intense vortex. Here, T_a is the atmospheric temperature at the height z=10 m. A PL should have ΔT>20^oC in accordance with the observations and numerical simulations. The conditions for a TC are not so straightforward but our diagram shows that the temperature difference of a few degrees, or possibly even a fraction of a degree, might be sufficient for TC development for a range of static stabilities and development times.     

Evaluation of Haney-Type Surface ThermalBoundary Conditions Using a CoupledAtmosphere and Ocean Model

A coupled atmosphere-ocean model developed at the Institute for Space Studies at NASA Goddard Space Flight Center (Russell et al., 1995) was used to verify the validity of Haney-type surface thermal boundary condition, which linearly connects net downward surface heat flux Q to air / sea temperature difference △T by a relaxation coefficient k. The model was initiated from the National Centers for Environmental Prediction (NCEP) atmospheric observations for 1 December 1977, and from the National Ocean Data Center (NODC) global climatological mean December temperature and salinity fields at 1° ×1° resolution. The time step is 7.5 minutes. We integrated the model for 450 days and obtained a complete model-generated global data set of daily mean downward net surface flux Q, surface air temperature TA,and sea surface temperature To. Then, we calculated the cross-correlation coefficients (CCC) between Q and △T. The ensemble mean CCC fields show (a) no correlation between Q and △T in the equatorial regions, and (b) evident correlation (CCC≥ 0.7) between Q and △T in the middle and high latitudes.Additionally, we did the variance analysis and found that when k= 120 W m-2K-1, the two standard deviations, σQ and σk△T, are quite close in the middle and high latitudes. These results agree quite well with a previous research (Chu et al., 1998) on analyzing the NCEP re-analyzed surface data, except that a smaller value of k (80 W m-2K-1) was found in the previous study.     

Measurement and investigation of chamber radical sources in the European Photoreactor (EUPHORE)

It is essential to quantify the background reactivity of smog-chambers, since this might be the major limitation of experiments carried out at low pollutant concentrations typical of the polluted atmosphere. Detailed investigation of three chamber experiments at zero-NO _x in the European Photoreactor (EUPHORE) were carried out by means of rate-of-production analysis and two uncertainty analysis tools: local uncertainty analysis and Monte Carlo simulations with Latin hypercube sampling. The chemical mechanism employed was that for methane plus the inorganic subset of the Master Chemical Mechanism (MCMv3.1). Newly installed instruments in EUPHORE allowed the measurement of nitrous acid and formaldehyde at sub-ppb concentrations with high sensitivity. The presence of HONO and HCHO during the experiments could be explained only by processes taking place on the FEP Teflon walls. The HONO production rate can be described by the empirical equation W(HONO)_EUPHORE ^dry = a × j _{NO 2}× exp (− T ₀/T) in the low relative humidity region (RH < 2%, a = 7.3×10²¹ cm⁻³, T ₀ = 8945K), and by the equation W(HONO)_EUPHORE ^humid = W(HONO)_EUPHORE ^dry+ j _{NO 2}× b × RH ^q in the higher relative humidity region (2% < RH < 15%, b = 5.8×10⁸ cm⁻³ and q = 0.36, and RH is the relative humidity in percentages). For HCHO the expression W(HCHO)_EUPHORE = c × j _{NO 2}exp (− T′₀/T) is applicable (c = 3.1×10¹⁷ cm⁻³ and T′₀ = 5686 K). In the 0–15% relative humidity range OH production from HONO generated at the wall is about a factor of two higher than that from the photolysis of 100 ppb ozone. Effect of added NO₂ was found to be consistent with the dark HONO formation rate coefficient of MCMv3.1.     

Variance Method to Determine Turbulent Fluxes of Momentum And Sensible Heat in The Stable Atmospheric Surface Layer

Evidence is presented that in the stable atmospheric surface layer turbulent fluxes of heat and momentum can be determined from the standard deviations of longitudinal wind velocity and temperature, σ _u and σ _T respectively, measured at a single level. An attractive aspect of this method is that it yields fluxes from measurements that can be obtained with two-dimensional sonic anemometers. These instruments are increasingly being used at official weather stations, where they replace the standard cup anemometer–wind vane system. With methods such as the one described in this note, a widespread, good quality, flux network can be established, which would greatly benefit the modelling community. It is shown that a ‘variance’ dimensionless height (ζ _σ) defined from σ _u and σ _T is highly related to the ‘conventional’ dimensionless stability parameter ζ=z/L, where z is height and L is the Obukhov length. Empirical functions for ζ _σ are proposed that allow direct calculation of heat and momentum fluxes from σ _u and σ _T. The method performs fairly well also during a night of intermittent turbulence.     

A Study of the Reaction of OH Radicals with N-Methyl Pyrrolidinone,N-Methyl Succinimide and N-Formyl Pyrrolidinone

A study of the oxidation mechanism of N-methyl pyrrolidinone (C₅H₉NO, NMP) initiated by hydroxyl radicals was made at EUPHORE at atmospheric pressure (1000 ± 10) mbar of air and ambient temperature (T = 300 ± 5 K). The main products were N-methyl succinimide (NMS) (52 ± 4)% and N-formyl pyrrolidinone (FP) (23 ± 9)%. The relative rate technique was used to determine the rate constants of OH with NMP, NMS and FP, the measured values were (in units of cm³ molecule ^{− 1} s^{− 1}): k_NMP = (2.2 ± 0.4) × 10^{− 11}, k_NMS = (1.4 ± 0.3) × 10^{− 12} and k_FP = (6 ± 1) × 10^{− 12}. The results are presented and discussed in terms of the atmospheric impact.     

Buoyancy and The Sensible Heat Flux Budget Within Dense Canopies

In contrast to atmospheric surface-layer (ASL) turbulence, a linear relationship between turbulent heat fluxes (F_T) and vertical gradients of mean air temperature within canopies is frustrated by numerous factors, including local variation in heat sources and sinks and large-scale eddy motion whose signature is often linked with the ejection-sweep cycle. Furthermore, how atmospheric stability modifies such a relationship remains poorly understood, especially in stable canopy flows. To date, no explicit model exists for relating F_T to the mean air temperature gradient, buoyancy, and the statistical properties of the ejection-sweep cycle within the canopy volume. Using third-order cumulant expansion methods (CEM) and the heat flux budget equation, a “diagnostic” analytical relationship that links ejections and sweeps and the sensible heat flux for a wide range of atmospheric stability classes is derived. Closure model assumptions that relate scalar dissipation rates with sensible heat flux, and the validity of CEM in linking ejections and sweeps with the triple scalar-velocity correlations, were tested for a mixed hardwood forest in Lavarone, Italy. We showed that when the heat sources (S_T) and F_T have the same sign (i.e. the canopy is heating and sensible heat flux is positive), sweeps dominate the sensible heat flux. Conversely, if S_T and F_T are opposite in sign, standard gradient-diffusion closure model predict that ejections must dominate the sensible heat flux.     

Skill as a function of time scale in ensembles of seasonal hindcasts

Forecast skill as a function of time lead and time averaging is examined in two 6-member ensembles of seasonal hindcasts. One ensemble is produced with the second generation general circulation model of the Canadian Centre for Climate Modelling and Analysis (GCM2) and the other with a reduced resolution version of the numerical weather prediction model of the Canadian Meteorological Centre (SEF). The integrations are initiated from the NCEP/NCAR reanalyzed data. Monthly sea surface temperature anomalies observed prior to the forecast period are maintained throughout the forecast season. A statistical forecast improvement technique, based on the singular value decomposition of forecast and reanalyzed fields, is discussed and evaluated. A simple analogue of the hindcast integrations is used to examine the behavior of two common skill scores, the correlation skill score and the explained variance skill score. The maximal skill score and the corresponding optimal forecast in this analogue are identified. The total skill of the optimal forecast is a sum of two terms, one associated with the initial conditions and the other with the lower boundary forcing. The two sources of skill operate on different time scales, with initial conditions being more important in the first one-two weeks and the atmospheric response to the boundary forcing becoming more dominant for longer time leads and time averages. This suggests that these sources of skill should be considered separately in forecast optimization. The statistical technique is moderately successful in improving the skill of monthly to seasonal forecasts of 500 hPa height (Z ₅₀₀) and 700 hPa temperature (T ₇₀₀) in the Northern Hemisphere and in the North Pacific/North America sector. The improvement is better when the forecasts for the first week and for the rest of the season are optimized separately. The SEF model produces better Z ₅₀₀ and T ₇₀₀ forecasts than GCM2 in the first one-two weeks whereas GCM2 performs slightly better at longer time leads. The skill of zero time lead forecast decays rapidly with averaging interval for time averages up to about 30–45 days and stabilizes, or even rises, for longer time averages. Excluding the first week from seasonal forecasts results in substantial degradation of predictive skill. Received: 1 November 1999 / Accepted: 24 May 2000     

梅汛期区域性暴雨的多尺度分析及临近预警

利用常规观测数据、NECP1°×1°逐6 h再分析资料、FY2E卫星黑体亮度温度TBB资料以及南京、常州多普勒天气雷达产品,对2012年江苏出梅之前的最后一场区域性暴雨过程进行多尺度分析。在此基础上,探索该类暴雨的临近预警线索,结果表明:(1)此次过程的雨带呈准纬向分布,属于典型的梅汛期静止锋降水。过程中主要有两次降水集中时段,两个阶段的降水性质存在差异,但都具备较高的降水效率。(2)中高纬大范围稳定的阻塞形势,为此次持续性暴雨过程的产生提供了有利的大尺度环流背景。而在此过程中,两段降水集中期的形成与地面触发系统的出现和维持有着较为密切的联系。(3)此次过程中两段降水集中期内的物理量特征以及TBB的演变情况和其对应的降水特征存在异同。(4)雷达特征分析表明,此次过程具有较高的降水效率和较长的持续时间。实际业务工作中可以通过判断回波的质心高度和边界层风速有无跃增来估计降水效率的潜势。当推断较高的降水效率潜势将持续较长时间时,应及时发布暴雨警报。     

An analysis of the turbulent structure in the unstable surface layer nearby a shelter belt

An analysis was performed of the turbulent data obtained from Yucheng experimental station in the Shandong Province in 1984. [t is shown that at variant wind speed, the spectra of streamwise velocity remain similar and the intensity of wind fluctuations is proportional to wind speed in the downwind area of shelter belt. Therefore, we may decide the similarity of wind fluctuations by a speed scale and a length scale which is not correlated with stability, σu /V0 = F(X / H). The -5/3 power range of temperature spectra extends to lower frequency. The variation of ratio σ0 /T. with stability becomes σ0 / T . = C(X / H)( - Z / L)-1/3 . There is not such an extension of -5 / 3 power range in the humidity spectra.     

Quantitative reconstruction of the last interglacial vegetation and climate based on the pollen record from Lake Baikal, Russia

Changes in mean temperature of the coldest (T _c) and warmest month (T _w), annual precipitation (P _ann) and moisture index (α) were reconstructed from a continuous pollen record from Lake Baikal, Russia. The pollen sequence CON01-603-2 (53°57′N, 108°54′E) was recovered from a 386 m water depth in the Continent Ridge and dated to ca. 130–114.8 ky BP. This time interval covers the complete last interglacial (LI), corresponding to MIS 5e. Results of pollen analysis and pollen-based quantitative biome reconstruction show pronounced changes in the regional vegetation throughout the record. Shrubby tundra covered the area at the beginning of MIS 5e (ca. 130–128 ky), consistent with the end of the Middle Pleistocene glaciation. The late glacial climate was characterised by low winter and summer temperatures (T _c ~ −38 to −35°C and T _w~11–13°C) and low annual precipitation (P _ann~300 mm). However, the wide spread of tundra vegetation suggests rather moist environments associated with low temperatures and evaporation (reconstructed α~1). Tundra was replaced by boreal conifer forest (taiga) by ca. 128 ky BP, suggesting a transition to the interglacial. Taiga-dominant phase lasted until ca. 117.4 ky BP, e.g. about 10 ky. The most favourable climate conditions occurred during the first half of the LI. P _ann reached 500 mm soon after 128 ky BP. However, temperature changed more gradually. Maximum values of T _c ~ −20°C and T _w~16–17°C are reconstructed from about 126 ky BP. Conditions became gradually colder after ca. 121 ky BP. T _c dropped to ~ −27°C and T _w to ~15°C by 119.5 ky BP. The reconstructed increase in continentality was accompanied by a decrease in P _ann to ~400–420 mm. However, the climate was still humid enough (α~0.9) to support growth of boreal evergreen conifers. A sharp turn towards a dry climate is reconstructed after ca. 118 ky BP, causing retreat of forest and spread of cool grass-shrub communities. Cool steppe dominated the vegetation in the area between ca. 117.5 ky and 114.8 ky BP, suggesting the end of the interglacial and transition to the last glacial (MIS 5d). Shift to the new glaciation was characterised by cool and very dry conditions with T _c ~ −28 to −30°C, T _w~14–15°C, P _ann~250 mm and α~0.5.     

Characteristics of the atmospheric heat source and moisture sink over the Qinghai-Tibetan Plateau during the second TIPEX of summer 1998 and their impact on surrounding monsoon

Summary ?Using the data of 6 automatic heat balance observation (AWS) stations and a data set of 52 surface observation stations over the Qinghai-Tibetan Plateau (“the Plateau”) and surroundings, the horizontal distribution is studied of “apparent atmospheric heat sources” 〈Q ₁〉 and of “apparent atmospheric moisture sinks” 〈Q ₂〉. The AWS stations were established during the period May to August 1998 of the Tibetan Plateau Meteorological Experiment (second TIPEX) by a cooperation of China and Japan. For this period the Plateau mean of 〈Q ₁〉 is positive. Its value of 74 W/m² is a little greater than a climate value and than values from MONEX and the first TIPEX in 1979, respectively. Also the corresponding 〈Q ₂〉 is positive. Hence during that time the Plateau is a heat source and a moisture sink. A day-to-day change of 〈Q ₁〉 and 〈Q ₂〉 is more pronounced over the middle and east part of the Plateau than over the west part. Diagnostics accompanied by numerical simulations are used to study the daily relationship between 〈Q ₁〉 over the Plateau and the weather over China and Asia for this summer. The results suggest that 〈Q ₁〉 may affect precipitation over northern China and position of the west Pacific subtropical high. Abnormal southward retreat of this Pacific high seems to have caused the second flood over the middle and lower Yangtse river basin in July. Received May 20, 2001; revised February 2, 2002     

Study on horizontal relative diffusion in the troposphere and lower stratosphere

The behaviour of relative diffusion theory and Gifford’s random-force theory for long-range atmospheric diffusion is examined. When a puff scale is smaller than the Lagrangian length scale, 2KTL, an accelerative relative diffusion region exists, i.e., σy∝t3/2. While the puff diffusion enters a two-dimensional turbulence region, in which the diffusion scale is larger than 500 km, or time scale is larger than 1 day, divergence and convergence are main cause of horizontal diffusion. Between the two above-mentioned regimes, diffusion deviation is given byσy=2KTL. The large-scale horizontal relative diffusion parameters were obtained by analyzing the data of radioactive cloud width collected in air nuclear tests.     

校正热带测雨卫星轨道抬升对微波成像仪亮温的影响

热带测雨卫星(Tropical Rainfall Measuring Mission,TRMM)在2001年8月轨道高度从350 km升高至402 km,搭载于其上的微波成像仪(TRMM Microwave Imager,TMI)的入射角随之发生变化,进而使得相应探测结果(亮温)发生改变,从而导致由此反演的大气参数出现虚假的突变。为保证轨道抬升前后TMI亮温资料的连致性,以便更好地用于气候研究,本研究首先分析了洋面轨道抬升前后亮温的差异及变化原因,然后结合微波辐射传输模式,分析了不同环境参数对亮温变化的影响,在此基础上用线性变换的方式对轨道抬升后的亮温进行了修正,并从不同角度检验修正效果。结果表明,轨道抬升前后亮温呈线性关系,低频垂直极化通道亮温轨道抬升后升高了0.8-1.6 K,其他通道亮温变化不大。经过修正,轨道抬升前后的亮温趋于一致,月平均亮温偏差明显减小,低频垂直极化通道亮温在轨道抬升期间的突变被消除,亮温变得连续平稳,可用于气候研究。     

The Theoretical Model of Atmospheric Turbulence Spectrum in Surface Layer

It is shown that the slope of energy spectrum obtained from the velocity solution of Kdv-Burgers equation lies between -5/3 and -2 in the dilogarithmic coordinates paper. The spectrum is very close to one of Kolmogorov’s isotropic turbulence and Frisch’s intermittent turbulence in inertial region. In this paper, the Kdv-Burgers equation to describe atmospheric boundary layer turbulence is obtained. In the equation, the 1 / Re, corresponds to dissipative coefficient v, to dispersive coefficient β, then (v/ 2β)2 corresponds to .We prove that the wave number corresponding to maximum energy spectrum decreases with the decrease of stability (i.e., the increase of in eddy-containing region. And the spectrim amplitude decreases with the increase of (i.e., the decrease of stability). These results are consistent with actual turbulence spectrum of atmospheric surface layer from turbulence data.     

Humidity Effect and Its Influence on the Seasonal Distribution of Precipitation δ18O in Monsoon Regions

The humidity effect, namely the markedly positive correlation between the stable isotopic ratio in precipitation and the dew-point deficit △Td in the atmosphere, is put forward firstly and the relationships between the δ18O in precipitation and △Td are analyzed for the Urumqi and Kunming stations, which have completely different climatic characteristics. Although the seasonal variations in δ18O and △Td exhibit differences between the two stations, their humidity effect is notable. The correlation coefficient and its confidence level of the humidity effect are higher than those of the amount effect at Kunming, showing the marked influence of the humidity conditions in the atmosphere on stable isotopes in precipitation.Using a kinetic model for stable isotopic fractionation, and according to the seasonal distribution of meanmonthly temperature at 500 hPa at Kunming, the variations of the δ18O in condensate in cloud aresimulated. A very good agreement between the seasonal variations of the simulated mean δ18O and themean monthly temperature at 500 hPa is obtained, showing that the oxygen stable isotope in condensateof cloud experiences a temperature effect. Such a result is markedly different from the amount effect atthe ground. Based on the simulations of seasonal variations of δ18O in falling raindrops, it can be foundthat, in the dry season from November to April, the increasing trend with falling distance of δ18O in fallingraindrops corresponds remarkably to the great ATd, showing a strong evaporation enrichment function infalling raindrops; however, in the wet season from May to October, the δ18O in falling raindrops displaysan unapparent increase corresponding to the small ATd, except in May. By comparing the simulated meanδ18O at the ground with the actual monthly δ18O in precipitation, we see distinctly that the two monthlyδ18O variations agree very well. On average, the δ18O values are relatively lower because of the highlymoist air, heavy rainfall, small △Td and weak evaporation enrichment function of stable isotopes in thefalling raindrops, under the influence of vapor from the oceans; but they are relatively higher because of the dry air, light rainfall, great △Td and strong evaporation enrichment function in falling raindrops, under the control of the continental air mass. Therefore, the δ18O in precipitation at Kunming can be used to indicate the humidity situation in the atmosphere to a certain degree, and thus indicate the intensity of the precipitation and the strength of the monsoon indirectly. The humidity effect changes not only the magnitude of the stable isotopic ratio in precipitation but also its seasonal distribution due to its influence on the strength of the evaporation enrichment of stable isotopes in falling raindrops and the direction of the net mass transfer of stable isotopes between the atmosphere and the raindrops. Consequently, it is inferred that the humidity effect is probably one of the foremost causes generating the amount effect.     

Atmospheric turbidity of urban and desert areas of the Nile Basin in the aftermath of Mt. Pinatubo’s eruption

Summary  The Linke T _L , ?ngstr?m β and Unsworth-Monteith δ _a turbidity parameters are investigated for two sites in Egypt: Cairo, a densely populated urban area, and Aswan, an arid unpolluted area. These three turbidity parameters are calculated from broadband pyrheliometric measurements recorded hourly over the period 1992–96. Monthly averages of T _L , β and δ _a show relatively flat and identical seasonal variations with a marked main maxima during spring at both sites, due to Khamsin depressions coming from the Great Sahara. A secondary maximum is observed at Aswan in summer, due to dust haze which prevails during that season, and at Cairo in autumn, due to the northern extension of the Sudan monsoon trough, which is accompanied by small scale depressions with dust particles. Annual mean values of T _L , β and δ _a (5.59, 0.250 and 0.372, respectively) at Cairo are larger than at Aswan (3.89, 0.139 and 0.213, respectively). In the same way, the seasonal mean values of T _L , β and δ _a at Cairo are larger than at Aswan. More generally, the monthly and yearly average turbidity values are significantly larger in Cairo than in Aswan for the whole period 1992–96, which is attributable in part to the urbanization/industrialization effect of Cairo. An estimate of the corresponding overburden is obtained by comparison between the present data and older T _L data from 1922–27. It is also shown that turbidity over both sites is largest during 1992, just after the eruption of Mount Pinatubo in 1991. The dependence of β on some meteorological parameters such as wind speed and direction, precipitable water, relative humidity, temperature and visibility, is also analyzed. This reveals in particular that visibility is not a good predictor of turbidity at either site. Conversely, the wind direction and speed have a definite effect on turbidity, and consequently, largest turbidities occur when the wind carries aerosols from the main industrial particle source areas around Cairo. For any season of the year, the average turbidity at the latter site is larger than that at other big cities such as Athens, Rome, and Toronto, but is lower than at Dhahran, Saudi Arabia. Received February 3, 2000 Revised August 25, 2000     

Simulated changes in vegetation distribution, land carbon storage, and atmospheric CO2 in response to a collapse of the North Atlantic thermohaline circulation

It is investigated how abrupt changes in the North Atlantic (NA) thermohaline circulation (THC) affect the terrestrial carbon cycle. The Lund–Potsdam–Jena Dynamic Global Vegetation Model is forced with climate perturbations from glacial freshwater experiments with the ECBILT-CLIO ocean–atmosphere–sea ice model. A reorganisation of the marine carbon cycle is not addressed. Modelled NA THC collapses and recovers after about a millennium in response to prescribed freshwater forcing. The initial cooling of several Kelvin over Eurasia causes a reduction of extant boreal and temperate forests and a decrease in carbon storage in high northern latitudes, whereas improved growing conditions and slower soil decomposition rates lead to enhanced storage in mid-latitudes. The magnitude and evolution of global terrestrial carbon storage in response to abrupt THC changes depends sensitively on the initial climate conditions. These were varied using results from time slice simulations with the Hadley Centre model HadSM3 for different periods over the past 21 kyr. Changes in terrestrial storage vary between −67 and +50 PgC for the range of experiments with different initial conditions. Simulated peak-to-peak differences in atmospheric CO₂ are 6 and 13 ppmv for glacial and late Holocene conditions. Simulated changes in δ¹³C are between 0.15 and 0.25‰. These simulated carbon storage anomalies during a NA THC collapse depend on their magnitude on the CO₂ fertilisation feedback mechanism. The CO₂ changes simulated for glacial conditions are compatible with available evidence from marine studies and the ice core CO₂ record. The latter shows multi-millennial CO₂ variations of up to 20 ppmv broadly in parallel with the Antarctic warm events A1 to A4 in the South and cooling in the North.     

热带太平洋表面水中CO2对El Niño事件响应的数值模拟

文中用一个带生物泵的三维全球海洋碳循环模式模拟了热带太平洋表面水中CO2总量(TCO2)在ElNi     

Determination of the distance between two adjacent stations,the observational vertical increment and the observational time interval in optimum sense

Considering the observational error, the truncation error and the requirements of numerical weather prediction, three formulas for determining the distance between two adjacent stations d1, the observational vertical increment △p1 and the observational time interval △t1 in optimum sense, have been derived. Since they depend on the shortest wavelength concerned and the ratio of maximum observational error to wave amplitude, the results are quite different for different scale systems.For the filtered model the values of d1, △p1,, and △t1 in general come near those required in the MANUAL on the GOS published in 1980 by WMO. But for the primitive equation model the estimated value of △t1 is much less than those required in the filtered model case.Therefore, it is improper to study the fast moving and developing processes of the atmospheric motion only on the basis of the conventional observations. It seems to be necessary to establish an optimum composite observational system including the surface-based system and the space-based system.