Graphics Processing Unit (GPU) Acceleration of the Goddard Earth Observing System Atmospheric Model

NASA Technical Reports Server (NTRS)

Putnam, Williama

2011-01-01

The Goddard Earth Observing System 5 (GEOS-5) is the atmospheric model used by the Global Modeling and Assimilation Office (GMAO) for a variety of applications, from long-term climate prediction at relatively coarse resolution, to data assimilation and numerical weather prediction, to very high-resolution cloud-resolving simulations. GEOS-5 is being ported to a graphics processing unit (GPU) cluster at the NASA Center for Climate Simulation (NCCS). By utilizing GPU co-processor technology, we expect to increase the throughput of GEOS-5 by at least an order of magnitude, and accelerate the process of scientific exploration across all scales of global modeling, including: The large-scale, high-end application of non-hydrostatic, global, cloud-resolving modeling at 10- to I-kilometer (km) global resolutions Intermediate-resolution seasonal climate and weather prediction at 50- to 25-km on small clusters of GPUs Long-range, coarse-resolution climate modeling, enabled on a small box of GPUs for the individual researcher After being ported to the GPU cluster, the primary physics components and the dynamical core of GEOS-5 have demonstrated a potential speedup of 15-40 times over conventional processor cores. Performance improvements of this magnitude reduce the required scalability of 1-km, global, cloud-resolving models from an unfathomable 6 million cores to an attainable 200,000 GPU-enabled cores.

Vertical resolution of baroclinic modes in global ocean models

NASA Astrophysics Data System (ADS)

Stewart, K. D.; Hogg, A. McC.; Griffies, S. M.; Heerdegen, A. P.; Ward, M. L.; Spence, P.; England, M. H.

2017-05-01

Improvements in the horizontal resolution of global ocean models, motivated by the horizontal resolution requirements for specific flow features, has advanced modelling capabilities into the dynamical regime dominated by mesoscale variability. In contrast, the choice of the vertical grid remains a subjective choice, and it is not clear that efforts to improve vertical resolution adequately support their horizontal counterparts. Indeed, considering that the bulk of the vertical ocean dynamics (including convection) are parameterized, it is not immediately obvious what the vertical grid is supposed to resolve. Here, we propose that the primary purpose of the vertical grid in a hydrostatic ocean model is to resolve the vertical structure of horizontal flows, rather than to resolve vertical motion. With this principle we construct vertical grids based on their abilities to represent baroclinic modal structures commensurate with the theoretical capabilities of a given horizontal grid. This approach is designed to ensure that the vertical grids of global ocean models complement (and, importantly, to not undermine) the resolution capabilities of the horizontal grid. We find that for z-coordinate global ocean models, at least 50 well-positioned vertical levels are required to resolve the first baroclinic mode, with an additional 25 levels per subsequent mode. High-resolution ocean-sea ice simulations are used to illustrate some of the dynamical enhancements gained by improving the vertical resolution of a 1/10° global ocean model. These enhancements include substantial increases in the sea surface height variance (∼30% increase south of 40°S), the barotropic and baroclinic eddy kinetic energies (up to 200% increase on and surrounding the Antarctic continental shelf and slopes), and the overturning streamfunction in potential density space (near-tripling of the Antarctic Bottom Water cell at 65°S).

High resolution simulations of aerosol microphysics in a global and regionally nested chemical transport model

NASA Astrophysics Data System (ADS)

Adams, P. J.; Marks, M.

2015-12-01

The aerosol indirect effect is the largest source of forcing uncertainty in current climate models. This effect arises from the influence of aerosols on the reflective properties and lifetimes of clouds, and its magnitude depends on how many particles can serve as cloud droplet formation sites. Assessing levels of this subset of particles (cloud condensation nuclei, or CCN) requires knowledge of aerosol levels and their global distribution, size distributions, and composition. A key tool necessary to advance our understanding of CCN is the use of global aerosol microphysical models, which simulate the processes that control aerosol size distributions: nucleation, condensation/evaporation, and coagulation. Previous studies have found important differences in CO (Chen, D. et al., 2009) and ozone (Jang, J., 1995) modeled at different spatial resolutions, and it is reasonable to believe that short-lived, spatially-variable aerosol species will be similarly - or more - susceptible to model resolution effects. The goal of this study is to determine how CCN levels and spatial distributions change as simulations are run at higher spatial resolution - specifically, to evaluate how sensitive the model is to grid size, and how this affects comparisons against observations. Higher resolution simulations are necessary supports for model/measurement synergy. Simulations were performed using the global chemical transport model GEOS-Chem (v9-02). The years 2008 and 2009 were simulated at 4ox5o and 2ox2.5o globally and at 0.5ox0.667o over Europe and North America. Results were evaluated against surface-based particle size distribution measurements from the European Supersites for Atmospheric Aerosol Research project. The fine-resolution model simulates more spatial and temporal variability in ultrafine levels, and better resolves topography. Results suggest that the coarse model predicts systematically lower ultrafine levels than does the fine-resolution model. Significant differences are also evident with respect to model-measurement comparisons, and will be discussed.

Global high-resolution simulations of tropospheric nitrogen dioxide using CHASER V4.0

NASA Astrophysics Data System (ADS)

Sekiya, Takashi; Miyazaki, Kazuyuki; Ogochi, Koji; Sudo, Kengo; Takigawa, Masayuki

2018-03-01

We evaluate global tropospheric nitrogen dioxide (NO2) simulations using the CHASER V4.0 global chemical transport model (CTM) at horizontal resolutions of 0.56, 1.1, and 2.8°. Model evaluation was conducted using satellite tropospheric NO2 retrievals from the Ozone Monitoring Instrument (OMI) and the Global Ozone Monitoring Experiment-2 (GOME-2) and aircraft observations from the 2014 Front Range Air Pollution and Photochemistry Experiment (FRAPPÉ). Agreement against satellite retrievals improved greatly at 1.1 and 0.56° resolutions (compared to 2.8° resolution) over polluted and biomass burning regions. The 1.1° simulation generally captured the regional distribution of the tropospheric NO2 column well, whereas 0.56° resolution was necessary to improve the model performance over areas with strong local sources, with mean bias reductions of 67 % over Beijing and 73 % over San Francisco in summer. Validation using aircraft observations indicated that high-resolution simulations reduced negative NO2 biases below 700 hPa over the Denver metropolitan area. These improvements in high-resolution simulations were attributable to (1) closer spatial representativeness between simulations and observations and (2) better representation of large-scale concentration fields (i.e., at 2.8°) through the consideration of small-scale processes. Model evaluations conducted at 0.5 and 2.8° bin grids indicated that the contributions of both these processes were comparable over most polluted regions, whereas the latter effect (2) made a larger contribution over eastern China and biomass burning areas. The evaluations presented in this paper demonstrate the potential of using a high-resolution global CTM for studying megacity-scale air pollutants across the entire globe, potentially also contributing to global satellite retrievals and chemical data assimilation.

Precipitation response to solar geoengineering in a high-resolution tropical-cyclone permitting coupled general circulation model

NASA Astrophysics Data System (ADS)

Irvine, P. J.; Keith, D.; Dykema, J. A.; Vecchi, G. A.; Horowitz, L. W.

2016-12-01

Solar geoengineering may limit or even halt the rise in global-average surface temperatures. Evidence from the geoMIP model intercomparison project shows that idealized geoengineering can greatly reduce temperature changes on a region-by-region basis. If solar geoengineering is used to hold radiative forcing or surface temperatures constant in the face of rising CO2, then the global evaporation and precipitation rates will be reduced below pre-industrial. The spartial and frequency distribution of the precipitation response is, however, much less well understood. There is limited evidence that solar geoengineering may reduce extreme precipitation events more that it reduces mean precipitation, but that evidence is based on relatively course resolution models that may to a poor job representing the distribution of extreme precipitation in the current climate. The response of global and regional climate, as well as tropical cyclone (TC) activity, to increasing solar geoengineering is explored through experiments with climate models spanning a broad range of atmospheric resolutions. Solar geoengineering is represented by an idealized adjustment of the solar constant that roughly halves the rate of increase in radiative forcing in a scenario with increasing CO2 concentration. The coarsest resolution model has approximately a 2-degree global resolution, representative of the typical resolution of past GCMs used to explore global response to CO2 increase, and its response is compared to that of two tropical cyclone permitting GCMs of approximately 0.5 and 0.25 degree resolution (FLOR and HiFLOR). The models have exactly the same ocean and sea-ice components, as well as the same parameterizations and parameter settings. These high-resolution models are used for real-time seasonal prediction, providing a unified framework for seasonal-to-multidecadal climate modeling. We assess the extreme precipitation response, comparing the frequency distribution of extreme events with and without solar geoengineering. We compare our results to two prior studies of the response of climate extremes to solar geoengineering.

The importance of vertical resolution in the free troposphere for modeling intercontinental plumes

NASA Astrophysics Data System (ADS)

Zhuang, Jiawei; Jacob, Daniel J.; Eastham, Sebastian D.

2018-05-01

Chemical plumes in the free troposphere can preserve their identity for more than a week as they are transported on intercontinental scales. Current global models cannot reproduce this transport. The plumes dilute far too rapidly due to numerical diffusion in sheared flow. We show how model accuracy can be limited by either horizontal resolution (Δx) or vertical resolution (Δz). Balancing horizontal and vertical numerical diffusion, and weighing computational cost, implies an optimal grid resolution ratio (Δx / Δz)opt ˜ 1000 for simulating the plumes. This is considerably higher than current global models (Δx / Δz ˜ 20) and explains the rapid plume dilution in the models as caused by insufficient vertical resolution. Plume simulations with the Geophysical Fluid Dynamics Laboratory Finite-Volume Cubed-Sphere Dynamical Core (GFDL-FV3) over a range of horizontal and vertical grid resolutions confirm this limiting behavior. Our highest-resolution simulation (Δx ≈ 25 km, Δz ≈ 80 m) preserves the maximum mixing ratio in the plume to within 35 % after 8 days in strongly sheared flow, a drastic improvement over current models. Adding free tropospheric vertical levels in global models is computationally inexpensive and would also improve the simulation of water vapor.

Prospects for improving the representation of coastal and shelf seas in global ocean models

NASA Astrophysics Data System (ADS)

Holt, Jason; Hyder, Patrick; Ashworth, Mike; Harle, James; Hewitt, Helene T.; Liu, Hedong; New, Adrian L.; Pickles, Stephen; Porter, Andrew; Popova, Ekaterina; Icarus Allen, J.; Siddorn, John; Wood, Richard

2017-02-01

Accurately representing coastal and shelf seas in global ocean models represents one of the grand challenges of Earth system science. They are regions of immense societal importance through the goods and services they provide, hazards they pose and their role in global-scale processes and cycles, e.g. carbon fluxes and dense water formation. However, they are poorly represented in the current generation of global ocean models. In this contribution, we aim to briefly characterise the problem, and then to identify the important physical processes, and their scales, needed to address this issue in the context of the options available to resolve these scales globally and the evolving computational landscape.We find barotropic and topographic scales are well resolved by the current state-of-the-art model resolutions, e.g. nominal 1/12°, and still reasonably well resolved at 1/4°; here, the focus is on process representation. We identify tides, vertical coordinates, river inflows and mixing schemes as four areas where modelling approaches can readily be transferred from regional to global modelling with substantial benefit. In terms of finer-scale processes, we find that a 1/12° global model resolves the first baroclinic Rossby radius for only ˜ 8 % of regions < 500 m deep, but this increases to ˜ 70 % for a 1/72° model, so resolving scales globally requires substantially finer resolution than the current state of the art.We quantify the benefit of improved resolution and process representation using 1/12° global- and basin-scale northern North Atlantic nucleus for a European model of the ocean (NEMO) simulations; the latter includes tides and a k-ɛ vertical mixing scheme. These are compared with global stratification observations and 19 models from CMIP5. In terms of correlation and basin-wide rms error, the high-resolution models outperform all these CMIP5 models. The model with tides shows improved seasonal cycles compared to the high-resolution model without tides. The benefits of resolution are particularly apparent in eastern boundary upwelling zones.To explore the balance between the size of a globally refined model and that of multiscale modelling options (e.g. finite element, finite volume or a two-way nesting approach), we consider a simple scale analysis and a conceptual grid refining approach. We put this analysis in the context of evolving computer systems, discussing model turnaround time, scalability and resource costs. Using a simple cost model compared to a reference configuration (taken to be a 1/4° global model in 2011) and the increasing performance of the UK Research Councils' computer facility, we estimate an unstructured mesh multiscale approach, resolving process scales down to 1.5 km, would use a comparable share of the computer resource by 2021, the two-way nested multiscale approach by 2022, and a 1/72° global model by 2026. However, we also note that a 1/12° global model would not have a comparable computational cost to a 1° global model in 2017 until 2027. Hence, we conclude that for computationally expensive models (e.g. for oceanographic research or operational oceanography), resolving scales to ˜ 1.5 km would be routinely practical in about a decade given substantial effort on numerical and computational development. For complex Earth system models, this extends to about 2 decades, suggesting the focus here needs to be on improved process parameterisation to meet these challenges.

Air Quality Forecasts Using the NASA GEOS Model: A Unified Tool from Local to Global Scales

NASA Technical Reports Server (NTRS)

Knowland, E. Emma; Keller, Christoph; Nielsen, J. Eric; Orbe, Clara; Ott, Lesley; Pawson, Steven; Saunders, Emily; Duncan, Bryan; Cook, Melanie; Liu, Junhua;

A framework for global river flood risk assessments

NASA Astrophysics Data System (ADS)

Winsemius, H. C.; Van Beek, L. P. H.; Jongman, B.; Ward, P. J.; Bouwman, A.

2012-08-01

There is an increasing need for strategic global assessments of flood risks in current and future conditions. In this paper, we propose a framework for global flood risk assessment for river floods, which can be applied in current conditions, as well as in future conditions due to climate and socio-economic changes. The framework's goal is to establish flood hazard and impact estimates at a high enough resolution to allow for their combination into a risk estimate. The framework estimates hazard at high resolution (~1 km2) using global forcing datasets of the current (or in scenario mode, future) climate, a global hydrological model, a global flood routing model, and importantly, a flood extent downscaling routine. The second component of the framework combines hazard with flood impact models at the same resolution (e.g. damage, affected GDP, and affected population) to establish indicators for flood risk (e.g. annual expected damage, affected GDP, and affected population). The framework has been applied using the global hydrological model PCR-GLOBWB, which includes an optional global flood routing model DynRout, combined with scenarios from the Integrated Model to Assess the Global Environment (IMAGE). We performed downscaling of the hazard probability distributions to 1 km2 resolution with a new downscaling algorithm, applied on Bangladesh as a first case-study application area. We demonstrate the risk assessment approach in Bangladesh based on GDP per capita data, population, and land use maps for 2010 and 2050. Validation of the hazard and damage estimates has been performed using the Dartmouth Flood Observatory database and damage estimates from the EM-DAT database and World Bank sources. We discuss and show sensitivities of the estimated risks with regard to the use of different climate input sets, decisions made in the downscaling algorithm, and different approaches to establish impact models.

Towards multi-resolution global climate modeling with ECHAM6-FESOM. Part II: climate variability

NASA Astrophysics Data System (ADS)

Rackow, T.; Goessling, H. F.; Jung, T.; Sidorenko, D.; Semmler, T.; Barbi, D.; Handorf, D.

2018-04-01

