FEM × DEM modelling of cohesive granular materials: Numerical homogenisation and multi-scale simulations

The article presents a multi-scale modelling approach of cohesive granular materials, its numerical implementation and its results. At microscopic level, Discrete Element Method (DEM) is used to model dense grains packing. At the macroscopic level, the numerical solution is obtained by a Finite Element Method (FEM). In order to bridge the micro- and macro-scales, the concept of Representative Elementary Volume (REV) is applied, in which the average REV stress and the consistent tangent operators are obtained in each macroscopic integration point as the results of DEM’s simulation. In this way, the numerical constitutive law is determined through the detailed modelling of the microstructure, taking into account the nature of granular materials. We first elaborate the principle of the computation homogenisation (FEM × DEM), then demonstrate the features of our multiscale computation in terms of a biaxial compression test. Macroscopic strain location is observed and discussed.     

Are satellite products good proxies for gauge precipitation over Singapore?

The uncertainties in two high-resolution satellite precipitation products (TRMM 3B42 v7.0 and GSMaP v5.222) were investigated by comparing them against rain gauge observations over Singapore on sub-daily scales. The satellite-borne precipitation products are assessed in terms of seasonal, monthly and daily variations, the diurnal cycle, and extreme precipitation over a 10-year period (2000–2010). Results indicate that the uncertainties in extreme precipitation is higher in GSMaP than in TRMM, possibly due to the issues such as satellite merging algorithm, the finer spatio-temporal scale of high intensity precipitation, and the swath time of satellite. Such discrepancies between satellite-borne and gauge-based precipitations at sub-daily scale can possibly lead to distorting analysis of precipitation characteristics and/or application model results. Overall, both satellite products are unable to capture the observed extremes and provide a good agreement with observations only at coarse time scales. Also, the satellite products agree well on the late afternoon maximum and heavier rainfall of gauge-based data in winter season when the Intertropical Convergence Zone (ITCZ) is located over Singapore. However, they do not reproduce the gauge-observed diurnal cycle in summer. The disagreement in summer could be attributed to the dominant satellite overpass time (about 14:00 SGT) later than the diurnal peak time (about 09:00 SGT) of gauge precipitation. From the analyses of extreme precipitation indices, it is inferred that both satellite datasets tend to overestimate the light rain and frequency but underestimate high intensity precipitation and the length of dry spells. This study on quantification of their uncertainty is useful in many aspects especially that these satellite products stand scrutiny over places where there are no good ground data to be compared against. This has serious implications on climate studies as in model evaluations and in particular, climate model simulated future projections, when information on precipitation extremes need to be reliable as they are highly crucial for adaptation and mitigation.     

Black carbon decomposition under varying water regimes

The stability of biomass-derived black carbon (BC) or biochar as a slow cycling pool in the global C cycle is an important property and is likely governed by environmental conditions. This study investigated the effects of water regimes (saturated, unsaturated and alternating saturated–unsaturated conditions) and differences in BC materials, produced by carbonizing corn residues and oak wood at two temperatures (350 °C and 600 °C) on BC degradation at 30 °C over 1 year in a full factorial experiment. Effects of water regime on C loss and potential cation exchange capacity (CECp at pH 7) significantly depended on biomass type. Corn BC was both mineralized (16% C loss for the first year) and was oxidized [1000 mmole(+) kg^?1 C] significantly faster under unsaturated conditions than under other water regimes, whereas oak BC mineralized most rapidly (12%) under alternating saturated–unsaturated conditions with similar oxidation, irrespective of water regime. Over 1 year of saturated incubation, the O/C ratio values did not significantly (P > 0.05) increase even though BC was mineralized by 9% and CECp increased by 170 mmole(+) kg^?1 C, in contrast to unsaturated and alternating saturated–unsaturated conditions. While mineralization and oxidation significantly decreased at higher charring temperature for corn, no difference was observed for oak (P > 0.05). Unsaturated and alternating conditions increased carboxylic and OH functional groups, while they decreased aliphatic groups. The pH increased by about one unit for corn BC, but decreased by 0.2 units for oak BC, indicating strong mineral dissolution of corn BC. Carbon loss strongly correlated with changes in O/C values of both corn BC and oak BC, indicating that oxidation of BC was most likely the major mechanism controlling its stability. However, under saturated conditions, additional mechanisms may govern BC degradation and require further investigation.     

Topography of the Moho and earth crust structure beneath the East Vietnam Sea from 3D inversion of gravity field data

The Moho depth, crustal thickness and fault systems of the East Vietnam Sea (EVS) are determined by 3D interpretation of satellite gravity. The Moho depth is calculated by 3D Parker inversion from residual gravity anomaly that is obtained by removing the gravity effects of seafloor and Pre-Cenozoic sediment basement topographies from the free air anomaly. The 3D inversion solution is constrained by power density spectrum of gravity anomaly and seismic data. The calculated Moho depths in the EVS vary from 30–31 km near the coast to 9 km in the Central Basin. A map of the lithosphere extension factor in the Cenozoic is constructed from Moho and Pre-Cenozoic sediment basement depths. The fault systems constructed by the maximum horizontal gradient approach include NE-SW, NW-SE, and N-S oriented faults. Based on the interpretation results, the EVS is sub-divided into five structural zones which demonstrated the different characteristics of the crustal structure.     

Influence of anthropogenic pressure and seasonal upwelling on coral reefs in Nha Trang Bay (Central Vietnam)

