Trends of sea level variations in the South China Sea from merged altimetry data

In this study, more than 13 yr of merged altimetry sea level anomalies (SLA) data were used to analyze the trends of sea level variations in the South China Sea (SCS). The result shows that the mean sea level over the SCS has a rise rate of 11.3 mm/yr during 1993–2000 and a fall rate of 11.8 mm/yr during 2001–2005. The geographical distribution of the sea level variations over the SCS is asymmetric with a pronounced variation existing in the deep water. The trends of thermosteric sea level variations were also examined using Ishii data and MITgcm assimilation data. The result indicates that the thermal change of the upper layer of the SCS has a significant contribution to the sea level variations. Heat budget analysis suggests that heat advection may be a key factor influencing the thermal change. Apart from thermal contribution, the effect of water exchange on the sea level variations was also studied.     

Reconstruction of past decades sea level using thermosteric sea level, tide gauge, satellite altimetry and ocean reanalysis data

This study investigates past sea level reconstruction (over 1950–2003) based on tide gauge records and EOF spatial patterns from different 2-D fields. In a first step, we test the influence on the reconstructed signal of the 2-D fields temporal coverage. For that purpose we use global grids of thermosteric sea level data, available over 1950–2003. Different time spans (in the range 10–50 yr) for the EOF spatial patterns, and different geographical distributions for the 1-D thermosteric sea level time series (interpolated at specific locations from the 2-D grids), are successively used to reconstruct the 54-year long thermosteric sea level signal. In each case we compare the reconstructed trend map with the reference. The simulation indicates that the longer the time span covered by the spatial EOFs, the closer to the reference the reconstructed thermosteric sea level trends. In a second step, we apply the method to reconstructing 2-D sea level data over 1950–2003, combining sparse tide gauge records available since 1950, with EOF spatial patterns from different sources: (1) thermosteric sea level grids over 1955–2003, (2) sea level grids from Topex/Poseidon satellite altimetry over 1993–2003, and (3) dynamic height grids from the SODA reanalysis over 1958–2001. The reconstructed global mean sea level trend based on thermosteric EOFs (case 1) is significantly lower than the observed trend, while the interannual/decadal sea level fluctuations are well reproduced. Case 2 (Topex/Poseidon EOFs over 1993–2003) leads to a global mean sea level trend over the 54-year time interval very close to the observed trend. But the spatial trends of the reconstruction over 1950–2003 are significantly different from those obtained with thermosteric EOFs. Case 3 (SODA EOFs over 1958–2001) provides a reconstruction trend map over 1950–2003 that differs significantly from the previous two cases. We discuss possible causes for such differences. For the three cases, on the other hand, reconstructed spatial trends over 1993–2003 agree well with the regional sea level trends observed by Topex/Poseidon.     

Steric and mass-induced Mediterranean sea level trends from 14 years of altimetry data

Long-term series of almost 14 years of altimetry data (1992–2005) have been analysed along with Sea Surface Temperature (SST) and temperature and salinity profiles to investigate sea level trends over the Mediterranean Sea. Although sea level variations are mainly driven by the steric contribution, the mass-induced component plays some role in modulating its oscillation. A spatially averaged positive trend of 2.1 ± 0.6 mm/year has been observed, but a change in sign in 2001 seems to appear. Steric effects (mainly on thermal origin) account for  55% of sea level trend. Although Mediterranean Sea is a semi-enclosed basin, this value is comparable to that reported for the global ocean. Sea level rise is particularly important in the Levantine basin south of Crete with values up to 10 ± 1 mm/year. Other areas of sea level rise are localised throughout the Levantine basin and in the Adriatic and Alboran Seas, with more moderate values. Sea level drop areas are localised in the Algerian basin, between the Balearic Islands and the African coasts and, particularly, in the Ionian basin. In this area, negative trends as high as − 10 ± 0.8 mm/year are detected mainly due to the mass-induced contribution, which suggests decadal changes of surface circulation. The inferred sea level trends have been correlated with North Atlantic Oscillation (NAO) indices and a low but significant correlation has been detected between sea level in the Levantine and Balearic basins and NAO index.     

Verification of a coupled climate-hydrological model against Holocene palaeohydrological records

