Ioannis Tsoukalas

Civil Engineer, MSc., PhD candidate
jtsoukalas@hotmail.com

Participation in research projects

Participation as Researcher

  1. Maintenance, upgrading and extension of the Decision Support System for the management of the Athens water resource system

Published work

Publications in scientific journals

  1. I. Tsoukalas, C. Makropoulos, and D. Koutsoyiannis, Simulation of stochastic processes exhibiting any-range dependence and arbitrary marginal distributions, Water Resources Research, 2018, (in review).
  2. I. Tsoukalas, S.M. Papalexiou, A. Efstratiadis, and C. Makropoulos, A cautionary note on the reproduction of dependencies through linear stochastic models with non-Gaussian white noise, Water, 10 (6), 771, doi:10.3390/w10060771, 2018.
  3. G. Papaioannou, A. Efstratiadis, L. Vasiliades, A. Loukas, S.M. Papalexiou, A. Koukouvinos, I. Tsoukalas, and P. Kossieris, An operational method for Floods Directive implementation in ungauged urban areas, Hydrology, 5 (2), 24, doi:10.3390/hydrology5020024, 2018.
  4. I. Tsoukalas, A. Efstratiadis, and C. Makropoulos, Stochastic periodic autoregressive to anything (SPARTA): Modelling and simulation of cyclostationary processes with arbitrary marginal distributions, Water Resources Research, 54 (1), 161–185, WRCR23047, doi:10.1002/2017WR021394, 2018.
  5. A. Tegos, N. Malamos, A. Efstratiadis, I. Tsoukalas, A. Karanasios, and D. Koutsoyiannis, Parametric modelling of potential evapotranspiration: a global survey, Water, 9 (10), 795, doi:10.3390/w9100795, 2017.
  6. C. Makropoulos, E. Rozos, I. Tsoukalas, and et al., Sewer-mining: A water reuse option supporting circular economy, public service provision and entrepreneurship, Journal of Environmental Management, doi:10.1016/j.jenvman.2017.07.026, 2017, (in press).
  7. E. Rozos, I. Tsoukalas, K. Ripis, E. Smeti, and C. Makropoulos, Turning black into green: Ecosystem services from treated wastewater, Desalination and Water Treatment, 91 (2017), 2017.
  8. I. Tsoukalas, P. Kossieris, A. Efstratiadis, and C. Makropoulos, Surrogate-enhanced evolutionary annealing simplex algorithm for effective and efficient optimization of water resources problems on a budget, Environmental Modelling and Software, 77, 122–142, doi:10.1016/j.envsoft.2015.12.008, 2016.
  9. I. Tsoukalas, and C. Makropoulos, A surrogate based optimization approach for the development of uncertainty-aware reservoir operational rules: the case of Nestos hydrosystem, Water Resources Management, 29 (13), 4719–4734, doi:10.1007/s11269-015-1086-8, 2015.
  10. I. Tsoukalas, and C. Makropoulos, Multiobjective optimisation on a budget: Exploring surrogate modelling for robust multi-reservoir rules generation under hydrological uncertainty, Environmental Modelling and Software, 69, 396–413, doi:10.1016/j.envsoft.2014.09.023, 2015.

Book chapters and fully evaluated conference publications

  1. I. Tsoukalas, C. Makropoulos, and A. Efstratiadis, Stochastic simulation of periodic processes with arbitrary marginal distributions, 15th International Conference on Environmental Science and Technology (CEST2017), Rhodes, Global Network on Environmental Science and Technology, 2017.
  2. E. Rozos, I. Tsoukalas, K. Ripis, E. Smeti, and C. Makropoulos, Turning black into green: ecosystem services from treated wastewater, 13th IWA Specialized Conference on Small Water and Wastewater Systems, Athens, Greece, National Technical University of Athens, 2016, (in press).
  3. I. Tsoukalas, P. Dimas, and C. Makropoulos, Hydrosystem optimization on a budget: Investigating the potential of surrogate based optimization techniques, 14th International Conference on Environmental Science and Technology (CEST2015), Global Network on Environmental Science and Technology, University of the Aegean, 2015.

Conference publications and presentations with evaluation of abstract

  1. I. Tsoukalas, P. Kossieris, A. Efstratiadis, C. Makropoulos, and D. Koutsoyiannis, CastaliaR: An R package for multivariate stochastic simulation at multiple temporal scales, European Geosciences Union General Assembly 2018, Geophysical Research Abstracts, Vol. 20, Vienna, EGU2018-18433, doi:10.13140/RG.2.2.20978.81605, European Geosciences Union, 2018.
  2. Y. Moustakis, P. Kossieris, I. Tsoukalas, and A. Efstratiadis, Quasi-continuous stochastic simulation framework for flood modelling, European Geosciences Union General Assembly 2017, Geophysical Research Abstracts, Vol. 19, Vienna, 19, EGU2017-534, European Geosciences Union, 2017.
  3. P. Kossieris, A. Efstratiadis, I. Tsoukalas, and D. Koutsoyiannis, Assessing the performance of Bartlett-Lewis model on the simulation of Athens rainfall, European Geosciences Union General Assembly 2015, Geophysical Research Abstracts, Vol. 17, Vienna, EGU2015-8983, doi:10.13140/RG.2.2.14371.25120, European Geosciences Union, 2015.
  4. A. Efstratiadis, I. Tsoukalas, P. Kossieris, G. Karavokiros, A. Christofides, A. Siskos, N. Mamassis, and D. Koutsoyiannis, Computational issues in complex water-energy optimization problems: Time scales, parameterizations, objectives and algorithms, European Geosciences Union General Assembly 2015, Geophysical Research Abstracts, Vol. 17, Vienna, EGU2015-5121, doi:10.13140/RG.2.2.11015.80802, European Geosciences Union, 2015.
  5. A. Drosou, P. Dimitriadis, A. Lykou, P. Kossieris, I. Tsoukalas, A. Efstratiadis, and N. Mamassis, Assessing and optimising flood control options along the Arachthos river floodplain (Epirus, Greece), European Geosciences Union General Assembly 2015, Geophysical Research Abstracts, Vol. 17, Vienna, EGU2015-9148, European Geosciences Union, 2015.
  6. I. Tsoukalas, P. Kossieris, A. Efstratiadis, and C. Makropoulos, Handling time-expensive global optimization problems through the surrogate-enhanced evolutionary annealing-simplex algorithm, European Geosciences Union General Assembly 2015, Geophysical Research Abstracts, Vol. 17, Vienna, EGU2015-5923, European Geosciences Union, 2015.

