Search NASA⌕ Search

SEARCH · Search NASA

Results for “global modeling”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS2004

Aerosol Absorption in the Atmosphere: Perspectives from Global Model, Ground-Based Measurements, and Field Observations

Aerosol absorption in the atmosphere poses a major uncertainty in assessing the aerosol climate effects. This uncertainty arises from the poorly quantified aerosol sources, especially black carbon emissions, and our limited knowledge of aerosol mixing state and optical properties. Here we use a global model GOCART to simulate atmospheric aerosols, including sulfate, black carbon, organic carbon, dust, and sea salt. We compare the model calculated total aerosol optical thickness, extinction, and absorption with those quantities from the ground-based sun photometer measurements from AERONET, satellite retrievals from the TOMS instrument, and field observations from ACE-Asia. We will examine the most sensitive factors in determining the aerosol absorption. and the consequences of assessing the aerosol radiative forcing and atmospheric heating associated with those factors.

Chin, Mian↗
NASA NTRSDate not supplied

Evaluation of Planetary Boundary Layer Structure from NASA Global Modeling and Assimilation Office’s Next Retrospective Analysis Product MERRA-21C

The Planetary Boundary Layer (PBL) is a complex interface between Earth’s surface and the atmosphere with high spatiotemporal variability in its characteristics, and accurate simulation and observation of the PBL has proven to be a challenge. In this study, we evaluate PBL structures from NASA Global Modeling and Assimilation Office (GMAO)’s next retrospective analysis product: the Goddard Earth Observing System Retrospective Analysis of the early 21st Century (GEOS-R21C), with 25-km horizontal resolution. The GEOS-R21C incorporates a wide range of observing systems and various improvements over previous GMAO reanalysis. We compare the PBL thermodynamic structure as well as PBL Height (PBLH) from GEOS-R21C with various PBL observations. For PBLH, definitions are different depending on each observing system. PBLHs derived from radiosonde observations and Global Navigation Satellite System Radio Occultation (GNSS RO) are based on bulk Richardson number and refractivity gradient, respectively. The results of the evaluation of PBL structure from GMAO’s next retrospective analysis product will be presented. In addition, comparison between near-real-time (NRT) commercial and NASA commercial GNSS RO PBLH will also be discussed.

Eun-Gyeong Yang↗
NASA NTRS2001

Global Modeling of Tropospheric Chemistry with Assimilated Meteorology: Model Description and Evaluation

We present a first description and evaluation of GEOS-CHEM, a global three-dimensional (3-D) model of tropospheric chemistry driven by assimilated meteorological observations from the Goddard Earth Observing System (GEOS) of the NASA Data Assimilation Office (DAO). The model is applied to a 1-year simulation of tropospheric ozone-NOx-hydrocarbon chemistry for 1994, and is evaluated with observations both for 1994 and for other years. It reproduces usually to within 10 ppb the concentrations of ozone observed from the worldwide ozonesonde data network. It simulates correctly the seasonal phases and amplitudes of ozone concentrations for different regions and altitudes, but tends to underestimate the seasonal amplitude at northern midlatitudes. Observed concentrations of NO and peroxyacetylnitrate (PAN) observed in aircraft campaigns are generally reproduced to within a factor of 2 and often much better. Concentrations of HNO3 in the remote troposphere are overestimated typically by a factor of 2-3, a common problem in global models that may reflect a combination of insufficient precipitation scavenging and gas-aerosol partitioning not resolved by the model. The model yields an atmospheric lifetime of methylchloroform (proxy for global OH) of 5.1 years, as compared to a best estimate from observations of 5.5 plus or minus 0.8 years, and simulates H2O2 concentrations observed from aircraft with significant regional disagreements but no global bias. The OH concentrations are approximately 20% higher than in our previous global 3-D model which included an UV-absorbing aerosol. Concentrations of CO tend to be underestimated by the model, often by 10-30 ppb, which could reflect a combination of excessive OH (a 20% decrease in model OH could be accommodated by the methylchloroform constraint) and an underestimate of CO sources (particularly biogenic). The model underestimates observed acetone concentrations over the South Pacific in fall by a factor of 3; a missing source from the ocean may be implicated.

Bey, Isabelle↗
NASA NTRS2019

MERRA-2 Ocean: The NASA Global Modeling and Assimilation Office's Weakly Coupled Atmosphere-Ocean Reanalysis Using GEOS-S2S Version 3

