Three-Dimensional Simulations of Bi-Directed Magnetohydrodynamic Jets Interacting with Cluster Environments

We report on a series of three-dimensional magnetohydrodynamic simulations of active galactic nucleus (AGN) jet propagation in realistic models of magnetized galaxy clusters. We are primarily interested in the details of energy transfer between jets and the intracluster medium (ICM) to help clarify what role such flows could have in the reheating of cluster cores. Our simulated jets feature a range of intermittency behaviors, including intermittent jets that periodically switch on and off and one model jet that shuts down completely, naturally creating a relic plume. The ICM into which these jets propagate incorporates tangled magnetic field geometries and density substructure designed to mimic some likely features of real galaxy clusters. We find that our jets are characteristically at least 60% efficient at transferring thermal energy to the ICM. Irreversible heat energy is not uniformly distributed, however, instead residing preferentially in regions very near the jet/cocoon boundaries. While intermittency affects the details of how, when, and where this energy is deposited, all of our models generically fail to heat the cluster cores uniformly. Both the detailed density structure and nominally weak magnetic fields in the ICM play interesting roles in perturbing the flows, particularly when the jets are non-steady. Still, this perturbation is never sufficient to isotropize the jet energy deposition, suggesting that some other ingredient is required for AGN jets to successfully reheat cluster cores.Comment: 19 pages, 18 figures, Accepted for publication in the Astrophysical Journa

Thin ice clouds in the Arctic: cloud optical depth and particle size retrieved from ground-based thermal infrared radiometry

Multiband downwelling thermal measurements of zenith sky radiance, along with cloud boundary heights, were used in a retrieval algorithm to estimate cloud optical depth and effective particle diameter of thin ice clouds in the Canadian High Arctic. Ground-based thermal infrared (IR) radiances for 150 semitransparent ice clouds cases were acquired at the Polar Environment Atmospheric Research Laboratory (PEARL) in Eureka, Nunavut, Canada (80° N, 86° W). We analyzed and quantified the sensitivity of downwelling thermal radiance to several cloud parameters including optical depth, effective particle diameter and shape, water vapor content, cloud geometric thickness and cloud base altitude. A lookup table retrieval method was used to successfully extract, through an optimal estimation method, cloud optical depth up to a maximum value of 2.6 and to separate thin ice clouds into two classes: (1) TIC1 clouds characterized by small crystals (effective particle diameter  ≤  30 µm), and (2) TIC2 clouds characterized by large ice crystals (effective particle diameter  >  30 µm). The retrieval technique was validated using data from the Arctic High Spectral Resolution Lidar (AHSRL) and Millimeter Wave Cloud Radar (MMCR). Inversions were performed over three polar winters and results showed a significant correlation (R2 =  0.95) for cloud optical depth retrievals and an overall accuracy of 83 % for the classification of TIC1 and TIC2 clouds. A partial validation relative to an algorithm based on high spectral resolution downwelling IR radiance measurements between 8 and 21 µm was also performed. It confirms the robustness of the optical depth retrieval and the fact that the broadband thermal radiometer retrieval was sensitive to small particle (TIC1) sizes

Inflating a chain of x-ray deficient bubbles by a single jet activity episode

We show that a continuous jet with time-independent launching properties can inflate a chain of close and overlapping X-ray deficient bubbles. Using the numerical code PLUTO we run 2.5D (i.e. spherical coordinate system with cylindrical symmetry) hydrodynamic simulations and study the interaction of the jets with the intra-cluster medium (ICM). A key process is vortex fragmentation due to several mechanisms, including vortex-shedding and Kelvin-Helmholtz (KH) instabilities. Our results can account for the structure of two opposite chains of close bubbles as observed in the galaxy cluster Hydra A. Our results imply that the presence of multiple pairs of bubbles does not necessarily imply several jet-launching episodes. This finding might have implications to feedback mechanisms operating by jets.Comment: Accepted by ApJ Letter

Aerosol Remote Sensing in Polar Regions

