skip to main content

Title: Data Assimilation of High‐Latitude Electric Fields: Extension of a Multi‐Resolution Gaussian Process Model (Lattice Kriging) to Vector Fields

We develop a new methodology for the multi‐resolution assimilation of electric fields by extending a Gaussian process model (Lattice Kriging) used for scalar field originally to vector field. This method takes the background empirical model as “a priori” knowledge and fuses real observations under the Gaussian process framework. The comparison of assimilated results under two different background models and three different resolutions suggests that (a) the new method significantly reduces fitting errors compared with the global spherical harmonic fitting (SHF) because it uses range‐limited basis functions ideal for the local fitting and (b) the fitting resolution, determined by the number of basis functions, is adjustable and higher resolution leads to smaller errors, indicating that more structures in the data are captured. We also test the sensitivity of the fitting results to the total amount of input data: (a) as the data amount increases, the fitting results deviate from the background model and become more determined by data and (b) the impacts of data can reach remote regions with no data available. The assimilation also better captures short‐period variations in local PFISR measurements than the SHF and maintains a coherent pattern with the surrounding. The multi‐resolution Lattice Kriging is examined via attributing basis functions into multiple levels with different resolutions (fine level is located in the region with observations). Such multi‐resolution fitting has the smallest error and shortest computation time, making the regional high‐resolution modeling efficient. Our method can be modified to achieve the multi‐resolution assimilation for other vector fields from unevenly distributed observations.