Details on research projects

Participation as Researcher

  1. Maintenance, upgrading and extension of the Decision Support System for the management of the Athens water resource system

    Duration: October 2008–November 2011

    Budget: €72 000

    Project director: N. Mamassis

    Principal investigator: D. Koutsoyiannis

    This research project includes the maintenance, upgrading and extension of the Decision Support System that developed by NTUA for EYDAP in the framework of the research project “Updating of the supervision and management of the water resources’ system for the water supply of the Athens’ metropolitan area”. The project is consisted of the following parts: (a) Upgrading of the Data Base, (b)Upgrading and extension of hydrometeorological network, (c) upgrading of the hydrometeorological data process software, (d) upgrading and extension of the Hydronomeas software, (e) hydrological data analysis and (f) support to the preparation of the annual master plans

Published work in detail

Publications in scientific journals

  1. I. Tsoukalas, C. Makropoulos, and D. Koutsoyiannis, Simulation of stochastic processes exhibiting any-range dependence and arbitrary marginal distributions, Water Resources Research, 2018, (in review).

    Additional material:

  1. I. Tsoukalas, S.M. Papalexiou, A. Efstratiadis, and C. Makropoulos, A cautionary note on the reproduction of dependencies through linear stochastic models with non-Gaussian white noise, Water, 10 (6), 771, doi:10.3390/w10060771, 2018.

    Since the prime days of stochastic hydrology back in 1960s, autoregressive (AR) and moving average (MA) models (as well as their extensions) have been widely used to simulate hydrometeorological processes. Initially, AR(1) or Markovian models with Gaussian noise prevailed due to their conceptual and mathematical simplicity. However, the ubiquitous skewed behavior of most hydrometeorological processes, particularly at fine time scales, necessitated the generation of synthetic time series to also reproduce higher-order moments. In this respect, the former schemes were enhanced to preserve skewness through the use of non-Gaussian white noise— a modification attributed to Thomas and Fiering (TF). Although preserving higher-order moments to approximate a distribution is a limited and potentially risky solution, the TF approach has become a common choice in operational practice. In this study, almost half a century after its introduction, we reveal an important flaw that spans over all popular linear stochastic models that employ non-Gaussian white noise. Focusing on the Markovian case, we prove mathematically that this generating scheme provides bounded dependence patterns, which are both unrealistic and inconsistent with the observed data. This so-called “envelope behavior” is amplified as the skewness and correlation increases, as demonstrated on the basis of real-world and hypothetical simulation examples.

    Full text: http://www.itia.ntua.gr/en/getfile/1848/1/documents/water-10-00771.pdf (14101 KB)

    See also: http://www.mdpi.com/2073-4441/10/6/771

    Other works that reference this work (this list might be obsolete):

    1. Papalexiou, S. M., Y. Markonis, F. Lombardo, A. AghaKouchak, and E. Foufoula‐Georgiou, Precise temporal Disaggregation Preserving Marginals and Correlations (DiPMaC) for stationary and non‐stationary processes, Water Resources Research, doi:10.1029/2018WR022726, 2018.

  1. G. Papaioannou, A. Efstratiadis, L. Vasiliades, A. Loukas, S.M. Papalexiou, A. Koukouvinos, I. Tsoukalas, and P. Kossieris, An operational method for Floods Directive implementation in ungauged urban areas, Hydrology, 5 (2), 24, doi:10.3390/hydrology5020024, 2018.

    An operational framework for flood risk assessment in ungauged urban areas is developed within the implementation of the EU Floods Directive in Greece, and demonstrated for Volos metropolitan area, central Greece, which is frequently affected by intense storms causing fluvial flash floods. A scenario-based approach is applied, accounting for uncertainties of key modeling aspects. This comprises extreme rainfall analysis, resulting to spatially-distributed Intensity-Duration-Frequency (IDF) relationships and their confidence intervals, and flood simulations, through the SCS-CN method and the unit hydrograph theory, producing design hydrographs at the sub-watershed scale, for several soil moisture conditions. The propagation of flood hydrographs and the mapping of inundated areas are employed by the HEC-RAS 2D model, with flexible mesh size, by representing the resistance caused by buildings through the local elevation rise method. For all hydrographs, upper and lower estimates on water depths, flow velocities and inundation areas are estimated, for varying roughness coefficient values. The methodology is validated against the flood event of the 9th October 2006, using observed flood inundation data. Our analyses indicate that although typical engineering practices for ungauged basins are subject to major uncertainties, the hydrological experience may counterbalance the missing information, thus ensuring quite realistic outcomes.

    Full text: http://www.itia.ntua.gr/en/getfile/1829/1/documents/hydrology-05-00024_Idnk8fW.pdf (5243 KB)

    Other works that reference this work (this list might be obsolete):

    1. Petroselli, A., M. Vojtek, and J. Vojteková, Flood mapping in small ungauged basins: A comparison of different approaches for two case studies in Slovakia, Hydrology Research, doi:10.2166/nh.2018.040, 2018.
    2. Manfreda, S., C. Samela, A. Refice, V. Tramutoli, and F. Nardi, Advances in large-scale flood monitoring and detection, Hydrology, 5(3), 49, doi:10.3390/hydrology5030049, 2018.

  1. I. Tsoukalas, A. Efstratiadis, and C. Makropoulos, Stochastic periodic autoregressive to anything (SPARTA): Modelling and simulation of cyclostationary processes with arbitrary marginal distributions, Water Resources Research, 54 (1), 161–185, WRCR23047, doi:10.1002/2017WR021394, 2018.

    Stochastic models in hydrology traditionally aim at reproducing the empirically derived statistical characteristics of the observed data rather than any specific distribution model that attempts to describe the usually non-Gaussian statistical behavior of the associated processes. SPARTA (Stochastic Periodic AutoRegressive To Anything) offers an alternative and novel approach which allows the explicit representation of each process of interest with any distribution model, while simultaneously establishes dependence patterns that cannot be fully captured by the typical linear stochastic schemes. Cornerstone of the proposed approach is the Nataf joint-distribution model, which is related with the Gaussian copula, combined with Gaussian periodic autoregressive processes. In order to obtain the target stochastic structure, we have also developed a computationally simple and efficient algorithm, based on a hybrid Monte-Carlo procedure that is used to approximate the required equivalent correlation coefficients. Theoretical and practical benefits of the proposed method, contrasted to outcomes from widely used stochastic models, are demonstrated by means of real-world as well as hypothetical monthly simulation examples involving both univariate and multivariate time series.

    Additional material:

    Other works that reference this work (this list might be obsolete):

    1. Papalexiou, S. M., Unified theory for stochastic modelling of hydroclimatic processes: Preserving marginal distributions, correlation structures, and intermittency, Advances in Water Resources, 115, 234-252, doi:10.1016/j.advwatres.2018.02.013, 2018.

  1. A. Tegos, N. Malamos, A. Efstratiadis, I. Tsoukalas, A. Karanasios, and D. Koutsoyiannis, Parametric modelling of potential evapotranspiration: a global survey, Water, 9 (10), 795, doi:10.3390/w9100795, 2017.

    We present and validate a global parametric model of potential evapotranspiration (PET) with two parameters which are estimated through calibration, using as explanatory variables temperature and extraterrestrial radiation. The model and the parameters estimation approach were tested over the globe, using the FAO CLIMWAT database that provides monthly averaged values of meteorological inputs at 4300 locations worldwide. A preliminary analysis of these data allowed explaining the major drivers of PET over the globe and across seasons. Next, we developed an automatic optimization software tool to calibrate the model and provide point PET estimations against the given Penman-Monteith values. We also employed extended analysis of model inputs and outputs, including the production of global maps of optimized model parameters and associated performance metrics. Also, we employed interpolated values of the optimized parameters to validate the predictive capacity of our model against monthly meteorological time series, at several stations worldwide. The results were very encouraging, since even with the use of abstract climatic information for model calibration and the use of interpolated parameters as local predictors, the model generally ensures reliable PET estimations. In few cases the model performs poorly in estimating the reference PET, due to irregular interactions between temperature and extraterrestrial radiation, as well as because the associated processes are influenced by additional drivers, e.g. relative humidity and wind speed. However, the analysis of the residuals showed that the model is consistent in terms of parameters estimation and model validation. The provided parameters maps allow the direct use of the parametric model wherever in the world, providing PET estimates in case of missing data, that can be further improved even with a short term acquisition of meteorological data.

    Full text: http://www.itia.ntua.gr/en/getfile/1738/2/documents/water-09-00795.pdf (6428 KB)

    Additional material:

    See also: http://www.mdpi.com/2073-4441/9/10/795

    Works that cite this document: View on Google Scholar or ResearchGate

    Other works that reference this work (this list might be obsolete):

    1. Elferchichi, A., G. A. Giorgio, N. Lamaddalena, M. Ragosta, and V. Telesca, Variability of temperature and its impact on reference evapotranspiration: the test case of the Apulia region (Southern Italy), Sustainability, 9(12), 2337, doi:10.3390/su9122337, 2017.
    2. Li, M., R. Chu, S. Shen, and A. R. T. Islam, Quantifying climatic impact on reference evapotranspiration trends in the Huai River Basin of Eastern China, Water, 10(2), 144, doi:10.3390/w10020144, 2018.
    3. Yan, N., F. Tian, B. Wu, W. Zhu, and M. Yu, Spatiotemporal analysis of actual evapotranspiration and its causes in the Hai basin, Remote Sensing, 10(2), 332; doi:10.3390/rs10020332, 2018.
    4. Li, M., R. Chu, A.R.M.T. Islam, and S. Shen, Reference evapotranspiration variation analysis and its approaches evaluation of 13 empirical models in sub-humid and humid regions: A case study of the Huai River Basin, Eastern China, Water, 10(4), 493, doi:10.3390/w10040493, 2018.
    5. Hao, X., S. Zhang, W. Li, W. Duan, G. Fang, Y. Zhang , and B. Guo, The uncertainty of Penman-Monteith method and the energy balance closure problem, Journal of Geophysical Research – Atmospheres, 123(14), 7433-7443, doi:10.1029/2018JD028371, 2018.
    6. Giménez, P. O., and S. G. García-Galiano, Assessing Regional Climate Models (RCMs) ensemble-driven reference evapotranspiration over Spain, Water, 10(9), 1181, doi:10.3390/w10091181, 2018.

  1. C. Makropoulos, E. Rozos, I. Tsoukalas, and et al., Sewer-mining: A water reuse option supporting circular economy, public service provision and entrepreneurship, Journal of Environmental Management, doi:10.1016/j.jenvman.2017.07.026, 2017, (in press).

    Water scarcity, either due to increased urbanisation or climatic variability, has motivated societies to reduce pressure on water resources mainly by reducing water demand. However, this practice alone is not sufficient to both protect resources and guarantee the quality of life water services underpin especially within a context of increased urbanisation. As such, the idea of water reuse has been gaining momentum for some time in the water sector and has recently found a more general context within the emerging concept of the Circular Economy. As a result of this growing trend, water recycling schemes at various scales have been applied worldwide. The most common scale of water reuse is reusing the effluent of a wastewater treatment plant for irrigation or industrial uses (e.g. cooling towers, or rinsing). This is favoured by economies of scale, but to be economically viable it requires that the recycled-water user is close enough to the treatment plant (and at a more or less similar or lower elevation), otherwise capital and operational costs for transmission getratherhigh. Another downside with this scale of (centralised) reuse is that this scheme does not break the monopoly of water supply, since it is again the water company that runs the treatment unit and provides the effluent for reuse and as such offers reduced benefits in terms of job creation, innovation drive and entrepreneurship. On the other side of the scale spectrum, at the level of the household, reuse options include mostly the reuse of grey water for non-potable uses (such as toilet flushing and garden irrigation). Although promising and with significant potential for demand reduction, this scale of reuse is not necessarily cost effective, with all costs borne by the end user, and usually relies on additional motivation, such as drought conditions or environmental attitudes to be implemented. This study argues for an intermediate scale of water reuse, termed sewer-mining, which is a water recycling scheme at the neighbourhood scale. We suggest it provides a feasible alternative reuse option when the geography of the wastewater treatment plant is problematic, it relies on mature treatment technologies and presents an excellent opportunity for Small Medium Enterprises (SME) to be involved in the water supply market, thus securing both environmental, social and economic benefits (including but not restricted to water for ecosystem services). To support this argument, we report on a pilot sewer mining application. The pilot, integrates to important subsystems: a packaged treatment unit and an Information and Communications Technology (ICT) infrastructure that would allow an operator to manage remotely several sewer mining units thus rendering the provided service economically viable even for SMEs. The paper reports on the pilot’s overall performance and critically evaluates the potential of the sewer mining idea to become a significant piece of the circular economy puzzle for water.

  1. E. Rozos, I. Tsoukalas, K. Ripis, E. Smeti, and C. Makropoulos, Turning black into green: Ecosystem services from treated wastewater, Desalination and Water Treatment, 91 (2017), 2017.

    To reduce the impact of urban effluents on the environment, strict regulatory requirements have been set up for the disposal of wastewater, in most parts of the western world, requiring treatment before disposal. At the same time, the urban environment requires water inflows to satisfy a range of urban water demands, and the corresponding water abstractions put pressure on (often scarce) water resources. A suggested synergistic solution is to use the effluents from treatment plants as an alternative resource for irrigation or for industrial uses. Despite the existence of numerous successful applications, this practice is not very common mainly because of increased capital and operational costs, usually exceeding the cost of fresh water. A possible response of the market to this drawback could be to introduce in-situ small scale treatment units to cover local water needs. In this study, we assess the benefits of such a compact wastewater treatment unit that is used to provide water for irrigating an urban green area. Apart from the aesthetic improvement, the evaporative cooling (latent heat), which reduces the air temperature, is expected to have a positive impact on thermal comfort. A pilot scheme was deployed in KEREFYT, the research centre of the Athens Water Supply and Sewerage Company (EYDAP). This scheme was simulated with the UWOT model to estimate heat fluxes and the results were fed into Energy2D (a model that simulates heat transfer) to estimate the expected temperature drop. The results are promising and suggest that these technologies could play an important role in a more sustainable, circular water economy.

    Full text: http://www.itia.ntua.gr/en/getfile/1715/1/documents/Manuscript_subm2_CM.pdf (636 KB)

  1. I. Tsoukalas, P. Kossieris, A. Efstratiadis, and C. Makropoulos, Surrogate-enhanced evolutionary annealing simplex algorithm for effective and efficient optimization of water resources problems on a budget, Environmental Modelling and Software, 77, 122–142, doi:10.1016/j.envsoft.2015.12.008, 2016.

    In water resources optimization problems, the objective function usually presumes to first run a simulation model and then evaluate its outputs. However, long simulation times may pose significant barriers to the procedure. Often, to obtain a solution within a reasonable time, the user has to substantially restrict the allowable number of function evaluations, thus terminating the search much earlier than required. A promising strategy to address these shortcomings is the use of surrogate modeling techniques. Here we introduce the Surrogate-Enhanced Evolutionary Annealing-Simplex (SEEAS) algorithm that couples the strengths of surrogate modeling with the effectiveness and efficiency of the evolutionary annealing-simplex method. SEEAS combines three different optimization approaches (evolutionary search, simulated annealing, downhill simplex). Its performance is benchmarked against other surrogate-assisted algorithms in several test functions and two water resources applications (model calibration, reservoir management). Results reveal the significant potential of using SEEAS in challenging optimization problems on a budget.

    Related works:

    • [19] Early presentation if EGU conference

    Full text: http://www.itia.ntua.gr/en/getfile/1587/2/documents/SEEAS_paper.pdf (4310 KB)

    Additional material:

    Other works that reference this work (this list might be obsolete):

    1. Dariane , A. B., and M. M. Javadianzadeh, Towards an efficient rainfall–runoff model through partitioning scheme, Water, 8, 63; doi:10.3390/w8020063, 2016.
    2. Yaseen, Z. M., O. Jaafar, R. C. Deo, O. Kisi, J. Adamowski, J. Quilty, and A. El-Shafie, Boost stream-flow forecasting model with extreme learning machine data-driven: A case study in a semi-arid region in Iraq, Journal of Hydrology, 542, 603-614, doi:10.1016/j.jhydrol.2016.09.035, 2016.
    3. Müller, R., and N. Schütze, Multi-objective optimization of multi-purpose multi-reservoir systems under high reliability constraints, Environmental Earth Sciences, 75:1278, doi:10.1007/s12665-016-6076-5, 2016.
    4. #Christelis, V., V. Bellos, and G. Tsakiris, Employing surrogate modelling for the calibration of a 2D flood simulation model, Sustainable Hydraulics in the Era of Global Change: Proceedings of the 4th IAHR Europe Congress (Liege, Belgium, 27-29 July 2016), A. S. Erpicum, M. Pirotton, B. Dewals, P. Archambeau (editors), CRC Press, 2016.
    5. Salazar, J. Z., P. M. Reed, J. D. Quinn, M. Giuliani, and A. Castelletti, Balancing exploration, uncertainty and computational demands in many objective reservoir optimization, Advances in Water Resources, 109, 196-210, doi:10.1016/j.advwatres.2017.09.014, 2017.
    6. Christelis, V., and A. Mantoglou, Physics-based and data-driven surrogate models for pumping optimization of coastal aquifers, European Water, 57, 481–488, 2017.
    7. #Thandayutham, K., E. Avital, N. Venkatesan, and A. Samad, Design and analysis of a marine current turbine, Proceedings of ASME 2017 Gas Turbine India Conference and Exhibition, GTINDIA2017-4912, V001T02A014, Bangalore, India, doi:10.1115/GTINDIA2017-4912, 2017.
    8. Christelis, V., R. G. Regis, and A. Mantoglou, Surrogate-based pumping optimization of coastal aquifers under limited computational budgets, Journal of Hydroinformatics, 20(1), 164-176, doi:10.2166/hydro.2017.063, 2018.
    9. Psarrou, E., I. Tsoukalas, and C. Makropoulos, A Monte-Carlo-based method for the optimal placement and operation scheduling of sewer mining units in urban wastewater networks, Water, 10(2), 200, doi:10.3390/w10020200, 2018.
    10. #Thandayutham, K., L. K. Mishra, and A. Samad, Optimal design of a marine current turbine using CFD & FEA, 4th International Conference in Ocean Engineering, Chennai, India, 2018.
    11. Christelis, V., and A. G. Hughes, Metamodel-assisted analysis of an integrated model composition: an example using linked surface water – groundwater models, Environmental Modelling and Software, 107, 298-306, doi:10.1016/j.envsoft.2018.05.004, 2018.
    12. Zischg, A. P., G. Felder, M. Mosimann, V. Röthlisberger, and R. Weingartner, Extending coupled hydrological-hydraulic model chains with a surrogate model for the estimation of flood losses, Environmental Modelling and Software, 108, 174-185, doi:10.1016/j.envsoft.2018.08.009, 2018.

  1. I. Tsoukalas, and C. Makropoulos, A surrogate based optimization approach for the development of uncertainty-aware reservoir operational rules: the case of Nestos hydrosystem, Water Resources Management, 29 (13), 4719–4734, doi:10.1007/s11269-015-1086-8, 2015.

    Operation of large-scale hydropower reservoirs is a complex problem that involves conflicting objectives, such as hydropower generation and water supply. Deriving optimal operational rules is a challenging task due to the non-linearity of the system dynamics and the uncertainty of future inflows and water demands. A common approach to derive optimal control policies is to couple simulation models with optimization algorithms. This paper in order to investigate the performance of a future reservoir and safely infer about its significance employs stochastic simulation, thus long synthetically generated time-series and a multi-objective version of the Parameterization-Simulation-Optimization (PSO) framework to develop uncertainty-aware operational rules. Furthermore, in order to handle the high computational effort that ensues from that coupling we investigate the potential of a surrogate-based multi-objective optimization algorithm, ParEGO. The PSO framework is deployed with WEAP21 water resources management model as simulation engine and MATLAB for the implementation of optimization algorithms. A comparison between NSGAII and ParEGO optimization algorithms is performed to assess the effectiveness of the proposed algorithm. The aforementioned comparison showed that ParEGO provides efficient approximations of the Pareto front while reducing the computational effort required. Finally, the potential benefit and the significance of the future reservoir is underlined.

    Full text: http://www.itia.ntua.gr/en/getfile/1569/1/documents/tsoukalas_WRM.pdf (2008 KB)

    Other works that reference this work (this list might be obsolete):

    1. Müller, R., and N. Schütze, Multi-objective optimization of multi-purpose multi-reservoir systems under high reliability constraints, Environmental Earth Sciences, 75:1278, doi:10.1007/s12665-016-6076-5, 2016.
    2. Christelis, V., and A. G. Hughes, Metamodel-assisted analysis of an integrated model composition: an example using linked surface water – groundwater models, Environmental Modelling and Software, doi:10.1016/j.envsoft.2018.05.004, 2018.

  1. I. Tsoukalas, and C. Makropoulos, Multiobjective optimisation on a budget: Exploring surrogate modelling for robust multi-reservoir rules generation under hydrological uncertainty, Environmental Modelling and Software, 69, 396–413, doi:10.1016/j.envsoft.2014.09.023, 2015.

    Developing long term operation rules for multi-reservoir systems is complicated due to the number of decision variables, the non-linearity of system dynamics and the hydrological uncertainty. This uncertainty can be addressed by coupling simulation models with multi-objective optimisation algorithms driven by stochastically generated hydrological timeseries but the computational effort required imposes barriers to the exploration of the solution space. The paper addresses this by (a) employing a parsimonious multi-objective parameterization-simulation-optimization (PSO) framework, which incorporates hydrological uncertainty through stochastic simulation and allows the use of probabilistic objective functions and (b) by investigating the potential of multi-objective surrogate based optimisation (MOSBO) to significantly reduce the resulting computational effort. Three MOSBO algorithms are compared against two multi-objective evolutionary algorithms. Results suggest that MOSBOs are indeed able to provide robust, uncertainty-aware operation rules much faster, without significant loss of neither the generality of evolutionary algorithms nor of the knowledge embedded in domain-specific models.

    Other works that reference this work (this list might be obsolete):

    1. Müller, R., and N. Schütze, Multi-objective optimization of multi-purpose multi-reservoir systems under high reliability constraints, Environmental Earth Sciences, 75:1278, doi:10.1007/s12665-016-6076-5, 2016.
    2. Christelis, V., and A. G. Hughes, Metamodel-assisted analysis of an integrated model composition: an example using linked surface water – groundwater models, Environmental Modelling and Software, doi:10.1016/j.envsoft.2018.05.004, 2018.