We have coupled a climate model (ECBilt-CLIO-VECODE) and a hydrological model (STREAM) offline to simulate palaeodischarge of nineteen rivers (Amazon, Congo, Danube, Ganges, Krishna, Lena, Mackenzie, Mekong, Meuse, Mississippi, Murray–Darling, Nile, Oder, Rhine, Sacramento–San Joaquin, Syr Darya, Volga, Volta, Zambezi) for three time-slices: Early Holocene (9000–8650 BP), Mid-Holocene (6200–5850 BP) and Recent (1750–2000 AD). To evaluate the model's skill in retrodicting broad changes in mean palaeodischarge we have compared the model results with palaeodischarge estimates from multi-proxy records. We have compared the general trends inferred from the proxy data with statistical differences in modelled discharge between the three periods, thereby developing a technique to assess the level of agreement between the model and proxy data. The quality of the proxy data for each basin has been classed as good, reasonable or low. Of the model runs for which the proxy data were good or reasonable, 72% were in good agreement with the proxy data, and 92% were in at least reasonable agreement. We conclude that the coupled climate-hydrological model performs well in simulating mean discharge in the time-slices studied. The discharge trends inferred from the proxy and model data closely follow latitudinal and seasonal variations in insolation over the Holocene. For a number of basins for which agreement was not good we have identified specific mechanisms which could be responsible for the discrepancy, primarily the absence of the Laurentide ice sheet in our model. In order to use the model in an operational sense within water management studies it would be useful to use a higher spatial resolution and a daily time-step.     

Sea level budget over 2003–2008: A reevaluation from GRACE space gravimetry,satellite altimetry and Argo

From the IPCC 4th Assessment Report published in 2007, ocean thermal expansion contributed by ～ 50% to the 3.1 mm/yr observed global mean sea level rise during the 1993–2003 decade, the remaining rate of rise being essentially explained by shrinking of land ice. Recently published results suggest that since about 2003, ocean thermal expansion change, based on the newly deployed Argo system, is showing a plateau while sea level is still rising, although at a reduced rate (～ 2.5 mm/yr). Using space gravimetry observations from GRACE, we show that recent years sea level rise can be mostly explained by an increase of the mass of the oceans. Estimating GRACE-based ice sheet mass balance and using published estimates for glaciers melting, we further show that ocean mass increase since 2003 results by about half from an enhanced contribution of the polar ice sheets – compared to the previous decade – and half from mountain glaciers melting. Taking also into account the small GRACE-based contribution from continental waters (< 0.2 mm/yr), we find a total ocean mass contribution of ～ 2 mm/yr over 2003–2008. Such a value represents ～ 80% of the altimetry-based rate of sea level rise over that period. We next estimate the steric sea level (i.e., ocean thermal expansion plus salinity effects) contribution from: (1) the difference between altimetry-based sea level and ocean mass change and (2) Argo data. Inferred steric sea level rate from (1) (～ 0.3 mm/yr over 2003–2008) agrees well with the Argo-based value also estimated here (0.37 mm/yr over 2004–2008). Furthermore, the sea level budget approach presented in this study allows us to constrain independent estimates of the Glacial Isostatic Adjustment (GIA) correction applied to GRACE-based ocean and ice sheet mass changes, as well as of glaciers melting. Values for the GIA correction and glacier contribution needed to close the sea level budget and explain GRACE-based mass estimates over the recent years agree well with totally independent determinations.     

Thermosteric sea level rise for the past 50 years; comparison with tide gauges and inference on water mass contribution