This study forms part II of two papers describing ECHAM6-FESOM, a newly established global climate model with a unique multi-resolution sea ice-ocean component. While part I deals with the model description and the mean climate state, here we examine the internal climate variability of the model under constant present-day (1990) conditions. We (1) assess the internal variations in the model in terms of objective variability performance indices, (2) analyze variations in global mean surface temperature and put them in context to variations in the observed record, with particular emphasis on the recent warming slowdown, (3) analyze and validate the most common atmospheric and oceanic variability patterns, (4) diagnose the potential predictability of various climate indices, and (5) put the multi-resolution approach to the test by comparing two setups that differ only in oceanic resolution in the equatorial belt, where one ocean mesh keeps the coarse 1° resolution applied in the adjacent open-ocean regions and the other mesh is gradually refined to 0.25°. Objective variability performance indices show that, in the considered setups, ECHAM6-FESOM performs overall favourably compared to five well-established climate models. Internal variations of the global mean surface temperature in the model are consistent with observed fluctuations and suggest that the recent warming slowdown can be explained as a once-in-one-hundred-years event caused by internal climate variability; periods of strong cooling in the model (`hiatus' analogs) are mainly associated with ENSO-related variability and to a lesser degree also to PDO shifts, with the AMO playing a minor role. Common atmospheric and oceanic variability patterns are simulated largely consistent with their real counterparts. Typical deficits also found in other models at similar resolutions remain, in particular too weak non-seasonal variability of SSTs over large parts of the ocean and episodic periods of almost absent deep-water formation in the Labrador Sea, resulting in overestimated North Atlantic SST variability. Concerning the influence of locally (isotropically) increased resolution, the ENSO pattern and index statistics improve significantly with higher resolution around the equator, illustrating the potential of the novel unstructured-mesh method for global climate modeling.

Towards Direct Simulation of Future Tropical Cyclone Statistics in a High-Resolution Global Atmospheric Model

DOE PAGES

Wehner, Michael F.; Bala, G.; Duffy, Phillip; ...

2010-01-01

We present a set of high-resolution global atmospheric general circulation model (AGCM) simulations focusing on the model's ability to represent tropical storms and their statistics. We find that the model produces storms of hurricane strength with realistic dynamical features. We also find that tropical storm statistics are reasonable, both globally and in the north Atlantic, when compared to recent observations. The sensitivity of simulated tropical storm statistics to increases in sea surface temperature (SST) is also investigated, revealing that a credible late 21st century SST increase produced increases in simulated tropical storm numbers and intensities in all ocean basins. Whilemore » this paper supports previous high-resolution model and theoretical findings that the frequency of very intense storms will increase in a warmer climate, it differs notably from previous medium and high-resolution model studies that show a global reduction in total tropical storm frequency. However, we are quick to point out that this particular model finding remains speculative due to a lack of radiative forcing changes in our time-slice experiments as well as a focus on the Northern hemisphere tropical storm seasons.« less

A 4.5 km resolution Arctic Ocean simulation with the global multi-resolution model FESOM 1.4

NASA Astrophysics Data System (ADS)

Wang, Qiang; Wekerle, Claudia; Danilov, Sergey; Wang, Xuezhu; Jung, Thomas

2018-04-01

In the framework of developing a global modeling system which can facilitate modeling studies on Arctic Ocean and high- to midlatitude linkage, we evaluate the Arctic Ocean simulated by the multi-resolution Finite Element Sea ice-Ocean Model (FESOM). To explore the value of using high horizontal resolution for Arctic Ocean modeling, we use two global meshes differing in the horizontal resolution only in the Arctic Ocean (24 km vs. 4.5 km). The high resolution significantly improves the model's representation of the Arctic Ocean. The most pronounced improvement is in the Arctic intermediate layer, in terms of both Atlantic Water (AW) mean state and variability. The deepening and thickening bias of the AW layer, a common issue found in coarse-resolution simulations, is significantly alleviated by using higher resolution. The topographic steering of the AW is stronger and the seasonal and interannual temperature variability along the ocean bottom topography is enhanced in the high-resolution simulation. The high resolution also improves the ocean surface circulation, mainly through a better representation of the narrow straits in the Canadian Arctic Archipelago (CAA). The representation of CAA throughflow not only influences the release of water masses through the other gateways but also the circulation pathways inside the Arctic Ocean. However, the mean state and variability of Arctic freshwater content and the variability of freshwater transport through the Arctic gateways appear not to be very sensitive to the increase in resolution employed here. By highlighting the issues that are independent of model resolution, we address that other efforts including the improvement of parameterizations are still required.

Air Quality Forecasts Using the NASA GEOS Model

NASA Technical Reports Server (NTRS)

Keller, Christoph A.; Knowland, K. Emma; Nielsen, Jon E.; Orbe, Clara; Ott, Lesley; Pawson, Steven; Saunders, Emily; Duncan, Bryan; Follette-Cook, Melanie; Liu, Junhua;

Spectral characteristics of mid-latitude continental convection from a global variable-resolution Voronoi-mesh atmospheric model

NASA Astrophysics Data System (ADS)

Wong, M.; Skamarock, W. C.

2015-12-01

Global numerical weather forecast tests were performed using the global nonhydrostatic atmospheric model, Model for Prediction Across Scales (MPAS), for the NOAA Storm Prediction Center 2015 Spring Forecast Experiment (May 2015) and the Plains Elevated Convection at Night (PECAN) field campaign (June to mid-July 2015). These two sets of forecasts were performed on 50-to-3 km and 15-to-3 km smoothly-varying horizontal meshes, respectively. Both variable-resolution meshes have nominal convection-permitting 3-km grid spacing over the entire continental US. Here we evaluate the limited-area (vs. global) spectra from these NWP simulations. We will show the simulated spectral characteristics of total kinetic energy, vertical velocity variance, and precipitation during these spring and summer periods when diurnal continental convection is most active over central US. Spectral characteristics of a high-resolution global 3-km simulation (essentially no nesting) from the 20 May 2013 Moore, OK tornado case are also shown. These characteristics include spectral scaling, shape, and anisotropy, as well as the effective resolution of continental convection representation in MPAS.

Hurricane Forecasting with the High-resolution NASA Finite-volume General Circulation Model

NASA Technical Reports Server (NTRS)

Atlas, R.; Reale, O.; Shen, B.-W.; Lin, S.-J.; Chern, J.-D.; Putman, W.; Lee, T.; Yeh, K.-S.; Bosilovich, M.; Radakovich, J.

2004-01-01

A high-resolution finite-volume General Circulation Model (fvGCM), resulting from a development effort of more than ten years, is now being run operationally at the NASA Goddard Space Flight Center and Ames Research Center. The model is based on a finite-volume dynamical core with terrain-following Lagrangian control-volume discretization and performs efficiently on massive parallel architectures. The computational efficiency allows simulations at a resolution of a quarter of a degree, which is double the resolution currently adopted by most global models in operational weather centers. Such fine global resolution brings us closer to overcoming a fundamental barrier in global atmospheric modeling for both weather and climate, because tropical cyclones and even tropical convective clusters can be more realistically represented. In this work, preliminary results of the fvGCM are shown. Fifteen simulations of four Atlantic tropical cyclones in 2002 and 2004 are chosen because of strong and varied difficulties presented to numerical weather forecasting. It is shown that the fvGCM, run at the resolution of a quarter of a degree, can produce very good forecasts of these tropical systems, adequately resolving problems like erratic track, abrupt recurvature, intense extratropical transition, multiple landfall and reintensification, and interaction among vortices.

Analyzing and leveraging self-similarity for variable resolution atmospheric models

NASA Astrophysics Data System (ADS)

O'Brien, Travis; Collins, William

2015-04-01

Variable resolution modeling techniques are rapidly becoming a popular strategy for achieving high resolution in a global atmospheric models without the computational cost of global high resolution. However, recent studies have demonstrated a variety of resolution-dependent, and seemingly artificial, features. We argue that the scaling properties of the atmosphere are key to understanding how the statistics of an atmospheric model should change with resolution. We provide two such examples. In the first example we show that the scaling properties of the cloud number distribution define how the ratio of resolved to unresolved clouds should increase with resolution. We show that the loss of resolved clouds, in the high resolution region of variable resolution simulations, with the Community Atmosphere Model version 4 (CAM4) is an artifact of the model's treatment of condensed water (this artifact is significantly reduced in CAM5). In the second example we show that the scaling properties of the horizontal velocity field, combined with the incompressibility assumption, necessarily result in an intensification of vertical mass flux as resolution increases. We show that such an increase is present in a wide variety of models, including CAM and the regional climate models of the ENSEMBLES intercomparision. We present theoretical arguments linking this increase to the intensification of precipitation with increasing resolution.

GIEMS-D3: A new long-term, dynamical, high-spatial resolution inundation extent dataset at global scale

NASA Astrophysics Data System (ADS)

Aires, Filipe; Miolane, Léo; Prigent, Catherine; Pham Duc, Binh; Papa, Fabrice; Fluet-Chouinard, Etienne; Lehner, Bernhard

2017-04-01

The Global Inundation Extent from Multi-Satellites (GIEMS) provides multi-year monthly variations of the global surface water extent at 25kmx25km resolution. It is derived from multiple satellite observations. Its spatial resolution is usually compatible with climate model outputs and with global land surface model grids but is clearly not adequate for local applications that require the characterization of small individual water bodies. There is today a strong demand for high-resolution inundation extent datasets, for a large variety of applications such as water management, regional hydrological modeling, or for the analysis of mosquitos-related diseases. A new procedure is introduced to downscale the GIEMS low spatial resolution inundations to a 3 arc second (90 m) dataset. The methodology is based on topography and hydrography information from the HydroSHEDS database. A new floodability index is adopted and an innovative smoothing procedure is developed to ensure the smooth transition, in the high-resolution maps, between the low-resolution boxes from GIEMS. Topography information is relevant for natural hydrology environments controlled by elevation, but is more limited in human-modified basins. However, the proposed downscaling approach is compatible with forthcoming fusion with other more pertinent satellite information in these difficult regions. The resulting GIEMS-D3 database is the only high spatial resolution inundation database available globally at the monthly time scale over the 1993-2007 period. GIEMS-D3 is assessed by analyzing its spatial and temporal variability, and evaluated by comparisons to other independent satellite observations from visible (Google Earth and Landsat), infrared (MODIS) and active microwave (SAR).

High Resolution Global Topography of Eros from NEAR Imaging and LIDAR Data

NASA Technical Reports Server (NTRS)

Gaskell, Robert W.; Konopliv, A.; Barnouin-Jha, O.; Scheeres, D.

2006-01-01

Principal Data Products: Ensemble of L-maps from SPC, Spacecraft state, Asteroid pole and rotation. Secondary Products: Global topography model, inertia tensor, gravity. Composite high resolution topography. Three dimensional image maps.

High-resolution Local Gravity Model of the South Pole of the Moon from GRAIL Extended Mission Data

NASA Technical Reports Server (NTRS)

Goossens, Sander Johannes; Sabaka, Terence J.; Nicholas, Joseph B.; Lemoine, Frank G.; Rowlands, David D.; Mazarico, Erwan; Neumann, Gregory A.; Smith, David E.; Zuber, Maria T.

2014-01-01

We estimated a high-resolution local gravity field model over the south pole of the Moon using data from the Gravity Recovery and Interior Laboratory's extended mission. Our solution consists of adjustments with respect to a global model expressed in spherical harmonics. The adjustments are expressed as gridded gravity anomalies with a resolution of 1/6deg by 1/6deg (equivalent to that of a degree and order 1080 model in spherical harmonics), covering a cap over the south pole with a radius of 40deg. The gravity anomalies have been estimated from a short-arc analysis using only Ka-band range-rate (KBRR) data over the area of interest. We apply a neighbor-smoothing constraint to our solution. Our local model removes striping present in the global model; it reduces the misfit to the KBRR data and improves correlations with topography to higher degrees than current global models.

Mars-solar wind interaction: LatHyS, an improved parallel 3-D multispecies hybrid model

NASA Astrophysics Data System (ADS)

Modolo, Ronan; Hess, Sebastien; Mancini, Marco; Leblanc, Francois; Chaufray, Jean-Yves; Brain, David; Leclercq, Ludivine; Esteban-Hernández, Rosa; Chanteur, Gerard; Weill, Philippe; González-Galindo, Francisco; Forget, Francois; Yagi, Manabu; Mazelle, Christian

2016-07-01

In order to better represent Mars-solar wind interaction, we present an unprecedented model achieving spatial resolution down to 50 km, a so far unexplored resolution for global kinetic models of the Martian ionized environment. Such resolution approaches the ionospheric plasma scale height. In practice, the model is derived from a first version described in Modolo et al. (2005). An important effort of parallelization has been conducted and is presented here. A better description of the ionosphere was also implemented including ionospheric chemistry, electrical conductivities, and a drag force modeling the ion-neutral collisions in the ionosphere. This new version of the code, named LatHyS (Latmos Hybrid Simulation), is here used to characterize the impact of various spatial resolutions on simulation results. In addition, and following a global model challenge effort, we present the results of simulation run for three cases which allow addressing the effect of the suprathermal corona and of the solar EUV activity on the magnetospheric plasma boundaries and on the global escape. Simulation results showed that global patterns are relatively similar for the different spatial resolution runs, but finest grid runs provide a better representation of the ionosphere and display more details of the planetary plasma dynamic. Simulation results suggest that a significant fraction of escaping O+ ions is originated from below 1200 km altitude.

Spatial Variability in Column CO2 Inferred from High Resolution GEOS-5 Global Model Simulations: Implications for Remote Sensing and Inversions

NASA Technical Reports Server (NTRS)

Ott, L.; Putman, B.; Collatz, J.; Gregg, W.

2012-01-01

Column CO2 observations from current and future remote sensing missions represent a major advancement in our understanding of the carbon cycle and are expected to help constrain source and sink distributions. However, data assimilation and inversion methods are challenged by the difference in scale of models and observations. OCO-2 footprints represent an area of several square kilometers while NASA s future ASCENDS lidar mission is likely to have an even smaller footprint. In contrast, the resolution of models used in global inversions are typically hundreds of kilometers wide and often cover areas that include combinations of land, ocean and coastal areas and areas of significant topographic, land cover, and population density variations. To improve understanding of scales of atmospheric CO2 variability and representativeness of satellite observations, we will present results from a global, 10-km simulation of meteorology and atmospheric CO2 distributions performed using NASA s GEOS-5 general circulation model. This resolution, typical of mesoscale atmospheric models, represents an order of magnitude increase in resolution over typical global simulations of atmospheric composition allowing new insight into small scale CO2 variations across a wide range of surface flux and meteorological conditions. The simulation includes high resolution flux datasets provided by NASA s Carbon Monitoring System Flux Pilot Project at half degree resolution that have been down-scaled to 10-km using remote sensing datasets. Probability distribution functions are calculated over larger areas more typical of global models (100-400 km) to characterize subgrid-scale variability in these models. Particular emphasis is placed on coastal regions and regions containing megacities and fires to evaluate the ability of coarse resolution models to represent these small scale features. Additionally, model output are sampled using averaging kernels characteristic of OCO-2 and ASCENDS measurement concepts to create realistic pseudo-datasets. Pseudo-data are averaged over coarse model grid cell areas to better understand the ability of measurements to characterize CO2 distributions and spatial gradients on both short (daily to weekly) and long (monthly to seasonal) time scales

Short-Range Prediction of Monsoon Precipitation by NCMRWF Regional Unified Model with Explicit Convection

NASA Astrophysics Data System (ADS)

Mamgain, Ashu; Rajagopal, E. N.; Mitra, A. K.; Webster, S.

2018-03-01

There are increasing efforts towards the prediction of high-impact weather systems and understanding of related dynamical and physical processes. High-resolution numerical model simulations can be used directly to model the impact at fine-scale details. Improvement in forecast accuracy can help in disaster management planning and execution. National Centre for Medium Range Weather Forecasting (NCMRWF) has implemented high-resolution regional unified modeling system with explicit convection embedded within coarser resolution global model with parameterized convection. The models configurations are based on UK Met Office unified seamless modeling system. Recent land use/land cover data (2012-2013) obtained from Indian Space Research Organisation (ISRO) are also used in model simulations. Results based on short-range forecast of both the global and regional models over India for a month indicate that convection-permitting simulations by the high-resolution regional model is able to reduce the dry bias over southern parts of West Coast and monsoon trough zone with more intense rainfall mainly towards northern parts of monsoon trough zone. Regional model with explicit convection has significantly improved the phase of the diurnal cycle of rainfall as compared to the global model. Results from two monsoon depression cases during study period show substantial improvement in details of rainfall pattern. Many categories in rainfall defined for operational forecast purposes by Indian forecasters are also well represented in case of convection-permitting high-resolution simulations. For the statistics of number of days within a range of rain categories between `No-Rain' and `Heavy Rain', the regional model is outperforming the global model in all the ranges. In the very heavy and extremely heavy categories, the regional simulations show overestimation of rainfall days. Global model with parameterized convection have tendency to overestimate the light rainfall days and underestimate the heavy rain days compared to the observation data.

GLOBATO: An enhanced global relief model at 30 arc-seconds resolution

NASA Astrophysics Data System (ADS)

O'Leary, V.; Amante, C.

2017-12-01

The National Centers for Environmental Information (NCEI), an office of the National Oceanic and Atmospheric Administration (NOAA), first developed a digital bathymetric and elevation model, ETOPO5, from publicly available data in 1993. For nearly 25 years, NCEI's ETOPO family of global relief models have supported research at a planetary scale, including tsunami forecasting, ocean circulation modeling, visualization of the seafloor, understanding geological phenomena, and aiding the development of other global and regional elevation models. GLOBATO (GLObal BAThymetry and TOpography) is now the most detailed version released by NCEI with a horizontal resolution of 30 arc-seconds and succeeds ETOPO1 with the inclusion of several new or updated data-sets for the seafloor as well as land areas. GLOBATO is a compilation of data derived from models of satellite measurements, ship depth soundings, and multibeam surveys, as well as regional models developed for Greenland and Antarctica. These data were converted from different formats, resolutions, spatial distributions, and projections into a single global model using GDAL v2.2 and MB-System v5.5. As with previous NCEI models, GLOBATO is available in two formats, "bedrock elevation" (measured as the base of major ice sheets) and "ice surface elevation" (measured as the surface of major ice sheets) which provides comprehensive topographic and bathymetric coverage between +- 90 degrees latitude and +- 180 degrees longitude. Adhering to best practices, GLOBATO, all related digital products, and any supporting documentation are available online through the NCEI data portal. These new, high resolution models will better support the variety of research ETOPO1 has made possible.

A framework for global river flood risk assessments

NASA Astrophysics Data System (ADS)

Winsemius, H. C.; Van Beek, L. P. H.; Jongman, B.; Ward, P. J.; Bouwman, A.

2013-05-01

There is an increasing need for strategic global assessments of flood risks in current and future conditions. In this paper, we propose a framework for global flood risk assessment for river floods, which can be applied in current conditions, as well as in future conditions due to climate and socio-economic changes. The framework's goal is to establish flood hazard and impact estimates at a high enough resolution to allow for their combination into a risk estimate, which can be used for strategic global flood risk assessments. The framework estimates hazard at a resolution of ~ 1 km2 using global forcing datasets of the current (or in scenario mode, future) climate, a global hydrological model, a global flood-routing model, and more importantly, an inundation downscaling routine. The second component of the framework combines hazard with flood impact models at the same resolution (e.g. damage, affected GDP, and affected population) to establish indicators for flood risk (e.g. annual expected damage, affected GDP, and affected population). The framework has been applied using the global hydrological model PCR-GLOBWB, which includes an optional global flood routing model DynRout, combined with scenarios from the Integrated Model to Assess the Global Environment (IMAGE). We performed downscaling of the hazard probability distributions to 1 km2 resolution with a new downscaling algorithm, applied on Bangladesh as a first case study application area. We demonstrate the risk assessment approach in Bangladesh based on GDP per capita data, population, and land use maps for 2010 and 2050. Validation of the hazard estimates has been performed using the Dartmouth Flood Observatory database. This was done by comparing a high return period flood with the maximum observed extent, as well as by comparing a time series of a single event with Dartmouth imagery of the event. Validation of modelled damage estimates was performed using observed damage estimates from the EM-DAT database and World Bank sources. We discuss and show sensitivities of the estimated risks with regard to the use of different climate input sets, decisions made in the downscaling algorithm, and different approaches to establish impact models.

Simulating PACE Global Ocean Radiances

NASA Technical Reports Server (NTRS)

Gregg, Watson W.; Rousseaux, Cecile S.

2017-01-01

The NASA PACE mission is a hyper-spectral radiometer planned for launch in the next decade. It is intended to provide new information on ocean biogeochemical constituents by parsing the details of high resolution spectral absorption and scattering. It is the first of its kind for global applications and as such, poses challenges for design and operation. To support pre-launch mission development and assess on-orbit capabilities, the NASA Global Modeling and Assimilation Office has developed a dynamic simulation of global water-leaving radiances, using an ocean model containing multiple ocean phytoplankton groups, particulate detritus, particulate inorganic carbon (PIC), and chromophoric dissolved organic carbon (CDOC) along with optical absorption and scattering processes at 1 nm spectral resolution. The purpose here is to assess the skill of the dynamic model and derived global radiances. Global bias, uncertainty, and correlation are derived using available modern satellite radiances at moderate spectral resolution. Total chlorophyll, PIC, and the absorption coefficient of CDOC (aCDOC), are simultaneously assimilated to improve the fidelity of the optical constituent fields. A 5-year simulation showed statistically significant (P < 0.05) comparisons of chlorophyll (r = 0.869), PIC (r = 0.868), and a CDOC (r =0.890) with satellite data. Additionally, diatoms (r = 0.890), cyanobacteria (r = 0.732), and coccolithophores (r = 0.716) were significantly correlated with in situ data. Global assimilated distributions of optical constituents were coupled with a radiative transfer model (Ocean-Atmosphere Spectral Irradiance Model, OASIM) to estimate normalized water-leaving radiances at 1 nm for the spectral range 250-800 nm. These unassimilated radiances were within 0.074 mW/sq cm/micron/sr of MODIS-Aqua radiances at 412, 443, 488, 531, 547, and 667 nm. This difference represented a bias of 10.4% (model low). A mean correlation of 0.706 (P < 0.05) was found with global distributions of MODIS radiances. These results suggest skill in the global assimilated model and resulting radiances. The reported error characterization suggests that the global dynamical simulation can support some aspects of mission design and analysis. For example, the high spectral resolution of the simulation supports investigations of band selection. The global nature of the radiance representations supports investigations of satellite observing scenarios. Global radiances at bands not available in current and past missions support investigations of mission capability. PACE, ocean color, water-leaving radiances, biogeochemical model, radiative transfer model

An Overview of Numerical Weather Prediction on Various Scales

NASA Astrophysics Data System (ADS)

Bao, J.-W.

2009-04-01

The increasing public need for detailed weather forecasts, along with the advances in computer technology, has motivated many research institutes and national weather forecasting centers to develop and run global as well as regional numerical weather prediction (NWP) models at high resolutions (i.e., with horizontal resolutions of ~10 km or higher for global models and 1 km or higher for regional models, and with ~60 vertical levels or higher). The need for running NWP models at high horizontal and vertical resolutions requires the implementation of non-hydrostatic dynamic core with a choice of horizontal grid configurations and vertical coordinates that are appropriate for high resolutions. Development of advanced numerics will also be needed for high resolution global and regional models, in particular, when the models are applied to transport problems and air quality applications. In addition to the challenges in numerics, the NWP community is also facing the challenges of developing physics parameterizations that are well suited for high-resolution NWP models. For example, when NWP models are run at resolutions of ~5 km or higher, the use of much more detailed microphysics parameterizations than those currently used in NWP model will become important. Another example is that regional NWP models at ~1 km or higher only partially resolve convective energy containing eddies in the lower troposphere. Parameterizations to account for the subgrid diffusion associated with unresolved turbulence still need to be developed. Further, physically sound parameterizations for air-sea interaction will be a critical component for tropical NWP models, particularly for hurricane predictions models. In this review presentation, the above issues will be elaborated on and the approaches to address them will be discussed.

Atmospheric Test Models and Numerical Experiments for the Simulation of the Global Distributions of Weather Data Transponders III. Horizontal Distributions

DOE Office of Scientific and Technical Information (OSTI.GOV)

Molenkamp, C.R.; Grossman, A.

1999-12-20

A network of small balloon-borne transponders which gather very high resolution wind and temperature data for use by modern numerical weather predication models has been proposed to improve the reliability of long-range weather forecasts. The global distribution of an array of such transponders is simulated using LLNL's atmospheric parcel transport model (GRANTOUR) with winds supplied by two different general circulation models. An initial study used winds from CCM3 with a horizontal resolution of about 3 degrees in latitude and longitude, and a second study used winds from NOGAPS with a 0.75 degree horizontal resolution. Results from both simulations show thatmore » reasonable global coverage can be attained by releasing balloons from an appropriate set of launch sites.« less

GFS Products

Science.gov Websites

Inventory Image of horizontal rule Global Products Updated: 7/28/2017 Global Forecast System (GFS) Model Global Data Assimilation System (GDAS) Model * Information about the GFS * Information about the GFS Name GFS GFS - Global longitude-latitude grid WCOSS File Name Inventory 0.25 degree resolution

P-CSI v1.0, an accelerated barotropic solver for the high-resolution ocean model component in the Community Earth System Model v2.0

NASA Astrophysics Data System (ADS)

Huang, Xiaomeng; Tang, Qiang; Tseng, Yuheng; Hu, Yong; Baker, Allison H.; Bryan, Frank O.; Dennis, John; Fu, Haohuan; Yang, Guangwen

2016-11-01

In the Community Earth System Model (CESM), the ocean model is computationally expensive for high-resolution grids and is often the least scalable component for high-resolution production experiments. The major bottleneck is that the barotropic solver scales poorly at high core counts. We design a new barotropic solver to accelerate the high-resolution ocean simulation. The novel solver adopts a Chebyshev-type iterative method to reduce the global communication cost in conjunction with an effective block preconditioner to further reduce the iterations. The algorithm and its computational complexity are theoretically analyzed and compared with other existing methods. We confirm the significant reduction of the global communication time with a competitive convergence rate using a series of idealized tests. Numerical experiments using the CESM 0.1° global ocean model show that the proposed approach results in a factor of 1.7 speed-up over the original method with no loss of accuracy, achieving 10.5 simulated years per wall-clock day on 16 875 cores.

Uncertainty of global summer precipitation in the CMIP5 models: a comparison between high-resolution and low-resolution models

NASA Astrophysics Data System (ADS)

Huang, Danqing; Yan, Peiwen; Zhu, Jian; Zhang, Yaocun; Kuang, Xueyuan; Cheng, Jing

2018-04-01

The uncertainty of global summer precipitation simulated by the 23 CMIP5 CGCMs and the possible impacts of model resolutions are investigated in this study. Large uncertainties exist over the tropical and subtropical regions, which can be mainly attributed to convective precipitation simulation. High-resolution models (HRMs) and low-resolution models (LRMs) are further investigated to demonstrate their different contributions to the uncertainties of the ensemble mean. It shows that the high-resolution model ensemble means (HMME) and low-resolution model ensemble mean (LMME) mitigate the biases between the MME and observation over most continents and oceans, respectively. The HMME simulates more precipitation than the LMME over most oceans, but less precipitation over some continents. The dominant precipitation category in the HRMs (LRMs) is the heavy precipitation (moderate precipitation) over the tropic regions. The combinations of convective and stratiform precipitation are also quite different: the HMME has much higher ratio of stratiform precipitation while the LMME has more convective precipitation. Finally, differences in precipitation between the HMME and LMME can be traced to their differences in the SST simulations via the local and remote air-sea interaction.

Observations and Modeling of Composition of Upper Troposphere/Lower Stratosphere (UTILS): Isentropic Mixing Events and Morphology of HNO3 as Observed by HIRDLS and Comparison with Results from Global Modeling Initiative

NASA Technical Reports Server (NTRS)

Rodriquez, J. M.; Douglass, A.R.; Yoshida, Y.; Strahan, S.; Duncan, B.; Olsen, M.; Gille, J.; Yudin, V.; Nardi, B.

2008-01-01

isentropic exchange of air masses between the tropical upper troposphere and mid-latitude lowermost stratosphere (the so-called "middle world") is an important pathway for stratospheric-tropospheric exchange. A seasonal, global view of this process has been difficult to obtain, in part due to the lack of the vertical resolution in satellite observations needed to capture the laminar character of these events. Ozone observations at a resolution of about 1 km from the High Resolution Dynamic Limb Sounder (HIRDLS) on NASA's Aura satellite show instances of these intrusions. Such intrusions should also be observable in HN03 observations; however, the abundances of nitric acid could be additionally controlled by chemical processes or incorporation and removal into ice clouds. We present a systematic examination of the HIRDLS data on O3 and HNO3 to determine the seasonal and spatial characteristics of the distribution of isentropic intrusions. At the same time, we compare the observed distributions with those calculated by the Global Modeling Initiative combined tropospheric-stratospheric model, which has a vertical resolution of about I km. This Chemical Transport Model (CTM) is driven by meteorological fields obtained from the GEOS-4 system of NASA/Goddard Global Modeling and Assimilation Office (GMAO), for the Aura time period, at a vertical resolution of about 1 km. Such comparison brings out the successes and limitations of the model in representing isentropic stratospheric-tropospheric exchange, and the different processes controlling HNO3 in the UTAS.

High resolution global flood hazard map from physically-based hydrologic and hydraulic models.

NASA Astrophysics Data System (ADS)

Begnudelli, L.; Kaheil, Y.; McCollum, J.

2017-12-01

The global flood map published online at http://www.fmglobal.com/research-and-resources/global-flood-map at 90m resolution is being used worldwide to understand flood risk exposure, exercise certain measures of mitigation, and/or transfer the residual risk financially through flood insurance programs. The modeling system is based on a physically-based hydrologic model to simulate river discharges, and 2D shallow-water hydrodynamic model to simulate inundation. The model can be applied to large-scale flood hazard mapping thanks to several solutions that maximize its efficiency and the use of parallel computing. The hydrologic component of the modeling system is the Hillslope River Routing (HRR) hydrologic model. HRR simulates hydrological processes using a Green-Ampt parameterization, and is calibrated against observed discharge data from several publicly-available datasets. For inundation mapping, we use a 2D Finite-Volume Shallow-Water model with wetting/drying. We introduce here a grid Up-Scaling Technique (UST) for hydraulic modeling to perform simulations at higher resolution at global scale with relatively short computational times. A 30m SRTM is now available worldwide along with higher accuracy and/or resolution local Digital Elevation Models (DEMs) in many countries and regions. UST consists of aggregating computational cells, thus forming a coarser grid, while retaining the topographic information from the original full-resolution mesh. The full-resolution topography is used for building relationships between volume and free surface elevation inside cells and computing inter-cell fluxes. This approach almost achieves computational speed typical of the coarse grids while preserving, to a significant extent, the accuracy offered by the much higher resolution available DEM. The simulations are carried out along each river of the network by forcing the hydraulic model with the streamflow hydrographs generated by HRR. Hydrographs are scaled so that the peak corresponds to the return period corresponding to the hazard map being produced (e.g. 100 years, 500 years). Each numerical simulation models one river reach, except for the longest reaches which are split in smaller parts. Here we show results for selected river basins worldwide.

Re-assessing Present Day Global Mass Transport and Glacial Isostatic Adjustment From a Data Driven Approach

NASA Astrophysics Data System (ADS)

Wu, X.; Jiang, Y.; Simonsen, S.; van den Broeke, M. R.; Ligtenberg, S.; Kuipers Munneke, P.; van der Wal, W.; Vermeersen, B. L. A.

2017-12-01

Determining present-day mass transport (PDMT) is complicated by the fact that most observations contain signals from both present day ice melting and Glacial Isostatic Adjustment (GIA). Despite decades of progress in geodynamic modeling and new observations, significant uncertainties remain in both. The key to separate present-day ice mass change and signals from GIA is to include data of different physical characteristics. We designed an approach to separate PDMT and GIA signatures by estimating them simultaneously using globally distributed interdisciplinary data with distinct physical information and a dynamically constructed a priori GIA model. We conducted a high-resolution global reappraisal of present-day ice mass balance with focus on Earth's polar regions and its contribution to global sea-level rise using a combination of ICESat, GRACE gravity, surface geodetic velocity data, and an ocean bottom pressure model. Adding ice altimetry supplies critically needed dual data types over the interiors of ice covered regions to enhance separation of PDMT and GIA signatures, and achieve half an order of magnitude expected higher accuracies for GIA and consequently ice mass balance estimates. The global data based approach can adequately address issues of PDMT and GIA induced geocenter motion and long-wavelength signatures important for large areas such as Antarctica and global mean sea level. In conjunction with the dense altimetry data, we solved for PDMT coefficients up to degree and order 180 by using a higher-resolution GRACE data set, and a high-resolution a priori PDMT model that includes detailed geographic boundaries. The high-resolution approach solves the problem of multiple resolutions in various data types, greatly reduces aliased errors from a low-degree truncation, and at the same time, enhances separation of signatures from adjacent regions such as Greenland and Canadian Arctic territories.

Validation of a Global Hydrodynamic Flood Inundation Model

NASA Astrophysics Data System (ADS)

Bates, P. D.; Smith, A.; Sampson, C. C.; Alfieri, L.; Neal, J. C.

2014-12-01

In this work we present first validation results for a hyper-resolution global flood inundation model. We use a true hydrodynamic model (LISFLOOD-FP) to simulate flood inundation at 1km resolution globally and then use downscaling algorithms to determine flood extent and depth at 90m spatial resolution. Terrain data are taken from a custom version of the SRTM data set that has been processed specifically for hydrodynamic modelling. Return periods of flood flows along the entire global river network are determined using: (1) empirical relationships between catchment characteristics and index flood magnitude in different hydroclimatic zones derived from global runoff data; and (2) an index flood growth curve, also empirically derived. Bankful return period flow is then used to set channel width and depth, and flood defence impacts are modelled using empirical relationships between GDP, urbanization and defence standard of protection. The results of these simulations are global flood hazard maps for a number of different return period events from 1 in 5 to 1 in 1000 years. We compare these predictions to flood hazard maps developed by national government agencies in the UK and Germany using similar methods but employing detailed local data, and to observed flood extent at a number of sites including St. Louis, USA and Bangkok in Thailand. Results show that global flood hazard models can have considerable skill given careful treatment to overcome errors in the publicly available data that are used as their input.

Effects of Model Resolution and Ocean Mixing on Forced Ice-Ocean Physical and Biogeochemical Simulations Using Global and Regional System Models

NASA Astrophysics Data System (ADS)

Jin, Meibing; Deal, Clara; Maslowski, Wieslaw; Matrai, Patricia; Roberts, Andrew; Osinski, Robert; Lee, Younjoo J.; Frants, Marina; Elliott, Scott; Jeffery, Nicole; Hunke, Elizabeth; Wang, Shanlin

2018-01-01

The current coarse-resolution global Community Earth System Model (CESM) can reproduce major and large-scale patterns but is still missing some key biogeochemical features in the Arctic Ocean, e.g., low surface nutrients in the Canada Basin. We incorporated the CESM Version 1 ocean biogeochemical code into the Regional Arctic System Model (RASM) and coupled it with a sea-ice algal module to investigate model limitations. Four ice-ocean hindcast cases are compared with various observations: two in a global 1° (40˜60 km in the Arctic) grid: G1deg and G1deg-OLD with/without new sea-ice processes incorporated; two on RASM's 1/12° (˜9 km) grid R9km and R9km-NB with/without a subgrid scale brine rejection parameterization which improves ocean vertical mixing under sea ice. Higher-resolution and new sea-ice processes contributed to lower model errors in sea-ice extent, ice thickness, and ice algae. In the Bering Sea shelf, only higher resolution contributed to lower model errors in salinity, nitrate (NO3), and chlorophyll-a (Chl-a). In the Arctic Basin, model errors in mixed layer depth (MLD) were reduced 36% by brine rejection parameterization, 20% by new sea-ice processes, and 6% by higher resolution. The NO3 concentration biases were caused by both MLD bias and coarse resolution, because of excessive horizontal mixing of high NO3 from the Chukchi Sea into the Canada Basin in coarse resolution models. R9km showed improvements over G1deg on NO3, but not on Chl-a, likely due to light limitation under snow and ice cover in the Arctic Basin.

Classification of Clouds and Deep Convection from GEOS-5 Using Satellite Observations

NASA Technical Reports Server (NTRS)

Putman, William; Suarez, Max

2010-01-01

With the increased resolution of global atmospheric models and the push toward global cloud resolving models, the resemblance of model output to satellite observations has become strikingly similar. As we progress with our adaptation of the Goddard Earth Observing System Model, Version 5 (GEOS-5) as a high resolution cloud system resolving model, evaluation of cloud properties and deep convection require in-depth analysis beyond a visual comparison. Outgoing long-wave radiation (OLR) provides a sufficient comparison with infrared (IR) satellite imagery to isolate areas of deep convection. We have adopted a binning technique to generate a series of histograms for OLR which classify the presence and fraction of clear sky versus deep convection in the tropics that can be compared with a similar analyses of IR imagery from composite Geostationary Operational Environmental Satellite (GOES) observations. We will present initial results that have been used to evaluate the amount of deep convective parameterization required within the model as we move toward cloud system resolving resolutions of 10- to 1-km globally.

Simulating PACE Global Ocean Radiances

PubMed Central

Gregg, Watson W.; Rousseaux, Cécile S.

2017-01-01

The NASA PACE mission is a hyper-spectral radiometer planned for launch in the next decade. It is intended to provide new information on ocean biogeochemical constituents by parsing the details of high resolution spectral absorption and scattering. It is the first of its kind for global applications and as such, poses challenges for design and operation. To support pre-launch mission development and assess on-orbit capabilities, the NASA Global Modeling and Assimilation Office has developed a dynamic simulation of global water-leaving radiances, using an ocean model containing multiple ocean phytoplankton groups, particulate detritus, particulate inorganic carbon (PIC), and chromophoric dissolved organic carbon (CDOC) along with optical absorption and scattering processes at 1 nm spectral resolution. The purpose here is to assess the skill of the dynamic model and derived global radiances. Global bias, uncertainty, and correlation are derived using available modern satellite radiances at moderate spectral resolution. Total chlorophyll, PIC, and the absorption coefficient of CDOC (aCDOC), are simultaneously assimilated to improve the fidelity of the optical constituent fields. A 5-year simulation showed statistically significant (P <0.05) comparisons of chlorophyll (r = 0.869), PIC (r = 0.868), and aCDOC (r = 0.890) with satellite data. Additionally, diatoms (r = 0.890), cyanobacteria (r = 0.732), and coccolithophores (r = 0.716) were significantly correlated with in situ data. Global assimilated distributions of optical constituents were coupled with a radiative transfer model (Ocean-Atmosphere Spectral Irradiance Model, OASIM) to estimate normalized water-leaving radiances at 1 nm for the spectral range 250–800 nm. These unassimilated radiances were within −0.074 mW cm−2 μm1 sr−1 of MODIS-Aqua radiances at 412, 443, 488, 531, 547, and 667 nm. This difference represented a bias of −10.4% (model low). A mean correlation of 0.706 (P < 0.05) was found with global distributions of MODIS radiances. These results suggest skill in the global assimilated model and resulting radiances. The reported error characterization suggests that the global dynamical simulation can support some aspects of mission design and analysis. For example, the high spectral resolution of the simulation supports investigations of band selection. The global nature of the radiance representations supports investigations of satellite observing scenarios. Global radiances at bands not available in current and past missions support investigations of mission capability. PMID:29292403

Vorticity-divergence semi-Lagrangian global atmospheric model SL-AV20: dynamical core

NASA Astrophysics Data System (ADS)

Tolstykh, Mikhail; Shashkin, Vladimir; Fadeev, Rostislav; Goyman, Gordey

2017-05-01

SL-AV (semi-Lagrangian, based on the absolute vorticity equation) is a global hydrostatic atmospheric model. Its latest version, SL-AV20, provides global operational medium-range weather forecast with 20 km resolution over Russia. The lower-resolution configurations of SL-AV20 are being tested for seasonal prediction and climate modeling. The article presents the model dynamical core. Its main features are a vorticity-divergence formulation at the unstaggered grid, high-order finite-difference approximations, semi-Lagrangian semi-implicit discretization and the reduced latitude-longitude grid with variable resolution in latitude. The accuracy of SL-AV20 numerical solutions using a reduced lat-lon grid and the variable resolution in latitude is tested with two idealized test cases. Accuracy and stability of SL-AV20 in the presence of the orography forcing are tested using the mountain-induced Rossby wave test case. The results of all three tests are in good agreement with other published model solutions. It is shown that the use of the reduced grid does not significantly affect the accuracy up to the 25 % reduction in the number of grid points with respect to the regular grid. Variable resolution in latitude allows us to improve the accuracy of a solution in the region of interest.

A high-resolution global flood hazard model

NASA Astrophysics Data System (ADS)

Sampson, Christopher C.; Smith, Andrew M.; Bates, Paul B.; Neal, Jeffrey C.; Alfieri, Lorenzo; Freer, Jim E.

2015-09-01

Floods are a natural hazard that affect communities worldwide, but to date the vast majority of flood hazard research and mapping has been undertaken by wealthy developed nations. As populations and economies have grown across the developing world, so too has demand from governments, businesses, and NGOs for modeled flood hazard data in these data-scarce regions. We identify six key challenges faced when developing a flood hazard model that can be applied globally and present a framework methodology that leverages recent cross-disciplinary advances to tackle each challenge. The model produces return period flood hazard maps at ˜90 m resolution for the whole terrestrial land surface between 56°S and 60°N, and results are validated against high-resolution government flood hazard data sets from the UK and Canada. The global model is shown to capture between two thirds and three quarters of the area determined to be at risk in the benchmark data without generating excessive false positive predictions. When aggregated to ˜1 km, mean absolute error in flooded fraction falls to ˜5%. The full complexity global model contains an automatically parameterized subgrid channel network, and comparison to both a simplified 2-D only variant and an independently developed pan-European model shows the explicit inclusion of channels to be a critical contributor to improved model performance. While careful processing of existing global terrain data sets enables reasonable model performance in urban areas, adoption of forthcoming next-generation global terrain data sets will offer the best prospect for a step-change improvement in model performance.

Enhancing SMAP Soil Moisture Retrievals via Superresolution Techniques

NASA Astrophysics Data System (ADS)

Beale, K. D.; Ebtehaj, A. M.; Romberg, J. K.; Bras, R. L.

2017-12-01

Soil moisture is a key state variable that modulates land-atmosphere interactions and its high-resolution global scale estimates are essential for improved weather forecasting, drought prediction, crop management, and the safety of troop mobility. Currently, NASA's Soil Moisture Active/Passive (SMAP) satellite provides a global picture of soil moisture variability at a resolution of 36 km, which is prohibitive for some hydrologic applications. The goal of this research is to enhance the resolution of SMAP passive microwave retrievals by a factor of 2 to 4 using modern superresolution techniques that rely on the knowledge of high-resolution land surface models. In this work, we explore several super-resolution techniques including an empirical dictionary method, a learned dictionary method, and a three-layer convolutional neural network. Using a year of global high-resolution land surface model simulations as training set, we found that we are able to produce high-resolution soil moisture maps that outperform the original low-resolution observations both qualitatively and quantitatively. In particular, on a patch-by-patch basis we are able to produce estimates of high-resolution soil moisture maps that improve on the original low-resolution patches by on average 6% in terms of mean-squared error, and 14% in terms of the structural similarity index.

Impact of recent Global Digital Bathymetry and Topography Models on geoid modelling: Results from two case studies in Balearic and Aegean Seas

NASA Astrophysics Data System (ADS)

Delikaraoglou, D.; Mintourakis, I.; Kallianou, F.

2009-04-01

With the realization of the Shuttle Radar Topographic Mission (SRTM) and the free distribution of its global elevation dataset with 3 arcsec (90 m) resolution and less than 16 m vertical accuracy, together with the availability of the higher resolution (30 m) and accuracy (10 m) Digital Terrain Models (DTM) from the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER), these two valuable sources of uniform DEM data represent a revolution in the world of terrain modelling. DEMs are an important source of data for the generation of high resolution geoids since they provide the high-frequency content of the gravity field spectrum and are suitable for the computation of terrain effects to gravity and indirect effects to the geoid, thus allowing the combination of global geopotential models, local gravity anomalies and information about the earth's topography (represented by a given DEM). However, although such models are available for land, there are no readily accessible Digital Bathymetry Models (DBMs) of equivalent quality for the coastal and oceanic regions. Most of the global DBM's (e.g. ETOPO1, SRTM30, and GEBCO global bathymetric grid) are compilations of heterogeneous data with medium resolution and accuracy. This prevents to exploit the potential of the recent high resolution (1 arcmin) marine free-air gravity anomalies datasets derived from satellite altimetry (such as the DNSC08, and the Sandwell & Smith v18.1 (S&Sv18.1) global solutions) in conjunction with such global DBM's. Fortunately, for some regions, recently have become available DBM's of much better accuracy and resolution, such as the DBM of 1 km resolution for many regions of the Mediterranean Sea which is distributed by IFREMER, the French Research Institute for Exploitation of the Sea. The scope of this study is to use this latest regional DBM in combination with the newly available DNSC08 and SSV18.1 global marine free-air gravity anomalies datasets for marine and near shore geoid modelling of archipelagic (island) areas. We have concentrated in two test regions: (a) the Catalano-Balearic Sea (South of Spain in the NW Meditteranean), where adequate marine and land gravity data allow a detailed evaluation of our processing methodologies and their results and, (b) the Aegean Sea where the presence of many islands in varying distances from the mainland Greece and located on the continental shelf and/or divided by steep sea floor topography present some unique challenges for any high resolution geoid modelling efforts. For both test regions, we generated a combined DEM (C-DEM) using the IFREMER and SRTM 30 arcsec bathymetric data for the sea areas and SRTM 3 arcsec data for the surrounding land areas. In this contribution, we discuss various computational aspects relating to the so-called "Direct Topographical Effect" (DTE) and the "Indirect Topographical Effect" (ITE), the two most significant topographical effects that have to be evaluated when a precise geoid is being compiled. In addition, we outline the evaluation and the impact of the results obtained, especially with regard to the differences in the geoid models when different elevation data are used, and point out the main limitations and possibilities for further improvements in the use of the aforementioned satellite and terrestrial data for regional and local geoid mapping in coastal and island regions. Keywords: IFREMER, SRTM, terrain effects, free-air gravity anomalies, geoid modelling,Digital Bathymetry Models.

Evaluation of GOCE-based Global Geoid Models in Finnish Territory

NASA Astrophysics Data System (ADS)

Saari, Timo; Bilker-Koivula, Mirjam

2015-04-01

The gravity satellite mission GOCE made its final observations in the fall of 2013. By then it had exceeded its expected lifespan of one year with more than three additional years. Thus, the mission collected more data from the Earth's gravitational field than expected, and more comprehensive global geoid models have been derived ever since. The GOCE High-level Processing Facility (HPF) by ESA has published GOCE global gravity field models annually. We compared all of the 12 HPF-models as well as 3 additional GOCE, 11 GRACE and 6 combined GOCE+GRACE models with GPS-levelling data and gravity observations in Finland. The most accurate models were compared against high resolution global geoid models EGM96 and EGM2008. The models were evaluated up to three different degrees and order: 150 (the common maximum for the GRACE models), 240 (the common maximum for the GOCE models) and maximum. When coefficients up to degree and order 150 are used, the results of the GOCE models are comparable with the results of the latest GRACE models. Generally, all of the latest GOCE and GOCE+GRACE models give standard deviations of the height anomaly differences of around 15 cm and of gravity anomaly differences of around 10 mgal over Finland. The best solutions were not always achieved with the highest maximum degree and order of the satellite gravity field models, since the highest coefficients (above 240) may be less accurately determined. Over Finland, the latest GOCE and GOCE+GRACE models give similar results as the high resolution models EGM96 and EGM2008 when coefficients up to degree and order 240 are used. This is mainly due to the high resolution terrestrial data available in the area of Finland, which was used in the high resolution models.

Downscaling global precipitation for local applications - a case for the Rhine basin

NASA Astrophysics Data System (ADS)

Sperna Weiland, Frederiek; van Verseveld, Willem; Schellekens, Jaap

2017-04-01

