ALBERT

Description: Efficient performance of the Met Office Unified Model v8.2 on Intel Xeon partially used nodes Geoscientific Model Development Discussions, 7, 7395-7425, 2014 Author(s): I. Bermous The atmospheric Unified Model (UM) developed at the UK Met Office is used for weather and climate prediction by forecast teams at a number of international meteorological centres and research institutes on a wide variety of hardware and software environments. Over its 25 year history the UM sources have been optimised for a better application performance on a number of HPC systems including NEC SX vector architecture systems and recently the IBM Power6/Power7 platforms. Understanding the influence of the compiler flags, MPI libraries and run configurations is crucial to achieving the shortest elapsed times for a UM application on any particular HPC system. These aspects are very important for applications that must run within operational time frames. Driving the current study is the HPC industry trend since 1980 for processor arithmetic performance to increase at a faster rate than memory bandwidth. This gap has been growing especially fast for multicore processors in the past 10 years and it can have significant implication for the performance and performance scaling of memory bandwidth intensive applications, such as the UM. Analysis of partially used nodes on Intel Xeon clusters is provided in this paper for short and medium range weather forecasting systems using global and limited-area configurations. It is shown that on the Intel Xeon based clusters the fastest elapsed times and the most efficient system usage can be achieved using partially committed nodes.

Description: Simulation of groundwater and surface water over the continental US using a hyperresolution, integrated hydrologic model Geoscientific Model Development Discussions, 7, 7317-7349, 2014 Author(s): R. M. Maxwell, L. E. Condon, and S. J. Kollet Interactions between surface and groundwater systems are well-established theoretically and observationally. While numerical models that solve both surface and subsurface flow equations in a single framework (matrix) are increasingly being applied, computational limitations have restricted their use to local and regional studies. Regional or watershed, scale simulations have been effective tools in understanding hydrologic processes, however there are still many questions, such as the adaptation of water resources to anthropogenic stressors and climate variability, that need to be answered across large spatial extents at high resolution. In response to this "grand challenge" in hydrology, we present the results of a parallel, integrated hydrologic model simulating surface and subsurface flow at high spatial resolution (1 km) over much of continental North America (~ 6 300 000 or 6.3 million km 2 ). These simulations provide predictions of hydrologic states and fluxes, namely water table depth and streamflow, at unprecedented scale and resolution. The physically-based modeling approach used here requires limited parameterizations and relies only on more fundamental inputs, such as topography, hydrogeologic properties and climate forcing. Results are compared to observations and provide mechanistic insight into hydrologic process interaction. This study demonstrates both the feasibility of continental scale integrated models and their utility for improving our understanding of large-scale hydrologic systems; the combination of high resolution and large spatial extent facilitates novel analysis of scaling relationships using model outputs.

Description: EDDA: integrated simulation of debris flow erosion, deposition and property changes Geoscientific Model Development Discussions, 7, 7267-7316, 2014 Author(s): H. X. Chen and L. M. Zhang Debris flow material properties change during the initiation, transportation and deposition processes, which influences the runout characteristics of the debris flow. A quasi-three-dimensional depth-integrated numerical model, EDDA, is presented in this paper to simulate debris flow erosion, deposition and induced material property changes. The model considers changes in debris flow density, yield stress and dynamic viscosity during the flow process. The yield stress of debris flow mixture is determined at limit equilibrium using the Mohr–Coulomb equation, which is applicable to clear water flow, hyper-concentrated flow and fully developed debris flow. To assure numerical stability and computational efficiency at the same time, a variable time stepping algorithm is developed to solve the governing differential equations. Four numerical tests are conducted to validate the model. The first two tests involve a one-dimensional dam-break water flow and a one-dimensional debris flow with constant properties. The last two tests involve erosion and deposition, and the movement of multi-directional debris flows. The changes in debris flow mass and properties due to either erosion or deposition are shown to affect the runout characteristics significantly. The model is also applied to simulate a large-scale debris flow in Xiaojiagou Ravine to test the performance of the model in catchment-scale simulations. The results suggest that the model estimates well the volume, inundated area, and runout distance of the debris flow. The model is intended for use as a module in a real-time debris flow warning system.

