The effects of frequency-dependent quasar variability on the celestial reference frame

We examine the relationship between source position stability and astrophysical properties of radio-loud quasars making up the International Celestial Reference Frame (ICRF2). Understanding this relationship is important for improving quasar selection and analysis strategies, and therefore reference frame stability. We construct flux density time series, known as light curves, for 95 of the most frequently observed ICRF2 quasars at both the 2.3 and 8.4 GHz geodetic very long baseline interferometry (VLBI) observing bands. Because the appearance of new quasar components corresponds to an increase in quasar flux density, these light curves alert us about potential changes in source structure before they appear in VLBI images. We test how source position stability depends on three astrophysical parameters: (1) flux density variability at X band; (2) time lag between flares in S and X bands; (3) spectral index root-mean-square (rms), defined as the variability in the ratio between S and X band flux densities. We find that the time lag between S and X band light curves provides a good indicator of position stability: sources with time lags $<$ 0.06 years are significantly more stable ( $>$ 20 % improvement in weighted rms) than sources with larger time lags. A similar improvement is obtained by observing sources with low $(<$ 0.12) spectral index variability. On the other hand, there is no strong dependence of source position stability on flux density variability in a single frequency band. These findings can be understood by interpreting the time lag between S and X band light curves as a measure of the size of the source structure. Monitoring of source flux density at multiple frequencies therefore appears to provide a useful probe of quasar structure on scales important to geodesy. The observed astrometric position of the brightest quasar component (the core) is known to depend on observing frequency. We show how multi-frequency flux density monitoring may allow the dependence on frequency of the relative core positions along the jet to be elucidated. Knowledge of the position–frequency relation has important implications for current and future geodetic VLBI programs, as well as the alignment between the radio and optical celestial reference frames.     

Appraisal of land use/land cover of mangrove forest ecosystem using support vector machine

Human activities in many parts of the world have greatly changed the natural land cover. This study has been conducted on Pichavaram forest, south east coast of India, famous for its unique mangrove bio-diversity. The main objectives of this study were focused on monitoring land cover changes particularly for the mangrove forest in the Pichavaram area using multi-temporal Landsat images captured in the 1991, 2000, and 2009. The land use/land cover (LULC) estimation was done by a unique hybrid classification approach consisting of unsupervised and support vector machine (SVM)-based supervised classification. Once the vegetation and non-vegetation classes were separated, training site-based classification technology i.e., SVM-based supervised classification technique was used. The agricultural area, forest/plantation, degraded mangrove and mangrove forest layers were separated from the vegetation layer. Mud flat, sand/beach, swamp, sea water/sea, aquaculture pond, and fallow land were separated from non-vegetation layer. Water logged areas were delineated from the area initially considered under swamp and sea water-drowned areas. In this study, the object-based post-classification comparison method was employed for detecting changes. In order to evaluate the performance, an accuracy assessment was carried out using the randomly stratified sampling method, assuring distribution in a rational pattern so that a specific number of observations were assigned to each category on the classified image. The Kappa accuracy of SVM classified image was highest (94.53 %) for the 2000 image and about 94.14 and 89.45 % for the 2009 and 1991 images, respectively. The results indicated that the increased anthropogenic activities in Pichavaram have caused an irreversible loss of forest vegetation. These findings can be used both as a strategic planning tool to address the broad-scale mangrove ecosystem conservation projects and also as a tactical guide to help managers in designing effective restoration measures.     

Overstepping the garnet isograd: a comparison of QuiG barometry and thermodynamic modeling