The NASA Modern Era Reanalysis for Research and Applications (MERRA2) has been a respected and widely used reanalysis that has so far been restricted to the atmosphere. Now a newly released version of the atmosphere/ocean coupled data assimilation system (AODAS) has been developed by the NASA/Goddard Global Modeling and Assimilation Office to perform a retrospective ocean reanalysis from 1982 to present. In addition to assimilating all available in situ data (e.g. Argo, mooring, XBT and CTD data) and altimetry information into the ocean, the new version (GEOS-S2S Version 3) model includes a higher resolution, eddy-permitting ocean model than previous versions, a more realistic implementation of the atmosphere-ocean interface layer, and an improved coupling between glacier and ocean (among other improvements). In addition, this ocean data assimilation was expanded to include the assimilation of satellite sea surface salinity. The MERRA-2 AODAS will be described, and preliminary results will be shown from the assimilation reanalysis and from retrospective forecasts issued using a new ensemble strategy. Following the Global Ocean Data Assimilation Experiment (GODAE) protocols, we will present Class 1 through Class 4 validation results from the ocean reanalysis. Results indicate an improved ocean mixed layer depth, improved salinity near Greenland, an improved diurnal cycle of the sea surface skin temperature, an improved estimate of ocean evaporation, and better representation of western boundary currents (e.g. Gulf Stream) from our new ocean reanalysis. One of the motivations of this project is to provide optimal initial states for ENSO forecasting. Therefore, we will also present some preliminary results of retrospective ENSO forecasts. After thorough testing, it is expected that the GEOS-S2S Version 3 will replace our contributions to North American Multi-Model Ensemble (NMME), WCRP Subseasonal to Seasonal (S2S), and IRI seasonal prediction forecast projects.

Molod, Andrea↗
NASA NTRS1989

The Southern Hemisphere lower stratosphere during August and September 1987 - Analyses based on the United Kingdom Meteorological Office Global Model

This paper presents a series of meteorological analyses used to aid the interpretation of the in situ Airborne Antarctic Ozone Experiment (AAOE) observations obtained aboard the ER-2 and DC-8 aircraft and examines the basis and accuracy of the analytical procedure. Maps and sections of meteorological variables derived from the UK Meteorological Office Global Model are presented for ER-2 and DC-8 flight days. It is found that analyzed temperatures and winds are generally in good agreement with AAOE observations at all levels; minor discrepancies were evident only at DC-8 altitudes. Maps of potential vorticity presented on the 428-K potential temperature surface show that the vortex is essentially circumpolar, although there are periods when major distortions are apparent.

Mckenna, D. S.↗
NASA NTRS2015

The Quick Fire Emissions Dataset (QFED): Documentation of Versions 2.1, 2.2 and 2.4: Technical Report Series on Global Modeling and Data Assimilation - Volume 38

Biomass burning is an important source of particulates and trace gases and a major element of the terrestrial carbon cycle. Well constrained emissions from vegetation fires are needed to model direct and indirect effects of biomass burning aerosols, to model homogeneous and heterogeneous chemistry in the atmosphere, and to perform credible Earth system analysis, and climate and air pollution studies. To improve the performance of NASA Goddard Earth Observing System Model (GEOS) in the areas of atmospheric constituent modeling with a focus on biomass burning we developed the Quick Fire Emissions Dataset (QFED). The QFED emissions are based on the fire radiative power (top-down) approach and draw on the cloud correction method developed in the Global Fire Assimilation System (GFAS). Location and fire radiative power of fires are obtained from the Moderate Resolution Imaging Spectroradiometer (MODIS) Level 2 fire products (MOD14 and MYD14) and the MODIS Geolocation products (MOD03 and MYD03). QFED strengths are high spatial and temporal resolutions and near-real time availability. Daily mean emissions are available at 0.3125 times 0.25 degrees and in recent versions also at 0.1 times 0.1 degrees. QFED provides emissions of black carbon, organic carbon, sulfur dioxide, carbon monoxide, carbon dioxide, PM2.5, ammonia, nitrogen oxides, methyl ethyl ketone, propylene, ethane, propane, n- and i-butane, acetaldehyde, formaldehyde, acetone and methane. Two QFED product systems are maintained by the NASA Global Modeling and Assimilation Office (GMAO): one that produces near real-time daily emissions used operationally in the GEOS-5 Data Assimilation System, and one that produces an extended historical dataset with daily emissions from March 2000 to the present. The historical dataset also provides monthly mean and monthly climatological emissions.

