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 271 records · Page 15

Exploring the Utility of Machine Learning-Based Passive Microwave Brightness Temperature Data Assimilation over Terrestrial Snow in High Mountain Asia

This study explores the use of a support vector machine (SVM) as the observation operator within a passive microwave brightness temperature data assimilation framework (herein SVM-DA) to enhance the characterization of snow water equivalent (SWE) over High Mountain Asia (HMA). A series of synthetic twin experiments were conducted with the NASA Land Information System (LIS) at a number of locations across HMA. Overall, the SVM-DA framework is effective at improving SWE estimates (~70% reduction in RMSE relative to the Open Loop) for SWE depths less than 200 mm during dry snowpack conditions. The SVM-DA framework also improves SWE estimates in deep, wet snow (~45% reduction in RMSE) when snow liquid water is well estimated by the land surface model, but can lead to model degradation when snow liquid water estimates diverge from values used during SVM training. In particular, two key challenges of using the SVM-DA framework were observed over deep, wet snowpacks. First, variations in snow liquid water content dominate the brightness temperature spectral difference (TB) signal associated with emission from a wet snowpack, which can lead to abrupt changes in SWE during the analysis update. Second, the ensemble of SVM-based predictions can collapse (i.e., yield a near-zero standard deviation across the ensemble) when prior estimates of snow are outside the range of snow inputs used during the SVM training procedure. Such a scenario can lead to the presence of spurious error correlations between SWE and TB, and as a consequence, can result in degraded SWE estimates from the analysis update. These degraded analysis updates can be largely mitigated by applying rule-based approaches. For example, restricting the SWE update when the standard deviation of the predicted TB is greater than 0.05 K helps prevent the occurrence of filter divergence. Similarly, adding a thin layer (i.e., 5 mm) of SWE when the synthetic TB is larger than 5 K can improve SVM-DA performance in the presence of a precipitation dry bias. The study demonstrates that a carefully constructed SVM-DA framework cognizant of the inherent limitations of passive microwave-based SWE estimation holds promise for snow mass data assimilation.

Kwon, Yonghwan↗

Satellite Altimetry, Ocean Circulation, and Data Assimilation

Ocean circulation is a critical factor in determining the Earth's climate. Satellite altimetry has been proven a powerful technique for measuring the height of the sea surface for the study of global ocean circulation dynamics. A major objective of my research is to investigate the utility of altimeter data for ocean circulation studies. The 6 years' data record of TOPEX/POSEIDON have been analyzed to study the spatial and temporal characteristics of large-scale ocean variability. A major result obtained in 1998 is the discovery of large-scale oscillations in sea level with a period of 25 days in the Argentine Basin of the South Atlantic Ocean (see diagram). They exhibit a dipole pattern with counterclockwise rotational propagation around the Zapiola Rise (centered at 45S and 317E), a small seamount in the abyssal plain of the basin. The peak-to-trough amplitude is about 10 cm over a distance of 500-1000 km. The amplitude of these oscillations has large seasonal-to-interannual variations. The period and rotational characteristics of these oscillations are remarkably similar to the observations made by two current meters deployed near the ocean bottom in the region. What TOPEX/POSEIDON has detected apparently are manifestations of the movement of the entire water column (barotropic motion). The resultant transport variation is estimated to be about 50 x 10(exp 6) cubic M/S, which is about 50% of the total water transport in the region. Preliminary calculations suggest that these oscillations are topographically trapped waves. A numerical model of the South Atlantic is used to investigate the nature of and causes for these waves. A very important property of sea surface height is that it is directly related to the surface geostrophic velocity, which is related to deep ocean circulation through the density field. Therefore altimetry observations are not only useful for determining the surface circulation but also for revealing information about the deep ocean. Another objective of my research is to investigate the extent of information one can infer about the deep ocean from altimetric observations. The approach is data assimilation by ocean models. The objective is to develop a modeling/data assimilation system that will produce estimates of the 3-dimensional state of the ocean for the entire duration of the TOPEX/POSEIDON Mission.

Fu, Lee-Lueng↗

Contributions of irrigation modeling, soil moisture and snow data assimilation to high-resolution water budget estimates over the Po basin: progress towards digital replicas

