A ground hydrologic model with inclusion of a layer of vegetation canopy that can interact with general-circulation model

In this paper, a ground hydrologic model(GHM) is presented in which the vapor, heat and momentum exchanges between ground surface covers (including vegetation canopy) and atmosphere is described more realistically. The model is used to simulate three sets of field data and results from the numerical simulation agree with the field data well. GHM has been tested using input data generated by general circulation model (GCM) runs for both the North American regions and the Chinese regions, The results from GHM are quite different from those of GHMs in GCMs. It shows that a more active concerted effort on the land surface process study to provide a physically realistic GHM for predicting the exchange between land and atmosphere is important and necessary.

未来气候变化对黄土高原黑河流域水资源的影响

气候变化对黄土高原的水资源有重要影响,对其影响进行评估可以为区域发展提供重要的决策依据。基于分布式水文模型SWAT和4种全球环流模式的各3种排放情景,评估了2010～2039年黄土高塬沟壑区黑河流域水资源对气候变化的潜在响应。结果表明,黑河流域2010～2039年的年均降水变化-2.3%～7.8%,年均最高和最低温度分别升高0.7～2.2℃和1.2～2.8℃,年均径流量变化-19.8%～37.0%,1.2m剖面年均土壤水分含量变化-5.5%～17.2%,年均蒸散量普遍增长0.1%～5.9%;水文气象变量变化趋势复杂,但T检验表明年降水、径流、土壤水分和蒸散增长的概率较大。对于季节变化,降水可能在12～7月份和9月份增长,8月份和10～11月份减少;径流在4～7月份和9～10月份增加,11～3月份和8月份减少;土壤水分在各月都增长;蒸散11～6月份普遍增长,7～10月份减少的可能性较大。未来气候将发生显著变化并对水资源有重要影响,需采取必要的措施来减缓其不利影响。

A global shallow water model using high order multi-moment constrained finite volume method and icosahedral grid

A novel accurate numerical model for shallow water equations on sphere have been developed by implementing the high order multi-moment constrained finite volume (MCV) method on the icosahedral geodesic grid. High order reconstructions are conducted cell-wisely by making use of the point values as the unknowns distributed within each triangular cell element. The time evolution equations to update the unknowns are derived from a set of constrained conditions for two types of moments, i.e. the point values on the cell boundary edges and the cell-integrated average. The numerical conservation is rigorously guaranteed. in the present model, all unknowns or computational variables are point values and no numerical quadrature is involved, which particularly benefits the computational accuracy and efficiency in handling the spherical geometry, such as coordinate transformation and curved surface. Numerical formulations of third and fourth order accuracy are presented in detail. The proposed numerical model has been validated by widely used benchmark tests and competitive results are obtained. The present numerical framework provides a promising and practical base for further development of atmospheric and oceanic general circulation models. (C) 2009 Elsevier Inc. All rights reserved.

Intraseasonal variability of Indian Ocean sea surface temperature during boreal winter: Madden-Julian Oscillation versus submonthly forcing and processes

[ 1] Intraseasonal variability of Indian Ocean sea surface temperature (SST) during boreal winter is investigated by analyzing available data and a suite of solutions to an ocean general circulation model for 1998 - 2004. This period covers the QuikSCAT and Tropical Rainfall Measuring Mission (TRMM) observations. Impacts of the 30 - 90 day and 10 - 30 day atmospheric intraseasonal oscillations (ISOs) are examined separately, with the former dominated by the Madden-Julian Oscillation (MJO) and the latter dominated by convectively coupled Rossby and Kelvin waves. The maximum variation of intraseasonal SST occurs at 10 degrees S - 2 degrees S in the wintertime Intertropical Convergence Zone (ITCZ), where the mixed layer is thin and intraseasonal wind speed reaches its maximum. The observed maximum warming ( cooling) averaged over ( 60 degrees E - 85 degrees E, 10 degrees S - 3 degrees S) is 1.13 degrees C ( - 0.97 degrees C) for the period of interest, with a standard deviation of 0.39 degrees C in winter. This SST change is forced predominantly by the MJO. While the MJO causes a basin-wide cooling ( warming) in the ITCZ region, submonthly ISOs cause a more complex SST structure that propagates southwestward in the western-central basin and southeastward in the eastern ocean. On both the MJO and submonthly timescales, winds are the deterministic factor for the SST variability. Short-wave radiation generally plays a secondary role, and effects of precipitation are negligible. The dominant role of winds results roughly equally from wind speed and stress forcing. Wind speed affects SST by altering turbulent heat fluxes and entrainment cooling. Wind stress affects SST via several local and remote oceanic processes.

