Atmospheric gravity waves with small vertical-to-horizotal wavelength ratios

NASA Astrophysics Data System (ADS)

Song, I. S.; Jee, G.; Kim, Y. H.; Chun, H. Y.

2017-12-01

Gravity wave modes with small vertical-to-horizontal wavelength ratios of an order of 10-3 are investigated through the systematic scale analysis of governing equations for gravity wave perturbations embedded in the quasi-geostrophic large-scale flow. These waves can be categorized as acoustic gravity wave modes because their total energy is given by the sum of kinetic, potential, and elastic parts. It is found that these waves can be forced by density fluctuations multiplied by the horizontal gradients of the large-scale pressure (geopotential) fields. These theoretical findings are evaluated using the results of a high-resolution global model (Specified Chemistry WACCM with horizontal resolution of 25 km and vertical resolution of 600 m) by computing the density-related gravity-wave forcing terms from the modeling results.

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.

HIGH-RESOLUTION DATASET OF URBAN CANOPY PARAMETERS FOR HOUSTON, TEXAS

EPA Science Inventory

Urban dispersion and air quality simulation models applied at various horizontal scales require different levels of fidelity for specifying the characteristics of the underlying surfaces. As the modeling scales approach the neighborhood level (~1 km horizontal grid spacing), the...

The Impact of Horizontal and Temporal Resolution on Convection and Precipitation with High-Resolution GEOS-5

NASA Technical Reports Server (NTRS)

Putman, William P.

2012-01-01

Using a high-resolution non-hydrostatic version of GEOS-5 with the cubed-sphere finite-volume dynamical core, the impact of spatial and temporal resolution on cloud properties will be evaluated. There are indications from examining convective cluster development in high resolution GEOS-5 forecasts that the temporal resolution within the model may playas significant a role as horizontal resolution. Comparing modeled convective cloud clusters versus satellite observations of brightness temperature, we have found that improved. temporal resolution in GEOS-S accounts for a significant portion of the improvements in the statistical distribution of convective cloud clusters. Using satellite simulators in GEOS-S we will compare the cloud optical properties of GEOS-S at various spatial and temporal resolutions with those observed from MODIS. The potential impact of these results on tropical cyclone formation and intensity will be examined as well.

Assessment of the effects of horizontal grid resolution on long-term air quality trends using coupled WRF-CMAQ simulations

EPA Science Inventory

The objective of this study is to determine the adequacy of using a relatively coarse horizontal resolution (i.e. 36 km) to simulate long-term trends of pollutant concentrations and radiation variables with the coupled WRF-CMAQ model. WRF-CMAQ simulations over the continental Uni...

Mesosacle eddies in a high resolution OGCM and coupled ocean-atmosphere GCM

NASA Astrophysics Data System (ADS)

Yu, Y.; Liu, H.; Lin, P.

2017-12-01

The present study described high-resolution climate modeling efforts including oceanic, atmospheric and coupled general circulation model (GCM) at the state key laboratory of numerical modeling for atmospheric sciences and geophysical fluid dynamics (LASG), Institute of Atmospheric Physics (IAP). The high-resolution OGCM is established based on the latest version of the LASG/IAP Climate system Ocean Model (LICOM2.1), but its horizontal resolution and vertical resolution are increased to 1/10° and 55 layers, respectively. Forced by the surface fluxes from the reanalysis and observed data, the model has been integrated for approximately more than 80 model years. Compared with the simulation of the coarse-resolution OGCM, the eddy-resolving OGCM not only better simulates the spatial-temporal features of mesoscale eddies and the paths and positions of western boundary currents but also reproduces the large meander of the Kuroshio Current and its interannual variability. Another aspect, namely, the complex structures of equatorial Pacific currents and currents in the coastal ocean of China, are better captured due to the increased horizontal and vertical resolution. Then we coupled the high resolution OGCM to NCAR CAM4 with 25km resolution, in which the mesoscale air-sea interaction processes are better captured.

Evaluation of a multi-scale WRF-CAM5 simulation during the 2010 East Asian Summer Monsoon

DOE PAGES

Campbell, Patrick; Zhang, Yang; Wang, Kai; ...

2017-09-08

The Weather Research and Forecasting model with Chemistry (WRF-Chem) with the physics package of the Community Atmosphere Model Version 5 (CAM5) has been applied at multiple scales over Eastern China (EC) and the Yangtze River Delta (YRD) to evaluate how increased horizontal resolution with physics designed for a coarser resolution climate model impacts aerosols and clouds, and the resulting precipitation characteristics and performance during the 2010 East Asian Summer Monsoon (EASM). Despite large underpredictions in surface aerosol concentrations and aerosol optical depth, there is good spatial agreement with surface observations of chemical predictions, and increasing spatial resolution tends to improvemore » performance. Model bias and normalized root mean square values for precipitation predictions are relatively small, but there are significant differences when comparing modeled and observed probability density functions for precipitation in EC and YRD. Increasing model horizontal resolution tends to reduce model bias and error for precipitation predictions. The surface and column aerosol loading is maximized between about 32°N and 42°N in early to mid-May during the 2010 EASM, and then shifts north while decreasing in magnitude during July and August. Changing model resolution moderately changes the spatiotemporal relationships between aerosols, cloud properties, and precipitation during the EASM, thus demonstrating the importance of model grid resolution in simulating EASM circulation and rainfall patterns over EC and the YRD. In conclusion, results from this work demonstrate the capability and limitations in the aerosol, cloud, and precipitation representation of WRF-CAM5 for regional-scale applications down to relatively fine horizontal resolutions. Further WRF-CAM5 model development and application in this area is needed.« less

Evaluation of a multi-scale WRF-CAM5 simulation during the 2010 East Asian Summer Monsoon

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

Campbell, Patrick; Zhang, Yang; Wang, Kai

The Weather Research and Forecasting model with Chemistry (WRF-Chem) with the physics package of the Community Atmosphere Model Version 5 (CAM5) has been applied at multiple scales over Eastern China (EC) and the Yangtze River Delta (YRD) to evaluate how increased horizontal resolution with physics designed for a coarser resolution climate model impacts aerosols and clouds, and the resulting precipitation characteristics and performance during the 2010 East Asian Summer Monsoon (EASM). Despite large underpredictions in surface aerosol concentrations and aerosol optical depth, there is good spatial agreement with surface observations of chemical predictions, and increasing spatial resolution tends to improvemore » performance. Model bias and normalized root mean square values for precipitation predictions are relatively small, but there are significant differences when comparing modeled and observed probability density functions for precipitation in EC and YRD. Increasing model horizontal resolution tends to reduce model bias and error for precipitation predictions. The surface and column aerosol loading is maximized between about 32N and 42N in early to mid-May during the 2010 EASM, and then shifts north while decreasing in magnitude during July and August. Changing model resolution moderately changes the spatiotemporal relationships between aerosols, cloud properties, and precipitation during the EASM, thus demonstrating the importance of model grid resolution in simulating EASM circulation and rainfall patterns over EC and the YRD. Results from this work demonstrate the capability and limitations in the aerosol, cloud, and precipitation representation of WRF-CAM5 for regional-scale applications down to relatively fine horizontal resolutions. Further WRF-CAM5 model development and application in this area is needed.« less

Evaluation of a multi-scale WRF-CAM5 simulation during the 2010 East Asian Summer Monsoon

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

Campbell, Patrick; Zhang, Yang; Wang, Kai

The Weather Research and Forecasting model with Chemistry (WRF-Chem) with the physics package of the Community Atmosphere Model Version 5 (CAM5) has been applied at multiple scales over Eastern China (EC) and the Yangtze River Delta (YRD) to evaluate how increased horizontal resolution with physics designed for a coarser resolution climate model impacts aerosols and clouds, and the resulting precipitation characteristics and performance during the 2010 East Asian Summer Monsoon (EASM). Despite large underpredictions in surface aerosol concentrations and aerosol optical depth, there is good spatial agreement with surface observations of chemical predictions, and increasing spatial resolution tends to improvemore » performance. Model bias and normalized root mean square values for precipitation predictions are relatively small, but there are significant differences when comparing modeled and observed probability density functions for precipitation in EC and YRD. Increasing model horizontal resolution tends to reduce model bias and error for precipitation predictions. The surface and column aerosol loading is maximized between about 32°N and 42°N in early to mid-May during the 2010 EASM, and then shifts north while decreasing in magnitude during July and August. Changing model resolution moderately changes the spatiotemporal relationships between aerosols, cloud properties, and precipitation during the EASM, thus demonstrating the importance of model grid resolution in simulating EASM circulation and rainfall patterns over EC and the YRD. In conclusion, results from this work demonstrate the capability and limitations in the aerosol, cloud, and precipitation representation of WRF-CAM5 for regional-scale applications down to relatively fine horizontal resolutions. Further WRF-CAM5 model development and application in this area is needed.« less

The influence of spatially and temporally high-resolution wind forcing on the power input to near-inertial waves in the ocean

NASA Astrophysics Data System (ADS)

Rimac, Antonija; von Storch, Jin-Song; Eden, Carsten

2013-04-01

The estimated power required to sustain global general circulation in the ocean is about 2 TW. This power is supplied with wind stress and tides. Energy spectrum shows pronounced maxima at near-inertial frequency. Near-inertial waves excited by high-frequency winds represent an important source for deep ocean mixing since they can propagate into the deep ocean and dissipate far away from the generation sites. The energy input by winds to near-inertial waves has been studied mostly using slab ocean models and wind stress forcing with coarse temporal resolution (e.g. 6-hourly). Slab ocean models lack the ability to reproduce fundamental aspects of kinetic energy balance and systematically overestimate the wind work. Also, slab ocean models do not account the energy used for the mixed layer deepening or the energy radiating downward into the deep ocean. Coarse temporal resolution of the wind forcing strongly underestimates the near-inertial energy. To overcome this difficulty we use an eddy permitting ocean model with high-frequency wind forcing. We establish the following model setup: We use the Max Planck Institute Ocean Model (MPIOM) on a tripolar grid with 45 km horizontal resolution and 40 vertical levels. We run the model with wind forcings that vary in horizontal and temporal resolution. We use high-resolution (1-hourly with 35 km horizontal resolution) and low-resolution winds (6-hourly with 250 km horizontal resolution). We address the following questions: Is the kinetic energy of near-inertial waves enhanced when high-resolution wind forcings are used? If so, is this due to higher level of overall wind variability or higher spatial or temporal resolution of wind forcing? How large is the power of near-inertial waves generated by winds? Our results show that near-inertial waves are enhanced and the near-inertial kinetic energy is two times higher (in the storm track regions 3.5 times higher) when high-resolution winds are used. Filtering high-resolution winds in space and time, the near-inertial kinetic energy reduces. The reduction is faster when a temporal filter is used suggesting that the high-frequency wind forcing is more efficient in generating near-inertial wave energy than the small-scale wind forcing. Using low-resolution wind forcing the wind generated power to near-inertial waves is 0.55 TW. When we use high-resolution wind forcing the result is 1.6 TW meaning that the result increases by 300%.

EXAMINATION OF MODEL PREDICTIONS AT DIFFERENT HORIZONTAL GRID RESOLUTIONS

EPA Science Inventory

While fluctuations in meteorological and air quality variables occur on a continuum of spatial scales, the horizontal grid spacing of coupled meteorological and photochemical models sets a lower limit on the spatial scales that they can resolve. However, both computational costs ...

Evaluation of the Operational Multi-scale Environment model with Grid Adaptivity (OMEGA) for use in Wind Energy Applications in the Great Basin of Nevada

NASA Astrophysics Data System (ADS)

King, Kristien C.

In order to further assess the wind energy potential for Nevada, the accuracy of a computational meteorological model, the Operational Multi-scale Environment model with Grid Adaptivity (OMEGA), was evaluated by comparing simulation results with data collected from a wind monitoring tower near Tonopah, NV. The state of Nevada is characterized by high mountains and low-lying valleys, therefore, in order to determine the wind potential for the state, meteorological models that predict the wind must be able to accurately represent and account for terrain features and simulate topographic forcing with accuracy. Topographic forcing has a dominant role in the development and modification of mesoscale flows in regions of complex terrain, like Tonopah, especially at the level of wind turbine blade heights (~80 m). Additionally, model factors such as horizontal resolution, terrain database resolution, model physics, time of model initialization, stability regime, and source of initial conditions may each affect the ability of a mesoscale model to forecast winds correctly. The observational tower used for comparison was located at Stone Cabin, Nevada. The tower had both sonic anemometers and cup anemometers installed at heights of 40 m, 60 m, and 80 m above the surface. During a previous experiment, tower data were collected for the period February 9 through March 10, 2007 and compared to model simulations using the MM5 and WRF models at a number of varying horizontal resolutions. In this previous research, neither the MM5 nor the WRF showed a significant improvement in ability to forecast wind speed with increasing horizontal grid resolution. The present research evaluated the ability of OMEGA to reproduce point winds as compared to the observational data from the Stone Cabin Tower at heights of 40 m, 60 m, and 80 m. Unlike other mesoscale atmospheric models, OMEGA incorporates an unstructured triangular adaptive grid which allows for increased flexibility and accuracy in characterizing areas of complex terrain. Model sensitivity to horizontal grid resolution, initial conditions, and time of initialization were tested. OMEGA was run over three different horizontal grid resolutions with minimum horizontal edge lengths of: 18 km, 6 km, and 2 km. For each resolution, the model was initialized using both the Global Forecasting System (GFS) and North American Regional Reanalysis (NARR) to determine model sensitivity to initial conditions. For both the NARR and GFS initializations, the model was started at both 0000 UTC and 1200 UTC to determine the effect of start time and stability regime on the performance of the model. An additional intensive study into the model's performance was also conducted by a detailed evaluation of model results during two separate 24-hour periods, the first a period where the model performed well and the second a period where the model performed poorly, to determine which atmospheric factors most affect the predictive ability of the OMEGA model. The statistical results were then compared with the results from the MM5 and WRF simulations to determine the most appropriate model for wind energy potential studies in complex terrain.

Resolution dependence of precipitation statistical fidelity in hindcast simulations

DOE PAGES

O'Brien, Travis A.; Collins, William D.; Kashinath, Karthik; ...

2016-06-19

This article is a U.S. Government work and is in the public domain in the USA. Numerous studies have shown that atmospheric models with high horizontal resolution better represent the physics and statistics of precipitation in climate models. While it is abundantly clear from these studies that high-resolution increases the rate of extreme precipitation, it is not clear whether these added extreme events are “realistic”; whether they occur in simulations in response to the same forcings that drive similar events in reality. In order to understand whether increasing horizontal resolution results in improved model fidelity, a hindcast-based, multiresolution experimental designmore » has been conceived and implemented: the InitiaLIzed-ensemble, Analyze, and Develop (ILIAD) framework. The ILIAD framework allows direct comparison between observed and simulated weather events across multiple resolutions and assessment of the degree to which increased resolution improves the fidelity of extremes. Analysis of 5 years of daily 5 day hindcasts with the Community Earth System Model at horizontal resolutions of 220, 110, and 28 km shows that: (1) these hindcasts reproduce the resolution-dependent increase of extreme precipitation that has been identified in longer-duration simulations, (2) the correspondence between simulated and observed extreme precipitation improves as resolution increases; and (3) this increase in extremes and precipitation fidelity comes entirely from resolved-scale precipitation. Evidence is presented that this resolution-dependent increase in precipitation intensity can be explained by the theory of Rauscher et al. (), which states that precipitation intensifies at high resolution due to an interaction between the emergent scaling (spectral) properties of the wind field and the constraint of fluid continuity.« less

Resolution dependence of precipitation statistical fidelity in hindcast simulations

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

O'Brien, Travis A.; Collins, William D.; Kashinath, Karthik

