Search NASA⌕ Search

SEARCH · Search NASA

Results for “reanalyses”

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 145 records · Page 8

Impact of Land Surface Initialization Approach on Subseasonal Forecast Skill: a Regional Analysis in the Southern Hemisphere

The authors use a sophisticated coupled land-atmosphere modeling system for a Southern Hemisphere subdomain centered over southeastern Australia to evaluate differences in simulation skill from two different land surface initialization approaches. The first approach uses equilibrated land surface states obtained from offline simulations of the land surface model, and the second uses land surface states obtained from reanalyses. The authors find that land surface initialization using prior offline simulations contribute to relative gains in subseasonal forecast skill. In particular, relative gains in forecast skill for temperature of 10%-20% within the first 30 days of the forecast can be attributed to the land surface initialization method using offline states. For precipitation there is no distinct preference for the land surface initialization method, with limited gains in forecast skill irrespective of the lead time. The authors evaluated the asymmetry between maximum and minimum temperatures and found that maximum temperatures had the largest gains in relative forecast skill, exceeding 20% in some regions. These results were statistically significant at the 98% confidence level at up to 60 days into the forecast period. For minimum temperature, using reanalyses to initialize the land surface contributed to relative gains in forecast skill, reaching 40% in parts of the domain that were statistically significant at the 98% confidence level. The contrasting impact of the land surface initialization method between maximum and minimum temperature was associated with different soil moisture coupling mechanisms. Therefore, land surface initialization from prior offline simulations does improve predictability for temperature, particularly maximum temperature, but with less obvious improvements for precipitation and minimum temperature over southeastern Australia.

Precipitation↗

Sensitivity of Simulated Global Ocean Carbon Flux Estimates to Forcing by Reanalysis Products

Reanalysis products from MERRA, NCEP2, NCEP1, and ECMWF were used to force an established ocean biogeochemical model to estimate air-sea carbon fluxes (FCO2) and partial pressure of carbon dioxide (pCO2) in the global oceans. Global air-sea carbon fluxes and pCO2 were relatively insensitive to the choice of forcing reanalysis. All global FCO2 estimates from the model forced by the four different reanalyses were within 20% of in situ estimates (MERRA and NCEP1 were within 7%), and all models exhibited statistically significant positive correlations with in situ estimates across the 12 major oceanographic basins. Global pCO2 estimates were within 1% of in situ estimates with ECMWF being the outlier at 0.6%. Basin correlations were similar to FCO2. There were, however, substantial departures among basin estimates from the different reanalysis forcings. The high latitudes and tropics had the largest ranges in estimated fluxes among the reanalyses. Regional pCO2 differences among the reanalysis forcings were muted relative to the FCO2 results. No individual reanalysis was uniformly better or worse in the major oceanographic basins. The results provide information on the characterization of uncertainty in ocean carbon models due to choice of reanalysis forcing.

MERRA↗

Using Reanalysis in Crop Monitoring and Forecasting Systems

Weather observations are essential for crop monitoring and forecasting but they are not always available and in some cases they have limited spatial representativeness. Thus, reanalyses represent an alternative source of information to be explored. In this study, we assess the feasibility of reanalysis-based crop monitoring and forecasting by using the system developed and maintained by the European Commission- Joint Research Centre, its gridded daily meteorological observations, the biased-corrected reanalysis AgMERRA and the ERA-Interim reanalysis. We focus on Europe and on two crops, wheat and maize, in the period 1980-2010 under potential and water-imited conditions. In terms of inter-annual yield correlation at the country scale, the reanalysis-driven systems show a very good performance for both wheat and maize (with correlation values higher than 0.6 in almost all EU28 countries) when compared to the observations-driven system. However, significant yield biases affect both crops. All simulations show similar correlations with respect to the FAO reported yield time series. These findings support the integration of reanalyses in current crop monitoring and forecasting systems and point to the emerging opportunities linked to the coming availability of higher-resolution reanalysis updated at near real time.

Reanalysis↗

Impact of Satellite Sea Surface Salinity Observations on ENSO Predictions from the GMAO S2S Forecast System

