Simulation techniques for cosmological simulations

Modern cosmological observations allow us to study in great detail the evolution and history of the large scale structure hierarchy. The fundamental problem of accurate constraints on the cosmological parameters, within a given cosmological model, requires precise modelling of the observed structure. In this paper we briefly review the current most effective techniques of large scale structure simulations, emphasising both their advantages and shortcomings. Starting with basics of the direct N-body simulations appropriate to modelling cold dark matter evolution, we then discuss the direct-sum technique GRAPE, particle-mesh (PM) and hybrid methods, combining the PM and the tree algorithms. Simulations of baryonic matter in the Universe often use hydrodynamic codes based on both particle methods that discretise mass, and grid-based methods. We briefly describe Eulerian grid methods, and also some variants of Lagrangian smoothed particle hydrodynamics (SPH) methods.Comment: 42 pages, 16 figures, accepted for publication in Space Science Reviews, special issue "Clusters of galaxies: beyond the thermal view", Editor J.S. Kaastra, Chapter 12; work done by an international team at the International Space Science Institute (ISSI), Bern, organised by J.S. Kaastra, A.M. Bykov, S. Schindler & J.A.M. Bleeke

Galaxy alignments: Observations and impact on cosmology

Galaxy shapes are not randomly oriented, rather they are statistically aligned in a way that can depend on formation environment, history and galaxy type. Studying the alignment of galaxies can therefore deliver important information about the physics of galaxy formation and evolution as well as the growth of structure in the Universe. In this review paper we summarise key measurements of galaxy alignments, divided by galaxy type, scale and environment. We also cover the statistics and formalism necessary to understand the observations in the literature. With the emergence of weak gravitational lensing as a precision probe of cosmology, galaxy alignments have taken on an added importance because they can mimic cosmic shear, the effect of gravitational lensing by large-scale structure on observed galaxy shapes. This makes galaxy alignments, commonly referred to as intrinsic alignments, an important systematic effect in weak lensing studies. We quantify the impact of intrinsic alignments on cosmic shear surveys and finish by reviewing practical mitigation techniques which attempt to remove contamination by intrinsic alignments.Comment: 52 pages excl. references, 16 figures; minor changes to match version published in Space Science Reviews; part of a topical volume on galaxy alignments, with companion papers arXiv:1504.05456 and arXiv:1504.0554

Land cover classification using fuzzy rules and aggregation of contextual information through evidence theory

Land cover classification using multispectral satellite image is a very challenging task with numerous practical applications. We propose a multi-stage classifier that involves fuzzy rule extraction from the training data and then generation of a possibilistic label vector for each pixel using the fuzzy rule base. To exploit the spatial correlation of land cover types we propose four different information aggregation methods which use the possibilistic class label of a pixel and those of its eight spatial neighbors for making the final classification decision. Three of the aggregation methods use Dempster-Shafer theory of evidence while the remaining one is modeled after the fuzzy k-NN rule. The proposed methods are tested with two benchmark seven channel satellite images and the results are found to be quite satisfactory. They are also compared with a Markov random field (MRF) model-based contextual classification method and found to perform consistently better.Comment: 14 pages, 2 figure

Polar mesoscale cyclones in the northeast Atlantic: Comparing climatologies from ERA-40 and satellite imagery

Polar mesoscale cyclones over the subarctic are thought to be an important component of the coupled atmosphere–ocean climate system. However, the relatively small scale of these features presents some concern as to their representation in the meteorological reanalysis datasets that are commonly used to drive ocean models. Here polar mesocyclones are detected in the 40-Year European Centre for Medium-Range Weather Forecasts (ECMWF) Re-Analysis dataset (ERA-40) in mean sea level pressure and 500-hPa geopotential height, using an automated cyclone detection algorithm. The results are compared to polar mesocyclones detected in satellite imagery over the northeast Atlantic, for the period October 1993–September 1995. Similar trends in monthly cyclone numbers and a similar spatial distribution are found. However, there is a bias in the size of cyclones detected in the reanalysis. Up to 80% of cyclones larger than 500 km are detected in MSL pressure, but this hit rate decreases, approximately linearly, to ∼40% for 250-km-scale cyclones and to ∼20% for 100-km-scale cyclones. Consequently a substantial component of the associated air–sea fluxes may be missing from the reanalysis, presenting a serious shortcoming when using such reanalysis data for ocean modeling simulations. Eight maxima in cyclone density are apparent in the mean sea level pressure, clustered around synoptic observing stations in the northeast Atlantic. They are likely spurious, and a result of unidentified shortcomings in the ERA-40 data assimilation procedure

Radial distribution and strong lensing statistics of satellite galaxies and substructure using high resolution LCDM hydrodynamical simulations

