ALBERT

Description: The Chemistry CATT-BRAMS model (CCATT-BRAMS 4.5): a regional atmospheric model system for integrated air quality and weather forecasting and research Geoscientific Model Development, 6, 1389-1405, 2013 Author(s): K. M. Longo, S. R. Freitas, M. Pirre, V. Marécal, L. F. Rodrigues, J. Panetta, M. F. Alonso, N. E. Rosário, D. S. Moreira, M. S. Gácita, J. Arteta, R. Fonseca, R. Stockler, D. M. Katsurayama, A. Fazenda, and M. Bela Coupled Chemistry Aerosol-Tracer Transport model to the Brazilian developments on the Regional Atmospheric Modeling System (CCATT-BRAMS, version 4.5) is an on-line regional chemical transport model designed for local and regional studies of atmospheric chemistry from the surface to the lower stratosphere suitable both for operational and research purposes. It includes gaseous/aqueous chemistry, photochemistry, scavenging and dry deposition. The CCATT-BRAMS model takes advantage of BRAMS-specific development for the tropics/subtropics as well as the recent availability of preprocessing tools for chemical mechanisms and fast codes for photolysis rates. BRAMS includes state-of-the-art physical parameterizations and dynamic formulations to simulate atmospheric circulations down to the meter. This on-line coupling of meteorology and chemistry allows the system to be used for simultaneous weather and chemical composition forecasts as well as potential feedback between the two. The entire system is made of three preprocessing software tools for user-defined chemical mechanisms, aerosol and trace gas emissions fields and the interpolation of initial and boundary conditions for meteorology and chemistry. In this paper, the model description is provided along with the evaluations performed by using observational data obtained from ground-based stations, instruments aboard aircrafts and retrieval from space remote sensing. The evaluation accounts for model applications at different scales from megacities and the Amazon Basin up to the intercontinental region of the Southern Hemisphere.

Description: Beyond multifractional Brownian motion: new stochastic models for geophysical modelling Nonlinear Processes in Geophysics, 20, 643-655, 2013 Author(s): J. Lévy Véhel Multifractional Brownian motion (mBm) has proved to be a useful tool in various areas of geophysical modelling. Although a versatile model, mBm is of course not always an adequate one. We present in this work several other stochastic processes which could potentially be useful in geophysics. The first alternative type is that of self-regulating processes : these are models where the local regularity is a function of the amplitude, in contrast to mBm where it is tuned exogenously. We demonstrate the relevance of such models for digital elevation maps and for temperature records. We also briefly describe two other types of alternative processes, which are the counterparts of mBm and of self-regulating processes when the intensity of local jumps is considered in lieu of local regularity: multistable processes allow one to prescribe the local intensity of jumps in space/time, while this intensity is governed by the amplitude for self-stabilizing processes .

Description: δ 18 O water isotope in the i LOVECLIM model (version 1.0) – Part 1: Implementation and verification Geoscientific Model Development, 6, 1481-1491, 2013 Author(s): D. M. Roche A new 18 O stable water isotope scheme is developed for three components of the i LOVECLIM coupled climate model: atmospheric, oceanic and land surface. The equations required to reproduce the fractionation of stable water isotopes in the simplified atmospheric model ECBilt are developed consistently with the moisture scheme. Simplifications in the processes are made to account for the simplified vertical structure including only one moist layer. Implementation of these equations together with a passive tracer scheme for the ocean and a equilibrium fractionation scheme for the land surface leads to the closure of the (isotopic-) water budget in our climate system. Following the implementation, verification of the existence of usual δ 18 O to climatic relationships are performed for the Rayleigh distillation, the Dansgaard relationship and the δ 18 O –salinity relationship. Advantages and caveats of the approach taken are outlined. The isotopic fields simulated are shown to reproduce most expected oxygen-18–climate relationships with the notable exception of the isotopic composition in Antarctica.

Description: δ 18 O water isotope in the i LOVECLIM model (version 1.0) – Part 3: A palaeo-perspective based on present-day data–model comparison for oxygen stable isotopes in carbonates Geoscientific Model Development, 6, 1505-1516, 2013 Author(s): T. Caley and D. M. Roche Oxygen stable isotopes (δ 18 O) are among the most useful tools in palaeoclimatology/palaeoceanography. Simulation of oxygen stable isotopes allows testing how the past variability of these isotopes in water can be interpreted. By modelling the proxy directly in the model, the results can also be directly compared with the data. Water isotopes have been implemented in the global three-dimensional model of intermediate complexity i LOVECLIM, allowing fully coupled atmosphere–ocean simulations. In this study, we present the validation of the model results for present-day climate against the global database for oxygen stable isotopes in carbonates. The limitation of the model together with the processes operating in the natural environment reveal the complexity of use the continental calcite-δ 18 O signal of speleothems for a global quantitative data–model comparison exercise. On the contrary, the reconstructed surface ocean calcite-δ 18 O signal in i LOVECLIM does show a very good agreement with the late Holocene database (foraminifers) at the global and regional scales. Our results indicate that temperature and the isotopic composition of the seawater are the main control on the fossil-δ 18 O signal recorded in foraminifer shells when all species are grouped together. Depth habitat, seasonality and other ecological effects play a more significant role when individual species are considered. We argue that a data–model comparison for surface ocean calcite δ 18 O in past climates, such as the Last Glacial Maximum (≈ 21 000 yr), could constitute an interesting tool for mapping the potential shifts of the frontal systems and circulation changes throughout time. Similarly, the potential changes in intermediate oceanic circulation systems in the past could be documented by a data (benthic foraminifers)-model comparison exercise whereas future investigations are necessary in order to quantitatively compare the results with data for the deep ocean.

Description: Assimilating water column and satellite data for marine export production estimation Geoscientific Model Development, 6, 1575-1590, 2013 Author(s): X. Yao and R. Schlitzer Recent advances in satellite retrieval methodology now allow for estimation of particular organic carbon (POC) concentration in ocean surface waters directly from satellite-based optical data. Because of the good coverage, these data reveal small-scale spatial and temporal concentration gradients and document the evolution of surface water POC as well as the underlying driving biogeochemical processes throughout the seasons. Water column nutrient data also reveal biogeochemical activity. However, because of the scarcity of data, the deduction of temporal changes of particle production and export is not possible in most parts of the ocean. Here we present first results from a new study combining both data streams, thereby exploiting the high spatio-temporal resolution of surface POC concentrations from satellite optical sensors with water column nutrient data having sparser coverage but providing information throughout the entire water column. We use a medium-resolution global model with steady-state 3-D circulation that has been optimized by fitting to a large number of hydrographic parameters and tracers, including CFCs and natural radiocarbon. Production and export of POC is allowed to vary monthly, and the magnitudes of the monthly export fluxes are determined by fitting the model to satellite POC data as well as water column nutrient data using the adjoint method. Two cases have been investigated: (1) the production rate of POC is set to be proportional to export production (EP) and the seasonal changes are assumed sinusoidal (meridionally varying amplitude and phase), and (2) the POC production rate is linked to primary production rates (literature). Both cases were run with the same initial state and model settings, and show total cost function decreases of 12 and 95%, respectively. The POC misfit term alone decreased by 75 and 99.8%. The integrated annual global POC exports of the two cases are 9.9 and 12.3 Gt C yr −1 , respectively. Overall, the remaining POC and phosphate misfits of both solutions are considered too large, and the difference fields still exhibit significant systematic geographical patterns. This indicates that the present model runs are too simplistic and do not fully explain the data. Further, more refined model setups are needed.

Description: Enhancing the representation of subgrid land surface characteristics in land surface models Geoscientific Model Development, 6, 1609-1622, 2013 Author(s): Y. Ke, L. R. Leung, M. Huang, and H. Li Land surface heterogeneity has long been recognized as important to represent in the land surface models. In most existing land surface models, the spatial variability of surface cover is represented as subgrid composition of multiple surface cover types, although subgrid topography also has major controls on surface processes. In this study, we developed a new subgrid classification method (SGC) that accounts for variability of both topography and vegetation cover. Each model grid cell was represented with a variable number of elevation classes and each elevation class was further described by a variable number of vegetation types optimized for each model grid given a predetermined total number of land response units (LRUs). The subgrid structure of the Community Land Model (CLM) was used to illustrate the newly developed method in this study. Although the new method increases the computational burden in the model simulation compared to the CLM subgrid vegetation representation, it greatly reduced the variations of elevation within each subgrid class and is able to explain at least 80% of the total subgrid plant functional types (PFTs). The new method was also evaluated against two other subgrid methods (SGC1 and SGC2) that assigned fixed numbers of elevation and vegetation classes for each model grid (SGC1: M elevation bands– N PFTs method; SGC2: N PFTs– M elevation bands method). Implemented at five model resolutions (0.1°, 0.25°, 0.5°, 1.0°and 2.0°) with three maximum-allowed total number of LRUs (i.e., N LRU of 24, 18 and 12) over North America (NA), the new method yielded more computationally efficient subgrid representation compared to SGC1 and SGC2, particularly at coarser model resolutions and moderate computational intensity ( N LRU = 18). It also explained the most PFTs and elevation variability that is more homogeneously distributed spatially. The SGC method will be implemented in CLM over the NA continent to assess its impacts on simulating land surface processes.

Description: Non-extensivity and long-range correlations in the earthquake activity at the West Corinth rift (Greece) Nonlinear Processes in Geophysics, 20, 713-724, 2013 Author(s): G. Michas, F. Vallianatos, and P. Sammonds In the present work the statistical properties of the earthquake activity in a highly seismic region, the West Corinth rift (Central Greece), are being studied by means of generalized statistical physics. By using a dataset that covers the period 2001–2008, we investigate the earthquake energy distribution and the distribution of the time intervals (interevent times) between the successive events. As has been reported previously, these distributions exhibit complex statistical properties and fractality. By using detrended fluctuation analysis (DFA), a well-established method for detection of long-range correlations in non-stationary signals, it is shown that long-range correlations are also present in the earthquake activity. The existence of these properties motivates us to use non-extensive statistical physics (NESP) to investigate the statistical properties of the frequency-magnitude and the interevent time distributions, along with other well-known relations in seismology, such as the gamma distribution for interevent times. The results of the analysis indicate that the statistical properties of the earthquake activity can be successfully reproduced by means of NESP and that the earthquake activity at the West Corinth rift is correlated at all-time scales.

Description: An efficient method to generate a perturbed parameter ensemble of a fully coupled AOGCM without flux-adjustment Geoscientific Model Development, 6, 1447-1462, 2013 Author(s): P. J. Irvine, L. J. Gregoire, D. J. Lunt, and P. J. Valdes We present a simple method to generate a perturbed parameter ensemble (PPE) of a fully-coupled atmosphere-ocean general circulation model (AOGCM), HadCM3, without requiring flux-adjustment. The aim was to produce an ensemble that samples parametric uncertainty in some key variables and gives a plausible representation of the climate. Six atmospheric parameters, a sea-ice parameter and an ocean parameter were jointly perturbed within a reasonable range to generate an initial group of 200 members. To screen out implausible ensemble members, 20 yr pre-industrial control simulations were run and members whose temperature responses to the parameter perturbations were projected to be outside the range of 13.6 ± 2 °C, i.e. near to the observed pre-industrial global mean, were discarded. Twenty-one members, including the standard unperturbed model, were accepted, covering almost the entire span of the eight parameters, challenging the argument that without flux-adjustment parameter ranges would be unduly restricted. This ensemble was used in 2 experiments; an 800 yr pre-industrial and a 150 yr quadrupled CO 2 simulation. The behaviour of the PPE for the pre-industrial control compared well to ERA-40 reanalysis data and the CMIP3 ensemble for a number of surface and atmospheric column variables with the exception of a few members in the Tropics. However, we find that members of the PPE with low values of the entrainment rate coefficient show very large increases in upper tropospheric and stratospheric water vapour concentrations in response to elevated CO 2 and one member showed an implausible nonlinear climate response, and as such will be excluded from future experiments with this ensemble. The outcome of this study is a PPE of a fully-coupled AOGCM which samples parametric uncertainty and a simple methodology which would be applicable to other GCMs.

Description: The impact of initial spread calibration on the RELO ensemble and its application to Lagrangian dynamics Nonlinear Processes in Geophysics, 20, 621-641, 2013 Author(s): M. Wei, G. Jacobs, C. Rowley, C. N. Barron, P. Hogan, P. Spence, O. M. Smedstad, P. Martin, P. Muscarella, and E. Coelho A number of real-time ocean model forecasts were carried out successfully at Naval Research Laboratory (NRL) to provide modeling support and numerical guidance to the CARTHE GLAD at-sea experiment during summer 2012. Two RELO ensembles and three single models using NCOM and HYCOM with different resolutions were carried out. A calibrated ensemble system with enhanced spread and reliability was developed to better support this experiment. The calibrated ensemble is found to outperform the un-calibrated ensemble in forecasting accuracy, skill, and reliability for all the variables and observation spaces evaluated. The metrics used in this paper include RMS error, anomaly correlation, PECA, Brier score, spread reliability, and Talagrand rank histogram. It is also found that even the un-calibrated ensemble outperforms the single forecast from the model with the same resolution. The advantages of the ensembles are further extended to the Lagrangian framework. In contrast to a single model forecast, the RELO ensemble provides not only the most likely Lagrangian trajectory for a particle in the ocean, but also an uncertainty estimate that directly reflects the complicated ocean dynamics, which is valuable for decision makers. The examples show that the calibrated ensemble with more reliability can capture trajectories in different, even opposite, directions, which would be missed by the un-calibrated ensemble. The ensembles are applied to compute the repelling and attracting Lagrangian coherent structures (LCSs), and the uncertainties of the LCSs, which are hard to obtain from a single model forecast, are estimated. It is found that the spatial scales of the LCSs depend on the model resolution. The model with the highest resolution produces the finest, small-scale, LCS structures, while the model with lowest resolution generates only large-scale LCSs. The repelling and attracting LCSs are found to intersect at many locations and create complex mesoscale eddies. The fluid particles and drifters in the middle of these tangles are subject to attraction and repulsion simultaneously from these two kinds of LCSs. As a result, the movements of particles near the Deepwater Horizon (DWH) location are severely limited. This is also confirmed by the Lagrangian trajectories predicted by the ensembles.

Description: Multifractal properties of embedded convective structures in orographic precipitation: toward subgrid-scale predictability Nonlinear Processes in Geophysics, 20, 605-620, 2013 Author(s): M. Nogueira, A. P. Barros, and P. M. A. Miranda Rain and cloud fields produced by fully nonlinear idealized cloud resolving numerical simulations of orographic convective precipitation display statistical multiscaling behavior, implying that multifractal diagnostics should provide a physically robust basis for the downscaling and sub-grid scale parameterizations of moist processes. Our results show that the horizontal scaling exponent function (and respective multiscaling parameters) of the simulated rainfall and cloud fields varies with atmospheric and terrain properties, particularly small-scale terrain spectra, atmospheric stability, and advective timescale. This implies that multifractal diagnostics of moist processes for these simulations are fundamentally transient, exhibiting complex nonlinear behavior depending on atmospheric conditions and terrain forcing at each location. A particularly robust behavior found here is the transition of the multifractal parameters between stable and unstable cases, which has a clear physical correspondence to the transition from stratiform to organized (banded and cellular) convective regime. This result is reinforced by a similar behavior in the horizontal spectral exponent. Finally, our results indicate that although nonlinearly coupled fields (such as rain and clouds) have different scaling exponent functions, there are robust relationships with physical underpinnings between the scaling parameters that can be explored for hybrid dynamical-statistical downscaling.

Description: GAPPARD: a computationally efficient method of approximating gap-scale disturbance in vegetation models Geoscientific Model Development, 6, 1517-1542, 2013 Author(s): M. Scherstjanoi, J. O. Kaplan, E. Thürig, and H. Lischke Models of vegetation dynamics that are designed for application at spatial scales larger than individual forest gaps suffer from several limitations. Typically, either a population average approximation is used that results in unrealistic tree allometry and forest stand structure, or models have a high computational demand because they need to simulate both a series of age-based cohorts and a number of replicate patches to account for stochastic gap-scale disturbances. The detail required by the latter method increases the number of calculations by two to three orders of magnitude compared to the less realistic population average approach. In an effort to increase the efficiency of dynamic vegetation models without sacrificing realism, we developed a new method for simulating stand-replacing disturbances that is both accurate and faster than approaches that use replicate patches. The GAPPARD (approximating GAP model results with a Probabilistic Approach to account for stand Replacing Disturbances) method works by postprocessing the output of deterministic, undisturbed simulations of a cohort-based vegetation model by deriving the distribution of patch ages at any point in time on the basis of a disturbance probability. With this distribution, the expected value of any output variable can be calculated from the output values of the deterministic undisturbed run at the time corresponding to the patch age. To account for temporal changes in model forcing (e.g., as a result of climate change), GAPPARD performs a series of deterministic simulations and interpolates between the results in the postprocessing step. We integrated the GAPPARD method in the vegetation model LPJ-GUESS, and evaluated it in a series of simulations along an altitudinal transect of an inner-Alpine valley. We obtained results very similar to the output of the original LPJ-GUESS model that uses 100 replicate patches, but simulation time was reduced by approximately the factor 10. Our new method is therefore highly suited for rapidly approximating LPJ-GUESS results, and provides the opportunity for future studies over large spatial domains, allows easier parameterization of tree species, faster identification of areas of interesting simulation results, and comparisons with large-scale datasets and results of other forest models.

