A standard test case suite for two-dimensional linear transport on the sphere

It is the purpose of this paper to propose a standard test case suite for two-dimensional transport schemes on the sphere intended to be used for model development and facilitating scheme intercomparison. The test cases are designed to assess important aspects of accuracy in geophysical fluid dynamics such as numerical order of convergence, "minimal" resolution, the ability of the transport scheme to preserve filaments, transport "rough" distributions, and to preserve pre-existing functional relations between species/tracers under challenging flow conditions. <br><br> The experiments are designed to be easy to set up. They are specified in terms of two analytical wind fields (one non-divergent and one divergent) and four analytical initial conditions (varying from smooth to discontinuous). Both conventional error norms as well as novel mixing and filament preservation diagnostics are used that are easy to implement. The experiments pose different challenges for the range of transport approaches from Lagrangian to Eulerian. The mixing and filament preservation diagnostics do not require an analytical/reference solution, which is in contrast to standard error norms where a "true" solution is needed. Results using the CSLAM (Conservative Semi-Lagrangian Multi-tracer) scheme on the cubed-sphere are presented for reference and illustrative purposes

Thunderstorms and upper troposphere chemistry during the early stages of the 2006 North American Monsoon

To study the meteorology and chemistry that is associated with the early stages of the North American Monsoon, the Weather Research and Forecasting model coupled with Chemistry (WRF-Chem) is applied for the first time at high resolution (4 km grid spacing, allowing for explicit representation of convection) over a large region (continental US and northern Mexico) for a multi-week (15 July to 7 August 2006) integration. Evaluation of model results shows that WRF-Chem reasonably represents the large-scale meteorology and strong convective storms, but tends to overestimate weak convection. In the upper troposphere, the WRF-Chem model predicts ozone (O<sub>3</sub>) and carbon monoxide (CO) to within 10–20% of aircraft and sonde measurements. Comparison of UT O<sub>3</sub> and CO frequency distributions between WRF-Chem and satellite data indicates that WRF-Chem is lofting CO too frequently from the boundary layer (BL). This excessive lofting should also cause biases in the WRF-Chem ozone frequency distribution; however it agrees well with satellite data suggesting that either the chemical production of O<sub>3</sub> in the model is overpredicted or there is too much stratosphere to troposphere transport in the model. Analysis of different geographic regions (West Coast, Rocky Mountains, Central Plains, Midwest, and Gulf Coast) reveals that much of the convective transport occurs in the Rocky Mountains, while much of the UT ozone chemical production occurs over the Gulf Coast and Midwest regions where both CO and volatile organic compounds (VOCs) are abundant in the upper troposphere and promote the production of peroxy radicals. In all regions most of the ozone chemical production occurs within 24 h of the air being lofted from the boundary layer. In addition, analysis of the anticyclone and adjacent air indicates that ozone mixing ratios within the anticyclone region associated with the North American Monsoon and just outside the anticyclone are similar. Increases of O<sub>3</sub> within the anticyclone are strongly coincident with entrainment of stratospheric air into the anticyclone, but also are from in situ O<sub>3</sub> chemical production. In situ O<sub>3</sub> production is up to 17% greater within the anticyclone than just outside the anticyclone when the anticyclone is over the southern US indicating that the enhancement of O<sub>3</sub> is most pronounced over regions with abundant VOCs

A meteorological overview of the MILAGRO field campaigns

International audienceWe describe the large-scale meteorological conditions that affected atmospheric chemistry over Mexico during March 2006 when several field campaigns were conducted in the region. In-situ and remote-sensing instrumentation was deployed to obtain measurements of wind, temperature, and humidity profiles in the boundary layer and free atmosphere at four primary sampling sites in central Mexico. Several models were run operationally during the field campaign to provide forecasts of the local, regional, and synoptic meteorology as well as the predicted location of the Mexico City pollutant plume for aircraft flight planning purposes. Field campaign measurements and large-scale analyses are used to define three regimes that characterize the overall meteorological conditions: the first regime prior to 14 March, the second regime between 14 and 23 March, and the third regime after 23 March. Mostly sunny and dry conditions with periods of cirrus and marine stratus along the coast occurred during the first regime. The beginning of the second regime was characterized by a sharp increase in humidity over the central plateau and the development of late afternoon convection associated with the passage of a weak cold surge on 14 March. Over the next several days, the atmosphere over the central plateau became drier so that deep convection gradually diminished. The third regime began with the passage of a strong cold surge that lead to humidity, afternoon convection, and precipitation over the central plateau that was higher than during the second regime. The frequency and intensity of fires, as determined by satellite measurements, also diminished significantly after the third cold surge. The synoptic-scale flow patterns that govern the transport of pollutants in the region are described and compared to previous March periods to put the transport into a climatological context. The complex terrain surrounding Mexico City produces local and regional circulations that govern short-range transport; however, the mean synoptic conditions modulate the thermally-driven circulations and on several days the near-surface flow is coupled to the ambient winds aloft