Koster, Randal D.↗
NASA NTRS2007

Simulations of Tropospheric NO2 by the Global Modeling Initiative (GMI) Model Utilizing Assimilated and Forecast Meteorological Fields: Comparison to Ozone Monitoring Instrument (OMI) Measurements

We present simulations of the tropospheric composition for the years 2004 and 2005, carried out by the GMI Combined Stratosphere-Troposphere (Combo) model, at a resolution of 2degx2.5deg. The model includes a new parameterization of lightning sources of NO(x) which is coupled to the cloud mass fluxes in the adopted meteorological fields. These simulations use two different sets of input meteorological fields: a)late-look assimilated fields from the Global Modeling and Assimilation Office (GMAO), GEOS-4 system and b) 12-hour forecast fields initialized with the assimilated data. Comparison of the forecast to the assimilated fields indicates that the forecast fields exhibit less vigorous convection, and yield tropical precipitation fields in better agreement with observations. Since these simulations include a complete representation of the stratosphere, they provide realistic stratosphere-tropospheric fluxes of O3 and NO(y). Furthermore, the stratospheric contribution to total columns of different troposheric species can be subtracted in a consistent fashion, and the lightning production of NO(y) will depend on the adopted meteorological field. We concentrate here on the simulated tropospheric columns of NO2, and compare them to observations by the OM1 instrument for the years 2004 and 2005. The comparison is used to address these questions: a) is there a significant difference in the agreement/disagreement between simulations for these two different meteorological fields, and if so, what causes these differences?; b) how do the simulations compare to OMI observations, and does this comparison indicate an improvement in simulations with the forecast fields? c) what are the implications of these simulations for our understanding of the NO2 emissions over continental polluted regions?

Rodriquez, J. M.↗
NASA NTRS2023

Simulations of 7 Be and 10 Be with the GEOS-Chem global model v14.0.2 using state-of-the-art production rates

The cosmogenic radionuclides 7 Be and 10 Be are useful aerosol tracers for atmospheric transport studies. Combining 7 Be and 10 Be measurements with an atmospheric transport model can not only improve our understanding of the radionuclide transport and deposition processes but also provide an evaluation of the transport process in the model. To simulate these aerosol tracers, it is critical to evaluate the influence of radionuclides production uncertainties on simulations. Here we use the GEOS-Chem chemical transport model driven by the MERRA-2 reanalysis to simulate 7 Be and 10 Be with different production scenarios: the default production rate in GEOS-Chem based on an empirical approach (Lal and Peters, 1967; denoted as LP67), and two production rates from the CRAC:Be (Cosmic Ray Atmospheric Cascade: Beryllium; Poluianov et al., 2016) model considering only geomagnetic cut-off rigidities for a geocentric axial dipole (denoted as P16) or realistic spatial geomagnetic cut-off rigidity variations due to non-dipole moments of the geomagnetic field (denoted as P16spa). The model results are comprehensively evaluated with a large number of measurements including surface air concentrations and deposition fluxes. The model with the P16spa production can reproduce the absolute values and temporal variability of 7 Be and 10 Be surface concentrations and deposition fluxes on annual and sub-annual scales, as well as the vertical profiles of air concentrations. Simulations with the LP67 production tend to overestimate the absolute values of 7 Be and 10 Be concentrations. The P16 simulations suggest less than 10% differences compared to P16spa but tend to produce a significant positive bias (>20%) in the 7 Be deposition fluxes over East Asia. We find that the deposition fluxes are more sensitive to the production in the troposphere and downward transport from the stratosphere. Independent of the production models, surface air concentrations and deposition fluxes from all simulations show similar seasonal variations, suggesting a dominant meteorological influence. The model can also simulate reasonably the stratosphere-troposphere exchange process of 7 Be and 10 Be by producing stratospheric contribution and 10 Be/ 7 Be ratio values that agree with measurements. Finally, we illustrate the importance of including the time-varying solar modulation in the production calculation, which can significantly improve the agreement between model results and measurements, especially at mid- and high- latitudes. Reduced uncertainties in the production rates, as demonstrated in this study, improve the utility of 7 Be and 10 Be as aerosol tracers for evaluating and testing transport and scavenging processes in global models.

Minjie Zheng↗
NASA NTRS2008

Magnetotail Current Sheet Thinning and Magnetic Reconnection Dynamics in Global Modeling of Substorms