Multi-year sets of ground-based sun-photometer measurements conducted at 12 Arctic sites and 9 Antarctic sites were examined to determine daily mean values of aerosol optical thickness tau(lambda) at visible and near-infrared wavelengths, from which best-fit values of ngstrm's exponent alpha were calculated. Analyzing these data, the monthly mean values of tau(0.50 micrometers) and alpha and the relative frequency histograms of the daily mean values of both parameters were determined for winter-spring and summer-autumn in the Arctic and for austral summer in Antarctica. The Arctic and Antarctic covariance plots of the seasonal median values of alpha versus tau(0.50 micrometers) showed: (i) a considerable increase in tau(0.50 micrometers) for the Arctic aerosol from summer to winter-spring, without marked changes in alpha; and (ii) a marked increase in tau(0.50 micrometer) passing from the Antarctic Plateau to coastal sites, whereas alpha decreased considerably due to the larger fraction of sea-salt aerosol. Good agreement was found when comparing ground-based sun-photometer measurements of tau(lambda) and alpha at Arctic and Antarctic coastal sites with Microtops measurements conducted during numerous AERONET/MAN cruises from 2006 to 2013 in three Arctic Ocean sectors and in coastal and off-shore regions of the Southern Atlantic, Pacific, and Indian Oceans, and the Antarctic Peninsula. Lidar measurements were also examined to characterize vertical profiles of the aerosol backscattering coefficient measured throughout the year at Ny-lesund. Satellite-based MODIS, MISR, and AATSR retrievals of tau(lambda) over large parts of the oceanic polar regions during spring and summer were in close agreement with ship-borne and coastal ground-based sun-photometer measurements. An overview of the chemical composition of mode particles is also presented, based on in-situ measurements at Arctic and Antarctic sites. Fourteen log-normal aerosol number size-distributions were defined to represent the average features of nuclei, accumulation and coarse mode particles for Arctic haze, summer background aerosol, Asian dust and boreal forest fire smoke, and for various background austral summer aerosol types at coastal and high-altitude Antarctic sites. The main columnar aerosol optical characteristics were determined for all 14 particle modes, based on in-situ measurements of the scattering and absorption coefficients. Diurnally averaged direct aerosol-induced radiative forcing and efficiency were calculated for a set of multimodal aerosol extinction models, using various Bidirectional Reflectance Distribution Function models over vegetation-covered, oceanic and snow-covered surfaces. These gave a reliable measure of the pronounced effects of aerosols on the radiation balance of the surface-atmosphere system over polar regions

Cosmological MHD simulations of cluster formation with anisotropic thermal conduction

(abridged) The ICM has been suggested to be buoyantly unstable in the presence of magnetic field and anisotropic thermal conduction. We perform first cosmological simulations of galaxy cluster formation that simultaneously include magnetic fields, radiative cooling and anisotropic thermal conduction. In isolated and idealized cluster models, the magnetothermal instability (MTI) tends to reorient the magnetic fields radially. Using cosmological simulations of the Santa Barbara cluster we detect radial bias in the velocity and magnetic fields. Such radial bias is consistent with either the inhomogeneous radial gas flows due to substructures or residual MTI-driven field rearangements that are expected even in the presence of turbulence. Although disentangling the two scenarios is challenging, we do not detect excess bias in the runs that include anisotropic thermal conduction. The anisotropy effect is potentially detectable via radio polarization measurements with LOFAR and SKA and future X-ray spectroscopic studies with the IXO. We demonstrate that radiative cooling boosts the amplification of the magnetic field by about two orders of magnitude beyond what is expected in the non-radiative cases. At z=0 the field is amplified by a factor of about 10^6 compared to the uniform magnetic field evolved due to the universal expansion alone. Interestingly, the runs that include both radiative cooling and anisotropic thermal conduction exhibit stronger magnetic field amplification than purely radiative runs at the off-center locations. In these runs, shallow temperature gradients away from the cluster center make the ICM neutrally buoyant. The ICM is more easily mixed in these regions and the winding up of the frozen-in magnetic field is more efficient resulting in stronger magnetic field amplification.Comment: submitted to ApJ, higher resolution figures available at: http://www.astro.lsa.umich.edu/~mateuszr

Multi year aerosol characterization in the tropical Andes and in adjacent Amazonia using AERONET measurements

This work focuses on the analysis of columnar aerosol properties in the complex geophysical tropical region of South America within 10-20 South and 50-70 West. The region is quite varied and encompasses a significant part of Amazonia (lowlands) as well as high mountains in the Andes (highlands,~4000 m a.s.l.). Several AERONET stations were included to study the aerosol optical characteristics of the lowlands (Rio Branco, Ji Parana and Cuiaba in Brazil and Santa Cruz in Bolivia) and the highlands (La Paz, Bolivia) during the 2000-2014 period. Biomass-burning is by far the most important source of aerosol in the lowlands, particularly during the dry season (August-October). Multi-annual variability was investigated and showed very strong burning activity in 2005, 2006, 2007 and 2010. This resulted in smoke characterized by correspondingly strong, above-average AODs (aerosol optical depths) and homogeneous single scattering albedo (SSA) across all the stations (~0.93). For other years, however, SSA differences arise between the northern stations (Rio Branco and Ji Parana) with SSAs of ~0.95 and the southern stations (Cuiaba and Santa Cruz) with lower SSAs of ~0.85.Such differences are explained by the different types of vegetation burned in the two different regions. In the highlands, however, the transport of biomass burning smoke is found to be sporadic in nature. This sporadicity results in highly variable indicators of aerosol load and type (Angstrom exponent and fine mode fraction) with moderately significant increases in both. Regional dust and local pollution are the background aerosol in this highland region, whose elevation places it close to the free troposphere. Transported smoke particles were generally found to be more optical absorbing than in the lowlands: the hypothesis to explain this is the significantly higher amount of water vapor in Amazonia relative to the high mountain areas. The air-mass transport to La Paz was investigated using the HYSPLIT air-concentration five-days back trajectories. Two different patterns were clearly differentiated: westerly winds from the Pacific that clean the atmosphere and easterly winds favoring the transport of particles from Amazonia.Marie Skłodowska-Curie Individual Fellowships (IF) ACE_GFAT (grant agreement No 659398).European Union's Horizon 2020 Research and Innovation Programme under grant agreement No. 654109, ACTRIS-2

