Search NASA⌕ Search

SEARCH · Search NASA

Results for “Variational data assimilation”

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 361 records · Page 20

Long-Period Tidal Variations in the Length of Day

A new model of long-period tidal variations in length of day is developed. The model comprises 80 spectral lines with periods between 18.6 years and 4.7 days, and it consistently includes effects of mantle anelasticity and dynamic ocean tides for all lines. The anelastic properties followWahr and Bergen; experimental confirmation for their results now exists at the fortnightly period, but there remains uncertainty when extrapolating to the longest periods. The ocean modeling builds on recent work with the fortnightly constituent, which suggests that oceanic tidal angular momentum can be reliably predicted at these periods without data assimilation. This is a critical property when modeling most long-period tides, for which little observational data exist. Dynamic ocean effects are quite pronounced at shortest periods as out-of-phase rotation components become nearly as large as in-phase components. The model is tested against a 20 year time series of space geodetic measurements of length of day. The current international standard model is shown to leave significant residual tidal energy, and the new model is found to mostly eliminate that energy, with especially large variance reduction for constituents Sa, Ssa, Mf, and Mt.

LOD tidal oscillations↗

Reducing Aerosol Forcing Uncertainty by Combining Models With Satellite and Within-The-Atmosphere Observations: A Three-Way Street

Aerosol forcing uncertainty represents the largest climate forcing uncertainty overall. Its magnitude has remained virtually undiminished over the past 20 years despite considerable advances in understanding most of the key contributing elements. Recent work has produced modest increases only in the confidence of the uncertainty estimate itself. This review summarizes the contributions toward reducing the uncertainty in the aerosol forcing of climate made by satellite observations, measurements taken within the atmosphere, as well as modeling and data assimilation. We adopt a more measurement-oriented perspective than most reviews of the subject in assessing the strengths and limitations of each; gaps and possible ways to fill them are considered. Currently planned programs supporting advanced, global-scale satellite and surface-based aerosol, cloud, and precursor gas observations, climate modeling, and intensive field campaigns aimed at characterizing the underlying physical and chemical processes involved, are all essential. But in addition, new efforts are needed: (a) to obtain systematic aircraft in situ measurements capturing the multi-variate probability distribution functions of particle optical, microphysical, and chemical properties (and associated uncertainty estimates), as well as co-variability with meteorology, for the major aerosol airmass types; (b) to conceive, develop, and implement a suborbital (aircraft plus surface-based) program aimed at systematically quantifying the cloud-scale microphysics, cloud optical properties, and cloud-related vertical velocities associated with aerosol-cloud interactions; and (c) to focus much more research on integrating the unique contributions of satellite observations, suborbital measurements, and modeling, to reduce the persistent uncertainty in aerosol climate forcing.

aerosol climate forcing↗

Shifts in bacterial traits under chronic nitrogen deposition align with soil processes in arbuscular, but not ectomycorrhizal-associated trees

Nitrogen (N) deposition increases soil carbon (C) storage by reducing microbial activity. These effects vary in soil beneath trees that associate with arbuscular (AM) and ectomycorrhizal (ECM) fungi. Variation in carbon C and N uptake traits among microbes may explain differences in soil nutrient cycling between mycorrhizal associations in response to high N loads, a mechanism not previously examined due to methodological limitations. Here, we used quantitative Stable Isotope Probing (qSIP) to measure bacterial C and N assimilation rates from an added organic compound, which we conceptualize as functional traits. As such, we applied a trait-based approach to explore whether variation in assimilation rates of bacterial taxa can inform shifts in soil function under chronic N deposition. We show taxon-specific and community-wide declines of bacterial C and N uptake under chronic N deposition in both AM and ECM soils. N deposition-induced reductions in microbial activity were mirrored by declines in soil organic matter mineralization rates in AM but not ECM soils. Our findings suggest C and N uptake traits of bacterial communities can predict C cycling feedbacks to N deposition in AM soils, but additional data, for instance on the traits of fungi, may be needed to connect microbial traits with soil C and N cycling in ECM systems. In conclusion, our study also highlights the potential of employing qSIP in conjunction with trait-based analytical approaches to inform how ecological processes of microbial communities influence soil functioning.

59 BASIC BIOLOGICAL SCIENCES↗

Determining the impact of a severe dry to wet transition on watershed hydrodynamics in California, USA with an integrated hydrologic model