High-resolution water budget estimates benefit from modeling of human water management and satellite data assimilation (DA) in river basins with a large human footprint. Utilizing the Noah-MP land surface model with dynamic vegetation growth and river routing, in combination with an irrigation module, Sentinel-1 backscatter and snow depth retrievals, we produce a set of 0.7-km2 water budget estimates of the Po river basin (Italy) for 2015–2023. The results demonstrate that irrigation modeling improves the seasonal soil moisture variation and summer streamflow at all gauges in the valley after withdrawal of irrigation water from the streamflow in postprocessing (12% error reduction relative to observed low summer streamflow), even if the basin-wide irrigation amount is underestimated. Sentinel-1 backscatter DA for soil moisture updating strongly interacts with irrigation modeling: when both are activated, the soil moisture updates are limited, and the simulated irrigation amounts are reduced. Backscatter DA systematically reduces soil moisture in the spring, which improves downstream spring streamflow. Assimilating Sentinel-1 snow depth retrievals over the surrounding Alps and Apennines further improves spring streamflow in a complementary way (2% error reduction relative to observed high spring streamflow). Despite the seasonal improvements, irrigation modeling and Sentinel-1 backscatter DA cannot significantly improve short-term or interannual variations in soil moisture, irrigation modeling causes a systematically prolonged high vegetation productivity, and snow depth DA only impacts the deep snowpacks. This study helps advancing the design of digital water budget replicas for river basins.

Gabrielle J M De Lannoy↗

Simulation experiments for testing the assimilation of geostationary satellite temperature retrievals into a numerical prediction model

A simulation experiment is presented which tests a procedure that assimilates temperature soundings from geostationary satellites using a method developed by Gal-Chen (1983). The simulation is used to decide whether the increased temporal resolution of geostationary satellite data can be used to improve numerical simulation of a baroclinically unstable wave, and to assess the impact of gaps in the satellite sounding fields on the model simulations. Background information on the simulation experiment is presented, including discussion of the variational approach utilized, the numerical model, updating procedures, and control run. The simulation results indicate that, for the case of a baroclinically unstable wave, considerable improvements in the short-range forecast may be realized if geostationary satellite data are inserted with a frequency near 1 hr during a 6 hr analysis cycle compared with a single insertion.

Gal-Chen, T.↗

Characterizing the Vertical Profile of Aerosol Particle Extinction and Linear Depolarization over Southeast Asia and the Maritime Continent: The 2007-2009 View from CALIOP

Vertical profiles of 0.532 μm aerosol particle extinction coefficient and linear volume depolarization ratio are described for Southeast Asia and the Maritime Continent. Quality-screened and cloud-cleared Version 3.01 Level 2 NASA Cloud Aerosol Lidar with Orthogonal Polarization (CALIOP) 5-km Aerosol Profile datasets are analyzed from 2007 to 2009. Numerical simulations from the U.S. Naval Aerosol Analysis and Predictive System (NAAPS), featuring two-dimensional variational assimilation of NASA Moderate Resolution Imaging Spectroradiometer and Multi-angle Imaging Spectro- Radiometer quality-assured datasets, combined with regional ground-based lidar measurements, are considered for assessing CALIOP retrieval performance, identifying bias, and evaluating regional representativeness. CALIOP retrievals of aerosol particle extinction coefficient and aerosol optical depth (AOD) are high over land and low over open waters relative to NAAPS (0.412/0.312 over land for all data points inclusive, 0.310/0.235 when the per bin average is used and each is treated as single data points; 0.102/0.151 and 0.086/0.124, respectively, over ocean). Regional means, however, are very similar (0.180/0.193 for all data points and 0.155/0.159 when averaged per normalized bin), as the two factors offset one another. The land/ocean offset is investigated, and discrepancies attributed to interpretation of particle composition and a-priori assignment of the extinction-to-backscatter ratio ("lidar ratio") necessary for retrieving the extinction coefficient from CALIOP signals. Over land, NAAPS indicates more dust present than CALIOP algorithms are identifying, indicating a likely assignment of a higher lidar ratio representative of more absorptive particles. NAAPS resolvesmore smoke overwater than identified with CALIOP, indicating likely usage of a lidar ratio characteristic of less absorptive particles to be applied that biases low AOD there. Over open waters except within the Bay of Bengal, aerosol particle scattering is largely capped below 1.5 km MSL, though ground-based lidar measurements at Singapore differ slightly from this finding. Significant aerosol particle presence over land is similarly capped near 3.0 km MSL over most regions. Particle presence at low levels regionally, except over India, is dominated by relatively non-depolarizing particles. Industrial haze, sea salt droplets and fresh smoke are thus most likely present.

