Partial ages: diagnosing transport processes by means of multiple clocks

The concept of age is widely used to quantify the transport rate of tracers - or pollutants - in the environment. The age focuses only on the time taken to reach a given location and disregards other aspects of the path followed by the tracer parcel. To keep track of the subregions visited by the tracer parcel along this path, partial ages are defined as the time spent in the different subregions. Partial ages can be computed in an Eulerian framework in much the same way as the usual age by extending the Constituent oriented Age and Residence Time theory (CART, www.climate.be/CART). In addition to the derivation of theoretical results and properties of partial ages, applications to a 1D model with lateral/transient storage, to the 1D advection-diffusion equation and to the diagnosis of the ventilation of the deep ocean are provided. They demonstrate the versatility of the concept of partial age and the potential new insights that can be gained with it.     

Probabilistic modelling of macroseismic attenuation and forecast of damage scenarios

According to the idea now widespread that macroseismic intensity should be expressed in probabilistic terms, a beta-binomial model has been proposed in the literature to estimate the probability of the intensity at site in the Bayesian framework and a clustering procedure has been adopted to define learning sets of macroseismic fields required to assign prior distributions of the model parameters. This article presents the results concerning the learning sets obtained by exploiting the large Italian macroseismic database DBM1I11 (Locati et al. in DBMI11, the 2011 version of the Italian Macroseismic Database, 2011. http://emidius.mi.ingv.it/DBMI11/) and discusses the problems related to their use in probabilistic modelling of the attenuation in seismic regions of the European countries partners of the UPStrat-MAFA project (2012), namely South Iceland, Portugal, SE Spain and Mt Etna volcano area (Italy). Anisotropy and the presence of offshore earthquakes are some of the problems faced. All the work has been carried out in the framework of the Task B of the project.     

On the 485-day Mode in the Atmospheric Angular Momentum: Spectral Analysis of IERS Data

The modification of spectral analysis especially intended for studying the disturbing functions of the atmosphere and ocean, as well as the observed polar motion (Wiener–Liouville spectrum), is used. The time series of the atmospheric disturbing functions obtained by the U.S. National Centers for Environmental Prediction (NCEP) of the International Earth Rotation and Reference Systems Service (IERS) for the period from January 1, 1980 to June 20, 2014 (http://www.iers.org/.cs1?pid=43-1100116) are analyzed. It is shown that the baric disturbing function contains a regular mode with a period of ~16 months; the contribution of this mode in the polar motion is estimated.     

Subsurface Hydrology of the Lake Chad Basin from Convection Modelling and Observations

In the Lake Chad basin, the quaternary phreatic aquifer (named hereafter QPA) presents large piezometric anomalies referred to as domes and depressions whose depths are ~15 and ~60 m, respectively. A previous study (Leblanc et al. in Geophys Res Lett, 2003, doi: 10.1029/2003GL018094) noticed that brightness temperatures from METEOSAT infrared images of the Lake Chad basin are correlated with the QPA piezometry. Indeed, at the same latitude, domes are ~4–5 K warmer than the depressions. Leblanc et al. (Geophys Res Lett, 2003, doi: 10.1029/2003GL018094) suggested that such a thermal behaviour results from an evapotranspiration excess above the piezometric depressions, an interpretation implicitly assuming that the QPA is separated from the other aquifers by the clay-rich Pliocene formation. Based on satellite visible images, here we find evidence of giant polygons, an observation that suggests instead a local vertical connectivity between the different aquifers. We developed a numerical water convective model giving an alternative explanation for the development of QPA depressions and domes. Beneath the depressions, a cold descending water convective current sucks down the overlying QPA, while, beneath the dome, a warm ascending current produces overpressure. Such a basin-wide circulation is consistent with the water geochemistry. We further propose that the thermal diurnal and evaporation/condensation cycles specific to the water ascending current explain why domes are warmer. We finally discuss the possible influence of the inferred convective circulation on the transient variations of the QPA reported from observations of piezometric levels and GRACE-based water mass change over the region.     

Intensity predictive attenuation models calibrated in Mw for metropolitan France

In the framework of the SIGMA project, a study was launched to develop a parametric earthquake catalog for the historical period, covering the metropolitan territory and calibrated in Mw. A set of candidate calibration events was selected corresponding to earthquakes felt over a part of the French metropolitan territory, which are fairly well documented both in terms of macroseismic intensity distributions (SisFrance BRGM-EDF-IRSN) and magnitude estimates. The detailed analysis of the macroseismic data led us to retain only 30 events out of 65 with Mw ranging from 3.6 to 5.8. In order to supplement the dataset with data from larger magnitude events, Italian earthquakes were also considered (11 events posterior to 1900 with Mw?≥?6.0 out of 15 in total), using both the DBMI11 macroseismic database (Locati et al. in Seismol Resour Lett 85(3):727–734, 2014) and the parametric information from the CPTI11 (Rovida et al. in CPTI11, la versione 2011 del Catalogo Parametrico dei Terremoti Italiani Istituto Nazionale di Geofisica et Vulcanologia, Milano, Bologna, 2011.  https://doi.org/10.6092/ingv.it-cpti11). To avoid introducing bias related to the differences in terms of intensity scales (MSK vs. MCS), only intensities smaller than or equal to VII were considered (Traversa et al. in On the use of cross-border macroseismic data to improve the estimation of past earthquakes seismological parameters, 2014). Mw and depth metadata were defined according to the Si-Hex catalogue (Cara et al. in Bull Soc Géol Fr 186:3–19, 2015.  https://doi.org/10.2113/qssqfbull.186.1.3), published information, and to the specific worked conducted within SIGMA related to early instrumental recordings (Benjumea et al. in Study of instrumented earthquakes that occurred during the first part of the 20th century (1905–1962), 2015). For the depth estimates, we also performed a macroseismic analysis to evaluate the range of plausible estimates and check the consistency of the solutions. Uncertainties on the metadata related to the calibration earthquakes were evaluated using the range of available alternative estimates. The intensity attenuation models were developed using a one-step maximum likelihood scheme. Several mathematical formulations and sub-datasets were considered to evaluate the robustness of the results (similarly to Baumont and Scotti in Accounting for data and modeling uncertainties in empirical macroseismic predictive equations (EMPEs). Towards “European” EMPEs based on SISFRANCE, DBMI, ECOS macroseismic database, 2008). In particular, as the region of interest may be characterized by significant laterally varying attenuation properties (Bakun and Scotti in Geophys J Int 164:596–610, 2006; Gasperini in Bull Seismol Soc Am 91:826–841, 2001), we introduced regional attenuation terms to account for this variability. Two zonation schemes were tested, one at the national scale (France/Italy), another at the regional scale based on the studies of Mayor et al. (Bull Earthq Eng, 2017.  https://doi.org/10.1007/s10518-017-0124-8) for France and Gasperini (2001) for Italy. Between and within event residuals were analyzed in detail to identify the best models, that is, the ones associated with the best misfit and most limited residual trends with intensity and distance. This analysis led us to select four sets of models for which no significant trend in the between- and within-event residuals is detected. These models are considered to be valid over a wide range of Mw covering?~?3.5–7.0.     

Cokriging for multivariate Hilbert space valued random fields: application to multi-fidelity computer code emulation

In this paper we propose Universal trace co-kriging, a novel methodology for interpolation of multivariate Hilbert space valued functional data. Such data commonly arises in multi-fidelity numerical modeling of the subsurface and it is a part of many modern uncertainty quantification studies. Besides theoretical developments we also present methodological evaluation and comparisons with the recently published projection based approach by Bohorquez et al. (Stoch Environ Res Risk Assess 31(1):53–70, 2016.  https://doi.org/10.1007/s00477-016-1266-y). Our evaluations and analyses were performed on synthetic (oil reservoir) and real field (uranium contamination) subsurface uncertainty quantification case studies. Monte Carlo analyses were conducted to draw important conclusions and to provide practical guidelines for all future practitioners.     

A perfect prognosis approach for daily precipitation series in consideration of space–time correlation structure

Downscaling techniques are the required tools to link the global climate model outputs provided at a coarse grid resolution to finer scale surface variables appropriate for climate change impact studies. Besides the at-site temporal persistence, the downscaled variables have to satisfy the spatial dependence naturally observed between the climate variables at different locations. Furthermore, the precipitation spatial intermittency should be fulfilled. Because of the complexity in describing these properties, they are often ignored, which can affect the effectiveness of the hydrologic process modeling. This study is a continuation of the work by Khalili and Nguyen (Clim Dyn 49(7–8):2261–2278.  https://doi.org/10.1007/s00382-016-3443-6, 2017) regarding the multi-site statistical downscaling of daily precipitation series. Different approach of multi-site statistical downscaling based on the concept of the spatial autocorrelation is presented in this paper. This approach has proven to give effective results for multi-site multivariate statistical downscaling of daily extreme temperature time series (Khalili et al. in Int J Climatol 33:15–32.  https://doi.org/10.1002/joc.3402, 2013). However, more challenges are presented by the precipitation variable because of the high spatio-temporal variability and intermittency. The proposed approach consists of logistic and multiple regression models, linking the global climate predictors to the precipitation occurrences and amounts respectively, and using the spatial autocorrelation concept to reproduce the spatial dependence observed between the precipitation series at different sites. An empirical technique has also been involved in this approach in order to fulfill the precipitation intermittency property. The proposed approach was performed using observed daily precipitation data from ten weather stations located in the southwest region of Quebec and southeast region of Ontario in Canada, and climate predictors from the NCEP/NCAR (National Centers for Environmental Prediction/National Centre for Atmospheric Research) reanalysis dataset. The results have proven the ability of the proposed approach to adequately reproduce the observed precipitation occurrence and amount characteristics, temporal and spatial dependence, spatial intermittency and temporal variability.     

The French seismic CATalogue (FCAT-17)

In regions that undergo low deformation rates, as is the case for metropolitan France (i.e. the part of France in Europe), the use of historical seismicity, in addition to instrumental data, is necessary when dealing with seismic hazard assessment. This paper presents the strategy adopted to develop a parametric earthquake catalogue using moment magnitude M_w, as the reference magnitude scale to cover both instrumental and historical periods for metropolitan France. Work performed within the framework of the SiHex (SIsmicité de l’HEXagone) (Cara et al. Bull Soc Géol Fr 186:3–19, 2015. doi: 10.2113/qssqfbull.186.1.3) and SIGMA (SeIsmic Ground Motion Assessment; EDF-CEA-AREVA-ENEL) projects, respectively on instrumental and historical earthquakes, have been combined to produce the French seismic CATalogue, version 2017 (FCAT-17). The SiHex catalogue is composed of ~40,000 natural earthquakes, for which the hypocentral location and M_w magnitude are given. In the frame of the SIGMA research program, an integrated study has been realized on historical seismicity from intensity prediction equations (IPE) calibration in M_w detailed in Baumont et al. (submitted) companion paper to their application to earthquakes of the SISFRANCE macroseismic database (BRGM, EDF, IRSN), through a dedicated strategy developed by Traversa et al. (Bull Earthq Eng, 2017. doi: 10.1007/s10518-017-0178-7) companion paper, to compute their M_w magnitude and depth. Macroseismic data and epicentral location and intensity used both in IPE calibration and inversion process, are those of SISFRANCE without any revision. The inversion process allows the main macroseismic field specificities reported by SISFRANCE to be taken into account with an exploration tree approach. It also allows capturing the epistemic uncertainties associated with macroseismic data and to IPEs selection. For events that exhibit a poorly constrained macroseismic field (mainly old, cross border or off-shore earthquakes), joint inversion of M_w and depth is not possible, and depth needs to be fixed to calculate M_w. Regional a priori depths have been defined for this purpose based on analysis of earthquakes with a well constrained macroseismic field where joint inversion of M_w and depth is possible. As a result, 27% of SISFRANCE earthquake seismological parameters have been jointly inverted and for the other 73% M_w has been calculated assuming a priori depths. The FCAT-17 catalogue is composed of the SIGMA historical parametric catalogue (magnitude range between 3.5 up to 7.0), covering from AD463 to 1965, and of the SiHex instrumental one, extending from 1965 to 2009. Historical part of the catalogue results from an automatic inversion of SISFRANCE data. A quality index is estimated for each historical earthquake according to the way the events are processed. All magnitudes are given in M_w which makes this catalogue directly usable as an input for probabilistic or deterministic seismic hazard studies. Uncertainties on magnitudes and depths are provided for historical earthquakes following calculation scheme presented in Traversa et al. (2017). Uncertainties on magnitudes for instrumental events are from Cara et al. (J Seismol 21:551–565, 2017. doi: 10.1007/s10950-016-9617-1).     

Proposal of an empirical site classification method based on target simulated horizontal over vertical spectral ratio

Nowadays, most of the site classifications schemes are based on the predominant period of the site as determined from the average horizontal to vertical spectral ratios of seismic motion or microtremor. However, the difficulty lies in the identification of the predominant period in particular if the observed average response spectral ratio does not present a clear peak but rather a broadband amplification or multiple peaks. In this work, based on the Eurocode-8 (2004) site classification, and assuming bounded random fields for both shear and compression waves-velocities, damping coefficient, natural period and depth of soil profile, one propose a new site-classification approach, based on “target” simulated average \( H/V \) spectral ratios, defined for each soil class. Taking advantage of the relationship of Kawase et al. (Bull Seismol Soc Am 101:2001–2014, 2011), which link the \( H/V \) spectral ratio to the horizontal (\( HTF \)) over the vertical (\( VTF \)) transfer functions, statistics of \( H/V \) spectral ratio via deterministic visco-elastic seismic analysis using the wave propagation theory are computed for the 4 soil classes. The obtained results show that \( H/V \) and \( HTF \) have amplitudes and shapes remarkably different among the four soil classes and exhibit fundamental peaks in the period ranges remarkably similar. Moreover, the “target” simulated average \( H/V \) spectral ratios for the 4 soil classes are in good agreement with the experimental ones obtained by Zhao et al. (Bull Seismol Soc Am 96:914–925, 2006) from the abundant and reliable Japanese strong motions database Kik-net, Ghasemi et al. (Soil Dyn Earthq Eng 29:121–132, 2009) from the Iranian strong motion data, and Di Alessandro et al. (Bull Sesismol Soc Am 106:2, 2011.  https://doi.org/10.1785/0120110084) from the Italian strong motion data. In addition to the 4 EC-8 standard soil classes (A, B, C and D), the superposition of the 4 target \( H/V \) reveals 3 new boundary site classes; AB, BC and CD, for overlapping \( V_{s,30} \) ranges when the predominant peak is not clearly consistent with any of the 4 proposed classes. Finally, one proposes a site classification index based on the ratio between the cross-correlation and the mean quadratic error between the in situ \( H/V \) spectral ratio and the “target” one. In order to test the reliability of the proposed approach, data from 139 sites were used, 132 collected from the Kik-net network database from Japan and 7 from Algeria. The site classification success rate per site class are around 93, 82, 89 and 100% for rock, hard soil, medium soil and soft soil, respectively. Zhao et al. (2006) found an average success for the 4 classes of soil close to 60%, similar to what one found in the present study (63%) without considering the new soil classes, but much smaller if one considers them (86%). In the absence of \( V_{s,30} \) data, the proposed approach can be an alternative to site classification.     

Decaying grid turbulence experiments in a stratified fluid: flow measurements and statistics

Laboratory experiments of decaying grid stratified turbulence were performed in a two-layer fluid and varying the stratification intensity. Turbulence was generated by towing an array of cylinders in a square vessel and the grid was moved at a constant velocity along the total vertical extent of the tank. In order to investigate the influence of the stratification intensity on the turbulence decay, both 2C-PIV and stereo PIV were used to provide time resolved velocity fields in the horizontal plane and the out-of-plane velocity. As expected, a faster decay of the turbulence level along the vertical axis and the collapse in a quasi-horizontal motion increased with the buoyancy frequency, N. In order to characterise the decay process we investigated the time evolution of the vortex statistics, the turbulence scales and the kinetic energy and enstrophy of the horizontal flow. The exponents recovered in the corresponding scaling laws were compared with the theoretical predictions and with reference values obtained in previous experimental studies. Both the spectral analysis and the evolution of characteristic length scales indicate that, in the examined range of N, the dynamics is substantially independent of the stratification intensity. The results obtained were explained in terms of the scaling analysis of decaying turbulence in strongly stratified fluids introduced by Brethouwer et al. (J Fluid Mech 585:343–368.  https://doi.org/10.1017/S0022112007006854, 2007).     

Reports from the Smithsonian’s Global Volcanism Network,May 2009

Information included in this summary is based on more detailed reports published in the Bulletin of the Global Volcanism Network, vol. 34, no. 5, May 2009 (on the Internet at http://www.volcano.si.edu/). Edited by scientists at the Smithsonian, this bulletin includes reports provided by a worldwide network of correspondents. The reports contain the names and contact information for all sources. Please note that these reports are preliminary and subject to change as events are studied in more detail. The Global Volcanism Program welcomes further reports of current volcanism, seismic unrest, monitoring data, and field observations.     

Reports from the Smithsonian’s Global Volcanism Network,March 2009

Information included in this summary is based on more detailed reports published in the Bulletin of the Global Volcanism Network, vol. 34, no. 3, March 2009 (on the Internet at http://www.volcano.si.edu/). Edited by scientists at the Smithsonian, this bulletin includes reports provided by a worldwide network of correspondents. The reports contain the names and contact information for all sources. Please note that these reports are preliminary and subject to change as events are studied in more detail. The Global Volcanism Program welcomes further reports of current volcanism, seismic unrest, monitoring data, and field observations.     

Ground-motion characterization for the probabilistic seismic hazard assessment in Turkey

This study describes the methodology implemented to establish the ground-motion logic-tree for national probabilistic seismic hazard map of Turkey for shallow active crustal regions. The presented procedure provides quantitative information to guide the hazard experts while establishing the logic tree to capture the epistemic uncertainty in ground-motion characterization. It uses non-data-driven and data-driven testing methods to identify and rank candidate ground-motion prediction equations (GMPEs) under a specific ground-motion database. The candidate GMPEs are subjected to visual inspection and are classified into center, body and range (CBR) spectral estimates for a proper consideration of epistemic uncertainty. The GMPEs classified into CBR are then used in a suite of seismic hazard sensitivity analysis to establish the most suitable GMPE logic-tree whose spectral estimates are not biased by any one of the GMPEs in the logic-tree structure. The sensitivity analysis considers normalized spectral ordinates and is not manipulated by the spectral amplitudes. The proposed procedure is inherited from the relevant studies of the Earthquake Model of the Middle East (EMME; www.efehr.org:8080/jetspeed/portal/emme.psml) regional seismic hazard project. This paper also highlights the similarities and differences in ground-motion characterization between EMME and our approach.     

Patuxent landscape model: 1. Hydrological model development

We developed a spatially explicit, process-based model of the 2352 km² Patuxent river watershed in Maryland, and its subwatersheds to integrate data and knowledge over several spatial, temporal and complexity scales, and to serve as an aid to regional management. The model was developed using the Library of Hydro-Ecological Modules (LHEM, http://giee.uvm.edu/LHEM), which was designed to create flexible landscape model structures that can be easily modified and extended to suit the requirements of a variety of goals and case studies. The LHEM includes modules that simulate various aspects of ecosystem dynamics. In this paper we consider modules that represent the physical conditions in the environment (climatic factors, geoporphology), and hydrologic processes, both locally and spatially. Where possible the modules are formulated as Stella(R) models, spatial transport processes are presented as C++ code.     

Water and Food in the Twenty-First Century

In 2000, the World population was 6.2 billion people; it reached 7 billion in 2012 and is expected to reach 9.5 billion (±0.4) in 2050 and 11 billion (±1.5) in 2100, according to the 2012 UN projections (Gerland et al. in Science 346:234–237, 2014). The trend after 2100 is still one of the global demographic growths, but after 2060, Africa is the only continent where the population would still increase. The amount of water consumed annually to produce the food necessary to meet the needs of the populations varies greatly between countries, from about 600 to 2500 m³/year per capita (Zimmer in L’empreinte eau. Les faces cachées d’une ressource vitale. Charles Léopold Meyer, Paris, 2013), depending on their wealth, their food habits, and the percentage of food waste they generate (on average, 30 % of the food produced is wasted). In 2000, the total food production was on the order of 3300 million tons (in cereal equivalents). In 2014, it is estimated that about 0.8 billion inhabitants of the planet suffer from hunger (FAO in World agriculture: towards 2030–2050. FAO, Rome, 2014. http://www.fao.org/docrep/004/Y3557E/y3557e00.HTM) and do not get the nutrition they need to be in good health or, in the case of children, to grow properly (both physically and intellectually). This food deficit was on the order of 40 million tons of cereal equivalents in 2014. The number of inhabitants with a food deficit was about 0.85 billion before the 2008 crisis and was decreasing annually, but it increased abruptly after 2008 up to 1 billion inhabitants and is slowly decreasing now. Assuming a World average water consumption for food of 1300 m³/year per capita in 2000, 1400 m³/year in 2050, and 1500 m³/year in 2100, a volume of water of around 8200 km³/year was needed in 2000, 13,000 km³/year will be needed in 2050, and 16,500 km³/year in 2100 (Marsily in L’eau, un trésor en partage. Dunod, Paris, 2009). Can bioenergy be added to food production? Will that much water be available on Earth, and where will it come from? Is climate change going to modify the answers to these questions? Can severe droughts occur? Can there be conflicts related to a food deficit? Some preliminary answers and scenarios for food production will be given in this paper from a hydrologist’s viewpoint.     

Burst regimes of long-period irregular pulsations at frequencies of 2.0–6.0 mHz and substorm activity in the nightside magnetosphere

The results of studying the simultaneous observations of burst regimes of long-period irregular pulsations at frequencies of 2.0–6.0 mHz (the series of ipcl bursts) in the region of the dayside polar cusp and magnetic field disturbances in the nightside auroral oval are presented. The data on the magnetic field at Mirny (MIR, Φ = 76.93°; Λ = 122.92°) and Yellowknife (YKC, Φ = 69.94°; Λ = 294.38°) antipodal observatories as well as the AE index values (http://www.cetp.ipsl.fr/~isgi/homepag1.htm) have been used in an analysis. It has been found out that 87% (group I) and 13% (group II) of events were registered against a back-ground of substorm activity and a quiet nightside magnetosphere, respectively. It has been revealed that several morphological characteristics of the group-I and -II ipcl bursts differ depending on the conditions in the nightside magnetosphere. It has been indicated that the intervals between peaks and the amplitudes of ipcl bursts of both types are distributed according to the exponential and power laws. The results indicate that magnetospheric plasma turbulence develops in the region where burst regimes are formed. It is assumed that the substorm processes in the magnetotail manifest themselves in plasma turbulence in the dayside cusp.     

Streamflow dynamics in a Mediterranean temporary river

Abstract

Mediterranean rivers are characterized by the irregularity of flow, harsh hydrological fluctuations and a profound transformation as the result of human activity. In this study, we investigate the streamflow response of a Mediterranean temporary river in which different groundwater, agriculture and urban contributions play an important role. Streamflow was measured at three nested gauging stations installed along Na Borges River, a lowland agricultural catchment (319 km²) on the island of Mallorca. Based on two hydrological years (2004/05 and 2005/06), potential evapotranspiration and surface water–groundwater interaction generated a succession of four different hydrological periods playing an important role in baseflow dynamics. The runoff coefficients were very low (<5%). At the event scale, groundwater also controlled runoff response, being very different according to hydrogeology, antecedent conditions and human impacts. During dry seasons, wastewater and karstic spring discharges maintain an influent regime into some streams. As a result, intense rainstorms in late summer generated water volumes over the impervious urban surfaces involved, with the result that quickflow was significant because the hydrological pathways were active. Citation Estrany, J., Garcia, C. & Alberich, R. (2010 Estrany, J., Garcia, C. and Batalla, R. J. 2010. Hydrological response of a small Mediterranean agricultural catchment. J. Hydrol., 380(1-2): 180–190. doi:10.1016/j.jhydrol.2009.10.035[Crossref] , [Google Scholar]) Streamflow dynamics in a Mediterranean temporary river. Hydrol. Sci. J. 55(5), 717–736.      

Numerical modeling of hydrodynamics and sediment transport—an integrated approach

Point measurement-based estimation of bedload transport in the coastal zone is very difficult. The only way to assess the magnitude and direction of bedload transport in larger areas, particularly those characterized by complex bottom topography and hydrodynamics, is to use a holistic approach. This requires modeling of waves, currents, and the critical bed shear stress and bedload transport magnitude, with a due consideration to the realistic bathymetry and distribution of surface sediment types. Such a holistic approach is presented in this paper which describes modeling of bedload transport in the Gulf of Gdańsk. Extreme storm conditions defined based on 138-year NOAA data were assumed. The SWAN model (Booij et al. 1999) was used to define wind–wave fields, whereas wave-induced currents were calculated using the Ko?odko and Gic-Grusza (2015) model, and the magnitude of bedload transport was estimated using the modified Meyer-Peter and Müller (1948) formula. The calculations were performed using a GIS model. The results obtained are innovative. The approach presented appears to be a valuable source of information on bedload transport in the coastal zone.     

The University of Athens Hellenic Macroseismic Database (HMDB.UoA): historical earthquakes

A systematic study of historical earthquakes leading to the quantification of earthquake effects in terms of macroseismic data points (MDPs) and, consequently, earthquake parameters has been carried out in the last decade at the Laboratory of Seismology of the University of Athens. For each earthquake, the available background information was evaluated and the corresponding macroseismic intensities assessed in terms of the European Macroseismic Scale 1998. A considerable amount of these MDPs contributed to the Archive of Historical Earthquake Data inventory through European initiatives (NERIES and SHARE). Based on the structure of the European Database, the local version of the Hellenic Macroseismic Database (HMDB.UoA) was designed incorporating historical earthquakes of the period 1000–1899 from the eastern Aegean area, central Greece and Ionian Islands. In its present form, the HMDB.UoA includes 90 events with I _max?≥?7 (868 MDPs) and 1,088 events with I _max?<?7 (1,273 MDPs). The database is hosted on the website http://macroseismology.geol.uoa.gr/.