This article is a U.S. Government work and is in the public domain in the USA. Numerous studies have shown that atmospheric models with high horizontal resolution better represent the physics and statistics of precipitation in climate models. While it is abundantly clear from these studies that high-resolution increases the rate of extreme precipitation, it is not clear whether these added extreme events are “realistic”; whether they occur in simulations in response to the same forcings that drive similar events in reality. In order to understand whether increasing horizontal resolution results in improved model fidelity, a hindcast-based, multiresolution experimental designmore » has been conceived and implemented: the InitiaLIzed-ensemble, Analyze, and Develop (ILIAD) framework. The ILIAD framework allows direct comparison between observed and simulated weather events across multiple resolutions and assessment of the degree to which increased resolution improves the fidelity of extremes. Analysis of 5 years of daily 5 day hindcasts with the Community Earth System Model at horizontal resolutions of 220, 110, and 28 km shows that: (1) these hindcasts reproduce the resolution-dependent increase of extreme precipitation that has been identified in longer-duration simulations, (2) the correspondence between simulated and observed extreme precipitation improves as resolution increases; and (3) this increase in extremes and precipitation fidelity comes entirely from resolved-scale precipitation. Evidence is presented that this resolution-dependent increase in precipitation intensity can be explained by the theory of Rauscher et al. (), which states that precipitation intensifies at high resolution due to an interaction between the emergent scaling (spectral) properties of the wind field and the constraint of fluid continuity.« less

Data for Figures and Tables in Journal Article Assessment of the Effects of Horizontal Grid Resolution on Long-Term Air Quality Trends using Coupled WRF-CMAQ Simulations, doi:10.1016/j.atmosenv.2016.02.036

EPA Pesticide Factsheets

The dataset represents the data depicted in the Figures and Tables of a Journal Manuscript with the following abstract: The objective of this study is to determine the adequacy of using a relatively coarse horizontal resolution (i.e. 36 km) to simulate long-term trends of pollutant concentrations and radiation variables with the coupled WRF-CMAQ model. WRF-CMAQ simulations over the continental United State are performed over the 2001 to 2010 time period at two different horizontal resolutions of 12 and 36 km. Both simulations used the same emission inventory and model configurations. Model results are compared both in space and time to assess the potential weaknesses and strengths of using coarse resolution in long-term air quality applications. The results show that the 36 km and 12 km simulations are comparable in terms of trends analysis for both pollutant concentrations and radiation variables. The advantage of using the coarser 36 km resolution is a significant reduction of computational cost, time and storage requirement which are key considerations when performing multiple years of simulations for trend analysis. However, if such simulations are to be used for local air quality analysis, finer horizontal resolution may be beneficial since it can provide information on local gradients. In particular, divergences between the two simulations are noticeable in urban, complex terrain and coastal regions.This dataset is associated with the following publication

Baroclinic mixing of potential vorticity as the principal sharpening mechanism for the extratropical Tropopause Inversion Layer

NASA Astrophysics Data System (ADS)

Wang, Shu Meir; Geller, Marvin A.

2016-09-01

Previous works have shown that a dry, idealized general circulation model could produce many features of the extratropical Tropopause Inversion Layer (TIL). In particular, the following have been shown, but no explanations were given for these results. (1) A sharper extratropical TIL resulted more from increased horizontal resolution than from increased vertical resolution. (2) If the Equator-to-Pole temperature gradient was varied, the annual variation of the extratropical TIL found in observations could be reproduced. (3) The extratropical TIL altitude showed excellent correlation with the upper tropospheric relative vorticity, as had been previously proposed. (4) Increased horizontal model resolutions led to extratropical TILs that were at lower altitudes. We show that these conclusions follow from baroclinic mixing of high stratospheric potential vorticity into the troposphere being the principal sharpening mechanism for the extratropical TIL and the increased baroclinic activity occurring in higher horizontal resolution models. We furthermore suggest that the distance from the jet exerts a greater influence on the height and sharpness of the extratropical TIL than does the upper tropospheric relative vorticity, and this accounts for the annual behavior of the extratropical TIL found in observations and reproduced with a dry, mechanistic, global model.

Impact of Variable-Resolution Meshes on Regional Climate Simulations

NASA Astrophysics Data System (ADS)

Fowler, L. D.; Skamarock, W. C.; Bruyere, C. L.

2014-12-01

The Model for Prediction Across Scales (MPAS) is currently being used for seasonal-scale simulations on globally-uniform and regionally-refined meshes. Our ongoing research aims at analyzing simulations of tropical convective activity and tropical cyclone development during one hurricane season over the North Atlantic Ocean, contrasting statistics obtained with a variable-resolution mesh against those obtained with a quasi-uniform mesh. Analyses focus on the spatial distribution, frequency, and intensity of convective and grid-scale precipitations, and their relative contributions to the total precipitation as a function of the horizontal scale. Multi-month simulations initialized on May 1st 2005 using ERA-Interim re-analyses indicate that MPAS performs satisfactorily as a regional climate model for different combinations of horizontal resolutions and transitions between the coarse and refined meshes. Results highlight seamless transitions for convection, cloud microphysics, radiation, and land-surface processes between the quasi-uniform and locally- refined meshes, despite the fact that the physics parameterizations were not developed for variable resolution meshes. Our goal of analyzing the performance of MPAS is twofold. First, we want to establish that MPAS can be successfully used as a regional climate model, bypassing the need for nesting and nudging techniques at the edges of the computational domain as done in traditional regional climate modeling. Second, we want to assess the performance of our convective and cloud microphysics parameterizations as the horizontal resolution varies between the lower-resolution quasi-uniform and higher-resolution locally-refined areas of the global domain.

Impact of Variable-Resolution Meshes on Regional Climate Simulations

NASA Astrophysics Data System (ADS)

Fowler, L. D.; Skamarock, W. C.; Bruyere, C. L.

2013-12-01

The Model for Prediction Across Scales (MPAS) is currently being used for seasonal-scale simulations on globally-uniform and regionally-refined meshes. Our ongoing research aims at analyzing simulations of tropical convective activity and tropical cyclone development during one hurricane season over the North Atlantic Ocean, contrasting statistics obtained with a variable-resolution mesh against those obtained with a quasi-uniform mesh. Analyses focus on the spatial distribution, frequency, and intensity of convective and grid-scale precipitations, and their relative contributions to the total precipitation as a function of the horizontal scale. Multi-month simulations initialized on May 1st 2005 using NCEP/NCAR re-analyses indicate that MPAS performs satisfactorily as a regional climate model for different combinations of horizontal resolutions and transitions between the coarse and refined meshes. Results highlight seamless transitions for convection, cloud microphysics, radiation, and land-surface processes between the quasi-uniform and locally-refined meshes, despite the fact that the physics parameterizations were not developed for variable resolution meshes. Our goal of analyzing the performance of MPAS is twofold. First, we want to establish that MPAS can be successfully used as a regional climate model, bypassing the need for nesting and nudging techniques at the edges of the computational domain as done in traditional regional climate modeling. Second, we want to assess the performance of our convective and cloud microphysics parameterizations as the horizontal resolution varies between the lower-resolution quasi-uniform and higher-resolution locally-refined areas of the global domain.

Alleviating tropical Atlantic sector biases in the Kiel climate model by enhancing horizontal and vertical atmosphere model resolution: climatology and interannual variability

NASA Astrophysics Data System (ADS)

Harlaß, Jan; Latif, Mojib; Park, Wonsun

2018-04-01

We investigate the quality of simulating tropical Atlantic (TA) sector climatology and interannual variability in integrations of the Kiel climate model (KCM) with varying atmosphere model resolution. The ocean model resolution is kept fixed. A reasonable simulation of TA sector annual-mean climate, seasonal cycle and interannual variability can only be achieved at sufficiently high horizontal and vertical atmospheric resolution. Two major reasons for the improvements are identified. First, the western equatorial Atlantic westerly surface wind bias in spring can be largely eliminated, which is explained by a better representation of meridional and especially vertical zonal momentum transport. The enhanced atmospheric circulation along the equator in turn greatly improves the thermal structure of the upper equatorial Atlantic with much reduced warm sea surface temperature (SST) biases. Second, the coastline in the southeastern TA and steep orography are better resolved at high resolution, which improves wind structure and in turn reduces warm SST biases in the Benguela upwelling region. The strongly diminished wind and SST biases at high atmosphere model resolution allow for a more realistic latitudinal position of the intertropical convergence zone. Resulting stronger cross-equatorial winds, in conjunction with a shallower thermocline, enable a rapid cold tongue development in the eastern TA in boreal spring. This enables simulation of realistic interannual SST variability and its seasonal phase locking in the KCM, which primarily is the result of a stronger thermocline feedback. Our findings suggest that enhanced atmospheric resolution, both vertical and horizontal, could be a key to achieving more realistic simulation of TA climatology and interannual variability in climate models.

High-resolution dust modelling over complex terrains in West Asia

NASA Astrophysics Data System (ADS)

Basart, S.; Vendrell, L.; Baldasano, J. M.

2016-12-01

The present work demonstrates the impact of model resolution in dust propagation in a complex terrain region such as West Asia. For this purpose, two simulations using the NMMB/BSC-Dust model are performed and analysed, one with a high horizontal resolution (at 0.03° × 0.03°) and one with a lower horizontal resolution (at 0.33° × 0.33°). Both model experiments cover two intense dust storms that occurred on 17-20 March 2012 as a consequence of strong northwesterly Shamal winds that spanned over thousands of kilometres in West Asia. The comparison with ground-based (surface weather stations and sunphotometers) and satellite aerosol observations (Aqua/MODIS and MSG/SEVIRI) shows that despite differences in the magnitude of the simulated dust concentrations, the model is able to reproduce these two dust outbreaks. Differences between both simulations on the dust spread rise on regional dust transport areas in south-western Saudi Arabia, Yemen and Oman. The complex orography in south-western Saudi Arabia, Yemen and Oman (with peaks higher than 3000 m) has an impact on the transported dust concentration fields over mountain regions. Differences between both model configurations are mainly associated to the channelization of the dust flow through valleys and the differences in the modelled altitude of the mountains that alters the meteorology and blocks the dust fronts limiting the dust transport. These results demonstrate how the dust prediction in the vicinity of complex terrains improves using high-horizontal resolution simulations.

Air quality modelling in the Berlin-Brandenburg region using WRF-Chem v3.7.1: sensitivity to resolution of model grid and input data

NASA Astrophysics Data System (ADS)

Kuik, Friderike; Lauer, Axel; Churkina, Galina; Denier van der Gon, Hugo A. C.; Fenner, Daniel; Mar, Kathleen A.; Butler, Tim M.

2016-12-01

Air pollution is the number one environmental cause of premature deaths in Europe. Despite extensive regulations, air pollution remains a challenge, especially in urban areas. For studying summertime air quality in the Berlin-Brandenburg region of Germany, the Weather Research and Forecasting Model with Chemistry (WRF-Chem) is set up and evaluated against meteorological and air quality observations from monitoring stations as well as from a field campaign conducted in 2014. The objective is to assess which resolution and level of detail in the input data is needed for simulating urban background air pollutant concentrations and their spatial distribution in the Berlin-Brandenburg area. The model setup includes three nested domains with horizontal resolutions of 15, 3 and 1 km and anthropogenic emissions from the TNO-MACC III inventory. We use RADM2 chemistry and the MADE/SORGAM aerosol scheme. Three sensitivity simulations are conducted updating input parameters to the single-layer urban canopy model based on structural data for Berlin, specifying land use classes on a sub-grid scale (mosaic option) and downscaling the original emissions to a resolution of ca. 1 km × 1 km for Berlin based on proxy data including traffic density and population density. The results show that the model simulates meteorology well, though urban 2 m temperature and urban wind speeds are biased high and nighttime mixing layer height is biased low in the base run with the settings described above. We show that the simulation of urban meteorology can be improved when specifying the input parameters to the urban model, and to a lesser extent when using the mosaic option. On average, ozone is simulated reasonably well, but maximum daily 8 h mean concentrations are underestimated, which is consistent with the results from previous modelling studies using the RADM2 chemical mechanism. Particulate matter is underestimated, which is partly due to an underestimation of secondary organic aerosols. NOx (NO + NO2) concentrations are simulated reasonably well on average, but nighttime concentrations are overestimated due to the model's underestimation of the mixing layer height, and urban daytime concentrations are underestimated. The daytime underestimation is improved when using downscaled, and thus locally higher emissions, suggesting that part of this bias is due to deficiencies in the emission input data and their resolution. The results further demonstrate that a horizontal resolution of 3 km improves the results and spatial representativeness of the model compared to a horizontal resolution of 15 km. With the input data (land use classes, emissions) at the level of detail of the base run of this study, we find that a horizontal resolution of 1 km does not improve the results compared to a resolution of 3 km. However, our results suggest that a 1 km horizontal model resolution could enable a detailed simulation of local pollution patterns in the Berlin-Brandenburg region if the urban land use classes, together with the respective input parameters to the urban canopy model, are specified with a higher level of detail and if urban emissions of higher spatial resolution are used.

Increasing horizontal resolution in numerical weather prediction and climate simulations: illusion or panacea?

PubMed

Wedi, Nils P

2014-06-28

The steady path of doubling the global horizontal resolution approximately every 8 years in numerical weather prediction (NWP) at the European Centre for Medium Range Weather Forecasts may be substantially altered with emerging novel computing architectures. It coincides with the need to appropriately address and determine forecast uncertainty with increasing resolution, in particular, when convective-scale motions start to be resolved. Blunt increases in the model resolution will quickly become unaffordable and may not lead to improved NWP forecasts. Consequently, there is a need to accordingly adjust proven numerical techniques. An informed decision on the modelling strategy for harnessing exascale, massively parallel computing power thus also requires a deeper understanding of the sensitivity to uncertainty--for each part of the model--and ultimately a deeper understanding of multi-scale interactions in the atmosphere and their numerical realization in ultra-high-resolution NWP and climate simulations. This paper explores opportunities for substantial increases in the forecast efficiency by judicious adjustment of the formal accuracy or relative resolution in the spectral and physical space. One path is to reduce the formal accuracy by which the spectral transforms are computed. The other pathway explores the importance of the ratio used for the horizontal resolution in gridpoint space versus wavenumbers in spectral space. This is relevant for both high-resolution simulations as well as ensemble-based uncertainty estimation. © 2014 The Author(s) Published by the Royal Society. All rights reserved.

National Centers for Environmental Prediction

Science.gov Websites

resolution at T574 becomes ~ 23 km T382 Spectral truncation equivalent to horizontal resolution ~37 km T254 Spectral truncation equivalent to horizontal resolution ~50-55 km T190 Spectral truncation equivalent to horizontal resolution ~70 km T126 Spectral truncation equivalent to horizontal resolution ~100 km UM Unified

Overflow Simulations using MPAS-Ocean in Idealized and Realistic Domains

NASA Astrophysics Data System (ADS)

Reckinger, S.; Petersen, M. R.; Reckinger, S. J.

2016-02-01

MPAS-Ocean is used to simulate an idealized, density-driven overflow using the dynamics of overflow mixing and entrainment (DOME) setup. Numerical simulations are benchmarked against other models, including the MITgcm's z-coordinate model and HIM's isopycnal coordinate model. A full parameter study is presented that looks at how sensitive overflow simulations are to vertical grid type, resolution, and viscosity. Horizontal resolutions with 50 km grid cells are under-resolved and produce poor results, regardless of other parameter settings. Vertical grids ranging in thickness from 15 m to 120 m were tested. A horizontal resolution of 10 km and a vertical resolution of 60 m are sufficient to resolve the mesoscale dynamics of the DOME configuration, which mimics real-world overflow parameters. Mixing and final buoyancy are least sensitive to horizontal viscosity, but strongly sensitive to vertical viscosity. This suggests that vertical viscosity could be adjusted in overflow water formation regions to influence mixing and product water characteristics. Also, the study shows that sigma coordinates produce much less mixing than z-type coordinates, resulting in heavier plumes that go further down slope. Sigma coordinates are less sensitive to changes in resolution but as sensitive to vertical viscosity compared to z-coordinates. Additionally, preliminary measurements of overflow diagnostics on global simulations using a realistic oceanic domain are presented.

Sensitivity of sea-level forecasting to the horizontal resolution and sea surface forcing for different configurations of an oceanographic model of the Adriatic Sea

NASA Astrophysics Data System (ADS)

Bressan, Lidia; Valentini, Andrea; Paccagnella, Tiziana; Montani, Andrea; Marsigli, Chiara; Stefania Tesini, Maria

2017-04-01