linear volume depolarization ratio↗

Scientific goals of the Cooperative Multiscale Experiment (CME)

Mesoscale Convective Systems (MCS) form the focus of CME. Recent developments in global climate models, the urgent need to improve the representation of the physics of convection, radiation, the boundary layer, and orography, and the surge of interest in coupling hydrologic, chemistry, and atmospheric models of various scales, have emphasized the need for a broad interdisciplinary and multi-scale approach to understanding and predicting MCS's and their interactions with processes at other scales. The role of mesoscale systems in the large-scale atmospheric circulation, the representation of organized convection and other mesoscale flux sources in terms of bulk properties, and the mutually consistent treatment of water vapor, clouds, radiation, and precipitation, are all key scientific issues concerning which CME will seek to increase understanding. The manner in which convective, mesoscale, and larger scale processes interact to produce and organize MCS's, the moisture cycling properties of MCS's, and the use of coupled cloud/mesoscale models to better understand these processes, are also major objectives of CME. Particular emphasis will be placed on the multi-scale role of MCS's in the hydrological cycle and in the production and transport of chemical trace constituents. The scientific goals of the CME consist of the following: understand how the large and small scales of motion influence the location, structure, intensity, and life cycles of MCS's; understand processes and conditions that determine the relative roles of balanced (slow manifold) and unbalanced (fast manifold) circulations in the dynamics of MCS's throughout their life cycles; assess the predictability of MCS's and improve the quantitative forecasting of precipitation and severe weather events; quantify the upscale feedback of MCS's to the large-scale environment and determine interrelationships between MCS occurrence and variations in the large-scale flow and surface forcing; provide a data base for initialization and verification of coupled regional, mesoscale/hydrologic, mesoscale/chemistry, and prototype mesoscale/cloud-resolving models for prediction of severe weather, ceilings, and visibility; provide a data base for initialization and validation of cloud-resolving models, and for assisting in the fabrication, calibration, and testing of cloud and MCS parameterization schemes; and provide a data base for validation of four dimensional data assimilation schemes and algorithms for retrieving cloud and state parameters from remote sensing instrumentation.

Cotton, William↗

The Yearly Variation in Fall-Winter Arctic Winter Vortex Descent

Using the change in HALOE methane profiles from early September to late March, we have estimated the minimum amount of diabatic descent within the polar which takes place during Arctic winter. The year to year variations are a result in the year to year variations in stratospheric wave activity which (1) modify the temperature of the vortex and thus the cooling rate; (2) reduce the apparent descent by mixing high amounts of methane into the vortex. The peak descent amounts from HALOE methane vary from l0km -14km near the arrival altitude of 25 km. Using a diabatic trajectory calculation, we compare forward and backward trajectories over the course of the winter using UKMO assimilated stratospheric data. The forward calculation agrees fairly well with the observed descent. The backward calculation appears to be unable to produce the observed amount of descent, but this is only an apparent effect due to the density decrease in parcels with altitude. Finally we show the results for unmixed descent experiments - where the parcels are fixed in latitude and longitude and allowed to descend based on the local cooling rate. Unmixed descent is found to always exceed mixed descent, because when normal parcel motion is included, the path average cooling is always less than the cooling at a fixed polar point.

Schoeberl, Mark R.↗

Geomagnetic Variability and Predictability

