Recent studies have reported shallow and deep seep areas offshore Mocha island. Gas hydrate occurrences along the Chilean margin could explain seeps presence. Gas phases (gas hydrate and free gas) and geothermal gradients were estimated analysing two seismic sections. Close to Mocha island (up to 20 km) were detected high (up to 1900 m/s) and low (1260 m/s) velocities associated with high gas hydrate (up to 20 % of total volume) and free gas (up to 1.1% of total volume) concentrations respectively. These values are in agreement with a variable and high geothermal gradient (65 to 110 °C/km) related to high supply deep fluids canalised by faults and fractures. Faraway from Mocha island (more than 60 km), free gas concentrations decrease to 0.3 % of total volume and low geothermal gradient (from 35 to 60 °C/km) are associated with low fluids supply. Finally, we propose gas hydrate dissociation processes as the main supply source for seeps in the vicinity of Mocha island. These processes can be triggered by ancient sliding reported in literature.

The use of salt rock for underground radioactive waste disposal facilities requires a comprehensive analysis of creep-damage process in salt rock. A computer-controlled creep setup is employed to carry out a creep test of salt rock lasted as long as 359 days under a constant uniaxial stress. The AE space-time evolution and energy releasing characteristics during creep test are studied in the meantime. A new creep-damage model is proposed on the basis of fractional derivative by combining the AE statistical regularity. It indicates that the AE data in non-decay creep process of salt rock can be divided into three stages. Furthermore, the parameters of new creep-damage model are determined by Quasi-Newton method. The fitting analysis suggests that the creep-damage model based on fractional derivative in this paper provides a precise description of full creep regions in salt rock.

This study presents an analysis of the spatial and temporal distribution of the two large destructive earthquakes that occurred in Lake Van area on October 23, and November 9, 2011, together with the azimuth-dependent distribution of the seismic activity and microseismicity clusters after the mainshocks, associated with the complex rupture processes of their aftershock sequence. The sequence began with the magnitude Mw 7.1 earthquake of 23 October and a second destructive earthquake of Mw 5.6. The aftershock sequences of the two mainshocks were linked to the local crustal faults beneath Lake Van area, followed successively and produced unusually intense activity and significant damage in the area. The main purposes of this study are to document the spatial and temporal distribution and evolution of the October 23, 2011 aftershock hypocenters and the azimuth-dependent distribution of seismic activity, and to understand the spatial and temporal character of the aftershock sequence using the distributional and evolutional patterns of the aftershock hypocenters. A total of 10,000 aftershocks were obtained from seismic data with a high signal-to-noise ratio over collected over three years from October 23, 2011 to March 2014. These aftershocks were plotted for the time periods from November 2011 through March 2012 to March 2014 and ≈ 5000 aftershocks were retained in the depth versus distance cross-sections to detect the clusters in the first step of study (November 2011–March 2012). The focal depth distribution of the aftershock clusters, the migration of hypocenter activity and microseismicity clusters were analyzed and the distributional patterns of the detected clusters were assessed using the geometric distribution of the aftershock hypocenters. The spatial and temporal distribution of aftershocks reveal interesting key features of the deep rupture complexity of the Van earthquake: (1) most prominent aftershocks have been located in the upper crust at depths shallower than 10 km beneath ruptured area, indicating that the upper crust is brittle and seismogenic; (2) two spatial clusters have been detected at 8-10 km depths and the upward extrapolation of these clusters intersects with faults; the main cluster (60 km wide) bounded by inferred reverse faults (f3 and f4) and the central cluster (25–30 km wide) bounded by faults (f1 and f2); (3) these spatial clusters form the largest volumetric pattern of the conical-shaped cluster at depths of about 25–30 km of the azimuth-dependent rotational projections, suggesting azimuthal distributions of deep rupture characteristics; and (4) the strongest temporal cluster of microseismicity derived from temporal distribution of aftershocks has been detected within an area of about 2.5–3.0 km² and it is spatially observed at 20 km depth within the central cluster, suggesting progressive failure of the adjacent patches of possible fault.

This paper deals with an analysis of cartographical studies, the real estate cadastre, and its practical implementation, as well as the introduction of cartography into different education modules in university-level studies in Montenegro. There is a discussion of the development, production, and creativity in the fields of cartography and real estate cadastre over time, cartographical projection, scientific results, and recent changes such as advanced computer- and satellite-based technologies, GIS, cartographical visualization, and digital cartography. The impact of these changes on cartographical studies at the University of Montenegro is considered. Particular attention is given to analyses of cartography and the cadastre in institutions, and their connection with the development of cartography teaching modules of Geography, Geodesy and Geoinformatics at the University of Montenegro. The integrated analysis also covers the results of the questionnaire and the significance of the geo-topographical and cartographical heritage of Montenegro, with the aim of carrying these out. It can be seen that the tasks solved by using maps have tended to become more complex and that the cartographical methods employed in this have been always directed towards Montenegro’s most prominent and most urgent problems, including those that appear in the area of education.

