38 research outputs found
Technical report series on global modeling and data assimilation. Volume 5: Documentation of the AIRES/GEOS dynamical core, version 2
A detailed description of the numerical formulation of Version 2 of the ARIES/GEOS 'dynamical core' is presented. This code is a nearly 'plug-compatible' dynamics for use in atmospheric general circulation models (GCMs). It is a finite difference model on a staggered latitude-longitude C-grid. It uses second-order differences for all terms except the advection of vorticity by the rotation part of the flow, which is done at fourth-order accuracy. This dynamical core is currently being used in the climate (ARIES) and data assimilation (GEOS) GCMs at Goddard
Maintaining Atmospheric Mass and Water Balance Within Reanalysis
This report describes the modifications implemented into the Goddard Earth Observing System Version-5 (GEOS-5) Atmospheric Data Assimilation System (ADAS) to maintain global conservation of dry atmospheric mass as well as to preserve the model balance of globally integrated precipitation and surface evaporation during reanalysis. Section 1 begins with a review of these global quantities from four current reanalysis efforts. Section 2 introduces the modifications necessary to preserve these constraints within the atmospheric general circulation model (AGCM), the Gridpoint Statistical Interpolation (GSI) analysis procedure, and the Incremental Analysis Update (IAU) algorithm. Section 3 presents experiments quantifying the impact of the new procedure. Section 4 shows preliminary results from its use within the GMAO MERRA-2 Reanalysis project. Section 5 concludes with a summary
Technical report series on global modeling and data assimilation. Volume 1: Documentation of the Goddard Earth Observing System (GEOS) General Circulation Model, version 1
This technical report documents Version 1 of the Goddard Earth Observing System (GEOS) General Circulation Model (GCM). The GEOS-1 GCM is being used by NASA's Data Assimilation Office (DAO) to produce multiyear data sets for climate research. This report provides a documentation of the model components used in the GEOS-1 GCM, a complete description of model diagnostics available, and a User's Guide to facilitate GEOS-1 GCM experiments
Large-Scale Atmospheric Transport in GEOS Replay Simulations
Offline chemical transport models (CTMs) have traditionally been used to perform studies of atmospheric chemistry in a fixed dynamical environment. An alternative to using CTMs is to constrain the flow in a general circulation model using winds from meteorological analyses. The Goddard Earth Observing System (GEOS) "replay" approach involves reading in analyzed fields every six hours and recomputing the analysis increments, which are applied as a forcing to the meteorology at every model time step. Unlike in CTM, all of the subgrid-scale processes are recalculated on-line so that they are consistent with the large-scale analysis fields, similar in spirit to "nudged" simulations, in which the online meteorology is relaxed to the analysis. Here we compare the transport of idealized tracers in different replay simulations constrained with meteorological fields taken from The Modern-Era Retrospective Analysis for Research and Applications, Version 2 (MERRA-2). We show that there are substantial differences in their large-scale stratospheric transport, depending on whether analysis fields or assimilated fields are used. Replay simulations constrained with the instantaneous analysis fields produce stratospheric mean age values that are up to 30% too young relative to observations; by comparison, simulations constrained with the time-averaged assimilated fields produce more credible stratospheric transport. Our study indicates that care should be taken to correctly configure the model when the replay technique is used to simulate stratospheric composition
Comparing Flood-Producing Catastrophic Hurricanes in GEOS: A Study of Juan 1985 and Harvey 2017
Hurricanes Juan (1985) and Harvey (2017) bear remarkable similarities because of the exceptional amount of accumulated precipitation over coastal regions in the Gulf, consequent to the peculiar and rare situation in which a hurricane stalls for several days on a coastal area, displaying minimal progress. This is particularly dangerous because: a) terrestrial drag affect only part of the circulation so that weakening is not as fast as when there is rapid motion inland; b) a part of the circulation remains over water and as such low-level moisture flux convergence is maintained; c) precipitation lasts for a very long time over the same area; d) moisture recycling from flooded areas further contributes to maintaining flood conditions. From the forecast perspective, hurricanes with track singularities such as stalling points, loops, or multiple landfalls over the same region are exceptionally difficult to predict because they result from small variations of opposing mechanisms and are therefore intrinsically unstable. Intensity fluctuations and track singularities are strongly coupled. The representation of fluxes undergoing differential drag (ocean vs. land) inside the same circulation is critically dependent upon resolution. Recent modeling development efforts within the NASA Global Modeling and Assimilation Office (GMAO) have brought the horizontal resolution of the global cubed-sphere Goddard Earth Observing System (GEOS) atmospheric model down to 12 km. This has led to innovative findings in the field of tropical cyclone (TC) modeling within a global framework. The GEOS atmospheric model and data assimilation system (DAS), used to produce the MERRA-2 (Modern Era Retrospective analysis for Research and Applications) reanalysis, has a "Replay" capability that adds a forcing term to the model equations to constrain it to follow a specified trajectory. This capability is used to produce global, downscaled versions of the 50 km MERRA-2 reanalysis at a 12.5 km resolution. The full DAS and replay modes both use the Incremental Analysis Update or IAU for applying an analysis increment to the model's background state four times each day at 00Z, 06Z, 12Z and 18Z. In a full DAS mode, a six-hour free run of the model is made (Predictor Step), writing states that will be used as backgrounds or first-guesses for the analysis centered on 00Z. The Analysis system then produces its best estimate of the 00Z state by blending backgrounds with observations. The difference between this best estimate and the model's original first guess at 00Z is treated as an analysis increment and added uniformly as a forcing term to the model during a second, now forced, run of the 12Z to 03Z window, which produces the assimilated results. Replay mode is very similar, except only the 00Z background is needed from the predictor segment, which is thus half as long, and instead of doing an analysis, a previously analyzed state at 00Z is read in. The increment is then computed and applied to the corrector step in the same way as in the original DAS. When downscaling horizontal spatial filters can be applied to the "analysis" increment used to force the Corrector step, so as to reproduce the full meteorological state or selectively choose or discard certain forcings. In this presentation, full global hindcast simulations of Juan and operational forecasts of Harvey are compared
Land Boundary Conditions for the Goddard Earth Observing System Model Version 5 (GEOS-5) Climate Modeling System: Recent Updates and Data File Descriptions
The Earths land surface boundary conditions in the Goddard Earth Observing System version 5 (GEOS-5) modeling system were updated using recent high spatial and temporal resolution global data products. The updates include: (i) construction of a global 10-arcsec land-ocean lakes-ice mask; (ii) incorporation of a 10-arcsec Globcover 2009 land cover dataset; (iii) implementation of Level 12 Pfafstetter hydrologic catchments; (iv) use of hybridized SRTM global topography data; (v) construction of the HWSDv1.21-STATSGO2 merged global 30 arc second soil mineral and carbon data in conjunction with a highly-refined soil classification system; (vi) production of diffuse visible and near-infrared 8-day MODIS albedo climatologies at 30-arcsec from the period 2001-2011; and (vii) production of the GEOLAND2 and MODIS merged 8-day LAI climatology at 30-arcsec for GEOS-5. The global data sets were preprocessed and used to construct global raster data files for the software (mkCatchParam) that computes parameters on catchment-tiles for various atmospheric grids. The updates also include a few bug fixes in mkCatchParam, as well as changes (improvements in algorithms, etc.) to mkCatchParam that allow it to produce tile-space parameters efficiently for high resolution AGCM grids. The update process also includes the construction of data files describing the vegetation type fractions, soil background albedo, nitrogen deposition and mean annual 2m air temperature to be used with the future Catchment CN model and the global stream channel network to be used with the future global runoff routing model. This report provides detailed descriptions of the data production process and data file format of each updated data set
Solar Occultation Satellite Data and Derived Meteorological Products: Sampling Issues and Comparisons with Aura MLS
Derived Meteorological Products (DMPs, including potential temperature (theta), potential vorticity, equivalent latitude (EqL), horizontal winds and tropopause locations) have been produced for the locations and times of measurements by several solar occultation (SO) instruments and the Aura Microwave Limb Sounder (MLS). DMPs are calculated from several meteorological analyses for the Atmospheric Chemistry Experiment-Fourier Transform Spectrometer, Stratospheric Aerosol and Gas Experiment II and III, Halogen Occultation Experiment, and Polar Ozone and Aerosol Measurement II and III SO instruments and MLS. Time-series comparisons of MLS version 1.5 and SO data using DMPs show good qualitative agreement in time evolution of O3, N2O, H20, CO, HNO3, HCl and temperature; quantitative agreement is good in most cases. EqL-coordinate comparisons of MLS version 2.2 and SO data show good quantitative agreement throughout the stratosphere for most of these species, with significant biases for a few species in localized regions. Comparisons in EqL coordinates of MLS and SO data, and of SO data with geographically coincident MLS data provide insight into where and how sampling effects are important in interpretation of the sparse SO data, thus assisting in fully utilizing the SO data in scientific studies and comparisons with other sparse datasets. The DMPs are valuable for scientific studies and to facilitate validation of non-coincident measurements
Effect of SGLT2 inhibitors on stroke and atrial fibrillation in diabetic kidney disease: Results from the CREDENCE trial and meta-analysis
BACKGROUND AND PURPOSE: Chronic kidney disease with reduced estimated glomerular filtration rate or elevated albuminuria increases risk for ischemic and hemorrhagic stroke. This study assessed the effects of sodium glucose cotransporter 2 inhibitors (SGLT2i) on stroke and atrial fibrillation/flutter (AF/AFL) from CREDENCE (Canagliflozin and Renal Events in Diabetes With Established Nephropathy Clinical Evaluation) and a meta-Analysis of large cardiovascular outcome trials (CVOTs) of SGLT2i in type 2 diabetes mellitus. METHODS: CREDENCE randomized 4401 participants with type 2 diabetes mellitus and chronic kidney disease to canagliflozin or placebo. Post hoc, we estimated effects on fatal or nonfatal stroke, stroke subtypes, and intermediate markers of stroke risk including AF/AFL. Stroke and AF/AFL data from 3 other completed large CVOTs and CREDENCE were pooled using random-effects meta-Analysis. RESULTS: In CREDENCE, 142 participants experienced a stroke during follow-up (10.9/1000 patient-years with canagliflozin, 14.2/1000 patient-years with placebo; hazard ratio [HR], 0.77 [95% CI, 0.55-1.08]). Effects by stroke subtypes were: ischemic (HR, 0.88 [95% CI, 0.61-1.28]; n=111), hemorrhagic (HR, 0.50 [95% CI, 0.19-1.32]; n=18), and undetermined (HR, 0.54 [95% CI, 0.20-1.46]; n=17). There was no clear effect on AF/AFL (HR, 0.76 [95% CI, 0.53-1.10]; n=115). The overall effects in the 4 CVOTs combined were: Total stroke (HRpooled, 0.96 [95% CI, 0.82-1.12]), ischemic stroke (HRpooled, 1.01 [95% CI, 0.89-1.14]), hemorrhagic stroke (HRpooled, 0.50 [95% CI, 0.30-0.83]), undetermined stroke (HRpooled, 0.86 [95% CI, 0.49-1.51]), and AF/AFL (HRpooled, 0.81 [95% CI, 0.71-0.93]). There was evidence that SGLT2i effects on total stroke varied by baseline estimated glomerular filtration rate (P=0.01), with protection in the lowest estimated glomerular filtration rate (45 mL/min/1.73 m2]) subgroup (HRpooled, 0.50 [95% CI, 0.31-0.79]). CONCLUSIONS: Although we found no clear effect of SGLT2i on total stroke in CREDENCE or across trials combined, there was some evidence of benefit in preventing hemorrhagic stroke and AF/AFL, as well as total stroke for those with lowest estimated glomerular filtration rate. Future research should focus on confirming these data and exploring potential mechanisms