Permo-Triassic high-pressure metamorphism in the central western Korean Peninsula,and its link to Paleo-Tethyan Ocean closure: Key issues revisited

Permo-Triassic high-pressure(HP) mafic granulites, together with the Bibong retrogressed eclogite,preserved along the central western Korean Peninsula provide important insights into the Late Permian to Triassic collisional orogeny in northeast Asia. The metamorphic pressureetemperatureetime(P-T-t)paths of these rocks, however, remain poorly constrained and even overestimated, owing to outdated geothermobarometers and inaccurate isopleth techniques. Here we evaluate the metamorphic Pe T conditions of Triassic HP mafic granulites including those in Baekdong, Sinri and Daepan and the Bibong Triassic retrogressed eclogite in the Hongseong area, and the Permo-Triassic Samgot mafic granulite in the Imjingang Belt of the central western Korean Peninsula through the application of modern phase equilibria techniques. The Baekdong and Samgot mafic granulites and the Bibong retrogressed eclogite yield a range of 12.0 -16.0 kbar and 800 -900℃, representing HP granulite facies conditions. The Sinri and Daepan granulites from the Hongseong area show relatively lower grade metamorphic conditions between HP granulite and normal granulite facies, and are characterized by sub-isothermal decompression during exhumation. The similarities in the metamorphic ages and the post-collisional igneous activity from the central western Korean Peninsula indicate that the Triassic ages represent the retrograde stage of the metamorphic Pe T paths. In contrast, the Late Permian metamorphic ages, which are older than protolith ages of the post-collisional igneous rocks, correspond to the possible prograde stage of metamorphism. The P-T-t paths presented in this paper, together with the metamorphic ages and post-orogenic igneous events reported from these areas suggest trace of the subduction, accretion and exhumation history, and indicate a tectonic linkage among the northeast Asian continents during the Paleo-Tethyan Ocean closure.     

Evaluation of rock fracture patterns based on the element-free Galerkin method for stability assessment of a highly pressurized gas storage cavern

Reasonable failure paths are crucial in the stability assessment of a highly pressurized gas storage cavern in a rock mass using limit equilibrium analysis. In this study, the fracture patterns in the rock mass around a gas storage cavern are evaluated based on the element-free Galerkin (EFG) method with a cohesive crack model. The proposed analysis method was first validated through scaled-down model testing results before extending to a full-scale problem. Using the current approach, a parametric sensitivity analysis was performed by taking into account the influences of the stress ratio, the depth of the cavern, and the tensile strength of the rock mass. Based on the analysis results obtained, this research showed that the representative failure paths, in which the uplift evaluation is necessary, can be established from the EFG results with careful consideration of crack initiation and the propagation direction. The site having an in situ stress ratio greater than 1.2 is preferred for operation of underground gas storage caverns. By performing the limit equilibrium analyses for stability assessment of the rock mass above the pressurized cavern, we demonstrated that the factor of safety against uplift computed from the failure paths established in this study is higher than or equal to those obtained from analyses with the simplified models in previous studies without considering the influence parameters on the fracture patterns.     

Large-scale rainfall–runoff–inundation modeling for upper Citarum River watershed,Indonesia

The Citarum River is one of the strategic rivers in West Java, Indonesia. Its total watershed area is approximately 1800 km². Almost every year, the overflow from the Citarum River causes the inundation of most of the upper Citarum River watershed. To prevent and mitigate flood damage, it is necessary to understand the flooding characteristics. The region, however, suffers from a lack of observational data. Therefore, to analyze the inundation caused by flooding in the upper Citarum River watershed, a rainfall–runoff–inundation (RRI) model was employed. It used the following multiple satellite-derived datasets as input data as well as for model verification: Global Satellite Mapping of Precipitation, Hydrological data and maps based on Shuttle elevation Derivatives at multiple scales, Global Mosaics of the standard MODIS land cover type data product, and Landsat 7 satellite images. Parameter calibration was performed using a Monte Carlo simulation. The simulation was performed for February 2010. The results of this study show that the RRI model identifies inundation areas in large-scale river watersheds more effectively when using multiple satellite-derived datasets compared with the observed inundation map obtained from JICA in 2010 and Landsat 7 images. The model results can be improved if high-quality observed rainfall data, topographic data, and river cross-sectional data are available.     

A finite difference method for earthquake sequences in poroelastic solids

Induced seismicity (earthquakes caused by injection or extraction of fluids in Earth’s subsurface) is a major, new hazard in the USA, the Netherlands, and other countries, with vast economic consequences if not properly managed. Addressing this problem requires development of predictive simulations of how fluid-saturated solids containing frictional faults respond to fluid injection/extraction. Here, we present a finite difference method for 2D linear poroelasticity with rate-and-state friction faults, accounting for spatially variable properties. Semi-discrete stability and accuracy are proven using the summation-by-parts, simultaneous-approximation-term (SBP-SAT) framework for discretization and boundary condition enforcement. Convergence rates are verified using the method of manufactured solutions and comparison to the analytical solution to Mandel’s problem. The method is then applied to study fault slip triggered by fluid injection and diffusion through high-permeability fault damage zones. We demonstrate that in response to the same, gradual forcing, fault slip can occur in either an unstable manner, as short-duration earthquakes that radiate seismic waves, or as stable, aseismic, slow slip that accumulates over much longer time scales. Finally, we use these simulation results to discuss the role of frictional and elastic properties in determining the stability and nature of slip.     

