ALBERT

Description: Regional and local patterns in depth to water table, hydrochemistry and peat properties of bogs and their laggs in coastal British Columbia Hydrology and Earth System Sciences, 17, 3421-3435, 2013 Author(s): S. A. Howie and H. J. van Meerveld In restoration planning for damaged raised bogs, the lagg at the bog margin is often not given considerable weight and is sometimes disregarded entirely. However, the lagg is critical for the proper functioning of the bog, as it supports the water mound in the bog. In order to include the lagg in a restoration plan for a raised bog, it is necessary to understand the hydrological characteristics and functions of this rarely studied transition zone. We studied 13 coastal British Columbia (BC) bogs and identified two different gradients in depth to water table, hydrochemistry and peat properties: (1) a local bog expanse–bog margin gradient, and (2) a regional gradient related to climate and proximity to the ocean. Depth to water table generally increased across the transition from bog expanse to bog margin. In the bog expanse, pH was above 4.2 in the Pacific Oceanic wetland region (cooler and wetter climate) and below 4.3 in the Pacific Temperate wetland region (warmer and drier climate). Both pH and pH-corrected electrical conductivity increased significantly across the transition from bog expanse to bog margin, though not in all cases. Na + and Mg 2+ concentrations were generally highest in exposed, oceanic bogs and lower in inland bogs. Ash content in peat samples increased across the bog expanse–bog margin transition, and appears to be a useful abiotic indicator of the location of the bog margin. The observed variation in the hydrological and hydrochemical gradients across the bog expanse–bog margin transition highlights both local and regional diversity of bogs and their associated laggs.

Description: Imperfect scaling in distributions of radar-derived rainfall fields Hydrology and Earth System Sciences Discussions, 10, 11385-11422, 2013 Author(s): M. J. van den Berg, L. Delobbe, and N. E. C. Verhoest Fine scale rainfall observations for modeling exercises are often not available, but rather coarser data derived from a variety of sources are used. Effectively using these data sources in models often requires the probability distribution of the data at the applicable scale. Although numerous models for scaling distributions exist, these are often based on theoretical developments, rather than on data. In this study, we develop a model based on the α-stable distribution of rainfall fields, and tested on 5 min radar data from a Belgian weather radar. We use these data to estimate functions that describe parameters of the distribution over various scales. Moreover, we study how the mean of the distribution and the intermittency change with scale, and validate and design functions to describe the shape parameter of the distribution. This information was combined into an effective model of the distribution. Finally, the model was fitted to data from numerous storms, and the resulting parameters were compared to investigate the change in scaling behavior through time.

Description: Spatially resolved information on karst conduit flow from in-cave dye-tracing Hydrology and Earth System Sciences Discussions, 10, 11311-11335, 2013 Author(s): U. Lauber, W. Ufrecht, and N. Goldscheider Artificial tracers are powerful tools to investigate karst systems. Tracers are commonly injected into sinking streams or dolines, while springs serve as monitoring sites. The obtained flow and transport parameters represent mixed information from the vadose, epiphreatic and phreatic zones, i.e., the aquifer remains a black box. Accessible active caves constitute valuable but underexploited natural laboratories to gain detailed insights into the hydrologic functioning of the aquifer. Two multi-tracer tests in the catchment of a major karst spring (Blautopf, Germany) with injections and monitoring in two associated water caves aimed at obtaining spatially and temporally resolved information on groundwater flow in different compartments of the system. Two tracers were injected in the caves to characterize the hydraulic connections between them and with the spring. Two injections at the land surface, far from the spring, aimed at resolving the aquifer's internal drainage structure. Tracer breakthrough curves were monitored by field fluorimeters in caves and at the spring. Results demonstrate the dendritic drainage structure of the aquifer. It was possible to obtain relevant flow and transport parameters for different sections of this system. The highest mean flow velocities (275 m h −1 ) were observed in the near-spring epiphreatic section (open-channel flow), while velocities in the phreatic zone (pressurized flow) were one order of magnitude lower. Determined conduit water volumes confirm results of water balances and hydrograph analyses. In conclusion, experiments and monitoring in caves can deliver spatially resolved information on karst aquifer heterogeneity and dynamics that cannot be obtained by traditional investigative methods.

Description: On the lack of robustness of hydrologic models regarding water balance simulation – a diagnostic approach on 20 mountainous catchments using three models of increasing complexity Hydrology and Earth System Sciences Discussions, 10, 11337-11383, 2013 Author(s): L. Coron, V. Andréassian, C. Perrin, M. Bourqui, and F. Hendrickx This paper investigates the robustness of rainfall–runoff models when their parameters are transferred in time. More specifically, we studied their ability to simulate water balance on periods with different hydroclimatic characteristics. The testing procedure consisted in a series of parameter transfers between 10-yr periods and the systematic analysis of mean-volume errors. This procedure was applied to three conceptual models of different structural complexity over 20 mountainous catchments in southern France. The results showed that robustness problems are common. Errors on 10-yr-mean flows were significant for all three models and calibration periods, even when the entire record was used for calibration. Various graphical and numerical tools were used to show strong similarities between the shapes of mean flow biases calculated on a 10-yr-long sliding window when various parameter sets are used. Unexpected behavioural similarities were observed between the three models tested, considering their large differences in structural complexity. While the actual causes for robustness problems in these models remain unclear, this work stresses the limited transferability in time of the water balance adjustments made through parameter optimization. Although absolute differences between simulations obtained with different calibrated parameter sets were sometimes substantial, relative differences in simulated mean flows between time periods remained similar regardless of the calibrated parameter sets.

Description: Bayesian networks modelling in support to cross-cutting analysis of water supply and sanitation in developing countries Hydrology and Earth System Sciences, 17, 3397-3419, 2013 Author(s): C. Dondeynaz, J. López Puga, and C. Carmona Moreno Despite the efforts made towards the Millennium Development Goals targets during the last decade, improved access to water supply or basic sanitation still remains unavailable for millions of people across the world. This paper proposes a set of models that use 25 key variables and country profiles from the WatSan4Dev data set involving water supply and sanitation (Dondeynaz et al., 2012). This paper suggests the use of Bayesian network modelling methods because they are more easily adapted to deal with non-normal distributions, and integrate a qualitative approach for data analysis. They also offer the advantage of integrating preliminary knowledge into the probabilistic models. The statistical performance of the proposed models ranges between 20 and 5% error rates, which are very satisfactory taking into account the strong heterogeneity of variables. Probabilistic scenarios run from the models allow an assessment of the relationships between human development, external support, governance aspects, economic activities and water supply and sanitation (WSS) access. According to models proposed in this paper, gaining a strong poverty reduction will require the WSS access to reach 75–76% through: (1) the management of ongoing urbanisation processes to avoid slums development; and (2) the improvement of health care, for instance for children. Improving governance, such as institutional efficiency, capacities to make and apply rules, or control of corruption is positively associated with WSS sustainable development. The first condition for an increment of the HDP (human development and poverty) remains of course an improvement of the economic conditions with higher household incomes. Moreover, a significant country commitment to the environment, associated with civil society freedom of expression constitutes a favourable setting for sustainable WSS services delivery. Intensive agriculture using irrigation practises also appears as a mean for sustainable WSS thanks to multi-uses and complementarities. With a WSS sector organised at national level, irrigation practices can support the structuring and efficiency of the agriculture sector. It may then induce rural development in areas where WSS access often is set back compared to urban areas 1 . External financial support, called Official Development Assistance (ODA CI), plays a role in WSS improvement but comes last in the sensitivity analyses of models. An overall 47% of the Official Development Assistance goes first to poor countries, and is associated to governance aspects: (1) political stability and (2) country commitment to the environment and civil society degree of freedom. These governance aspects constitute a good framework for aid implementation in recipient countries. Modelling is run with the five groups of countries as defined in Dondeynaz et al. (2012). Models for profile 4 (essential external support) and profile 5 (primary material consumption) are specifically detailed and analysed in this paper. For countries in profile 4, fighting against water scarcity and progressing desertification should be the priority. However, for countries in profile 5, efforts should first concentrate on consolidation of political stability while supporting diversification of the economic activities. Nevertheless, for both profiles, reduction of poverty should remain the first priority as previously indicated. 1 JMP statistics, 2004 http://www.wssinfo.org/data-estimates/table/ , last access: 22 July 2013.

Description: Technical Note: A comparison of model and empirical measures of catchment-scale effective energy and mass transfer Hydrology and Earth System Sciences, 17, 3389-3395, 2013 Author(s): C. Rasmussen and E. L. Gallo Recent work suggests that a coupled effective energy and mass transfer (EEMT) term, which includes the energy associated with effective precipitation and primary production, may serve as a robust prediction parameter of critical zone structure and function. However, the models used to estimate EEMT have been solely based on long-term climatological data with little validation using direct empirical measures of energy, water, and carbon balances. Here we compare catchment-scale EEMT estimates generated using two distinct approaches: (1) EEMT modeled using the established methodology based on estimates of monthly effective precipitation and net primary production derived from climatological data, and (2) empirical catchment-scale EEMT estimated using data from 86 catchments of the Model Parameter Estimation Experiment (MOPEX) and MOD17A3 annual net primary production (NPP) product derived from Moderate Resolution Imaging Spectroradiometer (MODIS). Results indicated positive and significant linear correspondence ( R 2 = 0.75; P 〈 0.001) between model and empirical measures with an average root mean square error (RMSE) of 4.86 MJ m −2 yr −1 . Modeled EEMT values were consistently greater than empirical measures of EEMT. Empirical catchment estimates of the energy associated with effective precipitation ( E PPT ) were calculated using a mass balance approach that accounts for water losses to quick surface runoff not accounted for in the climatologically modeled E PPT . Similarly, local controls on primary production such as solar radiation and nutrient limitation were not explicitly included in the climatologically based estimates of energy associated with primary production ( E BIO ), whereas these were captured in the remotely sensed MODIS NPP data. These differences likely explain the greater estimate of modeled EEMT relative to the empirical measures. There was significant positive correlation between catchment aridity and the fraction of EEMT partitioned into E BIO ( F BIO ), with an increase in F BIO as a fraction of the total as aridity increases and percentage of catchment woody plant cover decreases. In summary, the data indicated strong correspondence between model and empirical measures of EEMT with limited bias that agree well with other empirical measures of catchment energy and water partitioning and plant cover.

Description: Analysis of projected changes in the occurrence of heat waves in Hungary Advances in Geosciences, 35, 115-122, 2013 Author(s): R. Pongrácz, J. Bartholy, and E. B. Bartha Heat wave events are important temperature-related hazards due to their impacts on human health. In 2004, a Heat Health Warning System including three levels of heat wave warning was developed on the basis of a retrospective analysis of mortality and meteorological data in Hungary to anticipate heat waves that may result in a large excess of mortality. Projected changes in the frequency of different heat wave warning levels are analysed for the 21st century. For this purpose, outputs of regional climate model PRECIS (Providing REgional Climates for Impacts Studies) are used taking into account three different global emissions scenarios (A2, A1B, B2). The results clearly show an increase in occurrence and length of heat waves with respect to the underlying emissions scenarios and regional climate model used. Moreover, the potential season of heat wave occurrences is projected to be lengthened by two months in 2071–2100 compared to 1961–1990.

Description: Review Article: On the relation between the seismic activity and the Hurst exponent of the geomagnetic field at the time of the 2000 Izu swarm Natural Hazards and Earth System Science, 13, 2189-2194, 2013 Author(s): F. Masci and J. N. Thomas Many papers document the observation of earthquake-related precursory signatures in geomagnetic field data. However, the significance of these findings is ambiguous because the authors did not adequately take into account that these signals could have been generated by other sources, and the seismogenic origin of these signals have not been validated by comparison with independent datasets. Thus, they are not reliable examples of magnetic disturbances induced by the seismic activity. Hayakawa et al. (2004) claim that at the time of the 2000 Izu swarm the Hurst exponent of the Ultra-Low-Frequency (ULF: 0.001–10 Hz) band of the geomagnetic field varied in accord with the energy released by the seismicity. The present paper demonstrates that the behaviour of the Hurst exponent was insufficiently investigated and also misinterpreted by the authors. We clearly show that during the Izu swarm the changes of the Hurst exponent were strongly related to the level of global geomagnetic activity and not to the increase of the local seismic activity.

Description: Forest fire danger rating in complex topography – results from a case study in the Bavarian Alps in autumn 2011 Natural Hazards and Earth System Science, 13, 2157-2167, 2013 Author(s): C. Schunk, C. Wastl, M. Leuchner, C. Schuster, and A. Menzel Forest fire danger rating based on sparse meteorological stations is known to be potentially misleading when assigned to larger areas of complex topography. This case study examines several fire danger indices based on data from two meteorological stations at different elevations during a major drought period. This drought was caused by a persistent high pressure system, inducing a pronounced temperature inversion and its associated thermal belt with much warmer, dryer conditions in intermediate elevations. Thus, a massive drying of fuels, leading to higher fire danger levels, and multiple fire occurrences at mid-slope positions were contrasted by moderate fire danger especially in the valleys. The ability of fire danger indices to resolve this situation was studied based on a comparison with the actual fire danger as determined from expert observations, fire occurrences and fuel moisture measurements. The results revealed that, during temperature inversion, differences in daily cycles of meteorological parameters influence fire danger and that these are not resolved by standard meteorological stations and fire danger indices (calculated on a once-a-day basis). Additional stations in higher locations or high-resolution meteorological models combined with fire danger indices accepting at least hourly input data may allow reasonable fire danger calculations under these circumstances.

Description: The Mediterranean Sea system: a review and an introduction to the special issue Ocean Science, 9, 789-803, 2013 Author(s): T. Tanhua, D. Hainbucher, K. Schroeder, V. Cardin, M. Álvarez, and G. Civitarese The Mediterranean Sea is a semi-enclosed sea characterized by high salinities, temperatures and densities. The net evaporation exceeds the precipitation, driving an anti-estuarine circulation through the Strait of Gibraltar, contributing to very low nutrient concentrations. The Mediterranean Sea has an active overturning circulation, one shallow cell that communicates directly with the Atlantic Ocean, and two deep overturning cells, one in each of the two main basins. It is surrounded by populated areas and is thus sensitive to anthropogenic forcing. Several dramatic changes in the oceanographic and biogeochemical conditions have been observed during the past several decades, emphasizing the need to better monitor and understand the changing conditions and their drivers. During 2011 three oceanographic cruises were conducted in a coordinated fashion in order to produce baseline data of important physical and biogeochemical parameters that can be compared to historic data and be used as reference for future observational campaigns. In this article we provide information on the Mediterranean Sea oceanographic situation, and present a short review that will serve as background information for the special issue in Ocean Science on "Physical, chemical and biological oceanography of the Mediterranean Sea". An important contribution of this article is the set of figures showing the large-scale distributions of physical and chemical properties along the full length of the Mediterranean Sea.