At the Hydro-meteo-climate service of the Regional environmental agency of Emilia-Romagna, Italy (Arpae-SIMC), the oceanographic numerical model AdriaROMS is used in the operational forecasting suite to compute sea level, temperature, salinity and 3-D current fields of the Adriatic Sea (northern Mediterranean Sea). In order to evaluate the performance of the sea-level forecast and to study different configurations of the ROMS model, two marine storms occurred on the Emilia Romagna coast during the winter 2015-2016 are investigated. The main focus of this study is to analyse the sensitivity of the model to the horizontal resolution and to the meteorological forcing. To this end, the model is run with two different configurations and with two horizontal grids at 1 and 2 km resolution. To study the influence of the meteorological forcing, the two storms have been reproduced by running ROMS in ensemble mode, forced by the 16-members of the meteorological ensemble COSMO-LEPS system. Possible optimizations of the model set-up are deduced by the comparison of the different run outputs.

Simulation of the Summer Monsoon Rainfall over East Asia using the NCEP GFS Cumulus Parameterization at Different Horizontal Resolutions

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

Lim, Kyo-Sun; Hong, Song You; Yoon, Jin-Ho

2014-10-01

The most recent version of Simplified Arakawa-Schubert (SAS) cumulus scheme in National Center for Environmental Prediction (NCEP) Global Forecast System (GFS) (GFS SAS) has been implemented into the Weather and Research Forecasting (WRF) model with a modification of triggering condition and convective mass flux to become depending on model’s horizontal grid spacing. East Asian Summer Monsoon of 2006 from June to August is selected to evaluate the performance of the modified GFS SAS scheme. Simulated monsoon rainfall with the modified GFS SAS scheme shows better agreement with observation compared to the original GFS SAS scheme. The original GFS SAS schememore » simulates the similar ratio of subgrid-scale precipitation, which is calculated from a cumulus scheme, against total precipitation regardless of model’s horizontal grid spacing. This is counter-intuitive because the portion of resolved clouds in a grid box should be increased as the model grid spacing decreases. This counter-intuitive behavior of the original GFS SAS scheme is alleviated by the modified GFS SAS scheme. Further, three different cumulus schemes (Grell and Freitas, Kain and Fritsch, and Betts-Miller-Janjic) are chosen to investigate the role of a horizontal resolution on simulated monsoon rainfall. The performance of high-resolution modeling is not always enhanced as the spatial resolution becomes higher. Even though improvement of probability density function of rain rate and long wave fluxes by the higher-resolution simulation is robust regardless of a choice of cumulus parameterization scheme, the overall skill score of surface rainfall is not monotonically increasing with spatial resolution.« less

On the role of horizontal resolution over the Tibetan Plateau in the REMO regional climate model

NASA Astrophysics Data System (ADS)

Xu, Jingwei; Koldunov, Nikolay; Remedio, Armelle Reca C.; Sein, Dmitry V.; Zhi, Xiefei; Jiang, Xi; Xu, Min; Zhu, Xiuhua; Fraedrich, Klaus; Jacob, Daniela

2018-02-01

A number of studies have shown that added value is obtained by increasing the horizontal resolution of a regional climate model to capture additional fine-scale weather processes. However, the mechanisms leading to this added value are different over areas with complicated orographic features, such as the Tibetan Plateau (TP). To determine the role that horizontal resolution plays over the TP, a detailed comparison was made between the results from the REMO regional climate model at resolutions of 25 and 50 km for the period 1980-2007. The model was driven at the lateral boundaries by the European Centre for Medium-Range Weather Forecasts Interim Reanalysis data. The experiments differ only in representation of topography, all other land parameters (e.g., vegetation characteristics, soil texture) are the same. The results show that the high-resolution topography affects the regional air circulation near the ground surface around the edge of the TP, which leads to a redistribution of the transport of atmospheric water vapor, especially over the Brahmaputra and Irrawaddy valleys—the main water vapor paths for the southern TP—increasing the amount of atmospheric water vapor transported onto the TP by about 5%. This, in turn, significantly decreases the temperature at 2 m by > 1.5 °C in winter in the high-resolution simulation of the southern TP. The impact of topography on the 2 m temperature over the TP is therefore by influencing the transport of atmospheric water vapor in the main water vapor paths.

Land use and topography influence in a complex terrain area: A high resolution mesoscale modelling study over the Eastern Pyrenees using the WRF model

NASA Astrophysics Data System (ADS)

Jiménez-Esteve, B.; Udina, M.; Soler, M. R.; Pepin, N.; Miró, J. R.

2018-04-01

Different types of land use (LU) have different physical properties which can change local energy balance and hence vertical fluxes of moisture, heat and momentum. This in turn leads to changes in near-surface temperature and moisture fields. Simulating atmospheric flow over complex terrain requires accurate local-scale energy balance and therefore model grid spacing must be sufficient to represent both topography and land-use. In this study we use both the Corine Land Cover (CLC) and United States Geological Survey (USGS) land use databases for use with the Weather Research and Forecasting (WRF) model and evaluate the importance of both land-use classification and horizontal resolution in contributing to successful modelling of surface temperatures and humidities observed from a network of 39 sensors over a 9 day period in summer 2013. We examine case studies of the effects of thermal inertia and soil moisture availability at individual locations. The scale at which the LU classification is observed influences the success of the model in reproducing observed patterns of temperature and moisture. Statistical validation of model output demonstrates model sensitivity to both the choice of LU database used and the horizontal resolution. In general, results show that on average, by a) using CLC instead of USGS and/or b) increasing horizontal resolution, model performance is improved. We also show that the sensitivity to these changes in the model performance shows a daily cycle.

Evaluating hourly rainfall characteristics over the U.S. Great Plains in dynamically downscaled climate model simulations using NASA-Unified WRF

NASA Astrophysics Data System (ADS)

Lee, Huikyo; Waliser, Duane E.; Ferraro, Robert; Iguchi, Takamichi; Peters-Lidard, Christa D.; Tian, Baijun; Loikith, Paul C.; Wright, Daniel B.

2017-07-01

Accurate simulation of extreme precipitation events remains a challenge in climate models. This study utilizes hourly precipitation data from ground stations and satellite instruments to evaluate rainfall characteristics simulated by the NASA-Unified Weather Research and Forecasting (NU-WRF) regional climate model at horizontal resolutions of 4, 12, and 24 km over the Great Plains of the United States. We also examined the sensitivity of the simulated precipitation to different spectral nudging approaches and the cumulus parameterizations. The rainfall characteristics in the observations and simulations were defined as an hourly diurnal cycle of precipitation and a joint probability distribution function (JPDF) between duration and peak intensity of precipitation events over the Great Plains in summer. We calculated a JPDF for each data set and the overlapping area between observed and simulated JPDFs to measure the similarity between the two JPDFs. Comparison of the diurnal precipitation cycles between observations and simulations does not reveal the added value of high-resolution simulations. However, the performance of NU-WRF simulations measured by the JPDF metric strongly depends on horizontal resolution. The simulation with the highest resolution of 4 km shows the best agreement with the observations in simulating duration and intensity of wet spells. Spectral nudging does not affect the JPDF significantly. The effect of cumulus parameterizations on the JPDFs is considerable but smaller than that of horizontal resolution. The simulations with lower resolutions of 12 and 24 km show reasonable agreement but only with the high-resolution observational data that are aggregated into coarse resolution and spatially averaged.

A statistical analysis of the dependency of closure assumptions in cumulus parameterization on the horizontal resolution

NASA Technical Reports Server (NTRS)

Xu, Kuan-Man

1994-01-01

Simulated data from the UCLA cumulus ensemble model are used to investigate the quasi-universal validity of closure assumptions used in existing cumulus parameterizations. A closure assumption is quasi-universally valid if it is sensitive neither to convective cloud regimes nor to horizontal resolutions of large-scale/mesoscale models. The dependency of three types of closure assumptions, as classified by Arakawa and Chen, on the horizontal resolution is addressed in this study. Type I is the constraint on the coupling of the time tendencies of large-scale temperature and water vapor mixing ratio. Type II is the constraint on the coupling of cumulus heating and cumulus drying. Type III is a direct constraint on the intensity of a cumulus ensemble. The macroscopic behavior of simulated cumulus convection is first compared with the observed behavior in view of Type I and Type II closure assumptions using 'quick-look' and canonical correlation analyses. It is found that they are statistically similar to each other. The three types of closure assumptions are further examined with simulated data averaged over selected subdomain sizes ranging from 64 to 512 km. It is found that the dependency of Type I and Type II closure assumptions on the horizontal resolution is very weak and that Type III closure assumption is somewhat dependent upon the horizontal resolution. The influences of convective and mesoscale processes on the closure assumptions are also addressed by comparing the structures of canonical components with the corresponding vertical profiles in the convective and stratiform regions of cumulus ensembles analyzed directly from simulated data. The implication of these results for cumulus parameterization is discussed.

Satellite observed thermodynamics during FGGE

NASA Technical Reports Server (NTRS)

Smith, W. L.

1985-01-01

During the First Global Atmospheric Research Program (GARP) Global Experiment (FGGE), determinations of temperature and moisture were made from TIROS-N and NOAA-6 satellite infrared and microwave sounding radiance measurements. The data were processed by two methods differing principally in their horizontal resolution. At the National Earth Satellite Service (NESS) in Washington, D.C., the data were produced operationally with a horizontal resolution of 250 km for inclusion in the FGGE Level IIb data sets for application to large-scale numerical analysis and prediction models. High horizontal resolution (75 km) sounding data sets were produced using man-machine interactive methods for the special observing periods of FGGE at the NASA/Goddard Space Flight Center and archived as supplementary Level IIb. The procedures used for sounding retrieval and the characteristics and quality of these thermodynamic observations are given.

A study of overflow simulations using MPAS-Ocean: Vertical grids, resolution, and viscosity

NASA Astrophysics Data System (ADS)

Reckinger, Shanon M.; Petersen, Mark R.; Reckinger, Scott J.

2015-12-01

MPAS-Ocean is used to simulate an idealized, density-driven overflow using the dynamics of overflow mixing and entrainment (DOME) setup. Numerical simulations are carried out using three of the vertical coordinate types available in MPAS-Ocean, including z-star with partial bottom cells, z-star with full cells, and sigma coordinates. The results are first benchmarked against other models, including the MITgcm's z-coordinate model and HIM's isopycnal coordinate model, which are used to set the base case used for this work. A full parameter study is presented that looks at how sensitive overflow simulations are to vertical grid type, resolution, and viscosity. Horizontal resolutions with 50 km grid cells are under-resolved and produce poor results, regardless of other parameter settings. Vertical grids ranging in thickness from 15 m to 120 m were tested. A horizontal resolution of 10 km and a vertical resolution of 60 m are sufficient to resolve the mesoscale dynamics of the DOME configuration, which mimics real-world overflow parameters. Mixing and final buoyancy are least sensitive to horizontal viscosity, but strongly sensitive to vertical viscosity. This suggests that vertical viscosity could be adjusted in overflow water formation regions to influence mixing and product water characteristics. Lastly, the study shows that sigma coordinates produce much less mixing than z-type coordinates, resulting in heavier plumes that go further down slope. Sigma coordinates are less sensitive to changes in resolution but as sensitive to vertical viscosity compared to z-coordinates.

An Idealized Test of the Response of the Community Atmosphere Model to Near-Grid-Scale Forcing Across Hydrostatic Resolutions

NASA Astrophysics Data System (ADS)

Herrington, A. R.; Reed, K. A.

2018-02-01

A set of idealized experiments are developed using the Community Atmosphere Model (CAM) to understand the vertical velocity response to reductions in forcing scale that is known to occur when the horizontal resolution of the model is increased. The test consists of a set of rising bubble experiments, in which the horizontal radius of the bubble and the model grid spacing are simultaneously reduced. The test is performed with moisture, through incorporating moist physics routines of varying complexity, although convection schemes are not considered. Results confirm that the vertical velocity in CAM is to first-order, proportional to the inverse of the horizontal forcing scale, which is consistent with a scale analysis of the dry equations of motion. In contrast, experiments in which the coupling time step between the moist physics routines and the dynamical core (i.e., the "physics" time step) are relaxed back to more conventional values results in severely damped vertical motion at high resolution, degrading the scaling. A set of aqua-planet simulations using different physics time steps are found to be consistent with the results of the idealized experiments.

A meteorological distribution system for high-resolution terrestrial modeling (MicroMet)

Treesearch

Glen E. Liston; Kelly Elder

2006-01-01

An intermediate-complexity, quasi-physically based, meteorological model (MicroMet) has been developed to produce high-resolution (e.g., 30-m to 1-km horizontal grid increment) atmospheric forcings required to run spatially distributed terrestrial models over a wide variety of landscapes. The following eight variables, required to run most terrestrial models, are...

A two-dimensional kinematic dynamo model of the ionospheric magnetic field at Venus

NASA Technical Reports Server (NTRS)

Cravens, T. E.; Wu, D.; Shinagawa, H.

1990-01-01

The results of a high-resolution, two-dimensional, time dependent, kinematic dynamo model of the ionospheric magnetic field of Venus are presented. Various one-dimensional models are considered and the two-dimensional model is then detailed. In this model, the two-dimensional magnetic induction equation, the magnetic diffusion-convection equation, is numerically solved using specified plasma velocities. Origins of the vertical velocity profile and of the horizontal velocities are discussed. It is argued that the basic features of the vertical magnetic field profile remain unaltered by horizontal flow effects and also that horizontal plasma flow can strongly affect the magnetic field for altitudes above 300 km.

Challenges in the development of very high resolution Earth System Models for climate science

NASA Astrophysics Data System (ADS)

Rasch, Philip J.; Xie, Shaocheng; Ma, Po-Lun; Lin, Wuyin; Wan, Hui; Qian, Yun

2017-04-01

The authors represent the 20+ members of the ACME atmosphere development team. The US Department of Energy (DOE) has, like many other organizations around the world, identified the need for an Earth System Model capable of rapid completion of decade to century length simulations at very high (vertical and horizontal) resolution with good climate fidelity. Two years ago DOE initiated a multi-institution effort called ACME (Accelerated Climate Modeling for Energy) to meet this an extraordinary challenge, targeting a model eventually capable of running at 10-25km horizontal and 20-400m vertical resolution through the troposphere on exascale computational platforms at speeds sufficient to complete 5+ simulated years per day. I will outline the challenges our team has encountered in development of the atmosphere component of this model, and the strategies we have been using for tuning and debugging a model that we can barely afford to run on today's computational platforms. These strategies include: 1) evaluation at lower resolutions; 2) ensembles of short simulations to explore parameter space, and perform rough tuning and evaluation; 3) use of regionally refined versions of the model for probing high resolution model behavior at less expense; 4) use of "auto-tuning" methodologies for model tuning; and 5) brute force long climate simulations.

Increased horizontal viewing zone angle of a hologram by resolution redistribution of a spatial light modulator.

PubMed

Takaki, Yasuhiro; Hayashi, Yuki

2008-07-01

The narrow viewing zone angle is one of the problems associated with electronic holography. We propose a technique that enables the ratio of horizontal and vertical resolutions of a spatial light modulator (SLM) to be altered. This technique increases the horizontal resolution of a SLM several times, so that the horizontal viewing zone angle is also increased several times. A SLM illuminated by a slanted point light source array is imaged by a 4f imaging system in which a horizontal slit is located on the Fourier plane. We show that the horizontal resolution was increased four times and that the horizontal viewing zone angle was increased approximately four times.

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

Smagorinsky-type diffusion in a high-resolution GCM

NASA Astrophysics Data System (ADS)

Schaefer-Rolffs, Urs; Becker, Erich

2013-04-01

The parametrization of the (horizontal) momentum diffusion is a paramount component of a Global Circulation Model (GCM). Aside from friction in the boundary layer, a relevant fraction of kinetic energy is dissipated in the free atmosphere, and it is known that a linear harmonic turbulence model is not sufficient to obtain a reasonable simulation of the kinetic energy spectrum. Therefore, often empirical hyper-diffusion schemes are employed, regardless of disadvantages like the violation of energy conservation and the second law of thermodynamics. At IAP we have developed an improved parametrization of the horizontal diffusion that is based on Smagorinsky's nonlinear and energy conservation formulation. This approach is extended by the dynamic Smagorinsky model (DSM) of M. Germano. In this new scheme, the mixing length is no longer a prescribed parameter but calculated dynamically from the resolved flow such as to preserve scale invariance for the horizontal energy cascade. The so-called Germano identity is solved by a tensor norm ansatz which yields a positive definite frictional heating. We present results from an investigation using the DSM as a parametrization of horizontal diffusion in a high-resolution version of the Kühlungborn Mechanistic general Circulation Model (KMCM) with spectral truncation at horizontal wavenumber 330. The DSM calculates the Smagorinsky parameter cS independent from the resolution scale. We find that this method yields an energy spectrum that exhibits a pronounced transition from a synoptic -3 to a mesoscale -5-3 slope at wavenumbers around 50. At the highest wavenumber end, a behaviour similar to that often obtained by tuning the hyper-diffusion is achieved self-consistently. This result is very sensitive to the explicit choice of the test filter in the DSM.

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.