Magnetotail current sheet thinning and magnetic reconnection are key elements of magnetospheric substorms. We utilized the global MHD model BATS-R-US with Adaptive Mesh Refinement developed at the University of Michigan to investigate the formation and dynamic evolution of the magnetotail thin current sheet. The BATSRUS adaptive grid structure allows resolving magnetotail regions with increased current density up to ion kinetic scales. We investigated dynamics of magnetotail current sheet thinning in response to southwards IMF turning. Gradual slow current sheet thinning during the early growth phase become exponentially fast during the last few minutes prior to nightside reconnection onset. The later stage of current sheet thinning is accompanied by earthward flows and rapid suppression of normal magnetic field component $B-z$. Current sheet thinning set the stage for near-earth magnetic reconnection. In collisionless magnetospheric plasma, the primary mechanism controlling the dissipation in the vicinity of the reconnection site is non-gyrotropic effects with spatial scales comparable with the particle Larmor radius. One of the major challenges in global MHD modeling of the magnetotail magnetic reconnection is to reproduce fast reconnection rates typically observed in smallscale kinetic simulations. Bursts of fast reconnection cause fast magnetic field reconfiguration typical for magnetospheric substorms. To incorporate nongyritropic effects in diffusion regions we developed an algorithm to search for magnetotail reconnection sites, specifically where the magnetic field components perpendicular to the local current direction approaches zero and form an X-type configuration. Spatial scales of the diffusion region and magnitude of the reconnection electric field are calculated self-consistently using MHD plasma and field parameters in the vicinity of the reconnection site. The location of the reconnection sites and spatial scales of the diffusion region are updated during the simulations. Such an approach allows quantifying the interaction between large-scale global magnetospheric dynamics and microphysical processes in diffusion regions localized near reconnection sites. To clarify the role of smallscale non-MHD effects in diffusion region on the global magnetospheric dynamic and to test different models of dissipation we perform simulations with steady southward IMF driving of the magnetosphere.

Kuznetsova, M. M.↗
NASA NTRS1983

A global model of the neutral thermosphere in magnetic coordinates based on AE-C data

Molecular nitrogen, atomic oxygen, and helium densities obtained from the AE-C satellite are analyzed in magnetic latitude and magnetic local time coordinates and compared to OGO 6 data for various seasons and magnetic activity levels. A depletion region for atomic oxygen and helium with respect to molecular nitrogen at high magnetic latitudes in the postmidnight magnetic local time sector persists under both high and low solar activity conditions. A global model in magnetic coordinates, patterned after the mass spectrometer and incoherent scatter (MSIS) model, is developed for molecular nitrogen, atomic oxygen, and helium. It is shown to represent the data well, without any residual UT dependence, and with an accuracy comparable to that provided by the more complex MSIS model with longitude terms The advantage of using magnetic coordinates is that they are more directly related to the major energy inputs and momentum sources in the polar regions than are geographic coordinates and are more convenient to use in studies of high latitude energy deposition processes. This is important for comparison with theoretical models where the number of coordinates is limited.

Stehle, C. G.↗
NASA NTRS1984

Three dimensional global modeling of atmospheric CO2

The initial attempts to model the atmospheric CO2 distribution, including couplings to the ocean and biosphere as sources and sinks of atmospheric CO2, encourage the notion that this approach will lead to useful quantitative constraints on CO2 fluxes. Realization of this objective will require: (1) continued improvement in the realism of the global transport modeling; (2) extended timeline of atmospheric CO2 monitoring, which improved precision and improved definition of the uncertainties in the measured CO2 amounts; and (3) given an accurate knowledge of model capabilities and limitations and given a good understanding of CO2 observations and their limitations, there is a need for good ideas concerning what quantitative information on the carbon cycle can be inferred from global modeling.

Hanse, J.↗
NASA NTRS1988

Magnetohydrodynamic boundary conditions for global models

Boundary conditions in the ionosphere and the upstream solar wind are important in determining the dynamics of global magnetohydrodynamic models of the magnetosphere. It is generally recognized that the orientation of the magnetic field in the upstream solar wind strongly modulates the rate of energy input into the magnetosphere by magnetic reconnection. However, other aspects of the upstream boundary conditions may determine whether the reconnection occurs in a patchy manner, as in flux transfer events, or in a global manner, as in the Paschmann et al. (1979) events. Ionospheric boundary conditions should also affect the reconnection process. For example, ionospheric line-tying can cause x-line motion in the outer magnetosphere. If it is assumed that auroras occur on field lines mapping to x-lines, then auroral motions are different than the local convective motion of the plasma in which they occur. Global magnetohydrodynamic models which incorporate both magnetospheric reconnection and ionospheric convection could be used to investigate the effect of reconnection and convection upon dayside and nightside auroral motions during the course of a magnetic substorm.