In this paper we compare sea level trends observed at a few selected tide gauges of good quality records with thermosteric (i.e., due to ocean temperature change) sea level trends over 1950–1998 using different gridded ocean temperature data sets from Levitus et al. (2000) [Levitus, S., Stephens, C., Antonov, J.I., Boyer, T.P., 2000. Yearly and Year-Season Upper Ocean Temperature Anomaly Fields, 1948–1998. U.S. Gov. Printing Office, Washington, D.C. pp. 23.], Ishii et al. (2003) [Ishii, M., Kimoto, M., Kachi, M., 2003. Historical ocean subsurface temperature analysis with error estimates, Mon. Weather Rev., 131, 51–73.] and Levitus et al. (2005) [Levitus S., Antonov, J.I., Boyer, T.P., 2005. Warming of the world ocean, 1955–2003. Geophys. Res. Lett. 32, L02604. doi:10.1029/2004GL021592.]. When using the Levitus data, we observe very high thermosteric rates at sites located along the northeast coast of the US, north of 37°N. Such high rates are not observed with the Ishii data. Elsewhere, thermosteric rates agree reasonably well whatever the data set. Excluding the northeast US coastline sites north of 37°N, we compare tide gauge-based sea level trends with thermosteric trends and note that, in spite of a significant correlation, the latter are too small to explain the observed trends. After correcting for thermosteric sea level trends, residual (observed minus thermosteric) trends have an average value of 1.4 ± 0.5 mm/year, which should have an eustatic (i.e., due to ocean mass change) origin. This result supports the recent investigation by Miller and Douglas (2004) [Miller, L., Douglas, B.C., 2004. Mass and volume contributions to 20th century global sea level rise. Nature 428, 406–408.] which suggests that a dominant eustatic contribution is needed to explain the rate of sea level rise of the last decades observed by tide gauges, and shows that Cabanes et al. (2001) [Cabanes, C., Cazenave, A., Le Provost, C., 2001. Sea level rise during past 40 years determined from satellite and in situ observations. Science 294, 840–842.] arrived at an incorrect conclusion due to peculiarities in the gridded Levitus et al. (2000) [Levitus, S., Stephens, C., Antonov, J.I., and Boyer, T.P., 2000. Yearly and Year-Season Upper Ocean Temperature Anomaly Fields, 1948–1998. U.S. Gov. Printing Office, Washington, D.C. pp. 23.] data set.     

Effects of anthropogenic intervention in the land hydrologic cycle on global sea level rise

Recent studies suggest that anthropogenic modification of land hydrology (e.g. through groundwater mining, dam building, irrigation, deforestation, wetlands drainage, and urbanization) could significantly impact sea-level rise, although the magnitude and sign of this effect have been widely debated. This paper attempts a comprehensive overview of the effects of human activities on land hydrology. Estimates are provided for the volumes of water associated with each of the major anthropogenic processes and the corresponding equivalent in sea level.Groundwater mining; and runoff from paved and built-up areas are two major sources of water added to the ocean. In contrast, storage of water behind dams, losses through percolation, and evapotranspiration from irrigated fields withhold water that would otherwise flow to the sea. The net effect of these processes holds back the equivalent of 0.8 +- 0.4 mm/yr from sea-level rise. This is a magnitude comparable to, but in the opposite direction from the currently observed sea-level rise of 1–2 mm/yr. These estimates are still preliminary, awaiting better documentation. Coupling of improved land hydrology models with GCMs will help in analysis of feedbacks, especially the partitioning of water among runoff, infiltration, and evaporation.     

Geocentric sea-level trend estimates from GPS analyses at relevant tide gauges world-wide

The problem of correcting the tide gauge records for the vertical land motion upon which the gauges are settled has only been partially solved. At best, the analyses so far have included model corrections for one of the many processes that can affect the land stability, namely the Glacial-Isostatic Adjustment (GIA). An alternative approach is to measure (rather than to model) the rates of vertical land motion at the tide gauges by means of space geodesy. A dedicated GPS processing strategy is implemented to correct the tide gauges records, and thus to obtain a GPS-corrected set of ‘absolute’ or geocentric sea-level trends. The results show a reduced dispersion of the estimated sea-level trends after application of the GPS corrections. They reveal that the reference frame implementation is now achieved within the millimetre accuracy on a weekly basis. Regardless of the application, whether local or global, we have shown that GPS data analysis has reached the maturity to provide useful information to separate land motion from oceanic processes recorded by the tide gauges or to correct these latter. For comparison purposes, we computed the global average of sea-level change according to Douglas [Douglas, B.C., 2001. Sea level change in the era of the recording tide gauge. Int. Geophys. Ser., 75, pp. 37–64.] rules, whose estimate is 1.84 ± 0.35 mm/yr after correction for the GIA effect [Peltier, W.R., 2001. Global glacial isostatic adjustment and modern instrumental records of relative sea level history. Int. Geophys. Ser., 75, pp. 65–95.]. We obtain a value of 1.31 ± 0.30 mm/yr, a value which appears to resolve the ‘sea level enigma’ [Munk, W., 2002. Twentieth century sea level: an enigma. Proc. Natl. Acad. Sci. U.S.A., 99(10), pp. 6550–6555].     