Effects of Increased Horizontal Resolution on Simulation of the North American Monsoon in the NCAR CAM3: An Evaluation based on Surface, Satellite, and Reanalysis Data

NASA Astrophysics Data System (ADS)

Collier, J. C.; Zhang, G. J.

2006-05-01

Simulation of the North American monsoon system by the National Center for Atmospheric Research (NCAR) Community Atmosphere Model (CAM3) is evaluated in its sensitivity to increasing horizontal resolution. For two resolutions, T42 and T85, rainfall is compared to TRMM satellite-derived and surface gauge-based rainfall rates over the U.S. and northern Mexico as well as rainfall accumulations in gauges of the North American Monsoon Experiment (NAME) Enhanced Rain Gauge Network (NERN) in the Sierra Madre Occidental mountains. Simulated upper-tropospheric mass and wind fields are compared to those from NCEP-NCAR reanalyses. The comparison presented herein demonstrates that tropospheric motions associated with the North American monsoon system are sensitive to increasing the horizontal resolution of the model. An increase in resolution from T42 to T85 results in changes to a region of large-scale mid-tropospheric descent found north and east of the monsoon anticyclone. Relative to its simulation at T42, this region extends farther south and west at T85. Additionally, at T85, the subsidence is stronger. Consistent with the differences in large-scale descent, the T85 simulation of CAM3 is anomalously dry over Texas and northeastern Mexico during the peak monsoon months. Meanwhile, the geographic distribution of rainfall over the Sierra Madre Occidental region of Mexico is more satisfactorily simulated at T85 than at T42 for July and August. Moisture import into this region is greater at T85 than at T42 during these months. A focused study of the Sierra Madre Occidental region in particular shows that, in the regional average sense, the timing of the peak of the monsoon is relatively insensitive to the horizontal resolution of the model, while a phase bias in the diurnal cycle of monsoon-season precipitation is somewhat reduced in the higher-resolution run. At both resolutions, CAM3 poorly simulates the month-to-month evolution of monsoon rainfall over extreme northwestern Mexico and Arizona, though biases are considerably improved at T85.

Computer Simulation of Global Profiles of Carbon Dioxide Using a Pulsed, 2-Micron, Coherent-Detection, Column-Content DIAL System

NASA Technical Reports Server (NTRS)

Kavaya, Michael J.; Singh, Upendra N.; Koch, Grady J.; Yu, Jirong; Frehlich, Rod G.

2009-01-01

We present preliminary results of computer simulations of the error in measuring carbon dioxide mixing ratio profiles from earth orbit. The simulated sensor is a pulsed, 2-micron, coherent-detection lidar alternately operating on at least two wavelengths. The simulated geometry is a nadir viewing lidar measuring the column content signal. Atmospheric absorption is modeled using FASCODE3P software with the HITRAN 2004 absorption line data base. Lidar shot accumulation is employed up to the horizontal resolution limit. Horizontal resolutions of 50, 100, and 200 km are shown. Assuming a 400 km spacecraft orbit, the horizontal resolutions correspond to measurement times of about 7, 14, and 28 s. We simulate laser pulse-pair repetition frequencies from 1 Hz to 100 kHz. The range of shot accumulation is 7 to 2.8 million pulse-pairs. The resultant error is shown as a function of horizontal resolution, laser pulse-pair repetition frequency, and laser pulse energy. The effect of different on and off pulse energies is explored. The results are compared to simulation results of others and to demonstrated 2-micron operating points at NASA Langley.

High Resolution Model Intercomparison Project (HighResMIP v1.0) for CMIP6

NASA Astrophysics Data System (ADS)

Haarsma, Reindert J.; Roberts, Malcolm J.; Vidale, Pier Luigi; Senior, Catherine A.; Bellucci, Alessio; Bao, Qing; Chang, Ping; Corti, Susanna; Fučkar, Neven S.; Guemas, Virginie; von Hardenberg, Jost; Hazeleger, Wilco; Kodama, Chihiro; Koenigk, Torben; Leung, L. Ruby; Lu, Jian; Luo, Jing-Jia; Mao, Jiafu; Mizielinski, Matthew S.; Mizuta, Ryo; Nobre, Paulo; Satoh, Masaki; Scoccimarro, Enrico; Semmler, Tido; Small, Justin; von Storch, Jin-Song

2016-11-01

Robust projections and predictions of climate variability and change, particularly at regional scales, rely on the driving processes being represented with fidelity in model simulations. The role of enhanced horizontal resolution in improved process representation in all components of the climate system is of growing interest, particularly as some recent simulations suggest both the possibility of significant changes in large-scale aspects of circulation as well as improvements in small-scale processes and extremes. However, such high-resolution global simulations at climate timescales, with resolutions of at least 50 km in the atmosphere and 0.25° in the ocean, have been performed at relatively few research centres and generally without overall coordination, primarily due to their computational cost. Assessing the robustness of the response of simulated climate to model resolution requires a large multi-model ensemble using a coordinated set of experiments. The Coupled Model Intercomparison Project 6 (CMIP6) is the ideal framework within which to conduct such a study, due to the strong link to models being developed for the CMIP DECK experiments and other model intercomparison projects (MIPs). Increases in high-performance computing (HPC) resources, as well as the revised experimental design for CMIP6, now enable a detailed investigation of the impact of increased resolution up to synoptic weather scales on the simulated mean climate and its variability. The High Resolution Model Intercomparison Project (HighResMIP) presented in this paper applies, for the first time, a multi-model approach to the systematic investigation of the impact of horizontal resolution. A coordinated set of experiments has been designed to assess both a standard and an enhanced horizontal-resolution simulation in the atmosphere and ocean. The set of HighResMIP experiments is divided into three tiers consisting of atmosphere-only and coupled runs and spanning the period 1950-2050, with the possibility of extending to 2100, together with some additional targeted experiments. This paper describes the experimental set-up of HighResMIP, the analysis plan, the connection with the other CMIP6 endorsed MIPs, as well as the DECK and CMIP6 historical simulations. HighResMIP thereby focuses on one of the CMIP6 broad questions, "what are the origins and consequences of systematic model biases?", but we also discuss how it addresses the World Climate Research Program (WCRP) grand challenges.

Sensitivity of mesoscale-model forecast skill to some initial-data characteristics, data density, data position, analysis procedure and measurement error

NASA Technical Reports Server (NTRS)

Warner, Thomas T.; Key, Lawrence E.; Lario, Annette M.

1989-01-01

The effects of horizontal and vertical data resolution, data density, data location, different objective analysis algorithms, and measurement error on mesoscale-forecast accuracy are studied with observing-system simulation experiments. Domain-averaged errors are shown to generally decrease with time. It is found that the vertical distribution of error growth depends on the initial vertical distribution of the error itself. Larger gravity-inertia wave noise is produced in forecasts with coarser vertical data resolution. The use of a low vertical resolution observing system with three data levels leads to more forecast errors than moderate and high vertical resolution observing systems with 8 and 14 data levels. Also, with poor vertical resolution in soundings, the initial and forecast errors are not affected by the horizontal data resolution.

Impact of Resolution on the Representation of Precipitation Variability Associated With the ITCZ

NASA Astrophysics Data System (ADS)

De Benedetti, Marc; Moore, G. W. K.

2017-12-01

The Intertropical Convergence Zone (ITCZ) is responsible for most of the weather and climate in equatorial regions along with many tropical-midlatitude interactions. It is therefore important to understand how models represent its structure and variability. Most ITCZ-associated precipitation is convective, making it unclear how horizontal resolution impacts its representation. To assess this, we introduce a novel technique that involves calculation of the precipitation field's decorrelation length scale (DCLS) using model data sets that share a common lineage with horizontal resolutions from 16 to 160 km. All resolutions captured the ITCZ's mean structure; however, imprints of topography, such as Hawaii and sea surface temperature, including the variability associated with upwelling cold water off the coast of South America, are more clearly represented at higher resolutions. The DCLS analysis indicates that there are changes in the spatial variability of the ITCZ's precipitation that are not reflected in its mean structure, thus confirming its utility as a diagnostic.

Scale dependency of regional climate modeling of current and future climate extremes in Germany

NASA Astrophysics Data System (ADS)

Tölle, Merja H.; Schefczyk, Lukas; Gutjahr, Oliver

2017-11-01

A warmer climate is projected for mid-Europe, with less precipitation in summer, but with intensified extremes of precipitation and near-surface temperature. However, the extent and magnitude of such changes are associated with creditable uncertainty because of the limitations of model resolution and parameterizations. Here, we present the results of convection-permitting regional climate model simulations for Germany integrated with the COSMO-CLM using a horizontal grid spacing of 1.3 km, and additional 4.5- and 7-km simulations with convection parameterized. Of particular interest is how the temperature and precipitation fields and their extremes depend on the horizontal resolution for current and future climate conditions. The spatial variability of precipitation increases with resolution because of more realistic orography and physical parameterizations, but values are overestimated in summer and over mountain ridges in all simulations compared to observations. The spatial variability of temperature is improved at a resolution of 1.3 km, but the results are cold-biased, especially in summer. The increase in resolution from 7/4.5 km to 1.3 km is accompanied by less future warming in summer by 1 ∘C. Modeled future precipitation extremes will be more severe, and temperature extremes will not exclusively increase with higher resolution. Although the differences between the resolutions considered (7/4.5 km and 1.3 km) are small, we find that the differences in the changes in extremes are large. High-resolution simulations require further studies, with effective parameterizations and tunings for different topographic regions. Impact models and assessment studies may benefit from such high-resolution model results, but should account for the impact of model resolution on model processes and climate change.

Evaluation of predicted diurnal cycle of precipitation after tests with convection and microphysics schemes in the Eta Model

NASA Astrophysics Data System (ADS)

Gomes, J. L.; Chou, S. C.; Yaguchi, S. M.

2012-04-01

Physics parameterizations and the model vertical and horizontal resolutions, for example, can significantly contribute to the uncertainty in the numerical weather predictions, especially at regions with complex topography. The objective of this study is to assess the influences of model precipitation production schemes and horizontal resolution on the diurnal cycle of precipitation in the Eta Model . The model was run in hydrostatic mode at 3- and 5-km grid sizes, the vertical resolution was set to 50 layers, and the time steps to 6 and 10 s, respectively. The initial and boundary conditions were taken from ERA-Interim reanalysis. Over the sea the 0.25-deg sea surface temperature from NOAA was used. The model was setup to run for each resolution over Angra dos Reis, located in the Southeast region of Brazil, for the rainy period between 18 December 2009 and 01 de January 2010, the model simulation range was 48 hours. In one set of runs the cumulus parameterization was switched off, in this case the model precipitation was fully simulated by cloud microphysics scheme, and in the other set the model was run with weak cumulus convection. The results show that as the model horizontal resolution increases from 5 to 3 km, the spatial pattern of the precipitation hardly changed, although the maximum precipitation core increased in magnitude. Daily data from automatic station data was used to evaluate the runs and shows that the diurnal cycle of temperature and precipitation were better simulated for 3 km when compared against observations. The model configuration results without cumulus convection shows a small contraction in the precipitating area and an increase in the simulated maximum values. The diurnal cycle of precipitation was better simulated with some activity of the cumulus convection scheme. The skill scores for the period and for different forecast ranges are higher at weak and moderate precipitation rates.

Numerical simulations of significant orographic precipitation in Madeira island

NASA Astrophysics Data System (ADS)

Couto, Flavio Tiago; Ducrocq, Véronique; Salgado, Rui; Costa, Maria João

2016-03-01

High-resolution simulations of high precipitation events with the MESO-NH model are presented, and also used to verify that increasing horizontal resolution in zones of complex orography, such as in Madeira island, improve the simulation of the spatial distribution and total precipitation. The simulations succeeded in reproducing the general structure of the cloudy systems over the ocean in the four periods considered of significant accumulated precipitation. The accumulated precipitation over the Madeira was better represented with the 0.5 km horizontal resolution and occurred under four distinct synoptic situations. Different spatial patterns of the rainfall distribution over the Madeira have been identified.

Modelling Precipitation and Temperature Extremes: The Importance of Horizontal Resolution

NASA Astrophysics Data System (ADS)

Shields, C. A.; Kiehl, J. T.; Meehl, G. A.

2013-12-01

Understanding Earth's water cycle on a warming planet is of critical importance in society's ability to adapt to climate change. Extreme weather events, such as floods, heat waves, and drought will likely change with the water cycle as greenhouse gases continue to rise. Location, duration, and intensity of extreme events can be studied using complex earth system models. Here, we employ the fully coupled Community Earth System Model (CESM1.0) to evaluate extreme event impacts for different possible future forcing scenarios. Simulations applying the Representative Concentration Pathway (RCP) scenarios 2.6 and 8.5 were chosen to bracket the range of model responses. Because extreme weather events happen on a regional scale, there is a tendency to favor using higher resolution models, i.e. models that can represent regional features with greater accuracy. Within the CESM1.0 framework, we evaluate both the standard 1 degree resolution (1 degree atmosphere/land coupled to 1 degree ocean/sea ice), and the higher 0.5 degree resolution version (0.5 degree atmosphere/land coupled to 1 degree ocean/sea ice), focusing on extreme precipitation events, heat waves, and droughts. We analyze a variety of geographical regions, but generally find that benefits from increased horizontal resolution are most significant on the regional scale.

Generation of a pseudo-2D shear-wave velocity section by inversion of a series of 1D dispersion curves

USGS Publications Warehouse

Luo, Y.; Xia, J.; Liu, J.; Xu, Y.; Liu, Q.

2008-01-01

Multichannel Analysis of Surface Waves utilizes a multichannel recording system to estimate near-surface shear (S)-wave velocities from high-frequency Rayleigh waves. A pseudo-2D S-wave velocity (vS) section is constructed by aligning 1D models at the midpoint of each receiver spread and using a spatial interpolation scheme. The horizontal resolution of the section is therefore most influenced by the receiver spread length and the source interval. The receiver spread length sets the theoretical lower limit and any vS structure with its lateral dimension smaller than this length will not be properly resolved in the final vS section. A source interval smaller than the spread length will not improve the horizontal resolution because spatial smearing has already been introduced by the receiver spread. In this paper, we first analyze the horizontal resolution of a pair of synthetic traces. Resolution analysis shows that (1) a pair of traces with a smaller receiver spacing achieves higher horizontal resolution of inverted S-wave velocities but results in a larger relative error; (2) the relative error of the phase velocity at a high frequency is smaller than at a low frequency; and (3) a relative error of the inverted S-wave velocity is affected by the signal-to-noise ratio of data. These results provide us with a guideline to balance the trade-off between receiver spacing (horizontal resolution) and accuracy of the inverted S-wave velocity. We then present a scheme to generate a pseudo-2D S-wave velocity section with high horizontal resolution using multichannel records by inverting high-frequency surface-wave dispersion curves calculated through cross-correlation combined with a phase-shift scanning method. This method chooses only a pair of consecutive traces within a shot gather to calculate a dispersion curve. We finally invert surface-wave dispersion curves of synthetic and real-world data. Inversion results of both synthetic and real-world data demonstrate that inverting high-frequency surface-wave dispersion curves - by a pair of traces through cross-correlation with phase-shift scanning method and with the damped least-square method and the singular-value decomposition technique - can feasibly achieve a reliable pseudo-2D S-wave velocity section with relatively high horizontal resolution. ?? 2008 Elsevier B.V. All rights reserved.

Regional model simulations of New Zealand climate

NASA Astrophysics Data System (ADS)

Renwick, James A.; Katzfey, Jack J.; Nguyen, Kim C.; McGregor, John L.

1998-03-01

