Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reanalysis data”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 793 records · Page 44

Drivers of Cloud Condensation Nuclei in the Eastern North Atlantic as Observed at the ARM Site

Abstract Seven years of data collected at the Atmospheric Radiation Measurement's Eastern North Atlantic (ENA) site are analyzed to understand the controls of Cloud Condensation Nuclei (CCN) concentrations in the region. Day‐night differences in the aerosol data as segregated by wind direction demonstrate the aerosol observations to be impacted by local emissions when the wind direction (wdir) is between 90° and 310° (measured clockwise from the North where air is coming from). Data collected during marine conditions (wdir <90° or wdir >310°) show the CCN concentrations to be higher in the summer months as compared to the winter months. CCN budget analysis revealed advection and precipitation scavenging being primarily responsible for modulating the CCN concentrations at the site on monthly timescales, with rain rates driving the precipitation scavenging term. High (greater than 75th percentile) and low (lower than 25th percentile) CCN events were identified for each month to characterize the sub‐monthly variability of CCN concentrations. Low CCN events had thicker clouds, stronger rain rates, and lower reanalysis reported free‐tropospheric aerosol pseudo number concentration at the ENA site as compared to the high CCN events. Analysis of satellite data of air‐parcels 48 hr prior to their arrival at the ENA site demonstrated the air parcels during low CCN events to encounter higher cloudiness, stronger rain rates, and higher cloud top heights as compared to the high CCN events. The results presented herein provide key constraints for model evaluation studies and climatological studies conducted at the ENA site.

54 ENVIRONMENTAL SCIENCES↗

The role of plasma wave turbulence in the formation of shock waves in collisionless plasmas

A data set containing measurements for approximately 100 shock crossings was compiled, and plasma wave profiles were plotted in conjunction with magnetic field data. Three emissions analyzed include upstream electron plasma oscillations, ion acoustic noise, and low frequency electromagnetic noise. The data demonstrate that the observed waves are consistent with generation by the kinetic cross-field streaming instability. Reanalysis of shock overshoots prompted by the presence of plasma feet upstream of nominally subcritical shocks leads to the conclusion that no firm evidence exists for a sharp subcritical supercritical transition. Electron plasma oscillations persist up to the ramps of subcritical shocks and begin decaying at the front edge of the feet of supercritical shocks. Their intensity is anticorrelated with the amplitude of the foot ion acoustic noise. Wave measurements support the predicted predominance of lower hybrid-like modes over ion acoustic noise in the production of resistive (plasma) heating. Evidence for variable shock normal angles is seen upstream from nominally quasiparallel shocks. The first cases of very strong electron heating observed at the terrestrial bow shock were discovered.

Mellott, M. M.↗

Remote Sensing of Dipole Rings

Historical satellite-derived sea surface temperature (SST) data are reanalyzed with a zebra color palette and a thermal separatrix method. The new results from this reanalysis are as follows: (a) Thirteen observational sequences of six rings from the Gulf Stream and the Brazil Current, which have historically been interpreted as solitary vortices or monopoles are shown to have a dipolar character; (b) some of these dipole rings have been observed in the open ocean, thereby eliminating the possibility that they are sustained by topographic interactions with the continental slope; (c) whether interacting with other features or evolving as isolated circulations, dipoles are seen to rotate within a relatively narrow range of approximately 4-8 deg/day (interacting) and 10-11 deg/day (isolated); and (d) feature tracking delineates energetic fluid in both vortices and eliminates the possibility of interpreting dipole rings as transient features produced by active monopoles and patches of entrained fluid.

Hooker, Stanford B.↗

Impact of Satellite Viewing-Swath Width on Global and Regional Aerosol Optical Thickness Statistics and Trends

