Fabric evolution of granular assembly under K0 loading/unloading

This study attempted to investigate the fabric evolution in K₀ loading/unloading. The work made use of a field simulator to control K₀ loading/unloading in large specimens prepared by air‐pluviation. In each loading stage, wave velocities along various propagation directions were measured. On the basis of the theories of micro‐mechanics and wave propagation, the microscopic parameters of the granular assembly were back calculated to investigate the fabric evolution of granular soil during K₀ loading/unloading. In this study, the Geometric fabric was modelled by fabric tensors of ranks 2 and 4. The comparison of calibrated results using ranks 2 and 4 revealed the advantage of the usage of rank‐4 fabric tensor in modelling fabric evolution in spite of its complexity. By comparing relative magnitudes of vertical and horizontal components of geometric fabric, it was demonstrated that relative constraint in lateral directions increased during K₀‐unloading in order to maintain a K₀ condition. It revealed that fabric evolution was responsible for a higher K₀ in unloading than K₀ in loading. Copyright © 2003 John Wiley & Sons, Ltd.     

Observation of an Extremely Large-Density Heliospheric Plasma Sheet Compressed by an Interplanetary Shock at 1 AU

At 11:46 UT on 9 September 2011, the Wind spacecraft encountered an interplanetary (IP) fast-forward shock. The shock was followed almost immediately by a short-duration (～?35 minutes) extremely dense pulse (with a peak ～?94 cm^?3). The pulse induced an extremely large positive impulse (SYM-H = 74 nT and Dst = 48 nT) on the ground. A close examination of other in situ parameters from Wind shows that the density pulse was associated with i) a spike in the plasma \(\upbeta\) (ratio of thermal to magnetic pressure), ii) multiple sign changes in the azimuthal component of the magnetic field (\(B_{\phi}\)), iii) a depressed magnetic field magnitude, iv) a small radial component of the magnetic field, and v) a large (>?90°) change in the suprathermal (～?255 eV) electron pitch angle across the density pulse. We conclude that the density pulse is associated with the heliospheric plasma sheet (HPS). The thickness of the HPS is estimated to be \({\sim}\,8.2\times10^{5}\ \mbox{km}\). The HPS density peak is about five times the value of a medium-sized density peak inside the HPS (～?18 cm^?3) at 1 AU. Our global three-dimensional magnetohydrodynamic simulation results (Wu et al. in J. Geophys. Res. 212, 1839, 2016) suggest that the extremely large density pulse may be the result of the compression of the HPS by an IP shock crossing or an interaction between an interplanetary shock and a corotating interaction region.     

Ice microphysics and climatic temperature feedback

The potential effects of ice microphysics involving ice crystal size distribution and ice water path (IWP) on climatic temperature perturbations are investigated by using a one-dimensional radiative-turbulent climate model. We define a mean effective size, denoting the width of ice crystals weighted by the geometric cross section area, to represent ice crystal size distribution. Based on aircraft measurements, both the mean effective size and IWP are related to temperature and may be parameterized as functions of temperature. The radiative properties of cirrus clouds are further parameterized in terms of these two basic cloud physics parameters. Using CO₂ doubling as the radiative forcing, feedbacks among temperature, the mean effective size and IWP, and the radiative properties of clouds are analyzed from the model results. We show that overall, a positive feedback associated with ice microphysics and the coupled radiative transfer is produced by temperature increase.     

Protolith and metamorphic ages of the Haiyangsuo Complex, eastern China: a non-UHP exotic tectonic slab in the Sulu ultrahigh-pressure terrane