Wake response to an ocean-feedback mechanism: Madeira Island case study

This discussion focused on the numerical study of a wake episode. The Weather Research and Forecasting model was used in a downscale mode. The current literature focuses the discussion on the adiabatic dynamics of atmospheric wakes. Changes in mountain height and consequently on its relation to the atmospheric inversion layer should explain the shift in wake regimes: from a 'strong-wake' to a 'weak-wake' scenario. Nevertheless, changes in SST variability can also induce similar regime shifts. Increase in evaporation, contributes to increase convection and thus to an uplift of the stratified atmospheric layer, above the critical height, with subsequent internal gravity wave activity.Comment: Under review proces

ELPIS-JP: a dataset of local-scale daily climate change scenarios for Japan

We developed a dataset of local-scale daily climate change scenarios for Japan (called ELPIS-JP) using the stochastic weather generators (WGs) LARS-WG and, in part, WXGEN. The ELPIS-JP dataset is based on the observed (or estimated) daily weather data for seven climatic variables (daily mean, maximum and minimum temperatures; precipitation; solar radiation; relative humidity; and wind speed) at 938 sites in Japan and climate projections from the multi-model ensemble of global climate models (GCMs) used in the coupled model intercomparison project (CMIP3) and multi-model ensemble of regional climate models form the Japanese downscaling project (called S-5-3). The capability of the WGs to reproduce the statistical features of the observed data for the period 1981–2000 is assessed using several statistical tests and quantile–quantile plots. Overall performance of the WGs was good. The ELPIS-JP dataset consists of two types of daily data: (i) the transient scenarios throughout the twenty-first century using projections from 10 CMIP3 GCMs under three emission scenarios (A1B, A2 and B1) and (ii) the time-slice scenarios for the period 2081–2100 using projections from three S-5-3 regional climate models. The ELPIS-JP dataset is designed to be used in conjunction with process-based impact models (e.g. crop models) for assessment, not only the impacts of mean climate change but also the impacts of changes in climate variability, wet/dry spells and extreme events, as well as the uncertainty of future impacts associated with climate models and emission scenarios. The ELPIS-JP offers an excellent platform for probabilistic assessment of climate change impacts and potential adaptation at a local scale in Japan

Regional versus remote atmosphere‐ocean drivers of the rapid projected intensification of the East Australian Current

Like many western boundary currents, the East Australian Current extension is projected to get stronger and warmer in the future. The CMIP5 multi‐model mean (MMM) projection suggests up to 5°C of warming under an RCP85 scenario by 2100. Previous studies employed Sverdrup balance to associate a trend in basin wide zonally integrated wind stress curl (resulting from the multi‐decadal poleward intensification in the westerly winds over the Southern Ocean) with enhanced transport in the EAC extension. Possible regional drivers are yet to be considered. Here, we introduce the NEMO‐OASIS‐WRF coupled regional climate model as a framework to improve our understanding of CMIP5 projections. We analyse a hierarchy of simulations in which the regional atmosphere and ocean circulations are allowed to freely evolve subject to boundary conditions that represent present day and CMIP5 RCP8.5 climate change anomalies. Evaluation of the historical simulation shows an EAC extension that is stronger than similar ocean‐only models and observations. This bias is not explained by a linear response to differences in wind stress. The climate change simulations show that regional atmospheric CMIP5 MMM anomalies drive 73% of the projected 12 Sv increase in EAC extension transport whereas the remote ocean boundary conditions and regional radiative forcing (greenhouse gases within the domain) play a smaller role. The importance of regional changes in wind stress curl in driving the enhanced EAC extension is consistent with linear theory where the NEMO‐OASIS‐WRF response is closer to linear transport estimates compared to the CMIP5 MMM