We use the Moderate Resolution Imaging Spectroradiometer (MODIS) satellite aerosol optical thickness (AOT) product to assess the impact of reduced swath width on global and regional AOT statistics and trends. Alongtrack and across-track sampling strategies are employed, in which the full MODIS data set is sub-sampled with various narrow-swath (approximately 400-800 km) and single pixel width (approximately 10 km) configurations. Although view-angle artifacts in the MODIS AOT retrieval confound direct comparisons between averages derived from different sub-samples, careful analysis shows that with many portions of the Earth essentially unobserved, spatial sampling introduces uncertainty in the derived seasonal-regional mean AOT. These AOT spatial sampling artifacts comprise up to 60%of the full-swath AOT value under moderate aerosol loading, and can be as large as 0.1 in some regions under high aerosol loading. Compared to full-swath observations, narrower swath and single pixel width sampling exhibits a reduced ability to detect AOT trends with statistical significance. On the other hand, estimates of the global, annual mean AOT do not vary significantly from the full-swath values as spatial sampling is reduced. Aggregation of the MODIS data at coarse grid scales (10 deg) shows consistency in the aerosol trends across sampling strategies, with increased statistical confidence, but quantitative errors in the derived trends are found even for the full-swath data when compared to high spatial resolution (0.5 deg) aggregations. Using results of a model-derived aerosol reanalysis, we find consistency in our conclusions about a seasonal-regional spatial sampling artifact in AOT Furthermore, the model shows that reduced spatial sampling can amount to uncertainty in computed shortwave top-ofatmosphere aerosol radiative forcing of 2-3 W m(sup−2). These artifacts are lower bounds, as possibly other unconsidered sampling strategies would perform less well. These results suggest that future aerosol satellite missions having significantly less than full-swath viewing are unlikely to sample the true AOT distribution well enough to obtain the statistics needed to reduce uncertainty in aerosol direct forcing of climate.

remote sensing↗

Comparing Flood-Producing Catastrophic Hurricanes in GEOS: A Study of Juan 1985 and Harvey 2017

Hurricanes Juan (1985) and Harvey (2017) bear remarkable similarities because of the exceptional amount of accumulated precipitation over coastal regions in the Gulf, consequent to the peculiar and rare situation in which a hurricane stalls for several days on a coastal area, displaying minimal progress. This is particularly dangerous because: a) terrestrial drag affect only part of the circulation so that weakening is not as fast as when there is rapid motion inland; b) a part of the circulation remains over water and as such low-level moisture flux convergence is maintained; c) precipitation lasts for a very long time over the same area; d) moisture recycling from flooded areas further contributes to maintaining flood conditions. From the forecast perspective, hurricanes with track singularities such as stalling points, loops, or multiple landfalls over the same region are exceptionally difficult to predict because they result from small variations of opposing mechanisms and are therefore intrinsically unstable. Intensity fluctuations and track singularities are strongly coupled. The representation of fluxes undergoing differential drag (ocean vs. land) inside the same circulation is critically dependent upon resolution. Recent modeling development efforts within the NASA Global Modeling and Assimilation Office (GMAO) have brought the horizontal resolution of the global cubed-sphere Goddard Earth Observing System (GEOS) atmospheric model down to 12 km. This has led to innovative findings in the field of tropical cyclone (TC) modeling within a global framework. The GEOS atmospheric model and data assimilation system (DAS), used to produce the MERRA-2 (Modern Era Retrospective analysis for Research and Applications) reanalysis, has a "Replay" capability that adds a forcing term to the model equations to constrain it to follow a specified trajectory. This capability is used to produce global, downscaled versions of the 50 km MERRA-2 reanalysis at a 12.5 km resolution. The full DAS and replay modes both use the Incremental Analysis Update or IAU for applying an analysis increment to the model's background state four times each day at 00Z, 06Z, 12Z and 18Z. In a full DAS mode, a six-hour free run of the model is made (Predictor Step), writing states that will be used as backgrounds or first-guesses for the analysis centered on 00Z. The Analysis system then produces its best estimate of the 00Z state by blending backgrounds with observations. The difference between this best estimate and the model's original first guess at 00Z is treated as an analysis increment and added uniformly as a forcing term to the model during a second, now forced, run of the 12Z to 03Z window, which produces the assimilated results. Replay mode is very similar, except only the 00Z background is needed from the predictor segment, which is thus half as long, and instead of doing an analysis, a previously analyzed state at 00Z is read in. The increment is then computed and applied to the corrector step in the same way as in the original DAS. When downscaling horizontal spatial filters can be applied to the "analysis" increment used to force the Corrector step, so as to reproduce the full meteorological state or selectively choose or discard certain forcings. In this presentation, full global hindcast simulations of Juan and operational forecasts of Harvey are compared.

Flooding↗

A Modeling Study of Oceanic Response to Daily and Monthly Surface Forcing