El Nino/Southern Oscillation (ENSO) has far reaching global climatic impacts and so extending useful ENSO forecasts would be of great benefit for society. However, one key variable that has yet to be fully exploited within coupled forecast systems is accurate estimation of near-surface ocean density. Satellite Sea Surface Salinity (SSS), combined with temperature, help to identify ocean density changes and associated mixing near the ocean surface. We assess the impact of satellite SSS observations for improving near-surface dynamics within ocean analyses and how these impact dynamical ENSO forecasts using the NASA GMAO (Global Modeling and Assimilation Office) Sub-seasonal to Seasonal (S2S_v2.1) coupled forecast system (Molod et al. 2018 - i.e. NASA's contribution to the NMME (National Multi-Model Ensemble) project). For all initialization experiments, all available along-track absolute dynamic topography and in situ observations are assimilated using the LETKF ( Local Ensemble Transform Kalman Filter) scheme similar to Penny et al., 2013. A separate reanalysis additionally assimilates Aquarius V5 (September 2011 to June 2015) and SMAP (Soil Moisture Active Passive) V4 (March 2015 to present) along-track data.We highlight the impact of satellite SSS on ocean reanalyses by comparing validation statistics of experiments that assimilate SSS versus our current prediction system that withholds SSS. We find that near-surface validation versus observed statistics for salinity are slightly degraded when assimilating SSS. This is an expected result due to known biases between SSS (measured by the satellite at approximately 1 centimeter) and in situ measurements (typically measured by Argo floats at 3 meters). On the other hand, a very encouraging result is that both temperature, absolute dynamic topography, and mixed layer statistics are improved with SSS assimilation. Previous work has shown that correcting near-surface density structure via gridded SSS assimilation can improve coupled forecasts. Here we present results of coupled forecasts that are initialized from the GMAO S2S reanalyses that assimilates/withholds along-track (L2) SSS. In particular, we contrast forecasts of the overestimated 2014 El Nino, the big 2015 El Nino, and the minor 2016 La Nina. For each of these ENSO scenarios, assimilation of satellite SSS improves the forecast validation. Improved SSS and density upgrades the mixed layer depth leading to more accurate coupled air/sea interaction.

Hackert, E.↗

Consequences of Different Air-Sea Feedbacks on Ocean Using MITgcm and MERRA-2 Forcing: Implications for Coupled Data Assimilation Systems

Ocean surface flux estimates from atmospheric and oceanic reanalyses contain errors that compensate for inaccuracies in the respective atmosphere and ocean models used to generate these reanalyses. A conundrum for climate studies is the discrepancy between surface fluxes that minimize model-data differences for an atmosphere-only model vs surface fluxes that minimize model-data differences for an ocean model. As a first step towards a consistent coupled ocean-atmosphere data-assimilation (DA) system, we compare surface net heat flux from a state-of-the-art atmospheric reanalysis, the Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2), to net heat flux from a state-of-the-art ocean state estimate, the Estimating the Circulation and Climate of the Ocean Version 4 (ECCO-v4). The possible impacts of the MERRA-2 and ECCO-v4 air-sea net heat flux difference in a coupled DA system were assessed using a set of experiments designed to imitate different “flavors” of a coupled DA system in an ocean-only setup. This was done by forcing the ECCO-v4 underlying ocean model - the Massachusetts Institute of Technology general circulation model (MITgcm) - with different sets of MERRA-2 fields and utilizing different forcing methods. By doing so we were able to turn off different air-sea feedbacks which, in a coupled DA setup, are partially muted by the constraining observations. The set of experiments, therefore, represents a range of active feedbacks in different “flavors” of coupled data-assimilation systems. For the period 1992–2011, MERRA-2 net heat flux has a global mean difference of -4.9 Wm(exp -2) relative to ECCO-v4. When MERRA-2 surface fields are used to force MITgcm, imbalances in the energy and the hydrological cycles of MERRA-2, which are directly related to the fact that MERRA-2 was created without an interactive ocean, propagate to the ocean. The experiment in which MITgcm is forced with MERRA-2 fluxes (MERRA-2-flux experiment) results in a 2.5°C global mean Sea Surface Temperature (SST) cooling, a 1m reduction in global mean sea level, and other drastic changes in the large scale ocean circulation relative to those resulting when the MITgcm is forced with the optimized ECCO-v4 net heat flux (the ECCO-v4 experiment itself). When MITgcm is forced with MERRA-2 state variables (MERRA-2-state experiment), the SST is somewhat restored to the observed SST, but the errors are shifted to the water cycle, resulting in a global mean sea level increase of 2.7 m. To further explore the pros and cons of these two approaches, we introduce a new intermediate forcing method in which the ocean is forced with turbulent fluxes but has a long wave feedback. This method, unlike MERRA-2 state, preserves the MERRA-2 water and salinity cycles, and it reduces the SST error compared to the MERRA-2-flux experiment, but the SST is not as good as that in the MERRA-2-state experiment. Our results have implications for ocean-model forcing recipes and clearly reveal the undesirable consequences of limiting the feedbacks in either these types of experiments or in coupled DA.

