skip to main content


This content will become publicly available on September 1, 2024

Title: Log-Gaussian Cox process modeling of large spatial lightning data using spectral and Laplace approximations
Lightning is a destructive and highly visible product of severe storms, yet there is still much to be learned about the conditions under which lightning is most likely to occur. The GOES-16 and GOES-17 satellites, launched in 2016 and 2018 by NOAA and NASA, collect a wealth of data regarding individual lightning strike occurrence and potentially related atmospheric variables. The acute nature and inherent spatial correlation in lightning data renders standard regression analyses inappropriate. Further, computational considerations are foregrounded by the desire to analyze the immense and rapidly increasing volume of lightning data. We present a new computationally feasible method that combines spectral and Laplace approximations in an EM algorithm, denoted SLEM, to fit the widely popular log-Gaussian Cox process model to large spatial point pattern datasets. In simulations we find SLEM is competitive with contemporary techniques in terms of speed and accuracy. When applied to two lightning datasets, SLEM provides better out-of-sample prediction scores and quicker runtimes, suggesting its particular usefulness for analyzing lightning data which tend to have sparse signals.  more » « less
Award ID(s):
1916208
NSF-PAR ID:
10485004
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
Project Euclid
Date Published:
Journal Name:
The Annals of Applied Statistics
Volume:
17
Issue:
3
ISSN:
1932-6157
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    We relate the spatial and temporal distribution of lightning flash rates and cloud top brightness temperature (CTBT) to concentric atmospheric gravity wave (CGW) events observed at the Southern Space Observatory (SSO) in São Martinho da Serra (29.44°S, 53.82°W, 488.7 m) in southern Brazil. The selected identified cases from 2017 to 2018 were observed by a hydroxyl (OH) all‐sky imager. Backward ray tracing shows that the time of gravity wave excitation agrees with the highest values of lightning flash rates (indicating lightning jump) as well as the coldest brightness temperatures that indicate the time of convective overshoot. Radiosonde measurements show high convective available potential energy (CAPE), associated with a maximum updraft velocity just prior to the wave events. We find that these possible source locations correspond to the positions and times that convective plumes overshot the tropopause (seen in GOES‐16 CTBT images). We also show that higher spatial lightning density (i.e., number of lightning flashes at a given longitude and latitude) agree with the overshoot locations from the GOES satellite. We also find that the overshoot times from the GOES‐16 satellite agree with the times lightning jumps were observed in the lightning flash rate. Finally, we find that the periodicities in the lightning flash rate agree with the periods of the observed CGWs, which further strengthens the result that the CGWs were excited by the deep convective systems determined from backward ray tracing.

     
    more » « less
  2. Abstract

    Global satellite studies show a maximum in deep convection and lightning downstream of the Andes in subtropical South America. The Remote sensing of Electrification, Lightning, And Mesoscale/microscale Processes with Adaptive Ground Observations (RELAMPAGO) field campaign was designed to investigate the physical processes that contribute to the rapid development of deep convection and mesoscale convective systems (MCSs) in Argentina. A lightning mapping array (LMA) was deployed to Argentina as part of RELAMPAGO to collect lightning observations from extreme storms in the region. This study combines lightning data from the LMA and the Geostationary Lightning Mapper onboardGOES‐16with 1‐km gridded radar data to examine the electrical characteristics of a variety of convective storms throughout their life cycle observed during RELAMPAGO. Results from the full campaign show 48% of flashes are associated with deep convection that occurs along the eastern edge of the Sierras de Córdoba (SDC) overnight. These flashes are 65 km2smaller on average compared to stratiform flashes, which occur most frequently 50–100 km east of the SDC in the early morning hours, consistent with the upscale growth of MCSs off the terrain. Analysis of the 13–14 December MCS shows that sharp increases in flash rates correspond to deep and wide convective cores that have high graupel and hail mass, 35‐dBZ volume, and ice water path. This work validates previous satellite studies of lightning in the region, but also provides higher spatial and temporal resolution information across the convective life cycle that has not been available in previous studies.

     
    more » « less
  3. In this study, we present a nationwide machine learning model for hourly PM2.5 estimation for the continental United States (US) using high temporal resolution Geostationary Operational Environmental Satellites (GOES-16) Aerosol Optical Depth (AOD) data, meteorological variables from the European Center for Medium Range Weather Forecasting (ECMWF) and ancillary data collected between May 2017 and December 2020. A model sensitivity analysis was conducted on predictor variables to determine the optimal model. It turns out that GOES16 AOD, variables from ECMWF, and ancillary data are effective variables in PM2.5 estimation and historical reconstruction, which achieves an average mean absolute error (MAE) of 3.0 μg/m3, and a root mean square error (RMSE) of 5.8 μg/m3. This study also found that the model performance as well as the site measured PM2.5 concentrations demonstrate strong spatial and temporal patterns. Specifically, in the temporal scale, the model performed best between 8:00 p.m. and 11:00 p.m. (UTC TIME) and had the highest coefficient of determination (R2) in Autumn and the lowest MAE and RMSE in Spring. In the spatial scale, the analysis results based on ancillary data show that the R2 scores correlate positively with the mean measured PM2.5 concentration at monitoring sites. Mean measured PM2.5 concentrations are positively correlated with population density and negatively correlated with elevation. Water, forests, and wetlands are associated with low PM2.5 concentrations, whereas developed, cultivated crops, shrubs, and grass are associated with high PM2.5 concentrations. In addition, the reconstructed PM2.5 surfaces serve as an important data source for pollution event tracking and PM2.5 analysis. For this purpose, from May 2017 to December 2020, hourly PM2.5 estimates were made for 10 km by 10 km and the PM2.5 estimates from August through November 2020 during the period of California Santa Clara Unite (SCU) Lightning Complex fires are presented. Based on the quantitative and visualization results, this study reveals that a number of large wildfires in California had a profound impact on the value and spatial-temporal distributions of PM2.5 concentrations. 
    more » « less
  4. Abstract

    Overshooting tops (OTs) are a well‐known indicator of updrafts capable of transporting air from the troposphere to the stratosphere and generating hazardous weather conditions. Satellites and radars have long been used to identify OTs, but the results have not been entirely consistent due to differences in sensor and measurement characteristics. OT detection approaches based on satellite infrared (IR) imagery have often been validated using human‐expert OT identifications, but such datasets are time‐consuming to compile over broad geographic regions. Despite radar limitations to detect the true physical cloud top, OTs identified within multi‐radar composites can serve as a stable reference for comprehensive satellite OT analysis and detection validation. This study analyzes a large OT data set compiled from Geostationary Operational Environmental Satellites (GOES)‐13/16 geostationary IR data and gridded volumetric Next‐Generation Radar (NEXRAD) reflectivity to better understand radar and IR observations of OTs, quantify agreement between satellite and radar OT detections, and demonstrate how an increased spatial sampling from GOES‐13 to GOES‐16 impacts OT appearance and detection performance. For nearly time‐matched scenes and moderate OT probability, the GOES‐13 detection rate (∼60%) is ∼15% lower than GOES‐16 (∼75%), which is mostly attributed to coarser spatial resolution. NEXRAD column‐maximum reflectivity and tropopause‐relative echo‐top height as a function of GOES OT probability were quite consistent between the two satellites however, indicating that efforts to account for differing resolution were largely successful. GOES false detections are unavoidable because outflow from nearby or recently decayed OTs can be substantially colder than the tropopause and look like an OT to an automated algorithm.

     
    more » « less
  5. Abstract

    The decision to establish a network of researchers centers on identifying shared research goals. Ecologically specific regions, such as the USA’s National Ecological Observatory Network’s (NEON’s) eco-climatic domains, are ideal locations by which to assemble researchers with a diverse range of expertise but focused on the same set of ecological challenges. The recently established Great Lakes User Group (GLUG) is NEON’s first domain specific ensemble of researchers, whose goal is to address scientific and technical issues specific to the Great Lakes Domain 5 (D05) by using NEON data to enable advancement of ecosystem science. Here, we report on GLUG’s kick off workshop, which comprised lightning talks, keynote presentations, breakout brainstorming sessions and field site visits. Together, these activities created an environment to foster and strengthen GLUG and NEON user engagement. The tangible outcomes of the workshop exceeded initial expectations and include plans for (i) two journal articles (in addition to this one), (ii) two potential funding proposals, (iii) an assignable assets request and (iv) development of classroom activities using NEON datasets. The success of this 2.5-day event was due to a combination of factors, including establishment of clear objectives, adopting engaging activities and providing opportunities for active participation and inclusive collaboration with diverse participants. Given the success of this approach we encourage others, wanting to organize similar groups of researchers, to adopt the workshop framework presented here which will strengthen existing collaborations and foster new ones, together with raising greater awareness and promotion of use of NEON datasets. Establishing domain specific user groups will help bridge the scale gap between site level data collection and addressing regional and larger ecological challenges.

     
    more » « less