The Earth's magnetic field varies on time scales from seconds to millions of years. These variations can be measured directly by ground observatories and satellites, and indirectly through laboratories. They are from different electromagnetic processes in magnetosphere and ionosphere (external), ocean and crust (surface), and deep in the Earth's fluid outer core (internal). Of the measured magnetic signals, 95% is from the core, and is called the core field, or internal field. Variation of the core field, often called the geomagnetic secular variation (SV), is the manifestation of the magnetohydrodynamic processes in the core. Therefore, SV provides rich information on the core dynamical state that is then critical for understanding the geophysical mechanisms of the SV, the core-mantle interactions, and contribution of mass transport inside the core on Earth's gravity variation. Traditionally, the core state is probed via two independent methodologies: core flow inversion from the observed SV, and pure numerical dynamo simulation. The former is driven by observations, but does not include necessary interactions among physical variables; the latter is dynamically consistent, but is not constrained by any observation. Therefore, both are very limited in providing appropriate estimates of the core state. To avoid these limitations, a new approach, geomagnetic data assimilation, appears recently. In this approach, simulation results (model forecasts) are constantly corrected with observations, the corrected solutions (analysis), are then used as the initial conditions to make more accurate forecasts of future. Recent studies have demonstrated successfully that assimilation solutions are different from those of free running models. These new solutions are better estimates of the core state, and have been used to provide accurate prediction of SV for the period from 2010 to 2015. This prediction is part of the IGRF field model for international community applications.

Kuang, Weijia↗

Ionospheric convection response to slow, strong variations in a Northward interplanetary magnetic field: A case study for January 14, 1988

We analyze ionospheric convection patterns over the polar regions during the passage of an interplanetary magnetic cloud on January 14, 1988, when the interplanetary magnetic field (IMF) rotated slowly in direction and had a large amplitude. Using the assimilative mapping of ionospheric electrodynamics (AMIE) procedure, we combine simultaneous observations of ionspheric drifts and magnetic perturbations from many different instruments into consistent patterns of high-latitude electrodynamics, focusing on the period of northward IMF. By combining satellite data with ground-based observations, we have generated one of the most comprehensive data sets yet assembled and used it to produce convection maps for both hemispheres. We present evidence that a lobe convection cell was embedded within normal merging convection during a period when the IMF B(sub y) and B(sub z) components were large and positive. As the IMF became predominantly northward, a strong reversed convection pattern (afternoon-to-morning potential drop of around 100 kV) appeared in the southern (summer) polar cap, while convection in the northern (winter) hemisphere became weak and disordered with a dawn-to-dust potential drop of the order of 30 kV. These patterns persisted for about 3 hours, until the IMF rotated significantly toward the west. We interpret this behavior in terms of a recently proposed merging model for northward IMF under solstice conditions, for which lobe field lines from the hemisphere tilted toward the Sun (summer hemisphere) drape over the dayside magnetosphere, producing reverse convection in the summer hemisphere and impeding direct contact between the solar wind and field lines connected to the winter polar cap. The positive IMF B(sub x) component present at this time could have contributed to the observed hemispheric asymmetry. Reverse convection in the summer hemisphere broke down rapidly after the ratio absolute value of B(sub y)/B(sub z) exceeded unity, while convection in the winter hemisphere strengthened. A dominant dawn-to-dusk potential drop was established in both hemispheres when the magnitude of B(sub y) exceeded that of B(sub z) with potential drops of the order of 100 kV, even while B(sub z) remained northward. The latter transition to southward B(sub z) produced a gradual intensification of the convection, but a greater qualitative change occurred at the transition through absolute value of B(sub y)/B(sub z) = 1 than at at the transition through B(sub z) = 0. The various convection patterns we derive under northward IMF conditions illustrate all possibilities previously discussed in the literature: nearly single-cell and multicell, distorted and symmetric, ordered and unordered, and sunward and antisunward.

Knipp, D. J.↗

Evaluation of TAO Observation System on ENSO Predictions from the GMAO S2S Forecast System