Description: Corrigendum to "NEMO on the shelf: assessment of the Iberia–Biscay–Ireland configuration" published in Ocean Sci., 9, 745–771, 2013 Ocean Science, 9, 787-787, 2013 Author(s): C. Maraldi, J. Chanut, B. Levier, N. Ayoub, P. De Mey, G. Reffray, F. Lyard, S. Cailleau, M. Drévillon, E. A. Fanjul, M. G. Sotillo, P. Marsaleix, and the Mercator Research and Development Team No abstract available.

Description: Preface "Landslide hazard and risk assessment at different scales" Natural Hazards and Earth System Science, 13, 2169-2171, 2013 Author(s): P. Reichenbach, A. Günther, and T. Glade

Description: We present a methodology for infrasonic remote sensing of winds in the stratosphere that does not require discrete ground-truth events. Our method uses measured time delays between arrays of sensors to provide group velocities (referred to here as celerities) and then minimizes the difference between observed and predicted celerities by perturbing an initial atmospheric specification. Because we focus on interarray propagation effects, it is not necessary to simulate the full propagation path from source to receiver. This feature allows us to use a relatively simple forward model that is applicable over short-regional distances. By focusing on stratospheric returns, we show that our non-linear inversion scheme converges much better if the starting model contains a strong stratospheric duct. Using the Horizontal Wind Model (HWM)/Mass Spectrometer Incoherent Scatter (MSISE) empirical climatology as a starting model, we demonstrate that the inversion scheme is robust to large uncertainties in backazimuth, but that uncertainties in the measured trace velocity and celerity require the use of prior constraints to ensure suitable convergence. The inversion of synthetic data, using realistic estimates of measurement error, shows that our scheme will nevertheless improve upon a starting model under most scenarios. The inversion scheme is applied to infrasound data recorded from a large event on 2010 December 25, which is presumed to be a bolide, using data from a nine-element infrasound network in Utah. We show that our recorded data require a stronger zonal wind speed in the stratosphere than is present in the HWM profile, and are more consistent with the Ground-to-Space (G2S) profile.

Description: Coda- Q is a stochastic parameter reflecting the heterogeneities of medium that seismic waves travel through. We confirmed that coda- Q would vary with the stress loaded to an elastic medium using numerical simulations of seismic wave propagation. When the stress is loaded, cracks in the crust could either close or newly open. The closure and opening of the cracks are not random but depending on the magnitude and the direction of the stress and the crack aspect ratio. The cracks in the medium after loading stress could be aligned in a specific orientation, and elastic wave velocity field would become anisotropic due to the alignment of specific crack orientations. Elastic wave velocity is in general faster along the direction corresponding with the crack orientation while slower along the perpendicular direction. In the numerical simulation, the effect of anisotropy in elastic wave velocity field due to the selective closure and opening of the cracks is calculated using a 2-D finite difference method assuming elastic wave velocity to be a function of the magnitude of loaded stress. The coda- Q calculated from seismic waves simulated for a model varies when the averaged normal stress changes. Our simulation indicated that the sensitivity of coda- Q –1 , that is the reciprocal of the coda- Q , would be 1.0 10 –2 (1.0 MPa –1 ) against the magnitude of the confining pressure and 1.0 10 –3 (1.0 deg –1 ) against the direction of principal stress. We would like to conclude that coda- Q , a stochastic parameter reflecting heterogeneities of subsurface medium, could become a quantitative state indicator of the stress field of the medium where seismic waves propagate through. Spatiotemporal variation of coda- Q reflects change in the stress field in the crust.

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: Influence of microphysical schemes on atmospheric water in the Weather Research and Forecasting model Geoscientific Model Development Discussions, 6, 4563-4601, 2013 Author(s): F. Cossu and K. Hocke This study examines how different microphysical parameterization schemes influence orographically-induced precipitation and the distributions of hydrometeors and water vapour for mid-latitude summer conditions in the Weather Research and Forecasting (WRF) model. A high-resolution two-dimensional idealized simulation is used to assess the differences between the schemes in which a moist air flow is interacting with a bell-shaped 2 km high mountain. Periodic lateral boundary conditions are chosen to recirculate atmospheric water in the domain. It is found that the 13 selected microphysical schemes conserve the water in the model domain. The gain or loss of water is less than 0.81% over a simulation time interval of 61 days. The differences of the microphysical schemes in terms of the distributions of water vapour, hydrometeors and accumulated precipitation are presented and discussed. The Kessler scheme, the only scheme without ice-phase processes, shows final values of cloud liquid water 14 times greater than the other schemes. The differences among the other schemes are not as extreme, but still they differ up to 79% in water vapour, up to 10 times in hydrometeors and up to 64% in accumulated precipitation at the end of the simulation. The microphysical schemes also differ in the surface evaporation rate. The WRF single-moment 3-class scheme has the highest surface evaporation rate compensated by the highest precipitation rate. The different distributions of hydrometeors and water vapour of the microphysical schemes induce differences up to 49 W m −2 in the downwelling shortwave radiation and up to 33 W m −2 in the downwelling longwave radiation.

Description: How do body-wave traveltimes constrain the Earth's radial (1-D) seismic structure? Existing 1-D seismological models underpin 3-D seismic tomography and earthquake location algorithms. It is therefore crucial to assess the quality of such 1-D models, yet quantifying uncertainties in seismological models is challenging and thus often ignored. Ideally, quality assessment should be an integral part of the inverse method. Our aim in this study is twofold: (i) we show how to solve a general Bayesian non-linear inverse problem and quantify model uncertainties, and (ii) we investigate the constraint on spherically symmetric P -wave velocity ( V P ) structure provided by body-wave traveltimes from the EHB bulletin (phases Pn , P , PP and PKP ). Our approach is based on artificial neural networks, which are very common in pattern recognition problems and can be used to approximate an arbitrary function. We use a Mixture Density Network to obtain 1-D marginal posterior probability density functions (pdfs), which provide a quantitative description of our knowledge on the individual Earth parameters. No linearization or model damping is required, which allows us to infer a model which is constrained purely by the data. We present 1-D marginal posterior pdfs for the 22 V P parameters and seven discontinuity depths in our model. P -wave velocities in the inner core, outer core and lower mantle are resolved well, with standard deviations of ~0.2 to 1 per cent with respect to the mean of the posterior pdfs. The maximum likelihoods of V P are in general similar to the corresponding ak135 values, which lie within one or two standard deviations from the posterior means, thus providing an independent validation of ak135 in this part of the radial model. Conversely, the data contain little or no information on P -wave velocity in the D '' layer, the upper mantle and the homogeneous crustal layers. Further, the data do not constrain the depth of the discontinuities in our model. Using additional phases available in the ISC bulletin, such as PcP , PKKP and the converted phases SP and ScP , may enhance the resolvability of these parameters. Finally, we show how the method can be extended to obtain a posterior pdf for a multidimensional model space. This enables us to investigate correlations between model parameters.

Description: We have developed a network optimization method for regional-scale microseismic monitoring networks and applied it to optimize the densification of the existing seismic network in northeastern Switzerland. The new network will build the backbone of a 10-yr study on the neotectonic activity of this area that will help to better constrain the seismic hazard imposed on nuclear power plants and waste repository sites. This task defined the requirements regarding location precision (0.5 km in epicentre and 2 km in source depth) and detection capability [magnitude of completeness M c  = 1.0 ( M L )]. The goal of the optimization was to find the geometry and size of the network that met these requirements. Existing stations in Switzerland, Germany and Austria were considered in the optimization procedure. We based the optimization on the simulated annealing approach proposed by Hardt & Scherbaum, which aims to minimize the volume of the error ellipsoid of the linearized earthquake location problem ( D -criterion). We have extended their algorithm to: calculate traveltimes of seismic body waves using a finite difference ray tracer and the 3-D velocity model of Switzerland, calculate seismic body-wave amplitudes at arbitrary stations assuming the Brune source model and using scaling and attenuation relations recently derived for Switzerland, and estimate the noise level at arbitrary locations within Switzerland using a first-order ambient seismic noise model based on 14 land-use classes defined by the EU-project CORINE and open GIS data. We calculated optimized geometries for networks with 10–35 added stations and tested the stability of the optimization result by repeated runs with changing initial conditions. Further, we estimated the attainable magnitude of completeness ( M c ) for the different sized optimal networks using the Bayesian Magnitude of Completeness (BMC) method introduced by Mignan et al. The algorithm developed in this study is also applicable to smaller optimization problems, for example, small local monitoring networks. Possible applications are volcano monitoring, the surveillance of induced seismicity associated with geotechnical operations and many more. Our algorithm is especially useful to optimize networks in populated areas with heterogeneous noise conditions and if complex velocity structures or existing stations have to be considered.

Description: Fluid injection in and withdrawal from wells are basic procedures in mining activities and deep resources exploitation, such as oil and gas extraction, permeability enhancement for geothermal exploitation and waste fluid disposal. All of these activities have the potential to induce seismicity, as exemplified by the 2006 Basel earthquake ( M L 3.4). Despite several decades of experience, the mechanisms of induced seismicity are not known in detail, which prevents effective risk assessment and/or mitigation. In this study, we provide an interpretation of induced seismicity based on computation of Coulomb stress changes that result from fluid injection/withdrawal at depth, mainly focused on the interpretation of induced seismicity due to stimulation of a geothermal reservoir. Seismicity is, theoretically, more likely where Coulomb stress changes are larger. For modeling purposes, we simulate the thermodynamic evolution of a system after fluid injection/withdrawal. The associated changes in pressure and temperature are subsequently considered as sources of incremental stress changes, which are then converted to Coulomb stress changes on favourably oriented faults, taking into account the background regional stress. Numerical results are applied to the water injection that was performed to create the fractured reservoir at the enhanced-geothermal-system site, Soultz-sous-Forets (France). Our approach describes well the observed seismicity, and provides an explanation for the different behaviors of a system when fluids are injected or withdrawn.

Description: We report on a broad-band high-resolution attenuation model for the North China Craton and surrounding regions based on regional Lg -wave data. Vertical broad-band waveforms recorded at 39 stations from 176 crustal earthquakes are collected to extract the Lg -wave amplitude spectra between 0.05 and 10.0 Hz. We use the dual-station method to generate a preliminary Q Lg model and use it as the initial model. Then, we combine the dual- and single-station data together to jointly invert the Q Lg distribution and Lg source excitation functions. These inversions are conducted independently at individual frequencies without using any a priori assumption about the frequency dependences in Q Lg and source terms. The maximum spatial resolution is approximately 1° x 1° in well-covered areas for frequencies between 0.05 and 2.0 Hz. The Q Lg image is then used to determine the relationship between the attenuation and different geological structures. Results show an average Q 0 (1 Hz Q Lg ) of 374 for the entire North China Craton with an increasing trend from east to west. Average Q 0 values are 337, 361 and 421 for the east, central and west blocks, respectively. For the surrounding regions, the Eastern Tibetan plateau has a very low Q 0 of 188, while the Northeast China Plate and the Tianshan–Xingmeng fold belts are characterized by high Q 0 values of 506 and 424, respectively. We also investigate regional variations of the Lg attenuation in low-frequency band between 0.2 and 1.0 Hz.

Description: We present a systematic study on the influence of pressure (0.1–600 MPa), temperature (750–1200 °C), carbon dioxide fugacity (log f CO 2  = –4.41 to 3.60) and time (2–12 hr) on the chemical and physical properties of carbonate rock. Our experiments aim to reproduce the conditions at the periphery of magma chamber where carbonate host rock is influenced by, but not readily assimilated by, magma. This permits the investigation of the natural conditions at which circulating fluids/gases promote infiltration reactions typical of metasomatic skarns that can involve large volumes of subvolcanic carbonate basements. Results show that, providing that carbon dioxide is retained in the pore space, decarbonation does not proceed at any magmatic pressure and temperature. However, when the carbon dioxide is free to escape, decarbonation can occur rapidly and is not hindered by a low initial porosity or permeability. Together with carbon dioxide and lime, portlandite, a mineral commonly found in voluminous metasomatic skarns, readily forms during carbonate decomposition. Post-experimental analyses highlight that thermal microcracking, a result of the highly anisotropic thermal expansion of calcite, exerts a greater influence on rock physical properties (porosity, ultrasonic wave velocities and elastic moduli) than decarbonation. Our data suggest that this will be especially true at the margins of dykes or magma bodies, where temperatures can reach up to 1200 °C. However, rock compressive strength is significantly reduced by both thermal cracking and decarbonation, explained by the relative weakness of lime + portlandite compared to calcite, and an increase in grain size with increasing temperature. Metasomatic skarns, whose petrogenetic reactions may involve a few tens of cubic kilometres, could therefore represent an important source of volcanic instability.

Description: The scattering of plane SH waves incident on a circular sectorial canyon is considered. An accurate region-matching technique is applied to derive a rigorous series solution. Appropriate wavefunctions are employed to describe antiplane motions. Judicious basis functions, involving Gegenbauer polynomials, are well utilized to correctly capture the singular behaviour in stress fields near the canyon bottom. The enforcement of matching conditions on the auxiliary boundary leads to the determination of unknown coefficients. Plotted results demonstrate the influence of pertinent parameters on surface and subsurface motions. Both steady-state and transient results are included. The solution technique proposed achieves a considerable reduction in the computational effort, facilitating benchmark computations. The derived series solution enriches the limited list of series solutions presently known for canyon problems related to SH -wave scattering.

Description: We present 3-D models of the P- and S -wave velocity distributions in the crust and uppermost mantle beneath Sicily, Calabria (Southern Italy), and surrounding submerged areas, obtained by tomographic inversion of traveltimes of regional body waves phases. Our method combines double-difference tomographic inversion with a post-processing procedure [Weighted Average Model method (WAM)]. This procedure was applied to a set of models consistent with the experimental data. We tested the ability of the WAM procedure to mitigate the uncertainty associated with the arbitrary nature of the many input parameters required for each inversion. The local reliability and resolution of the obtained models have been assessed through: synthetic tests, experimental tests carried out with independent data sets and unconventional tests based on the analysis of the internal consistency of the P - and S -velocity models. The tomographic images provide a detailed sketch of P- and S- wave velocity anomalies. These clearly show the shape of the Sicilian-Maghrebian belt beneath Sicily and Calabrian Arc at different depths. Low V P and Vs bodies are imaged beneath Stromboli and Marsili volcanoes in the southern Tyrrhenian, whereas high and low seismic velocities alternate beneath the Etna giving inferences on the possible depth of the mantle melting feeding the volcano. In the upper crust, the main sedimentary basins and tectonic features are also well imaged. Finally, tomographic cross sections show the trend of the Moho in the study area, where its depth ranges between 35 and 40 km beneath the Sicilian belt and between 15 and 22 km in the southern Tyrrhenian basin and Ionian Sea.

