Search

Lockdowns and voluntary social distancing led to significant reduction in people’s mobility. Yet, there is scant evidence on the heterogeneous effects across segments of the population. Using unique mobility indicators based on anonymized and aggregate data provided by Vodafone for Italy, Portugal, and Spain, we find that lockdowns had a larger impact on the mobility of women and younger cohorts. Younger people also experienced a sharper drop in mobility in response to rising COVID-19 infections. Our findings, which are consistent across estimation methods and robust to a variety of tests, warn about a possible widening of gender and inter-generational inequality and provide important inputs for the formulation of targeted policies.

Reykjavik is almost entirely heated by geothermal energy. Yet, recent growth of the city significantly increased the heat demand. Past experiences in Iceland's capital region showed that hydraulic stimulation of existing geothermal wells is suited to improve hydraulic performance and energy supply. However, fluid injection may also trigger felt or even damaging earthquakes, which are of concern in populated areas and pose a significant risk to stimulation operations. Consequently, soft stimulation concepts have been developed to increase geothermal well performance while minimizing environmental effects such as induced seismicity. In a demonstration project of hydraulic soft stimulation in October 2019, more than 20.000 m³ of water were injected into well RV-43 in Reykjavik in multiple stages and with different injection schemes. The hydraulic performance of the well was improved without inducing felt seismicity. An a priori seismic risk assessment was conducted and for the first time the risk was continuously updated by an adaptive traffic light system supported by a sophisticated realtime microseismic monitoring. Our results confirm that it is possible to improve the performance of geothermal wells in Reykjavik and worldwide with acceptable technical, economic, and environmental risks. Here we provide an overview of the entire stimulation project including site description, stimulation design, zonal isolation, logging, seismic risk assessment and mitigation measures, realtime seismic, hydraulic and chemical monitoring, and stimulation results and challenges., Geothermics, 96, ISSN:0375-6505

Sound regional policies are essential for balanced and sustained economic growth. The interaction of federal and regional policies with cross-regional structural differences affects human and physical capital formation, the business climate, private investment, market depth, and competition. This paper summarizes the main elements of Russia’s fiscal federalism, describes the channels through which it operates, and assesses the effectiveness of regional transfers in reducing regional disparities. The results suggest that federal transfers to regions contributed to reducing disparities arising from heterogeneous regional tax bases and fiscal revenues. This allowed regions with initially lower per capita income to increase human and physical capital at higher rates. There is little evidence for transfers contributing to increased cross-regional growth synchronization. The results also suggest that federal transfers did not significantly improve regional fiscal sustainability, a conclusion that is supported by the lack of convergence in per capita real income across Russian regions in the last 15 years.

The rapid increase in energy demand in the city of Reykjavik has posed the need for an additional supply of deep geothermal energy. The deep-hydraulic (re-)stimulation of well RV-43 on the peninsula of Geldinganes (north of Reykjavik) is an essential component of the plan implemented by Reykjavik Energy to meet this energy target. Hydraulic stimulation is often associated with fluid-induced seismicity, most of which is not felt on the surface but which, in rare cases, can be a nuisance to the population and even damage the nearby building stock. This study presents a first-of-its-kind pre-drilling probabilistic induced seismic hazard and risk analysis for the site of interest. Specifically, we provide probabilistic estimates of peak ground acceleration, European microseismicity intensity, probability of light damage (damage risk), and individual risk. The results of the risk assessment indicate that the individual risk within a radius of 2 km around the injection point is below 0.1 micromorts, and damage risk is below 10−2, for the total duration of the project. However, these results are affected by several orders of magnitude of variability due to the deep uncertainties present at all levels of the analysis, indicating a critical need in updating this risk assessment with in situ data collected during the stimulation. Therefore, it is important to stress that this a priori study represents a baseline model and starting point to be updated and refined after the start of the project.