more » « less
Award ID(s):
1753214 2149695
Author(s) / Creator(s):
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Space Weather
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Mountain winds are the driving force behind snow accumulation patterns in mountainous catchments, making accurate wind fields a prerequisite to accurate simulations of snow depth for ecological or water resource applications. In this study, we examine the effect that wind fields derived from different coarse data sets and downscaling schemes has on simulations of modeled snow depth at resolutions suitable for basin‐scale modeling (>50 m). Simulations are run over the Tuolumne River Basin, CA for the accumulation season of Water Year 2017 using the distributed snow model, SnowModel. We derived wind fields using observations from either sensor networks, the North American Land Data Assimilation System (12.5 km resolution), or High‐Resolution Rapid Refresh (HRRR, 3 km resolution) data, and downscaled using terrain‐based multipliers (MicroMet), a mass‐conserving flow model (WindNinja), or bilinear interpolation. Two wind fields derived from 3 km HRRR data and downscaled with respect to terrain produced snow depth maps that best matched observations of snow depth from airborne LiDAR. We find that modeling these wind fields at 100 and 50 m resolutions do not produce improvements in simulated snow depth when compared to wind fields modeled at a 150 m resolution due to their inability to represent wind dynamics at these scales. For input to distributed snow models at the basin‐scale, we recommend deriving wind fields from high resolution numerical weather prediction model output and downscaling with respect to terrain. Future studies should compare suspension schemes used in blowing snow models and investigate wind downscaling schemes of complexity between statistical and fluid dynamic models.

    more » « less
  2. Context.3C 84 is a nearby radio source with a complex total intensity structure, showing linear polarisation and spectral patterns. A detailed investigation of the central engine region necessitates the use of very-long-baseline interferometry (VLBI) above the hitherto available maximum frequency of 86 GHz.

    Aims.Using ultrahigh resolution VLBI observations at the currently highest available frequency of 228 GHz, we aim to perform a direct detection of compact structures and understand the physical conditions in the compact region of 3C 84.

    Methods.We used Event Horizon Telescope (EHT) 228 GHz observations and, given the limited (u, v)-coverage, applied geometric model fitting to the data. Furthermore, we employed quasi-simultaneously observed, ancillary multi-frequency VLBI data for the source in order to carry out a comprehensive analysis of the core structure.

    Results.We report the detection of a highly ordered, strong magnetic field around the central, supermassive black hole of 3C 84. The brightness temperature analysis suggests that the system is in equipartition. We also determined a turnover frequency ofνm = (113 ± 4) GHz, a corresponding synchrotron self-absorbed magnetic field ofBSSA = (2.9 ± 1.6) G, and an equipartition magnetic field ofBeq = (5.2 ± 0.6) G. Three components are resolved with the highest fractional polarisation detected for this object (mnet = (17.0 ± 3.9)%). The positions of the components are compatible with those seen in low-frequency VLBI observations since 2017–2018. We report a steeply negative slope of the spectrum at 228 GHz. We used these findings to test existing models of jet formation, propagation, and Faraday rotation in 3C 84.

    Conclusions.The findings of our investigation into different flow geometries and black hole spins support an advection-dominated accretion flow in a magnetically arrested state around a rapidly rotating supermassive black hole as a model of the jet-launching system in the core of 3C 84. However, systematic uncertainties due to the limited (u, v)-coverage, however, cannot be ignored. Our upcoming work using new EHT data, which offer full imaging capabilities, will shed more light on the compact region of 3C 84.

    more » « less
  3. Abstract

    Remote sensing data have been widely used to study various geophysical processes. With the advances in remote sensing technology, massive amount of remote sensing data are collected in space over time. Different satellite instruments typically have different footprints, measurement‐error characteristics, and data coverages. To combine data sets from different satellite instruments, we propose a dynamic fused Gaussian process (DFGP) model that enables fast statistical inference such as filtering and smoothing for massive spatio‐temporal data sets in a data‐fusion context. Based upon a spatio‐temporal‐random‐effect model, the DFGP methodology represents the underlying true process with two components: a linear combination of a small number of basis functions and random coefficients with a general covariance matrix, together with a linear combination of a large number of basis functions and Markov random coefficients. To model the underlying geophysical process at different spatial resolutions, we rely on the change‐of‐support property, which also allows efficient computations in the DFGP model. To estimate model parameters, we devise a computationally efficient stochastic expectation‐maximization algorithm to ensure its scalability for massive data sets. The DFGP model is applied to a total of 3.7 million sea surface temperature data sets in the tropical Pacific Ocean for a one‐week time period in 2010 from Moderate Resolution Imaging Spectroradiometer (MODIS) and Advanced Microwave Scanning Radiometer‐Earth Observing System (AMSR‐E) instruments.

    more » « less
  4. Abstract

    Gridded monthly rainfall estimates can be used for a number of research applications, including hydrologic modeling and weather forecasting. Automated interpolation algorithms, such as the “autoKrige” function in R, can produce gridded rainfall estimates that validate well but produce unrealistic spatial patterns. In this work, an optimized geostatistical kriging approach is used to interpolate relative rainfall anomalies, which are then combined with long-term means to develop the gridded estimates. The optimization consists of the following: 1) determining the most appropriate offset (constant) to use when log-transforming data; 2) eliminating poor quality data prior to interpolation; 3) detecting erroneous maps using a machine learning algorithm; and 4) selecting the most appropriate parameterization scheme for fitting the model used in the interpolation. Results of this effort include a 30-yr (1990–2019), high-resolution (250-m) gridded monthly rainfall time series for the state of Hawai‘i. Leave-one-out cross validation (LOOCV) is performed using an extensive network of 622 observation stations. LOOCV results are in good agreement with observations (R2= 0.78; MAE = 55 mm month−1; 1.4%); however, predictions can underestimate high rainfall observations (bias = 34 mm month−1; −1%) due to a well-known smoothing effect that occurs with kriging. This research highlights the fact that validation statistics should not be the sole source of error assessment and that default parameterizations for automated interpolation may need to be modified to produce realistic gridded rainfall surfaces. Data products can be accessed through the Hawai‘i Data Climate Portal (HCDP;

    Significance Statement

    A new method is developed to map rainfall in Hawai‘i using an optimized geostatistical kriging approach. A machine learning technique is used to detect erroneous rainfall maps and several conditions are implemented to select the optimal parameterization scheme for fitting the model used in the kriging interpolation. A key finding is that optimization of the interpolation approach is necessary because maps may validate well but have unrealistic spatial patterns. This approach demonstrates how, with a moderate amount of data, a low-level machine learning algorithm can be trained to evaluate and classify an unrealistic map output.

    more » « less
  5. ArcticDEM provides the public with an unprecedented opportunity to access very high-spatial resolution digital elevation models (DEMs) covering the pan-Arctic surfaces. As it is generated from stereo-pairs of optical satellite imagery, ArcticDEM represents a mixture of a digital surface model (DSM) over a non-ground areas and digital terrain model (DTM) at bare grounds. Reconstructing DTM from ArcticDEM is thus needed in studies requiring bare ground elevation, such as modeling hydrological processes, tracking surface change dynamics, and estimating vegetation canopy height and associated forest attributes. Here we proposed an automated approach for estimating DTM from ArcticDEM in two steps: (1) identifying ground pixels from WorldView-2 imagery using a Gaussian mixture model (GMM) with local refinement by morphological operation, and (2) generating a continuous DTM surface using ArcticDEMs at ground locations and spatial interpolation methods (ordinary kriging (OK) and natural neighbor (NN)). We evaluated our method at three forested study sites characterized by different canopy cover and topographic conditions in Livengood, Alaska, where airborne lidar data is available for validation. Our results demonstrate that (1) the proposed ground identification method can effectively identify ground pixels with much lower root mean square errors (RMSEs) (<0.35 m) to the reference data than the comparative state-of-the-art approaches; (2) NN performs more robustly in DTM interpolation than OK; (3) the DTMs generated from NN interpolation with GMM-based ground masks decrease the RMSEs of ArcticDEM to 0.648 m, 1.677 m, and 0.521 m for Site-1, Site-2, and Site-3, respectively. This study provides a viable means of deriving high-resolution DTM from ArcticDEM that will be of great value to studies focusing on the Arctic ecosystems, forest change dynamics, and earth surface processes. 
    more » « less