Description: We used two tracks of ALOS PALSAR images to investigate the focal mechanism and slip distribution of the 2011 March 24, M W 6.8 Burma strike-slip earthquake. Three different SAR techniques, namely conventional interferometry, SAR pixel offsets (SPO) and multiple-aperture InSAR (MAI), were employed to obtain the coseismic surface deformation fields along the ~30 km length of the fault rupture. Along-track measurements from SPO and MAI techniques show a high correlation, and were subsequently used to precisely determine the location and extent of the surface fault trace. The best-fitting fault model geometry derived from an iterative inversion technique suggests that the rupture occurred on a near-vertical sinistral strike-slip fault west of the Nam Ma fault with a strike of 70°. A maximum slip of 4.2 m occurs at a depth of 2.5 km, with significant slip constrained only to the upper 10 km of the crust.

Description: Knowledge of the mantle reflectivity structure is highly dependent on our ability to efficiently extract, and properly interpret, small seismic arrivals. Among the various data types and techniques, long-period SS/PP precursors and high-frequency receiver functions are routinely utilized to increase the confidence of the recovered mantle stratifications at distinct spatial scales. However, low resolution and a complex Fresnel zone are glaring weaknesses of SS precursors, while over-reliance on receiver distribution is a formidable challenge for the analysis of converted waves from oceanic regions. A promising high frequency alternative to receiver functions is P ' P ' precursors, which are capable of resolving mantle structures at vertical and lateral resolution of ~5 and ~200 km, respectively, owing to their spectral content, shallow angle of incidence and near-symmetric Fresnel zones. This study presents a novel processing method for both SS (or PP) and P ' P ' precursors based on deconvolution, stacking, Radon transform and depth migration. A suite of synthetic tests is performed to quantify the fidelity and stability of this method under different data conditions. Our multiresolution survey of the mantle at targeted areas near Nazca-South America subduction zone reveal both olivine and garnet related transitions at depths below 400 km. We attribute a depressed 660 to thermal variations, whereas compositional variations atop the upper-mantle transition zone are needed to explain the diminished or highly complex reflected/scattered signals from the 410 km discontinuity. We also observe prominent P ' P ' reflections within the transition zone, and the anomalous amplitudes near the plate boundary zone indicate a sharp (~10 km thick) transition that likely resonates with the frequency content of P ' P ' precursors. The migration of SS precursors in this study shows no evidence of split 660 reflections, but potential majorite–ilmenite (590–640 km) and ilmenite–perovskite transitions (740–750 km) are identified based on similarly processed high-frequency P ' P ' precursors. Additional findings of severely scattered energy in the lithosphere and distinct lower mantle reflections at ~800 km could be potentially important but require further verifications. Overall, our improved imaging methods and the strong sensitivity of P ' P ' precursors to the existence, depth, sharpness and strength of reflective structures offer significant future promise for the understanding of mantle mineralogy and dynamics.

Description: Time-dependent probabilistic seismic hazard assessment requires a stochastic description of earthquake occurrences. While short-term seismicity models are well-constrained by observations, the recurrences of characteristic on-fault earthquakes are only derived from theoretical considerations, uncertain palaeo-events or proxy data. Despite the involved uncertainties and complexity, simple statistical models for a quasi-period recurrence of on-fault events are implemented in seismic hazard assessments. To test the applicability of statistical models, such as the Brownian relaxation oscillator or the stress release model, we perform a systematic comparison with deterministic simulations based on rate- and state-dependent friction, high-resolution representations of fault systems and quasi-dynamic rupture propagation. For the specific fault network of the Lower Rhine Embayment, Germany, we run both stochastic and deterministic model simulations based on the same fault geometries and stress interactions. Our results indicate that the stochastic simulators are able to reproduce the first-order characteristics of the major earthquakes on isolated faults as well as for coupled faults with moderate stress interactions. However, we find that all tested statistical models fail to reproduce the characteristics of strongly coupled faults, because multisegment rupturing resulting from a spatiotemporally correlated stress field is underestimated in the stochastic simulators. Our results suggest that stochastic models have to be extended by multirupture probability distributions to provide more reliable results.

Description: Stress waves, known as acoustic emissions (AEs), are released by localized inelastic deformation events during the progressive failure of brittle rocks. Although several numerical models have been developed to simulate the deformation and damage processes of rocks, such as non-linear stress–strain behaviour and localization of failure, only a limited number have been capable of providing quantitative information regarding the associated seismicity. Moreover, the majority of these studies have adopted a pseudo-static approach based on elastic strain energy dissipation that completely disregards elastodynamic effects. This paper describes a new AE modelling technique based on the combined finite-discrete element method (FEM/DEM), a numerical tool that simulates material failure by explicitly considering fracture nucleation and propagation in the modelling domain. Given the explicit time integration scheme of the solver, stress wave propagation and the effect of radiated seismic energy can be directly captured. Quasi-dynamic seismic information is extracted from a FEM/DEM model with a newly developed algorithm based on the monitoring of internal variables (e.g. relative displacements and kinetic energy) in proximity to propagating cracks. The AE of a wing crack propagation model based on this algorithm are cross-analysed by traveltime inversion and energy estimation from seismic recordings. Results indicate a good correlation of AE initiation times and locations, and scaling of energies, independently calculated with the two methods. Finally, the modelling technique is validated by simulating a laboratory compression test on a granite sample. The micromechanical parameters of the heterogeneous model are first calibrated to reproduce the macroscopic stress–strain response measured during standard laboratory tests. Subsequently, AE frequency–magnitude statistics, spatial clustering of source locations and the evolution of AE rate are investigated. The distribution of event magnitude tends to decay as power law while the spatial distribution of sources exhibits a fractal character, in agreement with experimental observations. Moreover, the model can capture the decrease of seismic b value associated with the macrorupture of the rock sample and the transition of AE spatial distribution from diffuse, in the pre-peak stage, to strongly localized at the peak and post-peak stages, as reported in a number of published laboratory studies. In future studies, the validated FEM/DEM-AE modelling technique will be used to obtain further insights into the micromechanics of rock failure with potential applications ranging from laboratory-scale microcracking to engineering-scale processes (e.g. excavations within mines, tunnels and caverns, petroleum and geothermal reservoirs) to tectonic earthquakes triggering.

Description: We used strong-motion records from the 2012 May 20 and 29 Emilia-Romagna earthquakes ( M w 6.1 and 5.9, respectively) and four aftershocks with magnitudes ranging between 4.9 and 5.5 to analyse the S -wave spectral amplitude decay with distance and estimate acceleration source functions and site effects. The data set consists of six earthquakes, 44 stations and 248 records with hypocentral distances in the range 10 〈 r  〈 100 km. We rotated the accelerograms to calculate transverse and radial components of the acceleration spectrum. We found non-parametric attenuation functions that describe the spectral amplitude decay of SH and SV waves with distance at 60 different frequencies between 0.1 and 40 Hz. These attenuation functions provide an estimate of the quality factor Q at each frequency analysed. Assuming that geometrical spreading is 1/ r for r  ≤ r x and 1/( r x r ) 0.5 for r  〉 r x with r x  = 60 km and normalizing at 15 km (the recording distance where the attenuation functions start to decay), we find that the average Q for SH waves can be approximated by Q SH  = 82 ± 1 f  1.2±0.02 and by Q SV  = 79 ± 1 f  1.24±0.03 for SV waves in the frequency range 0.10 ≤ f  ≤ 10.7 Hz. At higher frequencies, 11.8 ≤ f  ≤ 40 Hz, the frequency dependence of Q weakens and is approximated by Q SH  = 301 ± 1 f   0.36±0.04 and Q SV  = 384 ± 1 f  0.28±0.04 . These results indicate that the S -wave attenuation is radially isotropic at local distances in the epicentral area. Nevertheless, we used these attenuation parameters separately to correct the radial (with Q SV ) and transverse (with Q SH ) components of the acceleration spectra and to separate source and site effects using a non-parametric spectral inversion scheme. We found that the source function of the main event and the bigger aftershocks show enhanced low frequency radiation between 0.4 and 3.0 Hz. We converted the source functions into far-field source acceleration spectra and interpreted the resulting source spectra in terms of Brune's model. The stress drops obtained range between approximately 0.9 and 2.9 MPa. Although all the recording stations used are located in the Po Plain, the site functions obtained from the spectral inversion show important amplification variability between the sites. We compared these site functions with the average horizontal to vertical spectral ratios calculated for each station, and we found consistent results for most stations.

Description: The sensitivity of global climate to the episodicity of fire aerosol emissions Atmospheric Chemistry and Physics Discussions, 13, 23691-23717, 2013 Author(s): S. K. Clark, D. S. Ward, and N. M. Mahowald One of the major ways in which forest and grass fires have an impact on global climate is through the release of aerosols. Most studies focusing on calculating the radiative forcing and other climate impacts of fire aerosols use monthly mean emissions derived from the Global Fire Emissions Database that captures only the seasonal cycle of fire aerosol emissions. Here we present the results of a sensitivity study that investigates the climate response to the episodicity of the fires, based on the standard approach which releases emissions every day, and contrasts that to the response when fires are represented as intense pulses of emissions that occur only over 1–2 days on a monthly, yearly, or five-yearly basis. Overall we find that in the modified cases with increased levels of episodicity, the all sky direct effect radiative forcing increases, the clear sky direct effect radiative forcing remains relatively constant, and the magnitude of the indirect effect radiative forcing decreases by about 1 W m −2 (from −1.6 to −0.6 W m −2 ). In the long term, we find that an increase in aerosol emission episodicity leads to an asymmetric change in indirect radiative forcing in the Northern Hemisphere compared to the Southern Hemisphere contributes to a slight shift in the annual average position of the intertropical convergence zone (ITCZ). This shift is found to have a mixed effect on the overall performance of the model at predicting precipitation rates in the tropics. Given these results we conclude that future studies that look to assess the present day global climate impacts of fire aerosols should consider the need to accurately represent fire episodicity.

Description: Not all calcite ballast is created equal: differing effects of foraminiferan and coccolith calcite on the formation and sinking of aggregates Biogeosciences Discussions, 10, 14861-14885, 2013 Author(s): K. Schmidt, C. L. De La Rocha, M. Gallinari, and G. Cortese Correlation between particulate organic carbon (POC) and calcium carbonate sinking through the deep ocean has led to the idea that ballast provided by calcium carbonate is important for the export of POC from the surface ocean. While this idea is certainly to some extent true, it is worth considering in more nuance, for example, examining the different effects on the aggregation and sinking of POC of small, non-sinking calcite particles like coccoliths and large, rapidly sinking calcite like planktonic foraminiferan tests. We have done that here in a simple experiment carried out in roller tanks that allow particles to sink continuously without being impeded by container walls. Coccoliths were efficiently incorporated into aggregates that formed during the experiment, increasing their sinking speed compared to similarly sized aggregates lacking added calcite ballast. The foraminiferan tests, which sank as fast as 700 m d −1 , became associated with only very minor amounts of POC. In addition, when they collided with other, larger, foraminferan-less aggregates, they fragmented them into two smaller, more slowly sinking aggregates. While these effects were certainly exaggerated within the confines of the roller tanks, they clearly demonstrate that calcium carbonate ballast is not just calcium carbonate ballast- different forms of calcium carbonate ballast have notably different effects on POC aggregation, sinking, and export.

Description: Laboratory studies of immersion and deposition mode ice nucleation of ozone aged mineral dust particles Atmospheric Chemistry and Physics, 13, 9097-9118, 2013 Author(s): Z. A. Kanji, A. Welti, C. Chou, O. Stetzer, and U. Lohmann Ice nucleation in the atmosphere is central to the understanding the microphysical properties of mixed-phase and cirrus clouds. Ambient conditions such as temperature ( T ) and relative humidity (RH), as well as aerosol properties such as chemical composition and mixing state play an important role in predicting ice formation in the troposphere. Previous field studies have reported the absence of sulfate and organic compounds on mineral dust ice crystal residuals sampled at mountain top stations or aircraft based measurements despite the long-range transport mineral dust is subjected to. We present laboratory studies of ice nucleation for immersion and deposition mode on ozone aged mineral dust particles for 233 〈 T 〈 263 K. Heterogeneous ice nucleation of untreated kaolinite (Ka) and Arizona Test Dust (ATD) particles is compared to corresponding aged particles that are subjected to ozone concentrations of 0.4–4.3 ppmv in a stainless steel aerosol tank. The portable ice nucleation counter (PINC) and immersion chamber combined with the Zurich ice nucleation chamber (IMCA-ZINC) are used to conduct deposition and immersion mode measurements, respectively. Ice active fractions as well as ice active surface site densities ( n s ) are reported and observed to increase as a function of decreasing temperature. We present first results that demonstrate enhancement of the ice nucleation ability of aged mineral dust particles in both the deposition and immersion mode due to ageing. We also present the first results to show a suppression of heterogeneous ice nucleation activity without the condensation of a coating of (in)organic material. In immersion mode, low ozone exposed Ka particles showed enhanced ice activity requiring a median freezing temperature of 1.5 K warmer than that of untreated Ka, whereas high ozone exposed ATD particles showed suppressed ice nucleation requiring a median freezing temperature of 3 K colder than that of untreated ATD. In deposition mode, low exposure Ka had ice active fractions of an order of magnitude higher than untreated Ka, whereas high ozone exposed ATD had ice active fractions up to a factor of 4 lower than untreated ATD. From our results, we derive and present parameterizations in terms of n s ( T ) that can be used in models to predict ice nuclei concentrations based on available aerosol surface area.

Description: Antecedent flow conditions and nitrate concentrations in the Mississippi River Basin Hydrology and Earth System Sciences Discussions, 10, 11451-11484, 2013 Author(s): J. C. Murphy, R. M. Hirsch, and L. A. Sprague The influence of antecedent flow conditions on nitrate concentrations was explored at eight sites in the Mississippi River Basin, USA. Antecedent moisture conditions have been shown to influence nutrient export from small, relatively homogenous basins, but this influence has not been observed at a regional or continental scale. Antecedent flow conditions were quantified as the ratio between the mean daily flow of the previous year and the mean daily flow from the period of record ( Q ratio), and the Q ratio was statistically related to nitrate anomalies (the unexplained variability in nitrate concentration after filtering out season, long-term trend, and contemporaneous flow effects) at each site. Nitrate anomaly and Q ratio were negatively related at three of the four major tributary sites and upstream in the Mississippi River, indicating that when the previous year was drier than average, at these sites, nitrate concentrations were higher than expected. The strength of these relationships increased when data were subdivided by contemporaneous flow conditions. Five of the eight sites had significant negative relationships ( p ≤ 0.05) at high or moderately high contemporaneous flows, suggesting nitrate that accumulates in these basins during a drought is flushed during subsequent storm events. At half of the sites, when flow during the previous year was 50% drier than average, nitrate concentration can be from 9 and 27% higher than nitrate concentrations that follow a year with average daily flow. Conversely, nitrate concentration can be from 8 and 21% lower than expected when the previous year was 50% wetter than average. These relationships between nitrate concentration and Q ratio serve as the basis for future studies that can better define specific hydrologic processes occurring during and after a drought, which influence nitrate concentration, such as the duration or magnitude of low flows, and the timing of low and high flows.