We analyse the number density and radial distribution of substructures and satellite galaxies using cosmological simulations that follow the gas dynamics of a baryonic component, including shock heating, radiative cooling and star formation within the hierarchical concordance LCDM model. We find that the dissipation of the baryons greatly enhances the survival of subhaloes, expecially in the galaxy core, resulting in a radial distribution of satellite galaxies that closely follows the overall mass distribution in the inner part of the halo. Hydrodynamical simulations are necessary to resolve the adiabatic contraction and dense cores of galaxies, resulting in a total number of satellites a factor of two larger than found in pure dark matter simulation, in good agreement with the observed spatial distribution of satellite galaxies within galaxies and clusters. Convergence tests show that the cored distribution found by previous authors in pure N-body simulations was due to physical overmerging of dark matter only structures. We proceed to use a ray-shooting technique in order to study the impact of these additional substructures on the number of violations of the cusp caustic magnification relation. We develop a new approach to try to disentangle the effect of substructures from the intrinsic discreteness of N-Body simulations. Even with the increased number of substructures in the centres of galaxies, we are not able to reproduce the observed high numbers of discrepancies observed in the flux ratios of multiply lensed quasars.Comment: 11 pages, 15 figures, comparison with previous works updated, one more plot added, minor changes to match the accepted version by MNRA

Evaluation of Skylab (EREP) data for forest and rangeland surveys

The author has identified the following significant results. Four widely separated sites (near Augusta, Georgia; Lead, South Dakota; Manitou, Colorado; and Redding, California) were selected as typical sites for forest inventory, forest stress, rangeland inventory, and atmospheric and solar measurements, respectively. Results indicated that Skylab S190B color photography is good for classification of Level 1 forest and nonforest land (90 to 95 percent correct) and could be used as a data base for sampling by small and medium scale photography using regression techniques. The accuracy of Level 2 forest and nonforest classes, however, varied from fair to poor. Results of plant community classification tests indicate that both visual and microdensitometric techniques can separate deciduous, conifirous, and grassland classes to the region level in the Ecoclass hierarchical classification system. There was no consistency in classifying tree categories at the series level by visual photointerpretation. The relationship between ground measurements and large scale photo measurements of foliar cover had a correlation coefficient of greater than 0.75. Some of the relationships, however, were site dependent

The PAndAS view of the Andromeda satellite system - I. A Bayesian search for dwarf galaxies using spatial and color-magnitude information

We present a generic algorithm to search for dwarf galaxies in photometric catalogs and apply it to the Pan-Andromeda Archaeological Survey (PAndAS). The algorithm is developed in a Bayesian framework and, contrary to most dwarf-galaxy-search codes, makes use of both the spatial and color-magnitude information of sources in a probabilistic approach. Accounting for the significant contamination from the Milky Way foreground and from the structured stellar halo of the Andromeda galaxy, we recover all known dwarf galaxies in the PAndAS footprint with high significance, even for the least luminous ones. Some Andromeda globular clusters are also recovered and, in one case, discovered. We publish a list of the 143 most significant detections yielded by the algorithm. The combined properties of the 39 most significant isolated detections show hints that at least some of these trace genuine dwarf galaxies, too faint to be individually detected. Follow-up observations by the community are mandatory to establish which are real members of the Andromeda satellite system. The search technique presented here will be used in an upcoming contribution to determine the PAndAS completeness limits for dwarf galaxies. Although here tuned to the search of dwarf galaxies in the PAndAS data, the algorithm can easily be adapted to the search for any localised overdensity whose properties can be modeled reliably in the parameter space of any catalog.Comment: 20 pages, 16 figures, 1 table; accepted for publication in ApJ. High res pdf available at https://www.dropbox.com/s/7zk7pme2wunwkjv/PAndAS_dwarf_galaxies.pd

Spatial snow water equivalent estimation for mountainous areas using wireless-sensor networks and remote-sensing products

We developed an approach to estimate snow water equivalent (SWE) through interpolation of spatially representative point measurements using a k-nearest neighbors (k-NN) algorithm and historical spatial SWE data. It accurately reproduced measured SWE, using different data sources for training and evaluation. In the central-Sierra American River basin, we used a k-NN algorithm to interpolate data from continuous snow-depth measurements in 10 sensor clusters by fusing them with 14 years of daily 500-m resolution SWE-reconstruction maps. Accurate SWE estimation over the melt season shows the potential for providing daily, near real-time distributed snowmelt estimates. Further south, in the Merced-Tuolumne basins, we evaluated the potential of k-NN approach to improve real-time SWE estimates. Lacking dense ground-measurement networks, we simulated k-NN interpolation of sensor data using selected pixels of a bi-weekly Lidar-derived snow water equivalent product. k-NN extrapolations underestimate the Lidar-derived SWE, with a maximum bias of −10 cm at elevations below 3000 m and +15 cm above 3000 m. This bias was reduced by using a Gaussian-process regression model to spatially distribute residuals. Using as few as 10 scenes of Lidar-derived SWE from 2014 as training data in the k-NN to estimate the 2016 spatial SWE, both RMSEs and MAEs were reduced from around 20–25 cm to 10–15 cm comparing to using SWE reconstructions as training data. We found that the spatial accuracy of the historical data is more important for learning the spatial distribution of SWE than the number of historical scenes available. Blending continuous spatially representative ground-based sensors with a historical library of SWE reconstructions over the same basin can provide real-time spatial SWE maps that accurately represents Lidar-measured snow depth; and the estimates can be improved by using historical Lidar scans instead of SWE reconstructions