The 2014 April 1, M-w 8.1 Iquique (Chile) earthquake struck in the Northern Chile seismic gap. With a rupture length of less than 200 km, it left unbroken large segments of the former gap. Early studies were able to model the main rupture features but results are ambiguous with respect to the role of aseismic slip and left open questions on the remaining hazard at the Northern Chile gap. A striking observation of the 2014 earthquake has been its extensive preparation phase, with more than 1300 events with magnitude above M-L 3, occurring during the 15 months preceding the main shock. Increasing seismicity rates and observed peak magnitudes accompanied the last three weeks before the main shock. Thanks to the large data sets of regional recordings, we assess the precursor activity, compare foreshocks and aftershocks and model rupture preparation and rupture effects. To tackle inversion challenges for moderate events with an asymmetric network geometry, we use full waveforms techniques to locate events, map the seismicity rate and derive source parameters, obtaining moment tensors for more than 300 events (magnitudes M-w 4.0-8.1) in the period 2013 January 1-2014 April 30. This unique data set of fore- and aftershocks is investigated to distinguish rupture process models and models of strain and stress rotation during an earthquake. Results indicate that the spatial distributions of foreshocks delineated the shallower part of the rupture areas of the main shock and its largest aftershock, well matching the spatial extension of the aftershocks cloud. Most moment tensors correspond to almost pure double couple thrust mechanisms, consistent with the slab orientation. Whereas no significant differences are observed among thrust mechanisms in different areas, nor among thrust foreshocks and aftershocks, the early aftershock sequence is characterized by the presence of normal fault mechanisms, striking parallel to the trench but dipping westward. These events likely occurred in the shallow wedge structure close to the slab interface and are consequence of the increased extensional stress in this region after the largest events. The overall stress inversion result suggests a minor stress rotation after the main shock, but a significant release of the deviatoric stress. The temporal change in the distribution of focal mechanisms can also be explained in terms of the spatial heterogeneity of the stress field: under such interpretation, the potential of a large megathrust earthquake breaking a larger segment offshore Northern Chile remains high.

One of the most critical problems affecting geophysical data acquisition procedures is related to the misorientation of multicomponent sensors with respect to a common reference system (e.g., geographic north). In many applications, misoriented sensors affect data analysis procedures, leading to errors in results and interpretations. These problems generally occur in applications where the orientation of the sensor cannot be actively controlled and is not known a priori, e.g., geophysical sensors deployed in borehole installations or on the seafloor. We have developed a quaternion-based method for the optimal reorientation of multicomponent geophysical sensors. In contrast to other approaches, we took into account the full time-series record from all sensor components. Therefore, our method could be applied to all time-series data and was not restricted to a certain type of geophysical sensor. Our method allows the robust calculation of relative reorientations between two-component or three-component sensors. By using a reference sensor in an iterative process, this result can be extended to the estimation of absolute sensor orientations. In addition to finding an optimal solution for a full 3D sensor rotation, we have established a rigorous scheme for the estimation of uncertainties of the resulting orientation parameters. We tested the feasibility and applicability of our method using synthetic data examples for a vertical seismic profile and an ocean bottom seismometer array. We noted that the quaternion-based reorientation method is superior to the standard approach of a single-parameter estimation of rotation angles.