The goal of this study is to investigate the effect of high-frequency surface forcing (wind stresses and heat fluxes) on upper-ocean response. We use the reduced-gravity quasi-isopycnal ocean model by Schopf and Loughe (1995) for this study. Two experiments are performed: one with daily and the other with monthly surface forcing. The two experiments are referred to as DD and MM, respectively. The daily surface wind stress is produced from the SSM/I wind data (Atlas et al. 1991) using the drag coefficient of Large and Pond (1982). The surface latent and sensible heat fluxes are estimated using the atmospheric mixed layer model by Seager et al. (1995) with the time-varying air temperature and specific humidity from the NCEP-NCAR reanalysis (Kalnay et al. 1996). The radiation is based on climatological shortwave radiation from the Earth Radiation Budget Experiment (ERBE) [Harrison et al. 1993] and the daily GEWEX SRB data. The ocean model domain is restricted to the Pacific Ocean with realistic land boundaries. At the southern boundary the model temperature and salinity are relaxed to the Levitus (1994) climatology. The time-mean SST distribution from MM is close to the observed SST climatology while the mean SST field from DD is about 1.5 C cooler. To identify the responsible processes, we examined the mean heat budgets and the heat balance during the first year (when the difference developed) in the two experiments. The analysis reveals that this is contributed by two factors. One is the difference in latent heat flux. The other is the difference in mixing processes. To further evaluate the responsible processes, we repeated the DD experiment by reducing the based vertical diffusion from 1e-4 to 0.5e-5. The resultant SST field becomes quite closer to the observed SST field. SST variability from the two experiments is generally similar, but the equatorial SST differences between the two experiments show interannual variations. We are investigating the possible mechanisms responsible for the different responses.

Sui, Chung-Hsiung↗

Offshore reanalysis wind speed assessment across the wind turbine rotor layer off the United States Pacific coast

Abstract. The California Pacific coast is characterized by considerable wind resource and areas of dense population, propelling interest in offshore wind energy as the United States moves toward a sustainable and decarbonized energy future. Reanalysis models continue to serve the wind energy community in a multitude of ways, and the need for validation in locations where observations have been historically limited, such as offshore environments, is strong. The U.S. Department of Energy (DOE) owns two lidar buoys that collect wind speed observations across the wind turbine rotor layer along with meteorological and oceanographic data near the surface to characterize the wind resource. Lidar buoy data collected from recent deployments off the northern California coast near Humboldt County and the central California coast near Morro Bay allow for validation of commonly used reanalysis products. In this article, wind speeds from the Modern-Era Retrospective analysis for Research and Applications version 2 (MERRA-2), the Climate Forecast System version 2 (CFSv2), the North American Regional Reanalysis (NARR), the European Centre for Medium-Range Weather Forecasts Reanalysis version 5 (ERA5), and the analysis system of the Rapid Refresh (RAP) are validated at heights within the wind turbine rotor layer ranging from 50 to 100 m. The validation results offer guidance on the performance and uncertainty associated with utilizing reanalyses for offshore wind resource characterization, providing the offshore wind energy community with information on the conditions that lead to reanalysis error. At both California coast locations, the reanalyses tend to underestimate the observed rotor-level wind resource. Occasions of large reanalysis error occur in conjunction with stable atmospheric conditions, wind speeds associated with peak turbine power production (> 10 m s−1), and mischaracterization of the diurnal wind speed cycle in summer months.

17 WIND ENERGY↗

Climate Forcing Datasets for Agricultural Modeling: Merged Products for Gap-Filling and Historical Climate Series Estimation

The AgMERRA and AgCFSR climate forcing datasets provide daily, high-resolution, continuous, meteorological series over the 1980-2010 period designed for applications examining the agricultural impacts of climate variability and climate change. These datasets combine daily resolution data from retrospective analyses (the Modern-Era Retrospective Analysis for Research and Applications, MERRA, and the Climate Forecast System Reanalysis, CFSR) with in situ and remotely-sensed observational datasets for temperature, precipitation, and solar radiation, leading to substantial reductions in bias in comparison to a network of 2324 agricultural-region stations from the Hadley Integrated Surface Dataset (HadISD). Results compare favorably against the original reanalyses as well as the leading climate forcing datasets (Princeton, WFD, WFD-EI, and GRASP), and AgMERRA distinguishes itself with substantially improved representation of daily precipitation distributions and extreme events owing to its use of the MERRA-Land dataset. These datasets also peg relative humidity to the maximum temperature time of day, allowing for more accurate representation of the diurnal cycle of near-surface moisture in agricultural models. AgMERRA and AgCFSR enable a number of ongoing investigations in the Agricultural Model Intercomparison and Improvement Project (AgMIP) and related research networks, and may be used to fill gaps in historical observations as well as a basis for the generation of future climate scenarios.