Description: Northern Hemisphere storminess in the Norwegian Earth System Model (NorESM1-M) Geoscientific Model Development Discussions, 7, 8975-9015, 2014 Author(s): E. M. Knudsen and J. E. Walsh Metrics of storm activity in Northern Hemisphere high- and midlatitudes are evaluated from historical output and future projections by the Norwegian Earth System Model (NorESM1-M) coupled global climate model. The European Re-Analysis Interim (ERA-Interim) and the Community Climate System Model (CCSM4), a global climate model of the same vintage as NorESM1-M, provide benchmarks for comparison. The focus is on the autumn and early winter (September through December), the period when the ongoing and projected Arctic sea ice retreat is greatest. Storm tracks derived from a vorticity-based algorithm for storm identification are reproduced well by NorESM1-M, although the tracks are somewhat better resolved in the higher-resolution ERA-Interim and CCSM4. The tracks are projected to shift polewards in the future as climate changes under the Representative Concentration Pathway (RCP) forcing scenarios. Cyclones are projected to become generally more intense in the high-latitudes, especially over the Alaskan region, although in some other areas the intensity is projected to decrease. While projected changes in track density are less coherent, there is a general tendency towards less frequent storms in midlatitudes and more frequent storms in high-latitudes, especially the Baffin Bay/Davis Strait region. Autumn precipitation is projected to increase significantly across the entire high-latitudes. Together with the projected increases in storm intensity and sea level and the loss of sea ice, this increase in precipitation implies a greater vulnerability to coastal flooding and erosion, especially in the Alaskan region. The projected changes in storm intensity and precipitation (as well as sea ice and sea level pressure) scale generally linearly with the RCP value of the forcing and with time through the 21st century.

Description: A multi-layer land surface energy budget model for implicit coupling with global atmospheric simulations Geoscientific Model Development Discussions, 7, 8649-8701, 2014 Author(s): J. Ryder, J. Polcher, P. Peylin, C. Ottlé, Y. Chen, E. van Gorsel, V. Haverd, M. J. McGrath, K. Naudts, J. Otto, A. Valade, and S. Luyssaert In Earth system modelling, a description of the energy budget of the vegetated surface layer is fundamental as it determines the meteorological conditions in the planetary boundary layer and as such contributes to the atmospheric conditions and its circulation. The energy budget in most Earth system models has long been based on a "big-leaf approach", with averaging schemes that represent in-canopy processes. Such models have difficulties in reproducing consistently the energy balance in field observations. We here outline a newly developed numerical model for energy budget simulation, as a component of the land surface model ORCHIDEE-CAN (Organising Carbon and Hydrology In Dynamic Ecosystems – CANopy). This new model implements techniques from single-site canopy models in a practical way. It includes representation of in-canopy transport, a multilayer longwave radiation budget, height-specific calculation of aerodynamic and stomatal conductance, and interaction with the bare soil flux within the canopy space. Significantly, it avoids iterations over the height of tha canopy and so maintains implicit coupling to the atmospheric model LMDz. As a first test, the model is evaluated against data from both an intensive measurement campaign and longer term eddy covariance measurements for the intensively studied Eucalyptus stand at Tumbarumba, Australia. The model performs well in replicating both diurnal and annual cycles of fluxes, as well as the gradients of sensible heat fluxes. However, the model overestimates sensible heat flux against an underestimate of the radiation budget. Improved performance is expected through the implementation of a more detailed calculation of stand albedo and a more up-to-date stomatal conductance calculation.