Description: A top-down model to generate ensembles of runoff from a large number of hillslopes Nonlinear Processes in Geophysics, 20, 683-704, 2013 Author(s): P. R. Furey, V. K. Gupta, and B. M. Troutman We hypothesize that total hillslope water loss for a rainfall–runoff event is inversely related to a function of a lognormal random variable, based on basin- and point-scale observations taken from the 21 km 2 Goodwin Creek Experimental Watershed (GCEW) in Mississippi, USA. A top-down approach is used to develop a new runoff generation model both to test our physical-statistical hypothesis and to provide a method of generating ensembles of runoff from a large number of hillslopes in a basin. The model is based on the assumption that the probability distributions of a runoff/loss ratio have a space–time rescaling property. We test this assumption using streamflow and rainfall data from GCEW. For over 100 rainfall–runoff events, we find that the spatial probability distributions of a runoff/loss ratio can be rescaled to a new distribution that is common to all events. We interpret random within-event differences in runoff/loss ratios in the model to arise from soil moisture spatial variability. Observations of water loss during events in GCEW support this interpretation. Our model preserves water balance in a mean statistical sense and supports our hypothesis. As an example, we use the model to generate ensembles of runoff at a large number of hillslopes for a rainfall–runoff event in GCEW.

Description: A mechanism for catastrophic filter divergence in data assimilation for sparse observation networks Nonlinear Processes in Geophysics, 20, 705-712, 2013 Author(s): G. A. Gottwald and A. J. Majda We study catastrophic filter divergence in data assimilation procedures whereby the forecast model develops severe numerical instabilities leading to a blow-up of the solution. Catastrophic filter divergence can occur in sparse observational grids with small observational noise for intermediate observation intervals and finite ensemble sizes. Using a minimal five-dimensional model, we establish that catastrophic filter divergence is a numerical instability of the underlying forecast model caused by the filtering procedure producing analyses which are not consistent with the true dynamics, and stiffness caused by the fast attraction of the inconsistent analyses towards the attractor during the forecast step.

Description: Sub-inertial modulation of nonlinear Kelvin waves in the coastal zone Nonlinear Processes in Geophysics, 20, 357-364, 2013 Author(s): D. V. Stepanov and V. V. Novotryasov Observational evidence is presented for interaction between nonlinear internal Kelvin waves at the ω t,i (where the ω t is the semidiurnal frequency and the ω i is the inertial frequency) and random oscillations of the background coastal current at the sub-inertial Ω frequency in the Japan/East Sea. Enhanced coastal currents at the sum ω + and difference ω-frequencies ω ± =ω t,i ± Ω have properties of propagating Kelvin waves, which suggests permanent energy exchange from the sub-inertial band to the mesoscale ω ± band. This interaction may be responsible for a greater-than-predicted intensification, steepening and breaking of boundary-trapped Kelvin waves. The problem of interaction between the nonlinear Kelvin wave at the frequency ω and the low-frequency narrowband noise with representative frequency Ω≪ω is investigated using the theory of nonlinear weak dispersion waves.

Description: On the parallelization of atmospheric inversions of CO 2 surface fluxes within a variational framework Geoscientific Model Development, 6, 783-790, 2013 Author(s): F. Chevallier The variational formulation of Bayes' theorem allows inferring CO 2 sources and sinks from atmospheric concentrations at much higher time–space resolution than the ensemble or analytical approaches. However, it usually exhibits limited scalable parallelism. This limitation hinders global atmospheric inversions operated on decadal time scales and regional ones with kilometric spatial scales because of the computational cost of the underlying transport model that has to be run at each iteration of the variational minimization. Here, we introduce a physical parallelization (PP) of variational atmospheric inversions. In the PP, the inversion still manages a single physically and statistically consistent window, but the transport model is run in parallel overlapping sub-segments in order to massively reduce the computation wall-clock time of the inversion. For global inversions, a simplification of transport modelling is described to connect the output of all segments. We demonstrate the performance of the approach on a global inversion for CO 2 with a 32 yr inversion window (1979–2010) with atmospheric measurements from 81 sites of the NOAA global cooperative air sampling network. In this case, we show that the duration of the inversion is reduced by a seven-fold factor (from months to days), while still processing the three decades consistently and with improved numerical stability.

Description: Inclusion of ash and SO 2 emissions from volcanic eruptions in WRF-Chem: development and some applications Geoscientific Model Development, 6, 457-468, 2013 Author(s): M. Stuefer, S. R. Freitas, G. Grell, P. Webley, S. Peckham, S. A. McKeen, and S. D. Egan We describe a new functionality within the Weather Research and Forecasting (WRF) model with coupled Chemistry (WRF-Chem) that allows simulating emission, transport, dispersion, transformation and sedimentation of pollutants released during volcanic activities. Emissions from both an explosive eruption case and a relatively calm degassing situation are considered using the most recent volcanic emission databases. A preprocessor tool provides emission fields and additional information needed to establish the initial three-dimensional cloud umbrella/vertical distribution within the transport model grid, as well as the timing and duration of an eruption. From this source condition, the transport, dispersion and sedimentation of the ash cloud can be realistically simulated by WRF-Chem using its own dynamics and physical parameterization as well as data assimilation. Examples of model applications include a comparison of tephra fall deposits from the 1989 eruption of Mount Redoubt (Alaska) and the dispersion of ash from the 2010 Eyjafjallajökull eruption in Iceland. Both model applications show good coincidence between WRF-Chem and observations.

Description: PORT, a CESM tool for the diagnosis of radiative forcing Geoscientific Model Development, 6, 469-476, 2013 Author(s): A. J. Conley, J.-F. Lamarque, F. Vitt, W. D. Collins, and J. Kiehl The Parallel Offline Radiative Transfer (PORT) model is a stand-alone tool, driven by model-generated datasets, that can be used for any radiation calculation that the underlying radiative transfer schemes can perform, such as diagnosing radiative forcing. In its present distribution, PORT isolates the radiation code from the Community Atmosphere Model (CAM4) in the Community Earth System Model (CESM1). The current configuration focuses on CAM4 radiation with the constituents as represented in present-day conditions in CESM1, along with their optical properties. PORT includes an implementation of stratospheric temperature adjustment under the assumption of fixed dynamical heating, which is necessary to compute radiative forcing in addition to the more straightforward instantaneous radiative forcing. PORT can be extended to use radiative constituent distributions from other models or model simulations. Ultimately, PORT can be used with various radiative transfer models. As illustrations of the use of PORT, we perform the computation of radiative forcing from doubling of carbon dioxide, from the change of tropospheric ozone concentration from the year 1850 to 2000, and from present-day aerosols. The radiative forcing from tropospheric ozone (with respect to 1850) generated by a collection of model simulations under the Atmospheric Chemistry and Climate Model Intercomparison Project is found to be 0.34 (with an intermodel standard deviation of 0.07) W m −2 . Present-day aerosol direct forcing (relative to no aerosols) is found to be −1.3 W m −2 .

Description: Application of wavelet transform for evaluation of hydrocarbon reservoirs: example from Iranian oil fields in the north of the Persian Gulf Nonlinear Processes in Geophysics, 20, 231-238, 2013 Author(s): M. R. Saadatinejad and H. Hassani The Persian Gulf and its surrounding area are some of the biggest basins and have a very important role in producing huge amounts of hydrocarbon, and this potential was evaluated in order to explore the target for geoscientists and petroleum engineers. Wavelet transform is a useful and applicable technique to reveal frequency contents of various signals in different branches of science and especially in petroleum studies. We applied two major capacities of continuous mode of wavelet transform in seismic investigations. These investigations were operated to detect reservoir geological structures and some anomalies related to hydrocarbon to develop and explore new petroleum reservoirs in at least 4 oilfields in the southwest of Iran. It had been observed that continuous wavelet transform results show some discontinuities in the location of faults and are able to display them more clearly than other seismic methods. Moreover, continuous wavelet transform, utilizing Morlet wavelet, displays low-frequency shadows on 4 different iso-frequency vertical sections to identify reservoirs containing gas. By comparing these different figures, the presence of low-frequency shadows under the reservoir could be seen and we can relate these variations from anomalies at different frequencies as an indicator of the presence of hydrocarbons in the target reservoir.

Description: Downscaling a global climate model to simulate climate change over the US and the implication on regional and urban air quality Geoscientific Model Development, 6, 1429-1445, 2013 Author(s): M. Trail, A. P. Tsimpidi, P. Liu, K. Tsigaridis, Y. Hu, A. Nenes, and A. G. Russell Climate change can exacerbate future regional air pollution events by making conditions more favorable to form high levels of ozone. In this study, we use spectral nudging with the Weather Research and Forecasting (WRF) model to downscale NASA earth system GISS modelE2 results during the years 2006 to 2010 and 2048 to 2052 over the contiguous United States in order to compare the resulting meteorological fields from the air quality perspective during the four seasons of five-year historic and future climatological periods. GISS results are used as initial and boundary conditions by the WRF regional climate model (RCM) to produce hourly meteorological fields. The downscaling technique and choice of physics parameterizations used are evaluated by comparing them with in situ observations. This study investigates changes of similar regional climate conditions down to a 12 km by 12 km resolution, as well as the effect of evolving climate conditions on the air quality at major US cities. The high-resolution simulations produce somewhat different results than the coarse-resolution simulations in some regions. Also, through the analysis of the meteorological variables that most strongly influence air quality, we find consistent changes in regional climate that would enhance ozone levels in four regions of the US during fall (western US, Texas, northeastern, and southeastern US), one region during summer (Texas), and one region where changes potentially would lead to better air quality during spring (Northeast). Changes in regional climate that would enhance ozone levels are increased temperatures and stagnation along with decreased precipitation and ventilation. We also find that daily peak temperatures tend to increase in most major cities in the US, which would increase the risk of health problems associated with heat stress. Future work will address a more comprehensive assessment of emissions and chemistry involved in the formation and removal of air pollutants.

Description: The SOCOL version 3.0 chemistry–climate model: description, evaluation, and implications from an advanced transport algorithm Geoscientific Model Development, 6, 1407-1427, 2013 Author(s): A. Stenke, M. Schraner, E. Rozanov, T. Egorova, B. Luo, and T. Peter We present the third generation of the coupled chemistry–climate model (CCM) SOCOL (modeling tools for studies of SOlar Climate Ozone Links). The most notable modifications compared to the previous model version are (1) the dynamical core has been updated with the fifth generation of the middle-atmosphere general circulation model MA-ECHAM (European Centre/HAMburg climate model), and (2) the advection of the chemical species is now calculated by a mass-conserving and shape-preserving flux-form transport scheme instead of the previously used hybrid advection scheme. The whole chemistry code has been rewritten according to the ECHAM5 infrastructure and transferred to Fortran95. In contrast to its predecessors, SOCOLvs3 is now fully parallelized. The performance of the new SOCOL version is evaluated on the basis of transient model simulations (1975–2004) with different horizontal (T31 and T42) resolutions, following the approach of the CCMVal-1 model validation activity. The advanced advection scheme significantly reduces the artificial loss and accumulation of tracer mass in regions with strong gradients that was observed in previous model versions. Compared to its predecessors, SOCOLvs3 generally shows more realistic distributions of chemical trace species, especially of total inorganic chlorine, in terms of the mean state, but also of the annual and interannual variability. Advancements with respect to model dynamics are for example a better representation of the stratospheric mean state in spring, especially in the Southern Hemisphere, and a slowdown of the upward propagation in the tropical lower stratosphere. Despite a large number of improvements model deficiencies still remain. Examples include a too-fast vertical ascent and/or horizontal mixing in the tropical stratosphere, the cold temperature bias in the lowermost polar stratosphere, and the overestimation of polar total ozone loss during Antarctic springtime.

Description: Stable water isotopes in the coupled atmosphere–land surface model ECHAM5-JSBACH Geoscientific Model Development, 6, 1463-1480, 2013 Author(s): B. Haese, M. Werner, and G. Lohmann In this study we present first results of a new model development, ECHAM5-JSBACH-wiso, where we have incorporated the stable water isotopes H 2 18 O and HDO as tracers in the hydrological cycle of the coupled atmosphere–land surface model ECHAM5-JSBACH. The ECHAM5-JSBACH-wiso model was run under present-day climate conditions at two different resolutions (T31L19, T63L31). A comparison between ECHAM5-JSBACH-wiso and ECHAM5-wiso shows that the coupling has a strong impact on the simulated temperature and soil wetness. Caused by these changes of temperature and the hydrological cycle, the δ 18 O in precipitation also shows variations from −4‰ up to 4‰. One of the strongest anomalies is shown over northeast Asia where, due to an increase of temperature, the δ 18 O in precipitation increases as well. In order to analyze the sensitivity of the fractionation processes over land, we compare a set of simulations with various implementations of these processes over the land surface. The simulations allow us to distinguish between no fractionation, fractionation included in the evaporation flux (from bare soil) and also fractionation included in both evaporation and transpiration (from water transport through plants) fluxes. While the isotopic composition of the soil water may change for δ 18 O by up to +8&permil:, the simulated δ 18 O in precipitation shows only slight differences on the order of ±1‰. The simulated isotopic composition of precipitation fits well with the available observations from the GNIP (Global Network of Isotopes in Precipitation) database.

Description: δ 18 O water isotope in the i LOVECLIM model (version 1.0) – Part 2: Evaluation of model results against observed δ 18 O in water samples Geoscientific Model Development, 6, 1493-1504, 2013 Author(s): D. M. Roche and T. Caley The H 2 18 O stable isotope was previously introduced in the three coupled components of the earth system model i LOVECLIM: atmosphere, ocean and vegetation. The results of a long (5000 yr) pre-industrial equilibrium simulation are presented and evaluated against measurement of H 2 18 O abundance in present-day water for the atmospheric and oceanic components. For the atmosphere, it is found that the model reproduces the observed spatial distribution and relationships to climate variables with some merit, though limitations following our approach are highlighted. Indeed, we obtain the main gradients with a robust representation of the Rayleigh distillation but caveats appear in Antarctica and around the Mediterranean region due to model limitation. For the oceanic component, the agreement between the modelled and observed distribution of water δ 18 O is found to be very good. Mean ocean surface latitudinal gradients are faithfully reproduced as well as the mark of the main intermediate and deep water masses. This opens large prospects for the applications in palaeoclimatic context.

Description: The Rock Geochemical Model (RokGeM) v0.9 Geoscientific Model Development, 6, 1543-1573, 2013 Author(s): G. Colbourn, A. Ridgwell, and T. M. Lenton A new model of terrestrial rock weathering – the Rock Geochemical Model (RokGeM) – was developed for incorporation into the GENIE Earth System modelling framework. In this paper we describe the model. We consider a range of previously devised parameterizations, ranging from simple dependencies on global mean temperature following Berner et al. (1983), to spatially explicit dependencies on run-off and temperature (GKWM, Bluth and Kump, 1994; GEM-CO2, Amiotte-Suchet et al., 2003) – fields provided by the energy-moisture balance atmosphere model component in GENIE. Using long-term carbon cycle perturbation experiments, we test the effects of a wide range of model parameters, including whether or not the atmosphere was "short-circuited" in the carbon cycle; the sensitivity and feedback strength of temperature and run-off on carbonate and silicate weathering; different river-routing schemes; 0-D (global average) vs. 2-D (spatially explicit) weathering schemes; and the lithology dependence of weathering. Included are details of how to run the model and visualize the results.