With the onset of climate change, regions relied upon for water supply are increasingly subject to end-member fluctuations between periods of severe drought followed by extreme precipitation. The impacts of these extreme conditions on watershed hydrodynamics in water-resource sensitive regions such as California are unknown despite their great importance for resilience and water management purposes. Understanding these impacts requires high-resolution physically based models to capture sharp variations of topography, land use, wetting fronts, etc. Here, an integrated hydrologic model was used in a high-performance computing framework to study the complex nonlinear dynamics occurring at a representative Californian watershed. The Cosumnes Watershed, one of the last major rivers in California without a dam, offers a rare opportunity to isolate the effects of water management from climate extremes. Here, we show model validation with comparisons between model outputs and local measurements in addition to various satellite-based products including (1) Snow Water Equivalent (SWE) with Snow Data Assimilation System (SNODAS) and a reconstruction method by Bair and co-authors, (2) soil moisture with Soil Moisture Active Passive (SMAP), and (3) evapotranspiration (ET) with Mapping Evapotranspiration at high Resolution with Internalized Calibration (METRIC). To assess changes in hydrodynamic behavior during climate extremes and their transitions, a simulation spanning a recent drought followed by the highest precipitation year on record (2015–2017) is discussed. From these simulations, we are able to highlight regions that are the most sensitive to climate extremes, which depend on many factors including hydrologic connectivity, geology and topography. These analyses provide a better understanding of the physical phenomena occurring in the watershed, strengthening our knowledge of how the system may respond to extreme conditions which might become the “new normal.”

54 ENVIRONMENTAL SCIENCES↗

Assimilation of Gridded Terrestrial Water Storage Observations from GRACE into a Land Surface Model

Observations of terrestrial water storage (TWS) from the Gravity Recovery and Climate Experiment (GRACE) satellite mission have a coarse resolution in time (monthly) and space (roughly 150,000 km(sup 2) at midlatitudes) and vertically integrate all water storage components over land, including soil moisture and groundwater. Data assimilation can be used to horizontally downscale and vertically partition GRACE-TWS observations. This work proposes a variant of existing ensemble-based GRACE-TWS data assimilation schemes. The new algorithm differs in how the analysis increments are computed and applied. Existing schemes correlate the uncertainty in the modeled monthly TWS estimates with errors in the soil moisture profile state variables at a single instant in the month and then apply the increment either at the end of the month or gradually throughout the month. The proposed new scheme first computes increments for each day of the month and then applies the average of those increments at the beginning of the month. The new scheme therefore better reflects submonthly variations in TWS errors. The new and existing schemes are investigated here using gridded GRACE-TWS observations. The assimilation results are validated at the monthly time scale, using in situ measurements of groundwater depth and soil moisture across the U.S. The new assimilation scheme yields improved (although not in a statistically significant sense) skill metrics for groundwater compared to the open-loop (no assimilation) simulations and compared to the existing assimilation schemes. A smaller impact is seen for surface and root-zone soil moisture, which have a shorter memory and receive smaller increments from TWS assimilation than groundwater. These results motivate future efforts to combine GRACE-TWS observations with observations that are more sensitive to surface soil moisture, such as L-band brightness temperature observations from Soil Moisture Ocean Salinity (SMOS) or Soil Moisture Active Passive (SMAP). Finally, we demonstrate that the scaling parameters that are applied to the GRACE observations prior to assimilation should be consistent with the land surface model that is used within the assimilation system.

GRACE↗

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↗

Global Ocean Evaporation: How Well Can We Estimate Interannual to Decadal Variability?