The growing coastal development, dredging and dumping activities, overfishing and expansion of marine cage culture in Nha Trang Bay (NTB) of Central Vietnam since the beginning of the 2000s have resulted in a dramatic decrease of live coral cover. Surveys conducted in April–May 2013 and the same period in 2014 revealed that with an increase in distance from the outer part of the bay towards the mainland, the rivers’ influxes and dredged areas, coral cover decreased from 75% to 0.6% and species richness from 63 to 5, while the abundance of macroalgae increased from 0% to 56%. These changes correlate with differences in the concentration of suspended sediments on the same gradient. The abundance of the crown‐of‐thorns starfish Acanthaster planci and of the echinoid Diadema setosum significantly increased between the first estimation in 1998 and the survey in 2014, from 0 to 1.7 individuals (ind.) per 100 m^?2 and from 50.8 to 94.5 ind. per 100 m^?2, respectively, contributing to coral loss and intensive bioerosion of the reef framework in the bay. The large sizes of adult colonies of tabulate Acropora on the remote stations with negligible sedimentation and eutrophication loads were inconsistent with the assumptions that temperature‐induced coral bleaching or cyclones could be the major impacts in Nha Trang Bay. Analysis of the 16‐year thermal history of the bay did not reveal any instances in which the coral thermal bleaching threshold had been exceeded up to the present study. Seasonal upwelling, which occurs annually in the vicinity of Central Vietnam, may contribute to mitigation of thermal anomalies within NTB and to the maintenance of healthy coral communities on the remote reefs with relatively low anthropogenic impact.     

Beach profiles characteristics along Giao Thuy and Hai Hau coasts, Vietnam: A field study

Giao Thuy and Hai Hau coasts are located in Nam Dinh province, Vietnam, with a total coastline of 54.42 km in length. The sea-dike system has been seriously damaged and there have been many dike breaches which caused floods and losses. This situation is considered of a general representative for coastal area in the northern part of Vietnam. A variety of studies have shown that the gradient in the longshore sediment transport rate and the offshore ?ne sediment lost are the main mechanisms causing the beach erosion. This study presents a field investigation of the beach profiles at Giao Thuy and Hai Hau beaches. Three types of empirical functions for the equilibrium beach profile are applied and compared with the observations. Results show that all observed beach pro?les can be described by a single function. However, one specific equilibrium pro?le equation is not sufficient to assess all beach pro?les. In Section 1 of Giao Thuy and Section 3 of Hai Thinh beaches, beach profiles are consistent with the logarithmic function, while the exponential function fits well in Section 2. This difference is explained with respect to coastal morphology, sediment characteristics and hydrodynamic conditions which vary in site. An analysis of the validity of the beach profile functions is recommended for the numerical modeling and engineering designs in this area.     

EM and GPR investigations of contaminant spread around the Hoc Mon waste site,Vietnam

This paper is presenting the results from near-surface geophysical surveys near the waste site of Hoc Mon in southern Vietnam where leachate contamination has been recognized at the surface. Using EM and GPR surveys, we were able to determine the lateral extent of a contaminated area of high electrical conductivity and have identified channels that concentrate the contaminant flow. The simple relationship between the electrical resistivity and the leachate concentration is suggested and estimated the in situ leachate concentration from the inversion of the EM data; values are as high as 40%. Thanks to a permeability barrier leachate flow is confined to the shallow subsurface, making it easier to apply possible site remediation projects.     

Geophysical techniques to aquifer locating and monitoring for industrial zones in North Hanoi, Vietnam

Geophysical methods were applied for hydrogeological targets in many countries including Vietnam. This paper presents results of using complex geophysical techniques as well as 2D electrical resistivity imaging (ERI), vertical electrical sounding (VES), very low frequency (VLF), and seismic refraction for geological structure investigation for locating the aquifers and assessing the hydrogeological conditions for groundwater potential in industrial zones of North Hanoi, Vietnam. The locations of two aquifers are determined by their depth and thickness on the basis of resistivity and seismic velocity values which were proved by stratifications of three boreholes to 40–60 m of depth on the study area. There are connections from surface water to shallow aquifer by hydraulic windows, as follows from VLF data. The deeper aquifer can be considered as a potential groundwater supply, but the water level is descending in time, as shown by hydrological monitoring. However, with careful use and by reducing sources of pollution, groundwater can continue to be an important natural resource for future.     

Catchment storm runoff modelling using the geomorphologic instantaneous unit hydrograph

A lumped empirical model, the Geomorphologic Instantaneous Unit Hydrograph (GIUH) rainfall-runoff model, is developed for the Can Le catchment in the upstream region of the Sai Gon river Basin (Vietnam). This model can serve to simulate catchment runoff into the Dau Tieng Reservoir and can be used as a flood forecasting tool for the ungauged Can Le catchment. The GIUH couples geomorphology and hydrology quantitatively. The obtained Unit Hydrograph is based on Horton's morphometric parameters; bifurcation, length and area ratios. A new functionality within the ILWIS GIS-RS package, namely ‘DEM-hydro processing’, is applied to effectively process a Digital Elevation Model to extract these ratios from the drainage network. To supplement the limited field data available, various satellites images have been used such as ASTER, SRTM (Shuttle Radar Topography Mission) and METEOSAT 5. A short field campaign to collect missing ground data was executed between September and October 2005. The data collected included discharge (and stage – discharge curve), meteorological data, soil, land use information that are used for paramerisation, calibration and validation of the GIUH. The model was successfully applied for the Can Le catchment. Using the Horton's morphometric parameters derived from the DEM with estimated overland and stream flow velocities, the model is easy-to-use.