Simulation of New Zealand climate is examined through the use of a regional climate model nested within the output of the Commonwealth Scientific and Industrial Research Organisation nine-level general circulation model (GCM). R21 resolution GCM output is used to drive a regional model run at 125 km grid spacing over the Australasian region. The 125 km run is used in turn to drive a simulation at 50 km resolution over New Zealand. Simulations with a full seasonal cycle are performed for 10 model years. The focus is on the quality of the simulation of present-day climate, but results of a doubled-CO2 run are discussed briefly. Spatial patterns of mean simulated precipitation and surface temperatures improve markedly as horizontal resolution is increased, through the better resolution of the country's orography. However, increased horizontal resolution leads to a positive bias in precipitation. At 50 km resolution, simulated frequency distributions of daily maximum/minimum temperatures are statistically similar to those of observations at many stations, while frequency distributions of daily precipitation appear to be statistically different to those of observations at most stations. Modeled daily precipitation variability at 125 km resolution is considerably less than observed, but is comparable to, or exceeds, observed variability at 50 km resolution. The sensitivity of the simulated climate to changes in the specification of the land surface is discussed briefly. Spatial patterns of the frequency of extreme temperatures and precipitation are generally well modeled. Under a doubling of CO2, the frequency of precipitation extremes changes only slightly at most locations, while air frosts become virtually unknown except at high-elevation sites.

Impact of numerical choices on water conservation in the E3SM Atmosphere Model Version 1 (EAM V1)

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

Zhang, Kai; Rasch, Philip J.; Taylor, Mark A.

The conservation of total water is an important numerical feature for global Earth system models. Even small conservation problems in the water budget can lead to systematic errors in century-long simulations for sea level rise projection. This study quantifies and reduces various sources of water conservation error in the atmosphere component of the Energy Exascale Earth System Model. Several sources of water conservation error have been identified during the development of the version 1 (V1) model. The largest errors result from the numerical coupling between the resolved dynamics and the parameterized sub-grid physics. A hybrid coupling using different methods formore » fluid dynamics and tracer transport provides a reduction of water conservation error by a factor of 50 at 1° horizontal resolution as well as consistent improvements at other resolutions. The second largest error source is the use of an overly simplified relationship between the surface moisture flux and latent heat flux at the interface between the host model and the turbulence parameterization. This error can be prevented by applying the same (correct) relationship throughout the entire model. Two additional types of conservation error that result from correcting the surface moisture flux and clipping negative water concentrations can be avoided by using mass-conserving fixers. With all four error sources addressed, the water conservation error in the V1 model is negligible and insensitive to the horizontal resolution. The associated changes in the long-term statistics of the main atmospheric features are small. A sensitivity analysis is carried out to show that the magnitudes of the conservation errors decrease strongly with temporal resolution but increase with horizontal resolution. The increased vertical resolution in the new model results in a very thin model layer at the Earth’s surface, which amplifies the conservation error associated with the surface moisture flux correction. We note that for some of the identified error sources, the proposed fixers are remedies rather than solutions to the problems at their roots. Future improvements in time integration would be beneficial for this model.« less

Tests of high-resolution simulations over a region of complex terrain in Southeast coast of Brazil

NASA Astrophysics Data System (ADS)

Chou, Sin Chan; Luís Gomes, Jorge; Ristic, Ivan; Mesinger, Fedor; Sueiro, Gustavo; Andrade, Diego; Lima-e-Silva, Pedro Paulo

2013-04-01

The Eta Model is used operationally by INPE at the Centre for Weather Forecasts and Climate Studies (CPTEC) to produce weather forecasts over South America since 1997. The model has gone through upgrades along these years. In order to prepare the model for operational higher resolution forecasts, the model is configured and tested over a region of complex topography located near the coast of Southeast Brazil. The model domain includes the two Brazilians cities, Rio de Janeiro and Sao Paulo, urban areas, preserved tropical forest, pasture fields, and complex terrain where it can rise from sea level up to about 1000 m. Accurate near-surface wind direction and magnitude are needed for the power plant emergency plan. Besides, the region suffers from frequent events of floods and landslides, therefore accurate local forecasts are required for disaster warnings. The objective of this work is to carry out a series of numerical experiments to test and evaluate high resolution simulations in this complex area. Verification of model runs uses observations taken from the nuclear power plant and higher resolution reanalyses data. The runs were tested in a period when flow was predominately forced by local conditions and in a period forced by frontal passage. The Eta Model was configured initially with 2-km horizontal resolution and 50 layers. The Eta-2km is a second nesting, it is driven by Eta-15km, which in its turn is driven by Era-Interim reanalyses. The series of experiments consists of replacing surface layer stability function, adjusting cloud microphysics scheme parameters, further increasing vertical and horizontal resolutions. By replacing the stability function for the stable conditions substantially increased the katabatic winds and verified better against the tower wind data. Precipitation produced by the model was excessive in the region. Increasing vertical resolution to 60 layers caused a further increase in precipitation production. This excessive precipitation was reduced by adjusting some parameters in the cloud microphysics scheme. Precipitation overestimate still occurs and further tests are still necessary. The increase of horizontal resolution to 1 km required adjusting model diffusion parameters and refining divergence calculations. Available observations in the region for a thorough evaluation is a major constraint.

Global 7 km mesh nonhydrostatic Model Intercomparison Project for improving TYphoon forecast (TYMIP-G7): experimental design and preliminary results

NASA Astrophysics Data System (ADS)

Nakano, Masuo; Wada, Akiyoshi; Sawada, Masahiro; Yoshimura, Hiromasa; Onishi, Ryo; Kawahara, Shintaro; Sasaki, Wataru; Nasuno, Tomoe; Yamaguchi, Munehiko; Iriguchi, Takeshi; Sugi, Masato; Takeuchi, Yoshiaki

2017-03-01

Recent advances in high-performance computers facilitate operational numerical weather prediction by global hydrostatic atmospheric models with horizontal resolutions of ˜ 10 km. Given further advances in such computers and the fact that the hydrostatic balance approximation becomes invalid for spatial scales < 10 km, the development of global nonhydrostatic models with high accuracy is urgently required. The Global 7 km mesh nonhydrostatic Model Intercomparison Project for improving TYphoon forecast (TYMIP-G7) is designed to understand and statistically quantify the advantages of high-resolution nonhydrostatic global atmospheric models to improve tropical cyclone (TC) prediction. A total of 137 sets of 5-day simulations using three next-generation nonhydrostatic global models with horizontal resolutions of 7 km and a conventional hydrostatic global model with a horizontal resolution of 20 km were run on the Earth Simulator. The three 7 km mesh nonhydrostatic models are the nonhydrostatic global spectral atmospheric Double Fourier Series Model (DFSM), the Multi-Scale Simulator for the Geoenvironment (MSSG) and the Nonhydrostatic ICosahedral Atmospheric Model (NICAM). The 20 km mesh hydrostatic model is the operational Global Spectral Model (GSM) of the Japan Meteorological Agency. Compared with the 20 km mesh GSM, the 7 km mesh models reduce systematic errors in the TC track, intensity and wind radii predictions. The benefits of the multi-model ensemble method were confirmed for the 7 km mesh nonhydrostatic global models. While the three 7 km mesh models reproduce the typical axisymmetric mean inner-core structure, including the primary and secondary circulations, the simulated TC structures and their intensities in each case are very different for each model. In addition, the simulated track is not consistently better than that of the 20 km mesh GSM. These results suggest that the development of more sophisticated initialization techniques and model physics is needed to further improve the TC prediction.

Hyper-Resolution Groundwater Modeling using MODFLOW 6

NASA Astrophysics Data System (ADS)

Hughes, J. D.; Langevin, C.

2017-12-01

MODFLOW 6 is the latest version of the U.S. Geological Survey's modular hydrologic model. MODFLOW 6 was developed to synthesize many of the recent versions of MODFLOW into a single program, improve the way different process models are coupled, and to provide an object-oriented framework for adding new types of models and packages. The object-oriented framework and underlying numerical solver make it possible to tightly couple any number of hyper-resolution models within coarser regional models. The hyper-resolution models can be used to evaluate local-scale groundwater issues that may be affected by regional-scale forcings. In MODFLOW 6, hyper-resolution meshes can be maintained as separate model datasets, similar to MODFLOW-LGR, which simplifies the development of a coarse regional model with imbedded hyper-resolution models from a coarse regional model. For example, the South Atlantic Coastal Plain regional water availability model was converted from a MODFLOW-2000 model to a MODFLOW 6 model. The horizontal discretization of the original model is approximately 3,218 m x 3,218 m. Hyper-resolution models of the Aiken and Sumter County water budget areas in South Carolina with a horizontal discretization of approximately 322 m x 322 m were developed and were tightly coupled to a modified version of the original coarse regional model that excluded these areas. Hydraulic property and aquifer geometry data from the coarse model were mapped to the hyper-resolution models. The discretization of the hyper-resolution models is fine enough to make detailed analyses of the effect that changes in groundwater withdrawals in the production aquifers have on the water table and surface-water/groundwater interactions. The approach used in this analysis could be applied to other regional water availability models that have been developed by the U.S. Geological Survey to evaluate local scale groundwater issues.

High Resolution Model Intercomparison Project (HighResMIP v1.0) for CMIP6

DOE PAGES

Haarsma, Reindert J.; Roberts, Malcolm J.; Vidale, Pier Luigi; ...

2016-11-22

Robust projections and predictions of climate variability and change, particularly at regional scales, rely on the driving processes being represented with fidelity in model simulations. The role of enhanced horizontal resolution in improved process representation in all components of the climate system is of growing interest, particularly as some recent simulations suggest both the possibility of significant changes in large-scale aspects of circulation as well as improvements in small-scale processes and extremes. However, such high-resolution global simulations at climate timescales, with resolutions of at least 50 km in the atmosphere and 0.25° in the ocean, have been performed at relativelymore » few research centres and generally without overall coordination, primarily due to their computational cost. Assessing the robustness of the response of simulated climate to model resolution requires a large multi-model ensemble using a coordinated set of experiments. The Coupled Model Intercomparison Project 6 (CMIP6) is the ideal framework within which to conduct such a study, due to the strong link to models being developed for the CMIP DECK experiments and other model intercomparison projects (MIPs). Increases in high-performance computing (HPC) resources, as well as the revised experimental design for CMIP6, now enable a detailed investigation of the impact of increased resolution up to synoptic weather scales on the simulated mean climate and its variability. The High Resolution Model Intercomparison Project (HighResMIP) presented in this paper applies, for the first time, a multi-model approach to the systematic investigation of the impact of horizontal resolution. A coordinated set of experiments has been designed to assess both a standard and an enhanced horizontal-resolution simulation in the atmosphere and ocean. The set of HighResMIP experiments is divided into three tiers consisting of atmosphere-only and coupled runs and spanning the period 1950–2050, with the possibility of extending to 2100, together with some additional targeted experiments. This paper describes the experimental set-up of HighResMIP, the analysis plan, the connection with the other CMIP6 endorsed MIPs, as well as the DECK and CMIP6 historical simulations. Lastly, HighResMIP thereby focuses on one of the CMIP6 broad questions, “what are the origins and consequences of systematic model biases?”, but we also discuss how it addresses the World Climate Research Program (WCRP) grand challenges.« less

Impact of numerical choices on water conservation in the E3SM Atmosphere Model version 1 (EAMv1)

NASA Astrophysics Data System (ADS)

Zhang, Kai; Rasch, Philip J.; Taylor, Mark A.; Wan, Hui; Leung, Ruby; Ma, Po-Lun; Golaz, Jean-Christophe; Wolfe, Jon; Lin, Wuyin; Singh, Balwinder; Burrows, Susannah; Yoon, Jin-Ho; Wang, Hailong; Qian, Yun; Tang, Qi; Caldwell, Peter; Xie, Shaocheng

2018-06-01

The conservation of total water is an important numerical feature for global Earth system models. Even small conservation problems in the water budget can lead to systematic errors in century-long simulations. This study quantifies and reduces various sources of water conservation error in the atmosphere component of the Energy Exascale Earth System Model. Several sources of water conservation error have been identified during the development of the version 1 (V1) model. The largest errors result from the numerical coupling between the resolved dynamics and the parameterized sub-grid physics. A hybrid coupling using different methods for fluid dynamics and tracer transport provides a reduction of water conservation error by a factor of 50 at 1° horizontal resolution as well as consistent improvements at other resolutions. The second largest error source is the use of an overly simplified relationship between the surface moisture flux and latent heat flux at the interface between the host model and the turbulence parameterization. This error can be prevented by applying the same (correct) relationship throughout the entire model. Two additional types of conservation error that result from correcting the surface moisture flux and clipping negative water concentrations can be avoided by using mass-conserving fixers. With all four error sources addressed, the water conservation error in the V1 model becomes negligible and insensitive to the horizontal resolution. The associated changes in the long-term statistics of the main atmospheric features are small. A sensitivity analysis is carried out to show that the magnitudes of the conservation errors in early V1 versions decrease strongly with temporal resolution but increase with horizontal resolution. The increased vertical resolution in V1 results in a very thin model layer at the Earth's surface, which amplifies the conservation error associated with the surface moisture flux correction. We note that for some of the identified error sources, the proposed fixers are remedies rather than solutions to the problems at their roots. Future improvements in time integration would be beneficial for V1.

High Resolution Simulations of Arctic Sea Ice, 1979-1993

DTIC Science & Technology

2003-01-01