Summary The Haiyangsuo Complex in the NE Sulu ultrahigh-pressure (UHP) terrane has discontinuous, coastal exposures of Late Archean gneiss with amphibolitized granulite, amphibolite, Paleoproterozoic metagabbroic intrusives, and Cretaceous granitic dikes over an area of about 15 km². The U–Pb SHRIMP dating of zircons indicates that theprotolith age of a garnet-biotite gneiss is >2500 Ma, whereas the granulite-facie metamorphism occurred at around 1800 Ma. A gabbroic intrusion was dated at ∼1730 Ma, and the formation of amphibolite-facies assemblages in both metagabbro and granulite occurred at ∼340–460 Ma. Petrologic and geochronological data indicate that these various rocks show no evidence of Triassic eclogite-facies metamorphism and Neoproterozoic protolith ages that are characteristics of Sulu-Dabie HP-UHP rocks, except Neoproterozoic inherited ages from post-collisional Jurassic granitic dikes. Haiyangsuo retrograde granulites with amphibolite-facies assemblages within the gneiss preserve relict garnet formed during granulite-facies metamorphism at ∼1.85 Ga. The Paleoproterozoic metamorphic events are almost coeval with gabbroic intrusions. The granulite-bearing gneiss unit and gabbro-dominated unit of the Haiyangsuo Complex were intruded by thin granitic dikes at about 160 Ma, which is coeval with post-collisional granitic intrusions in the Sulu terrane. We suggest that the Haiyangsuo Complex may represent a fragment of the Jiao-Liao-Ji Paleoproterozoic terrane developed at the eastern margin of the Sino-Korean basement, which was juxtaposed with the Sulu terrane prior to Jurassic granitic activity and regional deformation.     

Transition of UHP eclogites to gneissic rocks of low-amphibolite facies during exhumation: evidence from the Dabie terrane, central China

The Shuanghe ultrahigh-pressure (UHP) slab in the Dabie Mountains consists of layered coesite-bearing eclogite, jadeite quartzite, marble and biotite gneiss, and is fault bounded against hosting orthogneiss. Representative assemblages of eclogite are Grt+Omp+Coe+Rt±Ky±Phn±Mgs; it formed at P>27 kbar and 680–720±50 °C. During exhumation, these UHP rocks experienced multistage retrograde metamorphism. Coesite was overprinted by quartz aggregates, phengite by biotite±muscovite and rutile by titanite. Garnet was successively replaced by a thin rim of Amp, Amp+Pl, and Amp+Ep±Bt+Pl (minor). Omphacite and kyanite were replaced by Amp+Pl±Cpx (or ±Bt) and by Zo+Pl+Ms±Mrg±Bt, respectively. Secondary calcite occurs as irregular pockets in some layers. An outcrop near the UHP slab border is composed of 20 thin, concordant layers of foliated eclogites, amphibolite and gneissic rocks of variable bulk composition. These layers exhibit mineral assemblages and textures transitional from less altered through extensively retrograded eclogite to gneissic rock of low-amphibolite facies through hydration, metasomatism and recrystallization. Retrograde metamorphism has caused oxygen and hydrogen isotope disequilibria between some of the minerals, but the fluid for retrograde reactions was internally buffered in the stable isotope compositions. Retrograde metamorphism of variable extent may be attributed to selective infiltration of retrograde fluids of CO₂-rich and low-salinity aqueous, intensity of deformation and mineral resistance to alteration. The fluid phase for retrogression may have occurred either as discontinuous flow along grain boundaries in completely retrograded eclogites, and/or as isolated pockets in extensive or less altered eclogite layers.     

Exploration of the remote sounding of infrared cooling rates due to water vapor

Summary We explore the feasibility of deriving atmospheric infrared cooling rates by direct inversion of radiances observed by satellites from space. In order to convert radiances to fluxes and achieve vertical profiling at the same time, we show that it is necessary to combine radiances from narrow channels with radiances averaged over spectral bands. We demonstrate that the vertical integral of the cooling rate in the spectral band, convolved with a kernel function associated with the narrow channel, can be related to a weighted sum of the channel and band radiances. The band radiance must be evaluated at a specific zenith angle, which is a result of use of the mean value theorem. With known kernel functions, the combined radiances may be inverted to obtain the cooling rate profile. These results are derived from use of a random model for the transmittance in its strong-and weak-line limits. The results are similar in the two limits leading us to conclude that there are expressions that are approximately valid over the entire range of transmittance. We show by numerical methods that this conclusion is correct and apply the retrieval technique successfully to get the cooling rate profile in the rotational band of water vapor.With 3 Figures     

Remote sounding of the cirrus optical depth and temperature from 3.7 and 11 micrometer windows