Within the EU FP7 project eartH2Observe a global Water Resources Re-analysis (WRR) is being developed. This re-analysis consists of meteorological and hydrological water balance variables with global coverage, spanning the period 1979-2014 at 0.25 degrees resolution (Schellekens et al., 2016). The dataset can be of special interest in regions with limited in-situ data availability, yet for local scale analysis particularly in mountainous regions, a resolution of 0.25 degrees may be too coarse and downscaling the data to a higher resolution may be required. A downscaling toolbox has been made that includes spatial downscaling of precipitation based on the global WorldClim dataset that is available at 1 km resolution as a monthly climatology (Hijmans et al., 2005). The input of the down-scaling tool are either the global eartH2Observe WRR1 and WRR2 datasets based on the WFDEI correction methodology (Weedon et al., 2014) or the global Multi-Source Weighted-Ensemble Precipitation (MSWEP) dataset (Beck et al., 2016). Here we present a validation of the datasets over the Rhine catchment by means of a distributed hydrological model (wflow, Schellekens et al., 2014) using a number of precipitation scenarios. (1) We start by running the model using the local reference dataset derived by spatial interpolation of gauge observations. Furthermore we use (2) the MSWEP dataset at the native 0.25-degree resolution followed by (3) MSWEP downscaled with the WorldClim dataset and final (4) MSWEP downscaled with the local reference dataset. The validation will be based on comparison of the modeled river discharges as well as rainfall statistics. We expect that down-scaling the MSWEP dataset with the WorldClim data to higher resolution will increase its performance. To test the performance of the down-scaling routine we have added a run with MSWEP data down-scaled with the local dataset and compare this with the run based on the local dataset itself. - Beck, H. E. et al., 2016. MSWEP: 3-hourly 0.25° global gridded precipitation (1979-2015) by merging gauge, satellite, and reanalysis data, Hydrol. Earth Syst. Sci. Discuss., doi:10.5194/hess-2016-236, accepted for final publication. - Hijmans, R.J. et al., 2005. Very high resolution interpolated climate surfaces for global land areas. International Journal of Climatology 25: 1965-1978. - Schellekens, J. et al., 2016. A global water resources ensemble of hydrological models: the eartH2Observe Tier-1 dataset, Earth Syst. Sci. Data Discuss., doi:10.5194/essd-2016-55, under review. - Schellekens, J. et al., 2014. Rapid setup of hydrological and hydraulic models using OpenStreetMap and the SRTM derived digital elevation model. Environmental Modelling&Software - Weedon, G.P. et al., 2014. The WFDEI meteorological forcing data set: WATCH Forcing Data methodology applied to ERA-Interim reanalysis data. Water Resources Research, 50, doi:10.1002/2014WR015638.

Regional Community Climate Simulations with variable resolution meshes in the Community Earth System Model

NASA Astrophysics Data System (ADS)

Zarzycki, C. M.; Gettelman, A.; Callaghan, P.

2017-12-01

Accurately predicting weather extremes such as precipitation (floods and droughts) and temperature (heat waves) requires high resolution to resolve mesoscale dynamics and topography at horizontal scales of 10-30km. Simulating such resolutions globally for climate scales (years to decades) remains computationally impractical. Simulating only a small region of the planet is more tractable at these scales for climate applications. This work describes global simulations using variable-resolution static meshes with multiple dynamical cores that target the continental United States using developmental versions of the Community Earth System Model version 2 (CESM2). CESM2 is tested in idealized, aquaplanet and full physics configurations to evaluate variable mesh simulations against uniform high and uniform low resolution simulations at resolutions down to 15km. Different physical parameterization suites are also evaluated to gauge their sensitivity to resolution. Idealized variable-resolution mesh cases compare well to high resolution tests. More recent versions of the atmospheric physics, including cloud schemes for CESM2, are more stable with respect to changes in horizontal resolution. Most of the sensitivity is due to sensitivity to timestep and interactions between deep convection and large scale condensation, expected from the closure methods. The resulting full physics model produces a comparable climate to the global low resolution mesh and similar high frequency statistics in the high resolution region. Some biases are reduced (orographic precipitation in the western United States), but biases do not necessarily go away at high resolution (e.g. summertime JJA surface Temp). The simulations are able to reproduce uniform high resolution results, making them an effective tool for regional climate studies and are available in CESM2.

A high‐resolution global flood hazard model†

PubMed Central

Smith, Andrew M.; Bates, Paul D.; Neal, Jeffrey C.; Alfieri, Lorenzo; Freer, Jim E.

2015-01-01

Abstract Floods are a natural hazard that affect communities worldwide, but to date the vast majority of flood hazard research and mapping has been undertaken by wealthy developed nations. As populations and economies have grown across the developing world, so too has demand from governments, businesses, and NGOs for modeled flood hazard data in these data‐scarce regions. We identify six key challenges faced when developing a flood hazard model that can be applied globally and present a framework methodology that leverages recent cross‐disciplinary advances to tackle each challenge. The model produces return period flood hazard maps at ∼90 m resolution for the whole terrestrial land surface between 56°S and 60°N, and results are validated against high‐resolution government flood hazard data sets from the UK and Canada. The global model is shown to capture between two thirds and three quarters of the area determined to be at risk in the benchmark data without generating excessive false positive predictions. When aggregated to ∼1 km, mean absolute error in flooded fraction falls to ∼5%. The full complexity global model contains an automatically parameterized subgrid channel network, and comparison to both a simplified 2‐D only variant and an independently developed pan‐European model shows the explicit inclusion of channels to be a critical contributor to improved model performance. While careful processing of existing global terrain data sets enables reasonable model performance in urban areas, adoption of forthcoming next‐generation global terrain data sets will offer the best prospect for a step‐change improvement in model performance. PMID:27594719

Global terrain classification using Multiple-Error-Removed Improved-Terrain (MERIT) to address susceptibility of landslides and other geohazards

NASA Astrophysics Data System (ADS)

Iwahashi, J.; Yamazaki, D.; Matsuoka, M.; Thamarux, P.; Herrick, J.; Yong, A.; Mital, U.

2017-12-01

A seamless model of landform classifications with regional accuracy will be a powerful platform for geophysical studies that forecast geologic hazards. Spatial variability as a function of landform on a global scale was captured in the automated classifications of Iwahashi and Pike (2007) and additional developments are presented here that incorporate more accurate depictions using higher-resolution elevation data than the original 1-km scale Shuttle Radar Topography Mission digital elevation model (DEM). We create polygon-based terrain classifications globally by using the 280-m DEM interpolated from the Multi-Error-Removed Improved-Terrain DEM (MERIT; Yamazaki et al., 2017). The multi-scale pixel-image analysis method, known as Multi-resolution Segmentation (Baatz and Schäpe, 2000), is first used to classify the terrains based on geometric signatures (slope and local convexity) calculated from the 280-m DEM. Next, we apply the machine learning method of "k-means clustering" to prepare the polygon-based classification at the globe-scale using slope, local convexity and surface texture. We then group the divisions with similar properties by hierarchical clustering and other statistical analyses using geological and geomorphological data of the area where landslides and earthquakes are frequent (e.g. Japan and California). We find the 280-m DEM resolution is only partially sufficient for classifying plains. We nevertheless observe that the categories correspond to reported landslide and liquefaction features at the global scale, suggesting that our model is an appropriate platform to forecast ground failure. To predict seismic amplification, we estimate site conditions using the time-averaged shear-wave velocity in the upper 30-m (VS30) measurements compiled by Yong et al. (2016) and the terrain model developed by Yong (2016; Y16). We plan to test our method on finer resolution DEMs and report our findings to obtain a more globally consistent terrain model as there are known errors in DEM derivatives at higher-resolutions. We expect the improvement in DEM resolution (4 times greater detail) and the combination of regional and global coverage will yield a consistent dataset of polygons that have the potential to improve relations to the Y16 estimates significantly.

High Resolution Land Surface Modeling with the next generation Land Data Assimilation Systems

NASA Astrophysics Data System (ADS)

Kumar, S. V.; Eylander, J.; Peters-Lidard, C.

2005-12-01

Knowledge of land surface processes is important to many real-world applications such as agricultural production, water resources management, and flood predication. The Air Force Weather Agency (AFWA) has provided the USDA and other customers global soil moisture and temperature data for the past 30 years using the agrometeorological data assimilation model (now called AGRMET), merging atmospheric data. Further, accurate initialization of land surface conditions has been shown to greatly influence and improve weather forecast model and seasonal-to-interannual climate predictions. The AFWA AGRMET model exploits real time precipitation observations and analyses, global forecast model and satellite data to generate global estimates of soil moisture, soil temperature and other land surface states at 48km spatial resolution. However, to truly address the land surface initialization and climate prediction problem, and to mitigate the errors introduced by the differences in spatial scales of models, representations of land surface conditions need to be developed at the same fine scales such as that of cloud resolving models. NASA's Goddard Space Flight Center has developed an offline land data assimilation system known as the Land Information System (LIS) capable of modeling land atmosphere interactions at spatial resolutions as fine as 1km. LIS provides a software architecture that integrates the use of the state of the art land surface models, data assimilation techniques, and high performance computing and data management tools. LIS also employs many high resolution surface parameters such as the NASA Earth Observing System (EOS)-era products. In this study we describe the development of a next generation high resolution land surface modeling and data assimilation system, combining the capabilities of LIS and AGRMET. We investigate the influence of high resolution land surface data and observations on the land surface conditions by comparing with the operational AGRMET outputs.

Coarsening of physics for biogeochemical model in NEMO

NASA Astrophysics Data System (ADS)

Bricaud, Clement; Le Sommer, Julien; Madec, Gurvan; Deshayes, Julie; Chanut, Jerome; Perruche, Coralie

2017-04-01

