A control-volume finite-element method for three-dimensional multiphase basin modeling

In this paper we describe a 3D control-volume finite-element method to solve numerically the coupled partial differential equations (PDEs) governing geological processes involved in the evolution of sedimentary basins. These processes include sediment deposition and deformation, hydrocarbon generation, multiphase fluid flow, and heat transfer in deforming porous media.     

An integrated procedure to evaluate hydrological models

The growing concern for health‐related problems deriving from pollutants leaching is driving national and international administrations to support the development of tools for evaluating the effects of alternate management scenarios and identifying vulnerable areas. Cropping systems models are powerful tools for evaluating leachates under different environmental, social, and management conditions. As percolating water is the transport vehicle for pollutants transport in soil, a reliable evaluation of water balance models is a fundamental prerequisite for investigating pesticides and nitrate fate. As specific approaches for the evaluation of multi‐layer evolution of state variables are missing, we propose a fuzzy‐based, integrated indicator (I_SWC: 0, best; 1, worst) for a comprehensive evaluation of soil water content (SWC) simulations. We aggregated error metrics with others quantifying the homogeneity of errors across different soil layers, the capability of models to reproduce complex dynamics function of both time and soil depth, and model complexity. We tested I_SWC on a sample dataset where the models CropSyst and CERES‐Wheat were used to simulate SWC for winter wheat systems. I_SWC revealed that, in the explored conditions, the global assessment of the two models' performances allowed identification of CropSyst as the best (average I_SWC = 0·441, with a value of 0·537 obtained by CERES‐Wheat), although each model prevailed for some of the metrics. CropSyst presented the highest accuracy (average agreement module = 0·400), whereas CERES‐Wheat's accuracy was slightly worse, although achieved with a simplified modelling approach (average Akaike Information Criterion = − 230·44), thereby favouring large‐area applicability. The non‐univocal scores achieved by the models for the different metrics support the use of multi‐metric evaluation approaches for quantifying the different aspects of water balance model performances. Copyright © 2010 John Wiley & Sons, Ltd.     

Recent Changes Occurred in the Terminus of the Debris-covered Bilafond Glacier in the Karakoram Himalayas Using Remotely Sensed Images and Digital Elevation Models (1978-2011)

Recent changes occurred in terminus of the debris-covered Bilafond Glacier in the Karakoram Range in the Himalayas, Northern Pakistan was investigated in this research. Landsat MSS, TM and ETM＋ images were used for this study. Digital elevation models derived from ASTER GDEM and SRTM were also utilized. Visible, infrared and thermal infrared channels were utilized in order to get accurate glacier change maps. Three methods were tried to map this debris-covered glacier in this research. The glacier has been mapped successfully and the changes in the glacier terminus from 1978 to 2011 have been calculated. Manual, semi-automatic and thermal methods were found to give similar results. It was found that the glacier has undergone serious ablation during this period despite of the fact that many of the larger glaciers in the Hindu Kush and Karakoram mountain regions in the Upper Indus Basin were reported to be expanding. The terminus has been moved back about 600 meters during this period and there was an abrupt change in the glacier terminus during 1990-2002. We propose that debris thickness is not the only factor that influences the glacier ablation but the altitude of the debris-covered glacier as well. Many glaciers in the Karakoram region reported to be expanding were having higher altitudes compared to the study area.     

Combined influence of PDO and ENSO on northern Andean glaciers: a case study on the Cotopaxi ice-covered volcano,Ecuador