On the basis of forward IR radiation transfer analyses for an atmosphere containing semi-transparent, non-black cirrus, parameterization equations are derived for the retrieval of the cloud optical depth and cloud temperature utilizing AVHRR 3.7 and 10.8 μum channels. The retrieval techniques developed involve the use of either dual-frequency or dual-scanning angle radiance observations. We show that the cloud optical depth and cloud temperature may be inferred successively from the observed brightness temperature differences using these two techniques. Numerical experiments anderror analyses demonstrate that the dual-frequency method is specifically appropriate for optically thin cirrus cases (τ <1). In case the optical depth of cirrus is close to 3, combination of dual-frequency and dual-scanning angle methods is shown to give reasonable accuracy for the cloud optical depth and temperature retrieval.     

Gamma random field simulation by a covariance matrix transformation method

In studies involving environmental risk assessment, Gaussian random field generators are often used to yield realizations of a Gaussian random field, and then realizations of the non-Gaussian target random field are obtained by an inverse-normal transformation. Such simulation process requires a set of observed data for estimation of the empirical cumulative distribution function (ECDF) and covariance function of the random field under investigation. However, if realizations of a non-Gaussian random field with specific probability density and covariance function are needed, such observed-data-based simulation process will not work when no observed data are available. In this paper we present details of a gamma random field simulation approach which does not require a set of observed data. A key element of the approach lies on the theoretical relationship between the covariance functions of a gamma random field and its corresponding standard normal random field. Through a set of devised simulation scenarios, the proposed technique is shown to be capable of generating realizations of the given gamma random fields.     

Theoretical connections between optimization algorithms based on an approximate gradient

Performing a line search method in the direction given by the simplex gradient is a well-known method in the mathematical optimization community. For reservoir engineering optimization problems, both a modification of the simultaneous perturbation stochastic approximation (SPSA) and ensemble-based optimization (EnOpt) have recently been applied for estimating optimal well controls in the production optimization step of closed-loop reservoir management. The modified SPSA algorithm has also been applied to assisted history-matching problems. A recent comparison of the performance of EnOpt and a SPSA-type algorithm (G-SPSA) for a set of production optimization test problems showed that the two algorithms resulted in similar estimates of the optimal net-present-value and required roughly the same amount of computational time to achieve these estimates. Here, we show that, theoretically, this result is not surprising. In fact, we show that both the simplex, preconditioned simplex, and EnOpt algorithms can be derived directly from a modified SPSA-type algorithm where the preconditioned simplex algorithm is presented for the first time in this paper. We also show that the expectation of all these preconditioned stochastic gradients is a first-order approximation of the preconditioning covariance matrix times the true gradient or a covariance matrix squared times the true gradient.     

On the radiative properties of ice clouds: Light scattering,remote sensing,and radiation parameterization

Presented is a review of the radiative properties of ice clouds from three perspectives: light scattering simulations, remote sensing applications, and broadband radiation parameterizations appropriate for numerical models. On the subject of light scattering simulations, several classical computational approaches are reviewed, including the conventional geometric-optics method and its improved forms, the finite-difference time domain technique, the pseudo-spectral time domain technique, the discrete dipole approximation method, and the T-matrix method, with specific applications to the computation of the singlescattering properties of individual ice crystals. The strengths and weaknesses associated with each approach are discussed.With reference to remote sensing, operational retrieval algorithms are reviewed for retrieving cloud optical depth and effective particle size based on solar or thermal infrared(IR) bands. To illustrate the performance of the current solar- and IR-based retrievals, two case studies are presented based on spaceborne observations. The need for a more realistic ice cloud optical model to obtain spectrally consistent retrievals is demonstrated. Furthermore, to complement ice cloud property studies based on passive radiometric measurements, the advantage of incorporating lidar and/or polarimetric measurements is discussed.The performance of ice cloud models based on the use of different ice habits to represent ice particles is illustrated by comparing model results with satellite observations. A summary is provided of a number of parameterization schemes for ice cloud radiative properties that were developed for application to broadband radiative transfer submodels within general circulation models(GCMs). The availability of the single-scattering properties of complex ice habits has led to more accurate radiation parameterizations. In conclusion, the importance of using nonspherical ice particle models in GCM simulations for climate studies is proven.