Forbes, T. G.↗
NASA NTRS2012

An Overview of Atmospheric Composition OSSE Activities at NASA's Global Modeling and Assimilation Office

A model-based Observing System Simulation Experiment (OSSE) is a framework for numerical experimentation in which observables are simulated from fields generated by an earth system model, including a parameterized description of observational error characteristics. Simulated observations can be used for sampling studies, quantifying errors in analysis or retrieval algorithms, and ultimately being a planning tool for designing new observing missions. While this framework has traditionally been used to assess the impact of observations on numerical weather prediction, it has a much broader applicability, in particular to aerosols and chemical constituents. In this talk we will give a general overview of Observing System Simulation Experiments (OSSE) activities at NASA's Global Modeling and Assimilation Office, with focus on its emerging atmospheric composition component.

daSilva, Arlinda↗
DOE OSTI2024

Insights into Sustainable Nitrogen Fixation by Gas-phase Spectroscopic Measurements and Global Modeling of Reaction Intermediates in Humid Nitrogen Plasma

In this article, there is growing interest in reacting molecular nitrogen and water to sustainably synthesize fixed forms of nitrogen such as ammonia. In particular, low-temperature plasmas can activate these relatively inert feedstocks at or near room temperature without a catalyst. However, because of the enhanced reactivity and nonequilibrium chemistry, a diverse range of products is formed, and the underlying reaction mechanisms are exceedingly complex. In this work, we studied a simplified reactor consisting of a gaseous plasma containing controlled mixtures of nitrogen gas and water vapor. Densities of key chemical species such as N, H, OH, NH, and NO were measured by emission and laser-based spectroscopy as a function of the relative humidity. A global model was constructed and the reaction network was validated by comparing calculated species densities with experiments. We discover that N, a key initial intermediate for ammonia, strongly decreases in the presence of water vapor, and as a result, ammonia formation becomes limited at high relative humidity. This decrease is surprisingly not because N itself reacts, but because one of its main sources, an excited molecular nitrogen state, is reacted away. In addition, oxidation pathways for nitrogen, which lead to NO and related products, are found to be favored over reduction pathways because the corresponding reverse reactions are less significant. Together, this understanding helps explain previously reported observations of selectivity toward nitrogen oxides over ammonia, particularly at higher relative humidities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2019

A Global Model for Estimating Atmospheric Phase Scintillation Statistics

Since 2007, the National Aeronautics and Space Administration (NASA) has been collecting atmospheric phase turbulence data from various NASA ground stations throughout the world. The goal of these measurement campaigns has been to generate statistics to characterize the local site turbulence conditions and their impact on widely distributed ground based antenna arrays. This is of critical importance for the situation of uplink arraying, in which a priori knowledge of the fast varying turbulent conditions of water vapor in the troposphere may not be known, and will impact the power combining efficiency of ground based transmitting arrays. Therefore, the design of these type of systems will be dependent on the local climatology of the particular ground station site. Based on the 30+ station years of data collected characterizing atmospheric phase scintillation statistics at various sites, a global model is presented which attempts to predict the average phase statistics of a generic site based on local surface weather data, such as surface pressure, temperature, relative humidity, wind speed, and median wind direction. A model is proposed and based on a standard log power distribution similar to amplitude scintillation models trained on the existing data sets and shows reasonable accuracy against existing data sets.

propagation↗
NASA NTRS2018

Synergistic and Collaborative Development Strategies for FV3 Powered Next Generation Unified Global Modeling System