Description: Distributed allocation of mobile sensing swarms in gyre flows Nonlinear Processes in Geophysics, 20, 657-668, 2013 Author(s): K. Mallory, M. A. Hsieh, E. Forgoston, and I. B. Schwartz We address the synthesis of distributed control policies to enable a swarm of homogeneous mobile sensors to maintain a desired spatial distribution in a geophysical flow environment, or workspace. In this article, we assume the mobile sensors (or robots) have a "map" of the environment denoting the locations of the Lagrangian coherent structures or LCS boundaries. Using this information, we design agent-level hybrid control policies that leverage the surrounding fluid dynamics and inherent environmental noise to enable the team to maintain a desired distribution in the workspace. We discuss the stability properties of the ensemble dynamics of the distributed control policies. Since realistic quasi-geostrophic ocean models predict double-gyre flow solutions, we use a wind-driven multi-gyre flow model to verify the feasibility of the proposed distributed control strategy and compare the proposed control strategy with a baseline deterministic allocation strategy. Lastly, we validate the control strategy using actual flow data obtained by our coherent structure experimental testbed.

Description: An approach to computing direction relations between separated object groups Geoscientific Model Development, 6, 1591-1599, 2013 Author(s): H. Yan, Z. Wang, and J. Li Direction relations between object groups play an important role in qualitative spatial reasoning, spatial computation and spatial recognition. However, none of existing models can be used to compute direction relations between object groups. To fill this gap, an approach to computing direction relations between separated object groups is proposed in this paper, which is theoretically based on gestalt principles and the idea of multi-directions. The approach firstly triangulates the two object groups, and then it constructs the Voronoi diagram between the two groups using the triangular network. After this, the normal of each Voronoi edge is calculated, and the quantitative expression of the direction relations is constructed. Finally, the quantitative direction relations are transformed into qualitative ones. The psychological experiments show that the proposed approach can obtain direction relations both between two single objects and between two object groups, and the results are correct from the point of view of spatial cognition.

Description: Four-dimensional ensemble-variational data assimilation for global deterministic weather prediction Nonlinear Processes in Geophysics, 20, 669-682, 2013 Author(s): M. Buehner, J. Morneau, and C. Charette The goal of this study is to evaluate a version of the ensemble-variational data assimilation approach (EnVar) for possible replacement of 4D-Var at Environment Canada for global deterministic weather prediction. This implementation of EnVar relies on 4-D ensemble covariances, obtained from an ensemble Kalman filter, that are combined in a vertically dependent weighted average with simple static covariances. Verification results are presented from a set of data assimilation experiments over two separate 6-week periods that used assimilated observations and model configuration very similar to the currently operational system. To help interpret the comparison of EnVar versus 4D-Var, additional experiments using 3D-Var and a version of EnVar with only 3-D ensemble covariances are also evaluated. To improve the rate of convergence for all approaches evaluated (including EnVar), an estimate of the cost function Hessian generated by the quasi-Newton minimization algorithm is cycled from one analysis to the next. Analyses from EnVar (with 4-D ensemble covariances) nearly always produce improved, and never degraded, forecasts when compared with 3D-Var. Comparisons with 4D-Var show that forecasts from EnVar analyses have either similar or better scores in the troposphere of the tropics and the winter extra-tropical region. However, in the summer extra-tropical region the medium-range forecasts from EnVar have either similar or worse scores than 4D-Var in the troposphere. In contrast, the 6 h forecasts from EnVar are significantly better than 4D-Var relative to radiosonde observations for both periods and in all regions. The use of 4-D versus 3-D ensemble covariances only results in small improvements in forecast quality. By contrast, the improvements from using 4D-Var versus 3D-Var are much larger. Measurement of the fit of the background and analyzed states to the observations suggests that EnVar and 4D-Var can both make better use of observations distributed over time than 3D-Var. In summary, the results from this study suggest that the EnVar approach is a viable alternative to 4D-Var, especially when the simplicity and computational efficiency of EnVar are considered. Additional research is required to understand the seasonal dependence of the difference in forecast quality between EnVar and 4D-Var in the extra-tropics.

Description: A method to represent ozone response to large changes in precursor emissions using high-order sensitivity analysis in photochemical models Geoscientific Model Development, 6, 1601-1608, 2013 Author(s): G. Yarwood, C. Emery, J. Jung, U. Nopmongcol, and T. Sakulyanontvittaya Photochemical grid models (PGMs) are used to simulate tropospheric ozone and quantify its response to emission changes. PGMs are often applied for annual simulations to provide both maximum concentrations for assessing compliance with air quality standards and frequency distributions for assessing human exposure. Efficient methods for computing ozone at different emission levels can improve the quality of ozone air quality management efforts. This study demonstrates the feasibility of using the decoupled direct method (DDM) to calculate first- and second-order sensitivity of ozone to anthropogenic NO x and VOC emissions in annual PGM simulations at continental scale. Algebraic models are developed that use Taylor series to produce complete annual frequency distributions of hourly ozone at any location and any anthropogenic emission level between zero and 100%, adjusted independently for NO x and VOC. We recommend computing the sensitivity coefficients at the midpoint of the emissions range over which they are intended to be applied, in this case with 50% anthropogenic emissions. The algebraic model predictions can be improved by combining sensitivity coefficients computed at 10 and 50% anthropogenic emissions. Compared to brute force simulations, algebraic model predictions tend to be more accurate in summer than winter, at rural than urban locations, and with 100% than zero anthropogenic emissions. Equations developed to combine sensitivity coefficients computed with 10 and 50% anthropogenic emissions are able to reproduce brute force simulation results with zero and 100% anthropogenic emissions with a mean bias of less than 2 ppb and mean error of less than 3 ppb averaged over 22 US cities.

Description: A Lagrangian advection scheme with shape matrix (LASM) for solving advection problems Geoscientific Model Development, 7, 2951-2968, 2014 Author(s): L. Dong, B. Wang, and L. Liu A new Lagrangian advection scheme with shape matrix (LASM) is proposed to take advantage of the extreme low numerical diffusion of the Lagrangian methods. The tracer is discretized into finite parcels, which move along the downstream trajectories. Different from other Lagrangian schemes, the parcel shape is simulated explicitly by a linear transformation matrix. By doing so, the aliasing error in the Lagrangian schemes is largely reduced without introducing substantial interparcel mixing in the pure advection stage, because the flow information will be respected when remapping tracer density onto the fixed model grids. An adaptive interparcel mixing algorithm is constructed to ensure the validity of the linear approximation of the parcel shape, where the mixing is only triggered when it is necessary and resembles the physical mixing. The total tracer mass on the parcels is conserved exactly. The new scheme is validated by using several test cases.

Description: Sensitivity of the Mediterranean sea level to atmospheric pressure and free surface elevation numerical formulation in NEMO Geoscientific Model Development, 7, 3001-3015, 2014 Author(s): P. Oddo, A. Bonaduce, N. Pinardi, and A. Guarnieri The sensitivity of the dynamics of the Mediterranean Sea to atmospheric pressure and free surface elevation formulation using NEMO (Nucleus for European Modelling of the Ocean) was evaluated. Four different experiments were carried out in the Mediterranean Sea using filtered or explicit free surface numerical schemes and accounting for the effect of atmospheric pressure in addition to wind and buoyancy fluxes. Model results were evaluated by coherency and power spectrum analysis with tide gauge data. We found that atmospheric pressure plays an important role for periods shorter than 100 days. The free surface formulation is important to obtain the correct ocean response for periods shorter than 30 days. At frequencies higher than 15 days −1 the Mediterranean basin's response to atmospheric pressure was not coherent and the performance of the model strongly depended on the specific area considered. A large-amplitude seasonal oscillation observed in the experiments using a filtered free surface was not evident in the corresponding explicit free surface formulation case, which was due to a phase shift between mass fluxes in the Gibraltar Strait and at the surface. The configuration with time splitting and atmospheric pressure always performed best; the differences were enhanced at very high frequencies.

Description: Gaseous chemistry and aerosol mechanism developments for version 3.5.1 of the online regional model, WRF-Chem Geoscientific Model Development, 7, 2557-2579, 2014 Author(s): S. Archer-Nicholls, D. Lowe, S. Utembe, J. Allan, R. A. Zaveri, J. D. Fast, Ø. Hodnebrog, H. Denier van der Gon, and G. McFiggans We have made a number of developments to the Weather, Research and Forecasting model coupled with Chemistry (WRF-Chem), with the aim of improving model prediction of trace atmospheric gas-phase chemical and aerosol composition, and of interactions between air quality and weather. A reduced form of the Common Reactive Intermediates gas-phase chemical mechanism (CRIv2-R5) has been added, using the Kinetic Pre-Processor (KPP) interface, to enable more explicit simulation of VOC degradation. N 2 O 5 heterogeneous chemistry has been added to the existing sectional MOSAIC aerosol module, and coupled to both the CRIv2-R5 and existing CBM-Z gas-phase schemes. Modifications have also been made to the sea-spray aerosol emission representation, allowing the inclusion of primary organic material in sea-spray aerosol. We have worked on the European domain, with a particular focus on making the model suitable for the study of nighttime chemistry and oxidation by the nitrate radical in the UK atmosphere. Driven by appropriate emissions, wind fields and chemical boundary conditions, implementation of the different developments are illustrated, using a modified version of WRF-Chem 3.4.1, in order to demonstrate the impact that these changes have in the Northwest European domain. These developments are publicly available in WRF-Chem from version 3.5.1 onwards.

Description: A strategy for GIS-based 3-D slope stability modelling over large areas Geoscientific Model Development, 7, 2969-2982, 2014 Author(s): M. Mergili, I. Marchesini, M. Alvioli, M. Metz, B. Schneider-Muntau, M. Rossi, and F. Guzzetti GIS-based deterministic models may be used for landslide susceptibility mapping over large areas. However, such efforts require specific strategies to (i) keep computing time at an acceptable level, and (ii) parameterize the geotechnical data. We test and optimize the performance of the GIS-based, 3-D slope stability model r.slope.stability in terms of computing time and model results. The model was developed as a C- and Python-based raster module of the open source software GRASS GIS and considers the 3-D geometry of the sliding surface. It calculates the factor of safety (FoS) and the probability of slope failure ( P f ) for a number of randomly selected potential slip surfaces, ellipsoidal or truncated in shape. Model input consists of a digital elevation model (DEM), ranges of geotechnical parameter values derived from laboratory tests, and a range of possible soil depths estimated in the field. Probability density functions are exploited to assign P f to each ellipsoid. The model calculates for each pixel multiple values of FoS and P f corresponding to different sliding surfaces. The minimum value of FoS and the maximum value of P f for each pixel give an estimate of the landslide susceptibility in the study area. Optionally, r.slope.stability is able to split the study area into a defined number of tiles, allowing parallel processing of the model on the given area. Focusing on shallow landslides, we show how multi-core processing makes it possible to reduce computing times by a factor larger than 20 in the study area. We further demonstrate how the number of random slip surfaces and the sampling of parameters influence the average value of P f and the capacity of r.slope.stability to predict the observed patterns of shallow landslides in the 89.5 km 2 Collazzone area in Umbria, central Italy.

Description: Corrigendum to "Breeding and predictability in the baroclinic rotating annulus using a perfect model" published in Nonlin. Processes Geophys., 15, 469–487, 2008 Nonlinear Processes in Geophysics, 18, 359-359, 2011 Author(s): R. M. B. Young and P. L. Read No abstract available.

Description: Strategies for coupling global and limited-area ensemble Kalman filter assimilation Nonlinear Processes in Geophysics, 18, 415-430, 2011 Author(s): D. Merkova, I. Szunyogh, and E. Ott This paper compares the forecast performance of four strategies for coupling global and limited area data assimilation: three strategies propagate information from the global to the limited area process, while the fourth strategy feeds back information from the limited area to the global process. All four strategies are formulated in the Local Ensemble Transform Kalman Filter (LETKF) framework. Numerical experiments are carried out with the model component of the National Centers for Environmental Prediction (NCEP) Global Forecast System (GFS) and the NCEP Regional Spectral Model (RSM). The limited area domain is an extended North-America region that includes part of the north-east Pacific. The GFS is integrated at horizontal resolution T62 (about 150 km in the mid-latitudes), while the RSM is integrated at horizontal resolution 48 km. Experiments are carried out both under the perfect model hypothesis and in a realistic setting. The coupling strategies are evaluated by comparing their deterministic forecast performance at 12-h and 48-h lead times. The results suggest that the limited area data assimilation system has the potential to enhance the forecasts at 12-h lead time in the limited area domain at the synoptic and sub-synoptic scales (in the global wave number range of about 10 to 40). There is a clear indication that between the forecast performance of the different coupling strategies those that cycle the limited area assimilation process produce the most accurate forecasts. In the realistic setting, at 12-h forecast time the limited area systems produce more modest improvements compared to the global system than under the perfect model hypothesis, and at 48-h forecast time the global forecasts are more accurate than the limited area forecasts.

Description: Spectral methods for internal waves: indistinguishable density profiles and double-humped solitary waves Nonlinear Processes in Geophysics, 18, 351-358, 2011 Author(s): M. Dunphy, C. Subich, and M. Stastna Internal solitary waves are widely observed in both the oceans and large lakes. They can be described by a variety of mathematical theories, covering the full spectrum from first order asymptotic theory (i.e. Korteweg-de Vries, or KdV, theory), through higher order extensions of weakly nonlinear-weakly nonhydrostatic theory, to fully nonlinear-weakly nonhydrostatic theories and finally exact theory based on the Dubreil-Jacotin-Long (DJL) equation that is formally equivalent to the full set of Euler equations. We discuss how spectral and pseudospectral methods allow for the computation of novel phenomena in both approximate and exact theories. In particular we construct markedly different density profiles for which the coefficients in the KdV theory are very nearly identical. These two density profiles yield qualitatively different behaviour for both exact, or fully nonlinear, waves computed using the DJL equation and in dynamic simulations of the time dependent Euler equations. For exact, DJL, theory we compute exact solitary waves with two-scales, or so-called double-humped waves.

Description: A nonlinear method of removing harmonic noise in geophysical data Nonlinear Processes in Geophysics, 18, 367-379, 2011 Author(s): Y. Jeng and C.-S. Chen A nonlinear, adaptive method to remove the harmonic noise that commonly resides in geophysical data is proposed in this study. This filtering method is based on the ensemble empirical mode decomposition algorithm in conjunction with the logarithmic transform. We present a synthetic model study to investigate the capability of signal reconstruction from the decomposed data, and compare the results with those derived from other 2-D adaptive filters. Applications to the real seismic data acquired by using an ocean bottom seismograph and to a shot gather of the ground penetrating radar demonstrate the robustness of this method. Our work proposes a concept that instead of Fourier-based approaches, the harmonic noise removal in geophysical data can be achieved effectively by using an alternative nonlinear adaptive data analysis method, which has been applied extensively in other scientific studies.

Description: Comparison of correlation analysis techniques for irregularly sampled time series Nonlinear Processes in Geophysics, 18, 389-404, 2011 Author(s): K. Rehfeld, N. Marwan, J. Heitzig, and J. Kurths Geoscientific measurements often provide time series with irregular time sampling, requiring either data reconstruction (interpolation) or sophisticated methods to handle irregular sampling. We compare the linear interpolation technique and different approaches for analyzing the correlation functions and persistence of irregularly sampled time series, as Lomb-Scargle Fourier transformation and kernel-based methods. In a thorough benchmark test we investigate the performance of these techniques. All methods have comparable root mean square errors (RMSEs) for low skewness of the inter-observation time distribution. For high skewness, very irregular data, interpolation bias and RMSE increase strongly. We find a 40 % lower RMSE for the lag-1 autocorrelation function (ACF) for the Gaussian kernel method vs. the linear interpolation scheme,in the analysis of highly irregular time series. For the cross correlation function (CCF) the RMSE is then lower by 60 %. The application of the Lomb-Scargle technique gave results comparable to the kernel methods for the univariate, but poorer results in the bivariate case. Especially the high-frequency components of the signal, where classical methods show a strong bias in ACF and CCF magnitude, are preserved when using the kernel methods. We illustrate the performances of interpolation vs. Gaussian kernel method by applying both to paleo-data from four locations, reflecting late Holocene Asian monsoon variability as derived from speleothem δ 18 O measurements. Cross correlation results are similar for both methods, which we attribute to the long time scales of the common variability. The persistence time (memory) is strongly overestimated when using the standard, interpolation-based, approach. Hence, the Gaussian kernel is a reliable and more robust estimator with significant advantages compared to other techniques and suitable for large scale application to paleo-data.