Book chapters and fully evaluated conference publications

  1. I. Tsoukalas, C. Makropoulos, and A. Efstratiadis, Stochastic simulation of periodic processes with arbitrary marginal distributions, 15th International Conference on Environmental Science and Technology (CEST2017), Rhodes, Global Network on Environmental Science and Technology, 2017.

    Stochastic simulation of hydrological processes has a key role in water resources planning and management due to its ability to incorporate hydrological uncertainty within decision-making. Due to seasonality, the statistical characteristics of such processes are considered periodic functions, thus implying the use of cyclostationary stochastic models, typically using a common statistical distribution. Yet, this may not be representative of the statistical structure of such processes across all seasons. In this context, we introduce a novel model suitable for the simulation of periodic processes with arbitrary marginal distributions, called Stochastic Periodic AutoRegressive To Anything (SPARTA). Apart from capturing the periodic correlation structure of the underlying processes, its major advantages are a) the accurate preservation of seasonally-varying marginal distributions; b) the explicit generation of non-negative values; and c) the parsimonious model structure. Finally, the performance of the model is demonstrated through a theoretical (artificial) case study.

    Full text: http://www.itia.ntua.gr/en/getfile/1731/1/documents/cest2017_00797_oral_paper_V2.pdf (655 KB)

    Additional material:

    See also: http://cest.gnest.org/sites/default/files/presentation_file_list/cest2017_00797_oral_paper.pdf

  1. E. Rozos, I. Tsoukalas, K. Ripis, E. Smeti, and C. Makropoulos, Turning black into green: ecosystem services from treated wastewater, 13th IWA Specialized Conference on Small Water and Wastewater Systems, Athens, Greece, National Technical University of Athens, 2016, (in press).

    In order to reduce the impact of the urban effluents on the environment, modern societies have imposed restrictions regarding the quality of the disposals. For this reason, in the majority of the western world cities, the wastewater is treated before disposal. However, on the other side of the urban water cycle, water abstractions keep putting an increasing pressure on the water resources. As a countermeasure, treated wastewater is used occasionally as an alternative resource by employing large scale infrastructure to treat and supply water for either irrigation or industrial uses. Despite the existence of numerous successful applications, this practice is not very common mainly because of the increased capital and operational costs, usually exceeding the cost of fresh water. The response of the market to this drawback was to introduce in-situ small scale treatment units to cover local water needs. In this study, we assess the benefits of a compact wastewater treatment unit that is used to provide water for irrigating a green area. Apart from the aesthetic improvement, benefits are expected because of the evaporative cooling (latent heat), which reduce the air temperature. A pilot scheme was set up in KEREFYT, the research centre of Athens water supply company. This scheme was simulated with UWOT model to estimate the heat fluxes and the results were fed into Energy2D (a model that simulates heat transfer) to estimate the expected temperature drop.

    Full text: http://www.itia.ntua.gr/en/getfile/1600/1/documents/Manuscript_QiNArbH.pdf (509 KB)

  1. I. Tsoukalas, P. Dimas, and C. Makropoulos, Hydrosystem optimization on a budget: Investigating the potential of surrogate based optimization techniques, 14th International Conference on Environmental Science and Technology (CEST2015), Global Network on Environmental Science and Technology, University of the Aegean, 2015.

    Development of uncertainty-aware operational rules for multi-reservoir systems is a demanding and challenging task due to the complexity of the system dynamics, the number of decision variables and the hydrological uncertainty. In order to overcome this issue the parsimonious parameterization-simulation-optimization (PSO) framework is employed coupled with stochastically generated hydrological time-series. However, when the simulation model requires long computational time this coupling imposes a computational barrier to the framework. The purpose of this paper is threefold: a) Investigate the potential of Efficient Global Optimization (EGO) algorithm (and its variants) which is capable of reaching global optima within a few simulation model evaluations (~500 or less). b) Extend the capabilities of WEAP21 water resources management model by using it within PSO framework (named WEAP21-PSO) and c) Validate and compare the results of WEAP21-PSO using the well-known hydrosystem management model Hydronomeas coupled with Evolutionary Annealing Simplex (EAS) optimization algorithm. Results confirm that EGO has the potential and the capabilities to handle computationally demanding problems and furthermore is capable of locating the optimal solution within few simulation model evaluations and that the WEAP21-PSO framework performs well at the task at hand.

    Full text: http://www.itia.ntua.gr/en/getfile/1574/1/documents/cest2015_00162_oral_paper.pdf (475 KB)

    See also: http://cest.gnest.org/cest15proceedings/public_html/papers/cest2015_00162_oral_paper.pdf

