skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Atmospheric GNSS RO 1D-Var in Use at UCAR: Description and Validation
This paper describes, along with some validation results, the one-dimensional variational method (1D-Var) that is in use at the University Corporation for Atmospheric Research (UCAR) to retrieve atmospheric profiles of temperature, pressure, and humidity from the observation of the Global Navigation Satellite System (GNSS) radio occultation (RO). The retrieved profiles are physically consistent among the variables and statistically optimal as regards to a priori error statistics. Tests with idealized data demonstrate that the 1D-Var is highly effective in spreading the observational information and confirm that the method works as designed and expected, provided that correct input data are given. Tests for real-world data sets show that the retrieved profiles agree remarkably well with global weather analyses and collocated high vertical resolution radiosonde observations, and that the 1D-Var can produce value-added retrievals with respect to a priori profiles. We also find that the retrieved profiles are of exceptional long-term stability, suggesting that the 1D-Var can provide an excellent climate data record.  more » « less
Award ID(s):
2054356
PAR ID:
10433016
Author(s) / Creator(s):
; ; ; ;
Date Published:
Journal Name:
Remote Sensing
Volume:
14
Issue:
21
ISSN:
2072-4292
Page Range / eLocation ID:
5614
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Comparisons of atmospheric retrievals can reveal powerful insights on the strengths and limitations of our data and modeling tools. In this paper, we examine a sample of five L dwarfs of similar effective temperature (Teff) or spectral type to compare their pressure–temperature (P-T) profiles. Additionally, we explore the impact of an object’s metallicity and the signal-to-noise ratio (S/N) of the observations on the parameters we can retrieve. We present the first atmospheric retrievals: 2MASS J15261405+2043414, 2MASS J05395200−0059019, 2MASS J15394189−0520428, and GD 165B increasing the small but growing number of L dwarfs retrieved. When compared to the atmospheric retrievals of SDSS J141624.08+134826.7, a low-metallicity d/sdL7 primary in a wide L+T binary, we find that similarTeffsources have similar P-T profiles with metallicity differences impacting the relative offset between their P-T profiles in the photosphere. We also find that for near-infrared spectra, when the S/N is ≳80 we are in a regime where model uncertainties dominate over data measurement uncertainties. As such, S/N does not play a role in the retrieval’s ability to distinguish between a cloud-free and cloudless model, but may impact the confidence of the retrieved parameters. Lastly, we also discuss how to break cloud model degeneracies and the impact of extraneous gases in a retrieval model. 
    more » « less
  2. Abstract An integral equation method is presented for the 1D steady-state Poisson-Nernst-Planck equations modeling ion transport through membrane channels. The differential equations are recast as integral equations using Green’s 3rd identity yielding a fixed-point problem for the electric potential gradient and ion concentrations. The integrals are discretized by a combination of midpoint and trapezoid rules, and the resulting algebraic equations are solved by Gummel iteration. Numerical tests for electroneutral and non-electroneutral systems demonstrate the method’s 2nd order accuracy and ability to resolve sharp boundary layers. The method is applied to a 1D model of the K$$^+$$ +  ion channel with a fixed charge density that ensures cation selectivity. In these tests, the proposed integral equation method yields potential and concentration profiles in good agreement with published results. 
    more » « less
  3. Abstract. The West Antarctic Peninsula (WAP) is a rapidly warming region, withsubstantial ecological and biogeochemical responses to the observed changeand variability for the past decades, revealed by multi-decadal observationsfrom the Palmer Antarctica Long-Term Ecological Research (LTER) program. Thewealth of these long-term observations provides an important resource forecosystem modeling, but there has been a lack of focus on the developmentof numerical models that simulate time-evolving plankton dynamics over theaustral growth season along the coastal WAP. Here, we introduce aone-dimensional variational data assimilation planktonic ecosystem model (i.e., theWAP-1D-VAR v1.0 model) equipped with a modelparameter optimization scheme. We first demonstrate the modified and newlyadded model schemes to the pre-existing food web and biogeochemicalcomponents of the other ecosystem models that WAP-1D-VAR model was adaptedfrom, including diagnostic sea-ice forcing and trophic interactions specificto the WAP region. We then present the results from model experiments wherewe assimilate 11 different data types from an example Palmer LTER growthseason (October 2002–March 2003) directly related to corresponding modelstate variables and flows between these variables. The iterative dataassimilation procedure reduces the misfits between observationsand model results by 58 %, compared to before optimization, via an optimized set of12 parameters out of a total of 72 free parameters. The optimized model resultscapture key WAP ecological features, such as blooms during seasonal sea-iceretreat, the lack of macronutrient limitation, and modeled variables andflows comparable to other studies in the WAP region, as well as severalimportant ecosystem metrics. One exception is that the model slightlyunderestimates particle export flux, for which we discuss potentialunderlying reasons. The data assimilation scheme of the WAP-1D-VAR modelenables the available observational data to constrain previously poorlyunderstood processes, including the partitioning of primary production bydifferent phytoplankton groups, the optimal chlorophyll-to-carbon ratio ofthe WAP phytoplankton community, and the partitioning of dissolved organiccarbon pools with different lability. The WAP-1D-VAR model can besuccessfully employed to link the snapshots collected by the available datasets together to explain and understand the observed dynamics along thecoastal WAP. 
    more » « less
  4. null (Ed.)
    Surface albedo is a fundamental radiative parameter as it controls the Earth’s energy budget and directly affects the Earth’s climate. Satellite observations have long been used to capture the temporal and spatial variations of surface albedo because of their continuous global coverage. However, space-based albedo products are often affected by errors in the atmospheric correction, multi-angular bi-directional reflectance distribution function (BRDF) modelling, as well as spectral conversions. To validate space-based albedo products, an in situ tower albedometer is often used to provide continuous “ground truth” measurements of surface albedo over an extended area. Since space-based albedo and tower-measured albedo are produced at different spatial scales, they can be directly compared only for specific homogeneous land surfaces. However, most land surfaces are inherently heterogeneous with surface properties that vary over a wide range of spatial scales. In this work, tower-measured albedo products, including both directional hemispherical reflectance (DHR) and bi-hemispherical reflectance (BHR), are upscaled to coarse satellite spatial resolutions using a new method. This strategy uses high-resolution satellite derived surface albedos to fill the gaps between the albedometer’s field-of-view (FoV) and coarse satellite scales. The high-resolution surface albedo is generated from a combination of surface reflectance retrieved from high-resolution Earth Observation (HR-EO) data and moderate resolution imaging spectroradiometer (MODIS) BRDF climatology over a larger area. We implemented a recently developed atmospheric correction method, the Sensor Invariant Atmospheric Correction (SIAC), to retrieve surface reflectance from HR-EO (e.g., Sentinel-2 and Landsat-8) top-of-atmosphere (TOA) reflectance measurements. This SIAC processing provides an estimated uncertainty for the retrieved surface spectral reflectance at the HR-EO pixel level and shows excellent agreement with the standard Landsat 8 Surface Reflectance Code (LaSRC) in retrieving Landsat-8 surface reflectance. Atmospheric correction of Sentinel-2 data is vastly improved by SIAC when compared against the use of in situ AErosol RObotic NETwork (AERONET) data. Based on this, we can trace the uncertainty of tower-measured albedo during its propagation through high-resolution EO measurements up to coarse satellite scales. These upscaled albedo products can then be compared with space-based albedo products over heterogeneous land surfaces. In this study, both tower-measured albedo and upscaled albedo products are examined at Ground Based Observation for Validation (GbOV) stations (https://land.copernicus.eu/global/gbov/), and used to compare with satellite observations, including Copernicus Global Land Service (CGLS) based on ProbaV and VEGETATION 2 data, MODIS and multi-angle imaging spectroradiometer (MISR). 
    more » « less
  5. This work extends the applicability of our recent convexification- based algorithm for constructing images of the dielectric constant of buried or occluded target. We are orientated towards the detection of explosive-like targets such as antipersonnel land mines and improvised explosive devices in the non-invasive inspections of buildings. In our previous work, the method is posed in the perspective that we use multiple source locations running along a line of source to get a 2D image of the dielectric function. Mathematically, we solve a 1D coefficient inverse problem for a hyperbolic equation for each source location. Different from any conventional Born approximation-based technique for synthetic-aperture radar, this method does not need any linearization. In this paper, we attempt to verify the method using several 3D numerical tests with simulated data. We revisit the global convergence of the gradient descent method of our computational approach. 
    more » « less