Land surface albedo determines the splitting of downwelling solar radiation into components which are either reflected back to the atmosphere or absorbed by the surface. Land surface albedo is an important variable for the climate community and therefore was defined by the Global Climate Observing System (GCOS) as an Essential Climate Variable (ECV). Within the scope of the Satellite Application Facility for Land Surface Analysis (LSA SAF) of EUMETSAT, a near-real time (NRT) daily albedo product was developed in the last decade from observations provided by the SEVIRI instrument on board the geostationary satellites of the Meteosat Second Generation (MSG) series. In this study we present a new collection of albedo satellite products based on the same satellite data. The MSG Ten-day ALbedo (MTAL) product incorporates MSG observations over 31 days with a frequency of NRT production of 10 days. The MTAL collection is more dedicated to climate analysis studies compared to the daily albedo that was initially designed for the weather prediction community. For this reason, a homogeneous reprocessing of MTAL was done in 2018 to generate a Climate Data Record (CDR). The resulting product is called MTAL-R and has been made available to the community in addition to the NRT version of the MTAL product which has been available for several years. The retrieval algorithm behind the MTAL products comprises three distinct modules: one for atmospheric correction, one for daily inversion of a semi-empirical model of the bidirectional reflectance distribution function, and one for monthly composition that also determines surface albedo values. In this study the MTAL-R CDR is compared to ground surface measurements and concomitant albedo products collected by sensors on-board polar-orbiting satellites (SPOT-VGT and MODIS). We show that MTAL-R meets the quality requirements if MODIS or SPOT-VGT are considered as reference. This work leads to 14 years of production of geostationary land surface albedo products with a guaranteed continuity in the LSA SAF for the future years with the forthcoming third generation of European geostationary satellites.

Landslide is a sliding movement of rock mass, debris and soil along the slope under the action of gravity. Small Baseline Subset (SBAS) is an established method for the investigation and monitoring of landslide moments. This study focuses on monitoring the long-temporal displacement of mountainous terrain in Danba County, Sichuan Province via SBAS technique, based on 31 scenes of L-band ALOS/PALSAR data from Feb. 2007 to Oct. 2010.The results show that the largest velocity rates in LOS direction are ±120 mm/yr and maximum accumulated displacement is up to -300, which indicates fast movement of the mountainous terrain during the observation time. These results get good consistency against the results of previous study. This demonstrates the strong potential of SBAS technique for monitoring the landslides geohazard.

In this study, regression analyses were used to find a relationship between the rain gauge rainfall rate R and radar reflectivity Z for the urban catchment of the Służewiecki Stream in Warsaw, Poland. Rainfall totals for 18 events which were measured at two rainfall stations were used for these analyses. Various methods for determining calculational values of radar reflectivity in reference to specific rainfall cells with 1-km resolution within an event duration were applied. The influence of each of these methods on the Z-R relationship was analyzed. The correction coefficient for data from the SRI (Surface Rainfall Intensity) product was established, in which the values of rainfall rate are calculated based on parameters a and b determined by Marshall and Palmer. Relatively good agreement between measured and estimated rainfall totals for the analyzed events was obtained using the Z-R relationships as well as the correction coefficient determined in this study. Rainfall depths estimated from radar data for two selected events were used to simulate flow hydrographs in the catchment using the SWMM (Storm Water Management Model) hydrodynamic model. Different scenarios were applied to investigate the stream response to changes in rainfall depths, in which the data both for 2 existing as well as 64 virtual rain gauges assigned to appropriate rainfall cells in the catchment were included.

The main features of the instrumental global mean surface temperature (GMST) are reasonably well described by a simple linear response model driven by anthropogenic, volcanic and solar forcing. This model acts as a linear long-memory filter of the forcing signal. The physical interpretation of this filtering is the delayed response due to the thermal inertia of the ocean. This description is considerably more accurate if El Ni\~{n}o Southern Oscillation (ENSO) and the Atlantic Multi-decadal Oscillation (AMO) are regarded as additional forcing of the global temperature and hence subject to the same filtering as the other forcing components. By considering these as predictors in a linear regression scheme, more than 92\% of the variance in the instrumental GMST over the period 1870-2017 is explained by this model, and in particular all features of the 1998 -- 2015 hiatus, including its death. While the more prominent pauses during 1870 -- 1915 and 1940 -- 1970 can be attributed to clustering in time of strong volcanic eruptions, the recent hiatus is an unremarkable phenomenon that is attributed to ENSO and a small contribution from solar activity.

