Thermo-hydro-mechanical modeling of unsaturated soils using isogeometric analysis: Model development and application to strain localization simulation

This study presents a thermo-hydro-mechanical (THM) model of unsaturated soils using isogeometric analysis (IGA). The framework employs Bézier extraction to connect IGA to the conventional finite element analysis (FEA), featuring the current study as one of the first attempts to develop an IGA-FEA framework for solving THM problems in unsaturated soils. IGA offers higher levels of interelement continuity making it an attractive method for solving highly nonlinear problems. The governing equations of linear momentum, mass, and energy balance are coupled based on the averaging procedure within the hybrid mixture theory. The Drucker-Prager yield surface is used to limit the modified effective stress where the model follows small strain, quasi-static loading conditions. Temperature dependency of the surface tension is implemented in the soil-water retention curve. Nonuniform rational B-splines (NURBS) basis functions are used in the standard Galerkin method and weak formulations of the balance equations. Displacement, capillary pressure, gas pressure, and temperature are four independent quantities that are approximated by NURBS in spatial discretization. The framework is used to simulate strain localization in an undrained dense sand subjected to plane strain biaxial compression under different temperatures and displacement velocities. Results show that an increase in the displacement rate leads to reduction in the equivalent plastic strain while an increase in the temperature leads to an increase in the equivalent plastic strain. The findings suggest that the proposed IGA-based framework offers a viable alternative for solving THM problems in unsaturated soils.     

The first peak ground motion attenuation relationships for North of Vietnam

The first attenuation relationships of peak ground acceleration (PGA) and peak ground velocity (PGV) for northern Vietnam are obtained in this study. Ground motion data are collected by a portable broadband seismic network in northern Vietnam as a part of cooperation between the Institute of Geophysics, Vietnamese Academy of Science and Technology, Vietnam and Institute of Earth Sciences, Academia Sinica, Taiwan. The database comprises a total of 330 amplitude records by 14 broadband stations from 53 shallow earthquakes, which were occurred in and around northern Vietnam in the period between 01/2006 and 12/2009. These earthquakes are of local magnitudes between 1.6 and 4.6, focal depths less than 30 km, and epicentral distances less than 500 km. The new attenuation relationships for PGA and PGV are:

log₁₀(PGA)=-0.987+0.7521M_L-log₁₀(R)-0.00475R,     

Projected 21st century climate change on snow conditions over Shasta Dam watershed by means of dynamical downscaling

Snow is an important component of the Earth's climate system and is particularly vulnerable to global warming. It has been suggested that warmer temperatures may cause significant declines in snow water content and snow cover duration. In this study, snowfall and snowmelt were projected by means of a regional climate model that was coupled to a physically based snow model over Shasta Dam watershed to assess changes in snow water content and snow cover duration during the 21st century. This physically based snow model requires both physical data and future climate projections. These physical data include topography, soils, vegetation, and land use/land cover, which were collected from associated organizations. The future climate projections were dynamically downscaled by means of the regional climate model under 4 emission scenarios simulated by 2 general circulation models (fifth‐generation of the ECHAM general circulation model and the third‐generation atmospheric general circulation model). The downscaled future projections were bias corrected before projecting snowfall and snowmelt processes over Shasta Dam watershed during 2010–2099. This study's results agree with those of previous studies that projected snow water equivalent is decreasing by 50–80% whereas the fraction of precipitation falling as snowfall is decreasing by 15% to 20%. The obtained projection results show that future snow water content will change in both time and space. Furthermore, the results confirm that physical data such as topography, land cover, and atmospheric–hydrologic data are instrumental in the studies on the impact of climate change on the water resources of a region.     

Physically based maximum precipitation estimation under future climate change conditions

