Imaging seismic wave-fields with AlpArray and neighboring European networks

The AlpArray experiment and the deployment of Swath-D together with the dense permanent network in Italy allow for detailed imaging of the spatio-temporal imaging complexity of seismic wave-fields within the greater Alpine region. The distance of any point within the area to the nearest station is less than 30 km, resulting in an average inter-station distance of about 45 km. With a much denser deployment in a smaller region of the Alps (320 km in length and 140 km wide), the Swath-D network possesses an average inter-station distance of about 15 km. We show that seismogram sections with a spatial sampling of less than 5 km can be obtained using recordings of these regional arrays for just a single event. Multiply reflected body waves can be observed for up to 2 h after source time. In addition, we provide and describe animations of long-period seismic wave-fields using recordings of about 1300–1600 broadband stations for six representative earthquakes. These illustrate the considerable spatio-temporal variability of the wave-field’s properties at a high lateral resolution. Within denser station distributions like those provided by Swath-D, even shorter period body and surface wave features can be recovered. The decrease of the horizontal wavelength from to to surface waves, deviations from spherically symmetric wavefronts, and the capability to detect multi-orbit arrivals are demonstrated qualitatively by the presented wave-field animations, which are a valuable tool for educational, quality control, and research purposes. We note that the information content of the acquired datasets can only be adequately explored by application of appropriate quantitative methods accounting for the considerable complexity of the seismic wave-fields as revealed by the now available station configuration.

Atmospheric tomography using the Nordic Meteor Radar Cluster and Chilean Observation Network De Meteor Radars: network details and 3D-Var retrieval

Ground-based remote sensing of atmospheric parameters is often limited to single station observations by vertical profiles at a certain geographic location. This is a limiting factor for investigating gravity wave dynamics as the spatial information is often missing, e.g., horizontal wavelength, propagation direction or intrinsic frequency. In this study, we present a new retrieval algorithm for multistatic meteor radar networks to obtain tomographic 3-D wind fields within a pre-defined domain area. The algorithm is part of the Agile Software for Gravity wAve Regional Dynamics (ASGARD) and called 3D-Var, and based on the optimal estimation technique and Bayesian statistics. The performance of the 3D-Var retrieval is demonstrated using two meteor radar networks: the Nordic Meteor Radar Cluster and the Chilean Observation Network De Meteor Radars (CONDOR). The optimal estimation implementation provide statistically sound solutions and diagnostics from the averaging kernels and measurement response. We present initial scientific results such as body forces of breaking gravity waves leading to two counter-rotating vortices and horizontal wavelength spectra indicating a transition between the rotational k−3 and divergent k-5/3 mode at scales of 80–120 km. In addition, we performed a keogram analysis over extended periods to reflect the latitudinal and temporal impact of a minor sudden stratospheric warming in December 2019. Finally, we demonstrate the applicability of the 3D-Var algorithm to perform large-scale retrievals to derive meteorological wind maps covering a latitude region from Svalbard, north of the European Arctic mainland, to central Norway.

Large Roll Vortices Exhibited by Post-Tropical Cyclone Sandy during Landfall

Roll vortices are frequent features of a hurricane’s boundary layer, with kilometer or sub-kilometer horizontal scale. In this study, we found that large roll vortices with O (10 km) horizontal wavelength occurred over land in Post-Tropical Cyclone Sandy (2012) during landfall on New Jersey. Various characteristics of roll vortices were corroborated by analyses of Doppler radar observations, a 500 m resolution Weather Research and Forecasting (WRF) simulation, and an idealized roll vortex model. The roll vortices were always linear-shaped, and their wavelengths of 5–14 km were generally larger than any previously published for a tropical cyclone over land. Based on surface wind observations and simulated WRF surface wind fields, we found that roll vortices significantly increased the probability of hazardous winds and likely caused the observed patchiness of treefall during Sandy’s landfall.

Tomographic imaging of a large-scale travelling ionospheric disturbance during the Halloween storm of 2003

The most intense ionospheric storm observed in recent times occurred between 29 and 31 October 2003. The disturbances to the high-latitude regions set off several large-scale travelling ionospheric disturbances (LSTIDs), wave-like perturbations in the ionospheric electron density. This paper investigates one particular TID on 31 October 2003 using North American Global Positioning System (GPS) receiver network data and a tomographic imaging technique. The TID has an estimated period of 30 min and an estimated horizontal wavelength of 700 km and propagates south-westward over North America. The tomographic reconstruction of the wave is validated using a simulation of the observations and with independent observations from ionosondes and the CHAMP planar Langmuir probe. The results are discussed in the context of the magnetic and ionospheric conditions that may have contributed to the launch of the wave. Large-scale TIDs are challenging to study over large regions of the Earth, and the GPS network here is shown to offer a unique perspective on the spatial and temporal variation of the TID. The experimental results are backed up by simulations that show a denser network of receivers, as is available in more recent years, would produce improved accuracy in the TID imaging.

Environment and Mechanisms of Severe Turbulence in a Midlatitude Cyclone

A large midlatitude cyclone occurred over the central United States from 0000 to 1800 UTC 30 April 2017. During this period, there were more than 1100 reports of moderate-or-greater turbulence at commercial aviation cruising altitudes east of the Rocky Mountains. Much of this turbulence was located above or, otherwise, outside the synoptic-scale cloud shield of the cyclone, thus complicating its avoidance. In this study we use two-way nesting in a numerical model with finest horizontal spacing of 370 m to investigate possible mechanisms producing turbulence in two distinct regions of the cyclone. In both regions, model-parameterized turbulence kinetic energy compares well to observed turbulence reports. Despite being outside of hazardous large radar reflectivity locations in deep convection, both regions experienced strong modification of the turbulence environment as a result of upper-tropospheric/lower-stratospheric (UTLS) convective outflow. For one region, where turbulence was isolated and short lived, simulations revealed breaking of ~100-km horizontal-wavelength lower-stratospheric gravity waves in the exit region of a UTLS jet streak as the most likely mechanism for the observed turbulence. Although similar waves occurred in a simulation without convection, the altitude at which wave breaking occurred in the control simulation was strongly affected by UTLS outflow from distant deep convection. In the other analyzed region, turbulence was more persistent and widespread. There, overturning waves of much shorter 5–10-km horizontal wavelengths occurred within layers of gradient Richardson number < 0.25, which promoted Kelvin–Helmholtz instability associated with strong vertical shear in different horizontal locations both above and beneath the convectively enhanced UTLS jet.

Occurrence and Development of an Extreme Precipitation Event in the Ili Valley, Xinjiang, China and Analysis of Gravity Waves

We used observational data and the results from a high-resolution numerical simulation model to analyze the occurrence and development of an extreme precipitation event in the Ili Valley, Xinjiang, China on 26 June 2015. We analyzed the horizontal wavelength, period, speed, ducting, energy propagation and feedback mechanism of inertial gravity waves. A low-level convergence line was formed in the valley by the northerly and westerly winds as a result of Central Asian vortices and the trumpet-shaped topography of the Ili Valley. There was sufficient water vapor in the valley for the precipitation event to develop. A mesoscale vortex formed and developed on the low-level convergence line and the rainfall was distributed either near the convergence line or the mesoscale vortex. The low-level convergence line and the uplift caused by the terrain triggered convection, and then the convection triggered waves at lower levels. The combination of ascending motion induced by the lower level waves and the mesoscale vortex led to the development of convection, causing the precipitation to intensify. When the convection moved eastward to Gongliu County, it was coupled with the ascending phase of upper level waves, causing both the convection and precipitation to intensify again. We applied spectral analysis methods to verify that the waves were inertial gravity waves. The upper level inertial gravity waves propagated westward at a mean speed of −12 m s−1 with periods of 73–179 min and horizontal wavelengths of 50–55 km. The lower level inertial gravity waves propagated eastward at a mean speed of 8 m s−1 with periods of 73–200 min and a horizontal wavelength of 85 km. The more (less) favorable waveguide conditions determined whether the gravity waves persisted for a long (short) time and propagated for a longer (shorter) distance. Based on the mesoscale Eliassen–Palm flux theory, the wave energy of inertial gravity waves had an important effect on the maintenance and development of convection and precipitation by affecting wind strength and wind divergence. Feedback was mainly through the meridional and vertical transport of zonal momentum and the meridional transport of heat.

Tomographic Imaging of a Large Scale TID during the Halloween Storm of 2003

The most intense ionospheric storm observed in recent times occurred between 29–31 October 2003. The disturbances to the high-latitude regions set off several Large-Scale Travelling Ionospheric Disturbances (LSTIDs), wavelike perturbations in the ionospheric electron density. This paper investigates one particular Travelling Ionospheric Disturbance (TID) on 31 October 2003 using North American Global Positioning System (GPS) receiver network data and a tomographic imaging technique. The TID has an estimated period of 30 min, an estimated horizontal wavelength of 700 km and propagates South-Westward over North America. The tomographic reconstruction of the wave is validated using a simulation of the observations and with independent observations from ionosondes and the CHAMP Planar Langmuir Probe. The results are discussed in the context of the magnetic and ionospheric conditions that may have contributed to the launch of the wave. Large-scale TIDs are challenging to study over large regions of the Earth, and the GPS network here is shown to offer a unique perspective on the spatial and temporal variation of the TID. The experimental results are backed up by simulations that show a denser network of receivers, as is available in more recent years, would produce improved accuracy in the TID imaging.