Sea level forcing in the Mediterranean Sea between 1960 and 2000

Sea level trends and inter-annual variability in the Mediterranean Sea for the period 1960–2000 is explored by comparing observations from tide gauges with sea level hindcasts from a barotropic 2D circulation model, and two full primitive equation 3D ocean circulation models, a regional one and the Mediterranean component of a global one,. In the 2D model, 50% of the sea level variance was found to result from the wind and atmospheric pressure forcing. In the 3D models, 20% of the sea level variance was explained by the steric effects. The sea level residuals at the tide gauges locations, calculated by subtraction of the 2D model output from the sea level observations are significantly correlated (r = 0.4) with the steric signals from the 3D models. After the removal of the atmospheric and the steric contributions the tide-gauge sea level records indicate a period where sea level was stable (1960–1975) and a period where sea level was rising (1975–2000) with rates in the range 1.1–1.8 mm/yr. A part of the residual trend can be explained by the contribution of local land movements (0.3 mm/yr) while its major part indicates a global signal, probably mass addition, appearing after 1975.     

The glacial–interglacial paradox: the relation between the global means of sea level and sea surface temperature

Sea level variability during the Quarternary is simulated using a stochastic climate model, and a sensitivity relation for the change in net oceanic evaporation due to a change in sea surface temperature. In the application of this relation, it is assumed that the greater part of the change in net oceanic evaporation causes changes in the land ice storage, rather than being directly returned to the ocean by rivers. The analysis suggests that the observed sea level changes can be interpreted as due to the transfer of heat to the deep ocean from the surface mixed layer, arising from random radiation perturbations of the same variance as would give rise to the interannual variability of the global temperature series. The paradox is that glacial conditions (increase in ice storage) are favoured by positive (temperate) sea surface temperature anomalies, and interglacial conditions (decrease in ice storage) by negative (temperate) sea surface temperature anomalies. The evolution of both these regimes, which are inherently unstable, appears to be controlled by the deep water formation process, while albedo feedback is of minor importance. Fluvial feedback, (in which as the ice storage increases the fluvial inflow decreases), however, is found to be an important process, and a small sensitivity of river inflow to storage is consistent with forcing by random variability or by astronomical forcing. A simple analytical model incorporating the key processes of oceanic evaporation and fluvial feedback is presented. The analysis points to the importance of an accurate river model for climate system modelling.     

20 years of mass balances on the Piloto glacier, Las Cuevas river basin, Mendoza, Argentina

Climatic changes of the 20th century have altered the water cycle in the Andean basins of central Argentina. The most visible change is seen in the mountain glaciers, with loss of part of their mass due to decreasing thickness and a substantial recession in the last 100 years. This paper briefly describes the results of glacier mass balance research since 1979 in the Piloto Glacier at the Cajón del Rubio, in the headwaters of Las Cuevas River, presenting new results for the period 1997–2003. Very large interannual variability of net annual specific balance is evident, due largely to variations in winter snow accumulation, with a maximum net annual value of + 151 cm w.e. and a minimum value of - 230 cm w.e. Wet El Niño years are normally associated with positive net annual balances, while dry La Niña years generally result in negative balances. Within the 24-year period, 67% of the years show negative net annual specific balances, with a cumulative mass balance loss of - 10.50 m water equivalent (w.e.). Except for exceptions normally related to El Niño events, a general decreasing trend of winter snow accumulation is evident in the record, particularly after 1992, which has a strong effect in the overall negative mass balance values. The glacier contribution to Las Cuevas River runoff is analysed based on the Punta de Vacas River gauge station for a hypothetical year without snow precipitation (YWSP), when the snowmelt component is zero. Extremely dry years similar to a YWSP have occurred in 1968–1969, 1969–1970 and 1996–1997. The Punta de Vacas gauge station is located 62 km downstream from Piloto Glacier, and the basin contains 3.0% of uncovered glacier ice and 3.7% of debris-covered ice. The total glacier contribution to Las Cuevas River discharge is calculated as 82 ± 8% during extremely dry years. If glacier wastage continues at the present trend as observed during the last 2 decades, it will severely affect the water resources in the arid central Andes of Argentina.     

GRACE重力计划在揭示地球系统质量重新分布中的应用