Automated location of seismic events is a very important task in microseismic monitoring operations as well for local and regional seismic monitoring. Since microseismic records are generally characterised by low signal-to-noise ratio, such methods are requested to be noise robust and sufficiently accurate. Most of the standard automated location routines are based on the automated picking, identification and association of the first arrivals of P and S waves and on the minimization of the residuals between theoretical and observed arrival times of the considered seismic phases. Although current methods can accurately pick P onsets, the automatic picking of the S onset is still problematic, especially when the P coda overlaps the S wave onset. In this thesis I developed a picking free automated method based on the Short-Term-Average/Long-Term-Average (STA/LTA) traces at different stations as observed data. I used the STA/LTA of several characteristic functions in order to increase the sensitiveness to the P wave and the S waves. For the P phases we use the STA/LTA traces of the vertical energy function, while for the S phases, we use the STA/LTA traces of the horizontal energy trace and then a more optimized characteristic function which is obtained using the principal component analysis technique. The orientation of the horizontal components can be retrieved by robust and linear approach of waveform comparison between stations within a network using seismic sources outside the network (chapter 2). To locate the seismic event, we scan the space of possible hypocentral locations and origin times, and stack the STA/LTA traces along the theoretical arrival time surface for both P and S phases. Iterating this procedure on a three-dimensional grid we retrieve a multidimensional matrix whose absolute maximum corresponds to the spatial and temporal coordinates of the seismic event. Location uncertainties are then estimated by perturbing the STA/LTA parameters (i.e the length of both long and short time windows) and relocating each event several times. In order to test the location method I firstly applied it to a set of 200 synthetic events. Then we applied it to two different real datasets. A first one related to mining induced microseismicity in a coal mine in the northern Germany (chapter 3). In this case we successfully located 391 microseismic event with magnitude range between 0.5 and 2.0 Ml. To further validate the location method I compared the retrieved locations with those obtained by manual picking procedure. The second dataset consist in a pilot application performed in the Campania-Lucania region (southern Italy) using a 33 stations seismic network (Irpinia Seismic Network) with an aperture of about 150 km (chapter 4). We located 196 crustal earthquakes (depth < 20 km) with magnitude range 1.1 < Ml < 2.7. A subset of these locations were compared with accurate locations retrieved by a manual location procedure based on the use of a double difference technique. In both cases results indicate good agreement with manual locations. Moreover, the waveform stacking location method results noise robust and performs better than classical location methods based on the automatic picking of the P and S waves first arrivals. Die automatische Lokalisierung seismischer Ereignisse ist eine wichtige Aufgabe, sowohl im Bereich des Mikroseismischen Monitorings im Bergbau und von Untegrund Aktivitäten, wie auch für die lokale und regionale Überwachung von natürlichen Erdbeben. Da mikroseismische Datensätze häufig ein schlechtes Signal-Rausch-Verhältnis haben müssen die Lokalisierungsmethoden robust gegen Rauschsignale und trotzdem hinreichend genau sein. Aufgrund der in der Regel sehr hochfrequent aufgezeichneten Messreihen und der dadurch sehr umfangreichen Datensätze sind automatische Auswertungen erstrebenswert. Solche Methoden benutzen in der Regel automatisch gepickte und den P und S Phasen zugeordnete Ersteinsätze und Minimieren die Summe der quadratischen Zeitdifferenz zwischen den beobachteten und theoretischen Einsatzzeiten. Obgleich das automatische Picken der P Phase in der Regel sehr genau möglich ist, hat man beim Picken der S Phasen häufig Probleme, z.B. wenn die Coda der P Phase sehr lang ist und in den Bereich der S Phase hineinreicht. In dieser Doktorarbeit wird eine Methode vollautomatische, Wellenform-basierte Lokalisierungsmethode entwickelt, die Funktionen des Verhältnisses "Short Term Average / Long Term Average" (STA/LTA) verwendet und keine Pickzeiten invertiert. Die STA/LTA charakteristische Funktion wurde für unterschiedliche Wellenform Attribute getestet, um die Empfindlichkeit für P und S Phasen zu erhöhen. Für die P Phase wird die STA/LTA Funktion für die Energie der Vertikalkomponente der Bodenbewegung benutzt, wohingegen für die S Phase entweder die Energie der horizontalen Partikelbewegung oder eine optimierte Funktion auf Basis der Eigenwertzerlegung benutzt wird. Um die Ereignisse zu lokalisieren wird eine Gittersuche über alle möglichen Untergrundlokalisierungen durchgeführt. Für jeden räumlichen und zeitlichen Gitterpunkt werden die charakteristischen Funktionen entlang der theoretischen Einsatzkurve aufsummiert. Als Ergebnis erhält man eine 4-dimensionale Matrix über Ort und Zeit des Ereignisses, deren Maxima die wahrscheinlichsten Lokalisierungen darstellen. Um die Unsicherheiten der Lokalisierung abzuschätzen wurden die Parameter der STA/LTA Funktionen willkürlich verändert und das Ereignis relokalisiert. Die Punktwolke aller möglichen Lokalisierungen gibt ein Maß für die Unsicherheit des Ergebnisses. Die neu entwickelte Methode wurde an einem synthetischen Datensatz von 200 Ereignissen getestet und für zwei beobachtete Datensätze demonstriert. Der erste davon betrifft induzierte Seismizität in einem Kohlebergbau in Norddeutschland. Es wurden 391 Mikrobeben mit Magnituden zwischen Ml 0.5 und 2.0 erfolgreich lokalisiert und durch Vergleich mit manuell ausgewerteten Lokalisierungen verifziert.Der zweite Datensatz stammt von einem Anwednung auf des Regionale Überwachungsnetz in der Region Campania-Lucania (Süditalien) mit 33 seismischen Stationen und einer Apertur von etwa 150 km. Wir konnten 196 Erdbeben mit Tiefen < 20 km und Magnituden zwischen Ml 1.1 und 2.7 lokalisieren. Eine Untergruppe der eigenen Lokalisierungen wurde mit den Lokalisierungen einer Standard Lokalisierung sowie einer hochgenauen Relativlokalisierung verglichen. In beiden Fällen ist die Übereinstimmung mit den manuellen Lokalisierungen groß. Außerdem finden wir, dass die Wellenform Summations Lokalisierung ronbust gegen Rauschen ist und bessere Ergebnisse liefert als die Standard Lokalisierung, die auf dem automatischen Picken von Ersteinsatzzeiten alleine basiert.