The consequences of overstepping the garnet isograd reaction have been investigated by comparing the composition of garnet formed at overstepped P–T conditions (the overstep or “OS” model) with the P–T conditions that would be inferred by assuming garnet nucleated in equilibrium with the matrix assemblage at the isograd (the equilibrium or “EQ” model). The garnet nucleus composition formed at overstepped conditions is calculated as the composition that produces the maximum decrease in Gibbs free energy from the equilibrated, garnet-absent, matrix assemblage for the bulk composition under study. Isopleths were then calculated for this garnet nucleus composition assuming equilibrium with the matrix assemblage (the EQ model). Comparison of the actual P–T conditions of nucleation (the OS model) with those inferred from the EQ model reveals considerable discrepancy between the two. In general, the inferred garnet nucleation P–T conditions (the EQ model) are at a lower temperature and higher or lower pressure (depending on the coexisting calcic phase(s)) than the actual (OS model) nucleation conditions. Moreover, the degree of discrepancy increases with the degree of overstepping. Independent estimates of the pressure of nucleation of garnet were made using the Raman shift of quartz inclusions in garnet (quartz-in-garnet or QuiG barometry). To test the validity of this method, an experimental synthesis of garnet containing quartz inclusions was made at 800 °C, 20 kbar, and the measured Raman shift reproduced the synthesis conditions to within 120 bars. Raman band shifts from three natural samples were then used to calculate an isochore along which garnet was presumed to have nucleated. Model calculations were made at several temperatures along this isochore (the OS model), and these P–T conditions were compared to those computed assuming equilibrium nucleation (the EQ model) to estimate the degree of overstepping displayed by these samples. A sample from the garnet isograd in eastern Vermont is consistent with overstepping of around 10 degrees and 0.6 kbar (affinities of around 2 kJ/mole garnet). A sample from the staurolite–kyanite zone in the same terrane requires overstepping of around 50 °C and 2–5 kbar (affinities of around 10–18 kJ/mole garnet). A similar amount of overstepping was inferred for a blueschist sample from Sifnos, Greece. These results indicate that overstepping of garnet nucleation reactions may be common and pronounced in regionally metamorphosed terranes, and that the P–T conditions and paths inferred from garnet zoning studies may be egregiously in error.     

Discrete element modelling of material non‐coaxiality in simple shear flows

We investigate the quasi‐static simple shear flow of a two‐dimensional assembly of cohesionless particles using discrete element method (DEM) simulations. We focus on the unsteady flow regime where the solid would experience significant evolution of stresses, mobilised shear strength and dilation. We construct the DEM model using a discretised‐wall confined granular cell where the apparent boundary is allowed to dilate or contract synchronously with the confined solid. A rather uniform simple shear field is achieved across the whole assembly, which benefits rheological studies in generalising constitutive laws for continuum methods. We examine two aspects of the simple shear behaviour: macroscopic stress and strain rate evolution, particularly the non‐coaxiality between the principal directions of the two; and micromechanics such as evolution of fabric. For an initially anisotropic specimen sheared under constant normal pressure, the direction of principal stress rotates towards that of the principal strain rate, gradually reducing the degree of non‐coaxiality from about 45° to fluctuating around 0°. The rate in approaching coaxiality is slower in samples with larger initial porosity, stress ratio and mean stress. Generally, a faster rate in approaching coaxiality in simple shear is observed in a more dilatant sample, which often shows a larger degree of mobilised fabric anisotropy, suggesting the possible important role of instantaneous internal friction angle. The evolution of principal fabric direction resembles that of the principal stress direction. © 2013 The Authors. International Journal for Numerical and Analytical Methods in Geomechanics published by John Wiley & Sons, Ltd.     

Geochemical features and origin of natural gas in heavy oil area of the Western Slope,Songliao Basin,China

The Western Slope of the Songliao Basin is rich in heavy oil resources (>70 × 10⁸ bbl), around which there are shallow gas reservoirs (∼1.0 × 10¹² m³). The gas is dominated by methane with a dryness over 0.99, and the non-hydrocarbon component being overwelmingly nitrogen. Carbon isotope composition of methane and its homologs is depleted in ¹³C, with δ¹³C₁ values being in the range of −55‰ to −75‰, δ¹³C₂ being in the range of −40‰ to −53‰ and δ¹³C₃ being in the range of −30‰ to −42‰, respectively. These values differ significantly from those solution gases source in the Daqing oilfield. This study concludes that heavy oils along the Western Slope were derived from mature source rocks in the Qijia-Gulong Depression, that were biodegraded. The low reservoir temperature (30–50 °C) and low salinity of formation water with neutral to alkaline pH (NaHCO₃) appeared ideal for microbial activity and thus biodegradation. Natural gas along the Western Slope appears mainly to have originated from biodegradation and the formation of heavy oil. This origin is suggested by the heavy δ¹³C of CO₂ (−18.78‰ to 0.95‰) which suggests that the methane was produced via fermentation as the terminal decomposition stage of the oil.     

Two methods for pore network of porous media