太阳辐射对热层和电离层变化性的影响

Solar ultraviolet (UV) radiation at wavelengths less than 400 nm is an important source of energy for aeronomic processes throughout the solar system. Solar UV photons are absorbed in planetary atmospheres, as well as throughout the heliosphere, via photodissociation of molecules, photoionization of molecules and atoms, and photoexcitation toexcitation including resonance scattering. In this paper, the solar irradiances data measured by TIMED SEE, as well as the solar proxies such as F10.7 and Mg II, thermosphere neutral density of CHAMP measurements and topside ionospheric plasmas densities from DMSP, are used to analyze solar irradiance effects on the variabilities of the thermosphere and the ionosphere. First, thermosphere densities near 410 km altitude are analyzed for solar irradiance variability effects during the period 2002-2004. Correlations between the densities and the solar irradiances for different spectral lines and wavelength ranges reveal significantly different characteristics. The density correlates remarkably well with all the selected solar irradiances except the lower chromospheric O I (130.4 nm) emission. Among the chosen solar proxies, the Mg II core-to-wing ratio index, EUV (30-120 nm) and F10.7 show the highest correlations with the density for short-term (< ~27 days) variations. For both long- (> ~27 days) and short-term variations, linear correlation coefficients exhibit a decreasing trend from low latitudes towards high latitudes. The density variability can be effectively modeled (capturing 71% of the variance) using multiple solar irradiance indices, including F10.7, SEUV (the EUV 30-120 nm index), and SFUV (the FUV 120-193 nm index), in which a lag time of 1 day was used for both F10.7 and SEUV, and 5 days for SFUV. In our regression formulation SEUV has the largest contribution to the density variation (40%), with the F10.7 having the next largest contribution (32%) and SFUV accounting for the rest (28%). Furthermore, a pronounced period of about 27.2 days (mean period of the Sun's rotation) is present in both density and solar irradiance data of 2003 and 2004, and a pronounced period of about 54.4 days (doubled period of the solar rotation) is also revealed in 2004. However, soft X-ray and FUV irradiances did not present a pronounced 54.4 day period in 2004, in spite of their high correlation with the densities. The Ap index also shows 54-day periodicities in 2004, and magnetic activity, together with solar irradiance, affects the 54-day variation in density significantly. In addition, NRLMSISE00, DTM-2000 and JB2006 model predictions are compared with density measurements from CHAMP to assess their accuracy, and the results show that these models underestimate the response of the thermosphere to variations induced by solar rotation. Next, the equatorial topside ionospheric plasmas densities Ni are analyzed for solar irradiance variability effects during the period 2002-2005. Linear correlations between Ni and the solar irradiances for different wavelength ranges reveal significantly different characteristics. XUV (0-35 nm) and EUV (115-130 nm) show higher correlation with Ni for the long-term variations, whereas EUV (35-115 nm) show higher correlation for the short-term variations. Moreover, partial correlation analysis shows that the long-term variations of Ni are affected by both XUV (0-35 nm) and EUV (35-115 nm), whereas XUV (0-35 nm) play a more important role; the short-term variations of Ni are mostly affected by EUV (35-115 nm). Furthermore, a pronounced period of about 27 days is present in both Ni and solar irradiance data of 2003 and 2004, and a pronounced period of about 54 days is also revealed in 2004. Finally, prompted by previous studies that have suggested solar EUV radiation as a means of driving the semiannual variation, we investigate the intra-annual variation in thermosphere neutral density near 400 km during 2002-2005. The intra-annual variation, commonly referred to as the ‘semiannual variation’, is characterized by significant latitude structure, hemispheric asymmetries, and inter-annual variability. The magnitude of the maximum yearly difference, from the yearly minimum to the yearly maximum, varies by as much as 60% from year to year, and the phases of the minima and maxima also change by 20-40 days from year to year. Each annual harmonic of the intra-annual variation, namely, annual, semiannual, ter-annual and quatra-annual, exhibits a decreasing trend from 2002 through 2005 that is correlated with the decline in solar activity. In addition, some variations in these harmonics are correlated with geomagnetic activity, as represented by the daily mean value of Kp. Recent empirical models of the thermosphere are found to be deficient in capturing most of the latitude dependencies discovered in our data. In addition, the solar flux and geomagnetic activity proxies that we have employed do not capture some latitude and inter-annual variations detected in our data. It is possible that these variations are partly due to other effects, such as seasonal-latitudinal variations in turbopause altitude (and hence O/N2 composition) and ionosphere coupling processes that remain to be discovered in the context of influencing the intra-annual variations depicted here. Our results provide a new dataset to challenge and validate thermosphere-ionosphere general circulation models that seek to delineate the thermosphere intra-annual variation and to understand the various competing mechanisms that may contribute to its existence and variability. We furthermore suggest that the term “intra-annual” variation be adopted to describe the variability in thermosphere and ionosphere parameters that is well-captured through a superposition of annual, semiannual, ter-annual, and quatra-annual harmonic terms, and that “semiannual’ be used strictly in reference to a pure 6-monthly sinusoidal variation. Moreover, we propose the term “intra-seasonal” to refer to those shorter-term variations that arise as residuals from the above Fourier representation.

Roles of equatorial waves and western boundary reflection in the seasonal circulation of the equatorial Indian Ocean

An ocean general circulation model (OGCM) is used to study the roles of equatorial waves and western boundary reflection in the seasonal circulation of the equatorial Indian Ocean. The western boundary reflection is defined as the total Kelvin waves leaving the western boundary, which include the reflection of the equatorial Rossby waves as well as the effects of alongshore winds, off-equatorial Rossby waves, and nonlinear processes near the western boundary. The evaluation of the reflection is based on a wave decomposition of the OGCM results and experiments with linear models. It is found that the alongshore winds along the east coast of Africa and the Rossby waves in the off-equatorial areas contribute significantly to the annual harmonics of the equatorial Kelvin waves at the western boundary. The semiannual harmonics of the Kelvin waves, on the other hand, originate primarily from a linear reflection of the equatorial Rossby waves. The dynamics of a dominant annual oscillation of sea level coexisting with the dominant semiannual oscillations of surface zonal currents in the central equatorial Indian Ocean are investigated. These sea level and zonal current patterns are found to be closely related to the linear reflections of the semiannual harmonics at the meridional boundaries. Because of the reflections, the second baroclinic mode resonates with the semiannual wind forcing; that is, the semiannual zonal currents carried by the reflected waves enhance the wind-forced currents at the central basin. Because of the different behavior of the zonal current and sea level during the reflections, the semiannual sea levels of the directly forced and reflected waves cancel each other significantly at the central basin. In the meantime, the annual harmonic of the sea level remains large, producing a dominant annual oscillation of sea level in the central equatorial Indian Ocean. The linear reflection causes the semiannual harmonics of the incoming and reflected sea levels to enhance each other at the meridional boundaries. In addition, the weak annual harmonics of sea level in the western basin, resulting from a combined effect of the western boundary reflection and the equatorial zonal wind forcing, facilitate the dominance by the semiannual harmonics near the western boundary despite the strong local wind forcing at the annual period. The Rossby waves are found to have a much larger contribution to the observed equatorial semiannual oscillations of surface zonal currents than the Kelvin waves. The westward progressive reversal of seasonal surface zonal currents along the equator in the observations is primarily due to the Rossby wave propagation.

Mean sea surface heights of the South and East China Seas from ocean circulation model and geodetic leveling

The mean sea surface heights (sea surface topography) of the South China, East China, Yellow and Bohai Seas are derived from an ocean general circulation model and surface air pressure. The circulation model covers the global oceans, with fine grid (1/6degrees) covering the East Asian marginal seas and coarse grid (31) covering the rest part of the global oceans. The result shows that the China 1985 National Altitude Datum is 24.7 cm above the me-an sea surface height of the world oceans. The mean sea surface in the coastal ocean adjacent to China is higher in the south than in the north. Intercomparison of the model results with the geodetic leveling measurements at 28 coastal tidal stations shows a standard deviation of 4.8 cm and a fitting coefficient of 95.3%. After correction through linear regression, the standard deviation is reduced to 4.5 cm. This indicates that the accuracy of model results is sufficient for practical application. Based on the model results, the mean sea surface heights for the study area with a resolution of 1/6 degree are given. This result also links the mean sea levels at islands with those on the mainland coast and gives the mean sea surface heights at tidal stations in the Taiwan Island, the Dongsha Islands, the Yisha Islands and the Nansha Islands relative to the China 1985 National Altitude Datum.

COASTAL SEA-LEVEL AND THE LARGE-SCALE CLIMATE STATE - A DOWNSCALING EXERCISE FOR THE JAPANESE ISLANDS

A major problem which is envisaged in the course of man-made climate change is sea-level rise. The global aspect of the thermal expansion of the sea water likely is reasonably well simulated by present day climate models; the variation of sea level, due to variations of the regional atmospheric forcing and of the large-scale oceanic circulation, is not adequately simulated by a global climate model because of insufficient spatial resolution. A method to infer the coastal aspects of sea level change is to use a statistical ''downscaling'' strategy: a linear statistical model is built upon a multi-year data set of local sea level data and of large-scale oceanic and/or atmospheric data such as sea-surface temperature or sea-level air-pressure. We apply this idea to sea level along the Japanese coast. The sea level is related to regional and North Pacific sea-surface temperature and sea-level air pressure. Two relevant processes are identified. One process is the local wind set-up of water due to regional low-frequency wind anomalies; the other is a planetary scale atmosphere-ocean interaction which takes place in the eastern North Pacific.

土壤侵蚀高维指数非线性模型及其参数辨识

模型验证和数据库组建基础上,用WinEPIC模型定量模拟研究了黄土高原半湿润区长武、半干旱区固原和半干旱偏旱区海原20～30年内苜蓿草地水分生产潜力、10m土层土壤有效含水量和土壤湿度剖面分布特征的动态变化.结果表明:长武、固原和海原苜蓿草地水分生产潜力模拟值随降水量变化而呈现波动性降低趋势,其平均值分别为8.81、3.83和2.48t.hm-2;长武、固原和海原苜蓿草地10m土层逐月土壤有效含水量模拟值均呈现明显的波动性降低趋势,模拟初期,4～8年生苜蓿草地土壤干燥化趋势十分强烈,此后,随降水量变化长期在较低水平上波动;随着苜蓿生长年限的延长,苜蓿草地土壤干层逐年加深、加厚,长武、固原和海原土壤干层分布深度达到10m所需时间依次为6、6和4年,此后苜蓿草地降水渗深以下土层长期维持较为稳定的干燥化状态;苜蓿草地水分持续利用的合理年限为半湿润区8～10年,半干旱区6～8年,半干旱偏旱区4～6年.

Dynamics of intraseasonal sea level and thermocline variability in the equatorial Atlantic during 2002-03

Satellite and in situ observations in the equatorial Atlantic Ocean during 2002-03 show dominant spectral peaks at 40-60 days and secondary peaks at 10-40 days in sea level and thermocline within the intraseasonal period band (10-80 days). A detailed investigation of the dynamics of the intraseasonal variations is carried out using an ocean general circulation model, namely, the Hybrid Coordinate Ocean Model (HYCOM). Two parallel experiments are performed in the tropical Atlantic Ocean basin for the period 2000-03: one is forced by daily scatterometer winds from the Quick Scatterometer (QuikSCAT) satellite together with other forcing fields, and the other is forced by the low-passed 80-day version of the above fields. To help in understanding the role played by the wind-driven equatorial waves, a linear continuously stratified ocean model is also used. Within 3 degrees S-3 degrees N of the equatorial region, the strong 40-60-day sea surface height anomaly (SSHA) and thermocline variability result mainly from the first and second baroclinic modes equatorial Kelvin waves that are forced by intraseasonal zonal winds, with the second baroclinic mode playing a more important role. Sharp 40-50-day peaks of zonal and meridional winds appear in both the QuikSCAT and Pilot Research Moored Array in the Tropical Atlantic (PIRATA) data for the period 2002-03, and they are especially strong in 2002. Zonal wind anomaly in the central-western equatorial basin for the period 2000-06 is significantly correlated with SSHA across the equatorial basin, with simultaneous/ lag correlation ranging from-0.62 to 0.74 above 95% significance. Away from the equator (3 degrees-5 degrees N), however, sea level and thermocline variations in the 40-60-day band are caused largely by tropical instability waves (TIWs). On 10-40-day time scales and west of 10 degrees W, the spectral power of sea level and thermocline appears to be dominated by TIWs within 5 degrees S-5 degrees N of the equatorial region. The wind-driven circulation, however, also provides a significant contribution. Interestingly, east of 10 W, SSHA and thermocline variations at 10 40- day periods result almost entirely from wind-driven equatorial waves. During the boreal spring of 2002 when TIWs are weak, Kelvin waves dominate the SSHA across the equatorial basin (2 degrees S-2 degrees N). The observed quasi-biweekly Yanai waves are excited mainly by the quasi-biweekly meridional winds, and they contribute significantly to the SSHA and thermocline variations in 1 degrees-5 degrees N and 1 degrees-5 degrees S regions.

Origin and Pathway of Equatorial 13 degrees C Water in the Pacific Identified by a Simulated Passive Tracer and Its Adjoint

The origin and pathway of the thermostad water in the eastern equatorial Pacific Ocean, often referred to as the equatorial 13 degrees C Water, are investigated using a simulated passive tracer and its adjoint, based on circulation estimates of a global general circulation model. Results demonstrate that the source region of the 13 degrees C Water lies well outside the tropics. In the South Pacific, some 13 degrees C Water is formed northeast of New Zealand, confirming an earlier hypothesis on the water's origin. The South Pacific origin of the 13 degrees C Water is also related to the formation of the Eastern Subtropical Mode Water (ESTMW) and the Sub-Antarctic Mode Water (SAMW). The portion of the ESTMW and SAMW that eventually enters the density range of the 13 degrees C Water (25.8 < sigma(theta) < 26.6 kg m(-3)) does so largely by mixing. Water formed in the subtropics enters the equatorial region predominantly through the western boundary, while its interior transport is relatively small. The fresher North Pacific ESTMW and Central Mode Water (CMW) are also important sources of the 13 degrees C Water. The ratio of the southern versus the northern origins of the water mass is about 2 to 1 and tends to increase with time elapsed from its origin. Of the total volume of initially tracer-tagged water in the eastern equatorial Pacific, approximately 47.5% originates from depths above sigma(theta) = 25.8 kg m(-3) and 34.6% from depths below sigma(theta) = 26.6 kg m(-3), indicative of a dramatic impact of mixing on the route of subtropical water to becoming the 13 degrees C Water. Still only a small portion of the water formed in the subtropics reaches the equatorial region, because most of the water is trapped and recirculates in the subtropical gyre.

Long-Wave Dynamics of Sea Level Variations during Indian Ocean Dipole Events

Long-wave dynamics of the interannual variations of the equatorial Indian Ocean circulation are studied using an ocean general circulation model forced by the assimilated surface winds and heat flux of the European Centre for Medium-Range Weather Forecasts. The simulation has reproduced the sea level anomalies of the Ocean Topography Experiment (TOPEX)/Poseidon altimeter observations well. The equatorial Kelvin and Rossby waves decomposed from the model simulation show that western boundary reflections provide important negative feedbacks to the evolution of the upwelling currents off the Java coast during Indian Ocean dipole (IOD) events. Two downwelling Kelvin wave pulses are generated at the western boundary during IOD events: the first is reflected from the equatorial Rossby waves and the second from the off-equatorial Rossby waves in the southern Indian Ocean. The upwelling in the eastern basin during the 1997-98 IOD event is weakened by the first Kelvin wave pulse and terminated by the second. In comparison, the upwelling during the 1994 IOD event is terminated by the first Kelvin wave pulse because the southeasterly winds off the Java coast are weak at the end of 1994. The atmospheric intraseasonal forcing, which plays an important role in inducing Java upwelling during the early stage of an IOD event, is found to play a minor role in terminating the upwelling off the Java coast because the intraseasonal winds are either weak or absent during the IOD mature phase. The equatorial wave analyses suggest that the upwelling off the Java coast during IOD events is terminated primarily by western boundary reflections.

Algorithm Studies on How to Obtain a Conditional Nonlinear Optimal Perturbation (CNOP)

The conditional nonlinear optimal perturbation (CNOP), which is a nonlinear generalization of the linear singular vector (LSV), is applied in important problems of atmospheric and oceanic sciences, including ENSO predictability, targeted observations, and ensemble forecast. In this study, we investigate the computational cost of obtaining the CNOP by several methods. Differences and similarities, in terms of the computational error and cost in obtaining the CNOP, are compared among the sequential quadratic programming (SQP) algorithm, the limited memory Broyden-Fletcher-Goldfarb-Shanno (L-BFGS) algorithm, and the spectral projected gradients (SPG2) algorithm. A theoretical grassland ecosystem model and the classical Lorenz model are used as examples. Numerical results demonstrate that the computational error is acceptable with all three algorithms. The computational cost to obtain the CNOP is reduced by using the SQP algorithm. The experimental results also reveal that the L-BFGS algorithm is the most effective algorithm among the three optimization algorithms for obtaining the CNOP. The numerical results suggest a new approach and algorithm for obtaining the CNOP for a large-scale optimization problem.

IMPACTS OF THE TROPICAL PACIFIC COUPLED PROCESS ON THE INTERANNUAL VARIABILITY IN THE INDIAN OCEAN

The basic features of climatology and interannual variations of tropical Pacific and Indian Oceans were analyzed using a coupled general circulation model (CGCM), which was constituted with an intermediate 2.5-layer ocean model and atmosphere model ECHAM4. The CGCM well captures the spatial and temporal structure of the Pacific El Nino-Southern Oscillation (ENSO) and the variability features in the tropical Indian Ocean. The influence of Pacific air-sea coupled process on the Indian Ocean variability was investigated carefully by conducting numerical experiments. Results show that the occurrence frequency of positive/negative Indian Ocean Dipole (IOD) event will decrease/increase with the presence/absence of the coupled process in the Pacific Ocean. Further analysis demonstrated that the air-sea coupled process in the Pacific Ocean affects the IOD variability mainly by influencing the zonal gradient of thermocline via modulating the background sea surface wind.

Discussion on the basic mathematical models of neurons - Double synaptic weight neuron in general-purpose neurocomputer, theory and application

Based on the introduction of the traditional mathematical models of neurons in general-purpose neurocomputer, a novel all-purpose mathematical model-Double synaptic weight neuron (DSWN) is presented, which can simulate all kinds of neuron architectures, including Radial-Basis-Function (RBF) and Back-propagation (BP) models, etc. At the same time, this new model is realized using hardware and implemented in the new CASSANN-II neurocomputer that can be used to form various types of neural networks with multiple mathematical models of neurons. In this paper, the flexibility of the new model has also been described in constructing neural networks and based on the theory of Biomimetic pattern recognition (BPR) and high-dimensional space covering, a recognition system of omni directionally oriented rigid objects on the horizontal surface and a face recognition system had been implemented on CASSANN-H neurocomputer. The result showed DSWN neural network has great potential in pattern recognition.