Description: Semi-Lagrangian methods in air pollution models Geoscientific Model Development, 4, 511-541, 2011 Author(s): A. B. Hansen, J. Brandt, J. H. Christensen, and E. Kaas Various semi-Lagrangian methods are tested with respect to advection in air pollution modeling. The aim is to find a method fulfilling as many of the desirable properties by Rasch andWilliamson (1990) and Machenhauer et al. (2008) as possible. The focus in this study is on accuracy and local mass conservation. The methods tested are, first, classical semi-Lagrangian cubic interpolation, see e.g. Durran (1999), second, semi-Lagrangian cubic cascade interpolation, by Nair et al. (2002), third, semi-Lagrangian cubic interpolation with the modified interpolation weights, Locally Mass Conserving Semi-Lagrangian (LMCSL), by Kaas (2008), and last, semi-Lagrangian cubic interpolation with a locally mass conserving monotonic filter by Kaas and Nielsen (2010). Semi-Lagrangian (SL) interpolation is a classical method for atmospheric modeling, cascade interpolation is more efficient computationally, modified interpolation weights assure mass conservation and the locally mass conserving monotonic filter imposes monotonicity. All schemes are tested with advection alone or with advection and chemistry together under both typical rural and urban conditions using different temporal and spatial resolution. The methods are compared with a current state-of-the-art scheme, Accurate Space Derivatives (ASD), see Frohn et al. (2002), presently used at the National Environmental Research Institute (NERI) in Denmark. To enable a consistent comparison only non-divergent flow configurations are tested. The test cases are based either on the traditional slotted cylinder or the rotating cone, where the schemes' ability to model both steep gradients and slopes are challenged. The tests showed that the locally mass conserving monotonic filter improved the results significantly for some of the test cases, however, not for all. It was found that the semi-Lagrangian schemes, in almost every case, were not able to outperform the current ASD scheme used in DEHM with respect to accuracy.

Description: GEOCLIM reloaded (v 1.0): a new coupled earth system model for past climate change Geoscientific Model Development, 4, 451-481, 2011 Author(s): S. Arndt, P. Regnier, Y. Goddéris, and Y. Donnadieu We present a new version of the coupled Earth system model GEOCLIM. The new release, GEOCLIM reloaded (v 1.0), links the existing atmosphere and weathering modules to a novel, temporally and spatially resolved model of the global ocean circulation, which provides a physical framework for a mechanistic description of the marine biogeochemical dynamics of carbon, nitrogen, phosphorus and oxygen. The ocean model is also coupled to a fully formulated, vertically resolved diagenetic model. GEOCLIM reloaded is thus a unique tool to investigate the short- and long-term feedbacks between climatic conditions, continental inputs, ocean biogeochemical dynamics and diagenesis. A complete and detailed description of the resulting Earth system model and its new features is first provided. The performance of GEOCLIM reloaded is then evaluated by comparing steady-state simulation under present-day conditions with a comprehensive set of oceanic data and existing global estimates of bio-element cycling in the pelagic and benthic compartments.

Description: Brief communication "On one mechanism of low frequency variability of the Antarctic Circumpolar Current" Nonlinear Processes in Geophysics, 18, 361-365, 2011 Author(s): O. G. Derzho and B. de Young In this paper we present a simple analytical model for low frequency and large scale variability of the Antarctic Circumpolar Current (ACC). The physical mechanism of the variability is related to temporal and spatial variations of the cyclonic mean flow (ACC) due to circularly propagating nonlinear barotropic Rossby wave trains. It is shown that the Rossby wave train is a fundamental mode, trapped between the major fronts in the ACC. The Rossby waves are predicted to rotate with a particular angular velocity that depends on the magnitude and width of the mean current. The spatial structure of the rotating pattern, including its zonal wave number, is defined by the specific form of the stream function-vorticity relation. The similarity between the simulated patterns and the Antarctic Circumpolar Wave (ACW) is highlighted. The model can predict the observed sequence of warm and cold patches in the ACW as well as its zonal number.

Description: The CSIRO Mk3L climate system model version 1.0 – Part 1: Description and evaluation Geoscientific Model Development, 4, 483-509, 2011 Author(s): S. J. Phipps, L. D. Rotstayn, H. B. Gordon, J. L. Roberts, A. C. Hirst, and W. F. Budd The CSIRO Mk3L climate system model is a coupled general circulation model, designed primarily for millennial-scale climate simulations and palaeoclimate research. Mk3L includes components which describe the atmosphere, ocean, sea ice and land surface, and combines computational efficiency with a stable and realistic control climatology. This paper describes the model physics and software, analyses the control climatology, and evaluates the ability of the model to simulate the modern climate. Mk3L incorporates a spectral atmospheric general circulation model, a z -coordinate ocean general circulation model, a dynamic-thermodynamic sea ice model and a land surface scheme with static vegetation. The source code is highly portable, and has no dependence upon proprietary software. The model distribution is freely available to the research community. A 1000-yr climate simulation can be completed in around one-and-a-half months on a typical desktop computer, with greater throughput being possible on high-performance computing facilities. Mk3L produces realistic simulations of the larger-scale features of the modern climate, although with some biases on the regional scale. The model also produces reasonable representations of the leading modes of internal climate variability in both the tropics and extratropics. The control state of the model exhibits a high degree of stability, with only a weak cooling trend on millennial timescales. Ongoing development work aims to improve the model climatology and transform Mk3L into a comprehensive earth system model.

Description: A semi-phenomenological approach to explain the event-size distribution of the Drossel-Schwabl forest-fire model Nonlinear Processes in Geophysics, 18, 381-388, 2011 Author(s): S. Hergarten and R. Krenn We present a novel approach to explain the complex scaling behavior of the Drossel-Schwabl forest-fire model in two dimensions. Clusters of trees are characterized by their size and perimeter only, whereas spatial correlations are neglected. Coalescence of clusters is restricted to clusters of similar sizes. Our approach derives the value of the scaling exponent τ of the event size distribution directly from the scaling of the accessible perimeter of percolation clusters. We obtain τ = 1.19 in the limit of infinite growth rate, in perfect agreement with numerical results. Furthermore, our approach predicts the unusual transition from a power law to an exponential decay even quantitatively, while the exponential decay at large event sizes itself is reproduced only qualitatively.

Description: Double rank-ordering technique of ROMA (Rank-Ordered Multifractal Analysis) for multifractal fluctuations featuring multiple regimes of scales Nonlinear Processes in Geophysics, 18, 405-414, 2011 Author(s): S. W. Y. Tam and T. Chang Rank-Ordered Multifractal Analysis (ROMA), a technique capable of deciphering the multifractal characteristics of intermittent fluctuations, was originally applied to the results of a magnetohydrodynamic (MHD) simulation. Application of ROMA to measured fluctuations in the auroral zone, due to the dominant physical effects changing from kinetic to MHD as the scale increases, requires an additional level of rank-ordering in order to divide the domain of scales into regimes. An algorithm for the additional step in this double rank-ordering technique is discussed, and is demonstrated in the application to the electric field fluctuations in the auroral zone as an example. As a result of the double rank-ordering, ROMA is able to take into account the nonlinear crossover behavior characterized by the multiple regimes of time scales by providing a scaling variable and a scaling function that are global to all the time scales.

Description: Medusa-1.0: a new intermediate complexity plankton ecosystem model for the global domain Geoscientific Model Development, 4, 381-417, 2011 Author(s): A. Yool, E. E. Popova, and T. R. Anderson The ongoing, anthropogenically-driven changes to the global ocean are expected to have significant consequences for plankton ecosystems in the future. Because of the role that plankton play in the ocean's "biological pump", changes in abundance, distribution and productivity will likely have additional consequences for the wider carbon cycle. Just as in the terrestrial biosphere, marine ecosystems exhibit marked diversity in species and functional types of organisms. Predicting potential change in plankton ecosystems therefore requires the use of models that are suited to this diversity, but whose parameterisation also permits robust and realistic functional behaviour. In the past decade, advances in model sophistication have attempted to address diversity, but have been criticised for doing so inaccurately or ahead of a requisite understanding of underlying processes. Here we introduce MEDUSA-1.0 ( M odel of E cosystem D ynamics, nutrient U tilisation, S equestration and A cidification), a new "intermediate complexity" plankton ecosystem model that expands on traditional nutrient-phytoplankton-zooplankton-detritus (NPZD) models, and remains amenable to global-scale evaluation. MEDUSA-1.0 includes the biogeochemical cycles of nitrogen, silicon and iron, broadly structured into "small" and "large" plankton size classes, of which the "large" phytoplankton class is representative of a key phytoplankton group, the diatoms. A full description of MEDUSA-1.0's state variables, differential equations, functional forms and parameter values is included, with particular attention focused on the submodel describing the export of organic carbon from the surface to the deep ocean. MEDUSA-1.0 is used here in a multi-decadal hindcast simulation, and its biogeochemical performance evaluated at the global scale.

Description: PREP-CHEM-SRC – 1.0: a preprocessor of trace gas and aerosol emission fields for regional and global atmospheric chemistry models Geoscientific Model Development, 4, 419-433, 2011 Author(s): S. R. Freitas, K. M. Longo, M. F. Alonso, M. Pirre, V. Marecal, G. Grell, R. Stockler, R. F. Mello, and M. Sánchez Gácita The preprocessor PREP-CHEM-SRC presented in the paper is a comprehensive tool aiming at preparing emission fields of trace gases and aerosols for use in atmospheric-chemistry transport models. The considered emissions are from the most recent databases of urban/industrial, biogenic, biomass burning, volcanic, biofuel use and burning from agricultural waste sources. For biomass burning, emissions can be also estimated directly from satellite fire detections using a fire emission model included in the tool. The preprocessor provides emission fields interpolated onto the transport model grid. Several map projections can be chosen. The inclusion of these emissions in transport models is also presented. The preprocessor is coded using Fortran90 and C and is driven by a namelist allowing the user to choose the type of emissions and the databases.

Description: Large eddy simulation of sediment transport over rippled beds Nonlinear Processes in Geophysics, 21, 1169-1184, 2014 Author(s): J. C. Harris and S. T. Grilli Wave-induced boundary layer (BL) flows over sandy rippled bottoms are studied using a numerical model that applies a one-way coupling of a "far-field" inviscid flow model to a "near-field" large eddy simulation (LES) Navier–Stokes (NS) model. The incident inviscid velocity and pressure fields force the LES, in which near-field, wave-induced, turbulent bottom BL flows are simulated. A sediment suspension and transport model is embedded within the coupled flow model. The numerical implementation of the various models has been reported elsewhere, where we showed that the LES was able to accurately simulate both mean flow and turbulent statistics for oscillatory BL flows over a flat, rough bed. Here we show that the model accurately predicts the mean velocity fields and suspended sediment concentration for oscillatory flows over full-scale vortex ripples. Tests show that surface roughness has a significant effect on the results. Beyond increasing our insight into wave-induced oscillatory bottom BL physics, sophisticated coupled models of sediment transport such as that presented have the potential to make quantitative predictions of sediment transport and erosion/accretion around partly buried objects in the bottom, which is important for a vast array of bottom deployed instrumentation and other practical ocean engineering problems.

Description: Tropical troposphere to stratosphere transport of carbon monoxide and long-lived trace species in the Chemical Lagrangian Model of the Stratosphere (CLaMS) Geoscientific Model Development, 7, 2895-2916, 2014 Author(s): R. Pommrich, R. Müller, J.-U. Grooß, P. Konopka, F. Ploeger, B. Vogel, M. Tao, C. M. Hoppe, G. Günther, N. Spelten, L. Hoffmann, H.-C. Pumphrey, S. Viciani, F. D'Amato, C. M. Volk, P. Hoor, H. Schlager, and M. Riese Variations in the mixing ratio of trace gases of tropospheric origin entering the stratosphere in the tropics are of interest for assessing both troposphere to stratosphere transport fluxes in the tropics and the impact of these transport fluxes on the composition of the tropical lower stratosphere. Anomaly patterns of carbon monoxide (CO) and long-lived tracers in the lower tropical stratosphere allow conclusions about the rate and the variability of tropical upwelling to be drawn. Here, we present a simplified chemistry scheme for the Chemical Lagrangian Model of the Stratosphere (CLaMS) for the simulation, at comparatively low numerical cost, of CO, ozone, and long-lived trace substances (CH 4 , N 2 O, CCl 3 F (CFC-11), CCl 2 F 2 (CFC-12), and CO 2 ) in the lower tropical stratosphere. For the long-lived trace substances, the boundary conditions at the surface are prescribed based on ground-based measurements in the lowest model level. The boundary condition for CO in the lower troposphere (below about 4 km) is deduced from MOPITT measurements. Due to the lack of a specific representation of mixing and convective uplift in the troposphere in this model version, enhanced CO values, in particular those resulting from convective outflow are underestimated. However, in the tropical tropopause layer and the lower tropical stratosphere, there is relatively good agreement of simulated CO with in situ measurements (with the exception of the TROCCINOX campaign, where CO in the simulation is biased low ≈10–15 ppbv). Further, the model results (and therefore also the ERA-Interim winds, on which the transport in the model is based) are of sufficient quality to describe large scale anomaly patterns of CO in the lower stratosphere. In particular, the zonally averaged tropical CO anomaly patterns (the so called "tape recorder" patterns) simulated by this model version of CLaMS are in good agreement with observations, although the simulations show a too rapid upwelling compared to observations as a consequence of the overestimated vertical velocities in the ERA-Interim reanalysis data set. Moreover, the simulated tropical anomaly patterns of N 2 O are in good agreement with observations. In the simulations, anomaly patterns of CH 4 and CFC-11 were found to be very similar to those of N 2 O; for all long-lived tracers, positive anomalies are simulated because of the enhanced tropical upwelling in the easterly shear phase of the quasi-biennial oscillation.

Description: Dependence of sandpile avalanche frequency–size distribution on coverage extent and compactness of embedded toppling threshold heterogeneity: implications for the variation of Gutenberg–Richter b value Nonlinear Processes in Geophysics, 21, 1185-1193, 2014 Author(s): L.-Y. Chiao and Q. Liu The effects of the spatiotemporal evolution of failure threshold heterogeneity on the dynamics of fault criticality, and thus on regional seismogenesis, have attracted strong interest in the field of regional seismotectonics. The heterogeneity might be a manifestation of the macroscopic distribution and multiscale strength variation of asperities, the distinct regional stress level, and (microscopically) heterogeneous fault surface roughness or friction regimes. In this study, rather than attempting to mimic the complex microscale slipping physics on a fault surface, sandpile cellular automata were implemented with a straightforward toppling rule. The objective is to examine the influence of distinct configurations of the embedded heterogeneous toppling threshold field on the global system avalanche event statistics. The examination results revealed that increasing the coverage extent and decreasing the compactness of the heterogeneous failure threshold, rather than the magnitude, range of contrast, diversity, or the geometric configuration of the threshold heterogeneity, leads to a systematic increase in the scaling exponent of the avalanche event power law statistics, implying the importance of mutual interaction among toppling sites with distinct thresholds. For tectonic provinces with differing stress regimes evolving spatio temporally, it is postulated that the distinct extent and compactness of the heterogeneous failure threshold are critical factors that manifest in the reported dynamic variations of seismicity scaling.

Description: Instability and change detection in exponential families and generalized linear models, with a study of Atlantic tropical storms Nonlinear Processes in Geophysics, 21, 1133-1143, 2014 Author(s): Y. Lu and S. Chatterjee Exponential family statistical distributions, including the well-known normal, binomial, Poisson, and exponential distributions, are overwhelmingly used in data analysis. In the presence of covariates, an exponential family distributional assumption for the response random variables results in a generalized linear model. However, it is rarely ensured that the parameters of the assumed distributions are stable through the entire duration of the data collection process. A failure of stability leads to nonsmoothness and nonlinearity in the physical processes that result in the data. In this paper, we propose testing for stability of parameters of exponential family distributions and generalized linear models. A rejection of the hypothesis of stable parameters leads to change detection. We derive the related likelihood ratio test statistic. We compare the performance of this test statistic to the popular normal distributional assumption dependent cumulative sum (Gaussian CUSUM) statistic in change detection problems. We study Atlantic tropical storms using the techniques developed here, so to understand whether the nature of these tropical storms has remained stable over the last few decades.