Estimation of the extreme precipitation over a target watershed under a changing climate would be necessary to design safe large hydraulic structures. For this purpose, the maximum precipitation (MP) estimation approach was applied to the American River Watershed (ARW) in Northern California under several future climate conditions over 90 water years (2010–2099). These future climate conditions were obtained using 13 future climate projections from two general circulation models (ECHAM5 and CCSM3) based on four future climate scenarios (Special Report on Emissions Scenarios: A1B, A1FI, A2, and B1). A total of 1,170 future projected severe storm events (90 years × 13 projections) were selected with respect to the 72‐hr basin‐average precipitation over the ARW. The 72‐hr basin‐average precipitation for each of the selected severe storm events was maximized over the ARW by horizontally shifting the atmospheric boundary conditions of a regional atmospheric model in order to optimize the path of the storm system that corresponded to the particular event. After maximization, the MP estimates, which are the largest precipitation depths among the maximized results, were obtained as 836.7 mm for the early half‐century period (2010–2054) and 1,056.5 mm for the late half‐century period (2055–2099).     

A new species Lobophora tsengii sp. nov. (Dictyotales; Phaeophyceae) from Bach Long Vy (Bailongwei) Island,Vietnam

Journal of Oceanology and Limnology - A new species Lobophora tsengii is described from Bach Long Vy (Bailongwei) Island, Haiphong Province, Vietnam. The plants inhabit the subtidal zone and have...     

3-D Shear Wave Velocity Structure of the Crust and Upper Mantle in South China Sea and its Surrounding Regions by Surface Wave Dispersion Analysis

In this study, we construct a 3-D shear wave velocity structure of the crust and upper mantle in South China Sea and its surrounding regions by surface wave dispersion analysis. We use the multiple filter technique to calculate the group velocity dispersion curves of fundamental mode Rayleigh and Love waves with periods from 14 s to 120 s for earthquakes occurred around the Southeast Asia. We divide the study region (80° E–140° E, 16° S–32° N) into 3° × 3° blocks and use the constrained block inversion method to get the regionalized dispersion curve for each block. At some chosen periods, we put together laterally the regionalized group velocities from different blocks at the same period to get group velocity image maps. These maps show that there is significant heterogeneity in the group velocity of the study region. The dispersion curve of each block was then processed by surface wave inversion method to obtain the shear wave velocity structure. Finally, we put the shear wave velocity structures of all the blocks together to obtain the three-dimensional shear wave velocity structure of crust and upper mantle. The three-dimensional shear wave velocity structure shows that the shear wave velocity distribution in the crust and upper mantle of the South China Sea and its surrounding regions displays significant heterogeneity. There are significant differences among the crustal thickness, the lithospheric thickness and the shear wave velocity of the lid in upper mantle of different structure units. This study shows that the South China Sea Basin, southeast Sulu Sea Basin and Celebes Sea Basin have thinner crust. The thickness of crust in South China Sea Basin is 5–10 km; in Indochina is 25–40 km; in Peninsular Malaysia is 30–35 km; in Borneo is 30–35 km; in Palawan is 35 km; in the Philippine Islands is 30–35 km, in Sunda Shelf is 30–35 km, in Southeast China is 30–40 km, in West Philippine Basin is 5–10 km. The South China Sea Basin has a lithosphere with thickness of about 45–50 km, and the shear wave velocity of its lid is about 4.3–4.7 km/s; Indochina has a lithosphere with thickness of about 55–70 km, and the shear wave velocity of its lid is about 4.3–4.5 km/s; Borneo has a lithosphere with thickness of about 55–60 km, and the shear wave velocity of its lid is about 4.1–4.3 km/s; the Philippine Islands has a lithosphere with thickness of about 55–60 km, and the shear wave velocity of its lid is about 4.2–4.3 km/s, West Philippine Basin has a lithosphere with thickness of about 50–55 km, and the shear wave velocity of its lid is about 4.7–4.8 km/s, Sunda Self has a lithosphere with thickness of about 55–65 km, and the shear wave velocity of its lid is about 4.3 km/s. The Red-River Fault Zone probably penetrates to a depth of at least 200 km and is plausibly the boundary between the South China Block and the Indosinia Block.     

Lumped surface and sub‐surface runoff for erosion modeling within a small hilly watershed in northern Vietnam