Secular evolution of galactic discs: constraints on phase-space density

It was argued in the past that bulges of galaxies cannot be formed through collisionless secular evolution because that would violate constraints on the phase-space density: the phase-space density in bulges is several times larger than in the inner parts of discs. We show that these arguments against secular evolution are not correct. Observations give estimates of the coarsely grained phase-space densities of galaxies, f'=rho_s/(sigma_R sigma_phi sigma_z), where rho_s is stellar density and sigma_R, sigma_phi, sigma_z are the radial, tangential, and vertical rms velocities of stars. Using high-resolution N-body simulations, we study the evolution of f' in discs of Galaxy-size models. During the secular evolution, the discs, which are embedded in live CDM haloes, form a bar and then a thick, dynamically hot, central mass concentration. In the course of evolution f' declines at all radii, not just in the central region. However, the decline is different in different parts of the disc. In the inner disc, f'(R) develops a valley with a minimum around the end of the central mass concentration. The final result is that the values of f' in the central regions are significantly larger than those in the inner disc. The minimum, which gets deeper with time, seems to be due to a large phase mixing produced by the outer bar. We find that the shape and the amplitude of f'(R) for different simulations agree qualitatively with the observed f'(R) in our Galaxy. Curiously enough, the fact that the coarsely grained phase-space density of the bulge is significantly larger than the one of the inner disc turns out to be an argument in favor of secular formation of bulges, not against it.Comment: 9 pages, 5 figures included. Accepted for publication in MNRAS. Minor changes after referee's report. Two figures added (possition-velocity diagrams) to show (i) the agreement in the mass distribution of one of our models with that of the Galaxy, and (ii) the (minor) influence of gas on this distributio

Retrieving the global distribution of the threshold of wind erosion from satellite data and implementing it into the Geophysical Fluid Dynamics Laboratory land–atmosphere model (GFDL AM4.0/LM4.0)

Dust emission is initiated when surface wind velocities exceed the threshold of wind erosion. Many dust models used constant threshold values globally. Here we use satellite products to characterize the frequency of dust events and land surface properties. By matching this frequency derived from Moderate Resolution Imaging Spectroradiometer (MODIS) Deep Blue aerosol products with surface winds, we are able to retrieve a climatological monthly global distribution of the wind erosion threshold (Vthreshold) over dry and sparsely vegetated surfaces. This monthly two-dimensional threshold velocity is then implemented into the Geophysical Fluid Dynamics Laboratory coupled land–atmosphere model (AM4.0/LM4.0). It is found that the climatology of dust optical depth (DOD) and total aerosol optical depth, surface PM10 dust concentrations, and the seasonal cycle of DOD are better captured over the “dust belt” (i.e., northern Africa and the Middle East) by simulations with the new wind erosion threshold than those using the default globally constant threshold. The most significant improvement is the frequency distribution of dust events, which is generally ignored in model evaluation. By using monthly rather than annual mean Vthreshold, all comparisons with observations are further improved. The monthly global threshold of wind erosion can be retrieved under different spatial resolutions to match the resolution of dust models and thus can help improve the simulations of dust climatology and seasonal cycles as well as dust forecasting

Trans-Pacific Transport of Saharan Dust to Western North America: A Case Study

The first documented case of long range transport of Saharan dust over a pathway spanning Asia and the Pacific to Western North America is described. Crustal material generated by North African dust storms during the period 28 February - 3 March 2005 reached western Canada on 13-14 March 2005 and was observed by lidar and sunphotometer in the Vancouver region and by high altitude aerosol instrumentation at Whistler Peak. Global chemical models (GEOS-CHEM and NRL NAAPS) confirm the transport pathway and suggest source attribution was simplified in this case by the distinct, and somewhat unusual, lack of dust activity over Eurasia (Gobi and Takla Makan deserts) at this time. Over western North America, the dust layer, although subsiding close to the boundary layer, did not appear to contribute to boundary layer particulate matter concentrations. Furthermore, sunphotometer observations (and associated inversion products) suggest that the dust layer had only subtle optical impact (Aerosol Optical Thickness (Tau(sub a500)) and Angstrom exponent (Alpha(sub 440-870) were 0.1 and 1.2 respectively) and was dominated by fine particulate matter (modes in aerodynamic diameter at 0.3 and 2.5microns). High Altitude observations at Whistler BC, confirm the crustal origin of the layer (rich in Ca(++) ions) and the bi-modal size distribution. Although a weak event compared to the Asian Trans-Pacific dust events of 1998 and 2001, this novel case highlights the possibility that Saharan sources may contribute episodically to the aerosol burden in western North America