Temporal and spatial variations of potential evaporation and the driving mechanism over Tibet during 1961–2001

Potential evaporation (PE) is the basic component of the global hydrological cycle and energy balance. This study detected the temporal and spatial variations of PE and related driving factors in Tibet, China, for the period 1961–2001, based on observed data recorded at 22 meteorological stations. The results showed that (1) Tibet experienced a statistically significant decrease of PE between 1961 and 2001, which started mainly in the 1980s, along with accelerated warming; (2) the mean annual PE in Tibet showed an east–west increasing trend, and the annual PE at most stations presented decreasing trends; (3) an inverse correlation of mean annual PE with elevation was detected (low–high decreasing trend), and the statistical equations to estimate PE were established based on longitude, latitude and elevation; and (4) PE in Tibet can be well expressed by related meteorological variables, with vapour pressure deficit the dominant factor in determining PE.

EDITOR Z. W. Kundzewicz ASSOCIATE EDITOR not assigned     

Temporal and spatial dynamical simulation of groundwater characteristics in Minqin Oasis

Application scope of geostatistics has been gradually extended from original geologic field to soil science and ecological field, etc. and its successful application results have been widely demonstrated. But little information is reported as to the direct use of geostatistical method to work out the distribu- tion map of groundwater characteristics. In this paper the semivariogram of geostatistics, in combina- tion with GIS, was used to quantitatively study the spatial variation characteristics of groundwater table depth and mineralization degree and their relation to the landuse changes. F test of the used spherical model reached a very significant level, and the theoretical model can well reflect the spatial structural characteristics of groundwater table depth and mineralization degree and achieve an ideal result. This shows that the application of the method in the dynamical simulation of groundwater is feasible. And this paper also provides useful reference for the application of geostatistics in the study of the dy- namical variations of groundwater resources in the oasis.     

Changing structures of summertime heatwaves over China during 1961–2017

Despite the prevalence of artificial separation of daytime and nighttime hot extremes, they may actually co-occur or occur sequentially. Considering their potential lead-lag configuration, this study identified an entire heatwave period as consecutive days with either daytime or nighttime hot extremes and investigated the changes of the prevalence and sequence of daytime and nighttime hot extremes during heatwaves over China from 1961 to 2017. It was found that the majority(82%) of heatwaves were compound heatwaves that had both daytime and nighttime hot extremes exceeding the 90 th percentile-based thresholds, while only 7%(11%) were purely daytime(nighttime) heatwaves that contained only daytime(nighttime) hot extremes. During the entire periods of compound heatwaves, daytime hot extremes usually occurred one day or a few days before nighttime hot extremes, which was in accordance with the daily variations in radiation and meteorological conditions, such as the increasing surface humidity and cloud cover, and decreasing solar radiation during the entire heatwave periods. From 1961 to2017, compound heatwave numbers exhibited the sharpest increase with a statistically significant trend of 0.44 times decade-1, in contrast to an insignificant trend of 0.00 times decade-1 for purely daytime heatwaves and a significant trend of 0.09 times decade-1 for purely nighttime heatwaves. Within the compound heatwave periods, hot nights were starting earlier and ending later, and numbers of concurrent daytime-nighttime hot extremes increased significantly at 0.20 days decade-1. In particular,urban area were not only subject to increasingly more frequent and longer compound heatwaves, but also to more occurrences of concurrent daytime-nighttime hot extremes with more serious impact. This study provides instructions for researchers to customize and select appropriate heatwave indices.     

Planetary waves observed by TIMED/SABER in coupling the stratosphere–mesosphere–lower thermosphere during the winter of 2003/2004: Part 2—Altitude and latitude planetary wave structure

Part 2 of the present paper is focused on the planetary wave coupling from the stratosphere to the lower thermosphere (30–120 km) during the Arctic winter of 2003/2004. The planetary waves seen in the TIMED/SABER temperature data in the latitudinal range 50°N–50°S are studied in detail. The altitude and latitude structures of the planetary wave (stationary and travelling) clearly indicate that the stratosphere and mesosphere (30–90 km) are coupled by direct vertical propagation of the planetary waves, while the lower thermosphere (above 90–95 km altitude) is only partly connected with the lower levels probably indirectly through in-situ generation of disturbances by the dissipation and breaking of gravity waves filtered by lower atmospheric planetary waves. A peculiar feature of the thermal regime in the lower thermosphere is that it is dominated by zonally symmetric planetary waves.     

Analysis of space–time patterns of rainfall events during 1996–2008 in Yilan County (Taiwan)