Developing models to predict on‐site soil erosion and off‐site sediment transport at the agricultural watershed scale represent an on‐going challenge in research today. This study attempts to simulate the daily discharge and sediment loss using a distributed model that combines surface and sub‐surface runoffs in a small hilly watershed (< 1 km²). The semi‐quantitative model, Predict and Localize Erosion and Runoff (PLER), integrates the Manning–Strickler equation to simulate runoff and the Griffith University Erosion System Template equation to simulate soil detachment, sediment storage and soil loss based on a map resolution of 30 m × 30 m and over a daily time interval. By using a basic input data set and only two calibration coefficients based, respectively, on water velocity and soil detachment, the PLER model is easily applicable to different agricultural scenarios. The results indicate appropriate model performance and a high correlation between measured and predicted data with both Nash–Sutcliffe efficiency (Ef) and correlation coefficient (r²) having values > 0.9. With the simple input data needs, PLER model is a useful tool for daily runoff and soil erosion modeling in small hilly watersheds in humid tropical areas. Copyright © 2013 John Wiley & Sons, Ltd.     

Using Geomorphological Indicators to Predict Earthquake Magnitude (MOb‒Max): A Case Study from Cao Bang Province and Adjasent Areas (Vietnam)

Geotectonics - This paper presents analysis of geomorphological indices for predicting the maximum observed earthquake (MOb?Max) in research region, conducted in three steps: (i)...     

Relating Radar Remote Sensing of Biomass to Modelling of Forest Carbon Budgets

This paper addresses the use of radar remote sensing to map forest above-ground biomass, and discusses the use of biomass maps to test a dynamic vegetation model that identifies carbon sources and sinks and predicts their variation over time. For current radar satellite data, only the biomass of young/sparse forests or regrowth after disturbances can be recovered. An example from central Siberia illustrates that biomass can be measured by radar at a continental scale, and that a significant proportion of the Siberian forests have biomass values less than 50 tonnes/ha. Comparison between the radar map and calculations by the Sheffield Dynamic Global Vegetation Model (SDGVM) indicates that the model considerably overestimates biomass; under-representation of managed areas, disturbed areas and areas of low site quality in the model are suggested reasons for this effect. A case study carried out at the Büdingen plantation forest in Germany supports the argument that inadequate representations of site quality and forest management may cause model overestimates of biomass. Comparison of the calculated biomass of stands planted after 1990 with biomass estimates by radar allows identification of forest stands where the growth conditions assumed by the model are not valid. This allows a quality check on model calculations of carbon fluxes: only calculations for stands where there is good agreement between the data and the model predictions should be accepted. Although the paper only uses the SDGVM model, similar effects are likely in other dynamic vegetation models, and the results show that model calculations attempting to quantify the role of forests as carbon sources or sinks could be qualified and potentially improved by exploiting remotely sensed measurements of biomass.     

Simultaneous Organic Matter and Nitrogen Removal from Piggery Slaughterhouse Wastewater Using Bioaugmented Gaslift Membrane Bioreactor

A single bioaugmentation reactor and a side-stream gaslift membrane bioreactor combined with bioaugmentation are conducted to treat real wastewater from a centralized piggery slaughterhouse in Vietnam. The bioaugmented reactor is inoculated with heterotrophic microorganisms (Bacillus sp.) isolated from piggery slaughterhouse wastewater. The results of a single bioaugmentation experiment show high removal efficiency of chemical oxygen demand (COD) (84.8%–97.5%) and total nitrogen (TN) (69.9%–87.2%) at loading rates of 1.28–3.89 and 0.14–0.37 kg m⁻³ d⁻¹, respectively. The combined system demonstrates a significantly higher TN removal efficiency (89.0%–96.1%) (p < 0.001), more stable flux (36.0–38.4 L m⁻² h⁻¹), and transmembrane pressure (0.95–1.05 bar), and better capacity of separation of solid–liquid phases compared to the single bioaugmentation. High COD and TN removal efficiency is possibly due to assimilation and simultaneous nitrification and denitrification processes. The results of this study also indicate the feasibility and propitious efficiency of the bioaugmented gaslift membrane bioreactor for piggery slaughterhouse wastewater treatment.