Since the early 1990’s, the Tropical Atmosphere Ocean (TAO) mooring array has been instrumental for observing Kelvin and Rossby wave propagation which is key for El Nino/Southern Oscillation (ENSO) prediction. However, recent funding and programmatical issues have highlighted the need to rigorously assess the impact of the TAO observing system on ENSO predictions. Therefore, we evaluate the TAO using data assimilation observation denial experiments (also known as Observing System Evaluation or OSE experiments). This presentation will evaluate the TAO observing impact on both reanalyses and prediction of the big 2015 El Nino. We have completed reanalyses for July 2014-Dec 2015 for both the CONTROL experiment (that assimilates all available data) and the NOTAO experiment (that is identical to the CONTROL but withholds all TAO observations) using the best available NASA GEOS-S2S V3 seasonal prediction system. Validation of these reanalyses shows that TAO assimilation generally improves comparisons of temperature and salinity versus gridded in situ observations. Temperature is universally improved above the thermocline near the equator, but the biggest improvement is found in the eastern Pacific, just below the thermocline, where the variation of the thermocline defines ENSO events. In addition, a surprising result is that even with relatively few observations, salinity is improved throughout the equatorial region except near 120oW near the surface. ENSO forecasts were performed that were initialized from these CONTROL and NOTAO reanalyses. For the 9-month forecasts which were initialized in January, July, and October 2015, the NINO3.4 SST shows that the CONTROL forecasts are warmer and closer to observations than the NOTAO forecasts. We will show that upwelling and shoaling of the mixed layer amplifies the ENSO signal due to TAO assimilation. Prior to (after) April 2015, this upwelling is caused by relatively stronger Rossby (Kelvin) waves in the CONTROL than in the NOTAO experiments.

E. Hackert↗

Use of three-cornered hat error estimates in MERRA-2 to guide an improved reanalysis-Part 1

The three-cornered hat (3CH) method estimates the uncertainties of three different co-located model or observational data sets (Anthes and Rieckh, 2018; Sjoberg et al., 2021). Rieckh et al. (2021) used the 3CH method to compare the random error statistics of different global forecast and reanalysis models, as well as radio occultation (RO) and radiosonde observations. That study showed that the MERRA-2 reanalysis, while having smaller errors in the stratosphere than its predecessor MERRA, had larger errors in the troposphere than many of the other data sets analyzed. The MERRA-2 errors were particularly large in the tropics. In a collaborative effort between UCAR’s COSMIC (Constellation Observing System for Meteorology, Ionosphere and Meteorology) program and NASA’s Global Modeling and Assimilation Office (GMAO), we carried out further 3CH error diagnostics to help isolate the causes of these larger errors and help guide the development of an improved reanalysis. This presentation summarizes random error statistics associated with MERRA-2, ECMWF’s ERA5 reanalysis, and COSMIC-2 (C2) RO observations. We compute 3CH error variance estimates of refractivity, as well as temperature and specific humidity using UCAR’s COSMIC Data Analysis and Archive Center (CDAAC) improved 1D-variational (1D-Var) retrieval (wetPf2) over 15 latitude bands from 45S to 45N. The 1D-Var retrievals of specific humidity and temperature for C2 use NCEP’s Global Forecast System (GFS) as the background. Anthes et al. (2021) showed that it gives accurate estimates of temperature and specific humidity in the tropics and subtropics, even in the challenging environment of intense Hurricane Dorian (2019). This presentation confirms the previous results that MERRA-2 has significantly larger errors in the tropics and subtropics than either C2 or ERA5. Its errors are larger between 30S and 30N compared to 30-45 N-S latitudes, and are also larger over land compared to oceans. Most of the MERRA-2 refractivity errors come from specific humidity, except over land below 3 km where temperature errors are large. These results suggest that moist convection and atmospheric boundary layer physics in MERRA-2 may be responsible for a significant part of the higher uncertainties. These results are being used to guide GMAO in developing an improved next-generation reanalysis, as shown in a companion presentation submitted to this conference (El Akkraoui et al., 2021), which extends this study and describes improvements to MERRA-2 leading to the next GMAO reanalysis.

Jeremiah Sjoberg↗

Effect of Cross-Correlation on Geomagnetic Forecast Accuracies