Evaporation from the world's oceans constitutes the largest component of the global water balance. It is important not only as the ultimate source of moisture that is tied to the radiative processes determining Earth's energy balance but also to freshwater availability over land, governing habitability of the planet. Here we focus on variability of ocean evaporation on scales from interannual to decadal by appealing to three sources of data: the new MERRA-2 (Modern-Era Retrospective analysis for Research and Applications -2); climate models run with historical sea-surface temperatures, ice and atmospheric constituents (so-called AMIP experiments); and state-of-the-art satellite retrievals from the Seaflux and HOAPS (Hamburg Ocean-Atmosphere Parameters and Fluxes from Satellite) projects. Each of these sources has distinct advantages as well as drawbacks. MERRA-2, like other reanalyses, synthesizes evaporation estimates consistent with observationally constrained physical and dynamical models-but data stream discontinuities are a major problem for interpreting multi-decadal records. The climate models used in data assimilation can also be run with lesser constraints such as with SSTs and sea-ice (i.e. AMIPs) or with additional, minimal observations of surface pressure and marine observations that have longer and less fragmentary observational records. We use the new ERA-20C reanalysis produced by ECMWF embodying the latter methodology. Still, the model physics biases in climate models and the lack of a predicted surface energy balance are of concern. Satellite retrievals and comparisons to ship-based measurements offer the most observationally-based estimates, but sensor inter-calibration, algorithm retrieval assumptions, and short records are dominant issues. Our strategy depends on maximizing the advantages of these combined records. The primary diagnostic tool used here is an analysis of bulk aerodynamic computations produced by these sources and uses a first-order Taylor series analysis of wind speed, SST, near-surface stability and relative humidity variations around climatology to gauge the importance of these components. We find that the MERRA-2 evaporation record is strongly influenced by the availability of wind speed and humidity from passive microwave imagers beginning in the late 1980s as well as by the SST record. The trend over the period 1980 to present is nearly 10%. AMIP or the ERA-20C trends are much smaller. We find that ENSO-related signals involving both wind speed and thermodynamic variability remain the primary signal in the latter and are confirmed by satellite retrievals. We present uncertainty estimates based on the various data sources and discuss the implications for GEWEX water and energy budget science challenges.

Robertson, Franklin R.↗

The Use of MERRA-2 Near Surface Meteorology to Understand the Behavior of Planetary Boundary Layer Heights Derived from Wind Profiler Data over the US Great Plains

The atmospheric general circulation model (GCM) that underlies the MERRA-2 reanalysis includes a suite of physical parameterizations that describe the processes that occur in the planetary boundary layer (PBL). The data assimilation system assures that the atmospheric state variables used as input to these parameterizations are constrained to the best fit to all of the available observations. Many studies, however, have shown that the GCM-based estimates of MERRA-2 PBL heights are biased high, and so are not reliable for boundary layer studies.A 20-year record of PBL heights was derived from Wind Profiler (WP) backscatter data measured at a wide network of stations throughout the US Great Plains and has been validated against independent estimates. The behavior of these PBL heights shows geographical and temporal variations that are difficult to attribute to particular physical processes without additional information that are not part of the observational record.In the present study, we use information on physical processes from MERRA-2 to understand the behavior of the WP derived PBL heights. The behavior of the annual cycle of both MERRA-2 and WP PBL heights shows four classes of behavior: (i) canonical, characterized by a monthly progression in PBL height that follows the solar insolation, (ii) double peak, characterized by canonical behavior that is interrupted by a minimum in July, (iii) late peak, characterized by a suppressed heights in May and June, and return to canonical in July and August, and (iv) early peak where the PBL height rises with solar insolation but is suppressed later in the summer. The explanation for these behaviors and the relationship to local precipitation, temperature, sensible and latent heat fluxes, net radiation and aerosol load is articulated using information from MERRA-2.

Salmun, Haydee↗

Leafweb: Dataset in Support of Coupled Modeling of Photophysics, Photochemistry, and Biochemistry of Photosynthesis, December 2022 Release