Stochastic Environmental Research and Risk Assessment - Understanding local precipitation patterns is essential to water resource management and flood mitigation. Precipitation patterns can vary in...     

The 20–30-day oscillation of the global circulation and heavy precipitation over the lower reaches of the Yangtze River valley

Based on the observational data in summer, the variations of intraseasonal oscillation (ISO) of the daily rainfall over the lower reaches of the Yangtze River valley (LYRV) were studied by using the non-integer spectrum analysis. The NCEP/NCAR reanalysis data for the period of 1979–2005 were analyzed by principal oscillation pattern analysis (POP) to investigate the spatial and temporal characteristics of principal ISO patterns of the global circulation. The relationships of these ISO patterns to the rainfall ISO and the heavy precipitation process over LYRV were also discussed. It is found that the rainfall over LYRV in May–August is mainly of periodic oscillations of 10–20, 20–30 and 60–70 days, and the interannual variation of the intensity of its 20–30-day oscillation has a strongly positive correlation with the number of the heavy precipitation process. Two modes (POP1, POP2) are revealed by POP for the 20–30-day oscillation of the global 850 hPa geopotential height. One is a circumglobal teleconnection wave train in the middle latitude of the Southern Hemisphere (SCGT) with an eastward propagation, and the other is the southward propagation pattern in the tropical western Pacific (TWP). The POP modes explain 7.72% and 7.66% of the variance, respectively. These two principal ISO patterns are closely linked to the low frequency rainfall and heavy precipitation process over LYRV, in which the probability for the heavy precipitation process over LYRV is 54.9% and 60.4% for the positive phase of the imaginary part of POP1 and real part of POP2, respectively. Furthermore, the models of the global atmospheric circulation for the 20–30-day oscillation in association with or without the heavy precipitation process over LYRV during the Northern Hemisphere summer are set up by means of the composite analysis method. Most of the heavy precipitation processes over LYRV appear in Phase 4 of SCGT or Phase 6 of TWP. When the positive phases of 20–30-day oscillations for the rainfall over LYRV are associated with (without) the heavy precipitation process, a strong westerly stream appears (disappears) from the Arabian Sea via India and Bay of Bengal (BOB) to southern China and LYRV for the global 850 hPa filtered wind field during Phase 4 of SCGT. This situation is favorable (unfavorable) for the forming of the heavy precipitation process over LYRV. Similarly, a strong (weak) western wind belt forms from India through BOB to southern China and LYRV and the subtropical northwestern Pacific and central and eastern equatorial Pacific during Phase 6 of TWP for the cases with (without) the heavy precipitation process. The evolutions of these ISO patterns related to the 20–30-day oscillation are excited by either the interaction of extratropical circulation in both hemispheres or the heat source forcing in Asia monsoon domain and internal interaction of circulation in East Asia. These two global circulation models might therefore provide valuable information for the extended-range forecastof the heavy precipitation process over LYRV during the 10–30 days.     

Study on spatial pattern of land-use change in China during 1995–2000

It is more and more acknowledged that land-use/cover dynamic change has become a key subject urgently to be dealt with in the study of global environmental change. Supported by the Landsat TM digital images, spatial patterns and temporal variation of land-use change during 1995 –2000 are studied in the paper. According to the land-use dynamic degree model, supported by the 1km GRID data of land-use change and the comprehensive characters of physical, economic and social features, a dynamic regionalization of land-use change is designed to disclose the spatial pattern of land-use change processes. Generally speaking, in the traditional agricultural zones, e.g., Huang-Huai-Hai Plains, Yangtze River Delta and Sichuan Basin, the built-up and residential areas occupy a great proportion of arable land, and in the interlock area of farming and pasturing of northern China and the oases agricultural zones, the reclamation of arable land is conspicuously driven by changes of production conditions, economic benefits and climatic conditions. The implementation of “returning arable land into woodland or grassland” policies has won initial success in some areas, but it is too early to say that the trend of deforestation has been effectively reversed across China. In this paper, the division of dynamic regionalization of land-use change is designed, for the sake of revealing the temporal and spatial features of land-use change and laying the foundation for the study of regional scale land-use changes. Moreover, an integrated study, including studies of spatial pattern and temporal process of land-use change, is carried out in this paper, which is an interesting try on the comparative studies of spatial pattern on change process and the change process of spatial pattern of land-use change.     

Comparison of model electron densities and temperatures with Millstone Hill observations during undisturbed periods and the geomagnetic storms of 16–23 March and 6–12 April 1990