Ocean mesoscale and submesoscale turbulence contribute to ocean tracer transport and to shaping ocean biogeochemical tracers distribution. Representing adequately tracer transport in ocean models therefore requires to increase model resolution so that the impact of ocean turbulence is adequately accounted for. But due to supercomputers power and storage limitations, global biogeochemical models are not yet run routinely at eddying resolution. Still, because the "effective resolution" of eddying ocean models is much coarser than the physical model grid resolution, tracer transport can be reconstructed to a large extent by computing tracer transport and diffusion with a model grid resolution close to the effective resolution of the physical model. This observation has motivated the implementation of a new capability in NEMO ocean model (http://www.nemo-ocean.eu/) that allows to run the physical model and the tracer transport model at different grid resolutions. In a first time, we present results obtained with this new capability applied to a synthetic age tracer in a global eddying model configuration. In this model configuration, ocean dynamic is computed at ¼° resolution but tracer transport is computed at 3/4° resolution. The solution obtained is compared to 2 reference setup ,one at ¼° resolution for both physics and passive tracer models and one at 3/4° resolution for both physics and passive tracer model. We discuss possible options for defining the vertical diffusivity coefficient for the tracer transport model based on information from the high resolution grid. We describe the impact of this choice on the distribution and one the penetration of the age tracer. In a second time we present results obtained by coupling the physics with the biogeochemical model PISCES. We look at the impact of this methodology on some tracers distribution and dynamic. The method described here can found applications in ocean forecasting, such as the Copernicus Marine service operated by Mercator-Ocean, and in Earth System Models for climate applications.

Nonhydrostatic icosahedral atmospheric model (NICAM) for global cloud resolving simulations

NASA Astrophysics Data System (ADS)

Satoh, M.; Matsuno, T.; Tomita, H.; Miura, H.; Nasuno, T.; Iga, S.

2008-03-01

A new type of ultra-high resolution atmospheric global circulation model is developed. The new model is designed to perform "cloud resolving simulations" by directly calculating deep convection and meso-scale circulations, which play key roles not only in the tropical circulations but in the global circulations of the atmosphere. Since cores of deep convection have a few km in horizontal size, they have not directly been resolved by existing atmospheric general circulation models (AGCMs). In order to drastically enhance horizontal resolution, a new framework of a global atmospheric model is required; we adopted nonhydrostatic governing equations and icosahedral grids to the new model, and call it Nonhydrostatic ICosahedral Atmospheric Model (NICAM). In this article, we review governing equations and numerical techniques employed, and present the results from the unique 3.5-km mesh global experiments—with O(10 9) computational nodes—using realistic topography and land/ocean surface thermal forcing. The results show realistic behaviors of multi-scale convective systems in the tropics, which have not been captured by AGCMs. We also argue future perspective of the roles of the new model in the next generation atmospheric sciences.

Clouds in ECMWF's 30 KM Resolution Global Atmospheric Forecast Model (TL639)

NASA Technical Reports Server (NTRS)

Cahalan, R. F.; Morcrette, J. J.

1999-01-01

Global models of the general circulation of the atmosphere resolve a wide range of length scales, and in particular cloud structures extend from planetary scales to the smallest scales resolvable, now down to 30 km in state-of-the-art models. Even the highest resolution models do not resolve small-scale cloud phenomena seen, for example, in Landsat and other high-resolution satellite images of clouds. Unresolved small-scale disturbances often grow into larger ones through non-linear processes that transfer energy upscale. Understanding upscale cascades is of crucial importance in predicting current weather, and in parameterizing cloud-radiative processes that control long term climate. Several movie animations provide examples of the temporal and spatial variation of cloud fields produced in 4-day runs of the forecast model at the European Centre for Medium-Range Weather Forecasts (ECMWF) in Reading, England, at particular times and locations of simultaneous measurement field campaigns. model resolution is approximately 30 km horizontally (triangular truncation TL639) with 31 vertical levels from surface to stratosphere. Timestep of the model is about 10 minutes, but animation frames are 3 hours apart, at timesteps when the radiation is computed. The animations were prepared from an archive of several 4-day runs at the highest available model resolution, and archived at ECMWF. Cloud, wind and temperature fields in an approximately 1000 km X 1000 km box were retrieved from the archive, then approximately 60 Mb Vis5d files were prepared with the help of Graeme Kelly of ECMWF, and were compressed into MPEG files each less than 3 Mb. We discuss the interaction of clouds and radiation in the model, and compare the variability of cloud liquid as a function of scale to that seen in cloud observations made in intensive field campaigns. Comparison of high-resolution global runs to cloud-resolving models, and to lower resolution climate models is leading to better understanding of the upscale cascade and suggesting new cloud-radiation parameterizations for climate models.

Orographic precipitation at global and regional scales: Observational uncertainty and evaluation of 25-km global model simulations

NASA Astrophysics Data System (ADS)

Schiemann, Reinhard; Roberts, Charles J.; Bush, Stephanie; Demory, Marie-Estelle; Strachan, Jane; Vidale, Pier Luigi; Mizielinski, Matthew S.; Roberts, Malcolm J.

2015-04-01

Precipitation over land exhibits a high degree of variability due to the complex interaction of the precipitation generating atmospheric processes with coastlines, the heterogeneous land surface, and orography. Global general circulation models (GCMs) have traditionally had very limited ability to capture this variability on the mesoscale (here ~50-500 km) due to their low resolution. This has changed with recent investments in resolution and ensembles of multidecadal climate simulations of atmospheric GCMs (AGCMs) with ~25 km grid spacing are becoming increasingly available. Here, we evaluate the mesoscale precipitation distribution in one such set of simulations obtained in the UPSCALE (UK on PrACE - weather-resolving Simulations of Climate for globAL Environmental risk) modelling campaign with the HadGEM-GA3 AGCM. Increased model resolution also poses new challenges to the observational datasets used to evaluate models. Global gridded data products such as those provided by the Global Precipitation Climatology Project (GPCP) are invaluable for assessing large-scale features of the precipitation distribution but may not sufficiently resolve mesoscale structures. In the absence of independent estimates, the intercomparison of different observational datasets may be the only way to get some insight into the uncertainties associated with these observations. Here, we focus on mid-latitude continental regions where observations based on higher-density gauge networks are available in addition to the global data sets: Europe/the Alps, South and East Asia, and the continental US. The ability of GCMs to represent mesoscale variability is of interest in its own right, as climate information on this scale is required by impact studies. An additional motivation for the research proposed here arises from continuing efforts to quantify the components of the global radiation budget and water cycle. Recent estimates based on radiation measurements suggest that the global mean precipitation/evaporation may be up to 10 Wm-2 (about 0.35 mm day-1) larger than the estimate obtained from GPCP. While the main part of this discrepancy is thought to be due to the underestimation of remotely-sensed ocean precipitation, there is also considerable uncertainty about 'unobserved' precipitation over land, in particular in the form of snow in regions of high latitude/altitude. We aim to contribute to this discussion, at least at a qualitative level, by considering case studies of how area-averaged mountain precipitation is represented in different observational datasets and by HadGEM3-GA3 at different resolutions. Our results show that the AGCM simulates considerably more orographic precipitation at higher resolution. We find this at the global scale both for the winter and summer hemispheres, as well as in several case studies in mid-latitude regions. Gridded observations based on gauge measurements generally capture the mesoscale spatial variability of precipitation, but differ strongly from one another in the magnitude of area-averaged precipitation, so that they are of very limited use for evaluating this aspect of the modelled climate. We are currently conducting a sensitivity experiment (coarse-grained orography in high-resolution HadGEM3) to further investigate the resolution sensitivity seen in the model.

The Resolution Sensitivity of Northern Hemisphere Blocking in Four 25-km Atmospheric Global Circulation Models

DOE Office of Scientific and Technical Information (OSTI.GOV)

Schiemann, Reinhard; Demory, Marie-Estelle; Shaffrey, Len C.

The aim of this study is to investigate if the representation of Northern Hemisphere blocking is sensitive to resolution in current-generation atmospheric global circulation models (AGCMs). An evaluation is thus conducted of how well atmospheric blocking is represented in four AGCMs whose horizontal resolution is increased from a grid spacing of more than 100 km to about 25 km. It is shown that Euro-Atlantic blocking is simulated overall more credibly at higher resolution (i.e., in better agreement with a 50-yr reference blocking climatology created from the reanalyses ERA-40 and ERA-Interim). The improvement seen with resolution depends on the season andmore » to some extent on the model considered. Euro-Atlantic blocking is simulated more realistically at higher resolution in winter, spring, and autumn, and robustly so across the model ensemble. The improvement in spring is larger than that in winter and autumn. Summer blocking is found to be better simulated at higher resolution by one model only, with little change seen in the other three models. The representation of Pacific blocking is not found to systematically depend on resolution. Despite the improvements seen with resolution, the 25-km models still exhibit large biases in Euro-Atlantic blocking. For example, three of the four 25-km models underestimate winter northern European blocking frequency by about one-third. The resolution sensitivity and biases in the simulated blocking are shown to be in part associated with the mean-state biases in the models' midlatitude circulation.« less

The Resolution Sensitivity of Northern Hemisphere Blocking in Four 25-km Atmospheric Global Circulation Models

DOE PAGES

Schiemann, Reinhard; Demory, Marie-Estelle; Shaffrey, Len C.; ...

2016-12-19

The aim of this study is to investigate if the representation of Northern Hemisphere blocking is sensitive to resolution in current-generation atmospheric global circulation models (AGCMs). An evaluation is thus conducted of how well atmospheric blocking is represented in four AGCMs whose horizontal resolution is increased from a grid spacing of more than 100 km to about 25 km. It is shown that Euro-Atlantic blocking is simulated overall more credibly at higher resolution (i.e., in better agreement with a 50-yr reference blocking climatology created from the reanalyses ERA-40 and ERA-Interim). The improvement seen with resolution depends on the season andmore » to some extent on the model considered. Euro-Atlantic blocking is simulated more realistically at higher resolution in winter, spring, and autumn, and robustly so across the model ensemble. The improvement in spring is larger than that in winter and autumn. Summer blocking is found to be better simulated at higher resolution by one model only, with little change seen in the other three models. The representation of Pacific blocking is not found to systematically depend on resolution. Despite the improvements seen with resolution, the 25-km models still exhibit large biases in Euro-Atlantic blocking. For example, three of the four 25-km models underestimate winter northern European blocking frequency by about one-third. The resolution sensitivity and biases in the simulated blocking are shown to be in part associated with the mean-state biases in the models' midlatitude circulation.« less

Satellite-enhanced dynamical downscaling for the analysis of extreme events

NASA Astrophysics Data System (ADS)

Nunes, Ana M. B.

2016-09-01

The use of regional models in the downscaling of general circulation models provides a strategy to generate more detailed climate information. In that case, boundary-forcing techniques can be useful to maintain the large-scale features from the coarse-resolution global models in agreement with the inner modes of the higher-resolution regional models. Although those procedures might improve dynamics, downscaling via regional modeling still aims for better representation of physical processes. With the purpose of improving dynamics and physical processes in regional downscaling of global reanalysis, the Regional Spectral Model—originally developed at the National Centers for Environmental Prediction—employs a newly reformulated scale-selective bias correction, together with the 3-hourly assimilation of the satellite-based precipitation estimates constructed from the Climate Prediction Center morphing technique. The two-scheme technique for the dynamical downscaling of global reanalysis can be applied in analyses of environmental disasters and risk assessment, with hourly outputs, and resolution of about 25 km. Here the satellite-enhanced dynamical downscaling added value is demonstrated in simulations of the first reported hurricane in the western South Atlantic Ocean basin through comparisons with global reanalyses and satellite products available in ocean areas.

Effects of Drake Passage on a strongly eddying global ocean

NASA Astrophysics Data System (ADS)

Viebahn, Jan; von der Heydt, Anna S.; Dijkstra, Henk A.

2015-04-01

During the past 65 Million (Ma) years, Earth's climate has undergone a major change from warm 'greenhouse' to colder 'icehouse' conditions with extensive ice sheets in the polar regions of both hemispheres. The Eocene-Oligocene (~34 Ma) and Oligocene-Miocene (~23 Ma) boundaries reflect major transitions in Cenozoic global climate change. Proposed mechanisms of these transitions include reorganization of ocean circulation due to critical gateway opening/deepening, changes in atmospheric CO2-concentration, and feedback mechanisms related to land-ice formation. Drake Passage (DP) is an intensively studied gateway because it plays a central role in closing the transport pathways of heat and chemicals in the ocean. The climate response to a closed DP has been explored with a variety of general circulation models, however, all of these models employ low model-grid resolutions such that the effects of subgrid-scale fluctuations ('eddies') are parameterized. We present results of the first high-resolution (0.1° horizontally) realistic global ocean model simulation with a closed DP in which the eddy field is largely resolved. The simulation extends over more than 200 years such that the strong transient adjustment process is passed and a near-equilibrium ocean state is reached. The effects of DP are diagnosed by comparing with both an open DP high-resolution control simulation (of same length) and corresponding low-resolution simulations. By focussing on the heat/tracer transports we demonstrate that the results are twofold: Considering spatially integrated transports the overall response to a closed DP is well captured by low-resolution simulations. However, looking at the actual spatial distributions drastic differences appear between far-scattered high-resolution and laminar-uniform low-resolution fields. We conclude that sparse and highly localized tracer proxy observations have to be interpreted carefully with the help of high-resolution model simulations.

New global fire emission estimates and evaluation of volatile organic compounds

Treesearch

C. Wiedinmyer; L. K. Emmons; S. K. Akagi; R. J. Yokelson; J. J. Orlando; J. A. Al-Saadi; A. J. Soja

2010-01-01

A daily, high-resolution, global fire emissions model has been built to estimate emissions from open burning for air quality modeling applications: The Fire INventory from NCAR (FINN version 1). The model framework uses daily fire detections from the MODIS instruments and updated emission factors, specifically for speciated non-methane organic compounds (NMOC). Global...

WaterWorld, a spatial hydrological model applied at scales from local to global: key challenges to local application

NASA Astrophysics Data System (ADS)

Burke, Sophia; Mulligan, Mark

2017-04-01

WaterWorld is a widely used spatial hydrological policy support system. The last user census indicates regular use by 1029 institutions across 141 countries. A key feature of WaterWorld since 2001 is that it comes pre-loaded with all of the required data for simulation anywhere in the world at a 1km or 1 ha resolution. This means that it can be easily used, without specialist technical ability, to examine baseline hydrology and the impacts of scenarios for change or management interventions to support policy formulation, hence its labelling as a policy support system. WaterWorld is parameterised by an extensive global gridded database of more than 600 variables, developed from many sources, since 1998, the so-called simTerra database. All of these data are available globally at 1km resolution and some variables (terrain, land cover, urban areas, water bodies) are available globally at 1ha resolution. If users have access to better data than is pre-loaded, they can upload their own data. WaterWorld is generally applied at the national or basin scale at 1km resolution, or locally (for areas of <10,000km2) at 1ha resolution, though continental (1km resolution) and global (10km resolution) applications are possible so it is a model with local to global applications. WaterWorld requires some 140 maps to run including monthly climate data, land cover and use, terrain, population, water bodies and more. Whilst publically-available terrain and land cover data are now well developed for local scale application, climate and land use data remain a challenge, with most global products being available at 1km or 10km resolution or worse, which is rather coarse for local application. As part of the EartH2Observe project we have used WFDEI (WATCH Forcing Data methodology applied to ERA-Interim data) at 1km resolution to provide an alternative input to WaterWorld's preloaded climate data. Here we examine the impacts of that on key hydrological outputs: water balance, water quality and outline the remaining challenges of using datasets like these for local scale application.

One technique for refining the global Earth gravity models

NASA Astrophysics Data System (ADS)

Koneshov, V. N.; Nepoklonov, V. B.; Polovnev, O. V.

2017-01-01

The results of the theoretical and experimental research on the technique for refining the global Earth geopotential models such as EGM2008 in the continental regions are presented. The discussed technique is based on the high-resolution satellite data for the Earth's surface topography which enables the allowance for the fine structure of the Earth's gravitational field without the additional gravimetry data. The experimental studies are conducted by the example of the new GGMplus global gravity model of the Earth with a resolution about 0.5 km, which is obtained by expanding the EGM2008 model to degree 2190 with the corrections for the topograohy calculated from the SRTM data. The GGMplus and EGM2008 models are compared with the regional geoid models in 21 regions of North America, Australia, Africa, and Europe. The obtained estimates largely support the possibility of refining the global geopotential models such as EGM2008 by the procedure implemented in GGMplus, particularly in the regions with relatively high elevation difference.

Near-global climate simulation at 1 km resolution: establishing a performance baseline on 4888 GPUs with COSMO 5.0

NASA Astrophysics Data System (ADS)

Fuhrer, Oliver; Chadha, Tarun; Hoefler, Torsten; Kwasniewski, Grzegorz; Lapillonne, Xavier; Leutwyler, David; Lüthi, Daniel; Osuna, Carlos; Schär, Christoph; Schulthess, Thomas C.; Vogt, Hannes

2018-05-01

The best hope for reducing long-standing global climate model biases is by increasing resolution to the kilometer scale. Here we present results from an ultrahigh-resolution non-hydrostatic climate model for a near-global setup running on the full Piz Daint supercomputer on 4888 GPUs (graphics processing units). The dynamical core of the model has been completely rewritten using a domain-specific language (DSL) for performance portability across different hardware architectures. Physical parameterizations and diagnostics have been ported using compiler directives. To our knowledge this represents the first complete atmospheric model being run entirely on accelerators on this scale. At a grid spacing of 930 m (1.9 km), we achieve a simulation throughput of 0.043 (0.23) simulated years per day and an energy consumption of 596 MWh per simulated year. Furthermore, we propose a new memory usage efficiency (MUE) metric that considers how efficiently the memory bandwidth - the dominant bottleneck of climate codes - is being used.

A New High Resolution Climate Dataset for Climate Change Impacts Assessments in New England

NASA Astrophysics Data System (ADS)

Komurcu, M.; Huber, M.

2016-12-01

Assessing regional impacts of climate change (such as changes in extreme events, land surface hydrology, water resources, energy, ecosystems and economy) requires much higher resolution climate variables than those available from global model projections. While it is possible to run global models in higher resolution, the high computational cost associated with these simulations prevent their use in such manner. To alleviate this problem, dynamical downscaling offers a method to deliver higher resolution climate variables. As part of an NSF EPSCoR funded interdisciplinary effort to assess climate change impacts on New Hampshire ecosystems, hydrology and economy (the New Hampshire Ecosystems and Society project), we create a unique high-resolution climate dataset for New England. We dynamically downscale global model projections under a high impact emissions scenario using the Weather Research and Forecasting model (WRF) with three nested grids of 27, 9 and 3 km horizontal resolution with the highest resolution innermost grid focusing over New England. We prefer dynamical downscaling over other methods such as statistical downscaling because it employs physical equations to progressively simulate climate variables as atmospheric processes interact with surface processes, emissions, radiation, clouds, precipitation and other model components, hence eliminates fix relationships between variables. In addition to simulating mean changes in regional climate, dynamical downscaling also allows for the simulation of climate extremes that significantly alter climate change impacts. We simulate three time slices: 2006-2015, 2040-2060 and 2080-2100. This new high-resolution climate dataset (with more than 200 variables saved in hourly (six hourly) intervals for the highest resolution domain (outer two domains)) along with model input and restart files used in our WRF simulations will be publicly available for use to the broader scientific community to support in-depth climate change impacts assessments for New England. We present results focusing on future changes in New England extreme events.

The impact of resolution on the dynamics of the martian global atmosphere: Varying resolution studies with the MarsWRF GCM

NASA Astrophysics Data System (ADS)

Toigo, Anthony D.; Lee, Christopher; Newman, Claire E.; Richardson, Mark I.

2012-09-01

We investigate the sensitivity of the circulation and thermal structure of the martian atmosphere to numerical model resolution in a general circulation model (GCM) using the martian implementation (MarsWRF) of the planetWRF atmospheric model. We provide a description of the MarsWRF GCM and use it to study the global atmosphere at horizontal resolutions from 7.5° × 9° to 0.5° × 0.5°, encompassing the range from standard Mars GCMs to global mesoscale modeling. We find that while most of the gross-scale features of the circulation (the rough location of jets, the qualitative thermal structure, and the major large-scale features of the surface level winds) are insensitive to horizontal resolution over this range, several major features of the circulation are sensitive in detail. The northern winter polar circulation shows the greatest sensitivity, showing a continuous transition from a smooth polar winter jet at low resolution, to a distinct vertically “split” jet as resolution increases. The separation of the lower and middle atmosphere polar jet occurs at roughly 10 Pa, with the split jet structure developing in concert with the intensification of meridional jets at roughly 10 Pa and above 0.1 Pa. These meridional jets appear to represent the separation of lower and middle atmosphere mean overturning circulations (with the former being consistent with the usual concept of the “Hadley cell”). Further, the transition in polar jet structure is more sensitive to changes in zonal than meridional horizontal resolution, suggesting that representation of small-scale wave-mean flow interactions is more important than fine-scale representation of the meridional thermal gradient across the polar front. Increasing the horizontal resolution improves the match between the modeled thermal structure and the Mars Climate Sounder retrievals for northern winter high latitudes. While increased horizontal resolution also improves the simulation of the northern high latitudes at equinox, even the lowest model resolution considered here appears to do a good job for the southern winter and southern equinoctial pole (although in detail some discrepancies remain). These results suggest that studies of the northern winter jet (e.g., transient waves and cyclogenesis) will be more sensitive to global model resolution that those of the south (e.g., the confining dynamics of the southern polar vortex relevant to studies of argon transport). For surface winds, the major effect of increased horizontal resolution is in the superposition of circulations forced by local-scale topography upon the large-scale surface wind patterns. While passive predictions of dust lifting are generally insensitive to model horizontal resolution when no lifting threshold is considered, increasing the stress threshold produces significantly more lifting in higher resolution simulations with the generation of finer-scale, higher-stress winds due primarily to better-resolved topography. Considering the positive feedbacks expected for radiatively active dust lifting, we expect this bias to increase when such feedbacks are permitted.

Global River Water Temperature Modelling at Hyper-Resolution

NASA Astrophysics Data System (ADS)

Wanders, N.; van Vliet, M. T. H.; Wada, Y.; Van Beek, L. P.

2017-12-01

The temperature of river water plays a crucial role in many physical, chemical and biological aquatic processes. The influence of changing water temperatures is not only felt locally, but also has regional and downstream impacts. Sectors that might be affected by sudden or gradual changes in the water temperature are: energy production, industry and recreation. Although it is very important to have detailed information on this environmental variable, high-resolution simulations of water temperature on a large scale are currently lacking. Here we present a novel hyper-resolution water temperature dataset at the global scale. We developed the 1-D energy routing model WARM, to simulate river temperature for the period 1980-2014 at 10 km and 50 km resolution. The WARM model accounts for surface water abstraction, reservoirs, riverine flooding and formation of ice, therefore enabling a realistic representation of the water temperature. The water temperature simulations have been validated against 358 river monitoring stations globally for the period 1980 to 2014. The results indicate the increase in resolution significantly improves the simulation performance with a decrease in the water temperature RMSE from 3.5°C to 3.0°C and an increase in the mean correlation of the daily discharge simulations, from R=0.4 to 0.6. We find an average global increase in water temperature of 0.22°C per decade between 1960-2014, with increasing trends towards the end of the simulations period. Strong increasing trends in maxima in the Northern Hemisphere (0.62°C per decade) and minima in the Southern Hemisphere (0.45°C per decade). Finally, we show the impact of major heatwaves and drought events on the water temperature and water availability. The high resolution not only improves the model performance; it also positively impacts the relevancy of the simulation for local and regional scale studies and impact assessments. This new global water temperature dataset could help to develop decision-support system related to water quality with increasing precision and accuracy.

Assessing the Global Climate Response to Freshwater Forcing from the Antarctic Ice Sheet Under Future Climate Scenarios

NASA Astrophysics Data System (ADS)

Rogstad, S.; Condron, A.; DeConto, R.; Pollard, D.

2017-12-01

Observational evidence indicates that the West Antarctic Ice Sheet (WAIS) is losing mass at an accelerating rate. Impacts to global climate resulting from changing ocean circulation patterns due to increased freshwater runoff from Antarctica in the future could have significant implications for global heat transport, but to-date this topic has not been investigated using complex numerical models with realistic freshwater forcing. Here, we present results from a high resolution fully coupled ocean-atmosphere model (CESM 1.2) forced with runoff from Antarctica prescribed from a high resolution regional ice sheet-ice shelf model. Results from the regional simulations indicate a potential freshwater contribution from Antarctica of up to 1 m equivalent sea level rise by the end of the century under RCP 8.5 indicating that a substantial input of freshwater into the Southern Ocean is possible. Our high resolution global simulations were performed under IPCC future climate scenarios RCP 4.5 and 8.5. We will present results showing the impact of WAIS collapse on global ocean circulation, sea ice, air temperature, and salinity in order to assess the potential for abrupt climate change triggered by WAIS collapse.

Variable-Resolution Ensemble Climatology Modeling of Sierra Nevada Snowpack within the Community Earth System Model (CESM)

NASA Astrophysics Data System (ADS)

Rhoades, A.; Ullrich, P. A.; Zarzycki, C. M.; Levy, M.; Taylor, M.

2014-12-01

Snowpack is crucial for the western USA, providing around 75% of the total fresh water supply (Cayan et al., 1996) and buffering against seasonal aridity impacts on agricultural, ecosystem, and urban water demands. The resilience of the California water system is largely dependent on natural stores provided by snowpack. This resilience has shown vulnerabilities due to anthropogenic global climate change. Historically, the northern Sierras showed a net decline of 50-75% in snow water equivalent (SWE) while the southern Sierras showed a net accumulation of 30% (Mote et al., 2005). Future trends of SWE highlight that western USA SWE may decline by 40-70% (Pierce and Cayan, 2013), snowfall may decrease by 25-40% (Pierce and Cayan, 2013), and more winter storms may tend towards rain rather than snow (Bales et al., 2006). The volatility of Sierran snowpack presents a need for scientific tools to help water managers and policy makers assess current and future trends. A burgeoning tool to analyze these trends comes in the form of variable-resolution global climate modeling (VRGCM). VRGCMs serve as a bridge between regional and global models and provide added resolution in areas of need, eliminate lateral boundary forcings, provide model runtime speed up, and utilize a common dynamical core, physics scheme and sub-grid scale parameterization package. A cubed-sphere variable-resolution grid with 25 km horizontal resolution over the western USA was developed for use in the Community Atmosphere Model (CAM) within the Community Earth System Model (CESM). A 25-year three-member ensemble climatology (1980-2005) is presented and major snowpack metrics such as SWE, snow depth, snow cover, and two-meter surface temperature are assessed. The ensemble simulation is also compared to observational, reanalysis, and WRF model datasets. The variable-resolution model provides a mechanism for reaching towards non-hydrostatic scales and simulations are currently being developed with refined nests of 12.5km resolution over California.

Real-time Assimilation of Altimeter Derived Synthetic Profiles Into a Global version of the Naval Research Laboratory's Coastal Ocean Model (NCOM)

NASA Astrophysics Data System (ADS)

Rhodes, R. C.; Barron, C. N.; Fox, D. N.; Smedstad, L. F.

2001-12-01

A global implementation of the Navy Coastal Ocean Model (NCOM), developed by the Naval Research Laboratory (NRL) at Stennis Space Center is currently running in real-time and is planned for transition to the Naval Oceanographic Office (NAVOCEANO) in 2002. The model encompasses the open ocean to 5 m depth on a curvilinear global model grid with 1/8 degree grid spacing at 45N, extending from 80 S to a complete arctic cap with grid singularities mapped into Canada and Russia. Vertically, the model employs 41 sigma-z levels with sigma in the upper-ocean and coastal regions and z in the deeper ocean. The Navy Operational Global Atmospheric Prediction System (NOGAPS) provides 6-hourly wind stresses and heat fluxes for forcing, while the operational Modular Ocean Data Assimilation System (MODAS) provides the background climatology and tools for data pre-processing. Operationally available sea surface temperature (SST) and altimetry (SSH) data are assimilated into the NAVOCEANO global 1/8 degree MODAS 2-D analysis and the 1/16 degree Navy Layered Ocean Model (NLOM) to provide analyses and forecasts of SSH and SST. The 2-D SSH and SST nowcast fields are used as input to the MODAS synthetic climatology database to yield three-dimensional fields of synthetic temperature and salinity for assimilation into global NCOM. The synthetic profiles are weighted higher at depth in the assimilation process to allow the numerical model to properly develop the mixed-layer structure driven by the real-time atmospheric forcing. Global NCOM nowcasts and forecasts provide a valuable resource for rapid response to the varied and often unpredictable operational requests for 3-dimensional fields of ocean temperature, salinity, and currents. In some cases, the resolution of the global product is sufficient for guidance. In cases requiring higher resolution, the global product offers a quick overview of local circulation and provides initial and boundary conditions for higher resolution coastal models that may be more specialized for a particular task or domain. Nowcast and forecast results are presented globally and in selected areas of interest and model results are compared with historical and concurrent observations and analyses.

A high resolution global scale groundwater model

NASA Astrophysics Data System (ADS)

de Graaf, I. E.; Sutanudjaja, E.; Van Beek, L. P.; Bierkens, M. F.

2013-12-01

As the world's largest accessible source of freshwater, groundwater plays a vital role in satisfying the basic needs of human society. It serves as a primary source of drinking water and also supplies water for agricultural and industrial activities. During times of drought, the large natural groundwater storage provides a buffer against water shortage and sustains flows to rivers and wetlands, supporting ecosystem habitats and biodiversity. Yet, the current generation of global scale hydrological models (GHMs) do not include a groundwater flow component, although it is a crucial part of the hydrological cycle. Thus, a realistic physical representation of the groundwater system that allows for the simulation of groundwater head dynamics and lateral flows is essential for GHMs that increasingly run at finer resolution. In this study we present a transient global groundwater model with a resolution of 5 arc-minutes (approximately 10 km at the equator) using MODFLOW (McDonald and Harbaugh, 1988). Aquifer schematization and properties of this groundwater model were developed from available global lithological maps and datasets (Dürr et al., 2005; Gleeson et al., 2010; Hartmann and Moosdorf, 2013) combined with information about e.g. aquifer thickness and presence of less permeable, impermeable, and semi-impermeable layers. For the parameterization, we relied entirely on available global datasets and did not calibrate the model so that it can equally be expanded to data poor environments. We forced the groundwater model with the output from the global hydrological model PCR-GLOBWB (van Beek et al., 2011), specifically the net groundwater recharge and average surface water levels derived from routed channel discharge. We validated simulated groundwater heads with observations, from North America and Australia, resulting in a coefficient of determination of 0.8 and 0.7 respectively. This shows that it is feasible to build a global groundwater model using best available global information, and estimated water table depths are within acceptable accuracy in many parts of the world.

Weather extremes in very large, high-resolution ensembles: the weatherathome experiment

NASA Astrophysics Data System (ADS)

Allen, M. R.; Rosier, S.; Massey, N.; Rye, C.; Bowery, A.; Miller, J.; Otto, F.; Jones, R.; Wilson, S.; Mote, P.; Stone, D. A.; Yamazaki, Y. H.; Carrington, D.

2011-12-01

Resolution and ensemble size are often seen as alternatives in climate modelling. Models with sufficient resolution to simulate many classes of extreme weather cannot normally be run often enough to assess the statistics of rare events, still less how these statistics may be changing. As a result, assessments of the impact of external forcing on regional climate extremes must be based either on statistical downscaling from relatively coarse-resolution models, or statistical extrapolation from 10-year to 100-year events. Under the weatherathome experiment, part of the climateprediction.net initiative, we have compiled the Met Office Regional Climate Model HadRM3P to run on personal computer volunteered by the general public at 25 and 50km resolution, embedded within the HadAM3P global atmosphere model. With a global network of about 50,000 volunteers, this allows us to run time-slice ensembles of essentially unlimited size, exploring the statistics of extreme weather under a range of scenarios for surface forcing and atmospheric composition, allowing for uncertainty in both boundary conditions and model parameters. Current experiments, developed with the support of Microsoft Research, focus on three regions, the Western USA, Europe and Southern Africa. We initially simulate the period 1959-2010 to establish which variables are realistically simulated by the model and on what scales. Our next experiments are focussing on the Event Attribution problem, exploring how the probability of various types of extreme weather would have been different over the recent past in a world unaffected by human influence, following the design of Pall et al (2011), but extended to a longer period and higher spatial resolution. We will present the first results of the unique, global, participatory experiment and discuss the implications for the attribution of recent weather events to anthropogenic influence on climate.

Hybrid Multiscale Finite Volume method for multiresolution simulations of flow and reactive transport in porous media

NASA Astrophysics Data System (ADS)

Barajas-Solano, D. A.; Tartakovsky, A. M.

2017-12-01

We present a multiresolution method for the numerical simulation of flow and reactive transport in porous, heterogeneous media, based on the hybrid Multiscale Finite Volume (h-MsFV) algorithm. The h-MsFV algorithm allows us to couple high-resolution (fine scale) flow and transport models with lower resolution (coarse) models to locally refine both spatial resolution and transport models. The fine scale problem is decomposed into various "local'' problems solved independently in parallel and coordinated via a "global'' problem. This global problem is then coupled with the coarse model to strictly ensure domain-wide coarse-scale mass conservation. The proposed method provides an alternative to adaptive mesh refinement (AMR), due to its capacity to rapidly refine spatial resolution beyond what's possible with state-of-the-art AMR techniques, and the capability to locally swap transport models. We illustrate our method by applying it to groundwater flow and reactive transport of multiple species.

Projections of Future Summertime Ozone over the U.S.

DOE Office of Scientific and Technical Information (OSTI.GOV)

Pfister, G. G.; Walters, Stacy; Lamarque, J. F.

This study uses a regional fully coupled chemistry-transport model to assess changes in surface ozone over the summertime U.S. between present and a 2050 future time period at high spatial resolution (12 km grid spacing) under the SRES A2 climate and RCP8.5 anthropogenic pre-cursor emission scenario. The impact of predicted changes in climate and global background ozone is estimated to increase surface ozone over most of the U.S; the 5th - 95th percentile range for daily 8-hour maximum surface ozone increases from 31-79 ppbV to 30-87 ppbV between the present and future time periods. The analysis of a set ofmore » meteorological drivers suggests that these mostly will add to increasing ozone, but the set of simulations conducted does not allow to separate this effect from that through enhanced global background ozone. Statistically the most robust positive feedbacks are through increased temperature, biogenic emissions and solar radiation. Stringent emission controls can counteract these feedbacks and if considered, we estimate large reductions in surface ozone with the 5th-95th percentile reduced to 27-55 ppbV. A comparison of the high-resolution projections to global model projections shows that even though the global model is biased high in surface ozone compared to the regional model and compared to observations, both the global and the regional model predict similar changes in ozone between the present and future time periods. However, on smaller spatial scales, the regional predictions show more pronounced changes between urban and rural regimes that cannot be resolved at the coarse resolution of global model. In addition, the sign of the changes in overall ozone mixing ratios can be different between the global and the regional predictions in certain regions, such as the Western U.S. This study confirms the key role of emission control strategies in future air quality predictions and demonstrates the need for considering degradation of air quality with future climate change in emission policy making. It also illustrates the need for high resolution modeling when the objective is to address regional and local air quality or establish links to human health and society.« less

A Parallel Icosahedral, Higher Order Discontinuous Galerkin, Global Shallow Water Model: Global Ocean Tides and Aquaplanet Benchmarks

NASA Astrophysics Data System (ADS)

Salehipour, H.; Stuhne, G.; Peltier, W. R.

2012-12-01

The development of models of the ocean tides with higher resolution near the coastlines and courser mesh offshore, has been required due to the significant impacts of coastline configuration and bathymetry (associated with sea level rise) on the amplitude and phase of tidal constituents, not only under present conditions but also in the deep past [Griffiths and Peltier GRL 2008, Griffiths and Peltier AMS 2009, Hill et al. JGR 2011]. A global tidal model with enhanced resolution at the poles has been developed by Griffiths and Peltier [2008, 2009], which, although capable of highly resolving polar ocean tides , is based upon a standard structured Arakawa C grid and hence is not capable of resolving coastlines locally. Furthermore the use of a nested modelling approach, although it may enable local spatial refinement [Hill et al. 2011], nevertheless suffers from its inherent dependence on the availability of a global tidal model with necessarily low spatial resolution to provide the open boundary conditions required for the local high resolution model. On the other hand, an unstructured triangulation of the global domain provides a standalone framework that may be employed to study highly resolved regions without relying on secondary models. The first step in the development of the structure we are employing was described in Stuhne and Peltier [Ocean Modeling, 2009]. In further extending this modelling structure we are employing a new discontinuous Galerkin (DG) discretization of the governing equations in order to provide very high order of accuracy while also ensuring that momentum transport is locally conserved [Giraldo et al. JCP 2002]. After validating the 2D shallow water model with several test suites appropriate to aquaplanets [Williamson et al. JCP 1992, Galewsky et al. Tellus 2004, Nair and Lauritzen JCP 2010], the governing equations are extended to include the influence of internal tide drag in the deep ocean as well as the drag in shallow marginal seas together with the influence of gravitational self-attraction and loading. In this paper, we will explain the mathematical and numerical framework employed in the development of the DG global tidal model and present the validation results obtained using the present-day satellite altimetry data-constrained TPXO 6.2 global tidal solutions of Egbert et al. [JGR 1994].igure 1. Barotropic Instability Test of Galewsky et al. (Tellus 2004), with 2nd order DG

File Specification for the 7-km GEOS-5 Nature Run, Ganymed Release Non-Hydrostatic 7-km Global Mesoscale Simulation

NASA Technical Reports Server (NTRS)

da Silva, Arlindo M.; Putman, William; Nattala, J.

2014-01-01

This document describes the gridded output files produced by a two-year global, non-hydrostatic mesoscale simulation for the period 2005-2006 produced with the non-hydrostatic version of GEOS-5 Atmospheric Global Climate Model (AGCM). In addition to standard meteorological parameters (wind, temperature, moisture, surface pressure), this simulation includes 15 aerosol tracers (dust, sea-salt, sulfate, black and organic carbon), O3, CO and CO2. This model simulation is driven by prescribed sea-surface temperature and sea-ice, daily volcanic and biomass burning emissions, as well as high-resolution inventories of anthropogenic sources. A description of the GEOS-5 model configuration used for this simulation can be found in Putman et al. (2014). The simulation is performed at a horizontal resolution of 7 km using a cubed-sphere horizontal grid with 72 vertical levels, extending up to to 0.01 hPa (approximately 80 km). For user convenience, all data products are generated on two logically rectangular longitude-latitude grids: a full-resolution 0.0625 deg grid that approximately matches the native cubed-sphere resolution, and another 0.5 deg reduced-resolution grid. The majority of the full-resolution data products are instantaneous with some fields being time-averaged. The reduced-resolution datasets are mostly time-averaged, with some fields being instantaneous. Hourly data intervals are used for the reduced-resolution datasets, while 30-minute intervals are used for the full-resolution products. All full-resolution output is on the model's native 72-layer hybrid sigma-pressure vertical grid, while the reduced-resolution output is given on native vertical levels and on 48 pressure surfaces extending up to 0.02 hPa. Section 4 presents additional details on horizontal and vertical grids. Information of the model surface representation can be found in Appendix B. The GEOS-5 product is organized into file collections that are described in detail in Appendix C. Additional details about variables listed in this file specification can be found in a separate document, the GEOS-5 File Specification Variable Definition Glossary. Documentation about the current access methods for products described in this document can be found on the GEOS-5 Nature Run portal: http://gmao.gsfc.nasa.gov/projects/G5NR. Information on the scientific quality of this simulation will appear in a forthcoming NASA Technical Report Series on Global Modeling and Data Assimilation to be available from http://gmao.gsfc.nasa.gov/pubs/tm/.

Evaluation of a 3D local multiresolution algorithm for the correction of partial volume effects in positron emission tomography

PubMed Central

Le Pogam, Adrien; Hatt, Mathieu; Descourt, Patrice; Boussion, Nicolas; Tsoumpas, Charalampos; Turkheimer, Federico E.; Prunier-Aesch, Caroline; Baulieu, Jean-Louis; Guilloteau, Denis; Visvikis, Dimitris

2011-01-01

Purpose Partial volume effects (PVE) are consequences of the limited spatial resolution in emission tomography leading to under-estimation of uptake in tissues of size similar to the point spread function (PSF) of the scanner as well as activity spillover between adjacent structures. Among PVE correction methodologies, a voxel-wise mutual multi-resolution analysis (MMA) was recently introduced. MMA is based on the extraction and transformation of high resolution details from an anatomical image (MR/CT) and their subsequent incorporation into a low resolution PET image using wavelet decompositions. Although this method allows creating PVE corrected images, it is based on a 2D global correlation model which may introduce artefacts in regions where no significant correlation exists between anatomical and functional details. Methods A new model was designed to overcome these two issues (2D only and global correlation) using a 3D wavelet decomposition process combined with a local analysis. The algorithm was evaluated on synthetic, simulated and patient images, and its performance was compared to the original approach as well as the geometric transfer matrix (GTM) method. Results Quantitative performance was similar to the 2D global model and GTM in correlated cases. In cases where mismatches between anatomical and functional information were present the new model outperformed the 2D global approach, avoiding artefacts and significantly improving quality of the corrected images and their quantitative accuracy. Conclusions A new 3D local model was proposed for a voxel-wise PVE correction based on the original mutual multi-resolution analysis approach. Its evaluation demonstrated an improved and more robust qualitative and quantitative accuracy compared to the original MMA methodology, particularly in the absence of full correlation between anatomical and functional information. PMID:21978037

DOE Office of Scientific and Technical Information (OSTI.GOV)

Sakaguchi, Koichi; Leung, Lai-Yung R.; Zhao, Chun

This study presents a diagnosis of a multi-resolution approach using the Model for Prediction Across Scales - Atmosphere (MPAS-A) for simulating regional climate. Four AMIP experiments are conducted for 1999-2009. In the first two experiments, MPAS-A is configured using global quasi-uniform grids at 120 km and 30 km grid spacing. In the other two experiments, MPAS-A is configured using variable-resolution (VR) mesh with local refinement at 30 km over North America and South America embedded inside a quasi-uniform domain at 120 km elsewhere. Precipitation and related fields in the four simulations are examined to determine how well the VR simulationsmore » reproduce the features simulated by the globally high-resolution model in the refined domain. In previous analyses of idealized aqua-planet simulations, the characteristics of the global high-resolution simulation in moist processes only developed near the boundary of the refined region. In contrast, the AMIP simulations with VR grids are able to reproduce the high-resolution characteristics across the refined domain, particularly in South America. This indicates the importance of finely resolved lower-boundary forcing such as topography and surface heterogeneity for the regional climate, and demonstrates the ability of the MPAS-A VR to replicate the large-scale moisture transport as simulated in the quasi-uniform high-resolution model. Outside of the refined domain, some upscale effects are detected through large-scale circulation but the overall climatic signals are not significant at regional scales. Our results provide support for the multi-resolution approach as a computationally efficient and physically consistent method for modeling regional climate.« less

Impact of Parameterized Lee Wave Drag on the Energy Budget of an Eddying Global Ocean Model

DTIC Science & Technology

2013-08-26

Teixeira, J., Peng, M., Hogan, T.F., Pauley, R., 2002. Navy Operational Global Atmospheric Prediction System (NOGAPS): Forcing for ocean models...Impact of parameterized lee wave drag on the energy budget of an eddying global ocean model David S. Trossman a,⇑, Brian K. Arbic a, Stephen T...input and output terms in the total mechanical energy budget of a hybrid coordinate high-resolution global ocean general circulation model forced by winds

Tethys – A Python Package for Spatial and Temporal Downscaling of Global Water Withdrawals

DOE Office of Scientific and Technical Information (OSTI.GOV)

Li, Xinya; Vernon, Chris R.; Hejazi, Mohamad I.

Downscaling of water withdrawals from regional/national to local scale is a fundamental step and also a common problem when integrating large scale economic and integrated assessment models with high-resolution detailed sectoral models. Tethys, an open-access software written in Python, is developed with statistical downscaling algorithms, to spatially and temporally downscale water withdrawal data to a finer scale. The spatial resolution will be downscaled from region/basin scale to grid (0.5 geographic degree) scale and the temporal resolution will be downscaled from year to month. Tethys is used to produce monthly global gridded water withdrawal products based on estimates from the Globalmore » Change Assessment Model (GCAM).« less

Tethys – A Python Package for Spatial and Temporal Downscaling of Global Water Withdrawals

DOE PAGES

Li, Xinya; Vernon, Chris R.; Hejazi, Mohamad I.; ...

2018-02-09

Downscaling of water withdrawals from regional/national to local scale is a fundamental step and also a common problem when integrating large scale economic and integrated assessment models with high-resolution detailed sectoral models. Tethys, an open-access software written in Python, is developed with statistical downscaling algorithms, to spatially and temporally downscale water withdrawal data to a finer scale. The spatial resolution will be downscaled from region/basin scale to grid (0.5 geographic degree) scale and the temporal resolution will be downscaled from year to month. Tethys is used to produce monthly global gridded water withdrawal products based on estimates from the Globalmore » Change Assessment Model (GCAM).« less

Downscaling Global Land Cover Projections from an Integrated Assessment Model for Use in Regional Analyses: Results and Evaluation for the US from 2005 to 2095

DOE Office of Scientific and Technical Information (OSTI.GOV)

West, Tristram O.; Le Page, Yannick LB; Huang, Maoyi

2014-06-05

Projections of land cover change generated from Integrated Assessment Models (IAM) and other economic-based models can be applied for analyses of environmental impacts at subregional and landscape scales. For those IAM and economic models that project land use at the sub-continental or regional scale, these projections must be downscaled and spatially distributed prior to use in climate or ecosystem models. Downscaling efforts to date have been conducted at the national extent with relatively high spatial resolution (30m) and at the global extent with relatively coarse spatial resolution (0.5 degree).

Realism of Indian Summer Monsoon Simulation in a Quarter Degree Global Climate Model

NASA Astrophysics Data System (ADS)

Salunke, P.; Mishra, S. K.; Sahany, S.; Gupta, K.

2017-12-01

This study assesses the fidelity of Indian Summer Monsoon (ISM) simulations using a global model at an ultra-high horizontal resolution (UHR) of 0.25°. The model used was the atmospheric component of the Community Earth System Model version 1.2.0 (CESM 1.2.0) developed at the National Center for Atmospheric Research (NCAR). Precipitation and temperature over the Indian region were analyzed for a wide range of space and time scales to evaluate the fidelity of the model under UHR, with special emphasis on the ISM simulations during the period of June-through-September (JJAS). Comparing the UHR simulations with observed data from the India Meteorological Department (IMD) over the Indian land, it was found that 0.25° resolution significantly improved spatial rainfall patterns over many regions, including the Western Ghats and the South-Eastern peninsula as compared to the standard model resolution. Convective and large-scale rainfall components were analyzed using the European Centre for Medium Range Weather Forecast (ECMWF) Re-Analysis (ERA)-Interim (ERA-I) data and it was found that at 0.25° resolution, there was an overall increase in the large-scale component and an associated decrease in the convective component of rainfall as compared to the standard model resolution. Analysis of the diurnal cycle of rainfall suggests a significant improvement in the phase characteristics simulated by the UHR model as compared to the standard model resolution. Analysis of the annual cycle of rainfall, however, failed to show any significant improvement in the UHR model as compared to the standard version. Surface temperature analysis showed small improvements in the UHR model simulations as compared to the standard version. Thus, one may conclude that there are some significant improvements in the ISM simulations using a 0.25° global model, although there is still plenty of scope for further improvement in certain aspects of the annual cycle of rainfall.

Century long observation constrained global dynamic downscaling and hydrologic implication

NASA Astrophysics Data System (ADS)

Kim, H.; Yoshimura, K.; Chang, E.; Famiglietti, J. S.; Oki, T.

2012-12-01

It has been suggested that greenhouse gas induced warming climate causes the acceleration of large scale hydrologic cycles, and, indeed, many regions on the Earth have been suffered by hydrologic extremes getting more frequent. However, historical observations are not able to provide enough information in comprehensive manner to understand their long-term variability and/or global distributions. In this study, a century long high resolution global climate data is developed in order to break through existing limitations. 20th Century Reanalysis (20CR) which has relatively low spatial resolution (~2.0°) and longer term availability (140 years) is dynamically downscaled into global T248 (~0.5°) resolution using Experimental Climate Prediction Center (ECPC) Global Spectral Model (GSM) by spectral nudging data assimilation technique. Also, Global Precipitation Climatology Centre (GPCC) and Climate Research Unit (CRU) observational data are adopted to reduce model dependent uncertainty. Downscaled product successfully represents realistic geographical detail keeping low frequency signal in mean state and spatiotemporal variability, while previous bias correction method fails to reproduce high frequency variability. Newly developed data is used to investigate how long-term large scale terrestrial hydrologic cycles have been changed globally and how they have been interacted with various climate modes, such as El-Niño Southern Oscillation (ENSO) and Atlantic Multidecadal Oscillation (AMO). As a further application, it will be used to provide atmospheric boundary condition of multiple land surface models in the Global Soil Wetness Project Phase 3 (GSWP3).

Global multi-resolution terrain elevation data 2010 (GMTED2010)

USGS Publications Warehouse

Danielson, Jeffrey J.; Gesch, Dean B.

2011-01-01

In 1996, the U.S. Geological Survey (USGS) developed a global topographic elevation model designated as GTOPO30 at a horizontal resolution of 30 arc-seconds for the entire Earth. Because no single source of topographic information covered the entire land surface, GTOPO30 was derived from eight raster and vector sources that included a substantial amount of U.S. Defense Mapping Agency data. The quality of the elevation data in GTOPO30 varies widely; there are no spatially-referenced metadata, and the major topographic features such as ridgelines and valleys are not well represented. Despite its coarse resolution and limited attributes, GTOPO30 has been widely used for a variety of hydrological, climatological, and geomorphological applications as well as military applications, where a regional, continental, or global scale topographic model is required. These applications have ranged from delineating drainage networks and watersheds to using digital elevation data for the extraction of topographic structure and three-dimensional (3D) visualization exercises (Jenson and Domingue, 1988; Verdin and Greenlee, 1996; Lehner and others, 2008). Many of the fundamental geophysical processes active at the Earth's surface are controlled or strongly influenced by topography, thus the critical need for high-quality terrain data (Gesch, 1994). U.S. Department of Defense requirements for mission planning, geographic registration of remotely sensed imagery, terrain visualization, and map production are similarly dependent on global topographic data. Since the time GTOPO30 was completed, the availability of higher-quality elevation data over large geographic areas has improved markedly. New data sources include global Digital Terrain Elevation Data (DTEDRegistered) from the Shuttle Radar Topography Mission (SRTM), Canadian elevation data, and data from the Ice, Cloud, and land Elevation Satellite (ICESat). Given the widespread use of GTOPO30 and the equivalent 30-arc-second DTEDRegistered level 0, the USGS and the National Geospatial-Intelligence Agency (NGA) have collaborated to produce an enhanced replacement for GTOPO30, the Global Land One-km Base Elevation (GLOBE) model and other comparable 30-arc-second-resolution global models, using the best available data. The new model is called the Global Multi-resolution Terrain Elevation Data 2010, or GMTED2010 for short. This suite of products at three different resolutions (approximately 1,000, 500, and 250 meters) is designed to support many applications directly by providing users with generic products (for example, maximum, minimum, and median elevations) that have been derived directly from the raw input data that would not be available to the general user or would be very costly and time-consuming to produce for individual applications. The source of all the elevation data is captured in metadata for reference purposes. It is also hoped that as better data become available in the future, the GMTED2010 model will be updated.

A new synoptic scale resolving global climate simulation using the Community Earth System Model

NASA Astrophysics Data System (ADS)

Small, R. Justin; Bacmeister, Julio; Bailey, David; Baker, Allison; Bishop, Stuart; Bryan, Frank; Caron, Julie; Dennis, John; Gent, Peter; Hsu, Hsiao-ming; Jochum, Markus; Lawrence, David; Muñoz, Ernesto; diNezio, Pedro; Scheitlin, Tim; Tomas, Robert; Tribbia, Joseph; Tseng, Yu-heng; Vertenstein, Mariana

2014-12-01

High-resolution global climate modeling holds the promise of capturing planetary-scale climate modes and small-scale (regional and sometimes extreme) features simultaneously, including their mutual interaction. This paper discusses a new state-of-the-art high-resolution Community Earth System Model (CESM) simulation that was performed with these goals in mind. The atmospheric component was at 0.25° grid spacing, and ocean component at 0.1°. One hundred years of "present-day" simulation were completed. Major results were that annual mean sea surface temperature (SST) in the equatorial Pacific and El-Niño Southern Oscillation variability were well simulated compared to standard resolution models. Tropical and southern Atlantic SST also had much reduced bias compared to previous versions of the model. In addition, the high resolution of the model enabled small-scale features of the climate system to be represented, such as air-sea interaction over ocean frontal zones, mesoscale systems generated by the Rockies, and Tropical Cyclones. Associated single component runs and standard resolution coupled runs are used to help attribute the strengths and weaknesses of the fully coupled run. The high-resolution run employed 23,404 cores, costing 250 thousand processor-hours per simulated year and made about two simulated years per day on the NCAR-Wyoming supercomputer "Yellowstone."

How does a High Resolution Global Model represent Mesoscale Convective Systems over the Amazon Basin?

NASA Astrophysics Data System (ADS)

Rehbein, A.; Ambrizzi, T.

2017-12-01

The mesoscale convective systems (MCSs) are very important meteorological systems, which can impact on the local, regional and global climate. Despite of their importance, the knowledge about their occurrence and behavior is still poor, mainly over the tropical region of South America where the data availability is scarce. Besides, few attentions are given to represent the MCSs in the numerical modeling in that region. The aim of the present work is to evaluate the representation of the MCSs by a global high resolution model over the Amazon basin. In this study, we will make a revision of the state of art involving the MCSs' over the Amazon basin and also how they are represented. For this last point, we will identify and track the MCSs using precipitation data from a high resolution nonhydrostatic global model, called Non-hydrostatic ICosahedral Atmospheric Model (NICAM). The spatial and temporal resolution of NICAM are 14 km and 1 hour, respectively. The MCSs identification and tracking will be performed by the algorithm Forecast and Tracking the evolution of Cloud Clusters (ForTraCC) for the period of 2000 to 2008. This will allow us evaluate the representation of the MCSs obtained by NICAM and compare them with those found using infrared satellite images. NICAM's precipitation was validated using Climate Hazards Group InfraRed Precipitation with Station data (CHIRPS), from 1981 to 2008. Once the model is validated, we will analyze the variability of the MCSs using the simulations of the NICAM for a future climate.

Global hydrodynamic modelling of flood inundation in continental rivers: How can we achieve it?

NASA Astrophysics Data System (ADS)

Yamazaki, D.

2016-12-01

Global-scale modelling of river hydrodynamics is essential for understanding global hydrological cycle, and is also required in interdisciplinary research fields . Global river models have been developed continuously for more than two decades, but modelling river flow at a global scale is still a challenging topic because surface water movement in continental rivers is a multi-spatial-scale phenomena. We have to consider the basin-wide water balance (>1000km scale), while hydrodynamics in river channels and floodplains is regulated by much smaller-scale topography (<100m scale). For example, heavy precipitation in upstream regions may later cause flooding in farthest downstream reaches. In order to realistically simulate the timing and amplitude of flood wave propagation for a long distance, consideration of detailed local topography is unavoidable. I have developed the global hydrodynamic model CaMa-Flood to overcome this scale-discrepancy of continental river flow. The CaMa-Flood divides river basins into multiple "unit-catchments", and assumes the water level is uniform within each unit-catchment. One unit-catchment is assigned to each grid-box defined at the typical spatial resolution of global climate models (10 100 km scale). Adopting a uniform water level in a >10km river segment seems to be a big assumption, but it is actually a good approximation for hydrodynamic modelling of continental rivers. The number of grid points required for global hydrodynamic simulations is largely reduced by this "unit-catchment assumption". Alternative to calculating 2-dimensional floodplain flows as in regional flood models, the CaMa-Flood treats floodplain inundation in a unit-catchment as a sub-grid physics. The water level and inundated area in each unit-catchment are diagnosed from water volume using topography parameters derived from high-resolution digital elevation models. Thus, the CaMa-Flood is at least 1000 times computationally more efficient compared to regional flood inundation models while the reality of simulated flood dynamics is kept. I will explain in detail how the CaMa-Flood model has been constructed from high-resolution topography datasets, and how the model can be used for various interdisciplinary applications.

Towards the impact of eddies on the response of the global ocean circulation to Southern Ocean gateway opening

NASA Astrophysics Data System (ADS)

Viebahn, Jan; von der Heydt, Anna S.; Dijkstra, Henk A.

2014-05-01

During the past 65 Million (Ma) years, Earth's climate has undergone a major change from warm 'greenhouse' to colder 'icehouse' conditions with extensive ice sheets in the polar regions of both hemispheres. The Eocene-Oligocene (~34 Ma) and Oligocene-Miocene (~23 Ma) boundaries reflect major transitions in Cenozoic global climate change. Proposed mechanisms of these transitions include reorganization of ocean circulation due to critical gateway opening/deepening, changes in atmospheric CO2-concentration, and feedback mechanisms related to land-ice formation. A long-standing hypothesis is that the formation of the Antarctic Circumpolar Current due to opening/deepening of Southern Ocean gateways led to glaciation of the Antarctic continent. However, while this hypothesis remains controversial, its assessment via coupled climate model simulations depends crucially on the spatial resolution in the ocean component. More precisely, only high-resolution modeling of the turbulent ocean circulation is capable of adequately describing reorganizations in the ocean flow field and related changes in turbulent heat transport. In this study, for the first time results of a high-resolution (0.1° horizontally) realistic global ocean model simulation with a closed Drake Passage are presented. Changes in global ocean temperatures, heat transport, and ocean circulation (e.g., Meridional Overturning Circulation and Antarctic Coastal Current) are established by comparison with an open Drake Passage high-resolution reference simulation. Finally, corresponding low-resolution simulations are also analyzed. The results highlight the essential impact of the ocean eddy field in palaeoclimatic change.

Earth System Modeling 2.0: A Blueprint for Models That Learn From Observations and Targeted High-Resolution Simulations

NASA Astrophysics Data System (ADS)

Schneider, Tapio; Lan, Shiwei; Stuart, Andrew; Teixeira, João.

2017-12-01

Climate projections continue to be marred by large uncertainties, which originate in processes that need to be parameterized, such as clouds, convection, and ecosystems. But rapid progress is now within reach. New computational tools and methods from data assimilation and machine learning make it possible to integrate global observations and local high-resolution simulations in an Earth system model (ESM) that systematically learns from both and quantifies uncertainties. Here we propose a blueprint for such an ESM. We outline how parameterization schemes can learn from global observations and targeted high-resolution simulations, for example, of clouds and convection, through matching low-order statistics between ESMs, observations, and high-resolution simulations. We illustrate learning algorithms for ESMs with a simple dynamical system that shares characteristics of the climate system; and we discuss the opportunities the proposed framework presents and the challenges that remain to realize it.

GEM-AQ, an On-line Global Multiscale Chemical Weather System: Model Description and Evaluation of Gas Phase Chemistry Processes

NASA Astrophysics Data System (ADS)

Neary, L.; Kaminski, J. W.; Struzewska, J.; Ainslie, B.; McConnell, J. C.

2007-12-01

Tropospheric chemistry and air quality processes were implemented on-line in the Global Environmental Multiscale model. The integrated model, GEM-AQ, has been developed as a platform to investigate chemical weather at scales from global to urban. On the global scale, the model was exercised for five years (2001-2005) to evaluate its ability to simulate seasonal variations and regional distributions of trace gases such as ozone, nitrogen dioxide and carbon monoxide. The model results are compared with observations from satellites, aircraft measurement campaigns and balloon sondes. The same model has also been evaluated on the regional (~15km resolution) and urban scale (~3km resolution). A simulation of the formation and transport of photooxidants during the European heat wave of 2006 was performed and compared with surface observations throughout central and eastern Europe. The complex topographic region of the Lower Fraser Valley in British Columbia was the focus of another model evaluation during the PACIFIC 2001 field campaign. Comparison of model results with observations during this period will be shown.

A description and evaluation of an air quality model nested within global and regional composition-climate models using MetUM

NASA Astrophysics Data System (ADS)

Neal, Lucy S.; Dalvi, Mohit; Folberth, Gerd; McInnes, Rachel N.; Agnew, Paul; O'Connor, Fiona M.; Savage, Nicholas H.; Tilbee, Marie

2017-11-01

There is a clear need for the development of modelling frameworks for both climate change and air quality to help inform policies for addressing these issues simultaneously. This paper presents an initial attempt to develop a single modelling framework, by introducing a greater degree of consistency in the meteorological modelling framework by using a two-step, one-way nested configuration of models, from a global composition-climate model (GCCM) (140 km resolution) to a regional composition-climate model covering Europe (RCCM) (50 km resolution) and finally to a high (12 km) resolution model over the UK (AQUM). The latter model is used to produce routine air quality forecasts for the UK. All three models are based on the Met Office's Unified Model (MetUM). In order to better understand the impact of resolution on the downscaling of projections of future climate and air quality, we have used this nest of models to simulate a 5-year period using present-day emissions and under present-day climate conditions. We also consider the impact of running the higher-resolution model with higher spatial resolution emissions, rather than simply regridding emissions from the RCCM. We present an evaluation of the models compared to in situ air quality observations over the UK, plus a comparison against an independent 1 km resolution gridded dataset, derived from a combination of modelling and observations, effectively producing an analysis of annual mean surface pollutant concentrations. We show that using a high-resolution model over the UK has some benefits in improving air quality modelling, but that the use of higher spatial resolution emissions is important to capture local variations in concentrations, particularly for primary pollutants such as nitrogen dioxide and sulfur dioxide. For secondary pollutants such as ozone and the secondary component of PM10, the benefits of a higher-resolution nested model are more limited and reasons for this are discussed. This study highlights the point that the resolution of models is not the only factor in determining model performance - consistency between nested models is also important.

Application of an Integrated Assessment Model with state-level resolution for examining strategies for addressing air, climate and energy goals

EPA Science Inventory

The Global Climate Assessment Model (GCAM) is a global integrated assessment model used for exploring future scenarios and examining strategies that address air pollution, climate change, and energy goals. GCAM includes technology-rich representations of the energy, transportati...

Estimation of Global 1km-grid Terrestrial Carbon Exchange Part II: Evaluations and Applications

NASA Astrophysics Data System (ADS)

Murakami, K.; Sasai, T.; Kato, S.; Niwa, Y.; Saito, M.; Takagi, H.; Matsunaga, T.; Hiraki, K.; Maksyutov, S. S.; Yokota, T.

2015-12-01

Global terrestrial carbon cycle largely depends on a spatial pattern in land cover type, which is heterogeneously-distributed over regional and global scales. Many studies have been trying to reveal distribution of carbon exchanges between terrestrial ecosystems and atmosphere for understanding global carbon cycle dynamics by using terrestrial biosphere models, satellite data, inventory data, and so on. However, most studies remained within several tens of kilometers grid spatial resolution, and the results have not been enough to understand the detailed pattern of carbon exchanges based on ecological community and to evaluate the carbon stocks by forest ecosystems in each countries. Improving the sophistication of spatial resolution is obviously necessary to enhance the accuracy of carbon exchanges. Moreover, the improvement may contribute to global warming awareness, policy makers and other social activities. We show global terrestrial carbon exchanges (net ecosystem production, net primary production, and gross primary production) with 1km-grid resolution. The methodology for these estimations are shown in the 2015 AGU FM poster "Estimation of Global 1km-grid Terrestrial Carbon Exchange Part I: Developing Inputs and Modelling". In this study, we evaluated the carbon exchanges in various regions with other approaches. We used the satellite-driven biosphere model (BEAMS) as our estimations, GOSAT L4A CO2 flux data, NEP retrieved by NICAM and CarbonTracer2013 flux data, for period from Jun 2001 to Dec 2012. The temporal patterns for this period were indicated similar trends between BEAMS, GOSAT, NICAM, and CT2013 in many sub-continental regions. Then, we estimated the terrestrial carbon exchanges in each countries, and could indicated the temporal patterns of the exchanges in large carbon stock regions.Global terrestrial carbon cycle largely depends on a spatial pattern of land cover type, which is heterogeneously-distributed over regional and global scales. Many studies have been trying to reveal distribution of carbon exchanges between terrestrial ecosystems and atmosphere for understanding global carbon cycle dynamics by using terrestrial biosphere models, satellite data, inventory data, and so on. However, most studies remained within several tens of kilometers grid spatial resolution, and the results have not been enough to understand the detailed pattern of carbon exchanges based on ecological community and to evaluate the carbon stocks by forest ecosystems in each countries. Improving the sophistication of spatial resolution is obviously necessary to enhance the accuracy of carbon exchanges. Moreover, the improvement may contribute to global warming awareness, policy makers and other social activities. We show global terrestrial carbon exchanges (net ecosystem production, net primary production, and gross primary production) with 1km-grid resolution. The methodology for these estimations are shown in the 2015 AGU FM poster "Estimation of Global 1km-grid Terrestrial Carbon Exchange Part I: Developing Inputs and Modelling". In this study, we evaluated the carbon exchanges in various regions with other approaches. We used the satellite-driven biosphere model (BEAMS) as our estimations, GOSAT L4A CO2 flux data, NEP retrieved by NICAM and CarbonTracer2013 flux data, for period from Jun 2001 to Dec 2012. The temporal patterns for this period were indicated similar trends between BEAMS, GOSAT, NICAM, and CT2013 in many sub-continental regions. Then, we estimated the terrestrial carbon exchanges in each countries, and could indicated the temporal patterns of the exchanges in large carbon stock regions.

Final Technical Report

DOE Office of Scientific and Technical Information (OSTI.GOV)

Held, Isaac; V. Balaji; Fueglistaler, Stephan

We have constructed and analyzed a series of idealized models of tropical convection interacting with large-scale circulations, with 25-50km resolution and with 1-2km cloud resolving resolution to set the stage for rigorous tests of convection closure schemes in high resolution global climate models. Much of the focus has been on the climatology of tropical cyclogenesis in rotating systems and the related problem of the spontaneous aggregation of convection in non-rotating systems. The PI (Held) will be delivering the honorary Bjerknes lecture at the Fall 2016 AGU meeting in December on this work. We have also provided new analyses of long-standingmore » issues related to the interaction between convection and the large-scale circulation: Kelvin waves in the upper troposphere and lower stratosphere, water vapor transport into the stratosphere, and upper tropospheric temperature trends. The results of these analyses help to improve our understanding of processes, and provide tests for future high resolution global modeling. Our final goal of testing new convections schemes in next-generation global atmospheric models at GFDL has been left for future work due to the complexity of the idealized model results meant as tests for these models uncovered in this work and to computational resource limitations. 11 papers have been published with support from this grant, 2 are in review, and another major summary paper is in preparation.« less

Spectral decomposition of internal gravity wave sea surface height in global models

NASA Astrophysics Data System (ADS)

Savage, Anna C.; Arbic, Brian K.; Alford, Matthew H.; Ansong, Joseph K.; Farrar, J. Thomas; Menemenlis, Dimitris; O'Rourke, Amanda K.; Richman, James G.; Shriver, Jay F.; Voet, Gunnar; Wallcraft, Alan J.; Zamudio, Luis

2017-10-01

Two global ocean models ranging in horizontal resolution from 1/12° to 1/48° are used to study the space and time scales of sea surface height (SSH) signals associated with internal gravity waves (IGWs). Frequency-horizontal wavenumber SSH spectral densities are computed over seven regions of the world ocean from two simulations of the HYbrid Coordinate Ocean Model (HYCOM) and three simulations of the Massachusetts Institute of Technology general circulation model (MITgcm). High wavenumber, high-frequency SSH variance follows the predicted IGW linear dispersion curves. The realism of high-frequency motions (>0.87 cpd) in the models is tested through comparison of the frequency spectral density of dynamic height variance computed from the highest-resolution runs of each model (1/25° HYCOM and 1/48° MITgcm) with dynamic height variance frequency spectral density computed from nine in situ profiling instruments. These high-frequency motions are of particular interest because of their contributions to the small-scale SSH variability that will be observed on a global scale in the upcoming Surface Water and Ocean Topography (SWOT) satellite altimetry mission. The variance at supertidal frequencies can be comparable to the tidal and low-frequency variance for high wavenumbers (length scales smaller than ˜50 km), especially in the higher-resolution simulations. In the highest-resolution simulations, the high-frequency variance can be greater than the low-frequency variance at these scales.

Hyperresolution Global Land Surface Modeling: Meeting a Grand Challenge for Monitoring Earth's Terrestrial Water

NASA Technical Reports Server (NTRS)

Wood, Eric F.; Roundy, Joshua K.; Troy, Tara J.; van Beek, L. P. H.; Bierkens, Marc F. P.; 4 Blyth, Eleanor; de Roo, Ad; Doell. Petra; Ek, Mike; Famiglietti, James;

Hyperresolution global land surface modeling: Meeting a grand challenge for monitoring Earth's terrestrial water

NASA Astrophysics Data System (ADS)

Wood, Eric F.; Roundy, Joshua K.; Troy, Tara J.; van Beek, L. P. H.; Bierkens, Marc F. P.; Blyth, Eleanor; de Roo, Ad; DöLl, Petra; Ek, Mike; Famiglietti, James; Gochis, David; van de Giesen, Nick; Houser, Paul; Jaffé, Peter R.; Kollet, Stefan; Lehner, Bernhard; Lettenmaier, Dennis P.; Peters-Lidard, Christa; Sivapalan, Murugesu; Sheffield, Justin; Wade, Andrew; Whitehead, Paul

2011-05-01

Monitoring Earth's terrestrial water conditions is critically important to many hydrological applications such as global food production; assessing water resources sustainability; and flood, drought, and climate change prediction. These needs have motivated the development of pilot monitoring and prediction systems for terrestrial hydrologic and vegetative states, but to date only at the rather coarse spatial resolutions (˜10-100 km) over continental to global domains. Adequately addressing critical water cycle science questions and applications requires systems that are implemented globally at much higher resolutions, on the order of 1 km, resolutions referred to as hyperresolution in the context of global land surface models. This opinion paper sets forth the needs and benefits for a system that would monitor and predict the Earth's terrestrial water, energy, and biogeochemical cycles. We discuss six major challenges in developing a system: improved representation of surface-subsurface interactions due to fine-scale topography and vegetation; improved representation of land-atmospheric interactions and resulting spatial information on soil moisture and evapotranspiration; inclusion of water quality as part of the biogeochemical cycle; representation of human impacts from water management; utilizing massively parallel computer systems and recent computational advances in solving hyperresolution models that will have up to 109 unknowns; and developing the required in situ and remote sensing global data sets. We deem the development of a global hyperresolution model for monitoring the terrestrial water, energy, and biogeochemical cycles a "grand challenge" to the community, and we call upon the international hydrologic community and the hydrological science support infrastructure to endorse the effort.

The Kain-Fritsch Scheme: Science Updates & Revisiting Gray-Scale Issues from the NWP & Regional Climatae Perspectives

EPA Science Inventory

It’s just a matter of time before we see global climate models increasing their spatial resolution to that now typical of regional models. This encroachment brings in an urgent need for making regional NWP and climate models applicable at certain finer resolutions. One of the hin...

Projecting 21st Century Snowpack Trends in the Western United States using Variable-Resolution CESM

NASA Astrophysics Data System (ADS)

Rhoades, A.; Huang, X.; Zarzycki, C. M.; Ullrich, P. A.

2015-12-01

The western USA is integrally reliant upon winter season snowpack, which supplies 3/4 of the region's fresh water and buffers against seasonal aridity on agricultural, ecosystem, and urban water demands. By the end of the 21st century, western USA snowpack (SWE) could decline by 40-70%, snowfall by 25-40%, more winter storms could tend towards rain rather than snow, and the peak timing of snowmelt will shift several weeks earlier in the season. Further, there has been evidence that mountain ranges could face more accelerated warming (elevational dependent warming) due to climate change. These future trends have largely been derived from global climate models (CMIP5) which can't resolve some of the more relatively narrow mountain ranges, like the California Sierra Nevada, in great detail. Therefore, due to the importance of orographic uplift on weather fronts, eastern Pacific sea-surface temperature anomalies, atmospheric river events, and mesoscale convective systems, high-resolution global scale modeling techniques are necessary to properly resolve western USA mountain range climatology. Variable-resolution global climate models (VRGCMs) are a promising next-generation technique to analyze both past and future hydroclimatic trends in the region. VRGCMs serve as a bridge between regional and global models by allowing for high-resolution in areas of interest, eliminate lateral boundary forcings (and resultant model biases), allow for more dynamically inclusive large-scale climate teleconnections, and require smaller simulation times and lower data storage demand (compared to conventional global models). This presentation focuses on validating these next-generation models as well as projecting future climate change scenario impacts on several of the western USA's key hydroclimate metrics (e.g., two-meter surface temperature, snow cover, snow water equivalent, and snowfall) to inform water managers and policy makers and offer resilience to climate change impacts facing the region.

An evaluation of the variable-resolution CESM for modeling California's climate: Evaluation of VR-CESM for Modeling California's Climate

DOE Office of Scientific and Technical Information (OSTI.GOV)

Huang, Xingying; Rhoades, Alan M.; Ullrich, Paul A.

In this paper, the recently developed variable-resolution option within the Community Earth System Model (VR-CESM) is assessed for long-term regional climate modeling of California at 0.25° (~ 28 km) and 0.125° (~ 14 km) horizontal resolutions. The mean climatology of near-surface temperature and precipitation is analyzed and contrasted with reanalysis, gridded observational data sets, and a traditional regional climate model (RCM)—the Weather Research and Forecasting (WRF) model. Statistical metrics for model evaluation and tests for differential significance have been extensively applied. VR-CESM tended to produce a warmer summer (by about 1–3°C) and overestimated overall winter precipitation (about 25%–35%) compared tomore » reference data sets when sea surface temperatures were prescribed. Increasing resolution from 0.25° to 0.125° did not produce a statistically significant improvement in the model results. By comparison, the analogous WRF climatology (constrained laterally and at the sea surface by ERA-Interim reanalysis) was ~1–3°C colder than the reference data sets, underestimated precipitation by ~20%–30% at 27 km resolution, and overestimated precipitation by ~ 65–85% at 9 km. Overall, VR-CESM produced comparable statistical biases to WRF in key climatological quantities. Moreover, this assessment highlights the value of variable-resolution global climate models (VRGCMs) in capturing fine-scale atmospheric processes, projecting future regional climate, and addressing the computational expense of uniform-resolution global climate models.« less

An evaluation of the variable-resolution CESM for modeling California's climate: Evaluation of VR-CESM for Modeling California's Climate

DOE PAGES

Huang, Xingying; Rhoades, Alan M.; Ullrich, Paul A.; ...

2016-03-01

In this paper, the recently developed variable-resolution option within the Community Earth System Model (VR-CESM) is assessed for long-term regional climate modeling of California at 0.25° (~ 28 km) and 0.125° (~ 14 km) horizontal resolutions. The mean climatology of near-surface temperature and precipitation is analyzed and contrasted with reanalysis, gridded observational data sets, and a traditional regional climate model (RCM)—the Weather Research and Forecasting (WRF) model. Statistical metrics for model evaluation and tests for differential significance have been extensively applied. VR-CESM tended to produce a warmer summer (by about 1–3°C) and overestimated overall winter precipitation (about 25%–35%) compared tomore » reference data sets when sea surface temperatures were prescribed. Increasing resolution from 0.25° to 0.125° did not produce a statistically significant improvement in the model results. By comparison, the analogous WRF climatology (constrained laterally and at the sea surface by ERA-Interim reanalysis) was ~1–3°C colder than the reference data sets, underestimated precipitation by ~20%–30% at 27 km resolution, and overestimated precipitation by ~ 65–85% at 9 km. Overall, VR-CESM produced comparable statistical biases to WRF in key climatological quantities. Moreover, this assessment highlights the value of variable-resolution global climate models (VRGCMs) in capturing fine-scale atmospheric processes, projecting future regional climate, and addressing the computational expense of uniform-resolution global climate models.« less

Estimation of Global 1km-grid Terrestrial Carbon Exchange Part I: Developing Inputs and Modelling

NASA Astrophysics Data System (ADS)

Sasai, T.; Murakami, K.; Kato, S.; Matsunaga, T.; Saigusa, N.; Hiraki, K.

2015-12-01

Global terrestrial carbon cycle largely depends on a spatial pattern in land cover type, which is heterogeneously-distributed over regional and global scales. However, most studies, which aimed at the estimation of carbon exchanges between ecosystem and atmosphere, remained within several tens of kilometers grid spatial resolution, and the results have not been enough to understand the detailed pattern of carbon exchanges based on ecological community. Improving the sophistication of spatial resolution is obviously necessary to enhance the accuracy of carbon exchanges. Moreover, the improvement may contribute to global warming awareness, policy makers and other social activities. In this study, we show global terrestrial carbon exchanges (net ecosystem production, net primary production, and gross primary production) with 1km-grid resolution. As methodology for computing the exchanges, we 1) developed a global 1km-grid climate and satellite dataset based on the approach in Setoyama and Sasai (2013); 2) used the satellite-driven biosphere model (Biosphere model integrating Eco-physiological And Mechanistic approaches using Satellite data: BEAMS) (Sasai et al., 2005, 2007, 2011); 3) simulated the carbon exchanges by using the new dataset and BEAMS by the use of a supercomputer that includes 1280 CPU and 320 GPGPU cores (GOSAT RCF of NIES). As a result, we could develop a global uniform system for realistically estimating terrestrial carbon exchange, and evaluate net ecosystem production in each community level; leading to obtain highly detailed understanding of terrestrial carbon exchanges.