Description: Tidally induced velocity variations of the Beardmore Glacier, Antarctica, and their representation in satellite measurements of ice velocity The Cryosphere, 7, 1375-1384, 2013 Author(s): O. J. Marsh, W. Rack, D. Floricioiu, N. R. Golledge, and W. Lawson Ocean tides close to the grounding line of outlet glaciers around Antarctica have been shown to directly influence ice velocity, both linearly and non-linearly. These fluctuations can be significant and have the potential to affect satellite measurements of ice discharge, which assume displacement between satellite passes to be consistent and representative of annual means. Satellite observations of horizontal velocity variation in the grounding zone are also contaminated by vertical tidal effects, the importance of which is highlighted here in speckle tracking measurements. Eight TerraSAR-X scenes from the grounding zone of the Beardmore Glacier are analysed in conjunction with GPS measurements to determine short-term and decadal trends in ice velocity. Diurnal tides produce horizontal velocity fluctuations of 〉50% on the ice shelf, recorded in the GPS data 4 km downstream of the grounding line. This variability decreases rapidly to 〈5% only 15 km upstream of the grounding line. Daily fluctuations are smoothed to 〈1% in the 11-day repeat pass TerraSAR-X imagery, but fortnightly variations over this period are still visible and show that satellite-velocity measurements can be affected by tides over longer periods. The measured tidal displacement observed in radar look direction over floating ice also allows the grounding line to be identified, using differential speckle tracking where phase information cannot be easily unwrapped.

Description: A new data set of soil mineralogy for dust-cycle modeling Atmospheric Chemistry and Physics Discussions, 13, 23943-23993, 2013 Author(s): E. Journet, Y. Balkanski, and S. P. Harrison The mineralogy of airborne dust affects the impact of dust particles on direct and indirect radiative forcing, on atmospheric chemistry and on biogeochemical cycling. It is determined partly by the mineralogy of the dust-source regions and partly by size-dependent fractionation during erosion and transport. Here we present a data set that characterizes the clay and silt sized fractions of global soil units in terms of the abundance of 12 minerals that are important for dust-climate interactions: quartz, feldspars, illite, smectite, kaolinite, chlorite, vermiculite, mica, calcite, gypsum, hematite and goethite. The basic mineralogical information is derived from the literature, and is then expanded following explicit rules, in order to characterize as many soil units as possible. We present three alternative realisations of the mineralogical maps that account for the uncertainties in the mineralogical data. We examine the implications of the new database for calculations of the single scattering albedo of airborne dust and thus for dust radiative forcing.

Description: Calibration of the passive cavity aerosol spectrometer probe for airborne determination of the size distribution Atmospheric Measurement Techniques, 6, 2349-2358, 2013 Author(s): Y. Cai, J. R. Snider, and P. Wechsler This work describes calibration methods for the particle sizing and particle concentration systems of the passive cavity aerosol spectrometer probe (PCASP). Laboratory calibrations conducted over six years, in support of the deployment of a PCASP on a cloud physics research aircraft, are analyzed. Instead of using the many calibration sizes recommended by the PCASP manufacturer, a relationship between particle diameter and scattered light intensity is established using three sizes of mobility-selected polystyrene latex particles, one for each amplifier gain stage. In addition, studies of two factors influencing the PCASP's determination of the particle size distribution – amplifier baseline and particle shape – are conducted. It is shown that the PCASP-derived size distribution is sensitive to adjustments of the sizing system's baseline voltage, and that for aggregates of spheres, a PCASP-derived particle size and a sphere-equivalent particle size agree within uncertainty dictated by the PCASP's sizing resolution. Robust determinations of aerosol concentration, and size distribution, also require calibration of the PCASP's aerosol flowrate sensor. Sensor calibrations, calibration drift, and the sensor's non-linear response are documented.