Conference publications and presentations with evaluation of abstract

  1. I. Tsoukalas, P. Kossieris, A. Efstratiadis, C. Makropoulos, and D. Koutsoyiannis, CastaliaR: An R package for multivariate stochastic simulation at multiple temporal scales, European Geosciences Union General Assembly 2018, Geophysical Research Abstracts, Vol. 20, Vienna, EGU2018-18433, doi:10.13140/RG.2.2.20978.81605, European Geosciences Union, 2018.

    In contrast to great advances on stochastic simulation techniques in hydrology and their importance on water management and uncertainty assessment studies, operational software packages for generating synthetic data are limited and hardly accessible. This limits their adoption to a narrow audience, excluding the vast majority of researchers and practitioners. In an effort to bridge this gap, we introduce CastaliaR package that constitutes the R-based, open-source implementation of a state-of-the art methodology for multivariate stochastic simulation. Its background builds upon the works of Koutsoyiannis and Manetas (1996), Koutsoyiannis (1999, 2000) and Efstratiadis et al. (2014). Briefly, the overall scheme reproduces the statistical characteristics of the historical data at three temporal scales (annual, monthly and daily). The generation procedure lies upon a symmetric moving average process for the annual scale and a periodic autoregressive process for the finer scales, while a Monte Carlo disaggregation approach re-establishes consistency across the three temporal scales.

    Full text:

  1. Y. Moustakis, P. Kossieris, I. Tsoukalas, and A. Efstratiadis, Quasi-continuous stochastic simulation framework for flood modelling, European Geosciences Union General Assembly 2017, Geophysical Research Abstracts, Vol. 19, Vienna, 19, EGU2017-534, European Geosciences Union, 2017.

    Typically, flood modelling in the context of everyday engineering practices is addressed through event-based deterministic tools, e.g., the well-known SCS-CN method. A major shortcoming of such approaches is the ignorance of uncertainty, which is associated with the variability of soil moisture conditions and the variability of rainfall during the storm event. In event-based modeling, the sole expression of uncertainty is the return period of the design storm, which is assumed to represent the acceptable risk of all output quantities (flood volume, peak discharge, etc.). On the other hand, the varying antecedent soil moisture conditions across the basin are represented by means of scenarios (e.g., the three AMC types by SCS), while the temporal distribution of rainfall is represented through standard deterministic patterns (e.g., the alternative blocks method). In order to address these major inconsistencies,simultaneously preserving the simplicity and parsimony of the SCS-CN method, we have developed a quasi-continuous stochastic simulation approach, comprising the following steps: (1) generation of synthetic daily rainfall time series; (2) update of potential maximum soil moisture retention, on the basis of accumulated five-day rainfall; (3) estimation of daily runoff through the SCS-CN formula, using as inputs the daily rainfall and the updated value of soil moisture retention;(4) selection of extreme events and application of the standard SCS-CN procedure for each specific event, on the basis of synthetic rainfall. This scheme requires the use of two stochastic modelling components, namely the CastaliaR model, for the generation of synthetic daily data, and the HyetosMinute model, for the disaggregation of daily rainfall to finer temporal scales. Outcomes of this approach are a large number of synthetic flood events, allowing for expressing the design variables in statistical terms and thus properly evaluating the flood risk.

    Full text: http://www.itia.ntua.gr/en/getfile/1680/2/documents/FINAL_Moustakis_EGU2017.pdf (1492 KB)

    Additional material:

  1. P. Kossieris, A. Efstratiadis, I. Tsoukalas, and D. Koutsoyiannis, Assessing the performance of Bartlett-Lewis model on the simulation of Athens rainfall, European Geosciences Union General Assembly 2015, Geophysical Research Abstracts, Vol. 17, Vienna, EGU2015-8983, doi:10.13140/RG.2.2.14371.25120, European Geosciences Union, 2015.

    Many hydrological applications require the use of long rainfall data across a wide range of fine time scales. To meet this necessity, stochastic approaches are usually employed for the generation of large number of rainfall events, following a Monte Carlo approach. In this framework, Bartlett-Lewis model (BL) is a key representative from the family of Poisson-cluster stochastic processes. Here, we examine the performance of three different versions of BL model, with number of parameters varying from 5 up to 7, in representing the characteristics of convective and frontal rainfall of Athens (Greece). Apart from the typical statistical characteristics that are explicitly preserved by the stochastic model (mean, variance, lag-1 autocorrelation, probability dry), we also attempt to preserve the statistical distribution of annual rainfall maxima, as well as two important temporal properties of the observed storm events, i.e. the duration of storms and the time distance between subsequent events. This task is not straightforward, given that these characteristics are not described in the theoretical equations of the model, but they should be empirically evaluated on the basis of synthetic data. The analysis is conducted on monthly basis and for multiple time scales, i.e. from hourly to daily. Further to that, we focus on the formulation of the calibration problem, by assessing the performance of the BL model against issues such as choice of statistics to preserve, time scales, distance metrics, etc.

    Full text:

    See also: http://dx.doi.org/10.13140/RG.2.2.14371.25120

    Other works that reference this work (this list might be obsolete):

    1. Li, X., A. Meshgi, X. Wang, J. Zhang, S. H. X. Tay, G. Pijcke, N. Manocha, M. Ong, M. T. Nguyen, and V. Babovic, Three resampling approaches based on method of fragments for daily-to-subdaily precipitation disaggregation, International Journal of Climatology, 38(Suppl.1), e1119-e1138, doi:10.1002/joc.5438, 2018.
    2. Park, J., C. Onof, and D. Kim, A hybrid stochastic rainfall model that reproduces rainfall characteristics at hourly through yearly time dcale, Hydrology and Earth System Science Discussions, doi:10.5194/hess-2018-267, 2018.

  1. A. Efstratiadis, I. Tsoukalas, P. Kossieris, G. Karavokiros, A. Christofides, A. Siskos, N. Mamassis, and D. Koutsoyiannis, Computational issues in complex water-energy optimization problems: Time scales, parameterizations, objectives and algorithms, European Geosciences Union General Assembly 2015, Geophysical Research Abstracts, Vol. 17, Vienna, EGU2015-5121, doi:10.13140/RG.2.2.11015.80802, European Geosciences Union, 2015.

    Modelling of large-scale hybrid renewable energy systems (HRES) is a challenging task, for which several open computational issues exist. HRES comprise typical components of hydrosystems (reservoirs, boreholes, conveyance networks, hydropower stations, pumps, water demand nodes, etc.), which are dynamically linked with renewables (e.g., wind turbines, solar parks) and energy demand nodes. In such systems, apart from the well-known shortcomings of water resources modelling (nonlinear dynamics, unknown future inflows, large number of variables and constraints, conflicting criteria, etc.), additional complexities and uncertainties arise due to the introduction of energy components and associated fluxes. A major difficulty is the need for coupling two different temporal scales, given that in hydrosystem modeling, monthly simulation steps are typically adopted, yet for a faithful representation of the energy balance (i.e. energy production vs. demand) a much finer resolution (e.g. hourly) is required. Another drawback is the increase of control variables, constraints and objectives, due to the simultaneous modelling of the two parallel fluxes (i.e. water and energy) and their interactions. Finally, since the driving hydrometeorological processes of the integrated system are inherently uncertain, it is often essential to use synthetically generated input time series of large length, in order to assess the system performance in terms of reliability and risk, with satisfactory accuracy. To address these issues, we propose an effective and efficient modeling framework, key objectives of which are: (a) the substantial reduction of control variables, through parsimonious yet consistent parameterizations; (b) the substantial decrease of computational burden of simulation, by linearizing the combined water and energy allocation problem of each individual time step, and solve each local sub-problem through very fast linear network programming algorithms, and (c) the substantial decrease of the required number of function evaluations for detecting the optimal management policy, using an innovative, surrogate-assisted global optimization approach.

    Full text:

    See also: http://dx.doi.org/10.13140/RG.2.2.11015.80802

  1. A. Drosou, P. Dimitriadis, A. Lykou, P. Kossieris, I. Tsoukalas, A. Efstratiadis, and N. Mamassis, Assessing and optimising flood control options along the Arachthos river floodplain (Epirus, Greece), European Geosciences Union General Assembly 2015, Geophysical Research Abstracts, Vol. 17, Vienna, EGU2015-9148, European Geosciences Union, 2015.

    We present a multi-criteria simulation-optimization framework for the optimal design and setting of flood protection structures along river banks. The methodology is tested in the lower course of the Arachthos River (Epirus, Greece), downstream of the hydroelectric dam of Pournari. The entire study area is very sensitive, particularly because the river crosses the urban area of Arta, which is located just after the dam. Moreover, extended agricultural areas that are crucial for the local economy are prone to floods. In the proposed methodology we investigate two conflicting criteria, i.e. the minimization of flood hazards (due to damages to urban infrastructures, crops, etc.) and the minimization of construction costs of the essential hydraulic structures (e.g. dikes). For the hydraulic simulation we examine two flood routing models, named 1D HEC-RAS and quasi-2D LISFLOOD, whereas the optimization is carried out through the Surrogate-Enhanced Evolutionary Annealing-Simplex (SE-EAS) algorithm that couples the strengths of surrogate modeling with the effectiveness and efficiency of the EAS method.

    Full text:

  1. I. Tsoukalas, P. Kossieris, A. Efstratiadis, and C. Makropoulos, Handling time-expensive global optimization problems through the surrogate-enhanced evolutionary annealing-simplex algorithm, European Geosciences Union General Assembly 2015, Geophysical Research Abstracts, Vol. 17, Vienna, EGU2015-5923, European Geosciences Union, 2015.

    In water resources optimization problems, the calculation of the objective function usually presumes to first run a simulation model and then evaluate its outputs. In several cases, however, long simulation times may pose significant barriers to the optimization procedure. Often, to obtain a solution within a reasonable time, the user has to substantially restrict the allowable number of function evaluations, thus terminating the search much earlier than required by the problem’s complexity. A promising novel strategy to address these shortcomings is the use of surrogate modelling techniques within global optimization algorithms. Here we introduce the Surrogate-Enhanced Evolutionary Annealing-Simplex (SE-EAS) algorithm that couples the strengths of surrogate modelling with the effectiveness and efficiency of the EAS method. The algorithm combines three different optimization approaches (evolutionary search, simulated annealing and the downhill simplex search scheme), in which key decisions are partially guided by numerical approximations of the objective function. The performance of the proposed algorithm is benchmarked against other surrogate-assisted algorithms, in both theoretical and practical applications (i.e. test functions and hydrological calibration problems, respectively), within a limited budget of trials (from 100 to 1000). Results reveal the significant potential of using SE-EAS in challenging optimization problems, involving time-consuming simulations.

    Full text: