the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Objective identification of pressure wave events from networks of 1 Hz, high-precision sensors
Matthew A. Miller
Laura M. Tomkins
Mesoscale pressure waves, including atmospheric gravity waves, outflow and frontal passages, and wake lows, are outputs of and can potentially modify clouds and precipitation. The vertical motions associated with these waves can modify the temperature and relative humidity of air parcels and thus yield potentially irreversible changes to the cloud and precipitation content of those parcels. A wavelet-based method for identifying and tracking these types of wave signals in time series data from networks of low-cost, high-precision (0.8 Pa noise floor, 1 Hz recording frequency) pressure sensors is demonstrated. Strong wavelet signals are identified using a wave-period-dependent (i.e., frequency-dependent) threshold, and then those signals are extracted by inverting the wavelet transform. Wave periods between 1 and 120 min were analyzed – a range which could capture acoustic, acoustic-gravity, and gravity wave modes. After extracting the signals from a network of pressure sensors, the cross-correlation function is used to estimate the time difference between the wave passage at each pressure sensor. From those time differences, the wave phase velocity vector is calculated using a least-squares fit. If the fitting error is sufficiently small (thresholds of RMSE < 90 s and NRMSE < 0.1 were used), then a wave event is considered robust and trackable. We present examples of tracked wave events, including a Lamb wave caused by the Hunga Tonga volcanic eruption in January 2020, a gravity wave train, an outflow boundary passage, a frontal passage, and a cold front passage. The data and processing techniques presented here can have research applications in wave climatology and testing associations between waves and atmospheric phenomena.
- Article
(11415 KB) - Full-text XML
- BibTeX
- EndNote
Gravity waves (i.e., buoyancy waves), which result from vertical perturbations of stably stratified fluid, are ubiquitous in the atmosphere and represent an important distributor of energy through the atmosphere (Fritts and Alexander, 2003; Nappo, 2013). The initial perturbations which generate gravity waves can have several sources, including but not limited to forced flow over topography, deep convection, shear instability, adjustment of unbalanced flow, and nonlinear interaction between waves (Fritts and Alexander, 2003). Within the troposphere, the vertical motions associated with gravity waves have been shown to influence cloud and precipitation processes. If parcels then reach saturation and produce precipitation which falls out, the changes to the parcel can be irreversible. Several studies have investigated the effects of gravity waves on marine stratocumulus. Allen et al. (2013) and Connolly et al. (2013) related gravity waves to changes in drizzle production within marine stratocumulus; specifically, enhanced condensation and collision–coalescence to form drizzle drops appeared to occur in the updrafts associated with gravity waves. Evidence also suggests a link between gravity waves and the rapid erosion of marine stratocumulus cloud decks (Yuter et al., 2018), perhaps because evaporation due to entrainment is enhanced as marine stratocumulus clouds are lifted by gravity waves (Connolly et al., 2013). Fovell et al. (2006) identified gravity waves as a potential trigger mechanism for deep moist convective cells. Case studies have linked gravity waves to snow bands, i.e., linear mesoscale enhancements in snowfall rate within winter storms (Bosart et al., 1998; Zhang et al., 2001; Gaffin et al., 2003), but it is unclear how often snow bands are associated with gravity waves.
To illustrate how the vertical motions associated with gravity waves could influence cloud microphysical properties, we use Fig. 1, which shows relative humidity with respect to ice (RHice) and with respect to liquid water (RHwater) as a function of temperature and water vapor mixing ratio for a standard atmosphere, as well as the temperature difference between 100 % RHice and 100 % RHwater for each water vapor mixing ratio. For temperatures below 0 ∘C, RHice>RHwater. A lifted parcel would be cooled at a constant water vapor mixing ratio (move upward in Fig. 1) until it intersects the 100 % RHice contour wherein vapor deposition would reduce the water vapor (further lifting would move upward and to the left in Fig. 1). If the parcel cools sufficiently to also intersect the 100 % RHwater contour, supercooled water droplets would form in the parcel and riming would likely occur, further depleting the available water vapor in the parcel. An up-and-down motion of an air parcel in a gravity wave that only crosses the 100 % RHice contour yields ice mass changes that are reversible; i.e., the ice mass added by vapor deposition in the upward motion will be lost to sublimation in the downward motion.
In contrast, if the lifting of the parcel within a gravity wave starts at RHice≥100 % and intersects the 100 % RHwater contour, the ice mass added by riming when ice particles collide with supercooled droplets is not reversible (i.e., there is no unriming process). For example, if an ice saturated parcel containing ice crystals (e.g., at −8.0 ∘C and 2.96 g kg−1 water vapor mixing ratio) is lifted and cooled to −9.7 ∘C in a gravity wave (which requires 190 m of lift assuming a 9 parcel lapse rate), it would become supersaturated with respect to water as well as with respect to ice. Water droplets would form in the parcel, and the ice crystals could potentially become rimed. If the parcel remains ice saturated when it descends in the gravity wave, the additional ice mass from the riming on the ice particle would not be removed. Ice mass added by riming can only be removed by sublimation in conditions where the parcel is subsaturated with respect to ice.
One way of distinguishing gravity waves from other wave phenomena such as Kelvin–Helmholtz waves is that gravity waves produce a surface pressure signal (Nappo, 2013, Sect. 8.2), given that the stable layer in which they occur is adjacent or nearly adjacent to the surface. That said, several different phenomena can produce surface pressure disturbances on similar spatiotemporal scales to gravity waves, including but not limited to outflow boundary passages, convective wake lows (Johnson and Hamilton, 1988), release of conditional symmetric instability (Gray et al., 2011), and Lamb waves generated by, e.g., distant volcanic eruptions (Matoza et al., 2022), large bolide impacts (ReVelle, 2008), and thermonuclear explosions (Pierce and Posey, 1971). Acoustic and acoustic-gravity waves can also produce pressure signals at shorter timescales. Such waves would include infrasound waves, which can be associated with, e.g., convective storms and strong flow over mountains (Coffer and Parker, 2022; Bedard, 1978).
Time series of surface pressure data have been analyzed to identify tropospheric wave events in previous studies (Kjelaas et al., 1974; Christie et al., 1978; Einaudi et al., 1989; Grivet-Talocia and Einaudi, 1998; Grivet-Talocia et al., 1999; Koch and Saleeby, 2001; de Groot-Hedlin et al., 2014). Kjelaas et al. (1974) and Christie et al. (1978) presented case studies of gravity wave events selected manually from time series pressure data. Grivet-Talocia and Einaudi (1998) and Grivet-Talocia et al. (1999) recorded data at Hz (i.e., every 2 min) to identify wave periods longer than 30 min. Einaudi et al. (1989) used a network of microbarographs recording at 0.1 Hz (i.e., every 10 s) placed within roughly 100 m of each other, which constrained the characteristics of disturbances which could be tracked through their network to waves with speeds up to 50 m s−1 and wave periods of 1–20 min. Koch and Saleeby (2001) used operational Automatic Surface Observing Systems (ASOS) data recorded at Hz (i.e., every 5 min) which resolved wavelengths ≥150 km. While strong pressure disturbances including outflows and wake lows can be detected by ASOS pressure sensors logging data at 1 min intervals, the large spatial separation between operational weather stations, which are mostly located at airports, precludes determination of the associated wave speed and direction for mesoscale disturbances. de Groot-Hedlin et al. (2014) used 337 barometers deployed with the USArray Transportable Array, recording at 1 Hz (i.e., every 1 s) frequency, to detect and track high-amplitude (roughly 3 hPa peak to trough) pressure waves associated with convective storms in the southern United States. The USArray Transportable Array barometers were spaced roughly 70 km apart, which might also preclude tracking of localized disturbances.
There is a scarcity of data for detecting and tracking pressure disturbances, including gravity waves, on the meso-β scale or meso-γ scale (ranging from 2 to 200 km). Pressure disturbances on those scales may be relevant to phenomena such as snow bands (e.g., McMurdie et al., 2022), trade wind cumulus (e.g., Seifert and Heus, 2013), and bow echoes (e.g., Adams-Selin and Johnson, 2010), which are active areas of research. We developed a measurement and analysis technique which will allow for questions regarding gravity waves on those scales to be addressed. To what degree and in what conditions information on gravity wave occurrence would be valuable in operational weather settings is yet to be determined.
This paper presents data from networks of internet appliance low-cost, high-precision air pressure sensors (i.e., microbarographs) and a methodology for objectively identifying mesoscale wave events and wave speed and direction. The methodology is intended to be used for post-processing in research applications rather than for real-time or near-real-time detection of wave events. Similar low-cost sensor networks have been used for detection of seismic waves (Anthony et al., 2018) and for detecting infrasound waves to monitor fan rotation speeds in nuclear-reactor cooling towers (Eaton et al., 2022). The former network covered the area of Oklahoma, and the latter used networks covering the area of a single nuclear reactor. Our networks of pressure sensors are on the scale of a medium-to-large-sized city or metropolitan area.
Section 2 of this paper describes the pressure sensors used in this study and the data they provide. Section 3 describes the methodology for objectively identifying pressure waves from the pressure time traces. Section 4 provides five examples of events captured by the wave identification method. Finally, a summary and avenues for future work are discussed in Sect. 5.
2.1 Networks of pressure sensors
Pressure sensors were placed in three separate networks: New York City metro area and Long Island, NY; Raleigh, NC; and Toronto, ON, Canada (Fig. 2). Each pressure sensor was either a Bosch BME280 (Bosch, 2022) or a Bosch BMP388 (Bosch, 2020) Adafruit breakout board connected to a Raspberry Pi Zero W single-board computer used to log the data. BME280 sensors measure pressure, temperature, and humidity. BMP388 sensors measure pressure and temperature. Each pressure sensor costs roughly USD 50–75, subject to changes in the cost of Raspberry Pi Zero W units. For context to another low-cost network concept, the Raspberry Shake 4D seismographs cost a few hundred dollars per unit (Anthony et al., 2018). The combined sensor and communications package is about the size of a deck of cards. Sensors are connected to the internet and sync their data to a server at North Carolina State University. Initial testing of the sensors outdoors on patios, in sheltered locations such as garages, and indoors revealed pressure waves were well resolved in all locations, and it was best practice to place the sensors indoors to minimize wind contamination in pressure measurements. When active, each sensor records pressure at 1 s intervals with a roughly 0.8 Pa noise floor depending on ambient conditions. The sensors synchronize to network time upon startup. The wave extraction method only depends on relative pressure variations and is not sensitive to absolute or relative calibration.
To examine the properties of gravity waves which are detectable by these pressure sensors, we consider an internal gravity wave occurring in an environment with constant background wind u0. The relationship between the pressure perturbation p′ and the horizontal velocity perturbation u′ associated with the wave is described by Nappo (2013):
where ρ0 is the environmental air density and c is the phase speed of the gravity wave. From Eq. (1), the maximum pressure perturbation can therefore be related to the maximum horizontal velocity perturbation by
Figure 3 shows values according to Eq. (2) at an air density of 1.225 kg m−3 (standard air density at sea level; American Meteorological Society, 2022) for and values up to 15 m s−1. In order to reliably detect a wave signal, the amplitude likely needs to be substantially larger than the noise floor.
2.2 Operational weather observations
For context, we compare extracted wave signals with available operational weather observations.
We use Automated Surface Observing Systems (ASOS; NOAA National Centers for Environmental Information, 2021a) data including surface temperature, dew point, wind speed and direction, and additional pressure measurements coincident with wave events. ASOS data are recorded each minute. For wave events detected in New York and Long Island, we examined ASOS data from John F. Kennedy International Airport (KJFK) and Long Island MacArthur Airport (KISP). For wave events detected in Toronto, we examined ASOS data from Buffalo Niagara International Airport (KBUF) and Niagara Falls International Airport (KIAG).
For one example case in Sect. 4.2, we show upper-air radiosonde data from a weather balloon launched in Buffalo, NY. We obtained data from the Integrated Global Radiosonde Archive (IGRA; NOAA National Centers for Environmental Information, 2021b) and interpolated to a constant 100 m resolution. The data include measurements of temperature, dew point, and winds, from which we calculated wet-bulb temperature, frost point, and saturation equivalent potential temperature (). Because radiosondes are typically launched every 12 h at a limited number of locations, representative radiosonde data are not available for every case.
We use horizontal maps of data from the US National Weather Service (NWS) WSR-88D radars (NOAA National Weather Service Radar Operations Center, 1991) to show storm features occurring coincident with wave events in Sects. 4.3 and 4.5. Radar reflectivities are processed following Tomkins et al. (2022) to indicate regions with mixed precipitation in the scan, by inferring that points with reflectivity above 20 dBZ and dual-polarization correlation coefficient below 0.97 have mixed precipitation. In maps of radar reflectivity, those regions with mixed precipitation are then shown in greyscale. Doppler velocity waves are extracted from radial velocity data following Miller et al. (2022), by calculating the difference in radial velocity from successive scans, converting those differences to a binary (positive/negative) field, and filtering out small objects in that binary field.
The methods outlined here for identifying wave events in the pressure time traces are adapted from the techniques used by Grivet-Talocia and Einaudi (1998) and Grivet-Talocia et al. (1999). The method uses wavelet transforms to identify wave events in time-wave period (or, equivalently, time-frequency) space. Wavelet transforms are preferable to Fourier transforms for the purpose of identifying transient waves which are localized in time (Torrence and Compo, 1998). To illustrate the step-by-step procedure, an example corresponding to a gravity wave event on 23 February 2023 in the Toronto pressure network is described in detail.
3.1 Identifying wave events in a single sensor
The full pressure time series for a gravity wave event on 23 February 2023 captured by sensor 25 in Toronto is shown in Fig. 4a. As an initial pre-processing step, 10 s samples of pressure (i.e., averages of 10 pressure measurements) are used to smooth out noise and pressure perturbations due to high-frequency turbulent eddies in the data (Fig. 4b). Hereafter, time series labeled as total pressure are the 10 s subsamples of the original pressure measurements.
A wavelet transform W of a finite energy signal f(t) (a pressure time series in this study) can be defined as in Grivet-Talocia and Einaudi (1998, their Eq. 1):
where a is the scale (related to the wave period), and b shifts the wavelet in time (t). ψ* represents the mother wavelet. An analytic Morse wavelet was used (e.g., Olhede and Walden, 2002; Lilly and Olhede, 2012) via the cwt function within MATLAB (Lilly, 2021). In this study, W(b,a) will always refer to the wavelet transform of a pressure time series. The resulting wavelet transform is an array of complex values in timescale space. The absolute value of the wavelet transform can be considered the wavelet power at a given time and scale. Figure 4c shows the wavelet power associated with the wave event at sensor 25 on 23 February 2023. In this study, wave periods between 1 and 120 min were analyzed, corresponding to expected periods for mesoscale disturbances.
To objectively identify wave event centers according to wavelet power, a scale-dependent (i.e., wave-period-dependent) threshold function A(a) is defined as the mean wavelet power across all available data for the sensor network by scale, multiplied by a constant K:
A scale-dependent threshold K is necessary because the “background” wavelet power for a pressure time series generally increases with scale (e.g., Canavero and Einaudi, 1987). Grivet-Talocia and Einaudi (1998) and Grivet-Talocia et al. (1999) used 2 as an appropriate value for K. Lower values of K lead to more wave events being detected, which can include potential artifacts in the pressure time trace. For the present study, a K value of 10 was used to ensure that only the strongest wave signals were identified (solid contour in Fig. 4d). This threshold can be adjusted, and different applications may warrant different values of K. The mean wavelet power as a function of wave period is shown for each regional sensor network and for all networks combined in Fig. 5. Event centers were identified as local maxima in wavelet amplitude which exceed A(a) and which are located at (bmax,amax). In Fig. 4d, an event center is located within the solid contour. From the identified event centers, the first iterations of event regions (Ω′) were identified in timescale space as connected regions where the wavelet power exceeds , i.e., half of the event center threshold. In Fig. 4d, Ω′ is represented by the region within a dashed contour which contains a solid contour.
The watershed transform (Meyer, 1994) was used to refine Ω′. Watersheds (i.e., catchment basins) were identified in the negative wavelet power array . Any watersheds within Ω′ whose period range was entirely outside the period range of the watershed containing the event center were removed from the event region Ω′. This step was included to correct cases where multiple “peaks” in wavelet power were present within Ω′ at different wave periods, with a “valley” in wavelet power in between where the wavelet power still exceeded , which likely represented distinct wave modes and should be considered separate wave events.
Ω′ was extended to define the final event region Ω for each event, first by taking the bounding box of Ω′ and then by extending the bounding box along the time axis in both directions until it reaches a local minimum in the wavelet amplitude to ensure that the entire signal of interest is contained in the event region. This could result in overlapping event regions. Figure 4d shows the wavelet power normalized by the mean wavelet power by scale () for the 23 February 2023 example in sensor 25, with the outline of the event region overlaid with the magenta box.
After defining the event region Ω, the wave event trace could then be extracted (i.e., reconstructed) by inverting the wavelet transform function over the event region Ω. Figure 4e shows the extracted wave event trace for the 23 February 2023 wave event in sensor 25. As in Grivet-Talocia and Einaudi (1998), wave events were identified and extracted one at a time, with the extracted wave event subtracted from the pressure trace and the wavelet transform recalculated at each iteration, until the absolute maximum of was less than K (i.e., until no more events are left to be found in the pressure time series).
We tested the method of detecting wave events in a single sensor using synthetic pressure data. The synthetic time series of pressure is created by an initial constant pressure value (which is randomly chosen from a normal distribution with mean 1000 hPa and standard deviation 2 hPa). We then added normally distributed random noise centered on 0 with standard deviation equal to the noise floor (0.008 hPa). Finally, we added 105 pre-defined wave events, which consist of sine waves of period ranging from 2 to 120 min and maximum amplitude ranging from ±0.01 to ±1 hPa. Each set of sine waves lasts for 2 h, with a 12 min ramp-up and ramp-down period at the start and end of those 2 h in which the amplitude increases and decreases linearly, respectively. One of these synthetic wave events is shown in Fig. 6a. Using the values shown in Fig. 5 and a K value of 10, 52 of the 105 synthetic wave events were detected, with no false-positive event detections (Fig. 6c). The weakest detected synthetic wave event had an amplitude of ±0.0464 hPa (or 0.0928 from peak to trough) and a wave period of 2 min (shown in Fig. 6a). Lower values of K do lead to more wave events being detected, with few false positives. The K=2 used by Grivet-Talocia and Einaudi (1998) results in 86 of the 105 wave events being detected, with only one false positive. However, this exercise likely fails to capture the full extent of noise and the interference of many signals present in real pressure data. Lower values of K can result in more weak pressure wave events being detected, which may be “real” at a single sensor location, but these may then be erroneously paired with other weak pressure wave events at other sensor locations using the methods described in the following sections. Therefore, we will use K=10 to detect and track pressure wave events across the networks of sensors.
3.2 Matching corresponding wave events between multiple sensors
Once wave events were identified for each sensor individually, the following steps were taken to identify coherent wave events across multiple sensors. For this purpose, the terms primary sensor (denoted by i) and secondary sensor (denoted by j) will be used to describe a pair of sensors for which events are identified and paired together.
For each event in the primary sensor pressure trace, events in the secondary sensor pressure trace that occurred within 2 h of that primary sensor event (i.e., with a gap between the end of the event in one sensor and the start of the event in the other sensor not exceeding 2 h) were considered “candidate” events to match with the primary sensor event. This 2 h threshold is subjective, and it affects the range of speeds of wave events which can be detected. The threshold can be altered depending on the distances between pressure sensors and the desired application. For example, the largest distance between any two sensors in the three networks is that between sensors 012 and 027 (136 km apart, Fig. 2). A wave feature propagating at 18.9 m s−1 would take 2 h to propagate that distance. However, wave features propagating at an angle could be slower and propagate over both sensors within 2 h, and as long as pressure sensors in between those two most distant sensors capture the event, the following processing technique will allow those sensors to bridge the gap even if the event takes longer than 2 h to propagate across the distance between those sensors. Thus, 18.9 m s−1 is a conservative estimate of the minimum phase speed required for this methodology to track a wave event. Candidate matching events in the secondary sensor trace had to have a center period which was within the primary sensor event period range, and vice versa. Then, for each candidate matching event in the secondary sensor trace, the event waveforms are reconstructed by inverting the wavelet transform over the event region for both sensors. Figure 7 shows the extracted waveform for the 23 February 2023 event in sensor 25 and the same wave passage in sensors 04, 23, 24, and 34. The time lag estimate for the wave passage between sensors is Δtopt, the time lag which maximizes the cross-correlation function Cij(Δt):
where pi(t) and pj(t) are the extracted waveforms for the events in the primary and secondary sensor, respectively. The black lines and subfigure titles in Fig. 7 show the optimal shift in the extracted waveforms for sensors 04, 23, 24, and 34 to maximize Cij to sensor 25 for the 23 February 2023 example. The match to the primary sensor event is the candidate event with the highest maximized cross-correlation to the primary sensor event. If the maximized cross-correlation exceeded 0.65 and the same pair of matched events results from switching the primary and secondary sensors (i.e., the event is matched two ways), the event from sensor i and the event from sensor j are paired together. Switching the primary and secondary sensors is necessary to avoid instances where multiple events in one sensor are matched with the same event in another sensor. This can occur, for example, when a set of waves manifests as one event in one sensor and multiple (separate) events in another sensor.
The process of matching events between sensors described above was repeated for each possible combination (within a sensor network) of primary and secondary sensors in order to obtain the full set of lag times between each pair of sensors which captured each event. In other words, N2 pairs of sensors, order-dependent, were analyzed, where N is the number of sensors in the network with data at a given time. Then, each event in each sensor was assigned an ID based on which other sensors had a matching event in order to track events across three or more sensors. This process required iterating through each sensor in a network. Each event in the first sensor was assigned a new (i.e., arbitrary) ID. For each subsequent sensor sc, events with no two-way matches in any prior sensor were also given new IDs. If there were two-way matches with an event in one or more prior sensor(s), the event in the sensor sc would share the ID assigned to the matched event in the prior sensor. If there were multiple prior sensors with matched events and those events had different IDs, the ID associated with the higher maximized cross-correlation between the event traces was assigned to the event in sensor sc. If this process results in multiple events in sensor sc sharing the same event ID D, the event in sensor sc associated with the highest maximized cross-correlation with any one prior sensor for an event with event ID D is assigned event ID D, and the previously outlined steps are repeated for the other event(s) in sensor sc, with event ID D and associated sensors excluded.
The result of this process is a set of events with associated ID numbers for every sensor in the network. For a single sensor, each event has a unique ID. For each ID number that appeared in at least three sensors, the wave phase velocity vector was calculated using the set of lag times between each pair of sensors which captured the event.
3.3 Estimating wave phase velocity vector
Once sets of matched events were identified, the wave propagation velocities (two-dimensional vectors) could be estimated for events which occurred in three or more sensors. It is hypothesized for each wave event that a plane wave crosses the sensor network with slowness vector , where sx and sy are the inverses of the x and y components of the wave propagation vector (in s m−1), respectively. s can be solved for from the following equation (Del Pezzo and Giudicepietro, 2002):
where t is the column vector of the Δtopt values for each possible pair of sensors which captured the event, and Δx is the two-column matrix of the x and y components of the distance vector between each pair of sensors which captured the event. t and Δx each have rows, where Ns is the number of sensors which captured the event. Equation (6) can be considered an overdetermined system of linear equations, as long as Ns≥3, and is solved for s by a least-squares approach represented by
where superscript T indicates the transpose of a matrix (Del Pezzo and Giudicepietro, 2002).
Once sx and sy are solved for, they can be inverted to obtain the wave phase velocity components, cx and cy, respectively. Additionally, the modeled delay times tm can be calculated by solving Eq. (6) for t. From tm, we estimate the model error using root mean square error (RMSE) and normalized root mean square error (NRMSE):
Events with sufficiently small RMSE and NRMSE in the modeled delay times can be considered “trackable” events in that there is higher confidence in the wave velocity estimates for those events. After processing multiple years of data from the Toronto and New York pressure sensor networks and analyzing the resulting RMSE and NRMSE distributions, it was found that a maximum RMSE threshold of 90 s and a maximum NRMSE threshold of 0.1 are reasonable to consider a wave event trackable.
Additionally, we require that wave events be captured by at least four sensors to be considered trackable. While the slowness vector and corresponding error metrics can be calculated for events captured by only three sensors, the calculation is less constrained because there are only 3 delay times Δtopt in the calculation (compared to 6 delay times for events captured by four sensors, 10 delay times for events captured by five sensors, etc.). The result is that events captured by only three sensors can have small RMSE and NRMSE by chance much more easily than events captured by four or more sensors. For each robust and trackable wave event, the mean amplitude was calculated by averaging the difference between the maximum and minimum values in the extracted event trace for each sensor which captured the event. The center period for the event was calculated as the mean of the wave period corresponding to bmax for each sensor which captured the event.
It is especially useful to have examples of wave events where other observational sources constrain the wave phase speed and direction. We discuss the following examples in this section, some of which have corroborating information on the wave phase speed and/or direction. The Lamb waves caused by the Hunga Tonga–Hunga Ha'apai volcanic eruption in January 2022 represent a case where the origin of the waves is known and the phase speed is known as a function of air temperature (Amores et al., 2022). A gravity wave train which passed over Toronto on 25 February 2022 occurred coincident with a surface cyclone 100 km distant but in local conditions of sparse radar echo. A wave event on 4 February 2022 is a case associated with an outflow boundary clearly captured by Doppler radar data from the WSR-88D radar located in Upton, NY. In another example, waves coincided with a cold front which passed over Toronto on 15 November 2020. The cold front's associated narrow rain band and Doppler velocity wave (Miller et al., 2022) can be identified in WSR-88D radar data from Buffalo, NY. We also describe a wake low associated with a long-lived mesoscale convective system (MCS) which passed over Long Island on 14 September 2021. These events were each manually chosen after roughly 40 months of pressure data were processed. The gravity wave train on 25 February 2022 and cold front example on 15 November 2020 are not unusual; several other gravity wave trains and pressure jumps due to cold front passages were found. The other three example cases are atypical events, however. It is extremely rare to detect a pressure signal due to a volcanic eruption several thousand kilometers away. While other outflow boundaries were found to produce detectable pressure waves, the case on 4 February 2022 was unusual in terms of the time of year when it occurred (most other pressure waves associated with outflow boundary passages occurred during the warm season) and the radar signature. The wake low detected on 14 September 2021 was the only wake low event detected by our pressure sensor network.
4.1 15–17 January 2022: Hunga Tonga–Hunga Ha'apai eruption and Lamb wave
On 14–15 January 2022, large eruptions occurred at the Hunga Tonga volcano in the south Pacific Ocean, which produced ash plumes reaching the mesosphere and a series of shock waves. A particularly violent, submarine eruption occurred at around 04:00 UTC on 15 January (Global Volcanism Program, 2022). Satellite data suggest that the ash plume associated with this eruption reached as high as 57 (Carr et al., 2022; Proud et al., 2022) and contained roughly 400×106 kg of sulfur dioxide. Damaging tsunami waves due to the eruption were observed as far as Peru (Global Volcanism Program, 2022). Subsequent analyses of the atmospheric pressure waves from the eruption have classified the pressure wave observed far from the source eruption as a Lamb wave (e.g., Amores et al., 2022). This Lamb wave has been the subject of several studies and media reports since the time of the eruption (e.g., Amores et al., 2022; Adam, 2022; Burt, 2022; Bhatia and Fountain, 2022; Vergoz et al., 2022). The pressure signal associated with the Lamb wave was observed to circle the Earth several times with estimated phase speeds exceeding 100 m s−1 (Adam, 2022; Burt, 2022; Vergoz et al., 2022).
To identify and characterize the pressure waves from this event, we combined the three regional sensor networks (Fig. 2) to effectively create an array of 18 sensors which were active at the time of the Lamb wave passages. Table 1 summarizes events identified during this period which meet the robust event criteria outlined in the methods (captured by at least four sensors, RMSE≤90 s, and NRMSE≤0.1), in addition to having a mean optimal cross-correlation between the extracted event traces exceeding 0.75. The initial outbound (traveling from Tonga to the antipode location in Algeria) waves manifested as three separate detected events between 15:09 and 17:14 UTC on 15 January, each with high wave frequencies (i.e., low wave periods barely over 1 min) and low amplitudes (up to 0.3 hPa). The earliest and strongest of these events had a phase speed of 326.6 m s−1 and direction of 64.2∘ (i.e., to the east-northeast). The subsequent rebound (traveling from the antipode location to Tonga) waves were detected as two separate events between 03:56 and 07:29 UTC on 16 January. The rebound waves had a much higher amplitude (roughly 2.4 hPa) and lower wave frequency (i.e., longer wave periods on the order of roughly 10 min) than the outbound waves (Fig. 8). The earliest rebound wave event had a phase speed of 292.1 m s−1 and direction of 261.9∘ (i.e., to the west-southwest). Our networks captured the initial outbound waves from Tonga to the antipode at Algeria (Outbound (i), (ii), and (iii) in Table 1) and the first rebound waves back from Algeria (Rebound (i) and (ii) in Table 1). Subsequent reverberations of the Lamb waves were not trackable with our sensors and methods by the above criteria due to a combination of the low amplitudes, high frequency, and large phase speeds confounding the process of approximating delay times between sensors (an issue also described by Grivet-Talocia and Einaudi, 1998).
4.2 25 February 2022: gravity wave train over Toronto
Four pressure sensors in Toronto were used to detect and track a series of pressure oscillations between 08:16 and 13:59 UTC on 25 February 2022, with a particularly large pressure peak near 10:00 UTC (Fig. 9). The mean amplitude of the event across the four sensors was 2.1 hPa, and the wave train was estimated to propagate at 45.4 m s−1 at 73∘ (i.e., to the east-northeast). The center wave period was 00:20:22. At this time a mature surface cyclone was located roughly 100 km to the south of Toronto. Linear bands of reflectivity were identified in WSR-88D radar data from Buffalo, NY, in the hours leading up to the detected pressure waves, but between 09:00 and 12:00 UTC there was only sparse radar echo over the Toronto area. Between 12:30 and 14:30 UTC there was radar echo across the Toronto area, and a set of Doppler velocity waves was identified from the WSR-88D data using the methods in Miller et al. (2022). Those Doppler velocity waves appeared to propagate toward the northeast, in roughly the same direction as the detected pressure waves (“Video supplement” Animation Figure S01).
The nearest available sounding during the duration of the wave event appeared to indicate adequate environmental conditions for gravity wave ducting (Lindzen and Tung, 1976; Koch and O'Handley, 1997). The sounding from Buffalo, NY, valid at 12:00 UTC on 25 February 2022, shows a temperature inversion roughly 2 km deep (Fig. 10), which serves as the “ducting layer” directly above a shallow boundary layer. A moist neutral or conditionally unstable layer (indicated by near-zero or negative values of the vertical gradient in equivalent saturation potential temperature) was above this inversion extending to around 4500 m a.m.s.l. and serves as the “trapping layer” in Fig. 10. The sharp change in stability between the inversion layer and conditionally unstable layer at roughly 2800 m a.m.s.l. could serve as a reflector of gravity wave energy. The apparent presence of a gravity wave duct during the detected pressure wave event raises confidence that the pressure waves were gravity waves.
4.3 4 February 2022: outflow pressure jump and subsequent oscillations over Long Island
Between 17:30 and 19:00 UTC on 4 February 2022, an event with amplitude of roughly 1.8 hPa was detected by five pressure sensors in the New York City metro area and Long Island. This event was a positive jump in pressure followed, to varying degrees in each sensor, by weak oscillations in the pressure trace (Fig. 11). Prior to the jump in pressure, there had been a decreasing trend in the pressure traces for several hours. At the same time, WSR-88D weather radar data from Upton, NY, showed widespread precipitation echo over Long Island. A wave feature was apparent in the Doppler radial velocity data, which could also be identified following the methods of Miller et al. (2022) (Fig. 12 and “Video supplement” Animation Figure S02). This wave event had a phase speed of 21.1 m s−1 and direction of 118.2∘ clockwise from north (i.e., southeastward). The values are consistent with the radar-detected Doppler velocity wave feature (“Video supplement” Animation Figure S02).
Operational 1 min Automated Surface Observing System (ASOS) data (Fig. 13a) also recorded a jump in the surface pressure of nearly 2 hPa. Near the time of this jump, there was also a peak in the wind speed and gusts, along with a brief shift in the wind direction from north-northeasterly to north-northwesterly (Fig. 13b). These features, along with the modest decrease in the temperature (Fig. 13a), are consistent with a convective outflow boundary (i.e., gust front). A fine line can be seen in WSR-88D reflectivity data at roughly the same location as the wave, which further suggests that a convective outflow was responsible for the pressure rise (Fig. 12 and “Video supplement” Animation Figure S02).
4.4 15 November 2020: cold front passage over Toronto
A robust and trackable wave event was detected by five pressure sensors in Toronto coincident with a cold front passage at roughly 20:00 UTC on 15 November 2020. The pressure steadily dropped in the hours leading up to the frontal passage before abruptly rising 1–2 hPa as the cold air mass arrived. The pressure then dropped roughly 1 hPa about 30 min later. Some sensors recorded oscillations in the pressure trace embedded within the gradual pressure rise in the following hours (Fig. 14). 1 min ASOS data from Buffalo, NY, also captured the pressure jump at roughly the same time as the temperature and dew point drop indicating the cold front passage (Fig. 15).
The pressure wave event for the cold front passage had an estimated phase speed of 27.5 m s−1 at 65∘ (i.e., to the east-northeast), a mean amplitude of 1.8 hPa, and a center wave period of 00:02:08. WSR-88D radar data from Buffalo, NY, show a narrow band of high reflectivity and a Doppler velocity wave (identified following the methods in Miller et al., 2022) associated with the cold front advancing over Toronto at roughly 20:00 UTC at a speed and direction consistent with the pressure wave (“Video supplement” Animation Figure S03).
4.5 14 September 2021: wake low associated with a mesoscale convective system
Between 03:00 and 04:00 UTC on 14 September 2021, a pressure drop of roughly 5 hPa and subsequent recovery occurred at four of the pressure sensors in the New York City metro area and Long Island network (Fig. 16). This was detected as a wave event with an estimated propagation speed of 20.6 m s−1 and propagation direction of 67.5∘ (i.e., to the east-northeast). ASOS data from Islip, NY (KISP; Fig. 17) and other stations in the area (not shown) also recorded the pressure minimum. This wave event occurred near the time of a mesoscale convective system (MCS) passage over Long Island as indicated by reflectivity data from the WSR-88D radar in Upton, NY (Fig. 18 and “Video supplement” Animation Figure S04). In addition to the precipitation echo associated with the MCS translating from northwest to southeast, there was also a stationary region of weak echo with low dual-polarization correlation coefficient (shown in greyscale following Tomkins et al., 2022) in the vicinity of the radar. The stationary weak echo was likely non-meteorological and due to either birds or insects. The precipitation echo appears to be entirely past KISP by 03:42 UTC (Fig. 18c and “Video supplement” Animation Figure S04), which is roughly the same time as the minimum pressure at KISP (Fig. 17a).
This pressure minimum appears to be consistent with a wake low, associated with subsidence heating in the rear inflow jet (Markowski and Richardson, 2010; Johnson and Hamilton, 1988). The subsidence heating does not necessarily lead to warming at the surface, which was not observed in the ASOS data (Fig. 17a), but decreased air density aloft due to warming will still lead to a surface pressure decrease. Markowski and Richardson (2010) also note that a property of wake lows associated with a translating squall line is that the center of convergence due to the wake low does not perfectly align with the center of the wake low. Rather, the convergence center slightly lags behind the wake low center. In the 14 September 2021 example, the ASOS time series data have a wind speed minimum co-occurring with a shift in the wind direction from near 100∘ (east-southeasterly) to near 280∘ (west-northwesterly), which can be interpreted as the convergence maximum (Fig. 17b). This convergence maximum occurs slightly after the pressure minimum associated with the wake low (Fig. 17a), consistent with the Markowski and Richardson (2010) description.
In this study, a wavelet-based method was used to identify wave events in time series pressure data from networks of high-precision sensors (0.8 Pa noise floor) recording the pressure every second. In addition to identifying wave events in each sensor individually, the delay times in wave passage among sensors in a given network were used to determine the direction of wave propagation and phase velocity. The methods shown are intended mainly for post-processing of pressure data for research applications and not for real-time, operational use. A benefit to this method is that it can be fully automated to detect wave events across many months of data, and we have made the processing code openly available (Allen, 2023).
Overall, the method was most successful at tracking pressure wave events with relatively large amplitudes (on the order of 0.3 hPa or more) and longer periods (i.e., lower frequencies; on the order of 5 min or more). Low-amplitude, high-frequency waves likely propagated across the sensor networks many times, but these waves were difficult to reconstruct, due to their wavelet signal being weaker, and to track, due to aliasing of the waveform conflating the time lag estimates between sensors. We use a rather strict criteria for detecting a wave event in a given sensor; the peak wavelet power must exceed 10 times the mean value for a given wave period (i.e., K=10; Eq. 4). Grivet-Talocia and Einaudi (1998) also used wavelet analysis to detect pressure waves with a scale-dependent threshold; their K value was only 2. Therefore, many wave events detected using our criteria will be relatively high amplitude (on the order of 1 hPa or more), including most examples shown in Sect. 4. Depending on the desired application, this threshold and the other thresholds used in the wave detection can be adjusted.
Environmental factors can influence whether or not a given gravity wave is detected by our surface pressure sensors. Gravity waves aloft will not always produce a detectable pressure signal at the surface, for example if the planetary boundary layer is neutral or unstable (e.g., Kjelaas et al., 1974). Another possible limitation is that in their current network deployments the pressure sensors are too far apart to track highly localized disturbances, particularly for the New York City and Long Island sensor network. Our method may not always properly track waves which are modified by local conditions (which may alter their amplitude, frequency, and/or phase velocity) as they propagate across the sensor network. Future work will examine data from networks of pressure sensors a few kilometers to a few meters apart and the degree to which signals associated with waves in shallow marine clouds are detectable with these sensors.
Deployment of networks of low-cost, high-precision sensors opens myriad opportunities for monitoring the direction and speed of gravity waves that have not been previously available with conventional pressure sensors on operational weather stations due to their longer measurement interval and larger station spacing. A forthcoming publication will describe a 3+ year climatology of wave events detected by the pressure sensors deployed in New York and Toronto and address hypotheses regarding the relationship between gravity waves and local enhancements in snowfall rate within winter storms (i.e., snow bands). There are observational case studies demonstrating this connection (e.g., Bosart et al., 1998; Gaffin et al., 2003), but a multi-year data set with continuously monitoring pressure sensors in the context of radar data will enable a more comprehensive examination of the co-occurrence, or lack thereof, of gravity waves with snow bands across many winter storms.
The pressure time series data used throughout this publication can be found at https://doi.org/10.5281/zenodo.8136536 (Miller and Allen, 2023). The NWS NEXRAD Level II data used in Figs. 12 and 18 can be accessed from the National Centers for Environmental Information (NCEI) at https://doi.org/10.7289/V5W9574V (NOAA National Weather Service Radar Operations Center, 1991). The NWS ASOS surface station data used to create Figs. 13 and 17 can be accessed from NCEI at https://www.ncei.noaa.gov/products/land-based-station/automated-surface-weather-observing-systems (NOAA National Centers for Environmental Information, 2021a). The radiosonde data used to create Fig. 10 can be accessed from NCEI at https://www.ncei.noaa.gov/access/metadata/landing-page/bin/iso?id=gov.noaa.ncdc:C01500 (NOAA National Centers for Environmental Information, 2021b).
The code used for processing the pressure time series data can be found at https://doi.org/10.5281/zenodo.10150876 (Allen, 2023).
List of animations with captions and filenames. All animations can be viewed at https://doi.org/10.5446/s_1476 (Allen et al., 2023a). Individual animations can be viewed by following the DOI URL.
Animation Figure S01. Animated maps of reflectivity and Doppler velocity wave detection for NWS WSR-88D radar data from Buffalo, NY, at 0.5∘ tilt, from 07:06 to 14:57 UTC on 25 February 2022. Reflectivity values are shown in greyscale when there is likely enhancement due to melting (Tomkins et al., 2022). Filled blue circles indicate locations of pressure sensors which captured the wave event described in Sect. 4.2, and unfilled blue circles indicate locations of pressure sensors which did not capture the wave event. Title: 25 Feb 2022 KBUF Reflectivity and Doppler Velocity Waves. https://doi.org/10.5446/62539 (Allen et al., 2023e).
Animation Figure S02. Animated maps of reflectivity and Doppler velocity wave detection for NWS WSR-88D radar data from Upton, NY, at 0.5∘ tilt, from 15:41 to 21:29 UTC on 4 February 2022. Reflectivity values are shown in greyscale when there is likely enhancement due to melting (Tomkins et al., 2022). Filled blue circles indicate locations of pressure sensors which captured the wave event described in Sect. 4.3, and unfilled blue circles indicate locations of pressure sensors which did not capture the wave event. This animation goes with Fig. 12. Title: 04 Feb 2022 KOKX Reflectivity and Doppler Velocity Waves. https://doi.org/10.5446/62540 (Allen et al., 2023d).
Animation Figure S03. Animated maps of reflectivity and Doppler velocity wave detection for NWS WSR-88D radar data from Buffalo, NY, at 0.5∘ tilt, from 18:05 to 23:24 UTC on 15 November 2020. Reflectivity values are shown in greyscale when there is likely enhancement due to melting (Tomkins et al., 2022). Filled blue circles indicate locations of pressure sensors which captured the wave event described in Sect. 4.4, and unfilled blue circles indicate locations of pressure sensors which did not capture the wave event. Title: 15 Nov 2020 KBUF Reflectivity and Doppler Velocity Waves. https://doi.org/10.5446/62541 (Allen et al., 2023c).
Animation Figure S04. Animated maps of reflectivity and Doppler velocity wave detection for NWS WSR-88D radar data from Upton, NY, at 0.5∘ tilt, from 00:03 to 07:27 UTC on 14 September 2021. Reflectivity values are shown in greyscale when there is likely enhancement due to melting (Tomkins et al., 2022). Filled blue circles indicate locations of pressure sensors which captured the wave event described in Sect. 4.5, and unfilled blue circles indicate locations of pressure sensors which did not capture the wave event. This animation goes with Fig. 18. Title: 14 Sep 2021 KOKX Reflectivity and Doppler Velocity Waves. https://doi.org/10.5446/62542 (Allen et al., 2023b).
LRA, SEY, and MAM conceptualized the project and designed the methodology. MAM designed and built the pressure sensors and managed the pressure sensor networks. LRA and MAM wrote the software. LRA and LMT created the visualizations with input from SEY. LRA prepared the manuscript, SEY edited the manuscript, and MAM and LMT contributed to the final stages of reviewing and editing.
The contact author has declared that none of the authors has any competing interests.
Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this paper. While Copernicus Publications makes every effort to include appropriate place names, the final responsibility lies with the authors.
The authors express their sincere appreciation to the pressure sensor hosts in Toronto, the New York City metro area and Long Island, and Raleigh, including colleagues from Environment Canada, Stony Brook University, Columbia University, and friends, including Jase Bernhardt, Drew Claybrook, Brian Colle, Daniel Horn, Daniel Michelson, Robert Pincus, Spencer Rhodes, Adam Sobel, David Stark, and Jeff Waldstreicher, who graciously agreed to plug in pressure sensors to their home internet. Figure 10 was made using code developed by Kevin Burris. The development and refinement of the methodology benefited from discussions with DelWayne Bohnenstiehl, Brian Colle, Brian Mapes, Matthew Parker, and Minghua Zhang.
This work was supported by the Center for Geospatial Analytics at North Carolina State University.
This research has been supported by the National Science Foundation Directorate for Geosciences (grant no. AGS-1905736), the National Aeronautics and Space Administration (grant no. 80NSSC19K0354), and the Office of Naval Research (grant no. N000142112116).
This paper was edited by Gerd Baumgarten and reviewed by four anonymous referees.
Adam, D.: Tonga Volcano created Puzzling Atmospheric Ripples, Nature, 602, 497, https://doi.org/10.1038/d41586-022-00127-1, 2022. a, b
Adams-Selin, R. D. and Johnson, R. H.: Mesoscale Surface Pressure and Temperature Features Associated with Bow Echoes, Mon. Weather Rev., 138, 212–227, https://doi.org/10.1175/2009MWR2892.1, 2010. a
Allen, G., Vaughan, G., Toniazzo, T., Coe, H., Connolly, P., Yuter, S. E., Burleyson, C. D., Minnis, P., and Ayers, J. K.: Gravity-wave-induced perturbations in marine stratocumulus: Gravity-Wave-Induced Perturbations in Marine Stratocumulus, Q. J. Roy. Meteor. Soc., 139, 32–45, https://doi.org/10.1002/qj.1952, 2013. a
Allen, L.: lrallen34/pressure-wave-detection-public: Code for Objective identification of pressure wave events from networks of 1 Hz, high-precision sensors, Zenodo [code], https://doi.org/10.5281/zenodo.10150876, 2023. a, b
Allen, L., Tomkins, L., and Yuter, S.: Supplemental videos of the paper “Objective identification of pressure wave events from networks of 1-Hz, high-precision sensors”, TIB AV-Portal [video], https://doi.org/10.5446/s_1476, 2023a. a
Allen, L., Tomkins, L., and Yuter, S.: 14 Sep 2021 KOKX Reflectivity and Doppler Velocity Waves. Supplemental videos of the paper “Objective identification of pressure wave events from networks of 1-Hz, high-precision sensors”, TIB AV-Portal [video], https://doi.org/10.5446/62542, 2023b. a
Allen, L., Tomkins, L., and Yuter, S.: 15 Nov 2020 KBUF Reflectivity and Doppler Velocity Waves. Supplemental videos of the paper “Objective identification of pressure wave events from networks of 1-Hz, high-precision sensors”, TIB AV-Portal [video], https://doi.org/10.5446/62541, 2023c. a
Allen, L., Tomkins, M., and Yuter, S.: 04 Feb 2022 KOKX Reflectivity and Doppler Velocity Waves. Supplemental videos of the paper “Objective identification of pressure wave events from networks of 1-Hz, high-precision sensors”, TIB AV-Portal [video], https://doi.org/10.5446/62540, 2023d. a
Allen, L., Tomkins, L., and Yuter, S.: 25 Feb 2022 KBUF Reflectivity and Doppler Velocity Waves. Supplemental videos of the paper “Objective identification of pressure wave events from networks of 1-Hz, high-precision sensors”, TIB AV-Portal [video], https://doi.org/10.5446/62539, 2023e. a
American Meteorological Society: Standard atmosphere – Glossary of Meteorology, https://glossary.ametsoc.org/wiki/Standard_atmosphere (last access: 18 December 2023), 2022. a
Amores, A., Monserrat, S., Marcos, M., Argüeso, D., Villalonga, J., Jordà, G., and Gomis, D.: Numerical Simulation of Atmospheric Lamb Waves Generated by the 2022 Hunga-Tonga Volcanic Eruption, Geophys. Res. Lett., 49, e2022GL098240, https://doi.org/10.1029/2022GL098240, 2022. a, b, c
Anthony, R. E., Ringler, A. T., Wilson, D. C., and Wolin, E.: Do Low-Cost Seismographs Perform Well Enough for Your Network? An Overview of Laboratory Tests and Field Observations of the OSOP Raspberry Shake 4D, Seismol. Res. Lett., 90, 219–228, https://doi.org/10.1785/0220180251, 2018. a, b
Bedard, A. J.: Infrasound Originating Near Mountainous Regions in Colorado, J. Appl. Meteorol. Clim., 17, 1014–1022, https://doi.org/10.1175/1520-0450(1978)017<1014:IONMRI>2.0.CO;2, 1978. a
Bhatia, A. and Fountain, H.: `It's Super Spectacular.' See How the Tonga Volcano Unleashed a Once-in-a-Century Shockwave, The New York Times, https://www.nytimes.com/interactive/2022/04/14/upshot/tonga-pressure-wave.html (last access: 18 December 2023), 2022. a
Bosart, L. F., Bracken, W. E., and Seimon, A.: A Study of Cyclone Mesoscale Structure with Emphasis on a Large-Amplitude Inertia–Gravity Wave, Mon. Weather Rev., 126, 1497–1527, https://doi.org/10.1175/1520-0493(1998)126<1497:ASOCMS>2.0.CO;2, 1998. a, b
Bosch: BMP388 Data Sheet, https://www.bosch-sensortec.com/products/environmental-sensors/pressure-sensors/bmp388/ (last access: 18 December 2023), 2020. a
Bosch: BME280 Data Sheet, https://www.bosch-sensortec.com/products/environmental-sensors/humidity-sensors-bme280/ (last access: 18 December 2023), 2022. a
Burt, S.: Multiple airwaves crossing Britain and Ireland following the eruption of Hunga Tonga–Hunga Ha'apai on 15 January 2022, Weather, 77, 76–81, https://doi.org/10.1002/wea.4182, 2022. a, b
Canavero, F. G. and Einaudi, F.: Time and Space Variability of Spectral Estimates of Atmospheric Pressure, J. Atmos. Sci., 44, 1589–1604, https://doi.org/10.1175/1520-0469(1987)044<1589:TASVOS>2.0.CO;2, 1987. a
Carr, J. L., Horváth, Á., Wu, D. L., and Friberg, M. D.: Stereo Plume Height and Motion Retrievals for the Record-Setting Hunga Tonga-Hunga Ha'apai Eruption of 15 January 2022, Geophys. Res. Lett., 49, e2022GL098131, https://doi.org/10.1029/2022GL098131, 2022. a
Christie, D. R., Muirhead, K. J., and Hales, A. L.: On Solitary Waves in the Atmosphere, J. Atmos. Sci., 35, 805–825, https://doi.org/10.1175/1520-0469(1978)035<0805:OSWITA>2.0.CO;2, 1978. a, b
Coffer, B. E. and Parker, M. D.: Infrasound signals in simulated nontornadic and pre-tornadic supercells, J. Acoust. Soc. Am., 151, 939–954, https://doi.org/10.1121/10.0009400, 2022. a
Connolly, P. J., Vaughan, G., Cook, P., Allen, G., Coe, H., Choularton, T. W., Dearden, C., and Hill, A.: Modelling the effects of gravity waves on stratocumulus clouds observed during VOCALS-UK, Atmos. Chem. Phys., 13, 7133–7152, https://doi.org/10.5194/acp-13-7133-2013, 2013. a, b
de Groot-Hedlin, C. D., Hedlin, M. A. H., and Walker, K. T.: Detection of gravity waves across the USArray: A case study, Earth Planet. Sc. Lett., 402, 346–352, https://doi.org/10.1016/j.epsl.2013.06.042, 2014. a, b
Del Pezzo, E. and Giudicepietro, F.: Plane wave fitting method for a plane, small aperture, short period seismic array: a MATHCAD program, Comput. Geosci., 28, 59–64, 2002. a, b
Eaton, S. W., Cárdenas, E. S., Hix, J. D., Johnson, J. T., Watson, S. M., Chichester, D. L., Garcés, M. A., Magaña-Zook, S. A., Maceira, M., Marcillo, O. E., Chai, C., d'Entremont, B. P., and Reichardt, T. A.: An algorithmic approach to predicting mechanical draft cooling tower fan speeds from infrasound signals, Appl. Acoust., 199, 109015, https://doi.org/10.1016/j.apacoust.2022.109015, 2022. a
Einaudi, F., Bedard, A. J., and Finnigan, J. J.: A Climatology of Gravity Waves and Other Coherent Disturbances at the Boulder Atmospheric Observatory during March–April 1984, J. Atmos. Sci., 46, 303–329, https://doi.org/10.1175/1520-0469(1989)046<0303:ACOGWA>2.0.CO;2, 1989. a, b
Fovell, R. G., Mullendore, G. L., and Kim, S.-H.: Discrete Propagation in Numerically Simulated Nocturnal Squall Lines, Mon. Weather Rev., 134, 3735–3752, https://doi.org/10.1175/MWR3268.1, 2006. a
Fritts, D. C. and Alexander, M. J.: Gravity wave dynamics and effects in the middle atmosphere, Rev. Geophys., 41, 1003, https://doi.org/10.1029/2001RG000106, 2003. a, b
Gaffin, D. M., Parker, S. S., and Kirkwood, P. D.: An Unexpectedly Heavy and Complex Snowfall Event across the Southern Appalachian Region, Weather Forecast., 18, 224–235, https://doi.org/10.1175/1520-0434(2003)018<0224:AUHACS>2.0.CO;2, 2003. a, b
Global Volcanism Program: Report on Hunga Tonga-Hunga Ha'apai (Tonga), in: Weekly Volcanic Activity Report, 12 January–18 January 2022, edited by: Sennert, S. K., Smithsonian Institution and US Geological Survey, https://volcano.si.edu/showreport.cfm?doi=GVP.WVAR20220112-243040 (last access: 18 December 2023), 2022. a, b
Gray, S. L., Martínez-Alvarado, O., Baker, L. H., and Clark, P. A.: Conditional symmetric instability in sting-jet storms, Q. J. Roy. Meteor. Soc., 137, 1482–1500, https://doi.org/10.1002/qj.859, 2011. a
Grivet-Talocia, S. and Einaudi, F.: Wavelet analysis of a microbarograph network, IEEE T. Geosci. Remote, 36, 418–432, https://doi.org/10.1109/36.662727, 1998. a, b, c, d, e, f, g, h, i
Grivet-Talocia, S., Einaudi, F., Clark, W. L., Dennett, R. D., Nastrom, G. D., and VanZandt, T. E.: A 4-yr Climatology of Pressure Disturbances Using a Barometer Network in Central Illinois, Mon. Weather Rev., 127, 1613–1629, https://doi.org/10.1175/1520-0493(1999)127<1613:AYCOPD>2.0.CO;2, 1999. a, b, c, d
Johnson, R. H. and Hamilton, P. J.: The Relationship of Surface Pressure Features to the Precipitation and Airflow Structure of an Intense Midlatitude Squall Line, Mon. Weather Rev., 116, 1444–1473, https://doi.org/10.1175/1520-0493(1988)116<1444:TROSPF>2.0.CO;2, 1988. a, b
Kjelaas, A. G., Beran, D. W., Hooke, W. H., and Bean, B. R.: Waves Observed in the Planetary Boundary Layer using an Array of Acoustic Sounders, J. Atmos. Sci., 31, 2040–2045, https://doi.org/10.1175/1520-0469(1974)031<2040:WOITPB>2.0.CO;2, 1974. a, b, c
Koch, S. E. and O'Handley, C.: Operational Forecasting and Detection of Mesoscale Gravity Waves, Weather Forecast., 12, 253–281, https://doi.org/10.1175/1520-0434(1997)012<0253:OFADOM>2.0.CO;2, 1997. a, b
Koch, S. E. and Saleeby, S.: An Automated System for the Analysis of Gravity Waves and Other Mesoscale Phenomena, Weather Forecast., 16, 661–679, https://doi.org/10.1175/1520-0434(2001)016<0661:AASFTA>2.0.CO;2, 2001. a, b
Lilly, J.: jLab: A data analysis package for Matlab, v.1.7.1, Zenodo, https://doi.org/10.5281/zenodo.4547006, 2021. a
Lilly, J. M. and Olhede, S. C.: Generalized Morse Wavelets as a Superfamily of Analytic Wavelets, IEEE T. Signal Proces., 60, 6036–6041, https://doi.org/10.1109/TSP.2012.2210890, 2012. a
Lindzen, R. S. and Tung, K.-K.: Banded Convective Activity and Ducted Gravity Waves, Mon. Weather Rev., 104, 1602–1617, https://doi.org/10.1175/1520-0493(1976)104<1602:BCAADG>2.0.CO;2, 1976. a, b
Markowski, P. and Richardson, Y.: Mesoscale Meteorology in Midlatitudes, John Wiley & Sons, Ltd, Chichester, UK, https://doi.org/10.1002/9780470682104, 2010. a, b, c
Matoza, R. S., Fee, D., Assink, J. D., Iezzi, A. M., Green, D. N., Kim, K., Toney, L., Lecocq, T., Krishnamoorthy, S., Lalande, J.-M., Nishida, K., Gee, K. L., Haney, M. M., Ortiz, H. D., Brissaud, Q., Martire, L., Rolland, L., Vergados, P., Nippress, A., Park, J., Shani-Kadmiel, S., Witsil, A., Arrowsmith, S., Caudron, C., Watada, S., Perttu, A. B., Taisne, B., Mialle, P., Le Pichon, A., Vergoz, J., Hupe, P., Blom, P. S., Waxler, R., De Angelis, S., Snively, J. B., Ringler, A. T., Anthony, R. E., Jolly, A. D., Kilgour, G., Averbuch, G., Ripepe, M., Ichihara, M., Arciniega-Ceballos, A., Astafyeva, E., Ceranna, L., Cevuard, S., Che, I.-Y., De Negri, R., Ebeling, C. W., Evers, L. G., Franco-Marin, L. E., Gabrielson, T. B., Hafner, K., Harrison, R. G., Komjathy, A., Lacanna, G., Lyons, J., Macpherson, K. A., Marchetti, E., McKee, K. F., Mellors, R. J., Mendo-Pérez, G., Mikesell, T. D., Munaibari, E., Oyola-Merced, M., Park, I., Pilger, C., Ramos, C., Ruiz, M. C., Sabatini, R., Schwaiger, H. F., Tailpied, D., Talmadge, C., Vidot, J., Webster, J., and Wilson, D. C.: Atmospheric waves and global seismoacoustic observations of the January 2022 Hunga eruption, Tonga, Science, 377, 95–100, https://doi.org/10.1126/science.abo7063, 2022. a
McMurdie, L. A., Heymsfield, G. M., Yorks, J. E., Braun, S. A., Skofronick-Jackson, G., Rauber, R. M., Yuter, S., Colle, B., McFarquhar, G. M., Poellot, M., Novak, D. R., Lang, T. J., Kroodsma, R., McLinden, M., Oue, M., Kollias, P., Kumjian, M. R., Greybush, S. J., Heymsfield, A. J., Finlon, J. A., McDonald, V. L., and Nicholls, S.: Chasing Snowstorms: The Investigation of Microphysics and Precipitation for Atlantic Coast-Threatening Snowstorms (IMPACTS) Campaign, B. Am. Meteorol. Soc., 103, E1243–E1269, https://doi.org/10.1175/BAMS-D-20-0246.1, 2022. a
Meyer, F.: Topographic distance and watershed lines, Signal Process., 38, 113–125, https://doi.org/10.1016/0165-1684(94)90060-4, 1994. a
Miller, M. and Allen, L.: Data for Objective identification of pressure wave events from networks of 1 Hz, high-precision sensors, Zenodo [data set], https://doi.org/10.5281/zenodo.8136536, 2023. a
Miller, M. A., Yuter, S. E., Hoban, N. P., Tomkins, L. M., and Colle, B. A.: Detecting wave features in Doppler radial velocity radar observations, Atmos. Meas. Tech., 15, 1689–1702, https://doi.org/10.5194/amt-15-1689-2022, 2022. a, b, c, d, e
Nappo, C. J.: An introduction to atmospheric gravity waves, Academic Press, USA, 276 pp., ISBN 0-12-514082-7, 2013. a, b, c
NOAA National Centers for Environmental Information: Automated Surface/Weather Observing Systems (ASOS/AWOS), NOAA National Centers for Environmental Information [data set], https://www.ncei.noaa.gov/products/land-based-station/automated-surface-weather-observing-systems (last access: 18 December 2023), 2021a. a, b
NOAA National Centers for Environmental Information: Global BUFR Data Stream: Upper Air Reports from the National Weather Service Telecommunications Gateway (NWS TG), NOAA National Centers for Environmental Information [data set], https://www.ncei.noaa.gov/access/metadata/landing-page/bin/iso?id=gov.noaa.ncdc:C01500 (last access: 18 December 2023), 2021b. a, b
NOAA National Weather Service Radar Operations Center: NOAA Next Generation Radar (NEXRAD) Level II Base Data, NOAA National Centers for Environmental Information [data set], https://doi.org/10.7289/V5W9574V, 1991. a, b
Olhede, S. and Walden, A.: Generalized Morse wavelets, IEEE T. Signal Proces., 50, 2661–2670, https://doi.org/10.1109/TSP.2002.804066, 2002. a
Pierce, A. D. and Posey, J. W.: Theory of the Excitation and Propagation of Lamb's Atmospheric Edge Mode from Nuclear Explosions, Geophys. J. Int., 26, 341–368, https://doi.org/10.1111/j.1365-246X.1971.tb03406.x, 1971. a
Proud, S. R., Prata, A. T., and Schmauß, S.: The January 2022 eruption of Hunga Tonga-Hunga Ha'apai volcano reached the mesosphere, Science, 378, 554–557, https://doi.org/10.1126/science.abo4076, 2022. a
ReVelle, D. O.: Acoustic-Gravity Waves from Bolide Sources, Earth Moon Planets, 102, 345–356, https://doi.org/10.1007/s11038-007-9181-3, 2008. a
Seifert, A. and Heus, T.: Large-eddy simulation of organized precipitating trade wind cumulus clouds, Atmos. Chem. Phys., 13, 5631–5645, https://doi.org/10.5194/acp-13-5631-2013, 2013. a
Tomkins, L. M., Yuter, S. E., Miller, M. A., and Allen, L. R.: Image muting of mixed precipitation to improve identification of regions of heavy snow in radar data, Atmos. Meas. Tech., 15, 5515–5525, https://doi.org/10.5194/amt-15-5515-2022, 2022. a, b, c, d, e, f, g, h
Torrence, C. and Compo, G. P.: A Practical Guide to Wavelet Analysis, B. Am. Meteorol. Soc., 79, 61–78, https://doi.org/10.1175/1520-0477(1998)079<0061:APGTWA>2.0.CO;2, 1998. a
Vergoz, J., Hupe, P., Listowski, C., Le Pichon, A., Garcés, M. A., Marchetti, E., Labazuy, P., Ceranna, L., Pilger, C., Gaebler, P., Näsholm, S. P., Brissaud, Q., Poli, P., Shapiro, N., De Negri, R., and Mialle, P.: IMS observations of infrasound and acoustic-gravity waves produced by the January 2022 volcanic eruption of Hunga, Tonga: A global analysis, Earth Planet. Sc. Lett., 591, 117639, https://doi.org/10.1016/j.epsl.2022.117639, 2022. a, b
Yuter, S. E., Hader, J. D., Miller, M. A., and Mechem, D. B.: Abrupt cloud clearing of marine stratocumulus in the subtropical southeast Atlantic, Science, 361, 697–701, https://doi.org/10.1126/science.aar5836, 2018. a
Zhang, F., Davis, C. A., Kaplan, M. L., and Koch, S. E.: Wavelet analysis and the governing dynamics of a large-amplitude mesoscale gravity-wave event along the East Coast of the United States, Q. J. Roy. Meteor. Soc., 127, 2209–2245, 2001. a