Measurements of F-region electron density and temperature at Millstone Hill are compared with results from the IZMIRAN time-dependent mathematical model of the Earths ionosphere and plasmasphere during the periods 16–23 March and 6–12 April 1990. Each of these two periods included geomagnetically quiet intervals followed by major storms. Satisfactory agreement between the model and the data is obtained during the quiet intervals, provided that the recombination rate of O⁺(⁴S) ions was decreased by a factor of 1.5 at all altitudes during the nighttime periods 17–18 March, 19–20 March, 6–8 April and 8–9 April in order to increase the NmF2 at night better to match observations. Good model/data agreement is also obtained during the storm periods when vibrationally excited N₂ brings about factor-of-2-4 reductions in daytime NmF2. Model calculations are carried out using different expressions for the O⁺ – O collision frequency for momentum transfer, and the best agreement between the electron-density measurements and the model results is obtained when the CEDAR interim standard formula for the O⁺ – O collision frequency is used. Deviations from the Boltzmann distribution for the first five vibrational levels of NI were calculated. The calculated distribution is highly non-Boltzmann at vibrational levels j > 2, and the Boltzmann distribution assumption results in the increase of 10–30% in calculated NmF2 during the storm-time periods. During the March storm at solar maximum the model results obtained using the EUVAC solar flux model agree a little better with the observations in comparison with the EUV94 solar flux model. For the April storm period of moderate solar activity the EUV94X model results agree better with the observations in comparison to the EUVAC model.     

Planetary waves observed by TIMED/SABER in coupling the stratosphere–mesosphere–lower thermosphere during the winter of 2003/2004: Part 1—Comparison with the UKMO temperature results

Part 1 of the present paper is focused on the types of planetary wave seen in the TIMED/SABER and UK Met Office (UKMO) temperature data in the Northern Hemisphere (NH) (0–50°N) stratosphere (30–60 km) during the Arctic winter of 2003/2004, as the emphasis is on their spatial structure (latitude and altitude) and temporal evolution particularly in relation to the stratospheric warmings. A new method for analysis of satellite data is presented in this study where the migrating and nonmigrating tides and planetary waves (stationary, zonally symmetric and travelling) are simultaneously extracted from the satellite data. The comparison between the altitude and latitude structure of the SABER and UKMO planetary waves in the temperature field of the NH stratosphere indicates a high degree of qualitative and quantitative resemblance and in this way the validity of the new data analysis method is verified as well.     

Variability of aerosol properties during the 2007–2010 spring seasons over central Europe

Aerosol optical properties have been studied for spring seasons when increased values of PM10 are registered. Measurements of aerosol optical properties were taken by collocated lidar and sun-photometers at Belsk, Poland, and Minsk, Belarus. A significant increase of registered aerosol optical thickness (AOT) was found during episodes with elevated PM10 concentrations. An increase of AOT at 1020 nm amounted to 50% in the case of Minsk and 18% in the case of Belsk, while an increase of AOT at 400 nm was 66% and 33%, respectively. We noted an increase of Ångström exponent by 6% at both stations and no significant increase of single scattering albedo. The LIDAR measurements together with NAAPS model results and backtrajectory analysis suggest that both the biomass burning products and the Saharan dust are responsible for increased PM10 concentrations and large AOT values during spring time. The smoke aerosol is transported over Central Europe mainly in the boundary layer, increasing both PM10 concentration and AOT. The dust aerosol transported in the free troposphere slightly affects the AOT values only. Statistically significant correlation between PM10 concentration and AOT was found during reporting period.     

Specific patterns of XCO2 observed by GOSAT during 2009–2016 and assessed with model simulations over China

Spatiotemporal patterns of column-averaged dry air mole fraction of CO2(XCO2)have not been well characterized on a regional scale due to limitations in data availability and precision.This paper addresses these issues by examining such patterns in China using the long-term mapping XCO2 dataset(2009-2016)derived from the Greenhouse gases Observing SATellite(GOSAT).XCO2 simulations are also constructed using the high-resolution nested-grid GEOS-Chem model.The following results are found:Firstly,the correlation coefficient between the anthropogenic emissions and XCO2 spatial distribution is nearly zero in summer but up to 0.32 in autumn.Secondly,on average,XCO2 increases by 2.08 ppm every year from2010 to 2015,with a sharp increase of 2.6 ppm in 2013.Lastly,in the analysis of three typical regions,the GOSAT XCO2 time series is inbetter agreement with the GEOS-Chem simulation of XCO2 in the Taklimakan Desert region(the least difference with bias 0.65±0.78 ppm),compared with the northern urban agglomerationregion(-1.3±1.2 ppm)and the northeastern forest region(-1.4±1.4 ppm).The results are likely attributable to uncertainty in both the satellite-retrieved XCO2 data and the model simulation data.     