A high-resolution global-scale groundwater model

NASA Astrophysics Data System (ADS)

de Graaf, I. E. M.; Sutanudjaja, E. H.; van Beek, L. P. H.; Bierkens, M. F. P.

2015-02-01

Groundwater is the world's largest accessible source of fresh water. It plays a vital role in satisfying basic needs for drinking water, agriculture and industrial activities. During times of drought groundwater sustains baseflow to rivers and wetlands, thereby supporting ecosystems. Most global-scale hydrological models (GHMs) do not include a groundwater flow component, mainly due to lack of geohydrological data at the global scale. For the simulation of lateral flow and groundwater head dynamics, a realistic physical representation of the groundwater system is needed, especially for GHMs that run at finer resolutions. In this study we present a global-scale groundwater model (run at 6' resolution) using MODFLOW to construct an equilibrium water table at its natural state as the result of long-term climatic forcing. The used aquifer schematization and properties are based on available global data sets of lithology and transmissivities combined with the estimated thickness of an upper, unconfined aquifer. This model is forced with outputs from the land-surface PCRaster Global Water Balance (PCR-GLOBWB) model, specifically net recharge and surface water levels. A sensitivity analysis, in which the model was run with various parameter settings, showed that variation in saturated conductivity has the largest impact on the groundwater levels simulated. Validation with observed groundwater heads showed that groundwater heads are reasonably well simulated for many regions of the world, especially for sediment basins (R2 = 0.95). The simulated regional-scale groundwater patterns and flow paths demonstrate the relevance of lateral groundwater flow in GHMs. Inter-basin groundwater flows can be a significant part of a basin's water budget and help to sustain river baseflows, especially during droughts. Also, water availability of larger aquifer systems can be positively affected by additional recharge from inter-basin groundwater flows.

Xanthos – A Global Hydrologic Model

DOE Office of Scientific and Technical Information (OSTI.GOV)

Li, Xinya; Vernon, Chris R.; Hejazi, Mohamad I.

Xanthos is an open-source hydrologic model, written in Python, designed to quantify and analyse global water availability. Xanthos simulates historical and future global water availability on a monthly time step at a spatial resolution of 0.5 geographic degrees. Xanthos was designed to be extensible and used by scientists that study global water supply and work with the Global Change Assessment Model (GCAM). Xanthos uses a user-defined configuration file to specify model inputs, outputs and parameters. Xanthos has been tested using actual global data sets and the model is able to provide historical observations and future estimates of renewable freshwater resourcesmore » in the form of total runoff.« less

Xanthos – A Global Hydrologic Model

DOE PAGES

Li, Xinya; Vernon, Chris R.; Hejazi, Mohamad I.; ...

2017-09-11

Xanthos is an open-source hydrologic model, written in Python, designed to quantify and analyse global water availability. Xanthos simulates historical and future global water availability on a monthly time step at a spatial resolution of 0.5 geographic degrees. Xanthos was designed to be extensible and used by scientists that study global water supply and work with the Global Change Assessment Model (GCAM). Xanthos uses a user-defined configuration file to specify model inputs, outputs and parameters. Xanthos has been tested using actual global data sets and the model is able to provide historical observations and future estimates of renewable freshwater resourcesmore » in the form of total runoff.« less

Technical Challenges and Solutions in Representing Lakes when using WRF in Downscaling Applications

EPA Science Inventory

The Weather Research and Forecasting (WRF) model is commonly used to make high resolution future projections of regional climate by downscaling global climate model (GCM) outputs. Because the GCM fields are typically at a much coarser spatial resolution than the target regional ...

Hydrologic downscaling of soil moisture using global data without site-specific calibration

USDA-ARS?s Scientific Manuscript database

Numerous applications require fine-resolution (10-30 m) soil moisture patterns, but most satellite remote sensing and land-surface models provide coarse-resolution (9-60 km) soil moisture estimates. The Equilibrium Moisture from Topography, Vegetation, and Soil (EMT+VS) model downscales soil moistu...

Lightning Forcing in Global Fire Models: The Importance of Temporal Resolution

NASA Astrophysics Data System (ADS)

Felsberg, A.; Kloster, S.; Wilkenskjeld, S.; Krause, A.; Lasslop, G.

2018-01-01

In global fire models, lightning is typically prescribed from observational data with monthly mean temporal resolution while meteorological forcings, such as precipitation or temperature, are prescribed in a daily resolution. In this study, we investigate the importance of the temporal resolution of the lightning forcing for the simulation of burned area by varying from daily to monthly and annual mean forcing. For this, we utilize the vegetation fire model JSBACH-SPITFIRE to simulate burned area, forced with meteorological and lightning data derived from the general circulation model ECHAM6. On a global scale, differences in burned area caused by lightning forcing applied in coarser temporal resolution stay below 0.55% compared to the use of daily mean forcing. Regionally, however, differences reach up to 100%, depending on the region and season. Monthly averaged lightning forcing as well as the monthly lightning climatology cause differences through an interaction between lightning ignitions and fire prone weather conditions, accounted for by the fire danger index. This interaction leads to decreased burned area in the boreal zone and increased burned area in the Tropics and Subtropics under the coarser temporal resolution. The exclusion of interannual variability, when forced with the lightning climatology, has only a minor impact on the simulated burned area. Annually averaged lightning forcing causes differences as a direct result of the eliminated seasonal characteristics of lightning. Burned area is decreased in summer and increased in winter where fuel is available. Regions with little seasonality, such as the Tropics and Subtropics, experience an increase in burned area.

Evaluation of a Mesoscale Convective System in Variable-Resolution CESM

NASA Astrophysics Data System (ADS)

Payne, A. E.; Jablonowski, C.

2017-12-01

Warm season precipitation over the Southern Great Plains (SGP) follows a well observed diurnal pattern of variability, peaking at night-time, due to the eastward propagation of mesoscale convection systems that develop over the eastern slopes of the Rockies in the late afternoon. While most climate models are unable to adequately capture the organization of convection and characteristic pattern of precipitation over this region, models with high enough resolution to explicitly resolve convection show improvement. However, high resolution simulations are computationally expensive and, in the case of regional climate models, are subject to boundary conditions. Newly developed variable resolution global climate models strike a balance between the benefits of high-resolution regional climate models and the large-scale dynamics of global climate models and low computational cost. Recently developed parameterizations that are insensitive to the model grid scale provide a way to improve model performance. Here, we present an evaluation of the newly available Cloud Layers Unified by Binormals (CLUBB) parameterization scheme in a suite of variable-resolution CESM simulations with resolutions ranging from 110 km to 7 km within a regionally refined region centered over the SGP Atmospheric Radiation Measurement (ARM) site. Simulations utilize the hindcast approach developed by the Department of Energy's Cloud-Associated Parameterizations Testbed (CAPT) for the assessment of climate models. We limit our evaluation to a single mesoscale convective system that passed over the region on May 24, 2008. The effects of grid-resolution on the timing and intensity of precipitation, as well as, on the transition from shallow to deep convection are assessed against ground-based observations from the SGP ARM site, satellite observations and ERA-Interim reanalysis.

A Synopsis of Global Mapping of Freshwater Habitats and Biodiversity: Implications for Conservation

DOE Office of Scientific and Technical Information (OSTI.GOV)

McManamay, Ryan A.; Griffiths, Natalie A.; DeRolph, Christopher R.

Accurately mapping freshwater habitats and biodiversity at high-resolutions across the globe is essential for assessing the vulnerability and threats to freshwater organisms and prioritizing conservation efforts. Since the 2000s, extensive efforts have been devoted to mapping global freshwater habitats (rivers, lakes, and wetlands), the spatial representation of which has changed dramatically over time with new geospatial data products and improved remote sensing technologies. Some of these mapping efforts, however, are still coarse representations of actual conditions. Likewise, the resolution and scope of global freshwater biodiversity compilation efforts have also increased, but are yet to mirror the spatial resolution and fidelitymore » of mapped freshwater environments. In our synopsis, we find that efforts to map freshwater habitats have been conducted independently of those for freshwater biodiversity; subsequently, there is little congruence in the spatial representation and resolution of the two efforts. We suggest that global species distribution models are needed to fill this information gap; however, limiting data on habitat characteristics at scales that complement freshwater habitats has prohibited global high-resolution biogeography efforts. Emerging research trends, such as mapping habitat alteration in freshwater ecosystems and trait biogeography, show great promise in mechanistically linking global anthropogenic stressors to freshwater biodiversity decline and extinction risk.« less