This paper describes the application of remote sensing in monitoring the fluctuations in one of the mountain glaciers in the Ecuadorean Andes during the past few decades using ASTER, EO-1 ALI, Landsat MSS, TM and ETM + images. Satellite images were used to calculate the snow line altitudes (SLAs) during the period 1979–2013. Cotopaxi ice covered volcano was studied as representative of Ecuadorian glaciers in the eastern cordillera. Precipitation and air temperature data from various gauging stations within the range of 30 km from the study site and monthly discharge and water level data from a gauging station were also utilized in this study. Anomalies in precipitation and temperature were found to be slightly different in the Cotopaxi region compared to nearby Antizana in the same cordillera and Chimborazo region in the western cordillera. An attempt to correlate the El Niño—southern oscillation phenomenon with the glacier fluctuations in Ecuadorian Andes was done successfully. Cold and warm regimes of Pacific Decadal Oscillation is also considered. The calculated glacier fluctuations obtained were similar to that performed on the nearby Antizana 15 in the eastern cordillera during 1995–2002. Precipitation and temperature anomalies were similar with Antizana 15. It is evident from the research that SLAs were highly fluctuated between the period of occurrence of El Niño and La Niña events. It is also seen that the glacier fluctuations show a negative mass balance trend in during the warm regime of Pacific Decadal Oscillation during the past three decades. Glaciated areas were advanced during the La Nina events in the cold regime of PDO during 1998–2002.     

Sensitivity Analysis and Investigation of the Behaviour of the UTOPIA Land-Surface Process Model: A Case Study for Vineyards in Northern Italy

We used sensitivity-analysis techniques to investigate the behaviour of the land-surface model UTOPIA while simulating the micrometeorology of a typical northern Italy vineyard (Vitis vinifera L.) under average climatic conditions. Sensitivity-analysis experiments were performed by sampling the vegetation parameter hyperspace using the Morris method and quantifying the parameter relevance across a wide range of soil conditions. This method was used since it proved its suitability for models with high computational time or with a large number of parameters, in a variety of studies performed on different types of biophysical models. The impact of input variability was estimated on reference model variables selected among energy (e.g. net radiation, sensible and latent heat fluxes) and hydrological (e.g. soil moisture, surface runoff, drainage) budget components. Maximum vegetation cover and maximum leaf area index were ranked as the most relevant parameters, with sensitivity indices exceeding the remaining parameters by about one order of magnitude. Soil variability had a high impact on the relevance of most of the vegetation parameters: coefficients of variation calculated on the sensitivity indices estimated for the different soils often exceeded 100 %. The only exceptions were represented by maximum vegetation cover and maximum leaf area index, which showed a low variability in sensitivity indices while changing soil type, and confirmed their key role in affecting model results.     

Monitoring and modelling of hydrological processes in the semiarid region of Brazil: The Cariri experimental basins

Catchments have highly variable yields of runoff and soil erosion. The size, land use and the surface cover play a significant role and influence the catchment response and parameter values of simulation models. Two experimental basins—the Cariri basins—were equipped in a semi-arid region of Brazil, for obtaining runoff and sediment yield at different catchment scales, as well as, to evaluate the influence of the land use and surface cover. In the first basin, located in the municipality of Sumé, the field studies were carried out at two different scales: four micro-catchments with an area of around 0.5 ha and nine standard Wischmeier-type erosion plots of 100 m². The experimental units had varied vegetation and management. They were subjected only to natural rainfall events, and were monitored from 1982 to 1991. The total runoff and total sediment yield were determined for each of the events. The installations in the second basin, in the municipality of São João do Cariri, from 1999, include two erosion plots, three micro-catchments, and two sub-catchments of a small basin. These basins are still being monitored for runoff and sediment production. Among the micro-catchments two are nested to detect any scale effect at the micro-catchment level. Nearly 600 events of precipitation, that produced runoff in at least one of the experimental units, have been registered. These data have been used to evaluate the influence of various factors, including cultivation practices and to calibrate hydrological models for plots and micro-catchments. Parameters have been tested by means of cross validations among micro-catchments and sub-catchments. The data sets are made available to all the catchment hydrology researchers and others at https://doi.org/10.5281/zenodo.4690886 .     

Midnight reversal of ionospheric plasma bubble eastward velocity to westward velocity during geomagnetically quiettime: Climatology and its model validation