Description: Evaluation of the ECHAM family radiation codes performance in the representation of the solar signal Geoscientific Model Development, 7, 2859-2866, 2014 Author(s): T. Sukhodolov, E. Rozanov, A. I. Shapiro, J. Anet, C. Cagnazzo, T. Peter, and W. Schmutz Solar radiation is the main source of energy for the Earth's atmosphere and in many respects defines its composition, photochemistry, temperature profile and dynamics. The magnitude of the solar irradiance variability strongly depends on the wavelength, making difficult its representation in climate models. Due to some deficiencies in the applied radiation codes, several models fail to show a clear response in middle stratospheric heating rates to solar spectral irradiance variability; therefore, it is important to evaluate model performance in this respect before doing multiple runs. In this work we evaluate the performance of three generations of ECHAM (4, 5 and 6) solar radiation schemes by a comparison with the reference high-resolution libRadtran code. We found that all original ECHAM radiation codes miss almost all solar signals in the heating rates in the mesosphere. In the stratosphere the two-band ECHAM4 code (E4) has an almost negligible radiative response to solar irradiance changes and the six-band ECHAM5 code (E5c) reproduces only about half of the reference signal, while representation in the ECHAM6 code (E6) is better – it misses a maximum of about 15% in the upper stratosphere. On the basis of the comparison results we suggest necessary improvements to the ECHAM family codes by the inclusion of available parameterizations of the heating rate due to absorption by oxygen (O 2 ) and ozone (O 3 ). Improvement is presented for E5c and E6, and both codes, with the introduced parameterizations, represent the heating rate response to the spectral solar irradiance variability simulated with libRadtran much better without a substantial increase in computer time. The suggested parameterizations are recommended to be applied in the middle-atmosphere version of the ECHAM-5 and 6 models for the study of the solar irradiance influence on climate.

Description: GEOtop 2.0: simulating the combined energy and water balance at and below the land surface accounting for soil freezing, snow cover and terrain effects Geoscientific Model Development, 7, 2831-2857, 2014 Author(s): S. Endrizzi, S. Gruber, M. Dall'Amico, and R. Rigon GEOtop is a fine-scale grid-based simulator that represents the heat and water budgets at and below the soil surface. It describes the three-dimensional water flow in the soil and the energy exchange with the atmosphere, considering the radiative and turbulent fluxes. Furthermore, it reproduces the highly non-linear interactions between the water and energy balance during soil freezing and thawing, and simulates the temporal evolution of the water and energy budgets in the snow cover and their effect on soil temperature. Here, we present the core components of GEOtop 2.0 and demonstrate its functioning. Based on a synthetic simulation, we show that the interaction of processes represented in GEOtop 2.0 can result in phenomena that are significant and relevant for applications involving permafrost and seasonally frozen soils, both in high altitude and latitude regions.