2002年3月成功发射的美德合作卫星重力计划GRACE(Gravity Recovery And ClimateExperiment),即将提供空间分辨率约为200 km而时间分辨率为1个月的时变地球重力场模型序列。GRACE计划的星座由两颗相距约220 km,高度保持为300-500 km、倾角保持约90°的近极轨卫星组成。由于采用星载GPS和非保守力加速度计等高精度定轨技术,以及高精度的星一星跟踪数据反演地球重力场,在几百公里和更大空间尺度上, GRACE重力场的精度大大超过此前的卫星重力计划。根据GRACE时变重力场反演的地球系统质量重新分布,将对固体地球物理、海洋物理、气候学以及大地测量等应用有重要的意义。虽然其设计寿命只有5 yr,但研究表明GRACE的结果可用于研究北极冰长期时间尺度的变化,并进而研究极冰融化对全球气候变化,特别是对海平面长期变化的影响。在季节性时间尺度上,利用GRACE重力场反演的质量重新分布足以揭示平均小于1 cm的地表水变化,或小于1 mbar。的海底压强变化。除了巨大的社会效益和经济效益外,这些变化对了解地球系统的物质循环(主要是水循环)和能量循环有非常重要的意义。介绍GRACE重力场揭示的地球系统质量重新分布,为理解其地球物理应用提供必需的准备;同时针对我国大陆和沿海地区的地球物理应用提出初步的设想。     

Climate variability, vegetation productivity and people at risk

Human domination of ecosystems has been pervasive over the last century, with nearly half of Earth's surface transformed by human actions. It is widely accepted that humans appropriate up to 50% of global net primary production (NPP), the energy base of all the trophic levels on the land surface. Yet, despite the important role of vegetation productivity for defining Earth habitability, the covariation of NPP and human population distribution has not been analyzed in depth. We used recently available satellite-based NPP estimates, along with gridded population at 0.5° resolution, first, to identify the global distribution of human population with reference to average NPP and to the various climatic constraints (temperature, water and cloud cover) that limit NPP, second, to analyze recent trends in global NPP in relation to population trends, and third, to identify populations that are vulnerable to changes in NPP due to interannual variability in climate. Our results indicate that over half of the global human population is presently living in areas with above the average NPP of 490 g C m⁻² year⁻¹. By 1998, nearly 56% of global population lived in regions where water availability strongly influences NPP. Per capita NPP declined over much of Africa between 1982 and 1998, in spite of the estimated increases in NPP over the same period. On average, NPP over 40% of the total vegetated land surface has shown significant correlations with ENSO-induced climate variability affecting over 2.8 billion people.     

Modeling the interannual variability and trends in gross and net primary productivity of tropical forests from 1982 to 1999

The role of tropical ecosystems in global carbon cycling is uncertain, at least partially due to an incomplete understanding of climatic forcings of carbon fluxes. To reduce this uncertainty, we simulated and analyzed 1982–1999 Amazonian, African, and Asian carbon fluxes using the Biome-BGC prognostic carbon cycle model driven by National Centers for Environmental Prediction reanalysis daily climate data. We first characterized the individual contribution of temperature, precipitation, radiation, and vapor pressure deficit to interannual variations in carbon fluxes and then calculated trends in gross primary productivity (GPP) and net primary productivity (NPP). In tropical ecosystems, variations in solar radiation and, to a lesser extent, temperature and precipitation, explained most interannual variation in GPP. On the other hand, temperature followed by solar radiation primarily determined variation in NPP. Tropical GPP gradually increased in response to increasing atmospheric CO₂. Confirming earlier studies, changes in solar radiation played a dominant role in CO₂ uptake over the Amazon relative to other tropical regions. Model results showed negligible impacts from variations and trends in precipitation or vapor pressure deficits on CO₂ uptake.     

Paleolakes, paleofloods, and depressions in Aurorae and Ophir Plana, Mars: Connectivity of surface and subsurface hydrological systems