Landslide susceptibility mapping using random forest and boosted tree models in Pyeong-Chang,Korea

Landslides susceptibility maps were constructed in the Pyeong-Chang area, Korea, using the Random Forest and Boosted Tree models. Landslide locations were randomly selected in a 50/50 ratio for training and validation of the models. Seventeen landslide-related factors were extracted and constructed in a spatial database. The relationships between the observed landslide locations and these factors were identified by using the two models. The models were used to generate a landslide susceptibility map and the importance of the factors was calculated. Finally, the landslide susceptibility maps were validated. Finally, landslide susceptibility maps were generated. For the Random Forest model, the validation accuracy in regression and classification algorithms showed 79.34 and 79.18%, respectively, and for the Boosted Tree model, these were 84.87 and 85.98%, respectively. The two models showed satisfactory accuracies, and the Boosted Tree model showed better results than the Random Forest model.     

Northern landscapes in transition: Evidence,approach and ways forward using the Krycklan Catchment Study

Improving our ability to detect changes in terrestrial and aquatic systems is a grand challenge in the environmental sciences. In a world experiencing increasingly rapid rates of climate change and ecosystem transformation, our ability to understand and predict how, when, where, and why changes occur is essential for adapting and mitigating human behaviours. In this context, long-term field research infrastructures have a fundamentally important role to play. For northern boreal landscapes, the Krycklan Catchment Study (KCS) has supported monitoring and research aimed at revealing these changes since it was initiated in 1980. Early studies focused on forest regeneration and microclimatic conditions, nutrient balances and forest hydrology, which included monitoring climate variables, water balance components, and stream water chemistry. The research infrastructure has expanded over the years to encompass a 6790 ha catchment, which currently includes 11 gauged streams, ca. 1000 soil lysimeters, 150 groundwater wells, >500 permanent forest inventory plots, and a 150 m tall tower (a combined ecosystem-atmosphere station of the ICOS, Integrated Carbon Observation System) for measurements of atmospheric gas concentrations and biosphere-atmosphere exchanges of carbon, water, and energy. In addition, the KCS has also been the focus of numerous high resolution multi-spectral LiDAR measurements and large scale experiments. This large collection of equipment and data generation supports a range of disciplinary studies, but more importantly fosters multi-, trans-, and interdisciplinary research opportunities. The KCS attracts a broad collection of scientists, including biogeochemists, ecologists, foresters, geologists, hydrologists, limnologists, soil scientists, and social scientists, all of whom bring their knowledge and experience to the site. The combination of long-term monitoring, shorter-term research projects, and large-scale experiments, including manipulations of climate and various forest management practices, has contributed much to our understanding of boreal landscape functioning, while also supporting the development of models and guidelines for research, policy, and management.     

Observations of neutral hydrogen in the corona

Measurements of the linear polarization degree (p) of the emission in the inner corona (<1.5 R _⊙) in the green and red continuum ranges are discussed. “Excess p” (Δp = pgreen — pred), which reaches its maximum at 1.2–1.3 R _⊙ (<5%) in the polar regions and 20% within large-scale streamers, is present everywhere in the corona. It has been mentioned that excess p can be explained by the presence of neutral hydrogen in the corona.     

Underwater Explosion (UWE) Analysis of the ROKS Cheonan Incident

The underwater explosion (UWE) resulting in the sinking of the South Korean warship, ROKS Cheonan occurred on March 26 2010. Raw data was analyzed from several 3-component stations—Baengyeong-do Korea Meteorological Administration (KMA) station (BAR), Ganghwa KMA station (GAHB), Incheon Incorporated Research Institutions for Seismology (IRIS) station (INCN), the short-period station—Deokjeok-do KMA station (DEI), as well as from the seismo-acoustic array Baengyeong-do Korea Institute of Geoscience and Mineral Resources (KIGAM) station (BRDAR). The ROKS Cheonan incident has been investigated by both the Multinational Civilian-Military Joint Investigation Group (Ministry of National Defense, 2010) and Hong (Bull Seism Soc Am 101:1554–1562, 2011). Their respective methods and conclusions are also presented in this study. One of the main differences between their findings and ours is that we deducted that the fundamental bubble frequency was 1.01 Hz with a subsequent oscillation of 1.72 Hz. Also, in contrast to findings by the MCMJIG and Hong, our analysis shows the first reverberation frequency to be 8.5 Hz and the subsequent one to be ≈25 Hz. The TNT-equivalent charge weight (seismic yield) and seismic magnitude were estimated from an observed bubble frequency of 1.01 Hz and the analytical model of a bubble pulse. From the data analyzed, we deducted that the seismic yield would be about 136 kg of TNT, which is equivalent to the individual yield of a large number of land control mines (LCM) which were abandoned in the vicinity of the ROKS Cheonan incident by the Republic of Korea (ROK) Navy in the 1970s (Ministry of National Defense 2010). Also, whereas both the MCMJIG and HONG estimated the local magnitude at 1.5, our findings came to the conclusion of a local magnitude of approximately 2.04 based on the bubble frequency of 1.01 Hz measured on the vertical component of BAR station data considering the empirical relationship between charge weight (TNT yield) and underwater explosion magnitude. Strong high-frequency signals collected at the 3-component BAR station approximately 30 s after P-wave arrivals and infrasound records at BRDAR clearly indicate powerful acoustic phases and N-waves caused by a relatively shallow UWE. T-phases are also observed on seismograms and spectra at 15–17 Hz on the DEI, GAHB, and INCN stations.