The morphology and the activity of a submarine caldera, Avyssos, at the northern part of Nisyros volcano in the South Aegean Sea (Greece), has been studied by means of remotely operated underwater vehicle dives. The recorded time series of temperature and conductivity over the submarine volcano have been analyzed in terms of the Generalized Moments Method. The findings of the mathematical analysis shed light on the volcanic activity, but also on the morphology (shape) of the submarine volcano. The conductivity time series indicates the volcano is at rest in agreement with other types of observations. On the other hand, temperature fluctuations, which in general describe a multifractal process, show that the submarine caldera operates as an open system that interacts with its surroundings. This type of analysis can be used as an indicator for the state of activity and the morphological structure (closed or open system) of a submarine volcano.

This paper examines the tropical cyclone (TC) activity in Solomon Islands (SI) using the best track data from Tropical Cyclone Warning Centre Brisbane and Regional Specialized Meteorological Centre Nadi. The long-term trend analysis showed that the frequency of TCs has been decreasing in this region while average TC intensity becomes strong. Then, the datasets were classified according to the phase of Madden-Julian Oscillation (MJO) and the index of El Nino Southern Oscillation (ENSO) provided by Bureau of Meteorology. The MJO has sufficiently influenced TC activity in the SI region with more genesis occurring in phases 6-8, in which the lower outgoing longwave radiation indicates enhanced convective activity. In contrast, TC genesis occurs less frequently in phases 1, 2, and 5. As for the influence of ENSO, more TCs are generated in El Nino period. The TC genesis locations during El Nino (La Nina) period were significantly displaced to the north (south) over SI region. TCs generated during El Nino condition tended to be strong. This paper also argues the modulation in terms of seasonal climatic variability of large-scale environmental conditions such as sea surface temperature, low level relative vorticity, vertical wind shear, and upper level divergence.

Fracture density, a critical parameter of unconventional reservoirs, can be used to evaluate potential of unconventional reservoirs and location of production wells. Many technologies, such as amplitude variation with offset and azimuth (AVOA) technology, vertical seismic profiling (VSP) technology, and multicomponent seismic technology, are generally used to predict fracture of reservoirs. they can qualitatively predict fracture by analyzing seismic attributes, including seismic wave amplitudes, seismic wave velocities, which are sensitive to fracture. However, it is important to quantitatively describe fracture of reservoirs. In this study, based on a double-layer model, the relationships between fracture density and the double-layer model’s physical parameters, such as velocity of fast shear-wave, velocity of slow shear-wave, and density, were established, and then a powerful quantitative prediction method for fracture density was proposed dramatically. Afterwards, the Hudson model for crack was used to test the applicability of the method. The result shown that the quantitative prediction method for fracture density can be applied suitable to the Hudson model for crack. Finally, the result of validation models indicated that the method can predict fracture density effective, in which absolute relative deviation (ARD) were less than 5% and root-mean-square error (RMSE) was 4.88×10-3.

We use earthquake ground motion modelling via Ground Motion Prediction Equations (GMPEs) and numerical simulation of seismic waves to consider the effects of site amplification and basin resonance in Jakarta, the capital city of Indonesia. While spectral accelerations at short periods are sensitive to near-surface conditions (i.e., Vs30), our results suggest that, for basins as deep as Jakarta’s, available GMPEs cannot be relied upon to accurately estimate the effect of basin depth on ground motions at long periods (>1 s). Amplitudes at such long periods are influenced by entrapment of seismic waves in the basin, resulting in longer duration of strong ground motion, and interference between incoming and reflected waves as well as focusing at basin edges may amplify seismic waves. In order to simulate such phenomena in detail, a basin model derived from a previous study is used as a computational domain for deterministic earthquake scenario modeling in a 2-dimensional cross-section. A M_w 9.0 megathrust, a M_w 6.5 crustal thrust and a M_w 7.0 instraslab earthquake are chosen as scenario events that pose credible threats to Jakarta, and the interactions with the basin of seismic waves generated by these events were simulated. The highest PGV amplifications are recorded at sites near the middle of the basin and near its southern edge, with maximum amplifications of PGV in the horizontal component of 200% for the crustal, 600% for the megathrust and 335% for the deep intraslab earthquake scenario, respectively. We find that the levels of ground motion response spectral acceleration fall below those of the 2012 Indonesian building Codes's design response spectrum for short periods (< 1 s), but closely approach or may even exceed these levels for longer periods.