The plains of Aurorae and Ophir in the equatorial region of Mars display geomorphic evidence indicative of extensive but generally short-lived paleohydrological processes. Elaver Vallis in Aurorae Planum south of Ganges Chasma is an outflow channel system >180 km long, and here inferred to have formed by cataclysmic spillover flooding from a paleolake(s) contained in the Morella crater basin. Ganges Cavus is an enormous 5-km-deep depression of probable collapse origin located in the Morella basin. The fluid responsible for the infilling of the Morella basin likely emerged at least partially through Ganges Cavus or its incipient depression, and it may have been supplied also from small-scale springs in the basin. Similar paleohydrological processes are inferred also in Ophir Planum. It is reasonable to assume that water, sometimes sediment-laden and/or mixed with gases, was the responsible fluid for these phenomena although some of the observed features could be explained by non-aqueous processes such as volcanism. Water emergence may have occurred as consequences of ground ice melting or breaching of cryosphere to release water from the underlying hydrosphere. Dike intrusion is considered to be an important cause of formation for the cavi and smaller depressions in Aurorae and Ophir Plana, explaining also melting of ground ice or breaching of cryosphere. Alternatively, the depressions and crater basins may have been filled by regional groundwater table rising during the period(s) when cryosphere was absent or considerably thin. The large quantities of water necessary for explaining the paleohydrological processes in Aurorae and Ophir Plana could have been derived through crustal migration from the crust of higher plains in western Ophir Planum where water existed in confined aquifers or was produced by melting of ground ice due to magmatic heating or climatic shift, or from a paleolake in Candor Chasma further west.     

Control of the Antarctic ice sheet by ocean–ice interaction

The Antarctic ice cap is the largest ice sheet of modern times. It is of considerable importance to predict the sea level variability due to the associated changes in ice volume. We present the results of a simple grounded ice sheet model, developed from Oerlemans [Oerlemans, J., 2002. Global dynamics of the Antarctic Ice Sheet, Climate Dynamics 19, 85–93.], in which the net oceanic evaporation influences the ice cap volume in two ways, through changes in: (i) the accumulation rate, and (ii) the mean sea level. The net evaporation changes are driven by the sea surface temperature (SST) anomaly time series of Howard [Howard, W.R., 1997. A warm future in the past, Nature, 388, 418–419.] for the subantarctic Southern Ocean over the period 220 kyr to the present. The effect of the waxing and waning of the northern hemisphere ice sheets is integrated into the model using an independent model, in which ice melting depends on the SST anomaly and ice calving depends on the sea level anomaly. A series of analytical expressions are derived for the related properties of the coupled ocean–ice system applicable over time scales of 100 kyr, which show, in particular, that the Antarctic ice cap volume changes are due mainly to the effects of the northern hemisphere ice sheets on sea level (which influences ice calving), rather than directly to changes in SST, and hence the ice cap volume is greatest during interglacial periods. This conclusion, which is independent of the specification of the ice melting regime for the northern hemisphere ice sheets, strongly suggests that the changes in accumulation flux estimated from the Vostok proxy temperature data and used in other studies of the Antarctic mass balance have been overestimated. A simple expression is also presented for the lag of ice cap volume to SST, and it is found that the predictions for the mean sea level variability are similar to observations for a melting flux of the northern hemisphere ice sheets about twice their accumulation flux due to the net oceanic evaporation, except during major deglaciations when these two fluxes appear to be of similar magnitude.     

Geo-informational simulation of possible changes in Central Asian water resources

This study simulates water resources in the Tien Shan alpine basins to forecast how global and regional climate changes would affect river runoff. The model employed annual mean values for the major characteristics of the water cycle: annual air temperature, precipitation, evapotranspiration and river runoff. The simulation was based on 304 hydro-meteorological stations, 23 precipitation sites, 328 high altitudinal points with glaciological measurements, 123 stream-gauges, and 54 evaporation sites, and it took into account topography. The findings were simulated over Tien Shan relief using a 1:500,000 scale 100 m grid resolution Digital Elevation Model. An applicable GIS-based distributed River Runoff Model was implemented in regional conditions and tested in the Tien Shan basins. The annual evapotranspiration exceeds the river runoff in the Tien Shan watersheds particularly up to 3700 m. Hypothetical climate-change scenarios in the Tien Shan predict that by 2100 river runoff will increase by 1.047 times with an increase in air temperature averaging 3 °C and an increase in precipitation averaging 1.2 times the current levels. Change in precipitation, rather than temperature, is the main parameter determining river runoff in the Tien Shan. The maximum ratio for predicted river runoff could reach up to 2.2 and the minimum is predicted to be 0.55 times current levels. This possibly dramatic change in river runoff indicates on non-linear system response caused mainly by the non-linear response of evapotranspiration from air temperature and precipitation changes. In the frame of forecasted possible climate change scenarios the probability of river runoff growth amounts 83–87% and probability of this decline is 17–13% by 2100 in the Tien Shan River basins.     