Ehud Strobach↗

Impact of Satellite Sea Surface Salinity Observations on ENSO Predictions from the GMAO Seasonal Forecast System

El Nino/Southern Oscillation (ENSO) has far reaching global climatic impacts and so extending useful ENSO forecasts would be of great benefit for society. However, one key variable that has yet to be fully exploited within coupled forecast systems is accurate estimation of near-surface ocean density. Satellite sea surface salinity (SSS), combined with temperature, help to identify ocean density changes and associated mixing near the ocean surface. We assess the impact of satellite SSS observations for improving near-surface dynamics within ocean analyses and how these impact dynamical ENSO forecasts using the NASA GMAO (Global Modeling and Assimilation Office) Sub-seasonal to Seasonal (S2S_v2.1) coupled forecast system (Molod et al. 2018 - i.e. NASA's contribution to the NMME (North American Multi-Model Ensemble) project). For all initialization experiments, all available along-track absolute dynamic topography and in situ observations are assimilated using the LETKF (Local Ensemble Transform Kalman Filter) scheme similar to Penny et al., 2013. A separate reanalysis additionally assimilates Aquarius V5 (September 2011 to June 2015) and SMAP (Soil Moisture Active Passive satellite) V4.1 (March 2015 to present) along-track data.We highlight the impact of satellite SSS on ocean reanalyses by comparing validation statistics of experiments that assimilate SSS versus our current prediction system that withholds SSS. We find that near-surface validation versus observed statistics for salinity are slightly degraded when assimilating SSS. This is an expected result due to known biases between SSS (measured by satellite at approximately 1-centimeter depth) and in situ measurements (typically measured by Argo floats at 3-meters depth). On the other hand, a very encouraging result is that both temperature, absolute dynamic topography, and mixed layer statistics are improved with SSS assimilation. Previous work has shown that correcting near-surface density structure via gridded SSS assimilation can improve coupled forecasts. Here we present results of coupled forecasts that are initialized from GMAO S2S spring reanalyses that assimilate/withhold along-track (L2) SSS. In particular, we contrast forecasts of the big 2015 El Nino, the 2017 La Nina and the 2018 weak El Nino. For each of these ENSO scenarios, assimilation of satellite SSS improves the forecast validation. Improved SSS and density upgrade the mixed layer depth leading to more accurate coupled air/sea interaction. From March to June 2015, the availability of two overlapping satellite SSS instruments, Aquarius and SMAP, allows a unique opportunity to compare and contrast forecasts initialized with the benefit of these two satellite SSS observation types. We assess the impact of gridded satellite sea surface salinity observations on dynamical ENSO forecasts for the big 2015 El Nino.

Hackert, E.↗

The Effects of a 1998 Observing System Change on MERRA-2-Based Ozone Profile Simulations

Model simulations of ozone (O3) driven by meteorological reanalyses are useful for filling observational gaps and interpreting observed O3 variability and trends. However, the transport circulation of reanalysis products is impacted by changes to the observing system (the data assimilated into the reanalyses). We examine the impacts of these changes on simulated O3 from two models, Global Modeling Initiative (GMI) Chemistry Transport Model (GMI CTM) and Modern‐Era Retrospective Analysis for Research Applications version 2 (MERRA‐2) GMI Replay (M2 GMI Replay) simulation, using observations from global ozonesondes (>50,000 profiles) and satellites from 1980 to 2016. Both models are constrained by meteorology from the NASA MERRA‐2 reanalysis, and both use versions of NASA's GMI chemical mechanism. We focus on an observing system change affecting simulated O3 after 1998, associated with the assimilation of temperature and humidity data from new microwave profiling satellites. A large post‐1998 O3 increase, mainly confined to 15–20‐km altitude, of ~10 Dobson units (DU) in midlatitudes occurs in the GMI CTM, worsening the bias compared to observations. In contrast, an increase in M2 GMI Replay simulation O3 of ~10 DU is observed only near −60° latitude, reducing the bias compared to observations. The GMI CTM O3 high biases display a Quasi‐Biennial Oscillation (QBO)‐like periodicity that result from excessive transport from the tropical stratosphere to the midlatitude lower stratosphere during the QBO westerly phase. We quantify O3 discontinuities caused by MERRA‐2 observing system changes and demonstrate how the MERRA‐2 Global Modeling Initiative Replay simulation dampens the effects of these changes and QBO‐driven artifacts on simulated lower stratospheric and total O3. We caution against using simulations driven by a reanalysis to derive multidecadal O3 trends, especially prior to 1998.

Replay↗

Observational Needs for Improving Ocean and Coupled Reanalysis, S2S Prediction, and Decadal Prediction

Developments in ocean data assimilation (DA) and observing system technologies are intertwined. New observation types lead to new DA methods, and new DA methods such as Coupled Data Assimilation can change the value of existing observations or indicate where new observations can have greater utility for monitoring and prediction. Practitioners are encouraged to make better use of observations that are already available, for example in strongly coupled data assimilation where ocean observations can be used to improve atmospheric analyses and vice versa. Ocean reanalyses are useful for the analysis of climate,as well as initializing operational long-range prediction models. There are remaining challenges for ocean reanalyses due to biases and abrupt changes in the ocean observing system throughout its history, the presence of biases and drifts in models, and simplifying assumptions made in the DA methods. From a governance point of view, more support is needed to interface the observing community and the ocean DA community. For prediction applications, the ocean DA community must work with the ocean observing community to establish protocols for rapid communication of ocean observing data on NWP timescales. There is potential for new observations to enhance the observing system by supporting prediction on multiple timescales, ranging from the typical timescale of numerical weather prediction covering hours to weeks, out to multiple decades. It is highly encouraged that communication be fostered between thesecommunities to allow operational prediction centers the ability to provide guidance to the design of a sustained and adaptive observing network.

Ocean reanalysis↗

Utilization of Airborne Observations to Assess Model Parametrizations of Critical RH Profiles in the Arctic Ocean

Observational paucity is a reality in the Arctic Ocean. This is especially true for near-surface variables such as temperature, moisture, heat fluxes and BL clouds. As a result, modeling has become one of the major avenues for understanding current and future Arctic trends. Reanalyses are frequently used to force global ocean circulation and sea-ice models. But in northern high latitudes, model integrations and reanalyses are known to have large uncertainties in temperature and humidity profiles, and in boundary layer cloudiness. These are common sources of error in the surface radiative budget terms.An important way to diagnose these biases spatially and temporally is by using satellite remote sensing data. However, remotely-sensed observations also have large uncertainties, especially in near-surface temperature and relative humidity profiles.In situ observational studies are important in bridging our knowledge gap in regions such as the Arctic Ocean. Here, we utilize airborne and ship observations from the ARISE, ACME-IV, and ASCOS campaigns to construct critical relative humidity (RH) profiles over the Beaufort Sea. Such profiles are used as parameterization inputs in the NASA GOES global model to derive the total water condensate in a model grid-box, which determines the cloud fraction. Currently, the critical RH profiles are derived by global AIRS data, relaying mostly on mid-latitude regions, which are not necessarily relevant to the Arctic.We derive campaign-wide mean, standard deviation, and critical RH values, for grid size of 50x50 km and altitude bins between 50 to 400 m, covering both open-ocean and sea-ice covered regions. We compare profiles over open ocean and sea-ice, and look at correlations between the observed critical RH values and water condensate (by cloud number concentration) from observations versus the modeled ones. We then input our calculated values of minimal critical RH values into a set of GEOS single column model (SCM) simulations over the ARISE and ASCOS regions and compare the differences between the predicted values of cloud liquid water path (LWP), ice water path (IWP) and surface fluxes with the observed ones under the range of input parameterizations. Finally, we discuss the implications on surface radiative budget predictions in this region.

Segal-Rozenhaimer, Michal↗

A Moments View of Climatology and Variability of the Asian Summer Monsoon Anticyclone

A comprehensive investigation of the climatology of and interannual variability and trends in the Asian summer monsoon anticyclone (ASMA) is presented, based on a novel area and moments analysis. Moments include centroid location, aspect ratio, angle, and “excess kurtosis” (measuring how far the shape is from elliptical) for an equivalent ellipse with the same area as the ASMA. Key results are robust among the three modern reanalyses studied. The climatological ASMA is nearly elliptical, with its major axis aligned along its centroid latitude and a typical aspect ratio of ∼5–8. The ASMA centroid shifts northward with height, northward and westward during development, and in the opposite direction as it weakens. New evidence finding no obvious climatological bimodality in the ASMA reinforces similar suggestions from previous studies using modern reanalyses. Most trends in ASMA moments are not statistically significant. ASMA area and duration, however, increased significantly during 1979–2018; the 1958–2018 record analyzed for one reanalysis suggests that these trends may have accelerated in recent decades. ASMA centroid latitude is significantly positively (negatively) correlated with subtropical jet core latitude (altitude), and significantly negatively correlated with concurrent ENSO; these results are consistent with and extend previous work relating monsoon intensity, ENSO, and jet shifts. ASMA area is significantly positively correlated with the MEI ENSO index two months previously. These results improve our understanding of the ASMA using consistently defined diagnostics of its size, geometry, interannual variability, and trends that have not previously been analyzed.

Anticyclone↗

Variations in radiative heating of humid biomass burning aerosols in the southeast Atlantic from airborne observations and reanalysis

The atmosphere over the southeast Atlantic Ocean (SEA) experiences consistent springtime biomass burning (BB) smoke from widespread agricultural fires on the African continent. This smoke layer is initially lofted high in a continental mixed layer (~5-6km) and is then transported westward into the free troposphere where it overlies and ultimately mixes into a stratocumulus-topped oceanic boundary layer. Coincident with this smoke is an elevated humidity signal which is present from the time a given airmass is over the continental source region. ORACLES (ObseRvations of Aerosols above CLouds and their intEractionS) was a NASA Earth Venture Suborbital mission with the goal of measuring aerosol, cloud, and atmospheric properties over this region over three deployments in September 2016, August 2017, and October 2018. Measurements included in situ profile measurements of aerosol physical and optical properties, trace gas concentrations, and cloud properties, as well as column measurements of aerosol optical depth and trace gas concentrations from an airborne sun photometer. Here we present results from each of the ORACLES deployments over the southeast Atlantic Ocean. Despite the differing locations and season of each deployment, we show that there is good agreement between the airborne ORACLES dataset and large-scale reanalyses, specifically the ECMWF ERA5 reanalyses (and, to a degree, the CAMS product) as compared to other available products. This agreement allows us to examine multi-year seasonal patterns and trends beyond the three months with available aircraft data. In both the data and reanalysis, we find distinct variations between each deployment in terms of vertical smoke distribution and correlation to atmospheric specific humidity, due to changing conditions over the BB season. We trace the origin and the history of these smoky airmasses and finally, we briefly discuss the broader radiative and dynamical implications of these results for conditions of aerosols overlying stratocumulus clouds.

K. Pistone↗

Effect of Rain-Adjusted Satellite Sea Surface Salinity on ENSO Predictions from the GMAO S2S Forecast System

ENSO has a significant impact on climate variability throughout the world and so has been the key focus for improving coupled ocean-atmosphere forecasts. Assimilation of altimetry and in situ data lead to improved short-term forecasts of the coupled system. However, few studies have focused on improving the near-surface density/mixing through assimilation of satellite sea surface salinity (SSS). For expediency, most projects that do assimilate SSS do so as if these data were observed at the top model layer instead of at the surface. In rainy regions where buoyant water sits as a fresh lens at the surface, this assumption is likely invalid. Therefore, we adjust SSS so that it more accurately represents the salinity at the first model layer. The Rain Impact Model (RIM) uses a simple diffusion model to determine the near surface salinity gradient (i.e., 1 cm to 5 m). Satellite SSS data are modified using this near-surface salinity gradient, so the salinity values are now valid at the first model layer (we call this SSS@5m). We assess the impact of satellite SSS observations for near-surface dynamics within ocean reanalyses and how these impact dynamical ENSO forecasts using the NASA GMAO S2S coupled system. For all reanalysis experiments, all available along-track sea level and in situ observations are assimilated for 2011 to 2020 using the LETKF scheme (Penny et al., 2013). One reanalysis assimilates Aquarius /SMAP SSS as before. An additional reanalysis is performed assimilating the SSS@5m data. Validation statistics are compared for experiments that assimilate SSS (sub-optimally as before) versus the SSS@5m. We also compare results of coupled forecasts that are initialized from these reanalyses in spring. For all but the big La Niña in 2017, all NINO3.4 forecasts were improved by using the RIM. We will show that improved SSS estimates upgrades density and near-surface mixing leading to more accurate coupled air/sea interaction and better forecasts.

ENSO↗

PolarMERRA: A Polar-Focused Global Reanalysis Project for Scientific and Stakeholder Needs

The adequate modeling of physical processes in the Arctic and Antarctic is vital for developing prediction capabilities and for obtaining an understanding of rapidly evolving polar conditions, including their potential global impacts. These processes are often poorly represented in global models and reanalyses, owing to a legacy modeling focus on midlatitude processes, as well as a scarcity of observations needed for model development in polar regions. The polarMERRA initiative, a joint effort between NASA’s Cryospheric Sciences and Modeling and Prediction programs, seeks to improve the representation of cryospheric and polar atmospheric processes and to develop an open-source framework for a quantitively evaluation of polar-relevant variables against current and future satellite and in-situ observations, models, and reanalyses. Here, we evaluate the impacts of spatial resolution and modifications to sea ice, ice sheet, and atmospheric parameterizations on the representation of high latitude conditions by using a quantitative scorecard approach that leverages NASA’s extensive satellite record. Investigations are conducted using the NASA Goddard Earth Observing System model (GEOS) and its data assimilation system (GEOS DAS). Through a quantitative identification of process deficiencies, bias reductions in key surface variables, including temperature and precipitation over cryospheric surfaces, may be achieved. The polarMERRA project additionally seeks to identify additional data sources for use in the GEOS DAS, and to incorporate data and parameterization improvements into future model and reanalysis products for scientific and stakeholder use.

Lauren C Andrews↗

Regional Impacts of Ice Sheet Surface Representation Improvements in an Earth System Model

Ice sheet surface conditions and their subsequent impact on surface mass balance play an important role in ice sheet dynamics and ice sheet interaction with the overlying atmosphere and surrounding ocean. The NASA Global Modeling and Assimilation Office’s (GMAO) Goddard Earth Observing System (GEOS) ecosystem of models and reanalyses model ice sheet and glacier surface mass balance with a moderately complex snow and ice module, including prognostic surface albedo evolution, fractional snow cover, and snowpack hydrology and meltwater retention. With this configuration, Greenland near surface air temperatures and ice sheet melt are well represented in MERRA-2, GMAO’s current atmospheric reanalysis; however, the similarity to observations in some regions is due, in part, to compensating biases in surface energy budget. Recent improvements to GEOS’ cloud microphysics and longwave radiation parameterizations reduced biases in surface net longwave radiation and exacerbated surface net shortwave radiation biases, effectively increasing near surface temperatures and melting, particularly during the summer months. Here, we demonstrate the Earth system response to an improved representation of the ice sheet surface focused on minimizing the surface energy biases and increasing spatial and temporal representativeness. New changes to the surface albedo parameterization, including the introduction of a prognostic aerosol darkening scheme, help mitigate shortwave biases, while modifications to surface roughness characteristics improve near surface sensible heat fluxes. Focusing on the Arctic and Greenland, we examine the combined impact on these changes on regional weather conditions and ice sheet surface mass balance using a high resolution atmospheric only experiment, and we use a coupled model configuration to assess the impact on oceanic conditions and subseasonal-to-seasonal prediction. Overall, these improvements increase modeled ice sheet surface realism and provide a strong basis for more accurate ice sheet mass balance in GMAO’s future reanalyses and forecasting systems.

Lauren C Andrews↗

Regional Impacts of Ice Sheet Surface Representation Improvements in an Earth System Model

Ice sheet surface conditions and their subsequent impact on surface mass balance play an important role in ice sheet dynamics and ice sheet interaction with the overlying atmosphere and surrounding ocean. The NASA Global Modeling and Assimilation Office’s (GMAO) Goddard Earth Observing System (GEOS) ecosystem of models and reanalyses model ice sheet and glacier surface mass balance with a moderately complex snow and ice module, including prognostic surface albedo evolution, fractional snow cover, and snowpack hydrology and meltwater retention. With this configuration, Greenland near surface air temperatures and ice sheet melt are well represented in MERRA-2, GMAO’s current atmospheric reanalysis; however, the similarity to observations in some regions is due, in part, to compensating biases in surface energy budget. Recent improvements to GEOS’ cloud microphysics and longwave radiation parameterizations reduced biases in surface net longwave radiation and exacerbated surface net shortwave radiation biases, effectively increasing near surface temperatures and melting, particularly during the summer months. Here, we demonstrate the Earth system response to an improved representation of the ice sheet surface focused on minimizing the surface energy biases and increasing spatial and temporal representativeness. New changes to the surface albedo parameterization, including the introduction of a prognostic aerosol darkening scheme, help mitigate shortwave biases, while modifications to surface roughness characteristics improve near surface sensible heat fluxes. Focusing on the Arctic and Greenland, we examine the combined impact on these changes on regional weather conditions and ice sheet surface mass balance using a high resolution atmospheric only experiment, and we use a coupled model configuration to assess the impact on oceanic conditions and subseasonal-to-seasonal prediction. Overall, these improvements increase modeled ice sheet surface realism and provide a strong basis for more accurate ice sheet mass balance in GMAO’s future reanalyses and forecasting systems.

Lauren C. Andrews↗

Variability of carbon monoxide in the Mars atmosphere

Recent measurements and reanalyses of past measurements are used to derive estimates of the variability of CO concentrations in the Martian atmosphere. The J = 1 - 2 rotational transition of CO in the microwave spectrum of Mars was observed in January, 1982, and abundances were calculated through surface and atmospheric temperature modeling. A column density of 4.6 + or - 2.0 x 10 to the 20th/sq cm is obtained, which is contrasted with values of 3.0 + or - 1.0 x 10 to the 20th/sq cm in 1980 and 1.7 x + or - 0.9 x 10 to the 20th/sq cm in 1975 obtained by reanalysis of the J = 0 - 1 transition observed by Good and Schloerb (1981) and Kakar et al. (1977), respectively. It is noted that due to the uncertainties in global atmospheric average temperatures, results are consistent with a CO variability of from 0 to 100 percent over a time scale of several years. It is suggested that simultaneous measurements of the J = 0 - 1 and J = 1 - 2 transitions may provide stronger constraints on the CO abundance and atmospheric temperature of Mars.

Clancy, R. T.↗

Regional scale simulation of atmospheric structure using VAS data

Investigations with both Subsynoptic Scale Model (SSM) and Limited Area Mesoscale Prediction System (LAMPS) are directed at model static initializations using Visible Infrared Spin Scan Radiometer Atmospheric Sounder (VAS) derived temperatures and moistures and also satellite derived winds. The general approach for each day's case study is similar. A 12-hour control forecast (no satellite data) is made 12Z-OZ, using a radiosonde analysis at 12Z as the initial field. For the satellite data experiments, reanalyses are performed at intermediate times, using VAS soundings or VAS soundings in combination with satellite winds. Forecasts are subsequently made from the reanalysis time(s) to the same validation hour (OZ) as the control forecasts. Comparisons are then made between the control forecast and satellite experiments.

Disk, G.↗

The cosmic ray spectrum above 10(17) eV

The final analysis of the data obtained by the Sydney University Giant Airshower Recorder (SUGAR) is presented. The data has been reanalysed to take into account the effects of afterpulsing in the photomultiplier tubes. Event data was used to produce a spectrum of equivalent vertical muon number and from this a model dependent primary energy spectrum was obtained. These spectra show good evidence for the Ankle: a flattening at 10(19) eV. There is no sign of the cut-off which would be expected from the effects of the universal black body radiation.

Winn, M. M.↗