William H. Lipscomb * PO[ARISSP To evaluate improvements in modelling Arctic sea ice, we compare results from two regional models at 1/120 horizontal...resolution. The first is a coupled ice-ocean model of the Arctic Ocean, consisting of an ocean model (adapted from the Parallel Ocean Program, Los...Alamos National Laboratory [LANL]) and the "old" sea ice model . The second model uses the same grid but consists of an improved "new" sea ice model (LANL

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/.

What model resolution is required in climatological downscaling over complex terrain?

NASA Astrophysics Data System (ADS)

El-Samra, Renalda; Bou-Zeid, Elie; El-Fadel, Mutasem

2018-05-01

This study presents results from the Weather Research and Forecasting (WRF) model applied for climatological downscaling simulations over highly complex terrain along the Eastern Mediterranean. We sequentially downscale general circulation model results, for a mild and wet year (2003) and a hot and dry year (2010), to three local horizontal resolutions of 9, 3 and 1 km. Simulated near-surface hydrometeorological variables are compared at different time scales against data from an observational network over the study area comprising rain gauges, anemometers, and thermometers. The overall performance of WRF at 1 and 3 km horizontal resolution was satisfactory, with significant improvement over the 9 km downscaling simulation. The total yearly precipitation from WRF's 1 km and 3 km domains exhibited < 10% bias with respect to observational data. The errors in minimum and maximum temperatures were reduced by the downscaling, along with a high-quality delineation of temperature variability and extremes for both the 1 and 3 km resolution runs. Wind speeds, on the other hand, are generally overestimated for all model resolutions, in comparison with observational data, particularly on the coast (up to 50%) compared to inland stations (up to 40%). The findings therefore indicate that a 3 km resolution is sufficient for the downscaling, especially that it would allow more years and scenarios to be investigated compared to the higher 1 km resolution at the same computational effort. In addition, the results provide a quantitative measure of the potential errors for various hydrometeorological variables.

Challenge toward the prediction of typhoon behaviour and down pour

NASA Astrophysics Data System (ADS)

Takahashi, K.; Onishi, R.; Baba, Y.; Kida, S.; Matsuda, K.; Goto, K.; Fuchigami, H.

2013-08-01

Mechanisms of interactions among different scale phenomena play important roles for forecasting of weather and climate. Multi-scale Simulator for the Geoenvironment (MSSG), which deals with multi-scale multi-physics phenomena, is a coupled non-hydrostatic atmosphere-ocean model designed to be run efficiently on the Earth Simulator. We present simulation results with the world-highest 1.9km horizontal resolution for the entire globe and regional heavy rain with 1km horizontal resolution and 5m horizontal/vertical resolution for urban area simulation. To gain high performance by exploiting the system capabilities, we propose novel performance evaluation metrics introduced in previous studies that incorporate the effects of the data caching mechanism between CPU and memory. With a useful code optimization guideline based on such metrics, we demonstrate that MSSG can achieve an excellent peak performance ratio of 32.2% on the Earth Simulator with the single-core performance found to be a key to a reduced time-to-solution.

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.

Effects of vertical shear in modelling horizontal oceanic dispersion

NASA Astrophysics Data System (ADS)

Lanotte, A. S.; Corrado, R.; Palatella, L.; Pizzigalli, C.; Schipa, I.; Santoleri, R.

2016-02-01

The effect of vertical shear on the horizontal dispersion properties of passive tracer particles on the continental shelf of the South Mediterranean is investigated by means of observation and model data. In situ current measurements reveal that vertical gradients of horizontal velocities in the upper mixing layer decorrelate quite fast ( ˜ 1 day), whereas an eddy-permitting ocean model, such as the Mediterranean Forecasting System, tends to overestimate such decorrelation time because of finite resolution effects. Horizontal dispersion, simulated by the Mediterranean sea Forecasting System, is mostly affected by: (1) unresolved scale motions, and mesoscale motions that are largely smoothed out at scales close to the grid spacing; (2) poorly resolved time variability in the profiles of the horizontal velocities in the upper layer. For the case study we have analysed, we show that a suitable use of deterministic kinematic parametrizations is helpful to implement realistic statistical features of tracer dispersion in two and three dimensions. The approach here suggested provides a functional tool to control the horizontal spreading of small organisms or substance concentrations, and is thus relevant for marine biology, pollutant dispersion as well as oil spill applications.

Some lessons and thoughts from development of an old-fashioned high-resolution atmospheric general circulation model

NASA Astrophysics Data System (ADS)

Ohfuchi, Wataru; Enomoto, Takeshi; Yoshioka, Mayumi K.; Takaya, Koutarou

2014-05-01

Some high-resolution simulations with a conventional atmospheric general circulation model (AGCM) were conducted right after the first Earth Simulator started operating in the spring of 2002. More simulations with various resolutions followed. The AGCM in this study, AFES (Agcm For the Earth Simulator), is a primitive equation spectral transform method model with a cumulus convection parameterization. In this presentation, some findings from comparisons between high and low-resolution simulations, and some future perspectives of old-fashioned AGCMs will be discussed. One obvious advantage of increasing resolution is capability of resolving the fine structures of topography and atmospheric flow. By increasing resolution from T39 (about 320 km horizontal grid interval) to T79 (160 km), to T159 (80 km) to T319 (40 km), topographic precipitation over Japan becomes increasingly realistic. This feature is necessary for climate and weather studies involving both global and local aspects. In order to resolve submesoscale (about 100 km horizontal scale) atmospheric circulation, about 10-km grid interval is necessary. Comparing T1279 (10 km) simulations with T319 ones, it is found that, for example, the intensity of heavy rain associated with Baiu front and the central pressure of typhoon become more realistic. These realistic submesoscale phenomena should have impact on larger-sale flow through dynamics and thermodynamics. An interesting finding by increasing horizontal resolution of a conventional AGCM is that some cumulus convection parameterizations, such as Arakawa-Schubert type scheme, gradually stop producing precipitation, while some others, such as Emanuel type, do not. With the former, the grid condensation increases with the model resolution to compensate. Which characteristics are more desirable is arguable but it is an important feature one has to consider when developing a high-resolution conventional AGCM. Many may think that conventional primitive equation spectral transform AGCMs, such as AFES, have no future. Developing globally homogeneous nonhydrostatic cloud resolving grid AGCMs is obviously a straightforward direction for the future. However these models will be very expensive for many users for a while, perhaps for the next some decades. On the other hand, old-fashioned AGCMs with a grid interval of 20-100 km will remain to be accurate and efficient tools for many users for many years to come. Also by coupling with a fine-resolution regional nonhydrostatic model, a conventional AGCM may overcome its limitation for use in climate and weather studies in the future.

Prediction of heavy rainfall over Chennai Metropolitan City, Tamil Nadu, India: Impact of microphysical parameterization schemes

NASA Astrophysics Data System (ADS)

Singh, K. S.; Bonthu, Subbareddy; Purvaja, R.; Robin, R. S.; Kannan, B. A. M.; Ramesh, R.

2018-04-01

This study attempts to investigate the real-time prediction of a heavy rainfall event over the Chennai Metropolitan City, Tamil Nadu, India that occurred on 01 December 2015 using Advanced Research Weather Research and Forecasting (WRF-ARW) model. The study evaluates the impact of six microphysical (Lin, WSM6, Goddard, Thompson, Morrison and WDM6) parameterization schemes of the model on prediction of heavy rainfall event. In addition, model sensitivity has also been evaluated with six Planetary Boundary Layer (PBL) and two Land Surface Model (LSM) schemes. Model forecast was carried out using nested domain and the impact of model horizontal grid resolutions were assessed at 9 km, 6 km and 3 km. Analysis of the synoptic features using National Center for Environmental Prediction Global Forecast System (NCEP-GFS) analysis data revealed strong upper-level divergence and high moisture content at lower level were favorable for the occurrence of heavy rainfall event over the northeast coast of Tamil Nadu. The study signified that forecasted rainfall was more sensitive to the microphysics and PBL schemes compared to the LSM schemes. The model provided better forecast of the heavy rainfall event using the logical combination of Goddard microphysics, YSU PBL and Noah LSM schemes, and it was mostly attributed to timely initiation and development of the convective system. The forecast with different horizontal resolutions using cumulus parameterization indicated that the rainfall prediction was not well represented at 9 km and 6 km. The forecast with 3 km horizontal resolution provided better prediction in terms of timely initiation and development of the event. The study highlights that forecast of heavy rainfall events using a high-resolution mesoscale model with suitable representations of physical parameterization schemes are useful for disaster management and planning to minimize the potential loss of life and property.

Eastern equatorial Pacific sea surface temperature annual cycle in the Kiel climate model: simulation benefits from enhancing atmospheric resolution

NASA Astrophysics Data System (ADS)

Wengel, C.; Latif, M.; Park, W.; Harlaß, J.; Bayr, T.

2018-05-01

A long-standing difficulty of climate models is to capture the annual cycle (AC) of eastern equatorial Pacific (EEP) sea surface temperature (SST). In this study, we first examine the EEP SST AC in a set of integrations of the coupled Kiel Climate Model, in which only atmosphere model resolution differs. When employing coarse horizontal and vertical atmospheric resolution, significant biases in the EEP SST AC are observed. These are reflected in an erroneous timing of the cold tongue's onset and termination as well as in an underestimation of the boreal spring warming amplitude. A large portion of these biases are linked to a wrong simulation of zonal surface winds, which can be traced back to precipitation biases on both sides of the equator and an erroneous low-level atmospheric circulation over land. Part of the SST biases also is related to shortwave radiation biases related to cloud cover biases. Both wind and cloud cover biases are inherent to the atmospheric component, as shown by companion uncoupled atmosphere model integrations forced by observed SSTs. Enhancing atmosphere model resolution, horizontal and vertical, markedly reduces zonal wind and cloud cover biases in coupled as well as uncoupled mode and generally improves simulation of the EEP SST AC. Enhanced atmospheric resolution reduces convection biases and improves simulation of surface winds over land. Analysis of a subset of models from the Coupled Model Intercomparison Project phase 5 (CMIP5) reveals that in these models, very similar mechanisms are at work in driving EEP SST AC biases.

Application of the LEPS technique for Quantitative Precipitation Forecasting (QPF) in Southern Italy: a preliminary study

NASA Astrophysics Data System (ADS)

Federico, S.; Avolio, E.; Bellecci, C.; Colacino, M.; Walko, R. L.

2006-03-01

This paper reports preliminary results for a Limited area model Ensemble Prediction System (LEPS), based on RAMS (Regional Atmospheric Modelling System), for eight case studies of moderate-intense precipitation over Calabria, the southernmost tip of the Italian peninsula. LEPS aims to transfer the benefits of a probabilistic forecast from global to regional scales in countries where local orographic forcing is a key factor to force convection. To accomplish this task and to limit computational time in an operational implementation of LEPS, we perform a cluster analysis of ECMWF-EPS runs. Starting from the 51 members that form the ECMWF-EPS we generate five clusters. For each cluster a representative member is selected and used to provide initial and dynamic boundary conditions to RAMS, whose integrations generate LEPS. RAMS runs have 12-km horizontal resolution. To analyze the impact of enhanced horizontal resolution on quantitative precipitation forecasts, LEPS forecasts are compared to a full Brute Force (BF) ensemble. This ensemble is based on RAMS, has 36 km horizontal resolution and is generated by 51 members, nested in each ECMWF-EPS member. LEPS and BF results are compared subjectively and by objective scores. Subjective analysis is based on precipitation and probability maps of case studies whereas objective analysis is made by deterministic and probabilistic scores. Scores and maps are calculated by comparing ensemble precipitation forecasts against reports from the Calabria regional raingauge network. Results show that LEPS provided better rainfall predictions than BF for all case studies selected. This strongly suggests the importance of the enhanced horizontal resolution, compared to ensemble population, for Calabria for these cases. To further explore the impact of local physiographic features on QPF (Quantitative Precipitation Forecasting), LEPS results are also compared with a 6-km horizontal resolution deterministic forecast. Due to local and mesoscale forcing, the high resolution forecast (Hi-Res) has better performance compared to the ensemble mean for rainfall thresholds larger than 10mm but it tends to overestimate precipitation for lower amounts. This yields larger false alarms that have a detrimental effect on objective scores for lower thresholds. To exploit the advantages of a probabilistic forecast compared to a deterministic one, the relation between the ECMWF-EPS 700 hPa geopotential height spread and LEPS performance is analyzed. Results are promising even if additional studies are required.

The influence of spatially and temporally high-resolution wind forcing on the power input to near-inertial waves in the ocean

NASA Astrophysics Data System (ADS)

Rimac, A.; Eden, C.; von Storch, J.

2012-12-01

Coexistence of stable stratification, the meridional overturning circulation and meso-scale eddies and their influence on the ocean's circulation still raise complex questions concerning the ocean energetics. Oceanic general circulation is mainly forced by the wind field and deep water tides. Its essential energetics are the conversion of kinetic energy of the winds and tides into oceanic potential and kinetic energy. Energy needed for the circulation is bound to internal wave fields. Direct internal wave generation by the wind at the sea surface is one of the sources of this energy. Previous studies using mixed-layer type of models and low frequency wind forcings (six-hourly and daily) left room for improvement. Using mixed-layer models it is not possible to assess the distribution of near-inertial energy into the deep ocean. Also, coarse temporal resolution of wind forcing strongly underestimates the near-inertial wave energy. To overcome this difficulty we use a high resolution ocean model with high frequency wind forcings. We establish the following model setup: We use the Max Planck Institute Ocean Model (MPIOM) on a tripolar grid with 45km horizontal resolution and 40 vertical levels. We run the model with wind forcings that vary in horizontal (250km versus 40km) and temporal resolution (six versus one-hourly). In our study we answer the following questions: How big is the wind kinetic energy input to the near-inertial waves? Is the kinetic energy of the near-inertial waves enhanced when high-frequency wind forcings are used? If so, by how much and why, due to higher level of temporal wind variability or due to better spatial representation of the near-inertial waves? How big is the total power of near-inertial waves generated by the wind at the surface of the ocean? We run the model for one year. Our model results show that the near-inertial waves are excited both using wind forcings of high and low horizontal and temporal resolution. Near-inertial energy is almost two times higher when we force the model with high frequency wind forcings. The influence on the energy mostly depends on the time difference between two forcing fields while the spatial difference has little influence.

The influence of model resolution on ozone in industrial volatile organic compound plumes.

PubMed

Henderson, Barron H; Jeffries, Harvey E; Kim, Byeong-Uk; Vizuete, William G

2010-09-01

Regions with concentrated petrochemical industrial activity (e.g., Houston or Baton Rouge) frequently experience large, localized releases of volatile organic compounds (VOCs). Aircraft measurements suggest these released VOCs create plumes with ozone (O3) production rates 2-5 times higher than typical urban conditions. Modeling studies found that simulating high O3 productions requires superfine (1-km) horizontal grid cell size. Compared with fine modeling (4-kmin), the superfine resolution increases the peak O3 concentration by as much as 46%. To understand this drastic O3 change, this study quantifies model processes for O3 and "odd oxygen" (Ox) in both resolutions. For the entire plume, the superfine resolution increases the maximum O3 concentration 3% but only decreases the maximum Ox concentration 0.2%. The two grid sizes produce approximately equal Ox mass but by different reaction pathways. Derived sensitivity to oxides of nitrogen (NOx) and VOC emissions suggests resolution-specific sensitivity to NOx and VOC emissions. Different sensitivity to emissions will result in different O3 responses to subsequently encountered emissions (within the city or downwind). Sensitivity of O3 to emission changes also results in different simulated O3 responses to the same control strategies. Sensitivity of O3 to NOx and VOC emission changes is attributed to finer resolved Eulerian grid and finer resolved NOx emissions. Urban NOx concentration gradients are often caused by roadway mobile sources that would not typically be addressed with Plume-in-Grid models. This study shows that grid cell size (an artifact of modeling) influences simulated control strategies and could bias regulatory decisions. Understanding the dynamics of VOC plume dependence on grid size is the first step toward providing more detailed guidance for resolution. These results underscore VOC and NOx resolution interdependencies best addressed by finer resolution. On the basis of these results, the authors suggest a need for quantitative metrics for horizontal grid resolution in future model guidance.

An intercomparison of GCM and RCM dynamical downscaling for characterizing the hydroclimatology of California and Nevada

NASA Astrophysics Data System (ADS)

Xu, Z.; Rhoades, A.; Johansen, H.; Ullrich, P. A.; Collins, W. D.

2017-12-01

Dynamical downscaling is widely used to properly characterize regional surface heterogeneities that shape the local hydroclimatology. However, the factors in dynamical downscaling, including the refinement of model horizontal resolution, large-scale forcing datasets and dynamical cores, have not been fully evaluated. Two cutting-edge global-to-regional downscaling methods are used to assess these, specifically the variable-resolution Community Earth System Model (VR-CESM) and the Weather Research & Forecasting (WRF) regional climate model, under different horizontal resolutions (28, 14, and 7 km). Two groups of WRF simulations are driven by either the NCEP reanalysis dataset (WRF_NCEP) or VR-CESM outputs (WRF_VRCESM) to evaluate the effects of the large-scale forcing datasets. The impacts of dynamical core are assessed by comparing the VR-CESM simulations to the coupled WRF_VRCESM simulations with the same physical parameterizations and similar grid domains. The simulated hydroclimatology (i.e., total precipitation, snow cover, snow water equivalent and surface temperature) are compared with the reference datasets. The large-scale forcing datasets are critical to the WRF simulations in more accurately simulating total precipitation, SWE and snow cover, but not surface temperature. Both the WRF and VR-CESM results highlight that no significant benefit is found in the simulated hydroclimatology by just increasing horizontal resolution refinement from 28 to 7 km. Simulated surface temperature is sensitive to the choice of dynamical core. WRF generally simulates higher temperatures than VR-CESM, alleviates the systematic cold bias of DJF temperatures over the California mountain region, but overestimates the JJA temperature in California's Central Valley.

Kinetic energy spectra, vertical resolution and dissipation in high-resolution atmospheric simulations.

NASA Astrophysics Data System (ADS)

Skamarock, W. C.

2017-12-01

We have performed week-long full-physics simulations with the MPAS global model at 15 km cell spacing using vertical mesh spacings of 800, 400, 200 and 100 meters in the mid-troposphere through the mid-stratosphere. We find that the horizontal kinetic energy spectra in the upper troposphere and stratosphere does not converge with increasing vertical resolution until we reach 200 meter level spacing. Examination of the solutions indicates that significant inertia-gravity waves are not vertically resolved at the lower vertical resolutions. Diagnostics from the simulations indicate that the primary kinetic energy dissipation results from the vertical mixing within the PBL parameterization and from the gravity-wave drag parameterization, with smaller but significant contributions from damping in the vertical transport scheme and from the horizontal filters in the dynamical core. Most of the kinetic energy dissipation in the free atmosphere occurs within breaking mid-latitude baroclinic waves. We will briefly review these results and their implications for atmospheric model configuration and for atmospheric dynamics, specifically that related to the dynamics associated with the mesoscale kinetic energy spectrum.

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.

Resolving vorticity-driven lateral fire spread using the WRF-Fire coupled atmosphere-fire numerical model

NASA Astrophysics Data System (ADS)

Simpson, C. C.; Sharples, J. J.; Evans, J. P.

2014-09-01

Vorticity-driven lateral fire spread (VLS) is a form of dynamic fire behaviour, during which a wildland fire spreads rapidly across a steep leeward slope in a direction approximately transverse to the background winds. VLS is often accompanied by a downwind extension of the active flaming region and intense pyro-convection. In this study, the WRF-Fire (WRF stands for Weather Research and Forecasting) coupled atmosphere-fire model is used to examine the sensitivity of resolving VLS to both the horizontal and vertical grid spacing, and the fire-to-atmosphere coupling from within the model framework. The atmospheric horizontal and vertical grid spacing are varied between 25 and 90 m, and the fire-to-atmosphere coupling is either enabled or disabled. At high spatial resolutions, the inclusion of fire-to-atmosphere coupling increases the upslope and lateral rate of spread by factors of up to 2.7 and 9.5, respectively. This increase in the upslope and lateral rate of spread diminishes at coarser spatial resolutions, and VLS is not modelled for a horizontal and vertical grid spacing of 90 m. The lateral fire spread is driven by fire whirls formed due to an interaction between the background winds and the vertical circulation generated at the flank of the fire front as part of the pyro-convective updraft. The laterally advancing fire fronts become the dominant contributors to the extreme pyro-convection. The results presented in this study demonstrate that both high spatial resolution and two-way atmosphere-fire coupling are required to model VLS with WRF-Fire.

Development of a High-Resolution Climate Model for Future Climate Change Projection on the Earth Simulator

NASA Astrophysics Data System (ADS)

Kanzawa, H.; Emori, S.; Nishimura, T.; Suzuki, T.; Inoue, T.; Hasumi, H.; Saito, F.; Abe-Ouchi, A.; Kimoto, M.; Sumi, A.

2002-12-01

The fastest supercomputer of the world, the Earth Simulator (total peak performance 40TFLOPS) has recently been available for climate researches in Yokohama, Japan. We are planning to conduct a series of future climate change projection experiments on the Earth Simulator with a high-resolution coupled ocean-atmosphere climate model. The main scientific aims for the experiments are to investigate 1) the change in global ocean circulation with an eddy-permitting ocean model, 2) the regional details of the climate change including Asian monsoon rainfall pattern, tropical cyclones and so on, and 3) the change in natural climate variability with a high-resolution model of the coupled ocean-atmosphere system. To meet these aims, an atmospheric GCM, CCSR/NIES AGCM, with T106(~1.1o) horizontal resolution and 56 vertical layers is to be coupled with an oceanic GCM, COCO, with ~ 0.28ox 0.19o horizontal resolution and 48 vertical layers. This coupled ocean-atmosphere climate model, named MIROC, also includes a land-surface model, a dynamic-thermodynamic seaice model, and a river routing model. The poles of the oceanic model grid system are rotated from the geographic poles so that they are placed in Greenland and Antarctic land masses to avoild the singularity of the grid system. Each of the atmospheric and the oceanic parts of the model is parallelized with the Message Passing Interface (MPI) technique. The coupling of the two is to be done with a Multi Program Multi Data (MPMD) fashion. A 100-model-year integration will be possible in one actual month with 720 vector processors (which is only 14% of the full resources of the Earth Simulator).