Reconstruction of Mediterranean sea level fields for the period 1945–2000

The distribution of sea level in the Mediterranean Sea is recovered for the period 1945–2000 by using a reduced space optimal interpolation analysis. The method involves estimating empirical orthogonal functions from satellite altimeter data spanning the period 1993–2005 that are then combined with tide gauge data to recover sea level fields over the period 1945–2000. The reconstruction technique is discussed and its robustness is checked through different tests. For the altimetric period (1993–2000) the prediction skill is quantified over the whole domain by comparing the reconstructed fields with satellite altimeter observations. For past times the skill can only be tested locally, by validating the reconstruction against independent tide gauge records. The reconstructed distribution of sea level trends for the period 1945–2000 shows a positive peak in the Ionian Sea (up to 1.5 mm yr^− 1) and a negative peak of − 0.5 mm yr^− 1 in a small area to the south-east of Crete. Positive trends are found nearly everywhere, being larger in the western Mediterranean (between 0.5 and 1 mm yr^− 1) than in the eastern Mediterranean (between 0 and 0.5 mm yr^− 1). The estimated rate of mean sea level rise for the period 1945–2000 is 0.7 ± 0.2 mm yr^− 1, i.e. about a half of the rate estimated for global mean sea level. These overall results do not appear to be very sensitive to the distribution of tide gauges. The poorest results are obtained in open-sea regions with intense mesoscale variability not correlated with any tide gauge station, such as the Algerian Basin.     

Temporal and spatial variability of annual extreme water level in the Pearl River Delta region, China

This paper is concerned with identifying the spatial and temporal patterns in the annual maximum and minimum water level in the Pearl River Delta (PRD) region. The Mann–Kendall test and Pettitt test are used to detect trends and abrupt change points, and the Trend Free Pre-Whitening (TFPW) approach then eliminates the effect of serial correlation in data series with significant autocorrelation. Approximately fifty years of the annual hydrological variables from 18 stations in the three major rivers (the West River, the North River, and the East River) are examined. The changing trends of the extremes in water level show different features in different parts of the PRD region. Generally speaking, in the upper part of the delta, the water levels show a decreasing trend while in the middle and lower part there is an increasing trend. This spatial pattern of the extreme water level variation is unlikely to be due to a long-term change in stream flow in the PRD region because the water level changes do not always coincide with the extreme stream flow variations. Sand excavation initiated in the 1980s and continuing for more than 20 years in almost all tributaries around the PRD region is one of the most serious intensive human activities affecting water levels. The result of the Pettitt test indicates that most abrupt change points occurred in 1980s–1990s, which reveals that sand excavation and channel regulation are likely to have been the most significant factors contributing to the change over this period. These anthropogenic activities modify the annual extreme water level dramatically in a way that affects the morphology of river channels and estuaries of the PRD and also the redistribution of discharge. However, there are differences in the geographic locations of significant trends for the water level investigated, which implies that these impacts are not spatially uniform.     

Sources of δ18O Concentration Variability in Greenland Ice Cores: Temperature, North Atlantic Oscillation and Solar Activity

This work concerns possible connection between variations of δ¹⁸O in Greenland ice, temperature, precipitation amount, NAO indices and relative sunspot numbers. The relationships between 10-year averaged values of this data were derived making use of regression analysis methods. It was discovered that February temperature yields the main contribution in variability of δ¹⁸O. The multiple correlation coefficient between weighted night temperature of the most important months and δ¹⁸O equals 0.95±0.01. Besides the local source of δ¹⁸O changes, defined by the amount of precipitation and temperature, an additional source correlated with the level of solar activity has been revealed. The relative contribution of this source in the δ¹⁸O variance is about 0.1. It has been shown that the multiple correlation coefficient between δ¹⁸O and NAO indices is equal to 0.85. Therefore the NAO yields more than 2/3 of the variances of climatic change in Greenland. An independent contribution of solar activity in climatic changes is about 1/3.