Missing pieces of the puzzle: understanding decadal variability of Sahel Rainfall

NASA Astrophysics Data System (ADS)

Vellinga, Michael; Roberts, Malcolm; Vidale, Pier-Luigi; Mizielinski, Matthew; Demory, Marie-Estelle; Schiemann, Reinhard; Strachan, Jane; Bain, Caroline

2015-04-01

The instrumental record shows that substantial decadal fluctuations affected Sahel rainfall from the West African monsoon throughout the 20th century. Climate models generally underestimate the magnitude of decadal Sahel rainfall changes compared to observations. This shows that the processes that control low-frequency Sahel rainfall change are misrepresented in most CMIP5-era climate models. Reliable climate information of future low-frequency rainfall changes thus remains elusive. Here we identify key processes that control the magnitude of the decadal rainfall recovery in the Sahel since the mid-1980s. We show its sensitivity to model resolution and physics in a suite of experiments with global HadGEM3 model configurations at resolutions between 130-25 km. The decadal rainfall trend increases with resolution and at 60-25 km falls within the observed range. Higher resolution models have stronger increases of moisture supply and of African Easterly wave activity. Easterly waves control the occurrence of strong organised rainfall events which carry most of the decadal trend. Weak rainfall events occur too frequently at all resolutions and at low resolution contribute substantially to the decadal trend. All of this behaviour is seen across CMIP5, including future scenarios. Additional simulations with a global 12km version of HadGEM3 show that treating convection explicitly dramatically improves the properties of Sahel rainfall systems. We conclude that interaction between convective scale and global scale processes is key to decadal rainfall changes in the Sahel. This work is distributed under the Creative Commons Attribution 3.0 Unported License together with an author copyright. This license does not conflict with the regulations of the Crown Copyright.Crown Copyright

Tropospheric carbon monoxide over the Pacific during HIPPO: two-way coupled simulation of GEOS-Chem and its multiple nested models

NASA Astrophysics Data System (ADS)

Yan, Y.-Y.; Lin, J.-T.; Kuang, Y.; Yang, D.; Zhang, L.

2014-07-01

Global chemical transport models (CTMs) are used extensively to study air pollution and transport at a global scale. These models are limited by coarse horizontal resolutions, not allowing for detailed representation of small-scale nonlinear processes over the pollutant source regions. Here we couple the global GEOS-Chem CTM and its three high-resolution nested models to simulate the tropospheric carbon monoxide (CO) over the Pacific Ocean during five HIAPER Pole-to-Pole Observations (HIPPO) campaigns between 2009 and 2011. We develop a two-way coupler, PKUCPL, to integrate simulation results for chemical constituents from the global model (at 2.5° long. × 2° lat.) and the three nested models (at 0.667° long. × 0.5° lat.) covering Asia, North America and Europe, respectively. The coupler obtains nested model results to modify the global model simulation within the respective nested domains, and simultaneously acquires global model results to provide lateral boundary conditions for the nested models. Compared to the global model alone, the two-way coupled simulation results in enhanced CO concentrations in the nested domains. Sensitivity tests suggest the enhancement to be a result of improved representation of the spatial distributions of CO, nitrogen oxides and non-methane volatile organic compounds, the meteorological dependence of natural emissions, and other resolution-dependent processes. The relatively long lifetime of CO allows for the enhancement to be accumulated and carried across the globe. We find that the two-way coupled simulation increases the global tropospheric mean CO concentrations in 2009 by 10.4%, with a greater enhancement at 13.3% in the Northern Hemisphere. Coincidently, the global tropospheric mean hydroxyl radical (OH) is reduced by 4.2% (as compared to the interannual variability of OH at 2.3%), resulting in a 4.2% enhancement in the methyl chloroform lifetime (MCF, via reaction with tropospheric OH). The resulting CO and OH contents and MCF lifetime are closer to observation-based estimates. Both the global and the two-way coupled models capture the general spatiotemporal patterns of HIPPO CO over the Pacific. The two-way coupled simulation is much closer to HIPPO CO, with a mean bias of 1.1 ppb (1.4%) below 9 km compared to the bias at -7.2 ppb (-9.2%) for the global model. The improvement is most apparent over the North Pacific. Our test simulations show that the global model could resemble the two-way coupled simulation (especially below 4 km) by increasing its global CO emissions by 15% for HIPPO-1 and HIPPO-3, by 25% for HIPPO-2 and HIPPO-4, and by 35% for HIPPO-5. This has important implications for using the global model to constrain CO emissions. Thus, the two-way coupled simulation is a significantly improved model tool to studying the global impacts of air pollutants from major anthropogenic source regions.

Development of JPSS VIIRS Global Gridded Vegetation Index products for NOAA NCEP Environmental Modeling Systems

NASA Astrophysics Data System (ADS)

Vargas, Marco; Miura, Tomoaki; Csiszar, Ivan; Zheng, Weizhong; Wu, Yihua; Ek, Michael

2017-04-01

The first Joint Polar Satellite System (JPSS) mission, the Suomi National Polar-orbiting Partnership (S-NPP) satellite, was successfully launched in October, 2011, and it will be followed by JPSS-1, slated for launch in 2017. JPSS provides operational continuity of satellite-based observations and products for NOAA's Polar Operational Environmental Satellites (POES). Vegetation products derived from satellite measurements are used for weather forecasting, land modeling, climate research, and monitoring the environment including drought, the health of ecosystems, crop monitoring and forest fires. The operationally produced S-NPP VIIRS Vegetation Index (VI) Environmental Data Record (EDR) includes two vegetation indices: the Top of the Atmosphere (TOA) Normalized Difference Vegetation Index (NDVI), and the Top of the Canopy (TOC) Enhanced Vegetation Index (EVI). For JPSS-1, the S-NPP Vegetation Index EDR algorithm has been updated to include the TOC NDV. The current JPSS operational VI products are generated in granule style at 375 meter resolution at nadir, but these products in granule format cannot be ingested into NOAA operational monitoring and decision making systems. For that reason, the NOAA JPSS Land Team is developing a new global gridded Vegetation Index (VI) product suite for operational use by the NOAA National Centers for Environmental Prediction (NCEP). The new global gridded VIs will be used in the Multi-Physics (MP) version of the Noah land surface model (Noah-MP) in NCEP NOAA Environmental Modeling System (NEMS) for plant growth and data assimilation and to describe vegetation coverage and density in order to model the correct surface energy partition. The new VI 4km resolution global gridded products (TOA NDVI, TOC NDVI and TOC EVI) are being designed to meet the needs of directly ingesting vegetation index variables without the need to develop local gridding and compositing procedures. These VI products will be consistent with the already operational SNPP VIIRS Green Vegetation Fraction (GVF) global gridded 4km resolution. The ultimate goal is a global consistent set of global gridded land products at 1-km resolution to enable consistent use of the products in the full suite of global and regional NCEP land models. The new JPSS vegetation products system is scheduled to transition to operations in the fall of 2017.

A high-resolution physically-based global flood hazard map

NASA Astrophysics Data System (ADS)

Kaheil, Y.; Begnudelli, L.; McCollum, J.

2016-12-01

We present the results from a physically-based global flood hazard model. The model uses a physically-based hydrologic model to simulate river discharges, and 2D hydrodynamic model to simulate inundation. The model is set up such that it allows the application of large-scale flood hazard through efficient use of parallel computing. For hydrology, we use the Hillslope River Routing (HRR) model. HRR accounts for surface hydrology using Green-Ampt parameterization. The model is calibrated against observed discharge data from the Global Runoff Data Centre (GRDC) network, among other publicly-available datasets. The parallel-computing framework takes advantage of the river network structure to minimize cross-processor messages, and thus significantly increases computational efficiency. For inundation, we implemented a computationally-efficient 2D finite-volume model with wetting/drying. The approach consists of simulating flood along the river network by forcing the hydraulic model with the streamflow hydrographs simulated by HRR, and scaled up to certain return levels, e.g. 100 years. The model is distributed such that each available processor takes the next simulation. Given an approximate criterion, the simulations are ordered from most-demanding to least-demanding to ensure that all processors finalize almost simultaneously. Upon completing all simulations, the maximum envelope of flood depth is taken to generate the final map. The model is applied globally, with selected results shown from different continents and regions. The maps shown depict flood depth and extent at different return periods. These maps, which are currently available at 3 arc-sec resolution ( 90m) can be made available at higher resolutions where high resolution DEMs are available. The maps can be utilized by flood risk managers at the national, regional, and even local levels to further understand their flood risk exposure, exercise certain measures of mitigation, and/or transfer the residual risk financially through flood insurance programs.

Tropical cyclones in a stabilized 1.5 and 2 degree warmer world.

NASA Astrophysics Data System (ADS)

Wehner, M. F.; Stone, D. A.; Loring, B.; Krishnan, H.

2017-12-01

We present an ensemble of very high resolution global climate model simulations of a stabilized 1.5oC and 2oC warmer climate as envisioned by the Paris COP21 agreement. The resolution of this global climate model (25km) permits simulated tropical cyclones up to Category Five on the Saffir-Simpson scale Projected changes in tropical cyclones are significant. Tropical cyclones in the two stabilization scenarios are less frequent but more intense than in simulations of the present. Output data from these simulations is freely available to all interested parties and should prove a useful resource to those interested in studying the impacts of stabilized global warming.

Tropical Cyclones, Hurricanes, and Climate: NASA's Global Cloud-Scale Simulations and New Observations that Characterize the Lifecycle of Hurricanes

NASA Technical Reports Server (NTRS)

Putman, William M.

2010-01-01

One of the primary interests of Global Change research is the impact of climate changes and climate variability on extreme weather events, such as intense tropical storms and hurricanes. Atmospheric climate models run at resolutions of global weather models have been used to study the impact of climate variability, as seen in sea surface temperatures, on the frequency and intensity of tropical cyclones. NASA's Goddard Earth Observing System Model, version 5 (GEOS-5) in ensembles run at 50 km resolution has been able to reproduce the interannual variations of tropical cyclone frequency seen in nature. This, and other global models, have found it much more difficult to reproduce the interannual changes in intensity, a result that reflects the inability of the models to simulate the intensities of the most extreme storms. Better representation of the structures of cyclones requires much higher resolution models. Such improved representation is also fundamental to making best use of satellite observations. In collaboration with NOAA's Geophysical Fluid Dynamics Laboratory, GEOS-5 now has the capability of running at much higher resolution to better represent cloud-scale resolutions. Global simulations at cloud-permitting resolutions (10- to 3.5-km) allows for the development of realistic tropical cyclones from tropical storm 119 km/hr winds) to category 5 (>249km1hr winds) intensities. GEOS-5 has produced realistic rain-band and eye-wall structures in tropical cyclones that can be directly analyzed against satellite observations. For the first time a global climate model is capable of representing realistic intensity and track variability on a seasonal scale across basins. GEOS-5 is also used in assimilation mode to test the impact of NASA's observations on tropical cyclone forecasts. One such test, for tropical cyclone Nargis in the Indian Ocean in May 2008, showed that observations from Atmospheric Infrared Sounder (AIRS) and the Advanced Microwave Sounding Unit (AMSU-A) on Aqua substantially reduced forecast track errors. Tropical cyclones in the northern Indian Ocean pose serious challenges to operational weather forecasting systems, partly due to their shorter lifespan and more erratic track, compared to those in the Atlantic and the Pacific. SA is also bringing several state of the art instruments in recent field campaigns to peer under the clouds and study the inner workings of the tropical storms. With the Genesis and Rapid Intensification Processes (GRIP) experiment, a NASA Earth science field experiment in 2010 that includes the Global Hawk Unmanned Airborne System (UAS) configured with a suite of in situ and remote sensing instruments that are observing and characterizing the lifecycle of hurricanes, we expect significant improvement in our understanding of the track and intensification processes with the assimilation of the satellite and field campaign observations of meteorological parameters in the numerical prediction models.

How does increasing horizontal resolution in a global climate model improve the simulation of aerosol-cloud interactions?

DOE PAGES

Ma, Po-Lun; Rasch, Philip J.; Wang, Minghuai; ...

2015-06-23

We report the Community Atmosphere Model Version 5 is run at horizontal grid spacing of 2, 1, 0.5, and 0.25°, with the meteorology nudged toward the Year Of Tropical Convection analysis, and cloud simulators and the collocated A-Train satellite observations are used to explore the resolution dependence of aerosol-cloud interactions. The higher-resolution model produces results that agree better with observations, showing an increase of susceptibility of cloud droplet size, indicating a stronger first aerosol indirect forcing (AIF), and a decrease of susceptibility of precipitation probability, suggesting a weaker second AIF. The resolution sensitivities of AIF are attributed to those ofmore » droplet nucleation and precipitation parameterizations. Finally, the annual average AIF in the Northern Hemisphere midlatitudes (where most anthropogenic emissions occur) in the 0.25° model is reduced by about 1 W m -2 (-30%) compared to the 2° model, leading to a 0.26 W m -2 reduction (-15%) in the global annual average AIF.« less

Toward GEOS-6, A Global Cloud System Resolving Atmospheric Model

NASA Technical Reports Server (NTRS)

Putman, William M.

2010-01-01

NASA is committed to observing and understanding the weather and climate of our home planet through the use of multi-scale modeling systems and space-based observations. Global climate models have evolved to take advantage of the influx of multi- and many-core computing technologies and the availability of large clusters of multi-core microprocessors. GEOS-6 is a next-generation cloud system resolving atmospheric model that will place NASA at the forefront of scientific exploration of our atmosphere and climate. Model simulations with GEOS-6 will produce a realistic representation of our atmosphere on the scale of typical satellite observations, bringing a visual comprehension of model results to a new level among the climate enthusiasts. In preparation for GEOS-6, the agency's flagship Earth System Modeling Framework [JDl] has been enhanced to support cutting-edge high-resolution global climate and weather simulations. Improvements include a cubed-sphere grid that exposes parallelism; a non-hydrostatic finite volume dynamical core, and algorithm designed for co-processor technologies, among others. GEOS-6 represents a fundamental advancement in the capability of global Earth system models. The ability to directly compare global simulations at the resolution of spaceborne satellite images will lead to algorithm improvements and better utilization of space-based observations within the GOES data assimilation system

Will high-resolution global ocean models benefit coupled predictions on short-range to climate timescales?

NASA Astrophysics Data System (ADS)

Hewitt, Helene T.; Bell, Michael J.; Chassignet, Eric P.; Czaja, Arnaud; Ferreira, David; Griffies, Stephen M.; Hyder, Pat; McClean, Julie L.; New, Adrian L.; Roberts, Malcolm J.

2017-12-01

As the importance of the ocean in the weather and climate system is increasingly recognised, operational systems are now moving towards coupled prediction not only for seasonal to climate timescales but also for short-range forecasts. A three-way tension exists between the allocation of computing resources to refine model resolution, the expansion of model complexity/capability, and the increase of ensemble size. Here we review evidence for the benefits of increased ocean resolution in global coupled models, where the ocean component explicitly represents transient mesoscale eddies and narrow boundary currents. We consider lessons learned from forced ocean/sea-ice simulations; from studies concerning the SST resolution required to impact atmospheric simulations; and from coupled predictions. Impacts of the mesoscale ocean in western boundary current regions on the large-scale atmospheric state have been identified. Understanding of air-sea feedback in western boundary currents is modifying our view of the dynamics in these key regions. It remains unclear whether variability associated with open ocean mesoscale eddies is equally important to the large-scale atmospheric state. We include a discussion of what processes can presently be parameterised in coupled models with coarse resolution non-eddying ocean models, and where parameterizations may fall short. We discuss the benefits of resolution and identify gaps in the current literature that leave important questions unanswered.

Outcomes and challenges of global high-resolution non-hydrostatic atmospheric simulations using the K computer

NASA Astrophysics Data System (ADS)

Satoh, Masaki; Tomita, Hirofumi; Yashiro, Hisashi; Kajikawa, Yoshiyuki; Miyamoto, Yoshiaki; Yamaura, Tsuyoshi; Miyakawa, Tomoki; Nakano, Masuo; Kodama, Chihiro; Noda, Akira T.; Nasuno, Tomoe; Yamada, Yohei; Fukutomi, Yoshiki

2017-12-01

This article reviews the major outcomes of a 5-year (2011-2016) project using the K computer to perform global numerical atmospheric simulations based on the non-hydrostatic icosahedral atmospheric model (NICAM). The K computer was made available to the public in September 2012 and was used as a primary resource for Japan's Strategic Programs for Innovative Research (SPIRE), an initiative to investigate five strategic research areas; the NICAM project fell under the research area of climate and weather simulation sciences. Combining NICAM with high-performance computing has created new opportunities in three areas of research: (1) higher resolution global simulations that produce more realistic representations of convective systems, (2) multi-member ensemble simulations that are able to perform extended-range forecasts 10-30 days in advance, and (3) multi-decadal simulations for climatology and variability. Before the K computer era, NICAM was used to demonstrate realistic simulations of intra-seasonal oscillations including the Madden-Julian oscillation (MJO), merely as a case study approach. Thanks to the big leap in computational performance of the K computer, we could greatly increase the number of cases of MJO events for numerical simulations, in addition to integrating time and horizontal resolution. We conclude that the high-resolution global non-hydrostatic model, as used in this five-year project, improves the ability to forecast intra-seasonal oscillations and associated tropical cyclogenesis compared with that of the relatively coarser operational models currently in use. The impacts of the sub-kilometer resolution simulation and the multi-decadal simulations using NICAM are also reviewed.

Using Global Plate Velocity Boundary Conditions for Embedded Regional Geodynamic Models

NASA Astrophysics Data System (ADS)

Taramon Gomez, Jorge; Morgan, Jason; Perez-Gussinye, Marta

2015-04-01

The treatment of far-field boundary conditions is one of the most poorly resolved issues for regional modeling of geodynamic processes. In viscous flow, the choice of far-field boundary conditions often strongly shapes the large-scale structure of a geosimulation. The mantle velocity field along the sidewalls and base of a modeling region is typically much more poorly known than the geometry of past global motions of the surface plates as constrained by global plate motion reconstructions. For regional rifting models it has become routine to apply highly simplified 'plate spreading' or 'uniform rifting' boundary conditions to a 3-D model that limits its ability to simulate the geodynamic evolution of a specific rifted margin. One way researchers are exploring the sensitivity of regional models to uncertain boundary conditions is to use a nested modeling approach in which a global model is used to determine a large-scale flow pattern that is imposed as a constraint along the boundaries of the region to be modeled. Here we explore the utility of a different approach that takes advantage of the ability of finite element models to use unstructured meshes than can embed much higher resolution sub-regions within a spherical global mesh. In our initial project to validate this approach, we create a global spherical mesh in which a higher resolution sub-region is created around the nascent South Atlantic Rifting Margin. Global Plate motion BCs and plate boundaries are applied for the time of the onset of rifting, continuing through several 10s of Ma of rifting. Thermal, compositional, and melt-related buoyancy forces are only non-zero within the high-resolution subregion, elsewhere, motions are constrained by surface plate-motion constraints. The total number of unknowns needed to solve an embedded regional model with this approach is less than 1/3 larger than that needed for a structured-mesh solution on a Cartesian or spherical cap sub-regional mesh. Here we illustrate the initial steps within this workflow for creating time-varying surface boundary conditions (using GPlates), and a time-variable unstructured 3-D spherical mesh.

Scaling a Convection-Resolving RCM to Near-Global Scales

NASA Astrophysics Data System (ADS)

Leutwyler, D.; Fuhrer, O.; Chadha, T.; Kwasniewski, G.; Hoefler, T.; Lapillonne, X.; Lüthi, D.; Osuna, C.; Schar, C.; Schulthess, T. C.; Vogt, H.

2017-12-01

In the recent years, first decade-long kilometer-scale resolution RCM simulations have been performed on continental-scale computational domains. However, the size of the planet Earth is still an order of magnitude larger and thus the computational implications of performing global climate simulations at this resolution are challenging. We explore the gap between the currently established RCM simulations and global simulations by scaling the GPU accelerated version of the COSMO model to a near-global computational domain. To this end, the evolution of an idealized moist baroclinic wave has been simulated over the course of 10 days with a grid spacing of up to 930 m. The computational mesh employs 36'000 x 16'001 x 60 grid points and covers 98.4% of the planet's surface. The code shows perfect weak scaling up to 4'888 Nodes of the Piz Daint supercomputer and yields 0.043 simulated years per day (SYPD) which is approximately one seventh of the 0.2-0.3 SYPD required to conduct AMIP-type simulations. However, at half the resolution (1.9 km) we've observed 0.23 SYPD. Besides formation of frontal precipitating systems containing embedded explicitly-resolved convective motions, the simulations reveal a secondary instability that leads to cut-off warm-core cyclonic vortices in the cyclone's core, once the grid spacing is refined to the kilometer scale. The explicit representation of embedded moist convection and the representation of the previously unresolved instabilities exhibit a physically different behavior in comparison to coarser-resolution simulations. The study demonstrates that global climate simulations using kilometer-scale resolution are imminent and serves as a baseline benchmark for global climate model applications and future exascale supercomputing systems.

NASA Downscaling Project: Final Report

NASA Technical Reports Server (NTRS)

Ferraro, Robert; Waliser, Duane; Peters-Lidard, Christa

2017-01-01

A team of researchers from NASA Ames Research Center, Goddard Space Flight Center, the Jet Propulsion Laboratory, and Marshall Space Flight Center, along with university partners at UCLA, conducted an investigation to explore whether downscaling coarse resolution global climate model (GCM) predictions might provide valid insights into the regional impacts sought by decision makers. Since the computational cost of running global models at high spatial resolution for any useful climate scale period is prohibitive, the hope for downscaling is that a coarse resolution GCM provides sufficiently accurate synoptic scale information for a regional climate model (RCM) to accurately develop fine scale features that represent the regional impacts of a changing climate. As a proxy for a prognostic climate forecast model, and so that ground truth in the form of satellite and in-situ observations could be used for evaluation, the MERRA and MERRA - 2 reanalyses were used to drive the NU - WRF regional climate model and a GEOS - 5 replay. This was performed at various resolutions that were at factors of 2 to 10 higher than the reanalysis forcing. A number of experiments were conducted that varied resolution, model parameterizations, and intermediate scale nudging, for simulations over the continental US during the period from 2000 - 2010. The results of these experiments were compared to observational datasets to evaluate the output.

NASA Downscaling Project

NASA Technical Reports Server (NTRS)

Ferraro, Robert; Waliser, Duane; Peters-Lidard, Christa

2017-01-01

A team of researchers from NASA Ames Research Center, Goddard Space Flight Center, the Jet Propulsion Laboratory, and Marshall Space Flight Center, along with university partners at UCLA, conducted an investigation to explore whether downscaling coarse resolution global climate model (GCM) predictions might provide valid insights into the regional impacts sought by decision makers. Since the computational cost of running global models at high spatial resolution for any useful climate scale period is prohibitive, the hope for downscaling is that a coarse resolution GCM provides sufficiently accurate synoptic scale information for a regional climate model (RCM) to accurately develop fine scale features that represent the regional impacts of a changing climate. As a proxy for a prognostic climate forecast model, and so that ground truth in the form of satellite and in-situ observations could be used for evaluation, the MERRA and MERRA-2 reanalyses were used to drive the NU-WRF regional climate model and a GEOS-5 replay. This was performed at various resolutions that were at factors of 2 to 10 higher than the reanalysis forcing. A number of experiments were conducted that varied resolution, model parameterizations, and intermediate scale nudging, for simulations over the continental US during the period from 2000-2010. The results of these experiments were compared to observational datasets to evaluate the output.

Impacts of spatial resolution and representation of flow connectivity on large-scale simulation of floods

NASA Astrophysics Data System (ADS)

Mateo, Cherry May R.; Yamazaki, Dai; Kim, Hyungjun; Champathong, Adisorn; Vaze, Jai; Oki, Taikan

2017-10-01

Global-scale river models (GRMs) are core tools for providing consistent estimates of global flood hazard, especially in data-scarce regions. Due to former limitations in computational power and input datasets, most GRMs have been developed to use simplified representations of flow physics and run at coarse spatial resolutions. With increasing computational power and improved datasets, the application of GRMs to finer resolutions is becoming a reality. To support development in this direction, the suitability of GRMs for application to finer resolutions needs to be assessed. This study investigates the impacts of spatial resolution and flow connectivity representation on the predictive capability of a GRM, CaMa-Flood, in simulating the 2011 extreme flood in Thailand. Analyses show that when single downstream connectivity (SDC) is assumed, simulation results deteriorate with finer spatial resolution; Nash-Sutcliffe efficiency coefficients decreased by more than 50 % between simulation results at 10 km resolution and 1 km resolution. When multiple downstream connectivity (MDC) is represented, simulation results slightly improve with finer spatial resolution. The SDC simulations result in excessive backflows on very flat floodplains due to the restrictive flow directions at finer resolutions. MDC channels attenuated these effects by maintaining flow connectivity and flow capacity between floodplains in varying spatial resolutions. While a regional-scale flood was chosen as a test case, these findings should be universal and may have significant impacts on large- to global-scale simulations, especially in regions where mega deltas exist.These results demonstrate that a GRM can be used for higher resolution simulations of large-scale floods, provided that MDC in rivers and floodplains is adequately represented in the model structure.

GEOS S2S-2_1: GMAO's New High Resolution Seasonal Prediction System

NASA Technical Reports Server (NTRS)

Molod, Andrea; Akella, Santha; Andrews, Lauren; Barahona, Donifan; Borovikov, Anna; Chang, Yehui; Cullather, Richard; Hackert, Eric; Kovach, Robin; Koster, Randal;

Fast Magnetotail Reconnection: Challenge to Global MHD Modeling

NASA Astrophysics Data System (ADS)

Kuznetsova, M. M.; Hesse, M.; Rastaetter, L.; Toth, G.; de Zeeuw, D.; Gombosi, T.

2005-05-01

Representation of fast magnetotail reconnection rates during substorm onset is one of the major challenges to global MHD modeling. Our previous comparative study of collisionless magnetic reconnection in GEM Challenge geometry demonstrated that the reconnection rate is controlled by ion nongyrotropic behavior near the reconnection site and that it can be described in terms of nongyrotropic corrections to the magnetic induction equation. To further test the approach we performed MHD simulations with nongyrotropic corrections of forced reconnection for the Newton Challenge setup. As a next step we employ the global MHD code BATSRUS and test different methods to model fast magnetotail reconnection rates by introducing non-ideal corrections to the induction equation in terms of nongyrotropic corrections, spatially localized resistivity, or current dependent resistivity. The BATSRUS adaptive grid structure allows to perform global simulations with spatial resolution near the reconnection site comparable with spatial resolution of local MHD simulations for the Newton Challenge. We select solar wind conditions which drive the accumulation of magnetic field in the tail lobes and subsequent magnetic reconnection and energy release. Testing the ability of global MHD models to describe magnetotail evolution during substroms is one of the elements of science based validation efforts at the Community Coordinated Modeling Center.

Tropical cyclones over the North Indian Ocean: experiments with the high-resolution global icosahedral grid point model GME

NASA Astrophysics Data System (ADS)

Kumkar, Yogesh V.; Sen, P. N.; Chaudhari, Hemankumar S.; Oh, Jai-Ho

2018-02-01

In this paper, an attempt has been made to conduct a numerical experiment with the high-resolution global model GME to predict the tropical storms in the North Indian Ocean during the year 2007. Numerical integrations using the icosahedral hexagonal grid point global model GME were performed to study the evolution of tropical cyclones, viz., Akash, Gonu, Yemyin and Sidr over North Indian Ocean during 2007. It has been seen that the GME model forecast underestimates cyclone's intensity, but the model can capture the evolution of cyclone's intensity especially its weakening during landfall, which is primarily due to the cutoff of the water vapor supply in the boundary layer as cyclones approach the coastal region. A series of numerical simulation of tropical cyclones have been performed with GME to examine model capability in prediction of intensity and track of the cyclones. The model performance is evaluated by calculating the root mean square errors as cyclone track errors.

On the evaluation of global sea-salt aerosol models at coastal/orographic sites

NASA Astrophysics Data System (ADS)

Spada, M.; Jorba, O.; Pérez García-Pando, C.; Janjic, Z.; Baldasano, J. M.

2015-01-01

Sea-salt aerosol global models are typically evaluated against concentration observations at coastal stations that are unaffected by local surf conditions and thus considered representative of open ocean conditions. Despite recent improvements in sea-salt source functions, studies still show significant model errors in specific regions. Using a multiscale model, we investigated the effect of high model resolution (0.1° × 0.1° vs. 1° × 1.4°) upon sea-salt patterns in four stations from the University of Miami Network: Baring Head, Chatam Island, and Invercargill in New Zealand, and Marion Island in the sub-antarctic Indian Ocean. Normalized biases improved from +63.7% to +3.3% and correlation increased from 0.52 to 0.84. The representation of sea/land interfaces, mesoscale circulations, and precipitation with the higher resolution model played a major role in the simulation of annual concentration trends. Our results recommend caution when comparing or constraining global models using surface concentration observations from coastal stations.

A high resolution global scale groundwater model

NASA Astrophysics Data System (ADS)

de Graaf, Inge; Sutanudjaja, Edwin; van Beek, Rens; Bierkens, Marc

2014-05-01

As the world's largest accessible source of freshwater, groundwater plays a vital role in satisfying the basic needs of human society. It serves as a primary source of drinking water and supplies water for agricultural and industrial activities. During times of drought, groundwater storage provides a large natural buffer against water shortage and sustains flows to rivers and wetlands, supporting ecosystem habitats and biodiversity. Yet, the current generation of global scale hydrological models (GHMs) do not include a groundwater flow component, although it is a crucial part of the hydrological cycle. Thus, a realistic physical representation of the groundwater system that allows for the simulation of groundwater head dynamics and lateral flows is essential for GHMs that increasingly run at finer resolution. In this study we present a global groundwater model with a resolution of 5 arc-minutes (approximately 10 km at the equator) using MODFLOW (McDonald and Harbaugh, 1988). With this global groundwater model we eventually intend to simulate the changes in the groundwater system over time that result from variations in recharge and abstraction. Aquifer schematization and properties of this groundwater model were developed from available global lithological maps and datasets (Dürr et al., 2005; Gleeson et al., 2010; Hartmann and Moosdorf, 2013), combined with our estimate of aquifer thickness for sedimentary basins. We forced the groundwater model with the output from the global hydrological model PCR-GLOBWB (van Beek et al., 2011), specifically the net groundwater recharge and average surface water levels derived from routed channel discharge. For the parameterization, we relied entirely on available global datasets and did not calibrate the model so that it can equally be expanded to data poor environments. Based on our sensitivity analysis, in which we run the model with various hydrogeological parameter settings, we observed that most variance in groundwater depth is explained by variation in saturated conductivity, and, for the sediment basins, also by variation in recharge. We validated simulated groundwater heads with piezometer heads (available from www.glowasis.eu), resulting in a coefficient of determination for sedimentary basins of 0.92 with regression constant of 0.8. This shows the used method is suitable to build a global groundwater model using best available global information, and estimated water table depths are within acceptable accuracy in many parts of the world.

Kawase & McDermott revisited with a proper ocean model.

NASA Astrophysics Data System (ADS)

Jochum, Markus; Poulsen, Mads; Nuterman, Roman

2017-04-01