In an effort to better understand the dynamics of westward velocities of the nocturnal F-region plasma, the climatology of the westward traveling plasma bubbles – WTB – occurring during quiettime is studied here. The climatology of the WTB is analyzed here based on airglow images obtained during 14 quiet days between 2001 and 2006 at the Brazilian station São João do Cariri (Geographic 7.45°S, 36.5°W, dip ～20°S). The frequency of occurrence of the WTB maximizes in the descending phase of the solar cycle. The WTB velocities ranged between ～20 and 40 ms^?1. The frequency of occurrence had a peak value of only 3.65% at 2345 LT. The maximum occurrence of the WTB was in July–September. No WTB have been observed from November until April in all years 2001–2006. We show for the first time theoretically that the WTB dominant forcing mechanisms during geomagnetically quiet days are westward thermospheric winds.     

Recent trends in annual snowline variations in the northern wet outer tropics: case studies from southern Cordillera Blanca,Peru

TL-moments approach has been used in an analysis to identify the best-fitting distributions to represent the annual series of maximum streamflow data over seven stations in Johor, Malaysia. The TL-moments with different trimming values are used to estimate the parameter of the selected distributions namely: Three-parameter lognormal (LN3) and Pearson Type III (P3) distribution. The main objective of this study is to derive the TL-moments (t ₁,0), t ₁?=?1,2,3,4 methods for LN3 and P3 distributions. The performance of TL-moments (t ₁,0), t ₁?=?1,2,3,4 was compared with L-moments through Monte Carlo simulation and streamflow data over a station in Johor, Malaysia. The absolute error is used to test the influence of TL-moments methods on estimated probability distribution functions. From the cases in this study, the results show that TL-moments with four trimmed smallest values from the conceptual sample (TL-moments [4, 0]) of LN3 distribution was the most appropriate in most of the stations of the annual maximum streamflow series in Johor, Malaysia.     

Recent variations of supraglacial lakes on the Baltoro Glacier in the central Karakoram Himalaya and its possible teleconnections with the pacific decadal oscillation

This study discusses the formation and variations of supraglacial lakes on the Baltoro glacier system in the central Karakoram Himalaya during the last four decades. We mapped supraglacial lakes on the Baltoro Glacier from 1978 to 2014 using Landsat MSS, TM, ETM + and LCDM images. Most of the glacial lakes were formed or expanded during the late 1970s–2008. After 2008, the total number and the area of glacial lakes were found to be lesser compared to previous years. We tried to find any teleconnections exists between the glacial lake changes in this region and the pacific decadal oscillation (PDO), which entered its prolonged warm regime in the late 1970s and again to its cold regime in 2008, based on observational investigation. The decrease in the number and area of the supraglacial lakes after 2008 is hypothesised to be linked with the recent cold phase of PDO.     

An integrated evaluation of thirteen modelling solutions for the generation of hourly values of air relative humidity

The availability of hourly air relative humidity (HARH) data is a key requirement for the estimation of epidemic dynamics of plant fungal pathogens, in particular for the simulation of both the germination of the spores and the infection process. Most of the existing epidemic forecasting models require these data as input directly or indirectly, in the latter case for the estimation of leaf wetness duration. In many cases, HARH must be generated because it is not available in historical series and when there is the need to simulate epidemics either on a wide scale or with different climate scenarios. Thirteen modelling solutions (MS) for the generation of this variable were evaluated, with different input requirements and alternative approaches, on a large dataset including several sites and years. A composite indicator was developed using fuzzy logic to compare and to evaluate the performances of the models. The indicator consists of four modules: Accuracy, Correlation, Pattern and Robustness. Results showed that when available, daily maximum and minimum air relative humidity data substantially improved the estimation of HARH. When such data are not available, the choice of the MS is crucial, given the difference in predicting skills obtained during the analysis, which allowed a clear detection of the best performing MS. This study represents the first step of the creation of a robust modelling chain coupling the MS for the generation of HARH and disease forecasting models, including the systematic validation of each step of the simulation.