Description: Contribution of snow and glacier melt to discharge for highly glacierised catchments in Norway Hydrology and Earth System Sciences Discussions, 10, 11485-11517, 2013 Author(s): M. Engelhardt, T. V. Schuler, and L. M. Andreassen Glacierised catchments significantly alter the streamflow regime due to snow and glacier meltwater contribution to discharge. In this study, we modelled the mass balance and discharge rates for three highly glacierised catchments (〉50% glacier cover) in western Norway over the period 1961–2012. The spatial pattern of the catchments follows a gradient in climate continentality from west to east. The model uses gridded temperature and precipitation values from seNorge ( http://senorge.no ) as input which are available at a daily resolution. It accounts for accumulation of snow, transformation of snow to firn and ice, evaporation and melt. The model was calibrated for each catchment based on measurements of seasonal glacier mass-balances and daily discharge rates. For validation, daily melt rates were compared with measurements from sonic rangers located in the ablation zones of two of the glaciers and an uncertainty analysis was performed for the third catchment. The discharge contributions from snowmelt, glacier melt and rain were analysed with respect to spatial variations and temporal evolution. The model simulations reveal an increase of the relative contribution from glacier melt for the three catchments from less than 10% in the early 1990s to 15–30% in the late 2000s. The decline in precipitation by 10–20% in the same period was therefore overcompensated resulting in an increase of the annual discharge by 5–20%. Annual discharge sums and annual glacier melt are strongest correlated with annual and winter precipitation at the most maritime glacier and, with increased climate continentality, variations in both glacier melt contribution and annual discharge are becoming stronger correlated with variations in summer temperatures.

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: Decadal trends in the Antarctic sea ice extent ultimately controlled by ice-ocean feedback The Cryosphere Discussions, 7, 4585-4632, 2013 Author(s): H. Goosse and V. Zunz The large natural variability of the Antarctic sea ice is a key characteristic of the system that might be responsible for the small positive trend in sea ice extent observed since 1979. In order to gain insight in the processes responsible for this variability, we have analysed in a control simulation performed with a coupled climate model a strong positive ice-ocean feedback that amplifies sea ice variations. When sea ice concentration increases in a region, in particular close to the ice edge, the mixed layer depth tends to decrease. This can be caused by a net inflow of ice and thus of freshwater that stabilizes the water column. Another stabilizing mechanism at interannual time scales that appears more widespread in our simulation is associated with the downward salt transport due to the seasonal cycle of ice formation: brine is released in winter when ice is formed and mixed over a deep layer while the freshwater flux caused by ice melting is included in a shallow layer, resulting in a net vertical transport of salt. Because of this stronger stratification due to the presence of sea ice, more heat is stored at depth in the ocean and the vertical oceanic heat flux is reduced, which contributes to maintain a higher ice extent. This positive feedback is not associated with a particular spatial pattern. Consequently, the spatial distribution of the trend in ice concentration is largely imposed by the wind changes that can provide the initial perturbation. A positive freshwater flux could alternatively be the initial trigger but the amplitude of the final response of the sea ice extent is finally set up by the amplification related to ice-ocean feedback. Initial conditions have also an influence as the chance to have a large increase in ice extent is higher if starting from a state characterized by a low value.

Description: Growth of sulphuric acid nanoparticles under wet and dry conditions Atmospheric Chemistry and Physics Discussions, 13, 24087-24125, 2013 Author(s): L. Škrabalová, D. Brus, T. Anttila, V. Ždímal, and H. Lihavainen New particle formation, which greatly influences the number concentrations and size distributions of an atmospheric aerosol, is often followed by a rapid growth of freshly formed particles. The initial growth of a newly formed aerosol is the crucial process determining the fraction of nucleated particles growing into cloud condensation nuclei sizes, which have a significant influence on climate. In this study, we report the laboratory observations of the growth of nanoparticles produced by nucleation of H 2 SO 4 and water in a laminar flow tube at temperatures of 283, 293 and 303 K, under dry (a relative humidity of 1%) and wet conditions (relative humidity of 30%) and residence times of 30, 45, 60 and 90 s. The initial H 2 SO 4 concentration spans the range from 2 × 10 8 to 1.4 × 10 10 molecule cm −3 and the calculated wall losses of H 2 SO 4 were assumed to be diffusion limited. The detected particle number concentrations, measured by the Ultrafine Condensation Particle Counter (UCPC) and Differential Mobility Particle Sizer (DMPS), were found to depend strongly on the residence time. Hygroscopic particle growth, presented by growth factors, was found to be in good agreement with the previously reported studies. The experimental growth rates ranged from 20 nm h −1 to 890 nm h −1 at RH 1% and from 7 nm h −1 to 980 nm h −1 at RH 30% and were found to increase significantly with the increasing concentration of H 2 SO 4 . Increases in the nucleation temperature had a slight enhancing effect on the growth rates under dry conditions. The influence of relative humidity on growth was not consistent – at lower H 2 SO 4 concentrations, the growth rates were higher under dry conditions while at H 2 SO 4 concentrations greater than 1×10 9 molecule cm −3 the growth rates were higher under wet conditions. The growth rates show only a weak dependence on the residence time. The experimental observations were compared with predictions made using a numerical model, which investigates the growth of particles with three different extents of neutralization by the ammonia NH 3 : (1) pure H 2 SO 4 – H 2 O particles (2) particles formed by ammonium bisulphate, (NH 4 )HSO 4 (3) particles formed by ammonium sulphate, (NH 4 ) 2 SO 4 . The highest growth rates were found for ammonium sulphate particles. Since the model accounting for the initial H 2 SO 4 concentration predicted the experimental growth rates correctly, our results suggest that the commonly presumed diffusional wall losses of H 2 SO 4 are not so significant. We therefore assume that there are not only losses of H 2 SO 4 on the wall but also a flux of H 2 SO 4 molecules from the wall into the flow tube, the effect being more profound under dry conditions and at higher temperatures of the tube wall. Based on a comparison with the atmospheric observations, our results indicate that sulphuric acid alone can not explain the growth rates of particles formed in the atmosphere.

Description: Experimental quantification of contact freezing in an electrodynamic balance Atmospheric Measurement Techniques, 6, 2373-2382, 2013 Author(s): N. Hoffmann, A. Kiselev, D. Rzesanke, D. Duft, and T. Leisner Heterogeneous nucleation of ice in a supercooled water droplet induced by external contact with a dry aerosol particle has long been known to be more effective than freezing induced by the same nucleus immersed in the droplet. However, the experimental quantification of contact freezing is challenging. Here we report an experimental method to determine the temperature-dependent ice nucleation probability of size-selected aerosol particles. The method is based on the suspension of supercooled charged water droplets in a laminar flow of air containing aerosol particles as contact freezing nuclei. The rate of droplet–particle collisions is calculated numerically with account for Coulomb attraction, drag force and induced dipole interaction between charged droplet and aerosol particles. The calculation is verified by direct counting of aerosol particles collected by a levitated droplet. By repeating the experiment on individual droplets for a sufficient number of times, we are able to reproduce the statistical freezing behavior of a large ensemble of supercooled droplets and measure the average rate of freezing events. The freezing rate is equal to the product of the droplet–particle collision rate and the probability of freezing on a single contact, the latter being a function of temperature, size and composition of the contact ice nuclei. Based on these observations, we show that for the types of particles investigated so far, contact freezing is the dominating freezing mechanism on the timescale of our experiment.

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: Contribution of directly connected and isolated impervious areas to urban drainage network hydrographs Hydrology and Earth System Sciences, 17, 3473-3483, 2013 Author(s): Y. Seo, N.-J. Choi, and A. R. Schmidt This paper addresses the mass balance error observed in runoff hydrographs in urban watersheds by introducing assumptions regarding the contribution of infiltrated rainfall from pervious areas and isolated impervious area (IIA) to the runoff hydrograph. Rainfall infiltrating into pervious areas has been assumed not to contribute to the runoff hydrograph until Hortonian excess rainfall occurs. However, mass balance analysis in an urban watershed indicates that rainfall infiltrated to pervious areas can contribute directly to the runoff hydrograph, thereby offering an explanation for the long hydrograph tail commonly observed in runoff from urban storm sewers. In this study, a hydrologic analysis based on the width function is introduced, with two types of width functions obtained from both pervious and impervious areas, respectively. The width function can be regarded as the direct interpretation of the network response. These two width functions are derived to obtain distinct response functions for directly connected impervious areas (DCIA), IIA, and pervious areas. The results show significant improvement in the estimation of runoff hydrographs and suggest the need to consider the flow contribution from pervious areas to the runoff hydrograph. It also implies that additional contribution from flow paths through joints and cracks in sewer pipes needs to be taken into account to improve the estimation of runoff hydrographs in urban catchments.

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: Phenology as a strategy for carbon optimality: a global model Biogeosciences Discussions, 10, 15107-15152, 2013 Author(s): S. Caldararu, D. W. Purves, and P. I. Palmer Phenology is essential to our understanding of biogeochemical cycles and the climate system. We develop a global mechanistic model of leaf phenology based on the hypothesis that phenology is a strategy for optimal carbon gain at the canopy level so that trees adjust leaf gains and losses in response to environmental factors such as light, temperature and soil moisture, to achieve maximum carbon assimilation. We fit this model to five years of satellite observations of leaf area index (LAI) using a Bayesian fitting algorithm. We show that our model is able to reproduce phenological patterns for all vegetation types and use it to explore variations in growing season length and the climate factors that limit leaf growth for different biomes. Phenology in wet tropical areas is limited by leaf age physiological constraints while at higher latitude leaf seasonality is limited by low temperature and light availability. Leaf growth in grassland regions is limited by water availability but often in combination with other factors. This model will advance the current understanding of phenology for ecosystem carbon models and our ability to predict future phenological behaviour.

Description: Bias correction can modify climate model-simulated precipitation changes without adverse affect on the ensemble mean Hydrology and Earth System Sciences Discussions, 10, 11585-11611, 2013 Author(s): E. P. Maurer and D. W. Pierce When applied to remove climate model biases in precipitation, quantile mapping can in some settings modify the simulated trends. This has important implications when the precipitation will be used to drive an impacts model that is sensitive to changes in precipitation. We use daily precipitation output from 12 general circulation models (GCMs) over the conterminous United States interpolated to a common 1° grid, and gridded observations aggregated to the same scale, to compare precipitation differences before and after quantile mapping bias correction. The change in seasonal mean (winter, DJF, and summer, JJA) precipitation between different 30-yr historical periods is compared to examine (1) the consensus among GCMs as to whether the bias correction tends to amplify or diminish their simulated precipitation trends, and (2) whether the modification of the change in precipitation tends to improve or degrade the correspondence to observed changes in precipitation for the same periods. In some cases, for a particular GCM, the trend modification can be as large as the original simulated change, though the areas where this occurs varies among GCMs so the ensemble median shows smaller trend modification. In specific locations and seasons the trend modification by quantile mapping improves correspondence with observed trends, and in others it degrades it. In the majority of the domain the ensemble median is for little effect on the correspondence of simulated precipitation trends with observed. This highlights the need to use an ensemble of GCMs rather than relying on a small number of models to estimate impacts.

Description: Linking ICT and society in early warning and adaptation to hydrological extremes in mountains Natural Hazards and Earth System Science, 13, 2253-2270, 2013 Author(s): C. de Jong The assessment of the societal impact of hydrological extremes is particularly important in mountain regions, since mountains can be considered both as the generators and victims of extreme events. ICT can provide a powerful tool for transmitting hydro-meteorological information to predict, prepare and adapt to such events. However, in remote regions, such as mountains, the poles, deserts and islands, preventive and adaptive measures are often restricted by data availability and lack and/or incoherence of data networks. This paper distinguishes between early warning of floods and droughts, emphasising the latter in particular in mountains and explores the possibilities of enhancing the role of society in data collection, the identification, activation and application of stakeholder knowledge and transferral of data from gauged to ungauged catchments.

Description: The dihedral angle formed at junctions between two plagioclase grains and a grain of augite is only very rarely in textural equilibrium in gabbros from kilometre-scale crustal layered intrusions. The median of a population of these disequilibrium angles, cpp , varies systematically within a single layered intrusion, remaining constant over large stretches of stratigraphy with significant increases and decreases associated with the addition or reduction respectively of the number of phases on the liquidus of the bulk magma. The stepwise changes in cpp are present in the Upper Zone of the Bushveld Complex, the Megacyclic Unit I of the Sept Iles Intrusion, and the Layered Series of the Skaergaard intrusion. The plagioclase-bearing cumulates of Rum have a bimodal distribution of cpp , dependent on whether the cumulus assemblage includes clinopyroxene. The presence of the stepwise changes is independent of the order of arrival of cumulus phases and of the composition of either the cumulus phases or the inferred composition of the interstitial liquid. The only parameter that behaves in an exactly analogous manner to cpp is the rate of change in enthalpy with temperature ( H / T ) during crystallization. Both H / T and cpp increase with the addition of a liquidus phase, and decrease with the removal of a liquidus phase. The replacement of one phase by another has little effect on H / T and no discernible effect on cpp . An increase of H / T results in an increase in the fraction of the total enthalpy budget that is the latent heat of crystallization (the fractional latent heat). It also increases the mass crystallized in each incremental temperature drop (the crystal productivity). These increases of both fractional latent heat and crystal productivity are likely to cause an increase in the time taken to form three-grain junctions in the mush via thermal buffering of a thickened mushy layer. We suggest these are the underlying causes of stepwise increases in cpp . Stepwise changes in the geometry of three-grain junctions in fully solidified gabbros thus provide a clear microstructural marker for the progress of fractionation down the liquid line of descent in layered intrusions.

Description: High-resolution sampling in monogenetic fields has the potential to reveal fine-scale heterogeneity of the mantle, a feature that may be overwhelmed by larger fluxes of magma, or missed by under-sampling. The Quaternary Auckland Volcanic Field (AVF) in northern New Zealand is a basaltic field of 51 small-volume volcanic centres, and is one of the best-sampled examples of a monogenetic volcanic field. We present data for 12 centres in the volcanic field. These show the large compositional variations between volcanoes as well as through single eruptive sequences. Whole-rock compositions range from subalkaline basalt in the larger centres, through alkali basalt to nephelinite in the smallest centres. Fractional crystallization has had a limited effect in many of the centres, but high-pressure clinopyroxene crystallization may have occurred in others. Three end-members are observed in Pb isotope space, indicating that distinct mantle source components are involved in the petrogenesis of the magmas. Whole-rock multi-element patterns show that the larger centres have prominent positive Sr anomalies and lack K anomalies, whereas the smaller centres have prominent negative K anomalies and lack Sr anomalies. The melting parameters and compositions of the sources involved are modelled using trace element ratios and multi-element patterns, and three components are characterized: (1) fertile peridotite with a Pb-isotope composition similar to Pacific mid-ocean ridge basalt; (2) eclogite domains with a HIMU-like isotope composition dispersed within the fertile peridotite; (3) slightly depleted subduction-metasomatized peridotitic lithospheric mantle (containing c . 3% subduction fluids). Modelling shows that melting in the AVF begins in garnet-bearing fertile asthenosphere (with preferential melting of eclogite domains) and that melts are variably diluted by melts of the lithospheric source. The U–Th isotope compositions of the end-members in the AVF show 230 Th excess [( 230 Th/ 232 Th) ratios of 1·11–1·38], with the samples of lower ( 230 Th/ 232 Th) exhibiting higher ( 238 U/ 232 Th), which we attribute to the dilution effect of the melts from the lithospheric mantle source. Modelling reveals a correlation between melting in the asthenosphere, the degree of melting and incorporation of the metasomatized lithospheric mantle source, and the resultant size of the volcanic centre. This suggests that the scale of the eruption may essentially be controlled by asthenospheric mantle dynamics.

Description: The origin of mafic and ultramafic sills exhibiting different whole-rock compositional profiles (e.g. I-, C-, D-, M- and S-shaped profiles) remains controversial. We have addressed this issue by revisiting three ~100 m thick Siberian dolerite sills (Vavukansky, Kuz’movsky and Vilyuysky) that display remarkable internal differentiation. The Vavukansky sill has an M-shaped profile with prominent basal and top reversals showing inward increases in whole-rock MgO, Mg-number [100Mg/(Mg + Fe)] and normative An content [100An/(An + Ab)], followed by the Layered and Upper Border Series with inward decreases in these indices. The Kuz’movsky and Vilyuysky sills both show S-shaped profiles similar to the Vavukansky sill, but lack a top reversal. These whole-rock M- and S-shaped profiles are accompanied by similar profiles in mineral compositions. Plagioclase and, to a lesser extent, olivine show systematic inward increases in An content and Mg-number, respectively, across basal and top reversals. These compositional trends are followed by inward decreases in these ratios in the interiors of the Vavukansky and Kuz’movsky sills. Currently accepted models attribute whole-rock M- and S-shaped compositional profiles to crystal settling, compositional convection or compaction operating in closed systems. Our observations challenge these traditional interpretations because variations in mineral compositions observed in marginal reversals cannot result from closed-system fractionation. We suggest instead that initially the sills evolved as open systems that were slowly inflated by magmas that became gradually more primitive with time. The inflation was accompanied by in situ crystallization that preserved the preceding fractionation history of the injected magmas by forming basal and top reversals with minerals becoming more primitive inwards. This process culminated with rapid inflation of the sills to their current size owing to a major influx of primitive magma. Subsequently, magma flow through the sills ceased and they evolved as closed systems by fractional crystallization. This resulted in the Layered and Upper Border Series with minerals becoming more evolved inwards. This model can be extended to explain other compositional profiles and petrological features in mafic and ultramafic sills.

Description: A plume-in-grid approach to characterize air quality impacts of aircraft emissions at the Hartsfield–Jackson Atlanta International Airport Atmospheric Chemistry and Physics, 13, 9285-9302, 2013 Author(s): J. Rissman, S. Arunachalam, M. Woody, J. J. West, T. BenDor, and F. S. Binkowski This study examined the impacts of aircraft emissions during the landing and takeoff cycle on PM 2.5 concentrations during the months of June and July 2002 at the Hartsfield–Jackson Atlanta International Airport. Primary and secondary pollutants were modeled using the Advanced Modeling System for Transport, Emissions, Reactions, and Deposition of Atmospheric Matter (AMSTERDAM). AMSTERDAM is a modified version of the Community Multiscale Air Quality (CMAQ) model that incorporates a plume-in-grid process to simulate emissions sources of interest at a finer scale than can be achieved using CMAQ's model grid. Three fundamental issues were investigated: the effects of aircraft on PM 2.5 concentrations throughout northern Georgia, the differences resulting from use of AMSTERDAM's plume-in-grid process rather than a traditional CMAQ simulation, and the concentrations observed in aircraft plumes at subgrid scales. Comparison of model results with an air quality monitor located in the vicinity of the airport found that normalized mean bias ranges from −77.5% to 6.2% and normalized mean error ranges from 40.4% to 77.5%, varying by species. Aircraft influence average PM 2.5 concentrations by up to 0.232 μg m −3 near the airport and by 0.001–0.007 μg m −3 throughout the Atlanta metro area. The plume-in-grid process increases concentrations of secondary PM pollutants by 0.005–0.020 μg m −3 (compared to the traditional grid-based treatment) but reduces the concentration of non-reactive primary PM pollutants by up to 0.010 μg m −3 , with changes concentrated near the airport. Examination of subgrid-scale results indicates that median aircraft contribution to grid cells is higher than median puff concentration in the airport's grid cell and outside of a 20 km × 20 km square area centered on the airport, while in a 12 km × 12 km square ring centered on the airport, puffs have median concentrations over an order of magnitude higher than aircraft contribution to the grid cells. Maximum puff impacts are seen within the 12 km × 12 km ring, not in the airport's own grid cell, while maximum grid cell impacts occur within the airport's grid cell. Twenty-one (21)% of all aircraft-related puffs from the Atlanta airport have at least 0.1 μg m −3 PM 2.5 concentrations. Near the airport, median daily puff concentrations vary between 0.017 and 0.134 μg m −3 (0.05 and 0.35 μg m −3 at ground level), while maximum daily puff concentrations vary between 6.1 and 42.1 μg m −3 (7.5 and 42.1 μg m −3 at ground level) during the 2-month period. In contrast, median daily aircraft contribution to grid concentrations varies between 0.015 and 0.091 μg m −3 (0.09 and 0.40 μg m −3 at ground level), while the maximum varies between 0.75 and 2.55 μg m −3 (0.75 and 2.0 μg m −3 at ground level). Future researchers may consider using a plume-in-grid process, such as the one used here, to understand the impacts of aircraft emissions at other airports, for proposed future airports, for airport expansion projects under various future scenarios, and for other national-scale studies specifically when the maximum impacts at fine scales are of interest.

Description: Corrigendum to "Atmospheric column-averaged mole fractions of carbon dioxide at 53 aircraft measurement sites" published in Atmos. Chem. Phys. 13, 5265–5275, 2013 Atmospheric Chemistry and Physics, 13, 9213-9216, 2013 Author(s): Y. Miyamoto, M. Inoue, I. Morino, O. Uchino, T. Yokota, T. Machida, Y. Sawa, H. Matsueda, C. Sweeney, P. P. Tans, A. E. Andrews, S. C. Biraud, and P. K. Patra No abstract available.

Description: Timing of sea ice retreat can alter phytoplankton community structure in the western Arctic Ocean Biogeosciences Discussions, 10, 15153-15180, 2013 Author(s): given_name prefix surname suffix, A. Fujiwara, T. Hirawake, K. Suzuki, I. Imai, and S.-I. Saitoh This study assesses the response of phytoplankton assemblages to recent climate change, especially with regard to the shrinking of sea ice in the northern Chukchi Sea of the western Arctic Ocean. Distribution patterns of phytoplankton groups in the late summers of 2008–2010 were analyzed based on HPLC pigment signatures and, the following four major algal groups were inferred via multiple regression and cluster analyses: prasinophytes, diatoms, haptophytes and dinoflagellates. A remarkable interannual difference in the distribution pattern of the groups was found in the northern basin area. Haptophytes dominated and dispersed widely in warm surface waters in 2008, whereas prasinophytes dominated in cold water in 2009 and 2010. A difference in the onset date of sea ice retreat was evident among years – the sea ice retreat in 2008 was 1–2 months earlier than in 2009 and 2010. The spatial distribution of early sea ice retreat matched the areas in which a shift in algal community composition was observed. Steel-Dwass's multiple comparison tests were used to assess the physical, chemical and biological parameters of the four clusters. We found a statistically significant difference in temperature between the haptophyte-dominated cluster and the other clusters, suggesting that the change in the phytoplankton communities was related to the earlier sea ice retreat in 2008 and the corollary increase in sea surface temperatures. Longer periods of open water during the summer, which are expected in the future, may affect food webs and biogeochemical cycles in the western Arctic due to shifts in phytoplankton community structure.

Description: Can we determine what controls the spatio-temporal distribution of d-excess and 17 O-excess in precipitation using the LMDZ general circulation model? Climate of the Past, 9, 2173-2193, 2013 Author(s): C. Risi, A. Landais, R. Winkler, and F. Vimeux Combined measurements of the H 2 18 O and HDO isotopic ratios in precipitation, leading to second-order parameter D-excess, have provided additional constraints on past climates compared to the H 2 18 O isotopic ratio alone. More recently, measurements of H 2 17 O have led to another second-order parameter: 17 O-excess. Recent studies suggest that 17 O-excess in polar ice may provide information on evaporative conditions at the moisture source. However, the processes controlling the spatio-temporal distribution of 17 O-excess are still far from being fully understood. We use the isotopic general circulation model (GCM) LMDZ to better understand what controls d-excess and 17 O-excess in precipitation at present-day (PD) and during the last glacial maximum (LGM). The simulation of D-excess and 17 O-excess is evaluated against measurements in meteoric water, water vapor and polar ice cores. A set of sensitivity tests and diagnostics are used to quantify the relative effects of evaporative conditions (sea surface temperature and relative humidity), Rayleigh distillation, mixing between vapors from different origins, precipitation re-evaporation and supersaturation during condensation at low temperature. In LMDZ, simulations suggest that in the tropics convective processes and rain re-evaporation are important controls on precipitation D-excess and 17 O-excess. In higher latitudes, the effect of distillation, mixing between vapors from different origins and supersaturation are the most important controls. For example, the lower d-excess and 17 O-excess at LGM simulated at LGM are mainly due to the supersaturation effect. The effect of supersaturation is however very sensitive to a parameter whose tuning would require more measurements and laboratory experiments. Evaporative conditions had previously been suggested to be key controlling factors of d-excess and 17 O-excess, but LMDZ underestimates their role. More generally, some shortcomings in the simulation of 17 O-excess by LMDZ suggest that general circulation models are not yet the perfect tool to quantify with confidence all processes controlling 17 O-excess.

Description: Similarity estimators for irregular and age uncertain time series Climate of the Past Discussions, 9, 5299-5346, 2013 Author(s): K. Rehfeld and J. Kurths Paleoclimate time series are often irregularly sampled and age uncertain, which is an important technical challenge to overcome for successful reconstruction of past climate variability and dynamics. Visual comparison and interpolation-based linear correlation approaches have been used to infer dependencies from such proxy time series. While the first is subjective, not measurable and not suitable for the comparison of many datasets at a time, the latter introduces interpolation bias, and both face difficulties if the underlying dependencies are nonlinear. In this paper we investigate similarity estimators that could be suitable for the quantitative investigation of dependencies in irregular and age uncertain time series. We compare the Gaussian-kernel based cross correlation (gXCF, Rehfeld et al., 2011) and mutual information (gMI, Rehfeld et al., 2013) against their interpolation-based counterparts and the new event synchronization function (ESF). We test the efficiency of the methods in estimating coupling strength and coupling lag numerically, using ensembles of synthetic stalagmites with short, autocorrelated, linear and nonlinearly coupled proxy time series, and in the application to real stalagmite time series. In the linear test case coupling strength increases are identified consistently for all estimators, while in the nonlinear test case the correlation-based approaches fail. The lag at which the time series are coupled is identified correctly as the maximum of the similarity functions in around 60–55% (in the linear case) to 53–42% (for the nonlinear processes) of the cases when the dating of the synthetic stalagmite is perfectly precise. If the age uncertainty increases beyond 5% of the time series length, however, the true coupling lag is not identified more often than the others for which the similarity function was estimated. Age uncertainty contributes up to half of the uncertainty in the similarity estimation process. Time series irregularity contributes less, particularly for the adapted Gaussian-kernel based estimators and the event synchronization function. The introduced link strength concept summarizes the hypothesis test results and balances the individual strengths of the estimators: while gXCF is particularly suitable for short and irregular time series, gMI and the ESF can identify nonlinear dependencies. ESF could, in particular, be suitable to study extreme event dynamics in paleoclimate records. Programs to analyze paleoclimatic time series for significant dependencies are included in a freely available software toolbox.

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: Simplified PRBMs of spatial compliant multi-beam modules for planar motion Mechanical Sciences, 4, 311-318, 2013 Author(s): G. Hao PRBMs (pseudo-rigid-body models) have been becoming important engineering technologies/methods in the field of compliant mechanisms to simplify the design and analysis through the use of the knowledge body of rigid-body mechanisms coupling with springs. This article addresses the PRBMs of spatial multi-beam modules for planar motion, which are composed of three or more symmetrical wire/slender beams parallel to each other where the planar twisting DOF (degree of freedom) is assumed to be very small for specific applications/loading conditions. Simplified PRBMs are firstly proposed through replacing each beam in spatial multi-beam module with a rigid-body link plus two identical spherical joints at its two ends. The characteristics factor, bending stiffness and twisting stiffness for the spherical joint are determined. Load-displacement equations are then derived for a class of spatial multi-beam modules and general spatial multi-beam modules using the virtual work principle and kinematic relationships. Finally, nonlinear FEA (finite element analysis) is employed with comparisons with the PRBMs. The present PRBMs have shown the ability to predict the primary nonlinear constraint characteristics such as load-stiffening effect, cross-axis coupling in the two primary translational directions and buckling load.

Description: We present measurements of higher order clustering of galaxies in the latest release of the Canada–France–Hawaii Telescope Legacy Survey (CFHTLS)-Wide. We construct a series of volume-limited sample of galaxies containing more than one million galaxies over the redshift range 0.2 〈  z  〈 1 in the four independent fields of the CFHTLS-Wide. Using a counts-in-cells technique we measure the variance ${\bar{\xi }}_2$ and the hierarchical moments $S_{n}= {{\bar{\xi }}_n / {\bar{\xi }}_2^{n-1}}$ (3 ≤  n  ≤ 5) as a function of redshift and angular scale. We find that the measured field-to-field scatter in our estimators is in excellent agreement with analytical predictions. At small scales, corresponding to the highly non-linear regime, we find tentative evidence at the 1 level that the hierarchical moments increase with redshift. At large scales, corresponding to the weakly non-linear regime, our measurements are marginally consistent with perturbation theory predictions for standard cold dark matter cosmology using a simple linear bias. The predictions of perturbation theory tend to slightly overestimate our measurements, which may be a signature of non-linear bias.

Description: We estimate cosmic microwave background (CMB) polarization power spectra and temperature–polarization cross-spectra, from the 9-year data of the Wilkinson Microwave Anisotropy Probe ( WMAP ). Foreground cleaning is implemented using minimum variance linear combinations of the coefficients of needlet decompositions of sky maps for all WMAP channels, to produce maps for CMB temperature anisotropies ( T -mode) and polarization ( E and B -modes), for nine different years of observation. The final power spectra are computed from averages of all possible cross-year power spectra obtained using foreground-cleaned maps for the different years. Our analysis technique yields a measurement of the EE spectrum that is in excellent agreement with theoretical expectations from the current cosmological model. By comparison, the publicly available WMAP EE power spectrum is higher on average (and significantly higher than the predicted EE spectrum from the current best fit) at scales larger than about a degree, an excess that is not confirmed by our analysis. Our TE and TB measurements are in good agreement overall with the WMAP ones and are compatible with the theoretical expectations, although a few data points are off by a few standard deviations, and yield a reduced 2 somewhat above expectation. As predicted for a standard cosmological model with low tensor-to-scalar ratio, the EB and BB power spectra obtained in our analysis are compatible with zero.

Description: We present optical integral field unit observations of two gas pillars surrounding the Galactic young massive star cluster NGC 3603. The high S/N and spectral resolution of these data have allowed us to accurately quantify the Hα, [N ii ] and [S ii ] emission line shapes, and we find a mixture of broad (FWHM ~ 70–100 km s –1 ) and narrow (〈50 km s –1 ) components. The broad components are found close to the edges of both pillars, suggesting that they originate in turbulent mixing layers (TMLs) driven by the effect of the star cluster wind. Both pillars exhibit surprisingly high ionized gas densities of 〉10 000 cm –3 . In one pillar we found that these high densities are only found in the narrow component, implying that they must originate from deeper within the pillar than the broad component. From this, together with our kinematical data, we conclude that the narrow component traces a photoevaporation flow, and that the TML forms at the interface with the hot wind. On the pillar surfaces, we find a consistent offset in radial velocity between the narrow (brighter) components of Hα and [N ii ] of ~5–8 km s –1 , for which we were unable to find a satisfactory explanation. We urge the theoretical community to simulate mechanical and radiative cloud interactions in more detail to address the many unanswered questions raised by this study.

Description: By exploiting the data base of early-type galaxy (ETG) members of the WINGS survey of nearby clusters, we address here the long debated question of the origin and shape of the Fundamental Plane (FP). Our data suggest that different physical mechanisms concur in shaping and ‘tilting’ the FP with respect to the virial plane (VP) expectation. In particular, a ‘hybrid solution’ in which the structure of galaxies and their stellar population are the main contributors to the FP tilt seems to be favoured. We find that the bulk of the tilt should be attributed to structural non-homology, while stellar population effects play an important but less crucial role. In addition, our data indicate that the differential FP tilt between the V and K band is due to a sort of entanglement between structural and stellar population effects, for which the inward steepening of colour profiles ( V  –  K ) tends to increase at increasing the stellar mass of ETGs. The same kind of analysis applied to the ATLAS 3 D and Sloan Digital Sky Survey (SDSS) data in common with WINGS ( WSDSS throughout the paper) confirms our results, the only remarkable difference being the less important role that our data attribute to the stellar mass-to-light-ratio (stellar populations) in determining the FP tilt . The ATLAS 3 D data also suggest that the FP tilt depends as well on the dark matter (DM) fraction and on the rotational contribution to the kinetic energy ( V rot /), thus again pointing towards the above-mentioned ‘hybrid solution’. We show that the global properties of the FP, i.e. its tilt and tightness, can be understood in terms of the underlying correlation among mass, structure and stellar population of ETGs, for which, at increasing the stellar mass, ETGs become (on average) ‘older’ and more centrally concentrated. Finally, we show that a Malmquist-like selection effect may mimic a differential evolution of the mass-to-light ratio for galaxies of different masses. This should be taken into account in the studies investigating the amount of the so-called ‘downsizing’ phenomenon.

Description: I present a new algorithm, Curved-sky grAvitational Lensing for Cosmological Light conE simulatioNS ( calclens ), for efficiently computing weak gravitational lensing shear signals from large N -body light cone simulations over a curved sky. This new algorithm properly accounts for the sky curvature and boundary conditions, is able to produce redshift-dependent shear signals including corrections to the Born approximation by using multiple-plane ray tracing and properly computes the lensed images of source galaxies in the light cone. The key feature of this algorithm is a new, computationally efficient Poisson solver for the sphere that combines spherical harmonic transform and multigrid methods. As a result, large areas of sky (~10 000 square degrees) can be ray traced efficiently at high resolution using only a few hundred cores. Using this new algorithm and curved-sky calculations that only use a slower but more accurate spherical harmonic transform Poisson solver, I study the convergence, shear E-mode, shear B-mode and rotation mode power spectra. Employing full-sky E/B-mode decompositions, I confirm that the numerically computed shear B-mode and rotation mode power spectra are equal at high accuracy (1 per cent) as expected from perturbation theory up to second order. Coupled with realistic galaxy populations placed in large N -body light cone simulations, this new algorithm is ideally suited for the construction of synthetic weak lensing shear catalogues to be used to test for systematic effects in data analysis procedures for upcoming large-area sky surveys. The implementation presented in this work, written in c and employing widely available software libraries to maintain portability, is publicly available at http://code.google.com/p/calclens .

Description: Quasar accretion disc winds observed via broad absorption lines (BALs) in the UV produce strong continuous absorption in X-rays. The X-ray absorber is believed to serve critically as a radiative shield to keep the outflow ionizations low enough for radiative driving. However, previous studies have shown that ‘mini-BAL’ and narrow absorption line (NAL) outflows have dramatically less X-ray absorption than BALs. Here, we examine X-ray and rest-frame UV spectra of eight mini-BAL quasars with outflow speeds in the range 0.1–0.2c to test the hypothesis that these extreme speeds require a strong shield. We find that the X-ray absorption is weak or moderate, with neutral-equivalent column densities N H  〈 few 10 22 cm –2 , consistent with mini-BALs at lower speeds. We use photoionization models to show that the amount of shielding consistent with our data is too weak to control the outflow ionizations and, therefore, it is not important for the acceleration. Shielding in complex geometries also seems unlikely because the alleged shield would need to extinguish the ionizing far-UV flux while avoiding detection in X-rays and the near-UV. We argue that the outflow ionizations are kept moderate, instead, by high gas densities in small clouds. If the mini-BALs form at radial distances of the order of R ~ 2 pc from the central quasar (broadly consistent with theoretical models and with the mini-BAL variabilities observed here and in previous work), and the total column densities in the mini-BAL gas are N H 10 21 cm –2 , then the total radial extent of outflow clouds is only R clouds 3 x 10 13 cm in cases of no/weak shielding or R clouds 3 x 10 14 cm behind the maximum shield allowed by our data. This implies radial filling factors R clouds / R 5 x 10 – 6 or 5 x 10 – 5 for the unshielded or maximally shielded cases, respectively. Compared to the transverse sizes 8 x 10 15 cm (based on measured line depths), the outflows have shapes like thin ‘pancakes’ viewed face-on, or they occupy larger volumes like a spray of many dense clouds with a small volume filling factor. These results favour models with magnetic confinement in magnetic disc winds. To the extent that BALs, mini-BALs and NALs probe the same general outflow phenomenon, our result for dense substructures should apply to all three outflow types.

Description: The estimated stellar masses of galaxies are widely used to characterize how the galaxy population evolves over cosmic time. If stellar masses can be estimated in a robust manner, free from any bias, global diagnostics such as the stellar mass function can be used to constrain the physics of galaxy formation. We explore how galaxy stellar masses, estimated by fitting broad-band spectral energy distributions (SEDs) with stellar population models, can be biased as a result of commonly adopted assumptions for the star formation and chemical enrichment histories, recycled fractions and dust attenuation curves of galaxies. We apply the observational technique of broad-band SED fitting to model galaxy SEDs calculated by the theoretical galaxy formation model GALFORM, isolating the effect of each of these assumptions. We find that, averaged over the entire galaxy population, the common assumption of exponentially declining star formation histories does not, by itself, adversely affect stellar mass estimation. However, we also show that this result does not hold when considering galaxies that have undergone a recent burst of star formation. We show that fixing the metallicity in SED fitting or using sparsely sampled metallicity grids can introduce mass-dependent systematics into stellar mass estimates. We find that the common assumption of a star–dust geometry corresponding to a uniform foreground dust screen can cause the stellar masses of dusty model galaxies to be significantly underestimated. Finally, we show that stellar mass functions recovered by applying SED fitting to model galaxies at high redshift can differ significantly in both shape and normalization from the intrinsic mass functions predicted by a given model. In particular, the effects of dust can reduce the normalization at the high-mass end by up to 0.6 dex in some cases. Given these differences, our methodology of using stellar masses estimated from model galaxy SEDs offers a new, self-consistent way to compare model predictions with observations. We conclude that great care should be taken when comparing theoretical galaxy formation models to observational results based on the estimated stellar masses of high-redshift galaxies.

Description: We analyse the 2D correlation function of the Sloan Digital Sky Survey-III Baryon Oscillation Spectroscopic Survey (BOSS) CMASS sample of massive galaxies of the ninth data release to measure cosmic expansion H and the angular diameter distance D A at a mean redshift of 〈 z 〉 = 0.57. We apply, for the first time, a new correlation function technique called clustering wedges μ ( s ). Using a physically motivated model, the anisotropic baryonic acoustic feature in the galaxy sample is detected at a significance level of 4.7 compared to a featureless model. The baryonic acoustic feature is used to obtain model-independent constraints cz / H / r s = 12.28 ± 0.82 (6.7 percent accuracy) and D A / r s =  9.05 ± 0.27 (3.0 per cent) with a correlation coefficient of –0.5, where r s is the sound horizon scale at the end of the baryonic drag era. We conduct thorough tests on the data and 600 simulated realizations, finding robustness of the results regardless of the details of the analysis method. Combining this with r s constraints from the cosmic microwave background, we obtain H (0.57) = 90.8 ± 6.2 km s –1 Mpc –1 and D A (0.57) = 1386 ± 45 Mpc. We use simulations to forecast results of the final BOSS CMASS data set. We apply the reconstruction technique on the simulations demonstrating that the sharpening of the anisotropic baryonic acoustic feature should improve the detection as well as tighten constraints of H and D A by ~30 per cent on average.

Description: We present an analysis of Chandra archival, pre-explosion data of the positions of three nearby (〈25 Mpc) Type Ia supernovae, SN2011iv, SN2012cu and SN2012fr. No sources corresponding to the progenitors were found in any of the observations. Combining all sources with well-defined backgrounds does not reveal any evidence for X-ray emission from the progenitors either. We calculated upper limits on the bolometric luminosities of the progenitors, under the assumption that they were black bodies with effective temperatures between 30 and 150 eV, corresponding to ‘canonical’ supersoft X-ray sources. The upper limits of SN2012fr straddle the Eddington luminosity of canonical supersoft sources, but fainter canonical supersoft sources cannot be ruled out by this study. We also compare our upper limits with known compact binary supersoft X-ray sources. This study is a continuation of the campaign to directly detect or constrain the X-ray characteristics of pre-explosion observations of nearby Type Ia supernova progenitors; with the results reported in Nielsen, Voss & Nelemans (see reference in Introduction), the number of nearby Type Ia supernovae for which pre-explosion images are available in the Chandra archive is now 13 and counting.

Description: In this work, we report the Photodetector Array Camera and Spectrometer (PACS) 100 μm/160 μm detections of a sample of 42 GALEX -selected and far-infrared (FIR)-detected Lyman break galaxies (LBGs) at z  ~ 1 located in the Cosmic Evolution Survey (COSMOS) field and analyse their ultraviolet (UV) to FIR properties. The detection of these LBGs in the FIR indicates that they have a dust content high enough so that its emission can be directly detected. According to a spectral energy distribution (SED) fitting with stellar population templates to their UV-to-near-IR observed photometry, PACS-detected LBGs tend to be bigger ( R eff  ~ 4.1 kpc), more massive [log ( M * /M ) ~ 10.7], dustier [ E s ( B  –  V ) ~ 0.40], redder in the UV continuum (β ~ –0.60) and UV-brighter [log ( L UV /L ) ~ 10.1] than PACS-undetected LBGs. PACS-detected LBGs at z  ~ 1 are mostly disc-like galaxies and are located over the green valley and red sequence of the colour–magnitude diagram of galaxies at their redshift. By using their UV and IR emission, we find that PACS-detected LBGs tend to be less dusty and have slightly higher total star formation rates (SFRs) than other PACS-detected UV-selected galaxies within the same redshift range. As a consequence of the selection effect due to the depth of the FIR observations employed, all our PACS-detected LBGs have total IR luminosities, L IR , higher than 10 11 L and thus are luminous IR galaxies. However, none of the PACS-detected LBGs are in the ultra-luminous IR galaxy (ULIRG) regime, L IR  ≥ 10 12 L , where the FIR observations are complete. The finding of ULIRGs-LBGs at higher redshifts ( z  ~ 3) suggests an evolution of the FIR emission of LBGs with cosmic time. In an IRX–β diagram, PACS-detected LBGs at z  ~ 1 tend to be located around the relation for local starburst similarly to other UV-selected PACS-detected galaxies at the same redshift. Consequently, the dust-correction factors obtained with their UV continuum slope allow us to determine their total SFR, unlike at higher redshifts. However, the dust attenuation derived from UV to NIR SED fitting overestimates the total SFR for most of our PACS-detected LBGs in an age-dependent way: the overestimation factor is higher in younger galaxies. This is likely due to the typical degeneracy between dust attenuation and age in the SED fitting with synthetic templates and highlights the importance of the FIR measurements in the analysis of star-forming galaxies at intermediate redshifts.

Description: We study 23 previously published Kepler targets to perform a consistent grid-based Bayesian asteroseismic analysis and compare our results to those obtained via the Asteroseismic Modelling Portal. We find differences in the derived stellar parameters of many targets and their uncertainties. While some of these differences can be attributed to systematic effects between stellar evolutionary models, we show that the different methodologies deliver incompatible uncertainties for some parameters. Using non-adiabatic models and our capability to measure surface effects, we also investigate the dependency of these surface effects on the stellar parameters. Our results suggest a dependence of the magnitude of the surface effect on the mixing length parameter which also, but only minimally, affects the determination of stellar parameters. While some stars in our sample show no surface effect at all, the most significant surface effects are found for stars that are close to the Sun's position in the Hertzsprung–Russell diagram.

Description: We present evidence for a strong relationship between galaxy size and environment for the quiescent population in the redshift range 1 〈  z  〈 2. Environments were measured using projected galaxy overdensities on a scale of 400 kpc, as determined from ~96 000 K -band-selected galaxies from the UKIDSS Ultra Deep Survey (UDS). Sizes were determined from ground-based K -band imaging, calibrated using space-based CANDELS HST observations in the centre of the UDS field, with photometric redshifts and stellar masses derived from 11-band photometric fitting. From the resulting size–mass relation, we confirm that quiescent galaxies at a given stellar mass were typically ~50 per cent smaller at z  ~ 1.4 compared to the present day. At a given epoch, however, we find that passive galaxies in denser environments are on average significantly larger at a given stellar mass. The most massive quiescent galaxies ( M *  〉 2 10 11  M ) at z  〉 1 are typically 50 per cent larger in the highest density environments compared to those in the lowest density environments. Using Monte Carlo simulations, we reject the null hypothesis that the size–mass relation is independent of environment at a significance 〉4.8 for the redshift range 1 〈  z  〈 2. In contrast, the evidence for a relationship between size and environment is much weaker for star-forming galaxies.

Description: We present the results from the Very Large Array radio continuum total power and polarized intensity observations of Stephan's Quintet at 1.43 and 4.86 GHz, along with complementary 4.85- and 8.35-GHz Effelsberg observations. Our study shows a large envelope of radio emission encompassing all the member galaxies and hence a large volume of intergalactic matter. Infall of the galaxy NGC 7318B produces a ridge of intergalactic, polarized emission, for which the magnetic field strength has been estimated as 11.0 ± 2.2 μG, with an ordered component of 2.6 ± 0.8 μG. The energy density of the field within the ridge area is of the same order as estimates of the thermal component, implying that the magnetic field has a significant role in the dynamics of the intergalactic matter. We also report that the tidal dwarf galaxy candidate SQ-B possesses a strong and highly anisotropic magnetic field, with a total strength equal to 6.5 ± 1.9 μG and an ordered component reaching 3.5 ± 1.2 μG, which is comparable to that found in normal-sized galaxies.

Description: We present a simple and efficient phenomenological model for the two-dimensional two-point galaxy correlation function that works well over a wide range of scales, from large scales down to scales as small as 25 h –1 Mpc. Our model incorporates non-linear effects and a scale-dependent galaxy bias on small scales, and it allows the redshift-space distortions to be scale and direction dependent. We validate our model using LasDamas mock catalogues and apply it to the Sloan Digital Sky Survey (SDSS) Data Release Seven (DR7) luminous red galaxies (LRGs). Using only the monopole and quadrupole of the correlation function measured from the SDSS DR7 LRGs, we obtain improved measurements H ( z ) r s ( z d )/ c  = 0.0433 ± 0.0042, D A ( z )/ r s ( z d ) = 6.59 ± 0.46 and f ( z ) 8 ( z ) = 0.429 ± 0.089 at z  = 0.35, using the scale range 25 〈  s  〈 120 h –1 Mpc. We expect our results and model to be useful in tightening dark energy and gravity constraints from the full analysis of current and future galaxy clustering data.

Description: We study the formation and evolution of filamentary configurations of dark matter haloes in voids. Our investigation uses the high-resolution cold dark matter simulation CosmoGrid to look for void systems resembling the VGS_31 elongated system of three interacting galaxies that was recently discovered by the Void Galaxy Survey inside a large void in the Sloan Digital Sky Survey galaxy redshift survey. H  i data revealed these galaxies to be embedded in a common elongated envelope, possibly embedded in intravoid filament. In the CosmoGrid simulation we look for systems similar to VGS_31 in mass, size and environment. We find a total of eight such systems. For these systems, we study the distribution of neighbour haloes, the assembly and evolution of the main haloes and the dynamical evolution of the haloes, as well as the evolution of the large-scale structure in which the systems are embedded. The spatial distribution of the haloes follows that of the dark matter environment. We find that VGS_31-like systems have a large variation in formation time, having formed between 10 Gyr ago and the present epoch. However, the environments in which the systems are embedded evolved to resemble each other substantially. Each of the VGS_31-like systems is embedded in an intravoid wall, that no later than z  = 0.5 became the only prominent feature in its environment. While part of the void walls retain a rather featureless character, we find that around half of them are marked by a pronounced and rapidly evolving substructure. Five haloes find themselves in a tenuous filament of a few h –1 Mpc long inside the intravoid wall. Finally, we compare the results to observed data from VGS_31. Our study implies that the VGS_31 galaxies formed in the same (proto)filament, and did not meet just recently. The diversity amongst the simulated halo systems indicates that VGS_31 may not be typical for groups of galaxies in voids.

Description: Using 3D numerical hydrodynamical simulations, we show that jets launched prior to Type Ia supernova (SN Ia) explosion in the core-degenerate (CD) scenario can account for the appearance of two opposite lobes (‘Ears’) along the symmetry axis of the SN remnant (SNR). In the double-degenerate (DD) and CD scenarios the merger of the two degenerate compact objects is very likely to lead to the formation of an accretion disc, that might launch two opposite jets. In the CD scenario, these jets interact with the envelope ejected during the preceding common envelope phase. If explosion occurs shortly after the merger process, the exploding gas and the jets will collide with the ejected nebula, leading to SNR with axisymmetric components including ‘Ears’. We also explore the possibility that the jets are launched by the companion white dwarf prior to its merger with the core. This last process is similar to the one where jets are launched in some pre-planetary nebulae. The SNR ‘Ears’ in this case are formed by a spherical SN Ia explosion inside an elliptical planetary nebula-like object. We compare our numerical results with two SNRs – Kepler and G299.2–2.9.

Description: Low-mass star-forming regions are more complex than the simple spherically symmetric approximation that is often assumed. We apply a more realistic infall/outflow physical model to molecular/continuum observations of three late Class 0 protostellar sources with the aims of (a) proving the applicability of a single physical model for all three sources and (b) deriving physical parameters for the molecular gas component in each of the sources. We have observed several molecular species in multiple rotational transitions. The observed line profiles were modelled in the context of a dynamical model which incorporates infall and bipolar outflows, using a three-dimensional radiative transfer code. This results in constraints on the physical parameters and chemical abundances in each source. Self-consistent fits to each source are obtained. We constrain the characteristics of the molecular gas in the envelopes as well as in the molecular outflows. We find that the molecular gas abundances in the infalling envelope are reduced, presumably due to freeze-out, whilst the abundances in the molecular outflows are enhanced, presumably due to dynamical activity. Despite the fact that the line profiles show significant source-to-source variation, which primarily derives from variations in the outflow viewing angle, the physical parameters of the gas are found to be similar in each core.

Description: We present predictions for hydrogen and helium emission line luminosities from circumstellar matter around Type Ia supernovae (SNe Ia) using time dependent photoionization modelling. Early high-resolution ESO/Very Large Telescope (VLT) optical echelle spectra of the SN Ia 2000cx were taken before and up to ~70 d after maximum to probe the existence of such narrow emission lines from the supernova. We detect no such lines, and from our modelling place an upper limit on the mass-loss rate for the putative wind from the progenitor system, $\skew4\dot{M}\lesssim 1.3\times 10^{-5} \,{\rm M}_{\odot }\,{\rm yr}^{-1}$ , assuming a speed of 10 km s –1 and solar abundances for the wind. If the wind would be helium-enriched and/or faster, the upper limit on $\skew4\dot{M}$ could be significantly higher. In the helium-enriched case, we show that the best line to constrain the mass-loss would be He i 10 830. In addition to confirming the details of interstellar Na i and Ca ii absorption towards SN 2000cx as discussed by Patat et al., we also find evidence for 6613.56 Å diffuse interstellar band absorption in the Milky Way. We also discuss measurements of the X-ray emission from the interaction between the supernova ejecta and the wind and we re-evaluate observations of SN 1992A obtained ~16 d after maximum by Schlegel & Petre. We find an upper limit of $\skew4\dot{M}\sim 1.3\times 10^{-5} \,{\rm M}_{\odot }\,{\rm yr}^{-1}$ which is significantly higher than that estimated by Schlegel & Petre. These results, together with the previous observational work on the normal SNe Ia 1994D and 2001el, disfavour a symbiotic star in the upper mass-loss rate regime (so-called Mira-type systems) from being the likely progenitor scenario for these SNe. Our model calculations are general, and can also be used for the subclass of SNe Ia that do show circumstellar interaction, e.g. the recent PTF 11kx. To constrain hydrogen in late-time spectra, we present ESO/VLT and ESO/New Technology Telescope optical and infrared observations of SNe Ia 1998bu and 2000cx in the nebular phase, 251-388 d after maximum. We see no signs of hydrogen line emission in SNe 1998bu and 2000cx at these epochs, and from the absence of Hα with a width of the order of ~10 3 km s –1 , we argue from modelling that the mass of such hydrogen-rich gas must be 0.03 M for both supernovae. Comparing similar upper limits with recent models of Pan et al., it seems that hydrogen-rich donors with a separation of 5 times the radius of the donor may be ruled out for the five SNe Ia 1998bu, 2000cx, 2001el, 2005am and 2005cf. Larger separation, helium-rich donors, or a double-degenerate origin for these supernovae seems more likely. Our models have also been used to put the limit on hydrogen-rich gas in the recent SN 2011fe, and for this supernova, a double-degenerate origin seems likely.

Description: We build differential-reddening maps for 66 Galactic globular clusters (GCs) with archival Hubble Space Telescope WFC/ACS F 606 W and F 814 W photometry. Because of the different GC sizes (characterized by the half-light radius R h ) and distances to the Sun, the WFC/ACS field of view (200 arcsec x 200 arcsec) coverage ( R obs ) lies in the range 1 R obs / R h 15 for about 85 per cent of the sample, with about 10 per cent covering only the inner ( R obs R h ) parts. We divide the WFC/ACS field of view across each cluster in a regular cell grid and extract the stellar-density Hess diagram from each cell, shifting it in colour and magnitude along the reddening vector until matching the mean diagram. Thus, the maps correspond to the internal dispersion of the reddening around the mean. Depending on the number of available stars (i.e. probable members with adequate photometric errors), the angular resolution of the maps range from 7 arcsec x 7 arcsec to 20 arcsec x 20 arcsec. We detect spatially variable extinction in the 66 GCs studied, with mean values ranging from 〈E(B-V)〉0.018 (NGC 6981) up to 〈E(B-V)〉0.016 (Palomar 2). Differential-reddening correction decreases the observed foreground reddening and the apparent distance modulus but, since they are related to the same value of E ( B  –  V ), the distance to the Sun is conserved. Fits to the mean-ridge lines of the highly extincted and photometrically scattered GC Palomar 2 show that age and metallicity also remain unchanged after the differential-reddening correction, but measurement uncertainties decrease because of the reduced scatter. The lack of systematic variations of 〈E(B-V)〉 with both the foreground reddening and the sampled cluster area indicates that the main source of differential reddening is interstellar.

Description: We present a physically motivated semi-analytic model to understand the clustering of high-redshift Lyman α emitters (LAEs). We show that the model parameters constrained by the observed luminosity functions can be used to predict large-scale bias and angular correlation function of LAEs. These predictions are shown to reproduce the observations remarkably well. We find that average masses of dark matter haloes hosting LAEs brighter than the threshold narrow-band magnitude ~25 are ~10 11 M . These are smaller than that of typical Lyman break galaxies (LBGs) brighter than a similar threshold continuum magnitude by a factor of ~10. This results in a smaller clustering strength of LAEs compared to LBGs. However, using the observed relationship between the UV continuum and Lyman α luminosity of LAEs, we show that both LAEs and LBGs belong to the same parent galaxy population with narrow-band techniques having greater efficiency in picking up galaxies with low UV luminosity. We also show that the lack of evidence for the presence of the one-halo term in the observed LAE angular correlation functions can be attributed to a sub-Poisson distribution of LAEs in dark matter haloes as a result of their low halo occupations.

Description: We present a detailed analysis of spectral line profiles in Type Ia supernova (SN Ia) spectra. We focus on the feature at ~3500–4000 Å, which is commonly thought to be caused by blueshifted absorption of Ca H&K. Unlike some other spectral features in SN Ia spectra, this feature often has two overlapping (blue and red) components. It is accepted that the red component comes from photospheric calcium. However, it has been proposed that the blue component is caused by either high-velocity calcium (from either abundance or density enhancements above the photosphere of the supernova, SN) or Si ii 3858. By looking at multiple data sets and model spectra, focusing on spectra near maximum brightness, we conclude that the blue component of the Ca H&K feature is caused by Si ii 3858 for most SNe Ia, although high-velocity calcium is likely important for some SNe. The strength of the Si ii 3858 feature varies strongly with the light-curve shape of an SN. As a result, the velocity measured from a single-Gaussian fit to the full line profile correlates with light-curve shape. The velocity of the Ca H&K component of the profile does not correlate with light-curve shape, contrary to previous claims. We detail the pitfalls of assuming that the blue component of the Ca H&K feature is caused by calcium, with implications for our understanding of SN Ia progenitors, explosions and cosmology.

Description: Subject of this paper is the statistical properties of ellipticity alignments between galaxies evoked by their coupled angular momenta. Starting from physical angular momentum models, we bridge the gap towards ellipticity correlations, ellipticity spectra and derived quantities such as aperture moments, comparing the intrinsic signals with those generated by gravitational lensing, with the projected galaxy sample of Euclid in mind. We investigate the dependence of intrinsic ellipticity correlations on cosmological parameters and show that intrinsic ellipticity correlations give rise to non-Gaussian likelihoods as a result of non-linear functional dependences. Comparing intrinsic ellipticity spectra to weak lensing spectra we quantify the magnitude of their contaminating effect on the estimation of cosmological parameters and find that biases on dark energy parameters are very small in an angular momentum-based model in contrast to the linear alignment model commonly used. Finally, we quantify whether intrinsic ellipticities can be measured in the presence of the much stronger weak lensing induced ellipticity correlations, if prior knowledge on a cosmological model is assumed.

Description: In this paper, I present a new method that has been developed for determining the brightness of a cloudless night sky, on the basis of widely available amateur observations of comets. The tests show the correctness of the method, which makes it possible to determine the level of light pollution, defined as the brightness of the artificial sky glow, through the use of the archival observations of comets. The use of data bases of comet observations in Poland in the period 1994–2009 has led to a positive verification of the known model map of the brightness of the night sky. Also, it has been possible to find changes in the level of light pollution in this period, at the selected observation sites.

Description: We present a complete sample of molecular clumps containing compact and ultracompact H ii (UC H ii ) regions between  = 10° and 60° and | b | 〈 1°, identified by combining the APEX Telescope Large Area Survey of the Galaxy submm and CORNISH radio continuum surveys with visual examination of archival infrared data. Our sample is complete to optically thin, compact and UC H ii regions driven by a zero-age main-sequence star of spectral type B0 or earlier embedded within a 1000 M clump. In total we identify 213 compact and UC H ii regions, associated with 170 clumps. Unambiguous kinematic distances are derived for these clumps and used to estimate their masses and physical sizes, as well as the Lyman continuum fluxes and sizes of their embedded H ii regions. We find a clear lower envelope for the surface density of molecular clumps hosting massive star formation of 0.05 g cm –2 , which is consistent with a similar sample of clumps associated with 6.7 GHz masers. The mass of the most massive embedded stars is closely correlated with the mass of their natal clump. Young B stars appear to be significantly more luminous in the ultraviolet than predicted by current stellar atmosphere models. The properties of clumps associated with compact and UC H ii regions are very similar to those associated with 6.7 GHz methanol masers and we speculate that there is little evolution in the structure of the molecular clumps between these two phases. Finally, we identify a significant peak in the surface density of compact and UC H ii -regions associated with the W49A star-forming complex, noting that this complex is truly one of the most massive and intense regions of star formation in the Galaxy.

Description: We present full Stokes radio polarization observations of the quasar PKS B2126–158 ( z  = 3.268) from 1 to 10 GHz using the Australia Telescope Compact Array. The source has large fractional circular polarization (CP), m c | V |/ I , detected at high significance across the entire band (from 15 to 90 per 128 MHz subband). This allows us to construct the most robust CP spectrum of an active galactic nucleus (AGN) jet to date. We find m c    +0.60 ± 0.03 from 1.5 to 6.5 GHz, with a peak of m c ~ 1 per cent before the spectrum turns over somewhere between 6.5 and 8 GHz, above which m c    –3.0 ± 0.4 . The fractional linear polarization (LP; p ) varies from 0.2 to ~1 per cent across our frequency range and is strongly anticorrelated with the fractional CP, with a best-fitting power law giving m c    p –0.24 ± 0.03 . This is the first clear relation between the observed LP and CP of an AGN jet, revealing the action of Faraday conversion of LP to CP within the jet. More detailed modelling in conjunction with high spatial resolution observations are required to determine the true driving force behind the conversion (i.e. magnetic twist or internal Faraday rotation). In particular determining whether the observed Faraday rotation is internal or entirely external to the jet is key to this goal. The simplest interpretation of our observations favours some internal Faraday rotation, implying that Faraday rotation-driven conversion of LP to CP is the dominant CP generation mechanism. In this case, a small amount of vector-ordered magnetic field along the jet axis is required, along with internal Faraday rotation from the low-energy end of the relativistic electron energy spectrum in an electron–proton-dominated jet.

Description: We present new H i observations of Hoag's Object (HO) obtained with the Westerbork Synthesis Radio Telescope. The data show that the luminous optical ring around the elliptical body has a bright H i counterpart that shares the kinematical properties of the optical ring. The entire H i structure is twice as large as the optical ring and shows a mild warp in its outer regions relative to the inner ring. We detect two additional H i sources close in redshift to that of HO, and report on a newly identified SDSS optical companion galaxy. The H i sources are ~0.3 and ~1 Mpc away in projected distance, and the companion galaxy is also ~1 Mpc away. Our main conclusion is that the H i detected in HO shows no indication that this galaxy has experienced a recent (less than ~1 Gyr ago) accretion event. At least one of the two additional H i detected objects does not have an optical counterpart. One possibility is that this object is an H i filament left over from an interaction shaping HO, in which case this interaction must also have occurred at least 1–2 Gyr ago.

Description: We report the first comparative study of strong Mg ii absorbers ( W r  ≥ 1.0 Å) seen towards radio-loud quasars of core-dominated (CDQ) and lobe-dominated (LDQ) types and normal quasars (QSOs). The CDQ and LDQ samples were derived from the Sloan Digital Sky Survey Data Release 7 after excluding known ‘broad-absorption-line’ quasars and blazars. The Mg ii associated absorption systems having a velocity offset v  〈 5000 km s –1 from the systemic velocity of the background quasar were also excluded. Existing spectroscopic data for redshift-matched sightlines of 3975 CDQs and 1583 LDQs, covering an emission redshift range 0.39–4.87, were analysed and 864 strong Mg ii absorbers were found, covering the redshift range 0.45–2.17. The conclusions reached using this well-defined large data set of strong Mg ii absorbers are (i) the number density, d N /d z , towards CDQs shows a small, marginally significant excess (~9 per cent at 1.5 significance) over the estimate available for QSOs; (ii) in the redshift space, this difference is reflected in terms of a 1.6 excess of d N /d z over the QSOs, within the narrow redshift interval 1.2–1.8; (iii) the d N /dβ distribution (with β =  v / c ) for CDQs shows a significant excess (at 3.75 level) over the distribution found for a redshift- and luminosity-matched sample of QSOs, at β in the range 0.05–0.1. This leads us to infer that a significant fraction of strong Mg ii absorption systems seen in this offset velocity range are probably associated with the CDQs and might be accelerated into the line of sight by their powerful jets and/or due to the accretion-disc outflows close to our direction. Support to this scenario comes from a consistency check in which we consider only the spectral range corresponding to β 〉 0.2. The computed redshift distribution for strong Mg ii absorbers towards CDQs now shows excellent agreement with that known for QSOs, as indeed is expected for purely intervening absorption systems. Thus, it appears that for CDQs (and blazars) the associated strong Mg ii absorbers can be seen at much larger velocities relative to the nucleus than the commonly adopted upper limit of 5000 km s –1 .

Description: Dynamics of a complete sample of small perihelion distance near-parabolic comets discovered in the years 2006–2010 are studied (i.e. of 22 comets of q osc  〈 3.1 au). First, osculating orbits are obtained after a very careful positional data inspection and processing, including where appropriate, the method of data partitioning for determination of pre- and post-perihelion orbit for tracking then its dynamical evolution. The non-gravitational acceleration in the motion is detected for 50 per cent of investigated comets, in a few cases for the first time. Different sets of non-gravitational parameters are determined from pre- and post-perihelion data for some of them. The influence of the positional data structure on the possibility of the detection of non-gravitational effects and the overall precision of orbit determination is widely discussed. Secondly, both original and future orbits were derived by means of numerical integration of swarms of virtual comets obtained using a Monte Carlo cloning method. This method allows us to follow the uncertainties of orbital elements at each step of dynamical evolution. The complete statistics of original and future orbits that includes significantly different uncertainties of 1/ a -values is presented, also in the light of our results obtained earlier. Basing on 108 comets examined by us so far, we conclude that only one of them, C/2007 W1 Boattini, seems to be a serious candidate for an interstellar comet. We also found that 53 per cent of 108 near-parabolic comets escaping in the future from the Solar system, and the number of comets leaving the Solar system as so called Oort spike comets (i.e. comets suffering very small planetary perturbations) is 14 per cent. A new method for cometary orbit quality assessment is also proposed by means of modifying the original method, introduced by Marsden, Sekanina & Everhart. This new method leads to a better diversification of orbit quality classes for contemporary comets.

Description: We present the chemical abundance measurements of the first large, medium-resolution, uniformly selected damped Lyman α system (DLA) survey. The sample contains 99 DLAs towards 89 quasars selected from the SDSS DR5 DLA sample in a uniform way. We analyse the metallicities and kinematic diagnostics, including the velocity width of 90 per cent of the optical depth, v 90 , and the equivalent widths of the Si ii 1526 ( W 1526 ), C iv 1548 and Mg ii 2796 transitions. To avoid strong line-saturation effects on the metallicities measured in medium-resolution spectra (FWHM ~ 71 km s –1 ), we derived metallicities from metal transitions which absorbed at most 35 per cent of the quasar continuum flux. We find the evolution in cosmic mean metallicity of the sample, 〈 Z 〉 = (–0.04 ± 0.13) z  – (1.06 ± 0.36), consistent with no evolution over the redshift range z ~ [2.2, 4.4], but note that the majority of our sample falls at z ~ [2.2, 3.5]. The apparent lack of metallicity evolution with redshift is also seen in a lack of evolution in the median v 90 and W 1526 values. While this result may seem to conflict with other large surveys that have detected significant metallicity evolution, such as Rafelski et al. who found 〈 Z 〉 = (–0.22 ± 0.03) z  – (0.65 ± 0.09) over z ~ [0, 5], several tests show that these surveys are not inconsistent with our new result. However, over the smaller redshift range covered by our uniformly selected sample, the true evolution of the cosmic mean metallicity in DLAs may be somewhat flatter than the Rafelski et al. estimate.

Description: Expulsion of neutron-rich matter following the merger of neutron star binaries is crucial to the radioactively powered electromagnetic counterparts of these events and to their relevance as sources of r -process nucleosynthesis. Here we explore the long-term (viscous) evolution of remnant black hole accretion discs formed in such mergers by means of two-dimensional, time-dependent hydrodynamical simulations. The evolution of the electron fraction due to charged-current weak interactions is included, and neutrino self-irradiation is modelled as a lightbulb that accounts for the disc geometry and moderate optical depth effects. Over several viscous times (~1 s), a fraction of ~10 per cent of the initial disc mass is ejected as a moderately neutron-rich wind ( Y e  ~ 0.2) powered by viscous heating and nuclear recombination, with neutrino self-irradiation playing a sub-dominant role. Although the properties of the outflow vary in time and direction, their mean values in the heavy-element production region are relatively robust to variations in the initial conditions of the disc and the magnitude of its viscosity. The outflow is sufficiently neutron-rich that most of the ejecta forms heavy r -process elements with mass number A   130, thus representing a new astrophysical source of r -process nucleosynthesis, distinct from that produced in the dynamical ejecta. Due to its moderately high entropy, disc outflows contain a small residual fraction ~1 per cent of helium, which could produce a unique spectroscopic signature.

Description: We have performed multiband UBVR C I C JHK S photometry of two young clusters located at large Galactocentric distances in the direction of the Perseus spiral arm. The obtained distances and colour excesses amount to 3.9 ± 0.11 kpc, E ( B  –  V ) = 0.62 ± 0.05 for Berkeley 94, and 4.3 ± 0.15 kpc, E ( B  –  V ) = 0.58 ± 0.06 for Berkeley 96. The respective ages, as measured from the comparison of the upper colour–magnitude diagrams to model isochrones, amount to log 10 Age(yr) = 7.5 ± 0.07 and 7.0 ± 0.07, respectively. A sequence of optical pre-main-sequence (PMS) members is proposed in both clusters. In addition, samples of objects showing ( H  –  K S ) excess are found. Part of these are suggested to be PMS cluster members of lower mass than the optical candidates. The spatial distribution of these sources, the comparison to Galactic models and to the expected number of contaminating distant red galaxies, and the spectral energy distribution in particular cases support this suggestion. The spatial distributions shown by members in different mass ranges can be interpreted in terms of the results from numerical simulations. According to these, different initial conditions and evolutionary dynamical paths are suggested for the clusters. Berkeley 94 would have formed under supervirial conditions, and followed the so-called warm collapse model in its evolution, whereas Berkeley 96 would have formed with a subvirial structure, and would have evolved following a cold collapse path. Both processes would be able to reproduce the suggested degree of mass segregation and their spatial distribution by mass range. Finally, the mass distributions of the clusters, from the most massive stars down to PMS stars around 1.3 M , are calculated. An acceptable general agreement with the Salpeter initial mass function slope is found.