Climate Forcing Data↗

Toward Coupled Data Assimilation in NASA’s GEOS: Developments in the Ocean Context

The Global Modeling & Assimilation Office (GMAO) at NASA GSFC produces analyses and predictions of the Earth system using various configurations of the Goddard Earth Observing System (GEOS) model and assimilation system. The current sub-seasonal-to-seasonal prediction system (GEOS-S2S) is based on a coupled atmosphere-ocean-land-ice configuration of GEOS which includes the Modular Ocean Model version 5 (MOM5) run at approximately 50-km resolution and a de-coupled OI-based ocean analysis that uses an initialization of MOM5 forced by the MERRA-2 reanalysis. GMAO will soon implement an updated GEOS-S2S system that will run at 25-km resolution and adopt aspects of the hybrid four-dimensional ensemble-variational (H4DEnVar) system already running in the production-version atmospheric analysis system, including a Local Ensemble Transform Kalman Filter (LETKF) to provide initial conditions for the oceanic state. This presentation will focus on developments to sustain the GMAO's systems on longer time horizons, where more radical transformations will be required to adapt to advanced computing environments, higher resolution and more diverse model components, and new observations for the Earth system. Results will describe progress toward a version of the GEOS coupled system that will be based around the Joint Effort for Data assimilation Integration (JEDI) framework being developed within Joint Center for Satellite Data Assimilation (JCSDA) and include an updated ocean model, MOM6. Discussion will focus specifically on the use of a Unified Forward Operator (UFO) for simulating observations and the Object Oriented Prediction System (OOPS) for providing the state estimate. These features are being developed as a multi-agency effort under the auspices of the JCSDA and are being adopted in the GMAO for all its applications of coupled data assimilation including S2S, numerical weather prediction, and reanalysis.

Mahajan, Rahul↗

Evaluating the Sensitivity of Agricultural Model Performance to Different Climate Inputs: Supplemental Material

Projections of future food production necessarily rely on models, which must themselves be validated through historical assessments comparing modeled and observed yields. Reliable historical validation requires both accurate agricultural models and accurate climate inputs. Problems with either may compromise the validation exercise. Previous studies have compared the effects of different climate inputs on agricultural projections but either incompletely or without a ground truth of observed yields that would allow distinguishing errors due to climate inputs from those intrinsic to the crop model. This study is a systematic evaluation of the reliability of a widely used crop model for simulating U.S. maize yields when driven by multiple observational data products. The parallelized Decision Support System for Agrotechnology Transfer (pDSSAT) is driven with climate inputs from multiple sources reanalysis, reanalysis that is bias corrected with observed climate, and a control dataset and compared with observed historical yields. The simulations show that model output is more accurate when driven by any observation-based precipitation product than when driven by non-bias-corrected reanalysis. The simulations also suggest, in contrast to previous studies, that biased precipitation distribution is significant for yields only in arid regions. Some issues persist for all choices of climate inputs: crop yields appear to be oversensitive to precipitation fluctuations but under sensitive to floods and heat waves. These results suggest that the most important issue for agricultural projections may be not climate inputs but structural limitations in the crop models themselves.

simulation↗

A High‐Resolution Tropical Mesoscale Convective System Reanalysis (TMeCSR)

Abstract Modern global reanalysis products have greatly accelerated meteorological research in synoptic‐to‐planetary‐scale phenomena. However, their use in studying tropical mesoscale convective systems (MCSs) and their regional‐to‐global impact has mostly been limited to supplying initial and boundary conditions for MCS‐resolving simulations and providing information about the large‐scale environments of MCSs. These limitations are due to difficulties in resolving tropical MCS dynamics in the relatively low‐resolution global models and that tropical MCSs often occur over poorly observed regions. In this work, a Tropical MCS‐resolving Reanalysis product (TMeCSR) was created over a region with frequent tropical MCSs. This region spans the tropical Indian Ocean, tropical continental Asia, Maritime Continent, and Western Pacific. TMeCSR is produced by assimilating all‐sky infrared radiances from geostationary satellites and other conventional observations into an MCS‐resolving regional model using the Ensemble Kalman Filter. The resulting observation‐constrained high‐resolution (9‐km grid spacing) data set is available hourly during the boreal summer (June‐August) of 2017, during which widespread severe flooding occurred. Comparisons of TMeCSR and European Center for Medium Range Weather Forecast Reanalysis version 5 (ERA5) against independent satellite retrievals indicate that TMeCSR's cloud and multiscale rain fields are better than those of ERA5. Furthermore, TMeCSR better captured the diurnal variability of rainfall and the statistical characteristics of MCSs. Forecasts initialized from TMeCSR also have more accurate rain and clouds than those initialized from ERA5. The TMeCSR and ERA5 forecasts have similar performances with respect to sounding and surface observations. These results indicate that TMeCSR is a promising MCS‐resolving data set for tropical MCS studies.