This data set contains measurements of leaf gas exchange and Pulse-Amplitude Modulated (PAM) fluorometry of light, carbon dioxide (CO2), oxygen (O2), and temperature responses from 26 C3 photosynthetic and four C4 photosynthetic species measured by independent researchers in Canada, China, Finland, Netherlands, and United States in the field, garden, or greenhouse. Data were collected between 1987 and 2021, however data for individual species only provide coverage over a few hours to one year. Species include three lianas, three shrubs, two boreal deciduous trees, two boreal evergreen needle-leaf tree, three temperate deciduous trees, four tropical deciduous trees, three tropical evergreen trees, one C3 grass, three C4 grasses, and six crop varieties. These measurements are conducted according to standard protocols in gas exchange (Long and Bernacchi 2003) and PAM fluorometry (Baker 2008). The Scots pine (Pinus sylvestris) dataset contains one-year continuous fluorometry observations made at intervals of 10 or 30 minutes in the field under natural environments using Walz monitoring PAM. Measurements from all other 29 species include simultaneous PAM fluorometry and gas exchange observations. Among these 29 species, seven species were measured with light response curves only (i.e., light intensity varied systematically with ambient CO2 concentration controlled at a constant level, e.g., 400 ppm). All the other 22 species were measured with both the light response and CO2 response (i.e., ambient CO2 concentration varied systematically with light intensity controlled at a constant level, e.g., 1200 µmolm-2s-1). For most species, measurements were made with temperature controlled at ~ 25 degrees Celsius (C) with the exceptions of Scots pine (natural diurnal and seasonal variations), tomato cultivar Basket Vee (~ 21 degrees C), and cotton which contained temperature stress experiments (9 to 40 degrees C). All measurements were made at ambient O2 concentration except for the tomato cultivar Basket Vee and cotton which used two O2 levels (2 and 21%) and rice and tomato cultivar Growdena which used five O2 levels (2, 10, 21, 35, and 50%). This dataset supports the complete modeling of photosynthesis at the leaf scale, including photophysical, photochemical, and biochemical reactions. It will support research in the applications of remotely sensed solar induced chlorophyll fluorescence and the responses of energy dissipation pathways to environmental variations across species and climates. This dataset contains 260 data files in comma separate (.csv) format contained inside 11 compressed (*.zip) files. Additional metadata are provided: one data dictionary, a species information table, and a file-level metadata file in comma separate (.csv) format and a user guide in PDF (*.pdf) format.

air temperature↗

The Reanalysis for Stratospheric Trace Gas Studies

The "Reanalysis for Stratospheric Trace Gas Studies" (ReSTS) project intends to provide a state-of-the-art meteorological dataset for use in chemistry-transport models (CTMs). The initial focus is on the period from May 1991 until April 1995, which encompasses the initial few years of observations from NASA's Upper Atmosphere Research Satellite (UARS) as well as the abrupt increase and slow decay of volcanic aerosol loading from Mount Pinatubo. The reanalysis is being performed using NASA's operational "GEOS-4" data assimilation system. The results presented in this paper will begin with an overview of the meteorology of the ReSTS period, validating the ReSTS data against independent analyses, followed by an evaluation of the system statistics (the "Observation minus Forecast" for different input data types). Some sensitivities to the assumptions made in the assimilation process will be presented, including aspects of the forecast model and the statistical model used in the observation-forecast merging. Aspects of trace-gas transport will be discussed, focussing on the on-line water vapor distribution (which is a quasi-inert trace gas in the stratosphere), off-line transport studies of ozone and selected long-lived species, as well as assimilated ozone. The results will be discussed in the context of transport studies and our ability to use reanalyzed datasets to examine long-term variations in the meteorology and composition of the middle atmosphere.

Pawson, Steven↗

Warm and Saline Events Embedded in the Meridional Circulation of the Northern North Atlantic

Ocean state estimates from 1958 to 2005 from the Simple Ocean Assimilation System (SODA) system are analyzed to understand circulation between subtropical and subpolar Atlantic and their connection with atmospheric forcing. This analysis shows three periods (1960s, around 1980, and 2000s) with enhanced warm, saline waters reaching high latitudes, alternating with freshwater events originating at high latitudes. It complements surface drifter and altimetry data showing the subtropical -subpolar exchange leading to a significant temperature and salinity increase in the northeast Atlantic after 2001. The warm water limb of the Atlantic meridional overturning cell represented by SODA expanded in density/salinity space during these warm events. Tracer simulations using SODA velocities also show decadal variation of the Gulf Stream waters reaching the subpolar gyre and Nordic seas. The negative phase of the North Atlantic Oscillation index, usually invoked in such variability, fails to predict the warming and salinization in the early 2000s, with salinities not seen since the 1960s. Wind stress curl variability provided a linkage to this subtropical/subpolar gyre exchange as illustrated using an idealized two ]layer circulation model. The ocean response to the modulation of the climatological wind stress curl pattern was found to be such that the northward penetration of subtropical tracers is enhanced when amplitude of the wind stress curl is weaker than normal. In this case both the subtropical and subpolar gyres weaken and the subpolar density surfaces relax; hence, the polar front moves westward, opening an enhanced northward access of the subtropical waters in the eastern boundary current.

Hakkinen, Sirpa↗

Pristine highland clasts in consortium breccia 14305 Petrology and geochemistry