Simulations of Madden-Julian Oscillation in High Resolution Atmospheric General Circulation Model

NASA Astrophysics Data System (ADS)

Deng, Liping; Stenchikov, Georgiy; McCabe, Matthew; Bangalath, HamzaKunhu; Raj, Jerry; Osipov, Sergey

2014-05-01

The simulation of tropical signals, especially the Madden-Julian Oscillation (MJO), is one of the major deficiencies in current numerical models. The unrealistic features in the MJO simulations include the weak amplitude, more power at higher frequencies, displacement of the temporal and spatial distributions, eastward propagation speed being too fast, and a lack of coherent structure for the eastward propagation from the Indian Ocean to the Pacific (e.g., Slingo et al. 1996). While some improvement in simulating MJO variance and coherent eastward propagation has been attributed to model physics, model mean background state and air-sea interaction, studies have shown that the model resolution, especially for higher horizontal resolution, may play an important role in producing a more realistic simulation of MJO (e.g., Sperber et al. 2005). In this study, we employ unique high-resolution (25-km) simulations conducted using the Geophysical Fluid Dynamics Laboratory global High Resolution Atmospheric Model (HIRAM) to evaluate the MJO simulation against the European Center for Medium-range Weather Forecasts (ECMWF) Interim re-analysis (ERAI) dataset. We specifically focus on the ability of the model to represent the MJO related amplitude, spatial distribution, eastward propagation, and horizontal and vertical structures. Additionally, as the HIRAM output covers not only an historic period (1979-2012) but also future period (2012-2050), the impact of future climate change related to the MJO is illustrated. The possible changes in intensity and frequency of extreme weather and climate events (e.g., strong wind and heavy rainfall) in the western Pacific, the Indian Ocean and the Middle East North Africa (MENA) region are highlighted.

A global multilevel atmospheric model using a vector semi-Lagrangian finite-difference scheme. I - Adiabatic formulation

NASA Technical Reports Server (NTRS)

Bates, J. R.; Moorthi, S.; Higgins, R. W.

1993-01-01

An adiabatic global multilevel primitive equation model using a two time-level, semi-Lagrangian semi-implicit finite-difference integration scheme is presented. A Lorenz grid is used for vertical discretization and a C grid for the horizontal discretization. The momentum equation is discretized in vector form, thus avoiding problems near the poles. The 3D model equations are reduced by a linear transformation to a set of 2D elliptic equations, whose solution is found by means of an efficient direct solver. The model (with minimal physics) is integrated for 10 days starting from an initialized state derived from real data. A resolution of 16 levels in the vertical is used, with various horizontal resolutions. The model is found to be stable and efficient, and to give realistic output fields. Integrations with time steps of 10 min, 30 min, and 1 h are compared, and the differences are found to be acceptable.

Allocating emissions to 4 km and 1 km horizontal spatial resolutions and its impact on simulated NOx and O3 in Houston, TX

NASA Astrophysics Data System (ADS)

Pan, Shuai; Choi, Yunsoo; Roy, Anirban; Jeon, Wonbae

2017-09-01

A WRF-SMOKE-CMAQ air quality modeling system was used to investigate the impact of horizontal spatial resolution on simulated nitrogen oxides (NOx) and ozone (O3) in the Greater Houston area (a non-attainment area for O3). We employed an approach recommended by the United States Environmental Protection Agency to allocate county-based emissions to model grid cells in 1 km and 4 km horizontal grid resolutions. The CMAQ Integrated Process Rate analyses showed a substantial difference in emissions contributions between 1 and 4 km grids but similar NOx and O3 concentrations over urban and industrial locations. For example, the peak NOx emissions at an industrial and urban site differed by a factor of 20 for the 1 km and 8 for the 4 km grid, but simulated NOx concentrations changed only by a factor of 1.2 in both cases. Hence, due to the interplay of the atmospheric processes, we cannot expect a similar level of reduction of the gas-phase air pollutants as the reduction of emissions. Both simulations reproduced the variability of NASA P-3B aircraft measurements of NOy and O3 in the lower atmosphere (from 90 m to 4.5 km). Both simulations provided similar reasonable predictions at surface, while 1 km case depicted more detailed features of emissions and concentrations in heavily polluted areas, such as highways, airports, and industrial regions, which are useful in understanding the major causes of O3 pollution in such regions, and to quantify transport of O3 to populated communities in urban areas. The Integrated Reaction Rate analyses indicated a distinctive difference of chemistry processes between the model surface layer and upper layers, implying that correcting the meteorological conditions at the surface may not help to enhance the O3 predictions. The model-observation O3 bias in our studies (e.g., large over-prediction during the nighttime or along Gulf of Mexico coastline), were due to uncertainties in meteorology, chemistry or other processes. Horizontal grid resolution is unlikely the major contributor to these biases.

A Vertically Flow-Following, Icosahedral Grid Model for Medium-Range and Seasonal Prediction. Part 1: Model Description

NASA Technical Reports Server (NTRS)

Bleck, Rainer; Bao, Jian-Wen; Benjamin, Stanley G.; Brown, John M.; Fiorino, Michael; Henderson, Thomas B.; Lee, Jin-Luen; MacDonald, Alexander E.; Madden, Paul; Middlecoff, Jacques;

Wind-Stress Dust Lifting in a Mars Global Circulation Model: Representation across Resolutions

NASA Astrophysics Data System (ADS)

Chapman, R.; Lewis, S.; Balme, M. R.; Steele, L.

2017-12-01

The formation of Martian dust storms is believed to be driven by dust lifting by near-surface wind stress (NSWS). Accurately representing this dust lifting within Mars Global Circulation Models (MGCMs) is important in order to gain a full understanding of the Martian dust storm cycle. Parameterisations of dust lifting by NSWS exist within several MGCMs; implementations differ but they all follow a similar design, so progress within one model is relevant to the entire field. Few studies have explored in detail how the results of these parameterisations can be affected by changing the horizontal resolution of the model. An accurate parameterisation of dust lifting by NSWS will lift a representative dust mass, reproducing characteristic dust optical depths in the atmosphere. The geographical distribution of the dust lifting by NSWS will also change throughout the year, affecting patterns of dust storm formation and development. Currently, suitable values for dust lifting parameters must be identified at every new model resolution. Resolutions of 5° latitude x 5° longitude are often used to model the Martian climate, as thermal tides and long-term weather patterns can be well represented at this resolution. However, smaller scale phenomena (such as near-surface winds driven by local topography) cannot be accurately depicted at this resolution. We use the LMD-UK MGCM to complete multi-year simulations across multiple model resolutions. Our experiments range from `low' resolution 5° lat x 5° lon to `high' resolution 1° lat x 1° lon. In experiments with fixed, constant lifting parameters, we find that higher resolution simulations lift more dust, but that this trend is asymptotic. At low resolutions, dust lifting increases proportionately with the increase in number of horizontal gridboxes. However, at high resolutions, doubling the number of gridboxes results only in a 30% increase in the total dust mass lifted. Geographical and temporal distributions of dust lifting are investigated, as well as the total dust lifted, in order to assess the optimum parameters for each resolution, and to develop a calibration scheme for this dust lifting across model resolutions. The scheme is verified through comparison with spacecraft observations of dust optical depths and dust storm locations.

Effects of different regional climate model resolution and forcing scales on projected hydrologic changes

NASA Astrophysics Data System (ADS)

Mendoza, Pablo A.; Mizukami, Naoki; Ikeda, Kyoko; Clark, Martyn P.; Gutmann, Ethan D.; Arnold, Jeffrey R.; Brekke, Levi D.; Rajagopalan, Balaji

2016-10-01

We examine the effects of regional climate model (RCM) horizontal resolution and forcing scaling (i.e., spatial aggregation of meteorological datasets) on the portrayal of climate change impacts. Specifically, we assess how the above decisions affect: (i) historical simulation of signature measures of hydrologic behavior, and (ii) projected changes in terms of annual water balance and hydrologic signature measures. To this end, we conduct our study in three catchments located in the headwaters of the Colorado River basin. Meteorological forcings for current and a future climate projection are obtained at three spatial resolutions (4-, 12- and 36-km) from dynamical downscaling with the Weather Research and Forecasting (WRF) regional climate model, and hydrologic changes are computed using four different hydrologic model structures. These projected changes are compared to those obtained from running hydrologic simulations with current and future 4-km WRF climate outputs re-scaled to 12- and 36-km. The results show that the horizontal resolution of WRF simulations heavily affects basin-averaged precipitation amounts, propagating into large differences in simulated signature measures across model structures. The implications of re-scaled forcing datasets on historical performance were primarily observed on simulated runoff seasonality. We also found that the effects of WRF grid resolution on projected changes in mean annual runoff and evapotranspiration may be larger than the effects of hydrologic model choice, which surpasses the effects from re-scaled forcings. Scaling effects on projected variations in hydrologic signature measures were found to be generally smaller than those coming from WRF resolution; however, forcing aggregation in many cases reversed the direction of projected changes in hydrologic behavior.

: “Developing Regional Modeling Techniques Applicable for Simulating Future Climate Conditions in the Carolinas”

EPA Science Inventory

Global climate models (GCMs) are currently used to obtain information about future changes in the large-scale climate. However, such simulations are typically done at coarse spatial resolutions, with model grid boxes on the order of 100 km on a horizontal side. Therefore, techniq...

Monsoon Variability in the Arabian Sea from Enhanced and Standard Horizontal Resolution Coupled Climate Models.

NASA Astrophysics Data System (ADS)

McClean, J.; Veneziani, C.; Maltrud, M. E.; Taylor, M.; Bader, D. C.; Branstetter, M. L.; Evans, K. J.; Mahajan, S.

2016-02-01

The circulation of the upper ocean in the Arabian Sea switches direction seasonally due to the change in direction of the prevailing winds associated with the Indian Monsoon. Predictability of the monsoon circulation, however, is uncertain due to incomplete understanding of the physical processes operating on the monsoon and other time scales, particularly interannual and intraseasonal. We use the Community Earth System Model (CESM) with enhanced horizontal resolution in each of its components relative to standard coupled climate model resolution, to better understand these time scale interactions. A standard resolution CESM counterpart is used to assess how horizontal resolution impacts the depiction of these processes. In the enhanced resolution case, 0.25° Community Atmosphere Model 5 (CAM5) is coupled to, among other components, the tripolar nominal 0.1° Parallel Ocean Program 2 (POP2). The fine resolution CESM simulation was run for 85 years; constant 1850 preindustrial forcing was used throughout the run, allowing us to isolate internal variability of the coupled system. Model parameters were adjusted ("tuned") to produce an acceptably small top of the atmosphere radiation imbalance. The reversal of the Somali Current (SC), the western boundary current off northeast Africa, has typically been associated with that of the monsoon. The SC reverses from southwestward in boreal winter to northeastward in summer; coastal upwelling is induced by the summer monsoonal winds. Recently it has been shown from new observations that the SC starts to reverse prior to the monsoon switch. Westward propagating Rossby waves have been implicated as responsible for the early SC reversal. We will discuss the sequencing of remote and local forcing on the timing of the spring inter-monsoonal switch in the direction of the SC and the appearance of the Great Whirl off the Oman Coast. Particularly, we consider how the Indian Ocean Dipole (IOD) acts to modify the seasonal strength and variability of the western boundary current system including upwelling. We look for a connection between interannual upwelling variability and that of rainfall off the west coast of India. As well, we examine changes due to the IOD in the upper ocean temperature and salinity structure along the Rossby wave propagation route in the Arabian Sea.

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.

Simulations of the transport and deposition of 137Cs over Europe after the Chernobyl NPP accident: influence of varying emission-altitude and model horizontal and vertical resolution

NASA Astrophysics Data System (ADS)

Evangeliou, N.; Balkanski, Y.; Cozic, A.; Møller, A. P.

2013-03-01

The coupled model LMDzORINCA has been used to simulate the transport, wet and dry deposition of the radioactive tracer 137Cs after accidental releases. For that reason, two horizontal resolutions were deployed and used in the model, a regular grid of 2.5°×1.25°, and the same grid stretched over Europe to reach a resolution of 0.45°×0.51°. The vertical dimension is represented with two different resolutions, 19 and 39 levels, respectively, extending up to mesopause. Four different simulations are presented in this work; the first uses the regular grid over 19 vertical levels assuming that the emissions took place at the surface (RG19L(S)), the second also uses the regular grid over 19 vertical levels but realistic source injection heights (RG19L); in the third resolution the grid is regular and the vertical resolution 39 vertical levels (RG39L) and finally, it is extended to the stretched grid with 19 vertical levels (Z19L). The best choice for the model validation was the Chernobyl accident which occurred in Ukraine (ex-USSR) on 26 May 1986. This accident has been widely studied since 1986, and a large database has been created containing measurements of atmospheric activity concentration and total cumulative deposition for 137Cs from most of the European countries. According to the results, the performance of the model to predict the transport and deposition of the radioactive tracer was efficient and accurate presenting low biases in activity concentrations and deposition inventories, despite the large uncertainties on the intensity of the source released. However, the best agreement with observations was obtained using the highest horizontal resolution of the model (Z19L run). The model managed to predict the radioactive contamination in most of the European regions (similar to Atlas), and also the arrival times of the radioactive fallout. As regards to the vertical resolution, the largest biases were obtained for the 39 layers run due to the increase of the levels in conjunction with the uncertainty of the source term. Moreover, the ecological half-life of 137Cs in the atmosphere after the accident ranged between 6 and 9 days, which is in good accordance to what previously reported and in the same range with the recent accident in Japan. The high response of LMDzORINCA model for 137Cs reinforces the importance of atmospheric modeling in emergency cases to gather information for protecting the population from the adverse effects of radiation.