Chan, Man‐Yau↗

Evaluation of Global Surface Water Temperature Data Sets for use in Passive Remote Sensing of Soil Moisture

Inland open water bodies often pose a systematic error source in the passive remote sensing retrievals of soil moisture. Water temperature is a necessary variable used to compute water emissions that is required to be subtracted from satellite observation to yield actual emissions from the land portion, which in turn generates accurate soil moisture retrievals. Therefore, overestimation of soil moisture can often be corrected using concurrent water temperature data in the overall mitigation procedure. In recent years, several data sets of lake water temperature have become available, but their specifications and accuracy have rarely been investigated in the context of passive soil moisture remote sensing on a global scale. For this reason, three lake temperature products were evaluated against in-situ measurements from 2007 to 2011. The data sets include the lake surface water temperature (LSWT) from Global Observatory of Lake Responses to Environmental Change (GloboLakes), the Copernicus Global Land Operations Cryosphere and Water (C-GLOPS), as well as the lake mix-layer temperature (LMLT) from the European Centers for Medium-Range Weather Forecast (ECMWF) ERA5 Land Reanalysis. GloboLakes, C-GLOPS, and ERA5 Land have overall comparable performance with Pearson correlations (R) of 0.87, 0.92 and 0.88 in comparison with in-situ measurements. LSWT products exhibit negative median biases of −0.27 K (GloboLakes) and −0.31 K (C-GLOPS), whereas the median bias of LMLT is 1.56 K. When mapped from their respective native resolutions to a common 9 km Equal-Area Scalable Earth (EASE) Grid 2.0 projection, similar relative performance was observed. LMLT and LSWT data are closer in performance over the 9 km grid cells that exhibit a small range of lake cover fractions (0.05–0.5). Despite comparable relative performance, ERA5 Land shows great advantages in spatial coverage and temporal resolution. In summary, an integrated evaluation on data accuracy, long-term availability, global coverage, temporal resolution, and regular forward processing with modest data latency led us to conclude that LMLT from the ERA5 Land Reanalysis product represents the most optimal path for use in the development of a long-term soil moisture product.

Runze Zhang↗

More Far-Side Deep Moonquake Nests Discovered

As reported last year, we started to reanalyze the seismic data acquired from 1969 to 1977 with a network of stations established on the Moon during the Apollo mission. The reason for the reanalysis was because recent advances in computer technology make it possible to employ much more sophisticated analysis techniques than was possible previously. The primary objective of the reanalysis was to search for deep moonquakes on the far side of the Moon and, if found, to use them to infer the structure of the Moon's deep interior, including a possible central core. The first step was to identify any new deep moonquakes that escaped our earlier search by applying a combination of waveform cross-correlation and single-link cluster analysis, and then to see if any of them are from previously unknown nests of deep moonquakes. We positively identified 7245 deep moonquakes, more than a five-fold increase from the previous 1360. We also found at least 88 previously unknown deep-moonquake nests. The question was whether any of these newly discovered nets were on the far side of the Moon, and we now report that our analysis of the data indicates that some of them are indeed on the far side.

Nakamura, Y.↗

Towards an Indian Land Data Assimilation System (ILDAS): A Coupled Hydrologic-Hydraulic System for Water Balance Assessments