Surface geomagnetic observation can determine up to degree L = 14 time-varying spherical harmonic coefficients of the poloidal magnetic field. Assimilation of these coefficients to numerical dynamo simulation could help us understand better the dynamical processes in the Earth's outer core, and to provide more accurate forecast of geomagnetic secular variations (SV). In our previous assimilation studies, only the poloidal magnetic field in the core is corrected by the observations in the analysis. Unobservable core state variables (the toroidal magnetic field and the core velocity field) are corrected via the dynamical equations of the geodynamo. Our assimilation experiments show that the assimilated core state converges near the CMB, implying that the dynamo state is strongly constrained by surface geomagnetic observations, and is pulled closer to the truth by the data. We are now carrying out an ensemble of assimilation runs with 1000 years of geomagnetic and archeo/paleo magnetic record. In these runs the cross correlation between the toroidal and the poloidal magnetic fields is incorporated into the analysis. This correlation is derived from the physical boundary conditions of the toroidal field at the core-mantle boundary (CMB). The assimilation results are then compared with those of the ensemble runs without the cross-correlation, aiming at understanding two fundamental issues: the effect of the crosscorrelation on (1) the convergence of the core state, and (2) the SV prediction accuracies. The constrained dynamo solutions will provide valuable insights on interpreting the observed SV, e.g. the near-equator magnetic flux patches, the core-mantle interactions, and possibly other geodynamic observables.

Kuang, Weijia↗

NASA GLDAS Evapotranspiration Data and Climatology

Evapotranspiration (ET) is the water lost to the atmosphere by evaporation and transpiration. ET is a shared component in the energy and water budget, therefore, a critical variable for global energy and water cycle and climate change studies. However, direct ET measurements and data acquisition are difficult and expensive, especially at the global level. Therefore, modeling is one common alternative for estimating ET. With the goal to generate optimal fields of land surface states and fluxes, the Global Land Data Assimilation System (GLDAS) has been generating quality-controlled, spatially and temporally consistent, terrestrial hydrologic data, including ET and other variables that affect evaporation and transpiration, such as temperature, precipitation, humidity, wind, soil moisture, heat flux, and solar radiation. This poster presents the long-term ET climatology (mean and monthly), derived from the 61-year GLDAS-2 monthly 1.0 deg x 1.0 deg. NOAH model Experiment-1 data, and describes the basic characteristics of spatial and seasonal variations of the climatology. The time series of GLDAS-2 precipitation and radiation, and ET are also discussed to show the improvement of GLDAS-2 forcing data and model output over those from GLDAS-1.

Rui, Hualan↗

Some Continuum Aspects of Data Assimilation

A long-sought goal in data assimilation is to build more fully the time dimension into the data assimilation process. Truly four-dimensional data assimilation requires evolving second-moment information, namely the estimation error covariance, along with the state estimate itself. The continuum evolution equations for the estimation error covariance constitute a system of partial differential equations in six space dimensions, forced by a model error covariance and by the observations themselves. The high dimensionality of this system poses a difficult computational problem. Numerous methods have been proposed for approximate, discrete solution of this system of equations. Four-dimensional variational schemes solve these equations indirectly, while schemes based on Kalman filtering solve the equations more directly, usually on a small subspace of the full six-dimensional space. In both cases, most proposed solution methods are not derived from the continuum covariance evolution equations themselves. Instead, most methods simply inherit the discretization used for evolving the state estimate, either in the form of a discrete tangent linear model or, in the case of ensemble schemes, in the form of a fully nonlinear discrete model. In this lecture we show that solutions of the continuum covariance evolution equations possess simple properties that are not readily satisfied upon discretization of the equations. we give numerical examples illustrating that "inherited" discretizations sometimes fail to produce meaningfully accurate solutions. Finally, we suggest alternative discretization methods that may yield more faithful approximate solutions of the continuum problem.

Cohn, Stephen E.↗

Investigating Meteorological Conditions Associated with Valley Fever Outbreaks in Southwestern United States

According to CDC, a sharp increase in reported Valley fever cases (Coccidioidomycosis) has been found between 2009 and 2012. Environmental conditions play an important role for Valley fever outbreaks. For example, dust storms can significantly increase the amount of the fungus Coccidioides in the air and associated strong winds can transport the fungus to other areas. In last year’s AGU session, we reported a preliminary study to investigate hydrometeorological conditions and their connection with dust storm activities in southwestern United States. We found wind is a major contributing factor for the seasonal variation of dust storm activities. Interannual variation of the regional hydrometeorological conditions are closely linked to the large-scale environment such as the Pacific Decadal Oscillation (PDO). Strong winds are linked with a number of weather events such cold front passages, thunderstorms that produce downbursts and strong winds, the Santa Ana winds, etc. In this presentation, we will report the results of our latest investigation on meteorological conditions associated with Valley fever outbreaks in southwestern United States, using NLDAS (North American Land Data Assimilation System), GPM (Global Precipitation Measurement), and MERRA-2 (Modern Era Retrospective-analysis for Research and Applications, Version 2) hourly datasets, from the NASA Goddard Earth Sciences (GES) Data and Information Services Center (DISC).