This paper provides a comparative study between microwave tomography and synthetic time-reversal imaging techniques as applied to ground penetrating radar (GPR) surveys. The comparison is carried out by processing experimental data collected at a controlled test site, with various types of buried targets at given subsurface depths and representative soil conditions. It is shown that the two techniques allow us to obtain complementary information about position, depth and size of the targets from a single GPR survey.

We use Matlab 3D finite element fluid flow/transport modelling to simulate localized wellbore temperature events of order 0.05-0.1oC logged in Fennoscandia basement rock at ~ 1.5km depths. The temperature events are approximated as steady-state heat transport due to fluid draining from the crust into the wellbore via naturally occurring fracture-connectivity structures. Flow simulation is based on the empirics of spatially-correlated fracture-connectivity fluid flow widely attested by well-log, well-core, and well-production data. Matching model wellbore-centric radial temperature profiles to a 2D analytic expression for steady-state radial heat transport with Peclet number Pe ≡ r0φv0/D (r0 = wellbore radius, v0 = Darcy velocity at r0, φ = ambient porosity, D = rock-water thermal diffusivity), gives Pe ~ 10-15 for fracture-connectivity flow intersecting the well, and Pe ~ 0 for ambient crust. Darcy flow for model Pe ~ 10 at radius ~ 10 meters from the wellbore gives permeability estimate κ ~ 0.02Darcy for flow driven by differential fluid pressure between least principal crustal stress pore pressure and hydrostatic wellbore pressure. Model temperature event flow permeability κm ~ 0.02Darcy is related to well-core ambient permeability κ ~ 1µDarcy by empirical poroperm relation κm ~ κ exp(αmφ) for φ ~ 0.01 and αm ~ 1000. Our modelling of wellbore temperature events calibrates the concept of reactivating fossilized fracture-connectivity flow for EGS permeability stimulation of basement rock.

Papua is one of part in Indonesia which is the geology research of that place isn’t developed and limited. It causes the seismotectonic of Papua hasn’t been known. WinITDB was used to determine the dip angle plate which was on the north part of Papua. The determination of angle was done through seismicity’s cross section analysis in the area. To show that seismicity, earthquake history data that ever occurred in the area is needed. The result on the seismicity’s cross section of plane A–A’, was confluence by two plates with angle 150° against horizontal on the depth up to ±68 km. On the seismicity’s cross section of plane B–B’ had angle 135° against horizontal on the ±82 km depth. On the plane C–C’ seismicity’s cross section, was confluence of two plate which located between -1,77°S until -4,97°S subducted until 171 km depth on 1,38°N - 4,97°S. It proved that subduction characteristic in the northern Papua which was Australia continent plate subducted to north, followed by collision and the Pacific plate subduction on New Guinea. It is also confirmed by focus mechanism analysis which showed the earthquake activities are controlled by the not really deep active fault.

Channel bends are one of the most important characteristic features of natural streams. These bends often create the conditions for a hyporheic zone, which has been recognized as a critical component of stream ecosystems. The streambed vertical hydraulic conductivity (K_v), vertical hydraulic gradient (VHG) and Darcian flux (DF) in the hyporheic zone were estimated at 61 locations along a channel bend of the Beiluo River during July 2015 and January 2016. All the streambed attributes showed great spatial variability along the channel bend. Both upward fluxes and downward fluxes occurred during the two test periods, most of studied stream sections were controlled by downwelling, indicating stream water discharge into the subsurface. The average downward flux was higher at the downstream side than at the upstream side of the channel bend, especially in July 2015. The distribution of streambed sediment grain size has a significant influence on the variability of K_v; high percentages of silt and clay sediments generally lead to low K_v values. Higher K_v at the depositional left bank at the upstream site shifted toward the erosional right bank at the downstream site, with K_v values positively correlated with the water depth. This study suggested that the variabilities of K_v and VHG were influenced by the stream geomorphology and that the distribution of K_v was inversely related, to a certain extent, to the distribution of VHG across the channel bend. K_v and VHG were found to have opposite effects on the DF, and the close relationship between K_v and DF indicated that the water fluxes were mainly controlled by K_v.