Internal tide energy flux over a ridge measured by a co-located ocean glider and moored acoustic Doppler current profiler

Internal tide energy flux is an important diagnostic for the study of energy pathways in the ocean, from large-scale input by the surface tide to small-scale dissipation by turbulent mixing. Accurate calculation of energy flux requires repeated full-depth measurements of both potential density (ρ) and horizontal current velocity (u) over at least a tidal cycle and over several weeks to resolve the internal spring–neap cycle. Typically, these observations are made using full-depth oceanographic moorings that are vulnerable to being “fished out” by commercial trawlers when deployed on continental shelves and slopes. Here we test an alternative approach to minimize these risks, with u measured by a low-frequency acoustic Doppler current profiler (ADCP) moored near the seabed and ρ measured by an autonomous ocean glider holding station by the ADCP. The method is used to measure the semidiurnal internal tide radiating from the Wyville Thomson Ridge in the North Atlantic. The observed energy flux (4.2±0.2 kW m−1) compares favourably with historic observations and a previous numerical model study. Error in the energy flux calculation due to imperfect co-location of the glider and ADCP is estimated by subsampling potential density in an idealized internal tide field along pseudorandomly distributed glider paths. The error is considered acceptable (<10 %) if all the glider data are contained within a “watch circle” with a diameter smaller than 1∕8 the mode-1 horizontal wavelength of the internal tide. Energy flux is biased low because the glider samples density with a broad range of phase shifts, resulting in underestimation of vertical isopycnal displacement and available potential energy. The negative bias increases with increasing watch circle diameter. If watch circle diameter is larger than 1∕8 the mode-1 horizontal wavelength, the negative bias is more than 3 % and all realizations within the 95 % confidence interval are underestimates. Over the Wyville Thomson Ridge, where the semidiurnal mode-1 horizontal wavelength is ≈100 km and all the glider dives are within a 5 km diameter watch circle, the observed energy flux is estimated to have a negative bias of only 0.4 % and an error of less than 3 % at the 95 % confidence limit. With typical glider performance, we expect energy flux error due to imperfect co-location to be <10 % in most mid-latitude shelf slope regions.

Assessing the accuracy of SRTM altitude data for the hilly area in northeastern Romania

SRTM data is still one of the most used data in geosciences for various purposes: geomorphometric analysis, environmental covariate modelling or geomorphic change detection. Although high resolution national/regional DEMs exist, very often accessing them is expensive, or their coverage is not complete over specific areas (only floodplains or cities are covered). Because of this SRTM still remains the best choice when elevation data is needed for regional/national or global areas. In order to assess the correctness of SRTM data to depict the real shape of Earth surface we used a regional high resolution DEM which cover a part of the hilly area of north-eastern Romanian. Both DEMs were converted to the same horizontal and vertical datum (Stereo 70 Romanian projection and the EGG97 geoid), interpolated to the same grid size and position and compared using raster algebra. The horizontal x and y components and the vertical component errors were assessed. The results show that the errors of the SRTM model are well consistent with its acquisition method (the presence of the trees and the topographic shadow) and does represent reasonably well the Earth’s surface in the study area. Anyhow, the resolution of the Earth features depicted on the SRTM model is limited by the acquisition method and does not incorporate landforms which have a vertical and horizontal wavelength under 100 m.

Assessing the accuracy of SRTM altitude data for the hilly area in northeastern Romania

SRTM data is still one of the most used data in geosciences for various purposes: geomorphometric analysis, environmental covariate modelling or geomorphic change detection. Although high resolution national/regional DEMs exist, very often accessing them is expensive, or their coverage is not complete over specific areas (only floodplains or cities are covered). Because of this SRTM still remains the best choice when elevation data is needed for regional/national or global areas. In order to assess the correctness of SRTM data to depict the real shape of Earth surface we used a regional high resolution DEM which cover a part of the hilly area of north-eastern Romanian. Both DEMs were converted to the same horizontal and vertical datum (Stereo 70 Romanian projection and the EGG97 geoid), interpolated to the same grid size and position and compared using raster algebra. The horizontal x and y components and the vertical component errors were assessed. The results show that the errors of the SRTM model are well consistent with its acquisition method (the presence of the trees and the topographic shadow) and does represent reasonably well the Earth’s surface in the study area. Anyhow, the resolution of the Earth features depicted on the SRTM model is limited by the acquisition method and does not incorporate landforms which have a vertical and horizontal wavelength under 100 m.