Data are presented on the petrography and mineral chemistry of six pristine highland clasts chipped from the polymict lunar breccia 14305. Major and trace elements in the clasts were determined by instrumental neutron activation analysis, and mineral analyses were performed by electron microprobe. Mg-suite clasts have 'eastern' geochemical affinities, reaffirming the importance of local variations in geochemistry. These local variations are superimposed on the moon-wide, longitudinal variations noted by Warren and Wasson (1980). Alkali anorthosites and Mg-suite troctolites and anorthosites are not comagmatic, and cannot be related to a single parent magma by either fractional crystallization or variable assimilation of KREEP. Both magma suites may have assimilated varied amounts of KREEP into distinct parent magmas. Alternatively, alkali anorthosites may have crystallized directly from a KREEP-basalt parent magma. A thick crust of ferroan anorthosite probably never existed on the western lunar nearside, or was removed by basin-forming impacts prior to intrusion of later plutonic suites.

Shervais, J. W.↗

Joint Assimilation of SMOS Brightness Temperature and GRACE Terrestrial Water Storage Observations for Improved Soil Moisture Estimation

Observations from recent soil moisture missions (e.g. SMOS) have been used in innovative data assimilation studies to provide global high spatial (i.e. 40 km) and temporal resolution (i.e. 3-days) soil moisture profile estimates from microwave brightness temperature observations. In contrast with microwave-based satellite missions that are only sensitive to near-surface soil moisture (0 - 5 cm), the Gravity Recovery and Climate Experiment (GRACE) mission provides accurate measurements of the entire vertically integrated terrestrial water storage column but, it is characterized by low spatial (i.e. 150,000 km2) and temporal (i.e. monthly) resolutions. Data assimilation studies have shown that GRACE-TWS primarily affects (in absolute terms) deeper moisture storages (i.e., groundwater). This work hypothesizes that unprecedented soil water profile accuracy can be obtained through the joint assimilation of GRACE terrestrial water storage and SMOS brightness temperature observations. A particular challenge of the joint assimilation is the use of the two different types of measurements that are relevant for hydrologic processes representing different temporal and spatial scales. The performance of the joint assimilation strongly depends on the chosen assimilation methods, measurement and model error spatial structures. The optimization of the assimilation technique constitutes a fundamental step toward a multi-variate multi-resolution integrative assimilation system aiming to improve our understanding of the global terrestrial water cycle.

Soil Moisture↗

Carbon-Water-Energy Relations for Selected River Basins

A biophysical process-based model was run using satellite, assimilated and ancillary data for four years (1987-1990) to calculate components of total evaporation (transpiration, interception, soil and snow evaporation), net radiation, absorbed photosynthetically active radiation and net primary productivity over the global land surface. Satellite observations provided fractional vegetation cover, solar and photosynthetically active radiation incident of the surface, surface albedo, fractional cloud cover, air temperature and vapor pressure. The friction velocity and surface air pressure are obtained from a four dimensional data assimilation results, while precipitation is either only surface observations or a blended product of surface and satellite observations. All surface and satellite data are monthly mean values; precipitation has been disaggregated into daily values. All biophysical parameters of the model are prescribed according to published records. From these global land surface calculations results for river basins are derived using digital templates of basin boundaries. Comparisons with field observations (micrometeorologic, catchment water balance, biomass production) and atmospheric water budget analysis for monthly evaporation from six river basins have been done to assess errors in the calculations. Comparisons are also made with previous estimates of zonal variations of evaporation and net primary productivity. Efficiencies of transpiration, total evaporation and radiation use, and evaporative fraction for selected river basins will be presented.

Choudhury, B. J.↗

Mission Simulation of Space Lidar Measurements for Seasonal and Regional CO2 Variations