Description: The North American Carbon Program Multi-scale Synthesis and Terrestrial Model Intercomparison Project – Part 2: Environmental driver data Geoscientific Model Development, 7, 2875-2893, 2014 Author(s): Y. Wei, S. Liu, D. N. Huntzinger, A. M. Michalak, N. Viovy, W. M. Post, C. R. Schwalm, K. Schaefer, A. R. Jacobson, C. Lu, H. Tian, D. M. Ricciuto, R. B. Cook, J. Mao, and X. Shi Ecosystems are important and dynamic components of the global carbon cycle, and terrestrial biospheric models (TBMs) are crucial tools in further understanding of how terrestrial carbon is stored and exchanged with the atmosphere across a variety of spatial and temporal scales. Improving TBM skills, and quantifying and reducing their estimation uncertainties, pose significant challenges. The Multi-scale Synthesis and Terrestrial Model Intercomparison Project (MsTMIP) is a formal multi-scale and multi-model intercomparison effort set up to tackle these challenges. The MsTMIP protocol prescribes standardized environmental driver data that are shared among model teams to facilitate model–model and model–observation comparisons. This paper describes the global and North American environmental driver data sets prepared for the MsTMIP activity to both support their use in MsTMIP and make these data, along with the processes used in selecting/processing these data, accessible to a broader audience. Based on project needs and lessons learned from past model intercomparison activities, we compiled climate, atmospheric CO 2 concentrations, nitrogen deposition, land use and land cover change (LULCC), C3 / C4 grasses fractions, major crops, phenology and soil data into a standard format for global (0.5° × 0.5° resolution) and regional (North American: 0.25° × 0.25° resolution) simulations. In order to meet the needs of MsTMIP, improvements were made to several of the original environmental data sets, by improving the quality, and/or changing their spatial and temporal coverage, and resolution. The resulting standardized model driver data sets are being used by over 20 different models participating in MsTMIP. The data are archived at the Oak Ridge National Laboratory Distributed Active Archive Center (ORNL DAAC, http://daac.ornl.gov ) to provide long-term data management and distribution.

Description: The atmospheric chemistry box model CAABA/MECCA-3.0 Geoscientific Model Development, 4, 373-380, 2011 Author(s): R. Sander, A. Baumgaertner, S. Gromov, H. Harder, P. Jöckel, A. Kerkweg, D. Kubistin, E. Regelin, H. Riede, A. Sandu, D. Taraborrelli, H. Tost, and Z.-Q. Xie We present version 3.0 of the atmospheric chemistry box model CAABA/MECCA. In addition to a complete update of the rate coefficients to the most recent recommendations, a number of new features have been added: chemistry in multiple aerosol size bins; automatic multiple simulations reaching steady-state conditions; Monte-Carlo simulations with randomly varied rate coefficients within their experimental uncertainties; calculations along Lagrangian trajectories; mercury chemistry; more detailed isoprene chemistry; tagging of isotopically labeled species. Further changes have been implemented to make the code more user-friendly and to facilitate the analysis of the model results. Like earlier versions, CAABA/MECCA-3.0 is a community model published under the GNU General Public License.

Description: Rank ordering multifractal analysis of the auroral electrojet index Nonlinear Processes in Geophysics, 18, 277-285, 2011 Author(s): G. Consolini and P. De Michelis In the second half of the 90s interest grew on the complex features of the magnetospheric dynamics in response to solar wind changes. An important series of papers were published on the occurrence of chaos, turbulence and complexity. Among them, particularly interesting was the study of the bursty and fractal/multifractal character of the high latitude energy release during geomagnetic storms, which was evidenced by analyzing the features of the Auroral Electrojet (AE) indices. Recently, the multifractal features of the small time-scale increments of AE-indices have been criticized in favor of a more simple fractal behavior. This is particularly true for the scaling features of the probability density functions (PDFs) of the AE index increments. Here, after a brief review of the nature of the fractal/multifractal features of the magnetospheric response to solar wind changes, we investigate the multifractal nature of the scaling features of the AE index increments PDFs using the Rank Ordering Multifractal Analysis (ROMA) technique. The ROMA results clearly demonstrate the existence of a hierarchy of scaling indices, depending on the increment amplitude, for the data collapsing of PDFs relative to increments at different time scales. Our results confirm the previous results by Consolini et al. (1996) and the more recent results by Rypdal and Rypdal (2010).

Description: A comprehensive set of benchmark tests for a land surface model of simultaneous fluxes of water and carbon at both the global and seasonal scale Geoscientific Model Development, 4, 255-269, 2011 Author(s): E. Blyth, D. B. Clark, R. Ellis, C. Huntingford, S. Los, M. Pryor, M. Best, and S. Sitch Evaluating the models we use in prediction is important as it allows us to identify uncertainties in prediction as well as guiding the priorities for model development. This paper describes a set of benchmark tests that is designed to quantify the performance of the land surface model that is used in the UK Hadley Centre General Circulation Model (JULES: Joint UK Land Environment Simulator). The tests are designed to assess the ability of the model to reproduce the observed fluxes of water and carbon at the global and regional spatial scale, and on a seasonal basis. Five datasets are used to test the model: water and carbon dioxide fluxes from ten FLUXNET sites covering the major global biomes, atmospheric carbon dioxide concentrations at four representative stations from the global network, river flow from seven catchments, the seasonal mean NDVI over the seven catchments and the potential land cover of the globe (after the estimated anthropogenic changes have been removed). The model is run in various configurations and results are compared with the data. A few examples are chosen to demonstrate the importance of using combined use of observations of carbon and water fluxes in essential in order to understand the causes of model errors. The benchmarking approach is suitable for application to other global models.

Description: A new version of the CNRM Chemistry-Climate Model, CNRM-CCM: description and improvements from the CCMVal-2 simulations Geoscientific Model Development, 4, 873-900, 2011 Author(s): M. Michou, D. Saint-Martin, H. Teyssèdre, A. Alias, F. Karcher, D. Olivié, A. Voldoire, B. Josse, V.-H. Peuch, H. Clark, J. N. Lee, and F. Chéroux This paper presents a new version of the Météo-France CNRM Chemistry-Climate Model, so-called CNRM-CCM. It includes some fundamental changes from the previous version (CNRM-ACM) which was extensively evaluated in the context of the CCMVal-2 validation activity. The most notable changes concern the radiative code of the GCM, and the inclusion of the detailed stratospheric chemistry of our Chemistry-Transport model MOCAGE on-line within the GCM. A 47-yr transient simulation (1960–2006) is the basis of our analysis. CNRM-CCM generates satisfactory dynamical and chemical fields in the stratosphere. Several shortcomings of CNRM-ACM simulations for CCMVal-2 that resulted from an erroneous representation of the impact of volcanic aerosols as well as from transport deficiencies have been eliminated. Remaining problems concern the upper stratosphere (5 to 1 hPa) where temperatures are too high, and where there are biases in the NO 2 , N 2 O 5 and O 3 mixing ratios. In contrast, temperatures at the tropical tropopause are too cold. These issues are addressed through the implementation of a more accurate radiation scheme at short wavelengths. Despite these problems we show that this new CNRM CCM is a useful tool to study chemistry-climate applications.

Description: MIROC-ESM 2010: model description and basic results of CMIP5-20c3m experiments Geoscientific Model Development, 4, 845-872, 2011 Author(s): S. Watanabe, T. Hajima, K. Sudo, T. Nagashima, T. Takemura, H. Okajima, T. Nozawa, H. Kawase, M. Abe, T. Yokohata, T. Ise, H. Sato, E. Kato, K. Takata, S. Emori, and M. Kawamiya An earth system model (MIROC-ESM 2010) is fully described in terms of each model component and their interactions. Results for the CMIP5 (Coupled Model Intercomparison Project phase 5) historical simulation are presented to demonstrate the model's performance from several perspectives: atmosphere, ocean, sea-ice, land-surface, ocean and terrestrial biogeochemistry, and atmospheric chemistry and aerosols. An atmospheric chemistry coupled version of MIROC-ESM (MIROC-ESM-CHEM 2010) reasonably reproduces transient variations in surface air temperatures for the period 1850–2005, as well as the present-day climatology for the zonal-mean zonal winds and temperatures from the surface to the mesosphere. The historical evolution and global distribution of column ozone and the amount of tropospheric aerosols are reasonably simulated in the model based on the Representative Concentration Pathways' (RCP) historical emissions of these precursors. The simulated distributions of the terrestrial and marine biogeochemistry parameters agree with recent observations, which is encouraging to use the model for future global change projections.

Description: Preface Large amplitude internal waves in the coastal ocean Nonlinear Processes in Geophysics, 18, 653-655, 2011 Author(s): R. Grimshaw, K. Helfrich, and A. Scotti No Abstract available.

Description: Size distribution and structure of Barchan dune fields Nonlinear Processes in Geophysics, 18, 455-467, 2011 Author(s): O. Durán, V. Schwämmle, P. G. Lind, and H. J. Herrmann Barchans are isolated mobile dunes often organized in large dune fields. Dune fields seem to present a characteristic dune size and spacing, which suggests a cooperative behavior based on dune interaction. In Duran et al. (2009), we propose that the redistribution of sand by collisions between dunes is a key element for the stability and size selection of barchan dune fields. This approach was based on a mean-field model ignoring the spatial distribution of dune fields. Here, we present a simplified dune field model that includes the spatial evolution of individual dunes as well as their interaction through sand exchange and binary collisions. As a result, the dune field evolves towards a steady state that depends on the boundary conditions. Comparing our results with measurements of Moroccan dune fields, we find that the simulated fields have the same dune size distribution as in real fields but fail to reproduce their homogeneity along the wind direction.

Description: Combining 2-m temperature nowcasting and short range ensemble forecasting Nonlinear Processes in Geophysics, 18, 903-910, 2011 Author(s): A. Kann, T. Haiden, and C. Wittmann During recent years, numerical ensemble prediction systems have become an important tool for estimating the uncertainties of dynamical and physical processes as represented in numerical weather models. The latest generation of limited area ensemble prediction systems (LAM-EPSs) allows for probabilistic forecasts at high resolution in both space and time. However, these systems still suffer from systematic deficiencies. Especially for nowcasting (0–6 h) applications the ensemble spread is smaller than the actual forecast error. This paper tries to generate probabilistic short range 2-m temperature forecasts by combining a state-of-the-art nowcasting method and a limited area ensemble system, and compares the results with statistical methods. The Integrated Nowcasting Through Comprehensive Analysis (INCA) system, which has been in operation at the Central Institute for Meteorology and Geodynamics (ZAMG) since 2006 (Haiden et al., 2011), provides short range deterministic forecasts at high temporal (15 min–60 min) and spatial (1 km) resolution. An INCA Ensemble (INCA-EPS) of 2-m temperature forecasts is constructed by applying a dynamical approach, a statistical approach, and a combined dynamic-statistical method. The dynamical method takes uncertainty information (i.e. ensemble variance) from the operational limited area ensemble system ALADIN-LAEF (Aire Limitée Adaptation Dynamique Développement InterNational Limited Area Ensemble Forecasting) which is running operationally at ZAMG (Wang et al., 2011). The purely statistical method assumes a well-calibrated spread-skill relation and applies ensemble spread according to the skill of the INCA forecast of the most recent past. The combined dynamic-statistical approach adapts the ensemble variance gained from ALADIN-LAEF with non-homogeneous Gaussian regression (NGR) which yields a statistical \mbox{correction} of the first and second moment (mean bias and dispersion) for Gaussian distributed continuous variables. Validation results indicate that all three methods produce sharp and reliable probabilistic 2-m temperature forecasts. However, the statistical and combined dynamic-statistical methods slightly outperform the pure dynamical approach, mainly due to the under-dispersive behavior of ALADIN-LAEF outside the nowcasting range. The training length does not have a pronounced impact on forecast skill, but a spread re-scaling improves the forecast skill substantially. Refinements of the statistical methods yield a slight further improvement.

Description: Towards an online-coupled chemistry-climate model: evaluation of trace gases and aerosols in COSMO-ART Geoscientific Model Development, 4, 1077-1102, 2011 Author(s): C. Knote, D. Brunner, H. Vogel, J. Allan, A. Asmi, M. Äijälä, S. Carbone, H. D. van der Gon, J. L. Jimenez, A. Kiendler-Scharr, C. Mohr, L. Poulain, A. S. H. Prévôt, E. Swietlicki, and B. Vogel The online-coupled, regional chemistry transport model COSMO-ART is evaluated for periods in all seasons against several measurement datasets to assess its ability to represent gaseous pollutants and ambient aerosol characteristics over the European domain. Measurements used in the comparison include long-term station observations, satellite and ground-based remote sensing products, and complex datasets of aerosol chemical composition and number size distribution from recent field campaigns. This is the first time these comprehensive measurements of aerosol characteristics in Europe are used to evaluate a regional chemistry transport model. We show a detailed analysis of the simulated size-resolved chemical composition under different meteorological conditions. Mean, variability and spatial distribution of the concentrations of O 3 and NO x are well reproduced. SO 2 is found to be overestimated, simulated PM 2.5 and PM 10 levels are on average underestimated, as is AOD. We find indications of an overestimation of shipping emissions. Time evolution of aerosol chemical composition is captured, although some biases are found in relative composition. Nitrate aerosol components are on average overestimated, and sulfates underestimated. The accuracy of simulated organics depends strongly on season and location. While strongly underestimated during summer, organic mass is comparable in spring and autumn. We see indications for an overestimated fractional contribution of primary organic matter in urban areas and an underestimation of SOA at many locations. Aerosol number concentrations compare well with measurements for larger size ranges, but overestimations of particle number concentration with factors of 2–5 are found for particles smaller than 50 nm. Size distribution characteristics are often close to measurements, but show discrepancies at polluted sites. Suggestions for further improvement of the modeling system consist of the inclusion of a revised secondary organic aerosols scheme, aqueous-phase chemistry and improved aerosol boundary conditions. Our work sets the basis for subsequent studies of aerosol characteristics and climate impacts with COSMO-ART, and highlights areas where improvements are necessary for current regional modeling systems in general.

Description: Identifying the causes of differences in ozone production from the CB05 and CBMIV chemical mechanisms Geoscientific Model Development, 5, 257-268, 2012 Author(s): R. D. Saylor and A. F. Stein An investigation was conducted to identify the mechanistic differences between two versions of the carbon bond gas-phase chemical mechanism (CB05 and CBMIV) which consistently lead to larger ground-level ozone concentrations being produced in the CB05 version of the National Air Quality Forecasting Capability (NAQFC) modeling system even though the two parallel forecast systems utilize the same meteorology and base emissions and similar initial and boundary conditions. Box models of each of the mechanisms as they are implemented in the NAQFC were created and a set of 12 sensitivity simulations was designed. The sensitivity simulations independently probed the conceptual mechanistic differences between CB05 and CBMIV and were exercised over a 45-scenario simulation suite designed to emulate the wide range of chemical regimes encountered in a continental-scale atmospheric chemistry model. Results of the sensitivity simulations indicate that two sets of reactions that were included in the CB05 mechanism, but which were absent from the CBMIV mechanism, are the primary causes of the greater ozone production in the CB05 version of the NAQFC. One set of reactions recycles the higher organic peroxide species of CB05 (ROOH), resulting in additional photochemically reactive products that act to produce additional ozone in some chemical regimes. The other set of reactions recycles reactive nitrogen from less reactive forms back to NO 2 , increasing the effective NO x concentration of the system. In particular, the organic nitrate species (NTR), which was a terminal product for reactive nitrogen in the CBMIV mechanism, acts as a reservoir species in CB05 to redistribute NO x from major source areas to potentially NO x -sensitive areas where additional ozone may be produced in areas remote from direct NO x sources.

Description: Influence of parallel computational uncertainty on simulations of the Coupled General Climate Model Geoscientific Model Development, 5, 313-319, 2012 Author(s): Z. Song, F. Qiao, X. Lei, and C. Wang This paper investigates the impact of the parallel computational uncertainty due to the round-off error on climate simulations using the Community Climate System Model Version 3 (CCSM3). A series of sensitivity experiments have been conducted and the analyses are focused on the Global and Nino3.4 average sea surface temperatures (SST). For the monthly time series, it is shown that the amplitude of the deviation induced by the parallel computational uncertainty is the same order as that of the climate system change. However, the ensemble mean method can reduce the influence and the ensemble member number of 15 is enough to ignore the uncertainty. For climatology, the influence can be ignored when the climatological mean is calculated by using more than 30-yr simulations. It is also found that the parallel computational uncertainty has no distinguishable effect on power spectrum analysis of climate variability such as ENSO. Finally, it is suggested that the influence of the parallel computational uncertainty on Coupled General Climate Models (CGCMs) can be a quality standard or a metric for developing CGCMs.

Description: Wavelet-based spatial comparison technique for analysing and evaluating two-dimensional geophysical model fields Geoscientific Model Development, 5, 223-230, 2012 Author(s): S. Saux Picart, M. Butenschön, and J. D. Shutler Complex numerical models of the Earth's environment, based around 3-D or 4-D time and space domains are routinely used for applications including climate predictions, weather forecasts, fishery management and environmental impact assessments. Quantitatively assessing the ability of these models to accurately reproduce geographical patterns at a range of spatial and temporal scales has always been a difficult problem to address. However, this is crucial if we are to rely on these models for decision making. Satellite data are potentially the only observational dataset able to cover the large spatial domains analysed by many types of geophysical models. Consequently optical wavelength satellite data is beginning to be used to evaluate model hindcast fields of terrestrial and marine environments. However, these satellite data invariably contain regions of occluded or missing data due to clouds, further complicating or impacting on any comparisons with the model. This work builds on a published methodology, that evaluates precipitation forecast using radar observations based on predefined absolute thresholds. It allows model skill to be evaluated at a range of spatial scales and rain intensities. Here we extend the original method to allow its generic application to a range of continuous and discontinuous geophysical data fields, and therefore allowing its use with optical satellite data. This is achieved through two major improvements to the original method: (i) all thresholds are determined based on the statistical distribution of the input data, so no a priori knowledge about the model fields being analysed is required and (ii) occluded data can be analysed without impacting on the metric results. The method can be used to assess a model's ability to simulate geographical patterns over a range of spatial scales. We illustrate how the method provides a compact and concise way of visualising the degree of agreement between spatial features in two datasets. The application of the new method, its handling of bias and occlusion and the advantages of the novel method are demonstrated through the analysis of model fields from a marine ecosystem model.

Description: Spatial patterns of linear and nonparametric long-term trends in Baltic sea-level variability Nonlinear Processes in Geophysics, 19, 95-111, 2012 Author(s): R. V. Donner, R. Ehrcke, S. M. Barbosa, J. Wagner, J. F. Donges, and J. Kurths The study of long-term trends in tide gauge data is important for understanding the present and future risk of changes in sea-level variability for coastal zones, particularly with respect to the ongoing debate on climate change impacts. Traditionally, most corresponding analyses have exclusively focused on trends in mean sea-level. However, such studies are not able to provide sufficient information about changes in the full probability distribution (especially in the more extreme quantiles). As an alternative, in this paper we apply quantile regression (QR) for studying changes in arbitrary quantiles of sea-level variability. For this purpose, we chose two different QR approaches and discuss the advantages and disadvantages of different settings. In particular, traditional linear QR poses very restrictive assumptions that are often not met in reality. For monthly data from 47 tide gauges from along the Baltic Sea coast, the spatial patterns of quantile trends obtained in linear and nonparametric (spline-based) frameworks display marked differences, which need to be understood in order to fully assess the impact of future changes in sea-level variability on coastal areas. In general, QR demonstrates that the general variability of Baltic sea-level has increased over the last decades. Linear quantile trends estimated for sliding windows in time reveal a wide-spread acceleration of trends in the median, but only localised changes in the rates of changes in the lower and upper quantiles.

Description: On closure parameter estimation in chaotic systems Nonlinear Processes in Geophysics, 19, 127-143, 2012 Author(s): J. Hakkarainen, A. Ilin, A. Solonen, M. Laine, H. Haario, J. Tamminen, E. Oja, and H. Järvinen Many dynamical models, such as numerical weather prediction and climate models, contain so called closure parameters. These parameters usually appear in physical parameterizations of sub-grid scale processes, and they act as "tuning handles" of the models. Currently, the values of these parameters are specified mostly manually, but the increasing complexity of the models calls for more algorithmic ways to perform the tuning. Traditionally, parameters of dynamical systems are estimated by directly comparing the model simulations to observed data using, for instance, a least squares approach. However, if the models are chaotic, the classical approach can be ineffective, since small errors in the initial conditions can lead to large, unpredictable deviations from the observations. In this paper, we study numerical methods available for estimating closure parameters in chaotic models. We discuss three techniques: off-line likelihood calculations using filtering methods, the state augmentation method, and the approach that utilizes summary statistics from long model simulations. The properties of the methods are studied using a modified version of the Lorenz 95 system, where the effect of fast variables are described using a simple parameterization.

Description: Comprehensive analysis of tornado statistics in comparison to earthquakes: intensity and temporal behaviour Nonlinear Processes in Geophysics, 20, 47-57, 2013 Author(s): L. Schielicke and P. Névir Tornadoes and earthquakes are characterised by a high variability in their properties concerning intensity, geometric properties and temporal behaviour. Earthquakes are known for power-law behaviour in their intensity (Gutenberg–Richter law) and temporal statistics (e.g. Omori law and interevent waiting times). The observed similarity of high variability of these two phenomena motivated us to compare the statistical behaviour of tornadoes using seismological methods and quest for power-law behaviour. In general, the statistics of tornadoes show power-law behaviour partly coextensive with characteristic scales when the temporal resolution is high (10 to 60 min). These characteristic scales match with the typical diurnal behaviour of tornadoes, which is characterised by a maximum of tornado occurrences in the late afternoon hours. Furthermore, the distributions support the observation that tornadoes cluster in time. Finally, we shortly discuss a possible similar underlying structure composed of heterogeneous, coupled, interactive threshold oscillators that possibly explains the observed behaviour.

Description: Quantifying soil complexity using network models of soil porous structure Nonlinear Processes in Geophysics, 20, 41-45, 2013 Author(s): M. Samec, A. Santiago, J. P. Cárdenas, R. M. Benito, A. M. Tarquis, S. J. Mooney, and D. Korošak This paper describes an investigation into the properties of spatially embedded complex networks representing the porous architecture of soil systems. We suggest an approach to quantify the complexity of soil pore structure based on the node-node link correlation properties of the networks. We show that the complexity depends on the strength of spatial embedding of the network and that this is related to the transition from a non-compact to compact phase of the network.

Description: Analysis of asymmetries in propagating mode-2 waves Nonlinear Processes in Geophysics, 20, 59-69, 2013 Author(s): J. Olsthoorn, A. Baglaenko, and M. Stastna Using numerical simulations performed with a pseudo-spectral incompressible Navier–Stokes solver, we describe the asymmetries that arise in the recirculating core of mode-2 internal, solitary-like waves. The waves are generated in a manner consistent with many laboratory studies, namely via the collapse of a region of mixed fluid. Analysis of the simulations reveals that asymmetries across both the wave crest and the pycnocline centre develop in the spatial distribution of density, kinetic energy and a passive tracer transported by the mode-2 waves. The simulations are extended to three-dimensions to allow for the formation of spanwise instabilities. We find that three-dimensionalization modifies the structure and energetics of the core, but that the majority of the results obtained from two dimensional simulations remain valid. Taken together, our simulations demonstrate that the cores of solitary-like mode-2 waves are different then their counterparts for mode-1 waves and that their accurate characterization on both lab and field scales should account for the core asymmetry across the pycnocline centre.

Description: Performance of McRAS-AC in the GEOS-5 AGCM: aerosol-cloud-microphysics, precipitation, cloud radiative effects, and circulation Geoscientific Model Development, 6, 57-79, 2013 Author(s): Y. C. Sud, D. Lee, L. Oreopoulos, D. Barahona, A. Nenes, and M. J. Suarez A revised version of the Microphysics of clouds with Relaxed Arakawa-Schubert and Aerosol-Cloud interaction scheme (McRAS-AC) including, among others, a new ice nucleation parameterization, is implemented in the GEOS-5 AGCM. Various fields from a 10-yr-long integration of the AGCM with McRAS-AC are compared with their counterparts from an integration of the baseline GEOS-5 AGCM, as well as satellite observations. Generally McRAS-AC simulations have smaller biases in cloud fields and cloud radiative effects over most of the regions of the Earth than the baseline GEOS-5 AGCM. Two systematic biases are identified in the McRAS-AC runs: one is underestimation of cloud particle numbers around 40° S–60° S, and one is overestimate of cloud water path during the Northern Hemisphere summer over the Gulf Stream and North Pacific. Sensitivity tests show that these biases potentially originate from biases in the aerosol input. The first bias is largely eliminated in a test run using 50% smaller radius of sea-salt aerosol particles, while the second bias is substantially reduced when interactive aerosol chemistry is turned on. The main weakness of McRAS-AC is the dearth of low-level marine stratus clouds, a probable outcome of lack of explicit dry-convection in the cloud scheme. Nevertheless, McRAS-AC largely simulates realistic clouds and their optical properties that can be improved further with better aerosol input. An assessment using the COSP simulator in a 1-yr integration provides additional perspectives for understanding cloud optical property differences between the baseline and McRAS-AC simulations and biases against satellite data. Overall, McRAS-AC physically couples aerosols, the microphysics and macrophysics of clouds, and their radiative effects and thereby has better potential to be a valuable tool for climate modeling research.

Description: Corrigendum to "CLM4-BeTR, a generic biogeochemical transport and reaction module for CLM4: model development, evaluation, and application" published in Geosci. Model Dev., 6, 127–140, 2013 Geoscientific Model Development, 6, 245-245, 2013 Author(s): J. Y. Tang, W. J. Riley, C. D. Koven, and Z. M. Subin No abstract available.

Description: Implementation of the chemistry module MECCA (v2.5) in the modal aerosol version of the Community Atmosphere Model component (v3.6.33) of the Community Earth System Model Geoscientific Model Development, 6, 255-262, 2013 Author(s): M. S. Long, W. C. Keene, R. Easter, R. Sander, A. Kerkweg, D. Erickson, X. Liu, and S. Ghan A coupled atmospheric chemistry and climate system model was developed using the modal aerosol version of the National Center for Atmospheric Research Community Atmosphere Model (modal-CAM; v3.6.33) and the Max Planck Institute for Chemistry's Module Efficiently Calculating the Chemistry of the Atmosphere (MECCA; v2.5) to provide enhanced resolution of multiphase processes, particularly those involving inorganic halogens, and associated impacts on atmospheric composition and climate. Three Rosenbrock solvers (Ros-2, Ros-3, RODAS-3) were tested in conjunction with the basic load-balancing options available to modal-CAM (1) to establish an optimal configuration of the implicitly-solved multiphase chemistry module that maximizes both computational speed and repeatability of Ros-2 and RODAS-3 results versus Ros-3, and (2) to identify potential implementation strategies for future versions of this and similar coupled systems. RODAS-3 was faster than Ros-2 and Ros-3 with good reproduction of Ros-3 results, while Ros-2 was both slower and substantially less reproducible relative to Ros-3 results. Modal-CAM with MECCA chemistry was a factor of 15 slower than modal-CAM using standard chemistry. MECCA chemistry integration times demonstrated a systematic frequency distribution for all three solvers, and revealed that the change in run-time performance was due to a change in the frequency distribution of chemical integration times; the peak frequency was similar for all solvers. This suggests that efficient chemistry-focused load-balancing schemes can be developed that rely on the parameters of this frequency distribution.

Description: Ion motion in the current sheet with sheared magnetic field – Part 1: Quasi-adiabatic theory Nonlinear Processes in Geophysics, 20, 163-178, 2013 Author(s): A. V. Artemyev, A. I. Neishtadt, and L. M. Zelenyi We present a theory of trapped ion motion in the magnetotail current sheet with a constant dawn–dusk component of the magnetic field. Particle trajectories are described analytically using the quasi-adiabatic invariant corresponding to averaging of fast oscillations around the tangential component of the magnetic field. We consider particle dynamics in the quasi-adiabatic approximation and demonstrate that the principal role is played by large (so called geometrical) jumps of the quasi-adiabatic invariant. These jumps appear due to the current sheet asymmetry related to the presence of the dawn–dusk magnetic field. The analytical description is compared with results of numerical integration. We show that there are four possible regimes of particle motion. Each regime is characterized by certain ranges of values of the dawn–dusk magnetic field and particle energy. We find the critical value of the dawn–dusk magnetic field, where jumps of the quasi-adiabatic invariant vanish.

Description: Implementation of multirate time integration methods for air pollution modelling Geoscientific Model Development, 5, 1395-1405, 2012 Author(s): M. Schlegel, O. Knoth, M. Arnold, and R. Wolke Explicit time integration methods are characterised by a small numerical effort per time step. In the application to multiscale problems in atmospheric modelling, this benefit is often more than compensated by stability problems and step size restrictions resulting from stiff chemical reaction terms and from a locally varying Courant-Friedrichs-Lewy (CFL) condition for the advection terms. Splitting methods may be applied to efficiently combine implicit and explicit methods (IMEX splitting). Complementarily multirate time integration schemes allow for a local adaptation of the time step size to the grid size. In combination, these approaches lead to schemes which are efficient in terms of evaluations of the right-hand side. Special challenges arise when these methods are to be implemented. For an efficient implementation, it is crucial to locate and exploit redundancies. Furthermore, the more complex programme flow may lead to computational overhead which, in the worst case, more than compensates the theoretical gain in efficiency. We present a general splitting approach which allows both for IMEX splittings and for local time step adaptation. The main focus is on an efficient implementation of this approach for parallel computation on computer clusters.

Description: Conditioning model output statistics of regional climate model precipitation on circulation patterns Nonlinear Processes in Geophysics, 19, 623-633, 2012 Author(s): F. Wetterhall, F. Pappenberger, Y. He, J. Freer, and H. L. Cloke Dynamical downscaling of Global Climate Models (GCMs) through regional climate models (RCMs) potentially improves the usability of the output for hydrological impact studies. However, a further downscaling or interpolation of precipitation from RCMs is often needed to match the precipitation characteristics at the local scale. This study analysed three Model Output Statistics (MOS) techniques to adjust RCM precipitation; (1) a simple direct method (DM), (2) quantile-quantile mapping (QM) and (3) a distribution-based scaling (DBS) approach. The modelled precipitation was daily means from 16 RCMs driven by ERA40 reanalysis data over the 1961–2000 provided by the ENSEMBLES (ENSEMBLE-based Predictions of Climate Changes and their Impacts) project over a small catchment located in the Midlands, UK. All methods were conditioned on the entire time series, separate months and using an objective classification of Lamb's weather types. The performance of the MOS techniques were assessed regarding temporal and spatial characteristics of the precipitation fields, as well as modelled runoff using the HBV rainfall-runoff model. The results indicate that the DBS conditioned on classification patterns performed better than the other methods, however an ensemble approach in terms of both climate models and downscaling methods is recommended to account for uncertainties in the MOS methods.

Description: Lidar signal simulation for the evaluation of aerosols in chemistry transport models Geoscientific Model Development, 5, 1543-1564, 2012 Author(s): S. Stromatas, S. Turquety, L. Menut, H. Chepfer, J. C. Péré, G. Cesana, and B. Bessagnet We present an adaptable tool, the OPTSIM (OPTical properties SIMulation) software, for the simulation of optical properties and lidar attenuated backscattered profiles ( β ') from aerosol concentrations calculated by chemistry transport models (CTM). It was developed to model both Level 1 observations and Level 2 aerosol lidar retrievals in order to compare model results to measurements: the level 2 enables to estimate the main properties of aerosols plume structures, but may be limited due to specific assumptions. The level 1, originally developed for this tool, gives access to more information about aerosols properties ( β ') requiring, at the same time, less hypothesis on aerosols types. In addition to an evaluation of the aerosol loading and optical properties, active remote sensing allows the analysis of aerosols' vertical structures. An academic case study for two different species (black carbon and dust) is presented and shows the consistency of the simulator. Illustrations are then given through the analysis of dust events in the Mediterranean region during the summer 2007. These are based on simulations by the CHIMERE regional CTM and observations from the CALIOP space-based lidar, and highlight the potential of this approach to evaluate the concentration, size and vertical structure of the aerosol plumes.

Description: A typical wave wake from high-speed vessels: its group structure and run-up Nonlinear Processes in Geophysics, 20, 179-188, 2013 Author(s): I. Didenkulova and A. Rodin High-amplitude water waves induced by high-speed vessels are regularly observed in Tallinn Bay, the Baltic Sea, causing intense beach erosion and disturbing marine habitants in the coastal zone. Such a strong impact on the coast may be a result of a certain group structure of the wave wake. In order to understand it, here we present an experimental study of the group structure of these wakes at Pikakari beach, Tallinn Bay. The most energetic vessel waves at this location (100 m from the coast at the water depth 2.7 m) have amplitudes of about 1 m and periods of 8–10 s and cause maximum run-up heights on a beach up to 1.4 m. These waves represent frequency modulated packets where the largest and longest waves propagate ahead of other smaller amplitude and period waves. Sometimes the groups of different heights and periods can be separated even within one wave wake event. The wave heights within a wake are well described by the Weibull distribution, which has different parameters for wakes from different vessels. Wave run-up heights can also be described by Weibull distribution and its parameters can be connected to the parameters of the distribution of wave heights 100 m from the coast. Finally, the run-up of individual waves within a packet is studied. It is shown that the specific structure of frequency modulated wave packets, induced by high-speed vessels, leads to a sequence of high wave run-ups at the coast, even when the original wave heights are rather moderate. This feature can be a key to understanding the significant impact on coasts caused by fast vessels.

Description: Modelling the Caspian Sea and its catchment area using a coupled regional atmosphere-ocean model (RegCM4-ROMS): model design and preliminary results Geoscientific Model Development, 6, 283-299, 2013 Author(s): U. U. Turuncoglu, G. Giuliani, N. Elguindi, and F. Giorgi We describe the development of a coupled regional atmosphere-ocean model (RegCM4-ROMS) and its implementation over the Caspian Sea basin. The coupled model is run for the period 1999–2008 (after a spin up of 4 yr) and it is compared to corresponding stand alone model simulations and a simulation in which a distributed 1d lake model is run for the Caspian Sea. All model versions show a good performance in reproducing the climatology of the Caspian Sea basin, with relatively minor differences across them. The coupled ROMS produces realistic, although somewhat overestimated, Caspian Sea Surface Temperature (SST), with a considerable improvement compared to the use of the simpler coupled lake model. Simulated near surface salinity and sea currents are also realistic, although the upwelling over the eastern coastal regions is underestimated. The sea ice extent over the shallow northern shelf of the Caspian Sea and its seasonal evolution are well reproduced, however, a significant negative bias in sea-ice fraction exists due to the relatively poor representation of the bathymetry. ROMS also calculates the Caspian Sea Level (CSL), showing that for the present experiment excessive evaporation over the lake area leads to a drift in estimated CSL. Despite this problem, which requires further analysis due to many uncertainties in the estimation of CSL, overall the coupled RegCM4-ROMS system shows encouraging results in reproducing both the climatology of the region and the basic characteristics of the Caspian Sea.

Description: A new marine ecosystem model for the University of Victoria Earth System Climate Model Geoscientific Model Development, 5, 1195-1220, 2012 Author(s): D. P. Keller, A. Oschlies, and M. Eby Earth System Climate Models (ESCMs) are valuable tools that can be used to gain a better understanding of the climate system, global biogeochemical cycles and how anthropogenically-driven changes may affect them. Here we describe improvements made to the marine biogeochemical ecosystem component of the University of Victoria's ESCM (version 2.9). Major changes include corrections to the code and equations describing phytoplankton light limitation and zooplankton grazing, the implementation of a more realistic zooplankton growth and grazing model, and the implementation of an iron limitation scheme to constrain phytoplankton growth. The new model is evaluated after a 10 000-yr spin-up and compared to both the previous version and observations. For the majority of biogeochemical tracers and ecosystem processes the new model shows significant improvements when compared to the previous version and evaluated against observations. Many of the improvements are due to better simulation of seasonal changes in higher latitude ecosystems and the effect that this has on ocean biogeochemistry. This improved model is intended to provide a basic new ESCM model component, which can be used as is or expanded upon (i.e., the addition of new tracers), for climate change and biogeochemical cycling research.

Description: Downscaling the climate change for oceans around Australia Geoscientific Model Development, 5, 1177-1194, 2012 Author(s): M. A. Chamberlain, C. Sun, R. J. Matear, M. Feng, and S. J. Phipps At present, global climate models used to project changes in climate poorly resolve mesoscale ocean features such as boundary currents and eddies. These missing features may be important to realistically project the marine impacts of climate change. Here we present a framework for dynamically downscaling coarse climate change projections utilising a near-global ocean model that resolves these features in the Australasian region, with coarser resolution elsewhere. A time-slice projection for a 2060s ocean was obtained by adding climate change anomalies to initial conditions and surface fluxes of a near-global eddy-resolving ocean model. Climate change anomalies are derived from the differences between present and projected climates from a coarse global climate model. These anomalies are added to observed fields, thereby reducing the effect of model bias from the climate model. The downscaling model used here is ocean-only and does not include the effects that changes in the ocean state will have on the atmosphere and air–sea fluxes. We use restoring of the sea surface temperature and salinity to approximate real-ocean feedback on heat flux and to keep the salinity stable. Extra experiments with different feedback parameterisations are run to test the sensitivity of the projection. Consistent spatial differences emerge in sea surface temperature, salinity, stratification and transport between the downscaled projections and those of the climate model. Also, the spatial differences become established rapidly ( 〈 3 yr), indicating the importance of mesoscale resolution. However, the differences in the magnitude of the difference between experiments show that feedback of the ocean onto the air–sea fluxes is still important in determining the state of the ocean in these projections. Until such a time when it is feasible to regularly run a global climate model with eddy resolution, our framework for ocean climate change downscaling provides an attractive way to explore the response of mesoscale ocean features with climate change and their effect on the broader ocean.

Description: Description of a hybrid ice sheet-shelf model, and application to Antarctica Geoscientific Model Development, 5, 1273-1295, 2012 Author(s): D. Pollard and R. M. DeConto The formulation of a 3-D ice sheet-shelf model is described. The model is designed for long-term continental-scale applications, and has been used mostly in paleoclimatic studies. It uses a hybrid combination of the scaled shallow ice and shallow shelf approximations for ice flow. Floating ice shelves and grounding-line migration are included, with parameterized ice fluxes at grounding lines that allows relatively coarse resolutions to be used. All significant components and parameterizations of the model are described in some detail. Basic results for modern Antarctica are compared with observations, and simulations over the last 5 million years are compared with previously published results. The sensitivity of ice volumes during the last deglaciation to basal sliding coefficients is discussed.

Description: The Nexus Land-Use model version 1.0, an approach articulating biophysical potentials and economic dynamics to model competition for land-use Geoscientific Model Development, 5, 1297-1322, 2012 Author(s): F. Souty, T. Brunelle, P. Dumas, B. Dorin, P. Ciais, R. Crassous, C. Müller, and A. Bondeau Interactions between food demand, biomass energy and forest preservation are driving both food prices and land-use changes, regionally and globally. This study presents a new model called Nexus Land-Use version 1.0 which describes these interactions through a generic representation of agricultural intensification mechanisms within agricultural lands. The Nexus Land-Use model equations combine biophysics and economics into a single coherent framework to calculate crop yields, food prices, and resulting pasture and cropland areas within 12 regions inter-connected with each other by international trade. The representation of cropland and livestock production systems in each region relies on three components: (i) a biomass production function derived from the crop yield response function to inputs such as industrial fertilisers; (ii) a detailed representation of the livestock production system subdivided into an intensive and an extensive component, and (iii) a spatially explicit distribution of potential (maximal) crop yields prescribed from the Lund-Postdam-Jena global vegetation model for managed Land (LPJmL). The economic principles governing decisions about land-use and intensification are adapted from the Ricardian rent theory, assuming cost minimisation for farmers. In contrast to the other land-use models linking economy and biophysics, crops are aggregated as a representative product in calories and intensification for the representative crop is a non-linear function of chemical inputs. The model equations and parameter values are first described in details. Then, idealised scenarios exploring the impact of forest preservation policies or rising energy price on agricultural intensification are described, and their impacts on pasture and cropland areas are investigated.

Description: Fractal analysis of ULF electromagnetic emissions in possible association with earthquakes in China Nonlinear Processes in Geophysics, 19, 577-583, 2012 Author(s): Y. Ida, D. Yang, Q. Li, H. Sun, and M. Hayakawa The long-term data (during the period of 1 March 2003 through 31 December 2006) of ULF geomagnetic variations observed at Kashi station (geographic coordinates: 39.5° N, 76.0° E) in China have been used to investigate the long-term variation of fractal dimension of ULF emissions. We have studied the changes in fractal dimension in association with several earthquakes around the observation station. It is then found that a significant change (or decrease) in the fractal dimension of the Z component took place before the 1 September 2003 earthquake, which lends a further support to our previous finding based on our improved polarization analysis for the same earthquake. The results obtained are discussed in the contexts of a few aspects (detectability of seismogenic emissions, comparison with previous results by other analysis methods, the importance of fractal analysis in the nonlinear process of the lithosphere and earthquake prediction).

Description: "Gtool5": a Fortran90 library of input/output interfaces for self-descriptive multi-dimensional numerical data Geoscientific Model Development, 5, 449-455, 2012 Author(s): M. Ishiwatari, E. Toyoda, Y. Morikawa, S. Takehiro, Y. Sasaki, S. Nishizawa, M. Odaka, N. Otobe, Y. O. Takahashi, K. Nakajima, T. Horinouchi, M. Shiotani, Y.-Y. Hayashi, and Gtool development group A Fortran90 input/output library, "gtool5", is developed for use with numerical simulation models in the fields of Earth and planetary sciences. The use of this library will simplify implementation of input/output operations into program code in a consolidated form independent of the size and complexity of the software and data. The library also enables simple specification of the metadata needed for post-processing and visualization of the data. These aspects improve the readability of simulation code, which facilitates the simultaneous performance of multiple numerical experiments with different software and efficiency in examining and comparing the numerical results. The library is expected to provide a common software platform to reinforce research on, for instance, the atmosphere and ocean, where a close combination of multiple simulation models with a wide variety of complexity of physics implementations from massive climate models to simple geophysical fluid dynamics models is required.

Description: Predictability of a low-order interactive ensemble Nonlinear Processes in Geophysics, 19, 273-282, 2012 Author(s): L. Siqueira and B. Kirtman In this paper, numerical and analytical studies were performed to uncover the mechanisms controlling the changes in ensemble spread of a low-order coupled model with multiple atmospheric realizations. An interactive ensemble approach was applied to a coupled dynamical system based on two versions of the Lorenz 63 model designed in order to imitate the behavior of a coupled system with different time scales. In the dynamic system used in this work the spread of ensemble members is highly dependent on the mean state corresponding to asymmetries in predictability. The slowness of the slow model and the intensity of the boundary forcing anomalies both contribute to the asymmetry and phase locking of both subsystems. The mechanisms controlling the fast model spread were uncovered revealing uncertainty dynamics depending on the location of ensemble members in the fast model phase space and implicitly on the slowness and magnitude of the slow model anomalies.

Description: Detection, tracking and event localization of jet stream features in 4-D atmospheric data Geoscientific Model Development, 5, 457-470, 2012 Author(s): S. Limbach, E. Schömer, and H. Wernli We introduce a novel algorithm for the efficient detection and tracking of features in spatiotemporal atmospheric data, as well as for the precise localization of the occurring genesis, lysis, merging and splitting events. The algorithm works on data given on a four-dimensional structured grid. Feature selection and clustering are based on adjustable local and global criteria, feature tracking is predominantly based on spatial overlaps of the feature's full volumes. The resulting 3-D features and the identified correspondences between features of consecutive time steps are represented as the nodes and edges of a directed acyclic graph, the event graph. Merging and splitting events appear in the event graph as nodes with multiple incoming or outgoing edges, respectively. The precise localization of the splitting events is based on a search for all grid points inside the initial 3-D feature that have a similar distance to two successive 3-D features of the next time step. The merging event is localized analogously, operating backward in time. As a first application of our method we present a climatology of upper-tropospheric jet streams and their events, based on four-dimensional wind speed data from European Centre for Medium-Range Weather Forecasts (ECMWF) analyses. We compare our results with a climatology from a previous study, investigate the statistical distribution of the merging and splitting events, and illustrate the meteorological significance of the jet splitting events with a case study. A brief outlook is given on additional potential applications of the 4-D data segmentation technique.

Description: Aerosol–climate interactions in the Norwegian Earth System Model – NorESM1-M Geoscientific Model Development, 6, 207-244, 2013 Author(s): A. Kirkevåg, T. Iversen, Ø. Seland, C. Hoose, J. E. Kristjánsson, H. Struthers, A. M. L. Ekman, S. Ghan, J. Griesfeller, E. D. Nilsson, and M. Schulz The objective of this study is to document and evaluate recent changes and updates to the module for aerosols and aerosol–cloud–radiation interactions in the atmospheric module CAM4-Oslo of the core version of the Norwegian Earth System Model (NorESM), NorESM1-M. Particular attention is paid to the role of natural organics, sea salt, and mineral dust in determining the gross aerosol properties as well as the anthropogenic contribution to these properties and the associated direct and indirect radiative forcing. The aerosol module is extended from earlier versions that have been published, and includes life-cycling of sea salt, mineral dust, particulate sulphate, black carbon, and primary and secondary organics. The impacts of most of the numerous changes since previous versions are thoroughly explored by sensitivity experiments. The most important changes are: modified prognostic sea salt emissions; updated treatment of precipitation scavenging and gravitational settling; inclusion of biogenic primary organics and methane sulphonic acid (MSA) from oceans; almost doubled production of land-based biogenic secondary organic aerosols (SOA); and increased ratio of organic matter to organic carbon (OM/OC) for biomass burning aerosols from 1.4 to 2.6. Compared with in situ measurements and remotely sensed data, the new treatments of sea salt and dust aerosols give smaller biases in near-surface mass concentrations and aerosol optical depth than in the earlier model version. The model biases for mass concentrations are approximately unchanged for sulphate and BC. The enhanced levels of modeled OM yield improved overall statistics, even though OM is still underestimated in Europe and overestimated in North America. The global anthropogenic aerosol direct radiative forcing (DRF) at the top of the atmosphere has changed from a small positive value to −0.08 W m −2 in CAM4-Oslo. The sensitivity tests suggest that this change can be attributed to the new treatment of biomass burning aerosols and gravitational settling. Although it has not been a goal in this study, the new DRF estimate is closer both to the median model estimate from the AeroCom intercomparison and the best estimate in IPCC AR4. Estimated DRF at the ground surface has increased by ca. 60%, to −1.89 W m −2 . We show that this can be explained by new emission data and omitted mixing of constituents between updrafts and downdrafts in convective clouds. The increased abundance of natural OM and the introduction of a cloud droplet spectral dispersion formulation are the most important contributions to a considerably decreased estimate of the indirect radiative forcing (IndRF). The IndRF is also found to be sensitive to assumptions about the coating of insoluble aerosols by sulphate and OM. The IndRF of −1.2 W m −2 , which is closer to the IPCC AR4 estimates than the previous estimate of −1.9 W m −2 , has thus been obtained without imposing unrealistic artificial lower bounds on cloud droplet number concentrations.

Description: Interaction of a monopole vortex with an isolated topographic feature in a three-layer geophysical flow Nonlinear Processes in Geophysics, 20, 107-119, 2013 Author(s): E. A. Ryzhov and K. V. Koshel In the frame of a three-layer, quasi-geostrophic analytical model of an f -plane geophysical flow, the Lagrangian advection induced by the interaction of a monopole vortex with an isolated topographic feature is addressed. Two different cases when the monopole is located either within the upper or the middle layer are of our interest. In the bottom layer, there is a delta-function topographic feature, which generates a closed recirculation region in its vicinity due to the background flow. This recirculation region extends to the middle and upper layers, and it plays the role of a topographic vortex. The interaction between the monopole and the topographic vortex causes a complex, including chaotic, advection of fluid particles. We show that the model's parameters, namely the monopole and topographic vortices' strengths and initial positions, and the layers' depths and densities, are responsible for the diverse advection patterns. While the patterns are rather complicated, one can single out two major processes, which mostly govern the fluid particle advection. The first one is the variation in time of the system's phase space structure, so that within the closed region of the topographic vortex, there appear periodically unclosed particle pathways by which the particles leave the topographic vortex. The second one is chaotic advection that arises from the nonstationarity of the monopole–topography interaction.

Description: The Atmospheric Chemistry and Climate Model Intercomparison Project (ACCMIP): overview and description of models, simulations and climate diagnostics Geoscientific Model Development, 6, 179-206, 2013 Author(s): J.-F. Lamarque, D. T. Shindell, B. Josse, P. J. Young, I. Cionni, V. Eyring, D. Bergmann, P. Cameron-Smith, W. J. Collins, R. Doherty, S. Dalsoren, G. Faluvegi, G. Folberth, S. J. Ghan, L. W. Horowitz, Y. H. Lee, I. A. MacKenzie, T. Nagashima, V. Naik, D. Plummer, M. Righi, S. T. Rumbold, M. Schulz, R. B. Skeie, D. S. Stevenson, S. Strode, K. Sudo, S. Szopa, A. Voulgarakis, and G. Zeng The Atmospheric Chemistry and Climate Model Intercomparison Project (ACCMIP) consists of a series of time slice experiments targeting the long-term changes in atmospheric composition between 1850 and 2100, with the goal of documenting composition changes and the associated radiative forcing. In this overview paper, we introduce the ACCMIP activity, the various simulations performed (with a requested set of 14) and the associated model output. The 16 ACCMIP models have a wide range of horizontal and vertical resolutions, vertical extent, chemistry schemes and interaction with radiation and clouds. While anthropogenic and biomass burning emissions were specified for all time slices in the ACCMIP protocol, it is found that the natural emissions are responsible for a significant range across models, mostly in the case of ozone precursors. The analysis of selected present-day climate diagnostics (precipitation, temperature, specific humidity and zonal wind) reveals biases consistent with state-of-the-art climate models. The model-to-model comparison of changes in temperature, specific humidity and zonal wind between 1850 and 2000 and between 2000 and 2100 indicates mostly consistent results. However, models that are clear outliers are different enough from the other models to significantly affect their simulation of atmospheric chemistry.

Description: Implementation of the Fast-JX Photolysis scheme (v6.4) into the UKCA component of the MetUM chemistry-climate model (v7.3) Geoscientific Model Development, 6, 161-177, 2013 Author(s): P. J. Telford, N. L. Abraham, A. T. Archibald, P. Braesicke, M. Dalvi, O. Morgenstern, F. M. O'Connor, N. A. D. Richards, and J. A. Pyle Atmospheric chemistry is driven by photolytic reactions, making their modelling a crucial component of atmospheric models. We describe the implementation and validation of Fast-JX, a state of the art model of interactive photolysis, into the MetUM chemistry-climate model. This allows for interactive photolysis rates to be calculated in the troposphere and augments the calculation of the rates in the stratosphere by accounting for clouds and aerosols in addition to ozone. In order to demonstrate the effectiveness of this new photolysis scheme we employ new methods of validating the model, including techniques for sampling the model to compare to flight track and satellite data.

Description: A generalized tagging method Geoscientific Model Development, 6, 247-253, 2013 Author(s): V. Grewe The understanding of causes of changes in climate-chemistry simulations is an important, but often challenging task. In atmospheric chemistry, one approach is to tag species according to their origin (e.g. emission categories) and to inherit these tags to other species during subsequent reactions. This concept was recently employed to calculate the contribution of atmospheric processes to temperature. Here a new concept for tagging any state variable is presented. This generalized tagging method results from a sensitivity analysis of the individual forcing terms of the right hand side of the governing differential equations. In a couple of examples, the consistency with previous approaches and the synergy by using different analysis techniques are shown. Since the method is based on a ratio describing relative sensitivities, singularities occur where the method is not applicable. For some applications, such as in atmospheric chemistry, these singularities can easily be removed. However, one theoretical example is given, where this method is not applicable at all.

Description: Characterization of turbulence stability through the identification of multifractional Brownian motions Nonlinear Processes in Geophysics, 20, 97-106, 2013 Author(s): K. C. Lee Multifractional Brownian motions have become popular as flexible models in describing real-life signals of high-frequency features in geoscience, microeconomics, and turbulence, to name a few. The time-changing Hurst exponent, which describes regularity levels depending on time measurements, and variance, which relates to an energy level, are two parameters that characterize multifractional Brownian motions. This research suggests a combined method of estimating the time-changing Hurst exponent and variance using the local variation of sampled paths of signals. The method consists of two phases: initially estimating global variance and then accurately estimating the time-changing Hurst exponent. A simulation study shows its performance in estimation of the parameters. The proposed method is applied to characterization of atmospheric stability in which descriptive statistics from the estimated time-changing Hurst exponent and variance classify stable atmosphere flows from unstable ones.

Description: A tri-stage cluster identification model for accurate analysis of seismic catalogs Nonlinear Processes in Geophysics, 20, 143-162, 2013 Author(s): S. J. Nanda, K. F. Tiampo, G. Panda, L. Mansinha, N. Cho, and A. Mignan In this paper we propose a tri-stage cluster identification model that is a combination of a simple single iteration distance algorithm and an iterative K-means algorithm. In this study of earthquake seismicity, the model considers event location, time and magnitude information from earthquake catalog data to efficiently classify events as either background or mainshock and aftershock sequences. Tests on a synthetic seismicity catalog demonstrate the efficiency of the proposed model in terms of accuracy percentage (94.81% for background and 89.46% for aftershocks). The close agreement between lambda and cumulative plots for the ideal synthetic catalog and that generated by the proposed model also supports the accuracy of the proposed technique. There is flexibility in the model design to allow for proper selection of location and magnitude ranges, depending upon the nature of the mainshocks present in the catalog. The effectiveness of the proposed model also is evaluated by the classification of events in three historic catalogs: California, Japan and Indonesia. As expected, for both synthetic and historic catalog analysis it is observed that the density of events classified as background is almost uniform throughout the region, whereas the density of aftershock events are higher near the mainshocks.

Description: Scaling for lobe and cleft patterns in particle-laden gravity currents Nonlinear Processes in Geophysics, 20, 121-130, 2013 Author(s): A. Jackson, B. Turnbull, and R. Munro Lobe and cleft patterns are frequently observed at the leading edge of gravity currents, including non-Boussinesq particle-laden currents such as powder snow avalanches. Despite the importance of the instability in driving air entrainment, little is known about its origin or the mechanisms behind its development. In this paper we seek to gain a better understanding of these mechanisms from a laboratory scale model of powder snow avalanches using lightweight granular material. The instability mechanisms in these flows appear to be a combination of those found in both homogeneous Boussinesq gravity currents and unsuspended granular flows, with the size of the granular particles playing a central role in determining the wavelength of the lobe and cleft pattern. When scaled by particle diameter a relationship between the Froude number and the wavelength of the lobe and cleft pattern is found, where the wavelength increases monotonically with the Froude number.

Description: Inversion of the western Pacific subtropical high dynamic model and analysis of dynamic characteristics for its abnormality Nonlinear Processes in Geophysics, 20, 131-142, 2013 Author(s): M. Hong, R. Zhang, J. X. Li, J. J. Ge, and K. F. Liu Based on time series data of 500 hPa potential field from NCEP/NCAR (National Center for Environmental Forecast of American/National Center for Atmospheric Research), a novel consideration of empirical orthogonal function (EOF) time–space separation and dynamic system reconstruction for time series is introduced. This method consists of two parts: first, the dynamical model inversion and model parameter optimization are carried out on the EOF time coefficient series using the genetic algorithm (GA), and, second, a nonlinear dynamic model representing the subtropical high (SH) activity and its abnormality is established. The SH activity and its abnormal mechanism is studied using the developed dynamical model. Results show that the configuration and diversification of the SH equilibriums have good correspondence with the actual short–medium term abnormal activity of the SH. Change of SH potential field brought by the combination of equilibriums is more complex than that by mutation, and their exhibition patterns are different. The mutation behavior from high-value to low-value equilibriums of the SH in summer corresponds with the southward drop of the SH in the observed weather process. The combination behavior of the two steady equilibriums corresponds with disappearance of the "double-ridge" phenomenon of the SH. Dynamical mechanisms of these phenomena are explained.

Description: Wavefield decomposition and phase space dynamics of the seismic noise at Volcàn de Colima, Mexico: evidence of a two-state source process Nonlinear Processes in Geophysics, 20, 71-84, 2013 Author(s): M. Palo and P. Cusano We analyse the seismic noise recorded at the Colima Volcano (Mexico) in the period December 2005–May 2006 by four broadband three-component seismic stations. Specifically, we characterize the spectral content of the signal and follow its time evolution along all the data set. Moreover, we infer the properties of the attractor in the phase space by false nearest neighbours analysis and Grassberger–Procaccia algorithm, and adopt a time-domain decomposition method (independent component analysis) to find the basic constituents (independent components) of the system. Constraints on the seismic wavefield are inferred by the polarization analysis. We find two states of the background seismicity visible in different time-intervals that are Phase A and Phase B. Phase A has a spectrum with two peaks at 0.15 Hz and 0.3 Hz, with the latter dominating, an attractor of correlation dimension close to 3, three quasi-monochromatic independent components, and a relevant fraction of crater-pointing polarization solutions in the near-field. In Phase B, the spectrum is preserved but with the highest peak at 0.15 Hz, the attractor has a correlation dimension close to 2, two independent components are extracted, and the polarization solutions are dominated by Rayleigh waves incoming from the southwest direction. We depict two sources acting on the background seismicity that are the microseismic noise loading on the Pacific coastline and a low-energy volcanic tremor. A change in the amplitude of the microseismic noise can induce the switching from a state of the system to the other.

Description: Optimising the FAMOUS climate model: inclusion of global carbon cycling Geoscientific Model Development, 6, 141-160, 2013 Author(s): J. H. T. Williams, R. S. Smith, P. J. Valdes, B. B. B. Booth, and A. Osprey FAMOUS fills an important role in the hierarchy of climate models, both explicitly resolving atmospheric and oceanic dynamics yet being sufficiently computationally efficient that either very long simulations or large ensembles are possible. An improved set of carbon cycle parameters for this model has been found using a perturbed physics ensemble technique. This is an important step towards building the "Earth System" modelling capability of FAMOUS, which is a reduced resolution, and hence faster running, version of the Hadley Centre Climate model, HadCM3. Two separate 100 member perturbed parameter ensembles were performed; one for the land surface and one for the ocean. The land surface scheme was tested against present-day and past representations of vegetation and the ocean ensemble was tested against observations of nitrate. An advantage of using a relatively fast climate model is that a large number of simulations can be run and hence the model parameter space (a large source of climate model uncertainty) can be more thoroughly sampled. This has the associated benefit of being able to assess the sensitivity of model results to changes in each parameter. The climatologies of surface and tropospheric air temperature and precipitation are improved relative to previous versions of FAMOUS. The improved representation of upper atmosphere temperatures is driven by improved ozone concentrations near the tropopause and better upper level winds.

Description: A community diagnostic tool for chemistry climate model validation Geoscientific Model Development, 5, 1061-1073, 2012 Author(s): A. Gettelman, V. Eyring, C. Fischer, H. Shiona, I. Cionni, M. Neish, O. Morgenstern, S. W. Wood, and Z. Li This technical note presents an overview of the Chemistry-Climate Model Validation Diagnostic (CCMVal-Diag) tool for model evaluation. The CCMVal-Diag tool is a flexible and extensible open source package that facilitates the complex evaluation of global models. Models can be compared to other models, ensemble members (simulations with the same model), and/or many types of observations. The initial construction and application is to coupled chemistry-climate models (CCMs) participating in CCMVal, but the evaluation of climate models that submitted output to the Coupled Model Intercomparison Project (CMIP) is also possible. The package has been used to assist with analysis of simulations for the 2010 WMO/UNEP Scientific Ozone Assessment and the SPARC Report on the Evaluation of CCMs. The CCMVal-Diag tool is described and examples of how it functions are presented, along with links to detailed descriptions, instructions and source code. The CCMVal-Diag tool supports model development as well as quantifies model changes, both for different versions of individual models and for different generations of community-wide collections of models used in international assessments. The code allows further extensions by different users for different applications and types, e.g. to other components of the Earth system. User modifications are encouraged and easy to perform with minimum coding.

Description: Models of soil organic matter decomposition: the SoilR package, version 1.0 Geoscientific Model Development, 5, 1045-1060, 2012 Author(s): C. A. Sierra, M. Müller, and S. E. Trumbore Soil organic matter decomposition is a very important process within the Earth system because it controls the rates of mineralization of carbon and other biogeochemical elements, determining their flux to the atmosphere and the hydrosphere. SoilR is a modeling framework that contains a library of functions and tools for modeling soil organic matter decomposition under the R environment for computing. It implements a variety of model structures and tools to represent carbon storage and release from soil organic matter. In SoilR, organic matter decomposition is represented as a linear system of ordinary differential equations that generalizes the structure of most compartment-based decomposition models. A variety of functions is also available to represent environmental effects on decomposition rates. This document presents the conceptual basis for the functions implemented in the package. It is complementary to the help pages released with the software.

Description: Pre-industrial and mid-Pliocene simulations with NorESM-L: AGCM simulations Geoscientific Model Development, 5, 1033-1043, 2012 Author(s): Z. Zhang and Q. Yan In the Pliocene Model Intercomparison Project (PlioMIP), two sets of experiments are suggested. One includes a reference and a mid-Pliocene experiment run with atmosphere general circulation models (AGCM experiments, referred to as Experiments I), the other includes a pre-industrial and a mid-Pliocene experiment run with coupled ocean-atmosphere general circulation models (AOGCM experiments, referred to as Experiments II). In this paper, we describe the AGCM experiments with the atmosphere component in the low-resolution version of the Norwegian Earth System Model (NorESM-L), and also assess the potential uncertainties in analyzing mid-Pliocene climate anomalies that might result from the choice of the sea surface temperature (SST) forcing for the reference experiment (pre-industrial or modern). We carry out a mid-Pliocene experiment, a control experiment forced by the modern SST fields, and a pre-industrial experiment forced by the monthly SST fields from HadISST averaged between 1879 and 1900. Our experiments illustrate that the simulated mid-Pliocene global mean annual surface air temperature (SAT) is 17.1 °C. It is 2.5 °C warmer than the control experiment, but 2.7 °C warmer than the pre-industrial experiment. We find that the uncertainties in analyses of mid-Pliocene climate anomalies are small on a global scale, but still large on a regional scale. On the regional scale, these uncertainties should be noted and assessed in future PlioMIP studies.

Description: Review Article: "The Lagrangian description of aperiodic flows: a case study of the Kuroshio Current" Nonlinear Processes in Geophysics, 19, 449-472, 2012 Author(s): C. Mendoza and A. M. Mancho This article reviews several recently developed Lagrangian tools and shows how their combined use succeeds in obtaining a detailed description of purely advective transport events in general aperiodic flows. In particular, because of the climate impact of ocean transport processes, we illustrate a 2-D application on altimeter data sets over the area of the Kuroshio Current, although the proposed techniques are general and applicable to arbitrary time dependent aperiodic flows. The first challenge for describing transport in aperiodical time dependent flows is obtaining a representation of the phase portrait where the most relevant dynamical features may be identified. areas that are related to confinement regions. This representation is accomplished by using global Lagrangian descriptors that when applied for instance to the altimeter data sets retrieve over the ocean surface a phase portrait where the geometry of interconnected dynamical systems is visible. The phase portrait picture is essential because it evinces which transport routes are acting on the whole flow. Once these routes are roughly recognised, it is possible to complete a detailed description by the direct computation of the finite time stable and unstable manifolds of special hyperbolic trajectories that act as organising centres of the flow.