Simulations of the transport and deposition of 137Cs over Europe after the Chernobyl Nuclear Power Plant accident: influence of varying emission-altitude and model horizontal and vertical resolution

NASA Astrophysics Data System (ADS)

Evangeliou, N.; Balkanski, Y.; Cozic, A.; Møller, A. P.

2013-07-01

The coupled model LMDZORINCA has been used to simulate the transport, wet and dry deposition of the radioactive tracer 137Cs after accidental releases. For that reason, two horizontal resolutions were deployed and used in the model, a regular grid of 2.5° × 1.27°, and the same grid stretched over Europe to reach a resolution of 0.66° × 0.51°. The vertical dimension is represented with two different resolutions, 19 and 39 levels respectively, extending up to the mesopause. Four different simulations are presented in this work; the first uses the regular grid over 19 vertical levels assuming that the emissions took place at the surface (RG19L(S)), the second also uses the regular grid over 19 vertical levels but realistic source injection heights (RG19L); in the third resolution the grid is regular and the vertical resolution 39 levels (RG39L) and finally, it is extended to the stretched grid with 19 vertical levels (Z19L). The model is validated with the Chernobyl accident which occurred in Ukraine (ex-USSR) on 26 May 1986 using the emission inventory from Brandt et al. (2002). This accident has been widely studied since 1986, and a large database has been created containing measurements of atmospheric activity concentration and total cumulative deposition for 137Cs from most of the European countries. According to the results, the performance of the model to predict the transport and deposition of the radioactive tracer was efficient and accurate presenting low biases in activity concentrations and deposition inventories, despite the large uncertainties on the intensity of the source released. The best agreement with observations was obtained using the highest horizontal resolution of the model (Z19L run). The model managed to predict the radioactive contamination in most of the European regions (similar to De Cort et al., 1998), and also the arrival times of the radioactive fallout. As regards to the vertical resolution, the largest biases were obtained for the 39 layers run due to the increase of the levels in conjunction with the uncertainty of the source term. Moreover, the ecological half-life of 137Cs in the atmosphere after the accident ranged between 6 and 9 days, which is in good accordance to what previously reported and in the same range with the recent accident in Japan. The high response of LMDZORINCA model for 137Cs reinforces the importance of atmospheric modelling in emergency cases to gather information for protecting the population from the adverse effects of radiation.

Mesoscale atmospheric modelling technology as a tool for the long-term meteorological dataset development

NASA Astrophysics Data System (ADS)

Platonov, Vladimir; Kislov, Alexander; Rivin, Gdaly; Varentsov, Mikhail; Rozinkina, Inna; Nikitin, Mikhail; Chumakov, Mikhail

2017-04-01

The detailed hydrodynamic modelling of meteorological parameters during the last 30 years (1985 - 2014) was performed for the Okhotsk Sea and the Sakhalin island regions. The regional non-hydrostatic atmospheric model COSMO-CLM used for this long-term simulation with 13.2, 6.6 and 2.2 km horizontal resolutions. The main objective of creation this dataset was the outlook of the investigation of statistical characteristics and the physical mechanisms of extreme weather events (primarily, wind speed extremes) on the small spatio-temporal scales. COSMO-CLM is the climate version of the well-known mesoscale COSMO model, including some modifications and extensions adapting to the long-term numerical experiments. The downscaling technique was realized and developed for the long-term simulations with three consequent nesting domains. ERA-Interim reanalysis ( 0.75 degrees resolution) used as global forcing data for the starting domain ( 13.2 km horizontal resolution), then these simulation data used as initial and boundary conditions for the next model runs over the domain with 6.6 km resolution, and similarly, for the next step to 2.2 km domain. Besides, the COSMO-CLM model configuration for 13.2 km run included the spectral nudging technique, i.e. an additional assimilation of reanalysis data not only at boundaries, but also inside the whole domain. Practically, this computational scheme realized on the SGI Altix 4700 supercomputer system in the Main Computer Center of Roshydromet and used 2,400 hours of CPU time total. According to modelling results, the verification of the obtained dataset was performed on the observation data. Estimations showed the mean error -0.5 0C, up to 2 - 3 0C RMSE in temperature, and overestimation in wind speed (RMSE is up to 2 m/s). Overall, analysis showed that the used downscaling technique with applying the COSMO-CLM model reproduced the meteorological conditions, spatial distribution, seasonal and synoptic variability of temperature and wind speed for the study area adequately. The dependences between reproduction quality of mesoscale atmospheric circulation features and the horizontal resolution of the model were revealed. In particular, it is shown that the use of 6 km resolution does not give any significant improvement comparing to 13 km resolution, whereas 2.2 km resolution provides an appreciable quality enhancement. Detailed synoptic analysis of extreme wind speed situations identified the main types of favorable to their genesis, associated with developing of cyclones over the Japan Islands or the Primorsky Kray of Russia, and penetration of intensified cyclones from Pacific Ocean through the Kamchatka peninsula, Kuril or Japan Islands. The obtained dataset will continue to be used for a full and comprehensive analysis of the reproduction quality of hydrometeorological fields, their statistical estimates, climatological trends and many other objectives.

Rainfall From Resolved Rather Than Parameterized Processes Better Represents the Present-Day and Climate Change Response of Moderate Rates in the Community Atmosphere Model

DOE PAGES

Kooperman, Gabriel J.; Pritchard, Michael S.; O'Brien, Travis A.; ...

2018-04-01

Deficiencies in the parameterizations of convection used in global climate models often lead to a distorted representation of the simulated rainfall intensity distribution (i.e., too much rainfall from weak rain rates). While encouraging improvements in high percentile rainfall intensity have been found as the horizontal resolution of the Community Atmosphere Model is increased to ~25 km, we demonstrate no corresponding improvement in the moderate rain rates that generate the majority of accumulated rainfall. Using a statistical framework designed to emphasize links between precipitation intensity and accumulated rainfall beyond just the frequency distribution, we show that CAM cannot realistically simulate moderatemore » rain rates, and cannot capture their intensification with climate change, even as resolution is increased. However, by separating the parameterized convective and large-scale resolved contributions to total rainfall, we find that the intensity, geographic pattern, and climate change response of CAM's large-scale rain rates are more consistent with observations (TRMM 3B42), superparameterization, and theoretical expectations, despite issues with parameterized convection. Increasing CAM's horizontal resolution does improve the representation of total rainfall intensity, but not due to changes in the intensity of large-scale rain rates, which are surprisingly insensitive to horizontal resolution. Rather, improvements occur through an increase in the relative contribution of the large-scale component to the total amount of accumulated rainfall. Analysis of sensitivities to convective timescale and entrainment rate confirm the importance of these parameters in the possible development of scale-aware parameterizations, but also reveal unrecognized trade-offs from the entanglement of precipitation frequency and total amount.« less

Rainfall From Resolved Rather Than Parameterized Processes Better Represents the Present-Day and Climate Change Response of Moderate Rates in the Community Atmosphere Model

NASA Astrophysics Data System (ADS)

Kooperman, Gabriel J.; Pritchard, Michael S.; O'Brien, Travis A.; Timmermans, Ben W.

2018-04-01

Deficiencies in the parameterizations of convection used in global climate models often lead to a distorted representation of the simulated rainfall intensity distribution (i.e., too much rainfall from weak rain rates). While encouraging improvements in high percentile rainfall intensity have been found as the horizontal resolution of the Community Atmosphere Model is increased to ˜25 km, we demonstrate no corresponding improvement in the moderate rain rates that generate the majority of accumulated rainfall. Using a statistical framework designed to emphasize links between precipitation intensity and accumulated rainfall beyond just the frequency distribution, we show that CAM cannot realistically simulate moderate rain rates, and cannot capture their intensification with climate change, even as resolution is increased. However, by separating the parameterized convective and large-scale resolved contributions to total rainfall, we find that the intensity, geographic pattern, and climate change response of CAM's large-scale rain rates are more consistent with observations (TRMM 3B42), superparameterization, and theoretical expectations, despite issues with parameterized convection. Increasing CAM's horizontal resolution does improve the representation of total rainfall intensity, but not due to changes in the intensity of large-scale rain rates, which are surprisingly insensitive to horizontal resolution. Rather, improvements occur through an increase in the relative contribution of the large-scale component to the total amount of accumulated rainfall. Analysis of sensitivities to convective timescale and entrainment rate confirm the importance of these parameters in the possible development of scale-aware parameterizations, but also reveal unrecognized trade-offs from the entanglement of precipitation frequency and total amount.

Rainfall From Resolved Rather Than Parameterized Processes Better Represents the Present‐Day and Climate Change Response of Moderate Rates in the Community Atmosphere Model

PubMed Central

Pritchard, Michael S.; O'Brien, Travis A.; Timmermans, Ben W.

2018-01-01

Abstract Deficiencies in the parameterizations of convection used in global climate models often lead to a distorted representation of the simulated rainfall intensity distribution (i.e., too much rainfall from weak rain rates). While encouraging improvements in high percentile rainfall intensity have been found as the horizontal resolution of the Community Atmosphere Model is increased to ∼25 km, we demonstrate no corresponding improvement in the moderate rain rates that generate the majority of accumulated rainfall. Using a statistical framework designed to emphasize links between precipitation intensity and accumulated rainfall beyond just the frequency distribution, we show that CAM cannot realistically simulate moderate rain rates, and cannot capture their intensification with climate change, even as resolution is increased. However, by separating the parameterized convective and large‐scale resolved contributions to total rainfall, we find that the intensity, geographic pattern, and climate change response of CAM's large‐scale rain rates are more consistent with observations (TRMM 3B42), superparameterization, and theoretical expectations, despite issues with parameterized convection. Increasing CAM's horizontal resolution does improve the representation of total rainfall intensity, but not due to changes in the intensity of large‐scale rain rates, which are surprisingly insensitive to horizontal resolution. Rather, improvements occur through an increase in the relative contribution of the large‐scale component to the total amount of accumulated rainfall. Analysis of sensitivities to convective timescale and entrainment rate confirm the importance of these parameters in the possible development of scale‐aware parameterizations, but also reveal unrecognized trade‐offs from the entanglement of precipitation frequency and total amount. PMID:29861837

Rainfall From Resolved Rather Than Parameterized Processes Better Represents the Present-Day and Climate Change Response of Moderate Rates in the Community Atmosphere Model

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

Kooperman, Gabriel J.; Pritchard, Michael S.; O'Brien, Travis A.

Deficiencies in the parameterizations of convection used in global climate models often lead to a distorted representation of the simulated rainfall intensity distribution (i.e., too much rainfall from weak rain rates). While encouraging improvements in high percentile rainfall intensity have been found as the horizontal resolution of the Community Atmosphere Model is increased to ~25 km, we demonstrate no corresponding improvement in the moderate rain rates that generate the majority of accumulated rainfall. Using a statistical framework designed to emphasize links between precipitation intensity and accumulated rainfall beyond just the frequency distribution, we show that CAM cannot realistically simulate moderatemore » rain rates, and cannot capture their intensification with climate change, even as resolution is increased. However, by separating the parameterized convective and large-scale resolved contributions to total rainfall, we find that the intensity, geographic pattern, and climate change response of CAM's large-scale rain rates are more consistent with observations (TRMM 3B42), superparameterization, and theoretical expectations, despite issues with parameterized convection. Increasing CAM's horizontal resolution does improve the representation of total rainfall intensity, but not due to changes in the intensity of large-scale rain rates, which are surprisingly insensitive to horizontal resolution. Rather, improvements occur through an increase in the relative contribution of the large-scale component to the total amount of accumulated rainfall. Analysis of sensitivities to convective timescale and entrainment rate confirm the importance of these parameters in the possible development of scale-aware parameterizations, but also reveal unrecognized trade-offs from the entanglement of precipitation frequency and total amount.« less

LFNet: A Novel Bidirectional Recurrent Convolutional Neural Network for Light-Field Image Super-Resolution.

PubMed

Wang, Yunlong; Liu, Fei; Zhang, Kunbo; Hou, Guangqi; Sun, Zhenan; Tan, Tieniu

2018-09-01

The low spatial resolution of light-field image poses significant difficulties in exploiting its advantage. To mitigate the dependency of accurate depth or disparity information as priors for light-field image super-resolution, we propose an implicitly multi-scale fusion scheme to accumulate contextual information from multiple scales for super-resolution reconstruction. The implicitly multi-scale fusion scheme is then incorporated into bidirectional recurrent convolutional neural network, which aims to iteratively model spatial relations between horizontally or vertically adjacent sub-aperture images of light-field data. Within the network, the recurrent convolutions are modified to be more effective and flexible in modeling the spatial correlations between neighboring views. A horizontal sub-network and a vertical sub-network of the same network structure are ensembled for final outputs via stacked generalization. Experimental results on synthetic and real-world data sets demonstrate that the proposed method outperforms other state-of-the-art methods by a large margin in peak signal-to-noise ratio and gray-scale structural similarity indexes, which also achieves superior quality for human visual systems. Furthermore, the proposed method can enhance the performance of light field applications such as depth estimation.

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

A Coupled Surface Nudging Scheme for use in Retrospective ...

EPA Pesticide Factsheets

A surface analysis nudging scheme coupling atmospheric and land surface thermodynamic parameters has been implemented into WRF v3.8 (latest version) for use with retrospective weather and climate simulations, as well as for applications in air quality, hydrology, and ecosystem modeling. This scheme is known as the flux-adjusting surface data assimilation system (FASDAS) developed by Alapaty et al. (2008). This scheme provides continuous adjustments for soil moisture and temperature (via indirect nudging) and for surface air temperature and water vapor mixing ratio (via direct nudging). The simultaneous application of indirect and direct nudging maintains greater consistency between the soil temperature–moisture and the atmospheric surface layer mass-field variables. The new method, FASDAS, consistently improved the accuracy of the model simulations at weather prediction scales for different horizontal grid resolutions, as well as for high resolution regional climate predictions. This new capability has been released in WRF Version 3.8 as option grid_sfdda = 2. This new capability increased the accuracy of atmospheric inputs for use air quality, hydrology, and ecosystem modeling research to improve the accuracy of respective end-point research outcome. IMPACT: A new method, FASDAS, was implemented into the WRF model to consistently improve the accuracy of the model simulations at weather prediction scales for different horizontal grid resolutions, as wel

Effects of Grid Resolution on Modeled Air Pollutant Concentrations Due to Emissions from Large Point Sources: Case Study during KORUS-AQ 2016 Campaign

NASA Astrophysics Data System (ADS)

Ju, H.; Bae, C.; Kim, B. U.; Kim, H. C.; Kim, S.

2017-12-01

Large point sources in the Chungnam area received a nation-wide attention in South Korea because the area is located southwest of the Seoul Metropolitan Area whose population is over 22 million and the summertime prevalent winds in the area is northeastward. Therefore, emissions from the large point sources in the Chungnam area were one of the major observation targets during the KORUS-AQ 2016 including aircraft measurements. In general, horizontal grid resolutions of eulerian photochemical models have profound effects on estimated air pollutant concentrations. It is due to the formulation of grid models; that is, emissions in a grid cell will be assumed to be mixed well under planetary boundary layers regardless of grid cell sizes. In this study, we performed series of simulations with the Comprehensive Air Quality Model with eXetension (CAMx). For 9-km and 3-km simulations, we used meteorological fields obtained from the Weather Research and Forecast model while utilizing the "Flexi-nesting" option in the CAMx for the 1-km simulation. In "Flexi-nesting" mode, CAMx interpolates or assigns model inputs from the immediate parent grid. We compared modeled concentrations with ground observation data as well as aircraft measurements to quantify variations of model bias and error depending on horizontal grid resolutions.

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.