Liu, Zhong↗

Observations of the Thermal Structure and Dynamics of the Martian Atmosphere

Information on the thermal structure of the Martian atmosphere has been obtained using several different techniques, including in situ sounding from decent probes, radio occultations and infrared thermal emission measurements from orbiters, and earth based full-disk microwave measurements. Knowledge of the 3- dimensional atmospheric temperature field and its temporal variation provides an important observational approach to the study of the dynamics of the atmosphere. Comparison of the observed thermal structure with results from a Mars General Circulation Model (MGCM) can be used to refine the model and provide greater insight into the physical processes controlling the dynamics of the atmosphere. Direct assimilation of the temperature observations into an MGCM can also be used to address similar goals. In addition, relatively simple diagnostic models can be applied to the data to directly obtain information on winds and other meteorological properties of the atmosphere. In this overview, we concentrate primarily on the structure retrieved from thermal emission measurements. Additional information is contained in the original extended abstract.

Conrath, B. J.↗

Interannual Variation in Phytoplankton Class-specific Primary Production at a Global Scale

Phytoplankton is responsible for over half of the net primary production on earth. The knowledge on the contribution of various phytoplankton groups to the total primary production is still poorly understood. Data from satellite observations suggest that for upwelling regions, photosynthetic rates by microplankton is higher than that of nanoplankton but that when the spatial extent is considered, the production by nanoplankton is comparable or even larger than microplankton. Here, we used the NASA Ocean Biogeochemical Model (NOBM) combined with remote sensing data via assimilation to evaluate the contribution of 4 phytoplankton groups to the total primary production. Globally, diatoms were the group that contributed the most to the total phytoplankton production (approx. 50%) followed by coccolithophores and chlorophytes. Primary production by diatoms was highest in high latitude (>45 deg) and in major upwelling systems (Equatorial Pacific and Benguela system). We assessed the effects of climate variability on the class-specific primary production using global (i.e. Multivariate El Nino Index, MEI) and 'regional' climate indices (e.g. Southern Annular Mode (SAM), Pacific Decadal Oscillation (PDO) and North Atlantic Oscillation (NAO)). Most interannual variability occurred in the Equatorial Pacific and was associated with climate variability. These results provide a modeling and data assimilation perspective to phytoplankton partitioning of primary production and contribute to our understanding of the dynamics of the carbon cycle in the oceans at a global scale.

OBM↗

The assimilation of satellite soundings, winds and satellite products in a mesoscale analysis/forecast system

Investigations in FY-85 were centered on three case study days in 1982. Two of these, March 6 and April 24, were Atmospheric Variability Experiment/Verical Atmospheric Sounder (AVE/VAS) days for which high spatial and temporal resolution RAOB and Vertical Atmospheric Sounder (VAS) data sets were available. The third investigation day, April 26, was a day of interesting severe weather. In the last part of FY-84 and early FY-85 we were able to demonstrate most importantly the complimentary nature of satellite soundings and winds in a forecast/analysis system. In our variational analysis scheme, cloud drift and water vapor winds enter into the height field as gradient information. The cloud drift winds especially, have the character of supplying information in cloudy areas where satellite soundings are not possible. In the April 26 experiments, analyses and forecasts using the combination satellite winds and soundings were superior to those using only soundings. Good consistency was shown between independent satellite forecasts from different initialization times run to the same verification time. A significant accomplishment in FY-85 was expanding experiments on April 26 to include quasi-continuous initialization inserting satellite soundings and winds from several different times into an analysis/forecast. Contrary to the first set of experiments on April 26, here forecast initialization fields were not independent, but contained satellite information from two data times.

Diak, G. R.↗