Effective management of water resources requires reliable estimates of land surface states and fluxes, including water balance components. But most land surface models run in uncoupled mode and do not produce river discharge at catchment scales to be useful for water resources management applications. Such integrated systems are also rare over India where hydrometeorological extremes have wreaked havoc on the economy and people. So, an Indian Land Data Assimilation System (ILDAS) with a coupled land surface and a hydrodynamic model has been developed and driven by multiple meteorological forcings (0.1°, daily) to estimate land surface states, channel discharge, and floodplain inundation. ILDAS benefits from an integrated framework as well as the largest suite of observation records collected over India and has been used to produce a reanalysis product for 1981–2021 using four forcing datasets, namely, Modern-Era Retrospective Analysis for Research and Applications, Version 2 (MERRA-2), Climate Hazards Group InfraRed Precipitation with Station data (CHIRPS), ECMWF’s ERA-5, and Indian Meteorological Department (IMD) gridded precipitation. We assessed the uncertainty and bias in these precipitation datasets and validated all major components of the terrestrial water balance, i.e., surface runoff, soil moisture, terrestrial water storage anomalies, evapotranspiration, and streamflow, against a combination of satellite and in situ observation datasets. Our assessment shows that ILDAS can represent the hydrological processes reasonably well over the Indian landmass with IMD precipitation showing the best relative performance. Evaluation against ESA-CCI soil moisture shows that MERRA-2 based estimates outperform the others, whereas ERA-5 performs best in simulating evapotranspiration when evaluated against MODIS ET. Evaluations against observed records show that CHIRPS-based estimates have the highest performance in reconstructing surface runoff and streamflow. Once operational, this system will be useful for supporting transboundary water management decision making in the region.

Indian Land Data Assimilation System (ILDAS)↗

Analyzing Error in Balloon-Based Wind Measurements Due to Spatial Separation Using the North American Regional Reanalysis

Weather balloons have been a longstanding asset to NASA and Aerospace meteorology. Balloons are used from launch vehicle design to day-of-launch operations. One of the most valuable assets from these balloons is wind data from the surface up to 30 kilometers. Due to aloft winds, the balloons may drift downrange of the launch site and vehicle flight path. A 2017 study found balloons at Kennedy Space Center (KSC) can drift as far as 200 kilometers from the launch site (Decker 2017). To obtain robust launch vehicle wind assessments, it is highly desirable to characterize the wind environment along the flight path. This study looks into the errors associated with spatial separation of wind measurements using the North American Regional Reanalysis (NARR).

Curtis, Nathan↗

Long-Term Ozone Variability and Trends from Reanalyses

Stratospheric ozone has a profound impact on radiation and chemistry over various spatial and temporal scales. The evolution of stratospheric ozone over the 21st century, however, is not well understood, especially in the lower stratosphere. Highly vertically resolved ozone data from satellite-borne limb sounders have proved to be invaluable for studying ozone in the middle and upper stratosphere but it was not until recently that these measurements were successfully incorporated in atmospheric reanalyses. Validation and comparison studies have demonstrated that the addition of observations from the Microwave Limb Sounder (MLS) on EOS (Earth Observing System) Aura greatly improved the quality of ozone fields in MERRA-2 (Modern-Era Retrospective analysis for Research and Applications, Version 2) making these assimilated data sets more useful for scientific research. In this presentation we demonstrate that multidecadal lower-stratospheric ozone variability and trends can be derived from NASA's MERRA-2 reanalysis ozone. In particular, the reanalysis ozone bias-corrected using a chemistry model simulation as a transfer function agrees very well with recently reprocessed long ozonesonde records. Ozone trends in the lower stratosphere will be discussed in the context of recent findings (Ball et al., 2018) and interpreted in connection with long-term circulation changes in the lower stratosphere. Next, we show that the use of ozone data retrieved from the next generation OMPS (Ozone Mapping Profiler Suite) instruments, including the OMPS Limb Profiler, can successfully extend the reanalyses into the future allowing comprehensive monitoring of global ozone and interpretation of its evolution during the critical period of expected ozone recovery and climate change from increasing concentration of greenhouse gases.

ozone↗

Ground Heat Flux Reconstruction Using Bayesian Uncertainty Quantification Machinery and Surrogate Modeling

Ground heat flux (G 0 ) is a key component of the land-surface energy balance of high-latitude regions. Despite its crucial role in controlling permafrost degradation due to global warming, G 0 is sparsely measured and not well represented in the outputs of global scale model simulation. In this study, an analytical heat transfer model is tested to reconstruct G 0 across seasons using soil temperature series from field measurements, Global Climate Model, and climate reanalysis outputs. The probability density functions of ground heat flux and of model parameters are inferred using available G 0 data (measured or modeled) for snow-free period as a reference. When observed G 0 is not available, a numerical model is applied using estimates of surface heat flux (dependent on parameters) as the top boundary condition. These estimates (and thus the corresponding parameters) are verified by comparing the distributions of simulated and measured soil temperature at several depths. Aided by state-of-the-art uncertainty quantification methods, the developed G 0 reconstruction approach provides novel means for assessing the probabilistic structure of the ground heat flux for regional permafrost change studies.

54 ENVIRONMENTAL SCIENCES↗