Two approaches of generating pore networks of porous media are presented to capture the pore fabric. The first methodology extracted pore structure from a computer simulated packing of spheres. The modified Delaunay tessellation was used to describe the porous media, and modified Nelder–Mead method in conjunction with three pore‐merging algorithms was used to generate the pore size and coordination number distributions of the randomly packed spheres. The Biconical Abscissa Asymmetric CONcentric bond was used to describe the connection between two adjacent voids. This algorithm was validated by predicting pore structure of a cubic array of spheres of equal radius with known pore sizes, throat sizes and coordination number distributions. The predicted distributions of pore structure agreed well with the measured. Then, the algorithm was used to predict pore structure and permeability of randomly packed spherical particles, and predicted permeability values were compared with published experimental data. The results showed that the predicted permeability values were in good agreement with those measured, confirming the proposed algorithm can capture the main flow paths of packed beds. The second methodology generated an equivalent pore network of porous media, of which the centers of voids were located in a regular lattice with constant pore center distance. However, this network allowed for matching both main geometrical and topological characteristics of the porous media. A comparison of the two approaches suggested that the second approach can also be used as a predictive tool to quantitatively study the microscopic properties of flow through porous media. Copyright © 2012 John Wiley & Sons, Ltd.     

A virtual experiment technique on the elementary behaviour of granular materials with discrete element method

Multi‐scale investigations aided by the discrete element method (DEM) play a vital role for current state‐of‐the‐art research on the elementary behaviour of granular materials. Similar to laboratory tests, there are three important aspects to be considered carefully, which are the proper stress/strain definition and measurement, the application of target loading paths and the designed experiment setup, to be addressed in the present paper. Considering the volume sensitive characteristics of granular materials, in the proposed technique, the deformation of the tested specimen is controlled and measured by deformation gradient tensor involving both the undeformed configuration and the current configuration. Definitions of Biot strain and Cauchy stress are adopted. The expressions of them in terms of contact forces and particle displacements, respectively, are derived. The boundary of the tested specimen consists of rigid massless planar units. It is suggested that the representative element uses a convex polyhedral (polygonal) shape to minimize possible boundary arching effects. General loading paths are described by directly specifying the changes in the stress/strain invariants or directions. Loading can be applied in the strain‐controlled mode by specifying the translations and rotations of the boundary units, or in the stress‐controlled mode by using a servo‐control mechanism, or in the combination of the two methods to realize mixed boundary conditions. Taking the simulation results as the natural consequences originated from a complex system, virtual experiments provide particle‐scale information database to conduct multi‐scale investigations for better understanding in granular material behaviours and possible development of the constitutive theories provided the qualitative similarity between the simulation results from virtual experiments and observations on real material behaviour. Copyright © 2011 John Wiley & Sons, Ltd.     

An experimental study of the effects of the upper boundary condition on the thermal convection in a rotating,differentially heated cylindrical annulus of water

Abstract

The effects of the upper boundary condition on the regime diagram and certain characteristics of the convection within a rotating, differentially heated cylindrical anulus of water were studied in the laboratory. It was found that the regime diagram was not substantially affected by the upper boundary condition. However, the thermal amplitude of the baroclinic waves, as a function of parameter space, and, as expected from previous work, the angular drift velocity of those waves were found to be strongly affected by the upper boundary condition. When the upper surface was free, the amplitude changes were explosive and highly non-linear (as discovered earlier by Kaiser, 1970). When the upper surface was rigid, they were smooth and quite linear. The baroclinic wave patterns drifted round the annulus at rates which were in direct response to the imposed “thermal wind”. However, (as previously observed by Koschmieder, 1972), when the upper surface was rigid they drifted approximately ten times more slowly than when it was free.     

Soil characterization based on land cover heterogeneity over a tropical landscape: an integrated approach using earth observation data-sets

Soil is a vital part of the natural environment and is always responding to changes in environmental factors, along with the influences of anthropogenic factors and land use changes. The long-term change in soil properties will result in change in soil health and fertility, and hence the soil productivity. Hence, the main aim of this paper focuses on the analysis of land use/land cover (LULC) change pattern in spatial and temporal perspective and to present its impact on soil properties in the Merawu catchment over the period of 18?years. Post classification change detection was performed to quantify the decadal changes in historical LULC over the periods of 1991, 2001 and 2009. The pixel to pixel comparison method was used to detect the LULC of the area. The key LULC types were selected for investigation of soil properties. Soil samples were analysed in situ to measure the physicochemical soil properties. The results of this study show remarkable changes in LULC in the period of 18?years. The effect of land cover change on soil properties, soil compaction and soil strength was found to be significant at a level of <0.05.     

The isolated non-circular ringlets of Saturn

We show that the combined effect of electrodynamic and gravitational forces can account for a number of features observed by Voyagers 1 and 2 in the isolated fine dust rings of Saturn. These include (a) the appearance and disappearnce of the braids in the F-ring, (b) the eccentricities of the F-ring and the ringlets within the Encke and Cassine divisions and a gap in the C-ring, and (c) the kinks in the eccentric Encke ring. They may also account for the very existence of these rings.