Description: The Secondary Organic Aerosol Processor (SOAP v1.0) model: a unified model with different ranges of complexity based on the molecular surrogate approach Geoscientific Model Development Discussions, 7, 379-429, 2014 Author(s): F. Couvidat and K. Sartelet The Secondary Organic Aerosol Processor (SOAP v1.0) model is presented. This model is designed to be modular with different user options depending on the computing time and the complexity required by the user. This model is based on the molecular surrogate approach, in which each surrogate compound is associated with a molecular structure to estimate some properties and parameters (hygroscopicity, absorption on the aqueous phase of particles, activity coefficients, phase separation). Each surrogate can be hydrophilic (condenses only on the aqueous phase of particles), hydrophobic (condenses only on the organic phase of particles) or both (condenses on both the aqueous and the organic phases of particles). Activity coefficients are computed with the UNIFAC thermodynamic model for short-range interactions and with the AIOMFAC parameterization for medium and long-range interactions between electrolytes and organic compounds. Phase separation is determined by Gibbs energy minimization. The user can choose between an equilibrium and a dynamic representation of the organic aerosol. In the equilibrium representation, compounds in the particle phase are assumed to be at equilibrium with the gas phase. However, recent studies show that the organic aerosol (OA) is not at equilibrium with the gas phase because the organic phase could be semi-solid (very viscous liquid phase). The condensation or evaporation of organic compounds could then be limited by the diffusion in the organic phase due to the high viscosity. A dynamic representation of secondary organic aerosols (SOA) is used with OA divided into layers, the first layer at the center of the particle (slowly reaches equilibrium) and the final layer near the interface with the gas phase (quickly reaches equilibrium).

Description: Comparison of the ensemble Kalman filter and 4D-Var assimilation methods using a stratospheric tracer transport model Geoscientific Model Development Discussions, 7, 339-377, 2014 Author(s): S. Skachko, Q. Errera, R. Ménard, Y. Christophe, and S. Chabrillat The Ensemble Kalman filter (EnKF) assimilation method is applied to the tracer transport using the same stratospheric transport model as in the 4D-Var assimilation system BASCOE. This EnKF version of BASCOE was built primarily to avoid the large costs associated with the maintenance of an adjoint model. The EnKF developed in BASCOE accounts for two adjustable parameters: a parameter α controlling the model error term and a parameter r controlling the observational error. The EnKF system is shown to be markedly sensitive to these two parameters, which are adjusted based on the monitoring of a χ 2 -test measuring the misfit between the control variable and the observations. The performance of the EnKF and 4D-Var versions was estimated through the assimilation of Aura-MLS ozone observations during an 8 month period which includes the formation of the 2008 Antarctic ozone hole. To ensure a proper comparison, despite the fundamental differences between the two assimilation methods, both systems use identical and carefully calibrated input error statistics. We provide the detailed procedure for these calibrations, and compare the two sets of analyses with a focus on the lower and middle stratosphere where the ozone lifetime is much larger than the observational update frequency. Based on the Observation-minus-Forecast statistics, we show that the analyses provided by the two systems are markedly similar, with biases smaller than 5% and standard deviation errors smaller than 10% in most of the stratosphere. Since the biases are markedly similar, they have most probably the same causes: these can be deficiencies in the model and in the observation dataset, but not in the assimilation algorithm nor in the error calibration. The remarkably similar performance also shows that in the context of stratospheric transport, the choice of the assimilation method can be based on application-dependent factors, such as CPU cost or the ability to generate an ensemble of forecasts.