A suite of experiments with global ocean models is used to test the hypothesis that Southern Ocean (SO) winds can modify the strength of the Atlantic Meridional Overturning Circulation (AMOC). It is found that for 3 and 1 degree resolution models the results are consistent with Toggweiler & Samuels (1995): stronger SO winds lead to a slight increase of the AMOC. In the simulations with 1/10 degree resolution, however, stronger SO winds weaken the AMOC. We show that these different outcomes are determined by the models' representation of topographic Rossby and Kelvin waves. Consistent with previous literature based on theory and idealized models, first baroclinic waves are slower in the coarse resolution models, but still manage to establish a pattern of global response that is similar to the one in the eddy-permitting model. Because of its different stratification, however, the Atlantic signal is transmitted by higher baroclinic modes. In the coarse resolution model these higher modes are dissipated before they reach 30N, whereas in the eddy-permitting model they reach the subpolar gyre undiminished. This inability of non-eddy-permitting ocean models to represent planetary waves with higher baroclinic modes casts doubt on the ability of climate models to represent non-local effects of climate change. Ideas on how to overcome these difficulties will be discussed.

Linking models and data on vegetation structure

NASA Astrophysics Data System (ADS)

Hurtt, G. C.; Fisk, J.; Thomas, R. Q.; Dubayah, R.; Moorcroft, P. R.; Shugart, H. H.

2010-06-01

For more than a century, scientists have recognized the importance of vegetation structure in understanding forest dynamics. Now future satellite missions such as Deformation, Ecosystem Structure, and Dynamics of Ice (DESDynI) hold the potential to provide unprecedented global data on vegetation structure needed to reduce uncertainties in terrestrial carbon dynamics. Here, we briefly review the uses of data on vegetation structure in ecosystem models, develop and analyze theoretical models to quantify model-data requirements, and describe recent progress using a mechanistic modeling approach utilizing a formal scaling method and data on vegetation structure to improve model predictions. Generally, both limited sampling and coarse resolution averaging lead to model initialization error, which in turn is propagated in subsequent model prediction uncertainty and error. In cases with representative sampling, sufficient resolution, and linear dynamics, errors in initialization tend to compensate at larger spatial scales. However, with inadequate sampling, overly coarse resolution data or models, and nonlinear dynamics, errors in initialization lead to prediction error. A robust model-data framework will require both models and data on vegetation structure sufficient to resolve important environmental gradients and tree-level heterogeneity in forest structure globally.

Spatial Representativeness Error in the Ground-Level Observation Networks for Black Carbon Radiation Absorption

NASA Astrophysics Data System (ADS)

Wang, Rong; Andrews, Elisabeth; Balkanski, Yves; Boucher, Olivier; Myhre, Gunnar; Samset, Bjørn Hallvard; Schulz, Michael; Schuster, Gregory L.; Valari, Myrto; Tao, Shu

2018-02-01

There is high uncertainty in the direct radiative forcing of black carbon (BC), an aerosol that strongly absorbs solar radiation. The observation-constrained estimate, which is several times larger than the bottom-up estimate, is influenced by the spatial representativeness error due to the mesoscale inhomogeneity of the aerosol fields and the relatively low resolution of global chemistry-transport models. Here we evaluated the spatial representativeness error for two widely used observational networks (AErosol RObotic NETwork and Global Atmosphere Watch) by downscaling the geospatial grid in a global model of BC aerosol absorption optical depth to 0.1° × 0.1°. Comparing the models at a spatial resolution of 2° × 2° with BC aerosol absorption at AErosol RObotic NETwork sites (which are commonly located near emission hot spots) tends to cause a global spatial representativeness error of 30%, as a positive bias for the current top-down estimate of global BC direct radiative forcing. By contrast, the global spatial representativeness error will be 7% for the Global Atmosphere Watch network, because the sites are located in such a way that there are almost an equal number of sites with positive or negative representativeness error.

For how long can we predict the weather? - Insights into atmospheric predictability from global convection-allowing simulations

NASA Astrophysics Data System (ADS)

Judt, Falko

2017-04-01

A tremendous increase in computing power has facilitated the advent of global convection-resolving numerical weather prediction (NWP) models. Although this technological breakthrough allows for the seamless prediction of weather from local to global scales, the predictability of multiscale weather phenomena in these models is not very well known. To address this issue, we conducted a global high-resolution (4-km) predictability experiment using the Model for Prediction Across Scales (MPAS), a state-of-the-art global NWP model developed at the National Center for Atmospheric Research. The goals of this experiment are to investigate error growth from convective to planetary scales and to quantify the intrinsic, scale-dependent predictability limits of atmospheric motions. The globally uniform resolution of 4 km allows for the explicit treatment of organized deep moist convection, alleviating grave limitations of previous predictability studies that either used high-resolution limited-area models or global simulations with coarser grids and cumulus parameterization. Error growth is analyzed within the context of an "identical twin" experiment setup: the error is defined as the difference between a 20-day long "nature run" and a simulation that was perturbed with small-amplitude noise, but is otherwise identical. It is found that in convectively active regions, errors grow by several orders of magnitude within the first 24 h ("super-exponential growth"). The errors then spread to larger scales and begin a phase of exponential growth after 2-3 days when contaminating the baroclinic zones. After 16 days, the globally averaged error saturates—suggesting that the intrinsic limit of atmospheric predictability (in a general sense) is about two weeks, which is in line with earlier estimates. However, error growth rates differ between the tropics and mid-latitudes as well as between the troposphere and stratosphere, highlighting that atmospheric predictability is a complex problem. The comparatively slower error growth in the tropics and in the stratosphere indicates that certain weather phenomena could potentially have longer predictability than currently thought.

Towards a Fine-Resolution Global Coupled Climate System for Prediction on Decadal/Centennial Scales

DOE Office of Scientific and Technical Information (OSTI.GOV)

McClean, Julie L.

The over-arching goal of this project was to contribute to the realization of a fully coupled fine resolution Earth System Model simulation in which a weather-scale atmosphere is coupled to an ocean in which mesoscale eddies are largely resolved. Both a prototype fine-resolution fully coupled ESM simulation and a first-ever multi-decadal forced fine-resolution global coupled ocean/ice simulation were configured, tested, run, and analyzed as part of this grant. Science questions focused on the gains from the use of high horizontal resolution, particularly in the ocean and sea-ice, with respect to climatically important processes. Both these fine resolution coupled ocean/sea icemore » and fully-coupled simulations and precedent stand-alone eddy-resolving ocean and eddy-permitting coupled ocean/ice simulations were used to explore the high resolution regime. Overall, these studies showed that the presence of mesoscale eddies significantly impacted mixing processes and the global meridional overturning circulation in the ocean simulations. Fourteen refereed publications and a Ph.D. dissertation resulted from this grant.« less

Multi-RCM ensemble downscaling of global seasonal forecasts (MRED)

NASA Astrophysics Data System (ADS)

Arritt, R. W.

2008-12-01

The Multi-RCM Ensemble Downscaling (MRED) project was recently initiated to address the question, Can regional climate models provide additional useful information from global seasonal forecasts? MRED will use a suite of regional climate models to downscale seasonal forecasts produced by the new National Centers for Environmental Prediction (NCEP) Climate Forecast System (CFS) seasonal forecast system and the NASA GEOS5 system. The initial focus will be on wintertime forecasts in order to evaluate topographic forcing, snowmelt, and the potential usefulness of higher resolution, especially for near-surface fields influenced by high resolution orography. Each regional model will cover the conterminous US (CONUS) at approximately 32 km resolution, and will perform an ensemble of 15 runs for each year 1982-2003 for the forecast period 1 December - 30 April. MRED will compare individual regional and global forecasts as well as ensemble mean precipitation and temperature forecasts, which are currently being used to drive macroscale land surface models (LSMs), as well as wind, humidity, radiation, turbulent heat fluxes, which are important for more advanced coupled macro-scale hydrologic models. Metrics of ensemble spread will also be evaluated. Extensive analysis will be performed to link improvements in downscaled forecast skill to regional forcings and physical mechanisms. Our overarching goal is to determine what additional skill can be provided by a community ensemble of high resolution regional models, which we believe will eventually define a strategy for more skillful and useful regional seasonal climate forecasts.

The Navy’s Application of Ocean Forecasting to Decision Support

DTIC Science & Technology

2014-09-01

Prediction Center (OPC) website for graphics or the National Operational Model Archive and Distribution System ( NOMADS ) for data files. Regional...inputs: » GLOBE = Global Land One-km Base Elevation » WVS = World Vector Shoreline » DBDB2 = Digital Bathymetry Data Base 2 minute resolution » DBDBV... Digital Bathymetry Data Base variable resolution Oceanography | Vol. 27, No.3130 Very High-Resolution Coastal Circulation Models Nearshore

Spatiotemporal models of global soil organic carbon stock to support land degradation assessments at regional and global scales: limitations, challenges and opportunities

NASA Astrophysics Data System (ADS)

Hengl, Tomislav; Heuvelink, Gerard; Sanderman, Jonathan; MacMillan, Robert

2017-04-01

There is an increasing interest in fitting and applying spatiotemporal models that can be used to assess and monitor soil organic carbon stocks (SOCS), for example, in support of the '4 pourmille' initiative aiming at soil carbon sequestration towards climate change adaptation and mitigation and UN's Land Degradation Neutrality indicators and similar degradation assessment projects at regional and global scales. The land cover mapping community has already produced several spatiotemporal data sets with global coverage and at relatively fine resolution e.g. USGS MODIS land cover annual maps for period 2000-2014; European Space Agency land cover maps at 300 m resolution for the year 2000, 2005 and 2010; Chinese GlobeLand30 dataset available for years 2000 and 2010; Columbia University's WRI GlobalForestWatch with deforestation maps at 30 m resolution for the period 2000-2016 (Hansen et al. 2013). These data sets can be used for land degradation assessment and scenario testing at global and regional scales (Wei et al 2014). Currently, however, no compatible global spatiotemporal data sets exist on status of soil quality and/or soil health (Powlson et al. 2013). This paper describes an initial effort to devise and evaluate a procedure for mapping spatio-temporal changes in SOC stocks using a complete stack of soil forming factors (climate, relief, land cover, land use, lithology and living organisms) represented mainly through remote sensing based time series of Earth images. For model building we used some 75,000 geo-referenced soil profiles and a stacks space-time covariates (land cover, land use, biomass, climate) at two standard resolutions: (1) 10 km resolution with data available for period 1920-2014 and (2) 1000 m resolution with data available for period 2000-2014. The initial results show that, although it is technically feasible to produce space time estimates of SOCS that demonstrate the procedure, the estimates are relatively uncertain (<45% of variation explained) and lead to obvious artifacts, especially in areas that have not be represented in time-dimension (temporal extrapolation). For some regions that possess somewhat adequate amounts of point data in space and time (e.g. USA) relatively credible space time estimates can be produced. By adding more training data (both legacy and newly collected points) these models can be gradually improved until they can become operational for decision making and scenario testing.

Atlantic Tropical Cyclogenetic Processes during SOP-3 NAMMA in the GEOS-5 Global Data Assimilation and Forecast System

NASA Technical Reports Server (NTRS)

Reale, Oreste; Lau, William K.; Kim, Kyu-Myong; Brin, Eugenia

2009-01-01

This article investigates the role of the Saharan Air Layer (SAL) in tropical cyclogenetic processes associated with a non-developing and a developing African easterly wave observed during the Special Observation Period (SOP-3) phase of the 2006 NASA African Monsoon Multidisciplinary Analyses (NAMMA). The two waves are chosen because both interact heavily with Saharan air. A global data assimilation and forecast system, the NASA GEOS-5, is being run to produce a set of high-quality global analyses, inclusive of all observations used operationally but with denser satellite information. In particular, following previous works by the same Authors, the quality-controlled data from the Atmospheric Infrared Sounder (AIRS) used to produce these analyses have a better coverage than the one adopted by operational centers. From these improved analyses, two sets of 31 5-day high resolution forecasts, at horizontal resolutions of both half and quarter degrees, are produced. Results show that very steep moisture gradients are associated with the SAL in forecasts and analyses even at great distance from the Sahara. In addition, a thermal dipole (warm above, cool below) is present in the non-developing case. Moderate Resolution Imaging Spectroradiometer (MODIS) show that aerosol optical thickness is higher in the non-developing case. Altogether, results suggest that radiative effect of dust may play some role in producing a thermal structure less favorable to cyclogenesis. Results also indicate that only global horizontal resolutions on the order of 20-30 kilometers can capture the large-scale transport and the fine thermal structure of the SAL, inclusive of the sharp moisture gradients, reproducing the effect of tropical cyclone suppression which has been hypothesized by previous authors from observational and regional modeling perspectives. These effects cannot be fully represented at lower resolutions. Global resolution of a quarter of a degree is a minimum critical threshold to investigate Atlantic tropical cyclogenesis from a global modeling perspective.

Mesoscale Numerical Simulations of the IAS Circulation

NASA Astrophysics Data System (ADS)

Mooers, C. N.; Ko, D.

2008-05-01

Real-time nowcasts and forecasts of the IAS circulation have been made for several years with mesoscale resolution using the Navy Coastal Ocean Model (NCOM) implemented for the IAS. It is commonly called IASNFS and is driven by the lower resolution Global NCOM on the open boundaries, synoptic atmospheric forcing obtained from the Navy Global Atmospheric Prediction System (NOGAPS), and assimilated satellite-derived sea surface height anomalies and sea surface temperature. Here, examples of the model output are demonstrated; e.g., Gulf of Mexico Loop Current eddy shedding events and the meandering Caribbean Current jet and associated eddies. Overall, IASNFS is ready for further analysis, application to a variety of studies, and downscaling to even higher resolution shelf models. Its output fields are available online through NOAA's National Coastal Data Development Center (NCDDC), located at the Stennis Space Center.

Global Weather Prediction and High-End Computing at NASA

NASA Technical Reports Server (NTRS)

Lin, Shian-Jiann; Atlas, Robert; Yeh, Kao-San

2003-01-01

We demonstrate current capabilities of the NASA finite-volume General Circulation Model an high-resolution global weather prediction, and discuss its development path in the foreseeable future. This model can be regarded as a prototype of a future NASA Earth modeling system intended to unify development activities cutting across various disciplines within the NASA Earth Science Enterprise.

Global distribution of urban parameters derived from high-resolution global datasets for weather modelling

NASA Astrophysics Data System (ADS)

Kawano, N.; Varquez, A. C. G.; Dong, Y.; Kanda, M.

2016-12-01

Numerical model such as Weather Research and Forecasting model coupled with single-layer Urban Canopy Model (WRF-UCM) is one of the powerful tools to investigate urban heat island. Urban parameters such as average building height (Have), plain area index (λp) and frontal area index (λf), are necessary inputs for the model. In general, these parameters are uniformly assumed in WRF-UCM but this leads to unrealistic urban representation. Distributed urban parameters can also be incorporated into WRF-UCM to consider a detail urban effect. The problem is that distributed building information is not readily available for most megacities especially in developing countries. Furthermore, acquiring real building parameters often require huge amount of time and money. In this study, we investigated the potential of using globally available satellite-captured datasets for the estimation of the parameters, Have, λp, and λf. Global datasets comprised of high spatial resolution population dataset (LandScan by Oak Ridge National Laboratory), nighttime lights (NOAA), and vegetation fraction (NASA). True samples of Have, λp, and λf were acquired from actual building footprints from satellite images and 3D building database of Tokyo, New York, Paris, Melbourne, Istanbul, Jakarta and so on. Regression equations were then derived from the block-averaging of spatial pairs of real parameters and global datasets. Results show that two regression curves to estimate Have and λf from the combination of population and nightlight are necessary depending on the city's level of development. An index which can be used to decide which equation to use for a city is the Gross Domestic Product (GDP). On the other hand, λphas less dependence on GDP but indicated a negative relationship to vegetation fraction. Finally, a simplified but precise approximation of urban parameters through readily-available, high-resolution global datasets and our derived regressions can be utilized to estimate a global distribution of urban parameters for later incorporation into a weather model, thus allowing us to acquire a global understanding of urban climate (Global Urban Climatology). Acknowledgment: This research was supported by the Environment Research and Technology Development Fund (S-14) of the Ministry of the Environment, Japan.

Climate Model Tests of the Early Anthropogenic Hypothesis

NASA Astrophysics Data System (ADS)

Vavrus, S.; Kutzbach, J.; Philippon, G.

2008-12-01

We test the hypothesis that greenhouse gas emissions produced by the combination of early and recent human activities, augmented by additional rises in greenhouse gases through ocean feedbacks, have kept the climate warmer than its natural level and offset an incipient glaciation. We use four different configurations of NCAR's Community Climate System Model to investigate the natural climate that should exist today if CO2 and CH4 concentrations had fallen to their average levels reached during previous interglaciations. The model simulations consist of three using a coupled atmosphere-slab ocean configuration---fixed land cover at moderate (T42) and high (T85) model resolution and interactive vegetation composition at T42 resolution--and one employing a coupled atmosphere-dynamical ocean configuration and fixed land cover at T42 resolution. With greenhouse gas concentrations lowered to their estimated natural levels, global mean temperature falls by 2.5-3.0 K in all four experiments. Of the total global cooling with fixed land cover and moderate model resolution, 38% (62%) is attributable to early agricultural activities (industrialization), while early agriculture accounts for approximately half of the expanded permanent snow cover area. The greenhouse cooling triggers widespread glacial inception in the Northern Hemisphere, where permanent snow cover expands by at least 80% and even more with the addition of enhanced model processes: 130% with the dynamical ocean, 150% with high (T85) model resolution, and 200% with vegetation feedbacks included. The regional pattern of incipient glaciation is strongly influenced by atmospheric and circulation changes, sea ice feedbacks, and model resolution. The simulation with a dynamical ocean produces a decrease in vertically integrated global ocean temperature of 1.25 K, a 20% weakening of the Atlantic meridional overturning cell, and an expansion of sea ice and reduced upwelling in the Southern Ocean. Viewed from the perspective of explaining the unusual late-Holocene increases of CO2 that occurred prior to the Industrial Revolution, these simulated changes in ocean temperature, sea ice cover, and circulation (with sign reversed) support the hypothesis that early agriculture played a role in initiating anomalous warming that thwarted incipient glaciation beginning several thousand years ago. Decreased ocean solubility globally and positive ocean/sea-ice feedbacks in the Southern Hemisphere probably augmented the initial CO2 increase and caused additional warming.

The Fire INventory from NCAR (FINN): a high resolution global model to estimate the emissions from open burning

Treesearch

C. Wiedinmyer; S. K. Akagi; R. J. Yokelson; L. K. Emmons; J. A. Al-Saadi; J. J. Orlando; A. J. Soja

2010-01-01

The Fire INventory from NCAR version 1.0 (FINNv1) provides daily, 1 km resolution, global estimates of the trace gas and particle emissions from open burning of biomass, which includes wildfire, agricultural fires, and prescribed burning and does not include 5 biofuel use and trash burning. Emission factors used in the calculations have been updated with recent data,...

High-resolution RCMs as pioneers for future GCMs

NASA Astrophysics Data System (ADS)

Schar, C.; Ban, N.; Arteaga, A.; Charpilloz, C.; Di Girolamo, S.; Fuhrer, O.; Hoefler, T.; Leutwyler, D.; Lüthi, D.; Piaget, N.; Ruedisuehli, S.; Schlemmer, L.; Schulthess, T. C.; Wernli, H.

2017-12-01

Currently large efforts are underway to refine the horizontal resolution of global and regional climate models to O(1 km), with the intent to represent convective clouds explicitly rather than using semi-empirical parameterizations. This refinement will move the governing equations closer to first principles and is expected to reduce the uncertainties of climate models. High resolution is particularly attractive in order to better represent critical cloud feedback processes (e.g. related to global climate sensitivity and extratropical summer convection) and extreme events (such as heavy precipitation events, floods, and hurricanes). The presentation will be illustrated using decade-long simulations at 2 km horizontal grid spacing, some of these covering the European continent on a computational mesh with 1536x1536x60 grid points. To accomplish such simulations, use is made of emerging heterogeneous supercomputing architectures, using a version of the COSMO limited-area weather and climate model that is able to run entirely on GPUs. Results show that kilometer-scale resolution dramatically improves the simulation of precipitation in terms of the diurnal cycle and short-term extremes. The modeling framework is used to address changes of precipitation scaling with climate change. It is argued that already today, modern supercomputers would in principle enable global atmospheric convection-resolving climate simulations, provided appropriately refactored codes were available, and provided solutions were found to cope with the rapidly growing output volume. A discussion will be provided of key challenges affecting the design of future high-resolution climate models. It is suggested that km-scale RCMs should be exploited to pioneer this terrain, at a time when GCMs are not yet available at such resolutions. Areas of interest include the development of new parameterization schemes adequate for km-scale resolution, the exploration of new validation methodologies and data sets, the assessment of regional-scale climate feedback processes, and the development of alternative output analysis methodologies.

High resolution global climate modelling; the UPSCALE project, a large simulation campaign

NASA Astrophysics Data System (ADS)

Mizielinski, M. S.; Roberts, M. J.; Vidale, P. L.; Schiemann, R.; Demory, M.-E.; Strachan, J.; Edwards, T.; Stephens, A.; Lawrence, B. N.; Pritchard, M.; Chiu, P.; Iwi, A.; Churchill, J.; del Cano Novales, C.; Kettleborough, J.; Roseblade, W.; Selwood, P.; Foster, M.; Glover, M.; Malcolm, A.

2014-01-01

The UPSCALE (UK on PRACE: weather-resolving Simulations of Climate for globAL Environmental risk) project constructed and ran an ensemble of HadGEM3 (Hadley centre Global Environment Model 3) atmosphere-only global climate simulations over the period 1985-2011, at resolutions of N512 (25 km), N216 (60 km) and N96 (130 km) as used in current global weather forecasting, seasonal prediction and climate modelling respectively. Alongside these present climate simulations a parallel ensemble looking at extremes of future climate was run, using a time-slice methodology to consider conditions at the end of this century. These simulations were primarily performed using a 144 million core hour, single year grant of computing time from PRACE (the Partnership for Advanced Computing in Europe) in 2012, with additional resources supplied by the Natural Environmental Research Council (NERC) and the Met Office. Almost 400 terabytes of simulation data were generated on the HERMIT supercomputer at the high performance computing center Stuttgart (HLRS), and transferred to the JASMIN super-data cluster provided by the Science and Technology Facilities Council Centre for Data Archival (STFC CEDA) for analysis and storage. In this paper we describe the implementation of the project, present the technical challenges in terms of optimisation, data output, transfer and storage that such a project involves and include details of the model configuration and the composition of the UPSCALE dataset. This dataset is available for scientific analysis to allow assessment of the value of model resolution in both present and potential future climate conditions.

High-resolution global climate modelling: the UPSCALE project, a large-simulation campaign

NASA Astrophysics Data System (ADS)

Mizielinski, M. S.; Roberts, M. J.; Vidale, P. L.; Schiemann, R.; Demory, M.-E.; Strachan, J.; Edwards, T.; Stephens, A.; Lawrence, B. N.; Pritchard, M.; Chiu, P.; Iwi, A.; Churchill, J.; del Cano Novales, C.; Kettleborough, J.; Roseblade, W.; Selwood, P.; Foster, M.; Glover, M.; Malcolm, A.

2014-08-01

The UPSCALE (UK on PRACE: weather-resolving Simulations of Climate for globAL Environmental risk) project constructed and ran an ensemble of HadGEM3 (Hadley Centre Global Environment Model 3) atmosphere-only global climate simulations over the period 1985-2011, at resolutions of N512 (25 km), N216 (60 km) and N96 (130 km) as used in current global weather forecasting, seasonal prediction and climate modelling respectively. Alongside these present climate simulations a parallel ensemble looking at extremes of future climate was run, using a time-slice methodology to consider conditions at the end of this century. These simulations were primarily performed using a 144 million core hour, single year grant of computing time from PRACE (the Partnership for Advanced Computing in Europe) in 2012, with additional resources supplied by the Natural Environment Research Council (NERC) and the Met Office. Almost 400 terabytes of simulation data were generated on the HERMIT supercomputer at the High Performance Computing Center Stuttgart (HLRS), and transferred to the JASMIN super-data cluster provided by the Science and Technology Facilities Council Centre for Data Archival (STFC CEDA) for analysis and storage. In this paper we describe the implementation of the project, present the technical challenges in terms of optimisation, data output, transfer and storage that such a project involves and include details of the model configuration and the composition of the UPSCALE data set. This data set is available for scientific analysis to allow assessment of the value of model resolution in both present and potential future climate conditions.

The contribution of future agricultural trends in the US Midwest to global climate change mitigation

DOE Office of Scientific and Technical Information (OSTI.GOV)

Thomson, Allison M.; Kyle, G. Page; Zhang, Xuesong

2014-01-19

Land use change is a complex response to changing environmental and socioeconomic systems. Historical drivers of land use change include changes in the natural resource availability of a region, changes in economic conditions for production of certain products and changing policies. Most recently, introduction of policy incentives for biofuel production have influenced land use change in the US Midwest, leading to concerns that bioenergy production systems may compete with food production and land conservation. Here we explore how land use may be impacted by future climate mitigation measures by nesting a high resolution agricultural model (EPIC – Environmental Policy Indicatormore » Climate) for the US Midwest within a global integrated assessment model (GCAM – Global Change Assessment Model). This approach is designed to provide greater spatial resolution and detailed agricultural practice information by focusing on the climate mitigation potential of agriculture and land use in a specific region, while retaining the global economic context necessary to understand the far ranging effects of climate mitigation targets. We find that until the simulated carbon prices are very high, the US Midwest has a comparative advantage in producing traditional food and feed crops over bioenergy crops. Overall, the model responds to multiple pressures by adopting a mix of future responses. We also find that the GCAM model is capable of simulations at multiple spatial scales and agricultural technology resolution, which provides the capability to examine regional response to global policy and economic conditions in the context of climate mitigation.« less

GEOS S2S-2_1: The GMAO new high resolution Seasonal Prediction System

NASA Astrophysics Data System (ADS)

Molod, A.; Vikhliaev, Y. V.; Hackert, E. C.; Kovach, R. M.; Zhao, B.; Cullather, R. I.; Marshak, J.; Borovikov, A.; Li, Z.; Barahona, D.; Andrews, L. C.; Chang, Y.; Schubert, S. D.; Koster, R. D.; Suarez, M.; Akella, S.

2017-12-01

A new version of the modeling and analysis system used to produce subseasonalto seasonal forecasts has just been released by the NASA/Goddard GlobalModeling and Assimilation Office. The new version runs at higher atmospheric resolution (approximately 1/2 degree globally), contains a subtantially improvedmodel description of the cryosphere, and includes additional interactive earth system model components (aerosol model). In addition, the Ocean data assimilationsystem has been replaced with a Local Ensemble Transform Kalman Filter.Here will describe the new system, along with the plans for the future (GEOS S2S-3_0) which will include a higher resolution ocean model and more interactive earth system model components (interactive vegetation, biomass burning from fires). We will alsopresent results from a free-running coupled simulation with the new system and resultsfrom a series of retrospective seasonal forecasts.Results from retrospective forecasts show significant improvements in surface temperaturesover much of the northern hemisphere and a much improved prediction of sea ice extent in bothhemispheres. The precipitation forecast skill is comparable to previous S2S systems, andthe only tradeoff is an increased "double ITCZ", which is expected as we go to higher atmospheric resolution.

The 4th order GISS model of the global atmosphere

NASA Technical Reports Server (NTRS)

Kalnay-Rivas, E.; Bayliss, A.; Storch, J.

1977-01-01

The new GISS 4th order model of the global atmosphere is described. It is based on 4th order quadratically conservative differences with the periodic application of a 16th order filter on the sea level pressure and potential temperature equations, a combination which is approximately enstrophy conserving. Several short range forecasts indicate a significant improvement over 2nd order forecasts with the same resolution (approximately 400 km). However the 4th order forecasts are somewhat inferior to 2nd order forecasts with double resolution. This is probably due to the presence of short waves in the range between 1000 km and 2000 km, which are computed more accurately by the 2nd order high resolution model. An operation count of the schemes indicates that with similar code optimization, the 4th order model will require approximately the same amount of computer time as the 2nd order model with the same resolution. It is estimated that the 4th order model with a grid size of 200 km provides enough accuracy to make horizontal truncation errors negligible over a period of a week for all synoptic scales (waves longer than 1000 km).

Concurrent Simulation of the Eddying General Circulation and Tides in a Global Ocean Model

DTIC Science & Technology

2010-01-01

Eddying General Circulation and Tides in a Global Ocean Model 5a. CONTRACT NUMBER 5b. GRANT NUMBER 5c. PROGRAM ELEMENT NUMBER 0602435N 6...STATEMENT Approved for public release, distribution is unlimited. 13. SUPPLEMENTARY NOTES 14. ABSTRACT This paper presents a five-year global ...running 25-h average to approximately separate tidal and non-tidal components of the near-bottom flow. In contrast to earlier high-resolution global

A Modeling Approach to Global Land Surface Monitoring with Low Resolution Satellite Imaging

NASA Technical Reports Server (NTRS)

Hlavka, Christine A.; Dungan, Jennifer; Livingston, Gerry P.; Gore, Warren J. (Technical Monitor)

1998-01-01

The effects of changing land use/land cover on global climate and ecosystems due to greenhouse gas emissions and changing energy and nutrient exchange rates are being addressed by federal programs such as NASA's Mission to Planet Earth (MTPE) and by international efforts such as the International Geosphere-Biosphere Program (IGBP). The quantification of these effects depends on accurate estimates of the global extent of critical land cover types such as fire scars in tropical savannas and ponds in Arctic tundra. To address the requirement for accurate areal estimates, methods for producing regional to global maps with satellite imagery are being developed. The only practical way to produce maps over large regions of the globe is with data of coarse spatial resolution, such as Advanced Very High Resolution Radiometer (AVHRR) weather satellite imagery at 1.1 km resolution or European Remote-Sensing Satellite (ERS) radar imagery at 100 m resolution. The accuracy of pixel counts as areal estimates is in doubt, especially for highly fragmented cover types such as fire scars and ponds. Efforts to improve areal estimates from coarse resolution maps have involved regression of apparent area from coarse data versus that from fine resolution in sample areas, but it has proven difficult to acquire sufficient fine scale data to develop the regression. A method for computing accurate estimates from coarse resolution maps using little or no fine data is therefore needed.

Geoscience Meets Social Science: A Flexible Data Driven Approach for Developing High Resolution Population Datasets at Global Scale

NASA Astrophysics Data System (ADS)

Rose, A.; McKee, J.; Weber, E.; Bhaduri, B. L.

2017-12-01

Leveraging decades of expertise in population modeling, and in response to growing demand for higher resolution population data, Oak Ridge National Laboratory is now generating LandScan HD at global scale. LandScan HD is conceived as a 90m resolution population distribution where modeling is tailored to the unique geography and data conditions of individual countries or regions by combining social, cultural, physiographic, and other information with novel geocomputation methods. Similarities among these areas are exploited in order to leverage existing training data and machine learning algorithms to rapidly scale development. Drawing on ORNL's unique set of capabilities, LandScan HD adapts highly mature population modeling methods developed for LandScan Global and LandScan USA, settlement mapping research and production in high-performance computing (HPC) environments, land use and neighborhood mapping through image segmentation, and facility-specific population density models. Adopting a flexible methodology to accommodate different geographic areas, LandScan HD accounts for the availability, completeness, and level of detail of relevant ancillary data. Beyond core population and mapped settlement inputs, these factors determine the model complexity for an area, requiring that for any given area, a data-driven model could support either a simple top-down approach, a more detailed bottom-up approach, or a hybrid approach.

Interactions between Arctic sea ice drift, concentration and thickness modeled by NEMO-LIM3 at different resolutions

NASA Astrophysics Data System (ADS)

Docquier, David; Massonnet, François; Raulier, Jonathan; Lecomte, Olivier; Fichefet, Thierry

2016-04-01