Characteristics of Changes in the Ozone Content in the Upper Stratosphere over Moscow during the Cold Half-Years of 2014–2015 and 2015–2016

Geomagnetism and Aeronomy - The results of ground-based microwave measurements of the stratospheric ozone profiles over Moscow during the cold half-years of 2014–2015 and 2015–2016 are...     

The role of vibrationally excited oxygen and nitrogen in the ionosphere during the undisturbed and geomagnetic storm period of 6–12 April 1990

We present a comparison of the observed behavior of the F-region ionosphere over Millstone Hill during the geomagnetically quiet and storm periods of 6–12 April 1990 with numerical model calculations from the IZMIRAN time-dependent mathematical model of the Earths ionosphere and plasmasphere. The major enhancement to the IZMIRAN model developed in this study is the use of a new loss rate of O⁺(⁴S) ions as a result of new high-temperature flowing afterglow measurements of the rate coefficients K₁ and K₂ for the reactions of O⁺(⁴S) with N₂ and O₂. The deviations from the Boltzmann distribution for the first five vibrational levels of O₂(v) were calculated, and the present study suggests that these deviations are not significant. It was found that the difference between the non-Boltzmann and Boltzmann distribution assumptions of O₂(v) and the difference between ion and neutral temperature can lead to an increase of up to about 3% or a decrease of up to about 4% of the calculated NmF2 as a result of a respective increase or a decrease in K₂. The IZMIRAN model reproduces major features of the data. We found that the inclusion of vibrationally excited N₂(v > 0) and O₂(v > 0) in the calculations improves the agreement between the calculated NmF2 and the data on 6, 9, and 10 April. However, both the daytime and nighttime densities are reproduced by the IZMIRAN model without the vibrationally excited nitrogen and oxygen on 8 and 11 April better than the IZMIRAN model with N₂(v > 0) and O₂(v > 0). This could be due to possible uncertainties in model neutral temperature and densities, EUV fluxes, rate coefficients, and the flow of ionization between the ionosphere and plasmasphere, and possible horizontal divergence of the flux of ionization above the station. Our calculations show that the increase in the O⁺ + N₂ rate factor due to N₂(v > 0) produces a 5–36% decrease in the calculated daytime peak density. The increase in the O⁺ + O₂ loss rate due to vibrationally excited O₂ produces 8–46% reductions in NmF2. The effects of vibrationally excited O₂ and N₂ on N_e and T_e are most pronounced during the daytime.     

Multi-decadal changes in the relationships between rainfall characteristics and debris-flow occurrences in response to gully evolution after the 1990–1995 Mount Unzen eruptions

Explosive volcanic eruptions can cause long-term landscape change, leading to increased sediment discharge that continues after the cessation of the eruptions. During the period 1990–1995, eruptions of Mount Unzen, Japan, generated large amounts of pyroclastic material, resulting in 57 debris-flow events during 1991–2018. To investigate changes in the relationships between rainfall characteristics and debris-flow occurrence, we conducted the following: geometric analysis of two gullies (i.e., debris-flow initiation zones) using LiDAR (light detection and ranging)-generated 1 m DEMs (digital elevation models); rainfall analysis, based on the relationship between rainfall duration and mean intensity (i.e., considering the intensity–duration, or ID, threshold); and debris-flow monitoring during 2016–2018. Since 1991, rainfall runoff has caused erosion of the supplied pyroclastic material, generating a channel network consisting of incised gullies. With sufficient rainfall, debris flows formed, accompanied by further gully erosion; this resulted in both vertical and lateral adjustments of the cross-sectional geometry. In the two decades since the eruptions ceased, readily mobilized pyroclastic material has become scarce as the gullies have adjusted to local hydrographic conditions. At the same time, the infiltration capacity of the volcanic flank has increased, reducing the capacity for overland flow. As a result, since 2000, rainfall events with intensities above the ID threshold have occurred; however, the lack of sediment supplied by the gullies appears to have hindered the occurrence and development of debris flows. This suggests that debris flows in volcanically perturbed landscapes may occur at lower rainfall thresholds as long as the corresponding upland channels are evolving as a result of intense overland flow. However, as such channels evolve towards equilibrium geometries, the frequency of debris flows decreases in response to the reduction in sediment availability.     