DCMIP2016: a review of non-hydrostatic dynamical core design and intercomparison of participating models

Atmospheric dynamical cores are a fundamental component of global atmospheric modeling systems and are responsible for capturing the dynamical behavior of the Earth's atmosphere via numerical integration of the Navier-Stokes equations. These systems have existed in one form or another for over half of a century, with the earliest discretizations having now evolved into a complex ecosystem of algorithms and computational strategies. In essence, no two dynamical cores are alike, and their individual successes suggest that no perfect model exists. To better understand modern dynamical cores, this paper aims to provide a comprehensive review of 11 non-hydrostatic dynamical cores, drawn from modeling centers and groups that participated in the 2016 Dynamical Core Model Intercomparison Project (DCMIP) workshop and summer school. This review includes a choice of model grid, variable placement, vertical coordinate, prognostic equations, temporal discretization, and the diffusion, stabilization, filters, and fixers employed by each syste

Recent Developments in Understanding Two-dimensional Turbulence and the Nastrom-Gage Spectrum

Two-dimensional turbulence appears to be a more formidable problem than three-dimensional turbulence despite the numerical advantage of working with one less dimension. In the present paper we review recent numerical investigations of the phenomenology of two-dimensional turbulence as well as recent theoretical breakthroughs by various leading researchers. We also review efforts to reconcile the observed energy spectrum of the atmosphere (the spectrum) with the predictions of two-dimensional turbulence and quasi-geostrophic turbulence.Comment: Invited review; accepted by J. Low Temp. Phys.; Proceedings for Warwick Turbulence Symposium Workshop on Universal features in turbulence: from quantum to cosmological scales, 200

DCMIP2016: the splitting supercell test case

This paper describes the splitting supercell idealized test case used in the 2016 Dynamical Core Model Intercomparison Project (DCMIP2016). These storms are useful test beds for global atmospheric models because the horizontal scale of convective plumes is O(1&thinsp;km), emphasizing non-hydrostatic dynamics. The test case simulates a supercell on a reduced-radius sphere with nominal resolutions ranging from 4 to 0.5&thinsp;km and is based on the work of Klemp et al. (2015). Models are initialized with an atmospheric environment conducive to supercell formation and forced with a small thermal perturbation. A simplified Kessler microphysics scheme is coupled to the dynamical core to represent moist processes. Reference solutions for DCMIP2016 models are presented. Storm evolution is broadly similar between models, although differences in the final solution exist. These differences are hypothesized to result from different numerical discretizations, physics–dynamics coupling, and numerical diffusion. Intramodel solutions generally converge as models approach 0.5&thinsp;km resolution, although exploratory simulations at 0.25&thinsp;km imply some dynamical cores require more refinement to fully converge. These results can be used as a reference for future dynamical core evaluation, particularly with the development of non-hydrostatic global models intended to be used in convective-permitting regimes.</p

Lateral terrestrial water flow contribution to summer precipitation at continental scale – A comparison between Europe and West Africa with WRF‐Hydro‐tag ensembles

It is well accepted that summer precipitation can be altered by soil moisture condition. Coupled land surface – atmospheric models have been routinely used to quantify soil moisture – precipitation feedback processes. However, most of the land surface models (LSMs) assume a vertical soil water transport and neglect lateral terrestrial water flow at the surface and in the subsurface, which potentially reduces the realism of the simulated soil moisture – precipitation feedback. In this study, the contribution of lateral terrestrial water flow to summer precipitation is assessed in two different climatic regions, Europe and West Africa, for the period June–September 2008. A version of the coupled atmospheric-hydrological model WRF-Hydro with an option to tag and trace land surface evaporation in the modelled atmosphere, named WRF-Hydro-tag, is employed. An ensemble of 30 simulations with terrestrial routing and 30 simulations without terrestrial routing is generated with random realizations of turbulent energy with the stochastic kinetic energy backscatter scheme, for both Europe and West Africa. The ensemble size allows to extract random noise from continental-scale averaged modelled precipitation. It is found that lateral terrestrial water flow increases the relative contribution of land surface evaporation to precipitation by 3.6% in Europe and 5.6% in West Africa, which enhances a positive soil moisture – precipitation feedback and generates more uncertainty in modelled precipitation, as diagnosed by a slight increase in normalized ensemble spread. This study demonstrates the small but non-negligible contribution of lateral terrestrial water flow to precipitation at continental scale