Sea ice concentration and thickness have substantially decreased in the Arctic since the beginning of the satellite era. As a result, mechanical strength has decreased allowing more fracturing and leading to increased sea ice drift. However, recent studies have highlighted that the interplay between sea ice thermodynamics and dynamics is poorly represented in contemporary global climate model (GCM) simulations. Thus, the considerable inter-model spread in terms of future sea ice extent projections could be reduced by better understanding the interactions between drift, concentration and thickness. This study focuses on the results coming from the global coupled ocean-sea ice model NEMO-LIM3 between 1979 and 2012. Three different simulations are forced by the Drakkar Forcing Set (DFS) 5.2 and run on the global tripolar ORCA grid at spatial resolutions of 0.25, 1° and 2°. The relation between modeled sea ice drift, concentration and thickness is further analyzed, compared to observations and discussed in the framework of the above-mentioned poor representation. It is proposed as a process-based metric for evaluating model performance. This study forms part of the EU Horizon 2020 PRIMAVERA project aiming at developing a new generation of advanced and well-evaluated high-resolution GCMs.

Global seasonal climate predictability in a two tiered forecast system: part I: boreal summer and fall seasons

NASA Astrophysics Data System (ADS)

Misra, Vasubandhu; Li, H.; Wu, Z.; DiNapoli, S.

2014-03-01

This paper shows demonstrable improvement in the global seasonal climate predictability of boreal summer (at zero lead) and fall (at one season lead) seasonal mean precipitation and surface temperature from a two-tiered seasonal hindcast forced with forecasted SST relative to two other contemporary operational coupled ocean-atmosphere climate models. The results from an extensive set of seasonal hindcasts are analyzed to come to this conclusion. This improvement is attributed to: (1) The multi-model bias corrected SST used to force the atmospheric model. (2) The global atmospheric model which is run at a relatively high resolution of 50 km grid resolution compared to the two other coupled ocean-atmosphere models. (3) The physics of the atmospheric model, especially that related to the convective parameterization scheme. The results of the seasonal hindcast are analyzed for both deterministic and probabilistic skill. The probabilistic skill analysis shows that significant forecast skill can be harvested from these seasonal hindcasts relative to the deterministic skill analysis. The paper concludes that the coupled ocean-atmosphere seasonal hindcasts have reached a reasonable fidelity to exploit their SST anomaly forecasts to force such relatively higher resolution two tier prediction experiments to glean further boreal summer and fall seasonal prediction skill.

Crustal density contrast detection by global gravity and topography models and in-situ gravity observations

NASA Astrophysics Data System (ADS)

Claessens, S. J.

2016-12-01

Mass density contrasts in the Earth's crust can be detected using an inversion of terrestrial or airborne gravity data. This contribution shows a technique to detect short-scale density contrasts using in-situ gravity observations in combination with a high-resolution global gravity model that includes variations in the gravity field due to topography. The technique is exemplified at various test sites using the Global Gravity Model Plus (GGMplus), which is a 7.2 arcsec resolution model of the Earth's gravitational field, covering all land masses and near-coastal areas within +/- 60° latitude. The model is a composite of GRACE and GOCE satellite observations, the EGM2008 global gravity model, and short-scale topographic gravity effects. Since variations in the Earth's gravity field due to topography are successfully modelled by GGMplus, any remaining differences with in-situ gravity observations are primarily due to mass density variations. It is shown that this technique effectively filters out large-scale density variations, and highlights short-scale near-surface density contrasts in the Earth's crust. Numerical results using recent high-density gravity surveys are presented, which indicate a strong correlation between density contrasts found and known lines of geological significance.

The effect of horizontal resolution on simulation quality in the Community Atmospheric Model, CAM5.1

DOE PAGES

Wehner, Michael F.; Reed, Kevin A.; Li, Fuyu; ...

2014-10-13

We present an analysis of version 5.1 of the Community Atmospheric Model (CAM5.1) at a high horizontal resolution. Intercomparison of this global model at approximately 0.25°, 1°, and 2° is presented for extreme daily precipitation as well as for a suite of seasonal mean fields. In general, extreme precipitation amounts are larger in high resolution than in lower-resolution configurations. In many but not all locations and/or seasons, extreme daily precipitation rates in the high-resolution configuration are higher and more realistic. The high-resolution configuration produces tropical cyclones up to category 5 on the Saffir-Simpson scale and a comparison to observations revealsmore » both realistic and unrealistic model behavior. In the absence of extensive model tuning at high resolution, simulation of many of the mean fields analyzed in this study is degraded compared to the tuned lower-resolution public released version of the model.« less

Large Scale Flood Risk Analysis using a New Hyper-resolution Population Dataset

NASA Astrophysics Data System (ADS)

Smith, A.; Neal, J. C.; Bates, P. D.; Quinn, N.; Wing, O.

2017-12-01

Here we present the first national scale flood risk analyses, using high resolution Facebook Connectivity Lab population data and data from a hyper resolution flood hazard model. In recent years the field of large scale hydraulic modelling has been transformed by new remotely sensed datasets, improved process representation, highly efficient flow algorithms and increases in computational power. These developments have allowed flood risk analysis to be undertaken in previously unmodeled territories and from continental to global scales. Flood risk analyses are typically conducted via the integration of modelled water depths with an exposure dataset. Over large scales and in data poor areas, these exposure data typically take the form of a gridded population dataset, estimating population density using remotely sensed data and/or locally available census data. The local nature of flooding dictates that for robust flood risk analysis to be undertaken both hazard and exposure data should sufficiently resolve local scale features. Global flood frameworks are enabling flood hazard data to produced at 90m resolution, resulting in a mis-match with available population datasets which are typically more coarsely resolved. Moreover, these exposure data are typically focused on urban areas and struggle to represent rural populations. In this study we integrate a new population dataset with a global flood hazard model. The population dataset was produced by the Connectivity Lab at Facebook, providing gridded population data at 5m resolution, representing a resolution increase over previous countrywide data sets of multiple orders of magnitude. Flood risk analysis undertaken over a number of developing countries are presented, along with a comparison of flood risk analyses undertaken using pre-existing population datasets.

ESiWACE: A Center of Excellence for HPC applications to support cloud resolving earth system modelling

NASA Astrophysics Data System (ADS)

Biercamp, Joachim; Adamidis, Panagiotis; Neumann, Philipp

2017-04-01

With the exa-scale era approaching, length and time scales used for climate research on one hand and numerical weather prediction on the other hand blend into each other. The Centre of Excellence in Simulation of Weather and Climate in Europe (ESiWACE) represents a European consortium comprising partners from climate, weather and HPC in their effort to address key scientific challenges that both communities have in common. A particular challenge is to reach global models with spatial resolutions that allow simulating convective clouds and small-scale ocean eddies. These simulations would produce better predictions of trends and provide much more fidelity in the representation of high-impact regional events. However, running such models in operational mode, i.e with sufficient throughput in ensemble mode clearly will require exa-scale computing and data handling capability. We will discuss the ESiWACE initiative and relate it to work-in-progress on high-resolution simulations in Europe. We present recent strong scalability measurements from ESiWACE to demonstrate current computability in weather and climate simulation. A special focus in this particular talk is on the Icosahedal Nonhydrostatic (ICON) model used for a comparison of high resolution regional and global simulations with high quality observation data. We demonstrate that close-to-optimal parallel efficiency can be achieved in strong scaling global resolution experiments on Mistral/DKRZ, e.g. 94% for 5km resolution simulations using 36k cores on Mistral/DKRZ. Based on our scalability and high-resolution experiments, we deduce and extrapolate future capabilities for ICON that are expected for weather and climate research at exascale.

Tropospheric carbon monoxide over the Pacific during HIPPO: two-way coupled simulation of GEOS-Chem and its multiple nested models

NASA Astrophysics Data System (ADS)

Yan, Y.-Y.; Lin, J.-T.; Kuang, Y.; Yang, D.; Zhang, L.

2014-12-01

Global chemical transport models (CTMs) are used extensively to study air pollution and transport at a global scale. These models are limited by coarse horizontal resolutions that do not allow for a detailed representation of small-scale nonlinear processes over the pollutant source regions. Here we couple the global GEOS-Chem CTM and its three high-resolution nested models to simulate the tropospheric carbon monoxide (CO) over the Pacific Ocean during five High-performance Instrumented Airborne Platform for Environmental Research (HIAPER) Pole-to-Pole Observations (HIPPO) campaigns between 2009 and 2011. We develop a two-way coupler, the PeKing University CouPLer (PKUCPL), allowing for the exchange and interaction of chemical constituents between the global model (at 2.5° long. × 2° lat.) and the three nested models (at 0.667° long. × 0.5° lat.) covering Asia, North America, and Europe. The coupler obtains nested model results to modify the global model simulation within the respective nested domains, and simultaneously acquires global model results to provide lateral boundary conditions (LBCs) for the nested models. Compared to the global model alone, the two-way coupled simulation results in enhanced CO concentrations in the nested domains. Sensitivity tests suggest the enhancement to be a result of improved representation of the spatial distributions of CO, nitrogen oxides, and non-methane volatile organic compounds, the meteorological dependence of natural emissions, and other resolution-dependent processes. The relatively long lifetime of CO allows for the enhancement to be accumulated and carried across the globe. We found that the two-way coupled simulation increased the global tropospheric mean CO concentrations in 2009 by 10.4%, with a greater enhancement at 13.3% in the Northern Hemisphere. Coincidently, the global tropospheric mean hydroxyl radical (OH) was reduced by 4.2%, resulting in a 4.2% enhancement in the methyl chloroform lifetime (MCF; via reaction with tropospheric OH). The resulting CO and OH contents and MCF lifetime are closer to observation-based estimates. Both the global and the two-way coupled models capture the general spatiotemporal patterns of HIPPO CO over the Pacific. The two-way coupled simulation is much closer to HIPPO CO, with a mean bias of 1.1 ppb (1.4%) below 9 km compared to the bias at -7.2 ppb (-9.2%) for the global model alone. The improvement is most apparent over the North Pacific. Our test simulations show that the global model alone could resemble the two-way coupled simulation (especially below 4 km) by increasing its global CO emissions by 15% for HIPPO-1 and HIPPO-3, by 25% for HIPPO-2 and HIPPO-4, and by 35% for HIPPO-5. This has important implications for using the global model alone to constrain CO emissions. Thus, the two-way coupled simulation is a significantly improved model tool for studying the global impacts of air pollutants from major anthropogenic source regions.

The Errors Sources Affect to the Results of One-Way Nested Ocean Regional Circulation Model

NASA Astrophysics Data System (ADS)

Pham, S. V.

2016-02-01

Pham-Van Sy1, Jin Hwan Hwang2 and Hyeyun Ku3 Dept. of Civil & Environmental Engineering, Seoul National University, KoreaEmail: 1phamsymt@gmail.com (Corresponding author) Email: 2jinhwang@snu.ac.krEmail: 3hyeyun.ku@gmail.comAbstractThe Oceanic Regional Circulation Model (ORCM) is an essential tool in resolving highly a regional scale through downscaling dynamically the results from the roughly revolved global model. However, when dynamic downscaling from a coarse resolution of the global model or observations to the small scale, errors are generated due to the different sizes of resolution and lateral updating frequency. This research evaluated the effect of four main sources on the results of the ocean regional circulation model (ORCMs) during downscaling and nesting the output data from the ocean global circulation model (OGCMs). Representative four error sources should be the way of the LBC formulation, the spatial resolution difference between driving and driven data, the frequency for up-dating LBCs and domain size. Errors which are contributed from each error source to the results of the ORCMs are investigated separately by applying the Big-Brother Experiment (BBE). Within resolution of 3km grid point of the ORCMs imposing in the BBE framework, it clearly exposes that the simulation results of the ORCMs significantly depend on the domain size and specially the spatial and temporal resolution of lateral boundary conditions (LBCs). The ratio resolution of spatial resolution between driving data and driven model could be up to 3, the updating frequency of the LBCs can be up to every 6 hours per day. The optimal domain size of the ORCMs could be smaller than the OGCMs' domain size around 2 to 10 times. Key words: ORCMs, error source, lateral boundary conditions, domain size Acknowledgement: This research was supported by grants from the Korean Ministry of Oceans and Fisheries entitled as "Developing total management system for the Keum river estuary and coast" and "Development of Technology for CO2 Marine Geological Storage". We also thank to the administrative supports of the Integrated Research Institute of Construction and Environmental Engineering of the Seoul National University.

DOE Office of Scientific and Technical Information (OSTI.GOV)

Li, Fuyu; Collins, William D.; Wehner, Michael F.

High-resolution climate models have been shown to improve the statistics of tropical storms and hurricanes compared to low-resolution models. The impact of increasing horizontal resolution in the tropical storm simulation is investigated exclusively using a series of Atmospheric Global Climate Model (AGCM) runs with idealized aquaplanet steady-state boundary conditions and a fixed operational storm-tracking algorithm. The results show that increasing horizontal resolution helps to detect more hurricanes, simulate stronger extreme rainfall, and emulate better storm structures in the models. However, increasing model resolution does not necessarily produce stronger hurricanes in terms of maximum wind speed, minimum sea level pressure, andmore » mean precipitation, as the increased number of storms simulated by high-resolution models is mainly associated with weaker storms. The spatial scale at which the analyses are conducted appears to have more important control on these meteorological statistics compared to horizontal resolution of the model grid. When the simulations are analyzed on common low-resolution grids, the statistics of the hurricanes, particularly the hurricane counts, show reduced sensitivity to the horizontal grid resolution and signs of scale invariant.« less

Detection and Attribution of Regional Climate Change

DOE Office of Scientific and Technical Information (OSTI.GOV)

Bala, G; Mirin, A

2007-01-19

We developed a high resolution global coupled modeling capability to perform breakthrough studies of the regional climate change. The atmospheric component in our simulation uses a 1{sup o} latitude x 1.25{sup o} longitude grid which is the finest resolution ever used for the NCAR coupled climate model CCSM3. Substantial testing and slight retuning was required to get an acceptable control simulation. The major accomplishment is the validation of this new high resolution configuration of CCSM3. There are major improvements in our simulation of the surface wind stress and sea ice thickness distribution in the Arctic. Surface wind stress and oceanmore » circulation in the Antarctic Circumpolar Current are also improved. Our results demonstrate that the FV version of the CCSM coupled model is a state of the art climate model whose simulation capabilities are in the class of those used for IPCC assessments. We have also provided 1000 years of model data to Scripps Institution of Oceanography to estimate the natural variability of stream flow in California. In the future, our global model simulations will provide boundary data to high-resolution mesoscale model that will be used at LLNL. The mesoscale model would dynamically downscale the GCM climate to regional scale on climate time scales.« less

Monitoring interannual variation in global crop yield using long-term AVHRR and MODIS observations

NASA Astrophysics Data System (ADS)

Zhang, Xiaoyang; Zhang, Qingyuan

2016-04-01

Advanced Very High Resolution Radiometer (AVHRR) and Moderate Resolution Imaging Spectroradiometer (MODIS) data have been extensively applied for crop yield prediction because of their daily temporal resolution and a global coverage. This study investigated global crop yield using daily two band Enhanced Vegetation Index (EVI2) derived from AVHRR (1981-1999) and MODIS (2000-2013) observations at a spatial resolution of 0.05° (∼5 km). Specifically, EVI2 temporal trajectory of crop growth was simulated using a hybrid piecewise logistic model (HPLM) for individual pixels, which was used to detect crop phenological metrics. The derived crop phenology was then applied to calculate crop greenness defined as EVI2 amplitude and EVI2 integration during annual crop growing seasons, which was further aggregated for croplands in each country, respectively. The interannual variations in EVI2 amplitude and EVI2 integration were combined to correlate to the variation in cereal yield from 1982-2012 for individual countries using a stepwise regression model, respectively. The results show that the confidence level of the established regression models was higher than 90% (P value < 0.1) in most countries in the northern hemisphere although it was relatively poor in the southern hemisphere (mainly in Africa). The error in the yield predication was relatively smaller in America, Europe and East Asia than that in Africa. In the 10 countries with largest cereal production across the world, the prediction error was less than 9% during past three decades. This suggests that crop phenology-controlled greenness from coarse resolution satellite data has the capability of predicting national crop yield across the world, which could provide timely and reliable crop information for global agricultural trade and policymakers.

A Variable-Resolution Stretched-Grid General Circulation Model and Data Assimilation System with Multiple Areas of Interest: Studying the Anomalous Regional Climate Events of 1998

NASA Technical Reports Server (NTRS)

Fox-Rabinovitz, Michael S.; Takacs, Lawrence; Govindaraju, Ravi C.; Atlas, Robert (Technical Monitor)

2002-01-01

The new stretched-grid design with multiple (four) areas of interest, one at each global quadrant, is implemented into both a stretched-grid GCM (general circulation model) and a stretched-grid data assimilation system (DAS). The four areas of interest include: the U.S./Northern Mexico, the El Nino area/Central South America, India/China, and the Eastern Indian Ocean/Australia. Both the stretched-grid GCM and DAS annual (November 1997 through December 1998) integrations are performed with 50 km regional resolution. The efficient regional down-scaling to mesoscales is obtained for each of the four areas of interest while the consistent interactions between regional and global scales and the high quality of global circulation, are preserved. This is the advantage of the stretched-grid approach. The global variable resolution DAS incorporating the stretched-grid GCM has been developed and tested as an efficient tool for producing regional analyses and diagnostics with enhanced mesoscale resolution. The anomalous regional climate events of 1998 that occurred over the U.S., Mexico, South America, China, India, African Sahel, and Australia are investigated in both simulation and data assimilation modes. Tree assimilated products are also used, along with gauge precipitation data, for validating the simulation results. The obtained results show that the stretched-grid GCM and DAS are capable of producing realistic high quality simulated and assimilated products at mesoscale resolution for regional climate studies and applications.

Quantifying the importance of spatial resolution and other factors through global sensitivity analysis of a flood inundation model

NASA Astrophysics Data System (ADS)

Thomas Steven Savage, James; Pianosi, Francesca; Bates, Paul; Freer, Jim; Wagener, Thorsten

2016-11-01

Where high-resolution topographic data are available, modelers are faced with the decision of whether it is better to spend computational resource on resolving topography at finer resolutions or on running more simulations to account for various uncertain input factors (e.g., model parameters). In this paper we apply global sensitivity analysis to explore how influential the choice of spatial resolution is when compared to uncertainties in the Manning's friction coefficient parameters, the inflow hydrograph, and those stemming from the coarsening of topographic data used to produce Digital Elevation Models (DEMs). We apply the hydraulic model LISFLOOD-FP to produce several temporally and spatially variable model outputs that represent different aspects of flood inundation processes, including flood extent, water depth, and time of inundation. We find that the most influential input factor for flood extent predictions changes during the flood event, starting with the inflow hydrograph during the rising limb before switching to the channel friction parameter during peak flood inundation, and finally to the floodplain friction parameter during the drying phase of the flood event. Spatial resolution and uncertainty introduced by resampling topographic data to coarser resolutions are much more important for water depth predictions, which are also sensitive to different input factors spatially and temporally. Our findings indicate that the sensitivity of LISFLOOD-FP predictions is more complex than previously thought. Consequently, the input factors that modelers should prioritize will differ depending on the model output assessed, and the location and time of when and where this output is most relevant.

High Resolution Modeling of Hurricanes in a Climate Context

NASA Astrophysics Data System (ADS)

Knutson, T. R.

2007-12-01

Modeling of tropical cyclone activity in a climate context initially focused on simulation of relatively weak tropical storm-like disturbances as resolved by coarse grid (200 km) global models. As computing power has increased, multi-year simulations with global models of grid spacing 20-30 km have become feasible. Increased resolution also allowed for simulation storms of increasing intensity, and some global models generate storms of hurricane strength, depending on their resolution and other factors, although detailed hurricane structure is not simulated realistically. Results from some recent high resolution global model studies are reviewed. An alternative for hurricane simulation is regional downscaling. An early approach was to embed an operational (GFDL) hurricane prediction model within a global model solution, either for 5-day case studies of particular model storm cases, or for "idealized experiments" where an initial vortex is inserted into an idealized environments derived from global model statistics. Using this approach, hurricanes up to category five intensity can be simulated, owing to the model's relatively high resolution (9 km grid) and refined physics. Variants on this approach have been used to provide modeling support for theoretical predictions that greenhouse warming will increase the maximum intensities of hurricanes. These modeling studies also simulate increased hurricane rainfall rates in a warmer climate. The studies do not address hurricane frequency issues, and vertical shear is neglected in the idealized studies. A recent development is the use of regional model dynamical downscaling for extended (e.g., season-length) integrations of hurricane activity. In a study for the Atlantic basin, a non-hydrostatic model with grid spacing of 18km is run without convective parameterization, but with internal spectral nudging toward observed large-scale (basin wavenumbers 0-2) atmospheric conditions from reanalyses. Using this approach, our model reproduces the observed increase in Atlantic hurricane activity (numbers, Accumulated Cyclone Energy (ACE), Power Dissipation Index (PDI), etc.) over the period 1980-2006 fairly realistically, and also simulates ENSO-related interannual variations in hurricane counts. Annual simulated hurricane counts from a two-member ensemble correlate with observed counts at r=0.86. However, the model does not simulate hurricanes as intense as those observed, with minimum central pressures of 937 hPa (category 4) and maximum surface winds of 47 m/s (category 2) being the most intense simulated so far in these experiments. To explore possible impacts of future climate warming on Atlantic hurricane activity, we are re-running the 1980- 2006 seasons, keeping the interannual to multidecadal variations unchanged, but altering the August-October mean climate according to changes simulated by an 18-member ensemble of AR4 climate models (years 2080- 2099, A1B emission scenario). The warmer climate state features higher Atlantic SSTs, and also increased vertical wind shear across the Caribbean (Vecchi and Soden, GRL 2007). A key assumption of this approach is that the 18-model ensemble-mean climate change is the best available projection of future climate change in the Atlantic. Some of the 18 global models show little increase in wind shear, or even a decrease, and thus there will be considerable uncertainty associated with the hurricane frequency results, which will require further exploration. Results from our simulations will be presented at the meeting.

TerraClimate, a high-resolution global dataset of monthly climate and climatic water balance from 1958-2015.

PubMed

Abatzoglou, John T; Dobrowski, Solomon Z; Parks, Sean A; Hegewisch, Katherine C

2018-01-09

We present TerraClimate, a dataset of high-spatial resolution (1/24°, ~4-km) monthly climate and climatic water balance for global terrestrial surfaces from 1958-2015. TerraClimate uses climatically aided interpolation, combining high-spatial resolution climatological normals from the WorldClim dataset, with coarser resolution time varying (i.e., monthly) data from other sources to produce a monthly dataset of precipitation, maximum and minimum temperature, wind speed, vapor pressure, and solar radiation. TerraClimate additionally produces monthly surface water balance datasets using a water balance model that incorporates reference evapotranspiration, precipitation, temperature, and interpolated plant extractable soil water capacity. These data provide important inputs for ecological and hydrological studies at global scales that require high spatial resolution and time varying climate and climatic water balance data. We validated spatiotemporal aspects of TerraClimate using annual temperature, precipitation, and calculated reference evapotranspiration from station data, as well as annual runoff from streamflow gauges. TerraClimate datasets showed noted improvement in overall mean absolute error and increased spatial realism relative to coarser resolution gridded datasets.

TerraClimate, a high-resolution global dataset of monthly climate and climatic water balance from 1958-2015

NASA Astrophysics Data System (ADS)

Abatzoglou, John T.; Dobrowski, Solomon Z.; Parks, Sean A.; Hegewisch, Katherine C.

2018-01-01

We present TerraClimate, a dataset of high-spatial resolution (1/24°, ~4-km) monthly climate and climatic water balance for global terrestrial surfaces from 1958-2015. TerraClimate uses climatically aided interpolation, combining high-spatial resolution climatological normals from the WorldClim dataset, with coarser resolution time varying (i.e., monthly) data from other sources to produce a monthly dataset of precipitation, maximum and minimum temperature, wind speed, vapor pressure, and solar radiation. TerraClimate additionally produces monthly surface water balance datasets using a water balance model that incorporates reference evapotranspiration, precipitation, temperature, and interpolated plant extractable soil water capacity. These data provide important inputs for ecological and hydrological studies at global scales that require high spatial resolution and time varying climate and climatic water balance data. We validated spatiotemporal aspects of TerraClimate using annual temperature, precipitation, and calculated reference evapotranspiration from station data, as well as annual runoff from streamflow gauges. TerraClimate datasets showed noted improvement in overall mean absolute error and increased spatial realism relative to coarser resolution gridded datasets.

A Prototype Nonhydrostatic Regional-to-Global Nested-Grid Atmosphere Model for Medium-range Weather Forecasting

NASA Astrophysics Data System (ADS)

Harris, L.; Lin, S. J.; Zhou, L.; Chen, J. H.; Benson, R.; Rees, S.

2016-12-01

Limited-area convection-permitting models have proven useful for short-range NWP, but are unable to interact with the larger scales needed for longer lead-time skill. A new global forecast model, fvGFS, has been designed combining a modern nonhydrostatic dynamical core, the GFDL Finite-Volume Cubed-Sphere dynamical core (FV3) with operational GFS physics and initial conditions, and has been shown to provide excellent global skill while improving representation of small-scale phenomena. The nested-grid capability of FV3 allows us to build a regional-to-global variable-resolution model to efficiently refine to 3-km grid spacing over the Continental US. The use of two-way grid nesting allows us to reach these resolutions very efficiently, with the operational requirement easily attainable on current supercomputing systems.Even without a boundary-layer or advanced microphysical scheme appropriate for convection-perrmitting resolutions, the effectiveness of fvGFS can be demonstrated for a variety of weather events. We demonstrate successful proof-of-concept simulations of a variety of phenomena. We show the capability to develop intense hurricanes with realistic fine-scale eyewalls and rainbands. The new model also produces skillful predictions of severe weather outbreaks and of organized mesoscale convective systems. Fine-scale orographic and boundary-layer phenomena are also simulated with excellent fidelity by fvGFS. Further expected improvements are discussed, including the introduction of more sophisticated microphysics and of scale-aware convection schemes.

High-resolution station-based diurnal ionospheric total electron content (TEC) from dual-frequency GPS observations

NASA Astrophysics Data System (ADS)

ćepni, Murat S.; Potts, Laramie V.; Miima, John B.

2013-09-01

electron content (TEC) estimates derived from Global Navigation Satellite System (GNSS) signal delays provide a rich source of information about the Earth's ionosphere. Networks of Global Positioning System (GPS) receivers data can be used to represent the ionosphere by a Global Ionospheric Map (GIM). Data input for GIMs is dual-frequency GNSS-only or a mixture of GNSS and altimetry observations. Parameterization of GNSS-only GIMs approaches the ionosphere as a single-layer model (SLM) to determine GPS TEC models over a region. Limitations in GNSS-only GIM TEC are due largely to the nonhomogenous global distribution of GPS tracking stations with large data gaps over the oceans. The utility of slant GPS ionospheric-induced path delays for high temporal resolution from a single-station data rate offers better representation of TEC over a small region. A station-based vertical TEC (TECV) approach modifies the traditional single-layer model (SLM) GPS TEC method by introducing a zenith angle weighting (ZAW) filter to capture signal delays from mostly near-zenith satellite passes. Comparison with GIMs shows the station-dependent TEC (SD-TEC) model exhibits robust performance under variable space weather conditions. The SD-TEC model was applied to investigate ionospheric TEC variability during the geomagnetic storm event of 9 March 2012 at midlatitude station NJJJ located in New Jersey, USA. The high temporal resolution TEC results suggest TEC production and loss rate differences before, during, and after the storm.

Global habitat suitability for framework-forming cold-water corals.

PubMed

Davies, Andrew J; Guinotte, John M

2011-04-15

Predictive habitat models are increasingly being used by conservationists, researchers and governmental bodies to identify vulnerable ecosystems and species' distributions in areas that have not been sampled. However, in the deep sea, several limitations have restricted the widespread utilisation of this approach. These range from issues with the accuracy of species presences, the lack of reliable absence data and the limited spatial resolution of environmental factors known or thought to control deep-sea species' distributions. To address these problems, global habitat suitability models have been generated for five species of framework-forming scleractinian corals by taking the best available data and using a novel approach to generate high resolution maps of seafloor conditions. High-resolution global bathymetry was used to resample gridded data from sources such as World Ocean Atlas to produce continuous 30-arc second (∼1 km(2)) global grids for environmental, chemical and physical data of the world's oceans. The increased area and resolution of the environmental variables resulted in a greater number of coral presence records being incorporated into habitat models and higher accuracy of model predictions. The most important factors in determining cold-water coral habitat suitability were depth, temperature, aragonite saturation state and salinity. Model outputs indicated the majority of suitable coral habitat is likely to occur on the continental shelves and slopes of the Atlantic, South Pacific and Indian Oceans. The North Pacific has very little suitable scleractinian coral habitat. Numerous small scale features (i.e., seamounts), which have not been sampled or identified as having a high probability of supporting cold-water coral habitat were identified in all ocean basins. Field validation of newly identified areas is needed to determine the accuracy of model results, assess the utility of modelling efforts to identify vulnerable marine ecosystems for inclusion in future marine protected areas and reduce coral bycatch by commercial fisheries.

Projecting the effects of climate change on Calanus finmarchicus distribution within the U.S. Northeast Continental Shelf.

PubMed

Grieve, Brian D; Hare, Jon A; Saba, Vincent S

2017-07-24

Calanus finmarchicus is vital to pelagic ecosystems in the North Atlantic Ocean. Previous studies suggest the species is vulnerable to the effects of global warming, particularly on the Northeast U.S. Shelf, which is in the southern portion of its range. In this study, we evaluate an ensemble of six different downscaled climate models and a high-resolution global climate model, and create a generalized additive model (GAM) to examine how future changes in temperature and salinity could affect the distribution and density of C. finmarchicus. By 2081-2100, we project average C. finmarchicus density will decrease by as much as 50% under a high greenhouse gas emissions scenario. These decreases are particularly pronounced in the spring and summer in the Gulf of Maine and Georges Bank. When compared to a high-resolution global climate model, the ensemble showed a more uniform change throughout the Northeast U.S. Shelf, while the high-resolution model showed larger decreases in the Northeast Channel, Shelf Break, and Central Gulf of Maine. C. finmarchicus is an important link between primary production and higher trophic levels, and the decrease projected here could be detrimental to the North Atlantic Right Whale and a host of important fishery species.

Sensitivities of the hydrologic cycle to model physics, grid resolution, and ocean type in the aquaplanet Community Atmosphere Model

NASA Astrophysics Data System (ADS)

Benedict, James J.; Medeiros, Brian; Clement, Amy C.; Pendergrass, Angeline G.

2017-06-01

Precipitation distributions and extremes play a fundamental role in shaping Earth's climate and yet are poorly represented in many global climate models. Here, a suite of idealized Community Atmosphere Model (CAM) aquaplanet simulations is examined to assess the aquaplanet's ability to reproduce hydroclimate statistics of real-Earth configurations and to investigate sensitivities of precipitation distributions and extremes to model physics, horizontal grid resolution, and ocean type. Little difference in precipitation statistics is found between aquaplanets using time-constant sea-surface temperatures and those implementing a slab ocean model with a 50 m mixed-layer depth. In contrast, CAM version 5.3 (CAM5.3) produces more time mean, zonally averaged precipitation than CAM version 4 (CAM4), while CAM4 generates significantly larger precipitation variance and frequencies of extremely intense precipitation events. The largest model configuration-based precipitation sensitivities relate to choice of horizontal grid resolution in the selected range 1-2°. Refining grid resolution has significant physics-dependent effects on tropical precipitation: for CAM4, time mean zonal mean precipitation increases along the Equator and the intertropical convergence zone (ITCZ) narrows, while for CAM5.3 precipitation decreases along the Equator and the twin branches of the ITCZ shift poleward. Increased grid resolution also reduces light precipitation frequencies and enhances extreme precipitation for both CAM4 and CAM5.3 resulting in better alignment with observational estimates. A discussion of the potential implications these hydrologic cycle sensitivities have on the interpretation of precipitation statistics in future climate projections is also presented.