Description: Regional climate modeling on European scales: a joint standard evaluation of the EURO-CORDEX RCM ensemble Geoscientific Model Development Discussions, 7, 217-293, 2014 Author(s): S. Kotlarski, K. Keuler, O. B. Christensen, A. Colette, M. Déqué, A. Gobiet, K. Goergen, D. Jacob, D. Lüthi, E. van Meijgaard, G. Nikulin, C. Schär, C. Teichmann, R. Vautard, K. Warrach-Sagi, and V. Wulfmeyer EURO-CORDEX is an international climate downscaling initiative that aims to provide high-resolution climate scenarios for Europe. Here an evaluation of the ERA-Interim-driven EURO-CORDEX regional climate model (RCM) ensemble is presented. The study documents the performance of the individual models in representing the basic spatio-temporal patterns of the European climate for the period 1989–2008. Model evaluation focuses on near-surface air temperature and precipitation, and uses the E-OBS dataset as observational reference. The ensemble consists of 17 simulations carried out by seven different models at grid resolutions of 12 km (nine experiments) and 50 km (eight experiments). Several performance metrics computed from monthly and seasonal mean values are used to assess model performance over eight sub-domains of the European continent. Results are compared to those for the ERA40-driven ENSEMBLES simulations. The analysis confirms the ability of RCMs to capture the basic features of the European climate, including its variability in space and time. But it also identifies non-negligible deficiencies of the simulations for selected metrics, regions and seasons. Seasonally and regionally averaged temperature biases are mostly smaller than 1.5 °C, while precipitation biases are typically located in the ±40% range. Some bias characteristics, such as a predominant cold and wet bias in most seasons and over most parts of Europe and a warm and dry summer bias over southern and south-eastern Europe reflect common model biases. For seasonal mean quantities averaged over large European sub-domains, no clear benefit of an increased spatial resolution (12 km vs. 50 km) can be identified. The bias ranges of the EURO-CORDEX ensemble mostly correspond to those of the ENSEMBLES simulations, but some improvements in model performance can be identified (e.g., a less pronounced southern European warm summer bias). The temperature bias spread across different configurations of one individual model can be of a similar magnitude as the spread across different models, demonstrating a strong influence of the specific choices in physical parameterizations and experimental setup on model performance. Based on a number of simply reproducible metrics, the present study quantifies the currently achievable accuracy of RCMs used for regional climate simulations over Europe and provides a quality standard for future model developments.

Description: Simulation of trace gases and aerosols over the Indian domain: evaluation of the WRF-Chem model Geoscientific Model Development Discussions, 7, 431-482, 2014 Author(s): M. Michael, A. Yadav, S. N. Tripathi, V. P. Kanawade, A. Gaur, P. Sadavarte, and C. Venkataraman The "online" meteorological and chemical transport Weather Research and Forecasting/Chemistry (WRF-Chem) model has been implemented over the Indian subcontinent for three consecutive summers in 2008, 2009 and 2010 to study the aerosol properties over the domain. The model simulated the meteorological parameters, trace gases and particulate matter. Predicted mixing ratios of trace gases (Ozone, carbon monoxide and sulfur dioxide) are compared with ground based observations over Kanpur. Simulated aerosol optical depth are compared with those observed at nine Aerosol Robotic Network stations (AERONET). The simulations show that the aerosol optical depth of the less polluted regions is better simulated compared to that of the locations where the aerosol loading is very high. The vertical profiles of extinction coefficient observed at the Kanpur Micropulse Lidar Network (MPLNET) station is underpredicted by the model by 10 to 50% for altitudes greater than 1.5 km and qualitatively simulate the elevated layers of aerosols. The simulated mass concentration of black carbon shows a correlation coefficient of 0.4 with observations. Vertical profiles of black carbon at various locations have also been compared with observations from an aircraft campaign held during pre-monsoon period of 2008 and 2009. This study shows that WRF-Chem model captures many important features of the observed atmospheric composition during the pre-monsoon season in India.

Description: Three-dimensional phase-field study of crack-seal microstructures – insights from innovative post-processing techniques Geoscientific Model Development Discussions, 7, 631-658, 2014 Author(s): K. Ankit, M. Selzer, and B. Nestler Numerical simulations of vein evolution contribute to a better understanding of processes involved in their formation and possess the potential to provide invaluable insights into the rock deformation history and fluid flow pathways. The primary aim of the present article is to investigate the influence of a "realistic" boundary condition, i.e. an algorithmically generated "fractal" surface, on the vein evolution in 3-D using a thermodynamically consistent approach, while explaining the benefits of accounting for an extra dimensionality. The 3-D simulation results are supplemented by innovative numerical post-processing and advanced visualization techniques. The new methodologies to measure the tracking efficiency demonstrate the importance of accounting the temporal evolution; no such information is usually accessible in field studies and notoriously difficult to obtain from laboratory experiments as well. The grain growth statistics obtained by numerically post-processing the 3-D computational microstructures explain the pinning mechanism which leads to arrest of grain boundaries/multi-junctions by crack peaks, thereby, enhancing the tracking behavior.