We introduce the Making Galaxies In a Cosmological Context (MAGICC) programme of smoothed particle hydrodynamics simulations. We describe a parameter study of galaxy formation simulations of an L* ...galaxy that uses early stellar feedback combined with supernova feedback to match the stellar mass-halo mass relationship. While supernova feedback alone can reduce star formation enough to match the stellar mass-halo mass relationship, the galaxy forms too many stars before z = 2 to match the evolution seen using abundance matching. Our early stellar feedback is purely thermal and thus operates like an ultraviolet ionization source as well as providing some additional pressure from the radiation of massive, young stars. The early feedback heats gas to >106 K before cooling it to 104 K. The pressure from this hot gas creates a more extended disc and prevents more star formation prior to z = 1 than supernova feedback alone. The resulting disc galaxy has a flat rotation curve, an exponential surface brightness profile, and matches a wide range of disc scaling relationships. The disc forms from the inside-out with an increasing exponential scale length as the galaxy evolves. Overall, early stellar feedback helps to simulate galaxies that match observational results at low and high redshifts.
A study of metal enrichment of the intergalactic medium (IGM) using a series of smooth particle hydrodynamic (SPH) simulations is presented, employing models for metal cooling and the turbulent ...diffusion of metals and thermal energy. An adiabatic feedback mechanism was adopted where gas cooling was prevented on the time-scale of supernova bubble expansion to generate galactic winds without explicit wind particles. The simulations produced a cosmic star formation history (SFH) that is broadly consistent with observations until z∼ 0.5, and a steady evolution of the universal neutral hydrogen fraction () that compares reasonably well with observations. The evolution of the mass and metallicities in stars and various gas phases was investigated. At z= 0, about 40 per cent of the baryons are in the warm–hot intergalactic medium (WHIM), but most metals (80–90 per cent) are locked in stars. At higher redshifts the proportion of metals in the IGM is higher due to more efficient loss from galaxies. The results also indicate that IGM metals primarily reside in the WHIM throughout cosmic history, which differs from simulations with hydrodynamically decoupled explicit winds. The metallicity of the WHIM lies between 0.01 and 0.1 solar with a slight decrease at lower redshifts. The metallicity evolution of the gas inside galaxies is broadly consistent with observations, but the diffuse IGM is under enriched at z∼ 2.5. Galactic winds most efficiently enrich the IGM for haloes in the intermediate mass range 1010–1011 M⊙. At the low-mass end gas is prevented from accreting on to haloes and has very low metallicities. At the high-mass end, the fraction of halo baryons escaped as winds declines along with the decline of stellar mass fraction of the galaxies. This is likely because of the decrease in star formation activity and decrease in wind escape efficiency. Metals enhance cooling which allows WHIM gas to cool on to galaxies and increases star formation. Metal diffusion allows winds to mix prior to escape, decreasing the IGM metal content in favour of gas within galactic haloes and star-forming gas. Diffusion significantly increases the amount of gas with low metallicities and changes the density–metallicity relation.
We present the McMaster Unbiased Galaxy Simulations (MUGS), the first nine galaxies of an unbiased selection ranging in total mass from 5 × 1011 M⊙ to 2 × 1012 M⊙ simulated using N-body smoothed ...particle hydrodynamics at high resolution. The simulations include a treatment of low-temperature metal cooling, UV background radiation, star formation and physically motivated stellar feedback. Mock images of the simulations show that the simulations lie within the observed range of relations such as that between colour and magnitude and that between brightness and circular velocity (Tully–Fisher). The greatest discrepancy between the simulated galaxies and observed galaxies is the high concentration of material at the centre of the galaxies as represented by the centrally peaked rotation curves and the high bulge-to-total ratios of the simulations determined both kinematically and photometrically. This central concentration represents the excess of low angular momentum material that long has plagued morphological studies of simulated galaxies and suggests that higher resolutions and a more accurate description of feedback will be required to simulate more realistic galaxies. Even with the excess central mass concentrations, the simulations suggest the important role merger history and halo spin play in the formation of discs.
The mountain pine beetle (Dendroctonus ponderosae Hopkins, Coleoptera: Curculionidae, Scolytinae) is a native insect of the pine forests of western North America, and its populations periodically ...erupt into large-scale outbreaks. During outbreaks, the resulting widespread tree mortality reduces forest carbon uptake and increases future emissions from the decay of killed trees. The impacts of insects on forest carbon dynamics, however, are generally ignored in large-scale modelling analyses. The current outbreak in British Columbia, Canada, is an order of magnitude larger in area and severity than all previous recorded outbreaks. Here we estimate that the cumulative impact of the beetle outbreak in the affected region during 2000-2020 will be 270 megatonnes (Mt) carbon (or 36 g carbon m-2 yr-1 on average over 374,000 km2 of forest). This impact converted the forest from a small net carbon sink to a large net carbon source both during and immediately after the outbreak. In the worst year, the impacts resulting from the beetle outbreak in British Columbia were equivalent to 75% of the average annual direct forest fire emissions from all of Canada during 1959-1999. The resulting reduction in net primary production was of similar magnitude to increases observed during the 1980s and 1990s as a result of global change. Climate change has contributed to the unprecedented extent and severity of this outbreak. Insect outbreaks such as this represent an important mechanism by which climate change may undermine the ability of northern forests to take up and store atmospheric carbon, and such impacts should be accounted for in large-scale modelling analyses.
Within a fully cosmological hydrodynamical simulation, we form a galaxy which rotates at 140 km s−1, and it is characterized by two loose spiral arms and a bar, indicative of a Hubble-type SBc/d ...galaxy. We show that our simulated galaxy has no classical bulge, with a pure disc profile at z = 1, well after the major merging activity has ended. A long-lived bar subsequently forms, resulting in the formation of a secularly formed 'pseudo-'bulge, with the final bulge-to-total light ratio of 0.21. We show that the majority of gas which loses angular momentum and falls to the central region of the galaxy during the merging epoch is blown back into the hot halo, with much of it returning later to form stars in the disc. We propose that this mechanism of redistribution of angular momentum via a galactic fountain, when coupled with the results from our previous study which showed why gas outflows are biased to have low angular momentum, can solve the angular momentum/bulgeless disc problem of the cold dark matter paradigm.
Aims. We examine the role of energy feedback in shaping the distribution of metals within cosmological hydrodynamical simulations of L∗ disc galaxies. While negative abundance gradients today provide ...a boundary condition for galaxy evolution models, in support of inside-out disc growth, empirical evidence as to whether abundance gradients steepen or flatten with time remains highly contradictory. Methods. We made use of a suite of L∗ discs, realised with and without “enhanced” feedback. All the simulations were produced using the smoothed particle hydrodynamics code Gasoline, and their in situ gas-phase metallicity gradients traced from redshift z ~ 2 to the present-day. Present-day age-metallicity relations and metallicity distribution functions were derived for each system. Results. The “enhanced” feedback models, which have been shown to be in agreement with a broad range of empirical scaling relations, distribute energy and re-cycled ISM material over large scales and predict the existence of relatively “flat” and temporally invariant abundance gradients. Enhanced feedback schemes reduce significantly the scatter in the local stellar age-metallicity relation and, especially, the O/Fe−Fe/H relation. The local O/Fe distribution functions for our L∗ discs show clear bimodality, with peaks at O/Fe = −0.05 and + 0.05 (for stars with Fe/H > −1), consistent with our earlier work on dwarf discs. Conclusions. Our results with “enhanced” feedback are inconsistent with our earlier generation of simulations realised with “conservative” feedback. We conclude that spatially-resolved metallicity distributions, particularly at high-redshift, offer a unique and under-utilised constraint on the uncertain nature of stellar feedback processes.
Abstract
We explore the circumgalactic medium (CGM) of two simulated star-forming galaxies with luminosities L ≈ 0.1 and 1 L
⋆ generated using the smooth particle hydrodynamic code gasoline. These ...simulations are part of the Making Galaxies In a Cosmological Context (magicc) program in which the stellar feedback is tuned to match the stellar mass-halo mass relationship. For comparison, each galaxy was also simulated using a 'lower feedback' (LF) model which has strength comparable to other implementations in the literature. The 'magicc feedback' (MF) model has a higher incidence of massive stars and an approximately two times higher energy input per supernova. Apart from the low-mass halo using LF, each galaxy exhibits a metal-enriched CGM that extends to approximately the virial radius. A significant fraction of this gas has been heated in supernova explosions in the disc and subsequently ejected into the CGM where it is predicted to give rise to substantial O vi absorption. The simulations do not yet address the question of what happens to the O vi when the galaxies stop forming stars. Our models also predict a reservoir of cool H i clouds that show strong Lyα absorption to several hundred kpc. Comparing these models to recent surveys with the Hubble Space Telescope, we find that only the MF models have sufficient O vi and H i gas in the CGM to reproduce the observed distributions. In separate analyses, these same MF models also show better agreement with other galaxy observables (e.g. rotation curves, surface brightness profiles and H i gas distribution). We infer that the CGM is the dominant reservoir of baryons for galaxy haloes.
Canada's managed boreal forest, 54% of the nation's total boreal forest area, stores 28 Pg carbon (C) in biomass, dead organic matter, and soil pools. The net C balance is dominated by the difference ...of two large continuous fluxes: C uptake (net primary production) and release during decomposition (heterotrophic respiration). Additional releases of C can be high in years, or in areas, that experience large anthropogenic or natural disturbances. From 1990 to 2008, Canada's managed boreal forest has acted as C sink of 28 Tg C year−1, removing CO2from the atmosphere to replace the 17 Tg of C annually harvested and store an additional 11 Tg of C year−1in ecosystem C pools. A large fraction (57%) of the C harvested since 1990 remains stored in wood products and solid waste disposal sites in Canada and abroad, replacing C emitted from the decay or burning of wood harvested prior to 1990 and contributing to net increases in product and landfill C pools. Wood product use has reduced emissions in other sectors by substituting for emission-intensive products (concrete, steel). The C balance of the unmanaged boreal forest is currently unknown. The future C balance of the Canadian boreal forest will affect the global atmospheric C budget and influence the mitigation efforts required to attain atmosphericCO2stabilization targets. The single biggest threat toCstocks is human-caused climate change. LargeCstocks have accumulated in the boreal because decomposition is limited by cold temperatures and often anoxic environments. Increases in temperatures and disturbance rates could result in a large net C source during the remainder of this century and beyond. Uncertainties about the impacts of global change remain high, but we emphasize the asymmetry of risk: sustained large-scale increases in productivity are unlikely to be of sufficient magnitude to offset higher emissions from increased disturbances and heterotrophic respiration. Reducing the uncertainties of the current and future C balance of Canada's 270 Mha of boreal forest requires addressing gaps in monitoring, observation, and quantification of forest C dynamics, with particular attention to 125 Mha of unmanaged boreal forest with extensive areas of deep organic soils, peatlands, and permafrost containing large quantities of C that are vulnerable to global warming.
Abstract
We use the same physical model to simulate four galaxies that match the relation between stellar and total mass, over a mass range that includes the vast majority of disc galaxies. The ...resultant galaxies, part of the Making Galaxies in a Cosmological Context (MaGICC) programme, also match observed relations between luminosity, rotation velocity, size, colour, star formation rate, H i mass, baryonic mass and metallicity. Radiation energy feedback from massive stars and supernova energy balance the complex interplay between cooling gas, regulated star formation, large-scale outflows and recycling of gas in a manner which correctly scales with the mass of the galaxy. Outflows, driven by the expansion of shells and superbubbles of overlapping supernova explosions, also play a key role in simulating galaxies with exponential surface brightness profiles, flat rotation curves and dark matter cores. Our study implies that large-scale outflows are the primary driver of the dependence of disc galaxy properties on mass. We show that the degree of outflows invoked in our model is required to meet the constraints provided by observations of O vi absorption lines in the circumgalactic media of nearby galaxies.
We compare the results of various cosmological gas-dynamical codes used to simulate the formation of a galaxy in the Λ cold dark matter structure formation paradigm. The various runs (13 in total) ...differ in their numerical hydrodynamical treatment smoothed particle hydrodynamics (SPH), moving mesh and adaptive mesh refinement but share the same initial conditions and adopt in each case their latest published model of gas cooling, star formation and feedback. Despite the common halo assembly history, we find large code-to-code variations in the stellar mass, size, morphology and gas content of the galaxy at z= 0, due mainly to the different implementations of star formation and feedback. Compared with observation, most codes tend to produce an overly massive galaxy, smaller and less gas rich than typical spirals, with a massive bulge and a declining rotation curve. A stellar disc is discernible in most simulations, although its prominence varies widely from code to code. There is a well-defined trend between the effects of feedback and the severity of the disagreement with observed spirals. In general, models that are more effective at limiting the baryonic mass of the galaxy come closer to matching observed galaxy scaling laws, but often to the detriment of the disc component. Although numerical convergence is not particularly good for any of the codes, our conclusions hold at two different numerical resolutions. Some differences can also be traced to the different numerical techniques; for example, more gas seems able to cool and become available for star formation in grid-based codes than in SPH. However, this effect is small compared to the variations induced by different feedback prescriptions. We conclude that state-of-the-art simulations cannot yet uniquely predict the properties of the baryonic component of a galaxy, even when the assembly history of its host halo is fully specified. Developing feedback algorithms that can effectively regulate the mass of a galaxy without hindering the formation of high angular momentum stellar discs remains a challenge.