In the paper the verification of forecasts of precipitation conditions measured by the standardized precipitation index SPI is presented. For the verification of categorical forecasts a contingency table was used. Standard verification measures were used for the SPI value forecast. The 30 day SPI moved every 10 days by 10 days was calculated in 2013-2015 from April to September on the basis of precipitation data from 35 meteorological stations in Poland. Predictions of the 30 day SPI were created in which precipitation was forecasted in the next 10 days (the SPI 10-day forecast) and 20 days (the SPI 20-day forecast). Both for the 10 and 20 days, the forecasts were skewed towards drier categories at the expense of wet categories. There was a good agreement between observed and 10-day forecast categories of precipitation. Less agreement is obtained for 20-day forecasts – these forecasts evidently “over-dry” the assessment of precipitation anomalies. The 10-day SPI value forecast accuracy is acceptable, whereas for the 20-day forecast is unsatisfactory. Both for the SPI categorical and the SPI value forecast, the 10-day SPI forecast is reliable and the 20-day forecast should be accepted with reservation and used with caution.

As the most active plateau on the Earth, the Qinghai-Tibet Plateau has a complex crust-mantle structure. Knowledge of the distribution of such a structure provides information for understanding the underlying geodynamic processes. We obtains a three-dimensional density model of crustal and upper mantle beneath Qinghai-Tibet plateau and its surrounding areas from the residual geoid anomalies using the Earth Gravitational Model (EGM) 2008. We estimate a refined density model by iterations, using an initial density contrast model. We confirm that the EGM2008 mission products can be used to constrain the crust-mantle density structures. Our major findings are: (1). At 300-400 km depth, high-D anomalies terminate around Jinsha River Suture (JRS) in the central TP, suggesting that the Indian plate has been reached over the Bangong Nujiang Suture (BNS) and almost reach to the JRS. (2). On the eastern TP, low-D anomalies at the depth of 0-300 km together with high-D anomalies at 400-670 km further verified the current eastward subduction of Indian plate. The ongoing subduction provides forces to the occurrences of frequent earthquakes and volcano. (3). At 600 km depth, low-D anomalies inside the TP illustrate the existence of hot weak material beneath there, contributing to the external material inward-thrusting.

Climate change is a global issue that draws widespread attention from the international society. As an important component of the climate system, the water cycle is directly affected by climate change. Thus, it is very important to study the influences of climate change on the basin water cycle with respect to maintenance of healthy rivers, sustainable use of water resources, and sustainable socioeconomic development in the basin. In this study, by assessing the suitability of multiple General Circulation Models (GCMs) recommended by the Intergovernmental Panel on Climate Change, Statistical Downscaling Model (SDSM) and Automated Statistical Downscaling model (ASD) were used to generate future climate change scenarios. These were then used to drive distributed hydrologic models (Variable Infiltration Capacity, Soil and Water Assessment Tool) for hydrological simulation of the Yangtze River and Yellow River basins, thereby quantifying the effects of climate change on the basin water cycle. The results showed that suitability assessment adopted in this study could effectively reduce the uncertainty of GCMs, and that statistical downscaling was able to greatly improve precipitation and temperature outputs in global climate mode. Compared to a baseline period (1961–1990), projected future periods (2046–2065 and 2081–2100) had a slightly decreasing tendency of runoff in the lower reaches of the Yangtze River basin. In particular, a significant increase in runoff was observed during flood seasons in the southeast part. However, runoff of the upper Yellow River basin decreased continuously. The results provide a reference for studying climate change in major river basins of China.

Eddy-driven jets are of importance in the ocean and atmosphere, and to a first approximation are governed by Rossby wave dynamics. This study addresses the time-dependent flux of fluid and potential vorticity between such a jet and an adjacent eddy, with specific regard to determining zonal and meridional wavenumber dependence. The flux amplitude in wavenumber space is obtained, which is easily computable for a given jet geometry, speed and latitude, and which provides instant information on the wavenumbers of the Rossby waves which maximize the flux. This new tool enables the quick determination of which modes are most influential in imparting fluid exchange, which in the long term will homogenize the potential vorticity between the eddy and the jet. The results are validated by computing backward- and forward-time finite-time Lyapunov exponent fields, and also stable and unstable manifolds; the intermingling of these entities defines the region of chaotic transport between the eddy and the jet. The relationship of all of these to the time-varying transport flux between the eddy and the jet is carefully elucidated. The flux quantification presented here works for general time-dependence, whether or not lobes (intersection regions between stable and unstable manifolds) are present in the mixing region, and is therefore also easily computable for wave packets consisting of infinitely many wavenumbers.