Refine
Year of publication
- 2020 (203) (remove)
Document Type
- Article (140)
- Postprint (32)
- Doctoral Thesis (24)
- Monograph/Edited Volume (3)
- Other (3)
- Review (1)
Is part of the Bibliography
- yes (203)
Keywords
- climate change (7)
- model (5)
- modelling (5)
- models (5)
- Chinese loess (4)
- climate-change (4)
- precipitation (4)
- remote sensing (4)
- tectonics (4)
- time series analysis (4)
Institute
- Institut für Geowissenschaften (203) (remove)
In recent years, nature-based solutions are receiving increasing attention in the field of disaster risk reduction and climate change adaptation as inclusive, no regret approaches. Ecosystem-based adaptation (EbA) can mitigate the impacts of climate change, build resilience and tackle environmental degradation thereby supporting the targets set by the 2030 Agenda, the Paris Agreement and the Sendai Framework. Despite these benefits, EbA is still rarely implemented in practice. To better understand the barriers to implementation, this research examines policy-makers' perceptions of EbA, using an extended version of Protection Motivation Theory as an analytical framework. Through semi-structured interviews with policy-makers at regional and provincial level in Central Vietnam, it was found that EbA is generally considered a promising response option, mainly due to its multiple ecosystem-service benefits. The demand for EbA measures was largely driven by the perceived consequences of natural hazards and climate change. Insufficient perceived response efficacy and time-lags in effectiveness for disaster risk reduction were identified as key impediments for implementation. Pilot projects and capacity building on EbA are important means to overcome these perceptual barriers. This paper contributes to bridging the knowledge-gap on political decision-making regarding EbA and can, thereby, promote its mainstreaming into policy plans.
Controversy over the plate tectonic affinity and evolution of the Saxon granulites in a two- or multi-plate setting during inter- or intracontinental collision makes the Saxon Granulite Massif a key area for the understanding of the Palaeozoic Variscan orogeny. The massif is a large dome structure in which tectonic slivers of metapelite and metaophiolite units occur along a shear zone separating a diapir-like body of high-Pgranulite below from low-Pmetasedimentary rocks above. Each of the upper structural units records a different metamorphic evolution until its assembly with the exhuming granulite body. New age and petrologic data suggest that the metaophiolites developed from early Cambrian protoliths during high-Pamphibolite facies metamorphism in the mid- to late-Devonian and thermal overprinting by the exhuming hot granulite body in the early Carboniferous. A correlation of new Ar-Ar biotite ages with publishedP-T-tdata for the granulites implies that exhumation and cooling of the granulite body occurred at average rates of similar to 8 mm/year and similar to 80 degrees C/Ma, with a drop in exhumation rate from similar to 20 to similar to 2.5 mm/year and a slight rise in cooling rate between early and late stages of exhumation. A time lag ofc. 2 Ma between cooling through the closure temperatures for argon diffusion in hornblende and biotite indicates a cooling rate of 90 degrees C/Ma when all units had assembled into the massif. A two-plate model of the Variscan orogeny in which the above evolution is related to a short-lived intra-Gondwana subduction zone conflicts with the oceanic affinity of the metaophiolites and the timescale ofc. 50 Ma for the metamorphism. Alternative models focusing on the internal Variscan belt assume distinctly different material paths through the lower or upper crust for strikingly similar granulite massifs. An earlier proposed model of bilateral subduction below the internal Variscan belt may solve this problem.
A silica undersaturated alkali-olivine basanitic magma intruded the late Paleocene/early Eocene Jafnayn Formation near Muscat. Geochemical analyses indicate that a significant amount of host rock (limestone) was assimilated into the magma. We dated the basanite as 42.7 +/- 1.0 Ma (2 sigma error; late Lutetian), using the whole rock Ar-40/Ar-39 step-wise heating technique. Intrusion occurred in the hanging wall of a major regional extensional shear zone (Frontal Range Fault, FRF) bounding the northern margin of two domes within the Oman Mountains (Jabal Akhdar and Saih Hatat domes). Two shear intervals along the FRF have been documented. The first interval lasted immediately after emplacement of the Semail Ophiolite (latest Cretaceous-early Eocene) while the second and poorly constrained interval was assumed to have occurred during the Oligocene.
The proximity of the basanite to the FRF suggests that magma used extensional faults for the upper part of its ascent path. Reactivated Permian rift faults of the Pangaea rift or other preexisting faults may have been used for the lower ascent part.
We conclude that the basanite intrusion coincided with the onset of the second deformation interval along the FRF, because (1) the position of the basanite is near a dextral releasing bend, associated with the second shear interval, (2) the overlap of our Ar-40/Ar-39 age with the cooling curves for rocks from the nearby Jabal Akhdar Dome, and (3) the basanite postdates the first FRF deformation episode by > 10 Ma. Thus, the second interval along the FRF had started already during the late Lutetian and probably lasted into the Miocene.
A ground motion logic tree for seismic hazard analysis in the stable cratonic region of Europe
(2020)
Regions of low seismicity present a particular challenge for probabilistic seismic hazard analysis when identifying suitable ground motion models (GMMs) and quantifying their epistemic uncertainty. The 2020 European Seismic Hazard Model adopts a scaled backbone approach to characterise this uncertainty for shallow seismicity in Europe, incorporating region-to-region source and attenuation variability based on European strong motion data. This approach, however, may not be suited to stable cratonic region of northeastern Europe (encompassing Finland, Sweden and the Baltic countries), where exploration of various global geophysical datasets reveals that its crustal properties are distinctly different from the rest of Europe, and are instead more closely represented by those of the Central and Eastern United States. Building upon the suite of models developed by the recent NGA East project, we construct a new scaled backbone ground motion model and calibrate its corresponding epistemic uncertainties. The resulting logic tree is shown to provide comparable hazard outcomes to the epistemic uncertainty modelling strategy adopted for the Eastern United States, despite the different approaches taken. Comparison with previous GMM selections for northeastern Europe, however, highlights key differences in short period accelerations resulting from new assumptions regarding the characteristics of the reference rock and its influence on site amplification.
Evaluation of a novel application of earthquake HVSR in site-specific amplification estimation
(2020)
Ground response analyses (GRA) model the vertical propagations of SH waves through flat-layered media (1DSH) and are widely carried out to evaluate local site effects in practice. Horizontal-to-vertical spectral ratio (HVSR) technique is a cost-effective approach to extract certain site-specific information, e.g., site fundamental frequency (f(0)), but HVSR values cannot be directly used to approximate the levels of S-wave amplifications. Motivated by the work of Kawase et al. (2019), we propose a procedure to correct earthquake HVSR amplitudes for direct amplification estimations. The empirical correction compensates HVSR by generic vertical amplification spectra categorized by the vertical fundamental frequency (f(0v)) via kappa-means clustering. In this investigation, we evaluate the effectiveness of the corrected HVSR in approximating observed linear amplifications in comparison with 1DSH modellings. We select a total of 90 KiK-net (Kiban Kyoshin network) surface-downhole sites which are found to have no velocity contrasts below their boreholes and thus of which surface-to-borehole spectral ratios (SBSRs) can be taken as their empirical transfer functions (ETFs). 1DSH-based theoretical transfer functions (TTFs) are computed in the linear domain considering uncertainties in Vs profiles through randomizations. Five goodness-of-fit metrics are adopted to gauge the closeness between observed (ETF) and predicted (i.e., TTF and corrected HVSR) amplifications in both amplitude and spectral shape over frequencies from f(0) to 25 Hz. We find that the empirical correction to HVSR is highly effective and achieves a "good match" in both spectral shape and amplitude at the majority of the 90 KiK-net sites, as opposed to less than one-third for the 1DSH modelling. In addition, the empirical correction does not require a velocity model, which GRAs require, and thus has great potentials in seismic hazard assessments.
The structural evolution of calderas in rifts helps to characterize the spatiotemporal relationships between magmatism, long wavelength crustal deformation and the formation of tectonic deformation zones along the rift axis. We document the structural characteristics of the c. 36 ka old Menengai Caldera located within a young zone of extension in the central Kenya Rift. Field mapping and high-resolution digital surface models show that NNE-striking Holocene normal faults perpendicular to the regional ESE-WNWextension direction dominate the interior sectors of the rift. Inside the caldera, these structures are overprinted by post-collapse doming and faulting of the magmatic centre, resulting in obliquely slipping normal faults bounding a resurgence horst. Radiocarbon dating of faulted units as young as 5 ka cal BP and the palaeo-shorelines of a lake formed during the African Humid Period in the Nakuru Basin indicate that volcanism and fault activity inside and in the vicinity of Menengai must have been sustained during the Holocene. Our analysis confirms that the caldera is located at the centre of an extending rift segment and suggests that other magmatic centres and young zones of faulting along the volcano-tectonic axis of the Kenya Rift may constitute nucleation points of faulting that ultimately foster future continental break-up.
Machine learning (ML) algorithms are being increasingly used in Earth and Environmental modeling studies owing to the ever-increasing availability of diverse data sets and computational resources as well as advancement in ML algorithms. Despite advances in their predictive accuracy, the usefulness of ML algorithms for inference remains elusive. In this study, we employ two popular ML algorithms, artificial neural networks and random forest, to analyze a large data set of flood events across Germany with the goals to analyze their predictive accuracy and their usability to provide insights to hydrologic system functioning. The results of the ML algorithms are contrasted against a parametric approach based on multiple linear regression. For analysis, we employ a model-agnostic framework named Permuted Feature Importance to derive the influence of models' predictors. This allows us to compare the results of different algorithms for the first time in the context of hydrology. Our main findings are that (1) the ML models achieve higher prediction accuracy than linear regression, (2) the results reflect basic hydrological principles, but (3) further inference is hindered by the heterogeneity of results across algorithms. Thus, we conclude that the problem of equifinality as known from classical hydrological modeling also exists for ML and severely hampers its potential for inference. To account for the observed problems, we propose that when employing ML for inference, this should be made by using multiple algorithms and multiple methods, of which the latter should be embedded in a cross-validation routine.
Underground coal gasification (UCG) is an in situ conversion technique that enables the production of high-calorific synthesis gas from resources that are economically not minable by conventional methods. A broad range of end-use options is available for the synthesis gas, including fuels and chemical feedstock production. Furthermore, UCG also offers a high potential for integration with Carbon Capture and Storage (CCS) to mitigate greenhouse gas emissions. In the present study, a stoichiometric equilibrium model, based on minimization of the Gibbs function has been used to estimate the equilibrium composition of the synthesis gas. Thereto, we further developed and applied a proven thermodynamic equilibrium model to simulate the relevant thermochemical coal conversion processes (pyrolysis and gasification). Our modeling approach has been validated against thermodynamic models, laboratory gasification experiments and UCG field trial data reported in the literature. The synthesis gas compositions have been found to be in good agreement under a wide range of different operating conditions. Consequently, the presented modeling approach enables an efficient quantification of synthesis gas quality resulting from UCG, considering varying coal and oxidizer compositions at deposit-specific pressures and temperatures.
Earth and environmental sciences rely on detailed information about subsurface processes. Whereas geophysical techniques typically provide highly resolved spatial images, monitoring subsurface processes is often associated with enormous effort and, therefore, is usually limited to point information in time or space. Thus, the development of spatial and temporal continuous field monitoring methods is a major challenge for the understanding of subsurface processes. We have developed a novel method for ground-penetrating-radar (GPR) reflection monitoring of subsurface flow processes under unsaturated conditions and applied it to a hydrological infiltration experiment performed across a periglacial slope deposit in northwest Luxembourg. Our approach relies on a spatial and temporal quasicontinuous data recording and processing, followed by an attribute analysis based on analyzing differences between individual time steps. The results demonstrate the ability of time-lapse GPR monitoring to visualize the spatial and temporal dynamics of preferential flow processes with a spatial resolution in the order of a few decimeters and temporal resolution in the order of a few minutes. We observe excellent agreement with water table information originating from different boreholes. This demonstrates the potential of surface-based GPR reflection monitoring to observe the spatiotemporal dynamics of water movements in the subsurface. It provides valuable, and so far not accessible, information for example in the field of hydrology and pedology that allows studying the actual subsurface processes rather than deducing them from point information.
Groningen is the largest onshore gas field under production in Europe. The pressure depletion of the gas field started in 1963. In 1991, the first induced micro-earthquakes have been located at reservoir level with increasing rates in the following decades. Most of these events are of magnitude less than 2.0 and cannot be felt. However, maximum observed magnitudes continuously increased over the years until the largest, significant event with ML=3.6 was recorded in 2014, which finally led to the decision to reduce the production. This causal sequence displays the crucial role of understanding and modeling the relation between production and induced seismicity for economic planing and hazard assessment. Here we test whether the induced seismicity related to gas exploration can be modeled by the statistical response of fault networks with rate-and-state-dependent frictional behavior. We use the long and complete local seismic catalog and additionally detailed information on production-induced changes at the reservoir level to test different seismicity models. Both the changes of the fluid pressure and of the reservoir compaction are tested as input to approximate the Coulomb stress changes. We find that the rate-and-state model with a constant tectonic background seismicity rate can reproduce the observed long delay of the seismicity onset. In contrast, so-called Coulomb failure models with instantaneous earthquake nucleation need to assume that all faults are initially far from a critical state of stress to explain the delay. Our rate-and-state model based on the fluid pore pressure fits the spatiotemporal pattern of the seismicity best, where the fit further improves by taking the fault density and orientation into account. Despite its simplicity with only three free parameters, the rate-and-state model can reproduce the main statistical features of the observed activity.
Receiver function approaches have proven to be valuable for the investigation of crustal and upper mantle discontinuities whose sharp changes in seismic velocities cause wave conversions. While the crustal and mantle transition zone discontinuities are largely understood, the X-discontinuity at 250-350 km depth is still an object of controversial debate. The origin and global distribution of this structure with a velocity jump of 1.5-4.8% for compressional and shear waves is still unexplained. Although the crustal and mantle transition zone discontinuities beneath SW Morocco and surroundings have been investigated, only a few studies observed the X-discontinuity and place the depth at 260-370 km beneath the region of western Morocco. In order to better locate and characterize the X-discontinuity beneath southwest Morocco, we create P-wave receiver functions using data recorded by the Morocco-Munster array and detect the X-discontinuity at apparent depths of 285-350 km. In the western part of our study region we find apparent depths of similar to 310-340 km. The eastern part of the study area appears more complex: we locate two velocity jumps at apparent depths of around 285-295 km and 330-350 km in the northeast, and in the southeast we find a discontinuity at apparent depths of 340-350 km. Due to the large depth range and the twofold appearance of the X-discontinuity, we suggest that two different phase transitions cause the X-discontinuity beneath SW Morocco. The velocity contrasts at larger depths likely point to the coesite-stishovite phase transition occurring in deep eclogitic pools. The shallower depths can be explained by the transition from orthoenstatite to high-pressure clinoenstatite which requires the reaction between eclogite and peridotite to form orthopyroxene-rich peridotite. This reaction is likely related to previously proposed small-scale mantle upwellings beneath SW Morocco. Since both phase transitions require eclogite occurrence, the location of the X-discontinuity in this region can be used to indicate the location of recycled oceanic crust.
Subduction of Neo-Tethys oceanic lithosphere beneath the Iranian plate during the Mesozoic formed several igneous bodies of ultramafic to intermediate and felsic composition. Intrusion of these magmas into a regional metamorphic sequence (the Sanandaj-Sirjan Zone) caused partial melting and formation of migmatites with meta-pelitic protoliths. The Alvand complex (west Iran) is a unique area comprising migmatites of both mafic and pelitic protoliths. In this area, the gabbroic rocks contain veins of leucosome at their contact with pyroxenite and olivine gabbro. These leucosomes are geochemically and mineralogically different from leucosomes of the meta-pelitic migmatites and clearly show properties of I-type granites. Microscopic observations and whole rock compositions of the mafic migmatite leucosomes show that migmatization occurred through partial melting of biotite, hornblende and plagioclase. Thermobarometric calculations indicate 800 degrees C and 3.7 kbar for partial melting, although phase diagram modeling demonstrates that the presence of water could decrease the solidus temperature by about 40 degrees C. Our results suggest an asthenospheric magma upwelling as the source of heat for partial melting of the gabbroic rock during subduction of Neo-Tethys oceanic crust under the western edge of the Iranian plate. The present study also reveals relationships between migmatization and formation of S- and I -type granites in the area.
High-resolution surface velocities and strain for Anatolia from Sentinel-1 InSAR and GNSS data
(2020)
Measurements of present-day surface deformation are essential for the assessment of long-term seismic hazard. The European Space Agency's Sentinel-1 satellites enable global, high-resolution observation of crustal motion from Interferometric Synthetic Aperture Radar (InSAR). We have developed automated InSAR processing systems that exploit the first similar to 5 years of Sentinel-1 data to measure surface motions for the similar to 800,000-km(2) Anatolian region. Our new 3-D velocity and strain rate fields illuminate deformation patterns dominated by westward motion of Anatolia relative to Eurasia, localized strain accumulation along the North and East Anatolian Faults, and rapid vertical signals associated with anthropogenic activities and to a lesser extent extension across the grabens of western Anatolia. We show that automatically processed Sentinel-1 InSAR data can characterize details of the velocity and strain rate fields with high resolution and accuracy over large regions. These results are important for assessing the relationship between strain accumulation and release in earthquakes. <br /> Plain Language Summary Satellite-based measurements of small rates of motion of the Earth's surface made at high spatial resolutions and over large areas are important for many geophysical applications including improving earthquake hazard models. We take advantage of recent advances in geodetic techniques in order to measure surface velocities and tectonic strain accumulation across the Anatolia region, including the highly seismogenic and often deadly North Anatolian Fault. We show that by combining Sentinel-1 Interferometric Synthetic Aperture Radar (InSAR) data with Global Navigation Satellite System (GNSS) measurements we can enhance our view of surface deformation associated with active tectonics, the earthquake cycle, and anthropogenic processes.
An effective strategy for combining variance- and distribution-based global sensitivity analysis
(2020)
We present a new strategy for performing global sensitivity analysis capable to estimate main and interaction effects from a generic sampling design. The new strategy is based on a meaningful combination of varianceand distribution-based approaches. The strategy is tested on four analytic functions and on a hydrological model. Results show that the analysis is consistent with the state-of-the-art Saltelli/Jansen formula but to better quantify the interaction effect between the input factors when the output distribution is skewed. Moreover, the estimation of the sensitivity indices is much more robust requiring a smaller number of simulations runs. Specific settings and alternative methods that can be integrated in the new strategy are also discussed. Overall, the strategy is considered as a new simple and effective tool for performing global sensitivity analysis that can be easily integrated in any environmental modelling framework.
The 10th edition of the International Congress on the Application of Raman Spectroscopy in Art and Archaeology (RAA2019) was held in Potsdam (Germany) from 3 to 7 September 2019, with eight keynote lectures, 35 oral presentations and 18 Poster Presentations. The number of active participants was 68 delegates from 20 countries among the 236 authors that presented at least one work.
Thermal erosion is a major mechanism of permafrost degradation, resulting in characteristic landforms. We inventory thermo-erosional valleys in ice-rich coastal lowlands adjacent to the Siberian Laptev Sea based on remote sensing, Geographic Information System (GIS), and field investigations for a first regional assessment of their spatial distribution and characteristics. Three study areas with similar geological (Yedoma Ice Complex) but diverse geomorphological conditions vary in valley areal extent, incision depth, and branching geometry. The most extensive valley networks are incised deeply (up to 35 m) into the broad inclined lowland around Mamontov Klyk. The flat, low-lying plain forming the Buor Khaya Peninsula is more degraded by thermokarst and characterized by long valleys of lower depth with short tributaries. Small, isolated Yedoma Ice Complex remnants in the Lena River Delta predominantly exhibit shorter but deep valleys. Based on these hydrographical network and topography assessments, we discuss geomorphological and hydrological connections to erosion processes. Relative catchment size along with regional slope interact with other Holocene relief-forming processes such as thermokarst and neotectonics. Our findings suggest that thermo-erosional valleys are prominent, hitherto overlooked permafrost degradation landforms that add to impacts on biogeochemical cycling, sediment transport, and hydrology in the degrading Siberian Yedoma Ice Complex.
Earthquakes often rupture across more than one fault segment. If such rupture segmentation occurs on a significant scale, a simple point-source or one-fault model may not represent the rupture process well. As a consequence earthquake characteristics inferred, based on one-source assumptions, may become systematically wrong. This might have effects on follow-up analyses, for example regional stress field inversions and seismic hazard assessments. While rupture segmentation is evident for most M-w > 7 earthquakes, also smaller ones with 5.5 < M-w < 7 can be segmented. We investigate the sensitivity of globally available data sets to rupture segmentation and their resolution to reliably estimate the mechanisms in presence of segmentation. We focus on the sensitivity of InSAR (Interferometric Synthetic Aperture Radar) data in the static near-field and seismic waveforms in the far-field of the rupture and carry out non-linear and Bayesian optimizations of single-source and two-sources kinematic models (double-couple point sources and finite, rectangular sources) using InSAR and teleseismic waveforms separately. Our case studies comprises of four M-w 6-7 earthquakes: the 2009 L'Aquila and 2016 Amatrice (Italy) and the 2005 and 2008 Zhongba (Tibet) earthquakes. We contrast the data misfits of different source complexity by using the Akaike informational criterion (AIC). We find that the AIC method is well suited for data-driven inferences on significant rupture segmentation for the given data sets. This is based on our observation that an AIC-stated significant improvement of data fit for two-segment models over one-segment models correlates with significantly different mechanisms of the two source segments and their average compared to the single-segment mechanism. We attribute these modelled differences to a sufficient sensitivity of the data to resolve rupture segmentation. Our results show that near-field data are generally more sensitive to rupture segmentation of shallow earthquakes than far-field data but that also teleseismic data can resolve rupture segmentation in the studied magnitude range. We further conclude that a significant difference in the modelled source mechanisms for different segmentations shows that an appropriate choice of model segmentation matters for a robust estimation of source mechanisms. It reduces systematic biases and trade-off and thereby improves the knowledge on the rupture. Our study presents a strategy and method to detect significant rupture segmentation such that an appropriate model complexity can be used in the source mechanism inference. A similar, systematic investigation of earthquakes in the range of M-w 5.5-7 could provide important hazard-relevant statistics on rupture segmentation. In these cases single-source models introduce a systematic bias. Consideration of rupture segmentation therefore matters for a robust estimation of source mechanisms of the studied earthquakes.
The simulation of broad-band (0.1 to 10 + Hz) ground-shaking over deep and spatially extended sedimentary basins at regional scales is challenging. We evaluate the ground-shaking of a potential M 6.5 earthquake in the southern Lower Rhine Embayment, one of the most important areas of earthquake recurrence north of the Alps, close to the city of Cologne in Germany. In a first step, information from geological investigations, seismic experiments and boreholes is combined for deriving a harmonized 3D velocity and attenuation model of the sedimentary layers. Three alternative approaches are then applied and compared to evaluate the impact of the sedimentary cover on ground-motion amplification. The first approach builds on existing response spectra ground-motion models whose amplification factors empirically take into account the influence of the sedimentary layers through a standard parameterization. In the second approach, site-specific 1D amplification functions are computed from the 3D basin model. Using a random vibration theory approach, we adjust the empirical response spectra predicted for soft rock conditions by local site amplification factors: amplifications and associated ground-motions are predicted both in the Fourier and in the response spectra domain. In the third approach, hybrid physics-based ground-motion simulations are used to predict time histories for soft rock conditions which are subsequently modified using the 1D site-specific amplification functions computed in method 2. For large distances and at short periods, the differences between the three approaches become less notable due to the significant attenuation of the sedimentary layers. At intermediate and long periods, generic empirical ground-motion models provide lower levels of amplification from sedimentary soils compared to methods taking into account site-specific 1D amplification functions. In the near-source region, hybrid physics-based ground-motions models illustrate the potentially large variability of ground-motion due to finite source effects.
Thirteen N-butylpyridinium salts, including three monometallic [C4Py](2)[MCl4], nine bimetallic [C4Py](2)[(M1-xMxCl4)-M-a-Cl-b] and one trimetallic compound [C4Py](2)[(M1-y-zMyMz (c) Cl4)-M-a-M-b] (M=Co, Cu, Mn; x=0.25, 0.50 or 0.75 and y=z=0.33), were synthesized and their structure and thermal and electrochemical properties were studied. All compounds are ionic liquids (ILs) with melting points between 69 and 93 degrees C. X-ray diffraction proves that all ILs are isostructural. The conductivity at room temperature is between 10(-4) and 10(-8) S cm(-1). Some Cu-based ILs reach conductivities of 10(-2) S cm(-1), which is, however, probably due to IL dec. This correlates with the optical bandgap measurements indicating the formation of large bandgap semiconductors. At elevated temperatures approaching the melting points, the conductivities reach up to 1.47x10(-1) S cm(-1) at 70 degrees C. The electrochemical stability windows of the ILs are between 2.5 and 3.0 V.