Climatological Ertel’s potential-vorticity flux and mean meridional circulation in the extratropical troposphere – lower stratosphere

We investigated to what extent the isentropic, non-geostrophic formulation of zonally averaged circulation derived for stratospheric conditions is applicable to climatological transport in the extratropical troposphere and lower stratosphere. The study is based on 10 years of daily data of ECMWF analysis and on the ECHAM3 climate model of the German Climate Computing Centre. The main result is a scalar isentropic mixing coefficient, K_yy, and a mean meridional transport circulation consistently derived from the same data base. For both data sources, isentropic mean meridional circulation is derived from horizontal mass flow rate for 4 representative months. Alternatively, a mean meridional circulation is calculated from total diabatic heating rates of the ECHAM3 model. It is shown that only the latter is in good agreement with the ECMWF mean meridional circulation. Isentropic analysis also comprises the seasonal cycle of the climatological meridional gradient and flux of Ertels potential vorticity (PV). Application of Tungs flux-gradient relation yields that for all seasons K_yy is positive in height-latitude regions where statistical significance is reached. Large K_yy values, marking regions of more efficient mixing, have been found in the subtropical vertical band of weak westerly wind and in mid-latitudes in regions of upward-propagating baroclinic wave activity in the middle and upper troposphere. Based on the ECMWF data and results of baroclinic-wave behaviour, strong indications are presented that positive zonally averaged PV flux polewards of the jet core in the NH is strengthened by stationary waves and nonlinear effects. Reduced eddy transport is apparent in winter and spring slightly below the subtropical tropopause jet. The seasonal cycle of K_yy from ECHAM3 data is to a great extent in agreement with the result based on ECMWF analysis. In the model, reduced interannual variability enlarges the height-latitude range where sign of K_yy is significant.     

Equatorial spread and dynamics in the F layer over West Africa from ionogram analysis,during the declining solar flux year 1994–1995

Equatorial spread-F (ESF) dynamics is studied based on quarter-hourly series of ionograms recorded in West Africa during the year December 1994 to November 1995. The time evolution of F layer at the typical equatorial site OUAGADOUGOU (Burkina Faso) differs slightly from that at DAKAR (Senegal), located at the northern edge of the magnetic equatorial belt. Comparing ESF occurrences at both stations with the pre-sunrise ones at FORTALEZA (Brazil), we find generally antisymmetric variations on both sides of the Atlantic, except for a common December solstice maximum. We stress the convergence of the multiple-scale morphologies drawn from new radars with those from high-quality ionosonde density (h′,t) plots. This revives the analysis of ESF dynamics and points toward long-term network studies of ESF coupling phenomena.     

Spatio-temporal characteristics of acoustic emission during the deformation of rock samples with compressional and extensional en-echelon faults

Introduction By the experimental techniques of acoustic emission (AE) and surface strain observation, a lot of achievements on the study of deformation and fracture process of samples with en-echelon faults has been obtained (Bomblakis, 1973; LIU, et al, 1986; DU, et al, 1989; LI, et al, 1989; ZHANG, FANG, 1989; ZHANG, LI, 1989; MA, et al, 1995a, b; LI, et al, 2000). However, there still exist many problems. For example, most of discussion on the spatio-temporal evolution of microf…     

Vertical motion and diabatic heating over the Indian monsoon region during the Bay of Bengal depression of 5–8 July, 1979

The vertical velocity, , and the diabatic heating were computed at 800, 600, 400 and 200 mb surfaces using the Omega equation. The highest contribution to is from the diabatic heating produced by condensation associated with the precipitations appearing to be the main source of diabatic heating. The net radiative cooling and the thermal advection in the upper troposphere over the warm anticyclone result in diabatic cooling over the eastern part of the Bay of Bengal and adjoining northern and eastern regions.List of Symbols Used C _p Heat capacity at constant pressure - f Coriolis parameter - g Acceleration due to gravity - P Atmospheric pressure - Q Diabatic heating rate per unit mass - R Gas constant of air - S Static stability parameter - t Time - U, V Zonal and meridional wind components - Specific volume - Relative vorticity - Absolute vorticity - Potential temperature - Geopotential - Vertical velocity (dP/dt) - ₁ Adiabatic vertical velocity - ₂ Vertical velocity due to certain forcing - ₃ Diabatic vertical velocity - Isobaric gradient operator - ² Laplacian operator - J(A, B) Jacobian operator