SUMMARYEarthquake location is one of the oldest problems in seismology, yet remains an active research topic. With dense seismic monitoring networks, it is possible to obtain reliable locations for microearthquakes; however, in many cases dense networks are lacking, limiting the location accuracy, or preventing location when there are too few observations. For small events in all settings, recording may be sparse and location may be difficult due to low signal-to-noise ratio. In this work, we introduce a new, distance-geometry-based method to locate seismicity clusters using only one or two seismic stations. A distance geometry problem consists in determining the location of sets of points based only on the distances between member pairs. Applied to seismology, our approach allows earthquake location using the interevent distance between earthquake pairs, which can be estimated using only one or two seismic stations. We first validate the method with synthetic data that resemble common cluster shapes, and then test the method with two seismic sequences in California: the August 2014 Mw 6.0 Napa earthquake and the July 2019 Mw 6.4 Ridgecrest earthquake sequence. We demonstrate that our approach provides robust and reliable results even for a single station. When using two seismic stations, the results capture the same structures recovered with high-resolution double-difference locations based on multiple stations. The proposed method is particularly useful for poorly monitored areas, where only a limited number of stations are available.

Since the end of 2013, the region around the two volcanoes Masaya and Momotombo, which includes the Nicaraguan capital Managua, has shown an unusually high seismic and volcanic activity. In December 2015, the Momotombo volcano erupted after 110 years of quiescence. Since mid-December 2015, the Masaya volcano has also shown gradually increasing activity, including the formation of a lava lake in its main crater. By adding 30 broadband stations, we had temporarily (December 2016–March 2017) densified the permanent Nicaraguan seismic network around these volcanoes to study the local seismicity and image the subsurface structure. During the observation period, we observed an overall low level of seismicity. Recorded events around Momotombo likely consist of aftershocks of the M5.5 earthquake that struck this area on September, 26th, 2016. At Masaya, we did not observe volcano-tectonic events. Using the continuous waveform recordings, we perform a 3D ambient seismic noise tomography that reveals a first image of the subsurface velocity structure below the Masaya and Momotombo volcanoes. While Momotombo shows a typical elongated low shear-wave velocity anomaly that reaches depths of about 8 km, Masaya does not show indications of a deep plumbing system. At Masaya, we have indications of a shallow (0–3 km) magmatic chamber, offset to the west and not directly below the active Santiago vent, where the crater lake is located At greater depth (3–8 km) a low velocity anomaly towards the northeast coincides in location with a modelled positive gravity anomaly and could indicate the presence of a former intrusive body. With this study we want to trigger further interest in the diverse tectonic and volcanic features of Nicaragua. Future, long-term seismic imaging and monitoring projects are of critical interest for the estimation of seismic and volcanic risks in Managua and the surroundings.