The GFDL (Geophysical Fluid Dynamics Laboratory, NOAA - National Oceanic and Atmospheric Administration) Finite-­Volume Cubed-Sphere Dynamical Core (FV3) is a scalable and flexible dynamical core capable of both hydrostatic and non-hydrostatic atmospheric simulations. FV3 has been chosen as the dynamical core for the Next Generation Global Prediction System project (NGGPS), designed to upgrade the current NCEP (National Centers for Environmental Prediction, NOAA) operational Global Forecast System (GFS) to run as a unified, fully-coupled system in NOAA's Environmental Modeling System infrastructure. FV3 dynamic core has a long history of serving as the main engine for global atmospheric models at various Government and Academic Research Laboratories including NOAA: GFDL Climate Modeling Suite (AM4 (Atmosphere Modeling 4), CM4 (Climate Model 4) ESM4 (Earth System Model 4), Hiram (HIgh Resolution Atmospheric Model)); NASA: GMAO (Global Modeling and Assimilation Office) Goddard Earth Observing System Model (GEOS); and NCAR (National Center for Atmospheric Research) Community Earth System Model (CESM). The three primary stakeholders in FV3 (GFDL, GMAO, EMC (Environmental Modeling Center - NOAA)) have embarked on synergistic and collaborative strategies with focus on the advancement of non-hydrostratic dynamic core, physics, chemistry, and data assimilation efforts, leveraging collective strengths of each of the partnering agencies. This talk presents the ongoing plans for model advancements at GFDL, GMAO and EMC, with special emphasis on contributions towards developing community based unified global modeling system for operational and research applications at respective organizations. Future plans include extending the collaborations to the development of earth system components including GMAO's aerosol chemistry models (GOCART/MAM (Goddard Chemistry Aerosol Radiation and Transport / Modal Aerosol Model)), Land Information System (LIS), and advanced data assimilation techniques; and GFDL's Modular Ocean Model (MOM) and Sea Ice Simulator (SIS) for transition to operations at NCEP.

Putman, William↗
NASA NTRS2021

Precipitable water vapor over oceans from the Maritime Aerosol Network: Evaluation of global models and satellite products under clear sky conditions

We present results from an evaluation of precipitable water vapor (W) over remote oceanic areas as derived from global reanalysis models and from satellites against observations from the Maritime Aerosol Network (MAN) for cloudless skies during the period of 2004–2017. They cover polar, mid latitude and tropical oceanic regions and represent a first effort to use MAN observations for such evaluation. The global reanalysis model products evaluated in this study are from the Modern-Era Retrospective analysis for Research and Applications Version 2 (MERRA-2), the European Centre for Medium-Range Weather Forecasts (ECMWF) Interim Reanalysis (ERA I), and the Climate Forecast System Reanalysis (CFSR) model. The satellite products evaluated are from the Moderate Resolution Imaging Spectroradiometer (MODIS), the Polarization and Directionality of the Earth's Reflectances (POLDER), the Global Ozone Monitoring Experiment (GOME-2), the Scanning Imaging Absorption Spectrometer for Atmospheric Chartography (SCIAMACHY), and the Atmospheric Infra-red Sounder (AIRS). Satellite retrievals of W are based on the attenuation of solar reflected light by water vapor absorption bands, except those from AIRS that rely on brightness temperature measurements. A very good agreement is observed between the model estimates and MAN, with mean differences of ~5% and standard deviations of ~15%. These results are within the uncertainties associated with the models and the measurements, indicating the skill of the reanalysis models to estimate W over oceans under clear sky conditions. Mean differences of W between the satellite and MAN products are ~11, 6.7, 12, −7, and 3% for MODIS, POLDER, GOME-2, SCIAMACHY and AIRS respectively, while their standard deviations are 31, 29, 28, 20 and 17%. These differences reveal the need to address inconsistencies among different satellite sensors and ground-based measurements to reduce the uncertainties associated with the retrievals.

Daniel Pérez-Ramírez↗
NASA NTRS1994

Technical report series on global modeling and data assimilation. Volume 2: Direct solution of the implicit formulation of fourth order horizontal diffusion for gridpoint models on the sphere

High order horizontal diffusion of the form K Delta(exp 2m) is widely used in spectral models as a means of preventing energy accumulation at the shortest resolved scales. In the spectral context, an implicit formation of such diffusion is trivial to implement. The present note describes an efficient method of implementing implicit high order diffusion in global finite difference models. The method expresses the high order diffusion equation as a sequence of equations involving Delta(exp 2). The solution is obtained by combining fast Fourier transforms in longitude with a finite difference solver for the second order ordinary differential equation in latitude. The implicit diffusion routine is suitable for use in any finite difference global model that uses a regular latitude/longitude grid. The absence of a restriction on the timestep makes it particularly suitable for use in semi-Lagrangian models. The scale selectivity of the high order diffusion gives it an advantage over the uncentering method that has been used to control computational noise in two-time-level semi-Lagrangian models.

Max J. Suarez↗