Results of mission simulation studies are presented for a laser-based atmospheric [82 sounder. The simulations are based on real-time carbon cycle process modeling and data analysis. The mission concept corresponds to the Active Sensing of [82 over Nights, Days, and Seasons (ASCENDS) recommended by the US National Academy of Sciences Decadal Survey of Earth Science and Applications from Space. One prerequisite for meaningful quantitative sensor evaluation is realistic CO2 process modeling across a wide range of scales, i.e., does the model have representative spatial and temporal gradients? Examples of model comparison with data will be shown. Another requirement is a relatively complete description of the atmospheric and surface state, which we have obtained from meteorological data assimilation and satellite measurements from MODIS and [ALIPS0. We use radiative transfer model calculations, an instrument model with representative errors ' and a simple retrieval approach to complete the cycle from "nature" run to "pseudo-data" CO2, Several mission and instrument configuration options are examined/ and the sensitivity to key design variables is shown. We use the simulation framework to demonstrate that within reasonable technological assumptions for the system performance, relatively high measurement precision can be obtained, but errors depend strongly on environmental conditions as well as instrument specifications. Examples are also shown of how the resulting pseudo - measurements might be used to address key carbon cycle science questions.

Kawa, Stephan↗

SHIELDS Final Technical Report

Predicting variations in the near-Earth space environment that can lead to spacecraft damage and failure, i.e. “space weather”, remains a big space physics challenge. A new capability was developed at Los Alamos National Laboratory (LANL) to understand, model, and predict Space Hazards Induced near Earth by Large Dynamic Storms, the SHIELDS framework. This framework simulates the dynamics of the Surface Charging Environment (SCE), the hot (keV) electrons representing the source and seed populations for the radiation belts, on both macro- and micro-scale. In addition to using physics-based models (like RAM-SCB, BATS-R-US, and iPIC3D), new data assimilation techniques employing data from LANL instruments on the Van Allen Probes and geosynchronous satellites were developed. An order of magnitude improvement in the accuracy in the simulation of the spacecraft surface charging environment was thus obtained. SHIELDS also includes a post-processing tool designed to calculate the surface charging for specific spacecraft geometry using the Curvilinear Particle-In-Cell (CPIC) code and to evaluate anomalies' relation to SCE dynamics. Such diagnostics is critically important when performing forensic analyses of space-system failures.

79 ASTRONOMY AND ASTROPHYSICS↗

Diagnostic Studies with GLA Fields

Assessments of the NASA Goddard Earth Observing System-1 Data Assimilation System(GEOS-1 DAS), regarding heating rates, energetics, and angular momentum quantities were made. These diagnostics can be viewed as measures of climate variability. Comparisons with the NOAA/NCEP reanalysis system of momentum and energetics diagnostics are included. Water vapor and angular momentum are diagnosed in many models, including those of NASA, as part of the Atmospheric Model Intercomparison Project. 'Me GEOS-I and NOAA/NCEP global atmospheric angular momentum values are coherent on time scales down to about three days. Furthermore, they agree with the series of Earth angular momentum, as measured by tiny fluctuations in the rotation rate of the Earth, as variations in the length of day. The torques that effect such changes in atmospheric and Earth momentum are dominated by the influence of particular mountain systems, including the Rockies, Himalayas, and Andes, upon mountain torques on time scales shorter than about two weeks. Other project areas included collaboration with Goddard Space Flight Center to examine the impact of mountainous areas and the treatments of parameterizations on diagnoses of the atmosphere. Relevant preprints are included herein.

Salstein, David A.↗

The Use of MERRA-2 Near Surface Meteorology to Understand the Behavior of Planetary Boundary Layer Heights Derived from Wind Profiler Data over the US Great Plains

The atmospheric general circulation model (GCM) that underlies the MERRA-2 reanalysis includes a suite of physical parameterizations that describe the processes that occur in the planetary boundary layer (PBL). The data assimilation system assures that the atmospheric state variables used as input to these parameterizations are constrained to the best fit to all of the available observations. Many studies, however, have shown that the GCM-based estimates of MERRA-2 PBL heights are biased high, and so are not reliable for application related to constituent transport or the carbon cycle. A new 20-year record of PBL heights was derived from Wind Profiler (WP) backscatter data measured at a wide network of stations throughout the US Great Plains and has been validated against independent estimates. The behavior of these PBL heights shows geographical and temporal variations that are difficult to attribute to particular physical processes without additional information that are not part of the observational record. In the present study, we use information on physical processes from MERRA-2 to understand the behavior of the WP derived PBL heights. The behavior of the annual cycle of both MERRA-2 and WP PBL heights shows three classes of behavior: (i) canonical, where the annual cycle follows the annual cycle of the sun, (ii) delayed, where the PBL height reaches its annual maximum after the annual maximum of the solar insolation, and (iii) double maxima, where the PBL height begins to rise with the solar insolation but falls sometimes during the summer and then rises again. Although the magnitude of these types of variations is described by the WP PBL record, the explanation for these behaviors and the relationship to local precipitation, temperature, hydrology and sensible and latent heat fluxes is articulated using information from MERRA-2.

MERRA↗