Bayesian Markov-chain-Monte-Carlo inversion of time-lapse cross hole ground-penetrating radar data to characterize the vadose zone at the Arrenaes field site, Denmark

The ground-penetrating radar (GPR) geophysical method has the potential to provide valuable information on the hydraulic properties of the vadose zone because of its strong sensitivity to soil water content. In particular, recent evidence has suggested that the stochastic inversion of crosshole GPR traveltime data can allow for a significant reduction in uncertainty regarding subsurface van Genuchten-Mualem (VGM) parameters. Much of the previous work on the stochastic estimation of VGM parameters from crosshole GPR data has considered the case of steady-state infiltration conditions, which represent only a small fraction of practically relevant scenarios. We explored in detail the dynamic infiltration case, specifically examining to what extent time-lapse crosshole GPR traveltimes, measured during a forced infiltration experiment at the Arreneas field site in Denmark, could help to quantify VGM parameters and their uncertainties in a layered medium, as well as the corresponding soil hydraulic properties. We used a Bayesian Markov-chain-Monte-Carlo inversion approach. We first explored the advantages and limitations of this approach with regard to a realistic synthetic example before applying it to field measurements. In our analysis, we also considered different degrees of prior information. Our findings indicate that the stochastic inversion of the time-lapse GPR data does indeed allow for a substantial refinement in the inferred posterior VGM parameter distributions compared with the corresponding priors, which in turn significantly improves knowledge of soil hydraulic properties. Overall, the results obtained clearly demonstrate the value of the information contained in time-lapse GPR data for characterizing vadose zone dynamics

When threat matters: Self-regulation, threat salience, and stereotyping

The final publication is available at Elsevier via http://dx.doi.org/10.1016/j.jesp.2015.03.003 © 2015. This manuscript version is made available under the CC-BY-NC-ND 4.0 license http://creativecommons.org/licenses/by-nc-nd/4.0/Four experiments examined whether information implying imminent threat to safety would interact with regulatory focus (Higgins, 1997) to affect the utilization of threat-relevant stereotypes. Because information suggesting imminent danger is more relevant to the safety goals of prevention-focused individuals than the advancement goals of promotion-focused individuals, utilization of threat-relevant stereotypes was expected to increase under such conditions only under prevention focus. Support for this prediction was obtained in four distinct and socially important domains. Using scenarios describing a violent crime committed by an African-American male (Experiment 1) or a petty crime committed by an undocumented immigrant (Experiment 2), prevention-focused individuals made judgments consistent with stereotypes when threat was perceived to be high rather than low. In studies that manipulated the stereotypicality of the target in a terrorism scenario (Experiments 3 & 4), prevention-focused individuals were more likely to endorse scrutinizing a stereotypical compared with a non-stereotypical target when terrorism was described as an increasing problem. Implications for models of stereotyping, self-regulation, and responding to threat are discussed.NSF Grant [1147779

Fast magnetic reconnection in free space: self-similar evolution process

We present a new model for time evolution of fast magnetic reconnection in free space, which is characterized by self-similarity. Reconnection triggered by locally enhanced resistivity assumed at the center of the current sheet can self-similarly and unlimitedly evolve until external factors affect the evolution. The possibility and stability of this type of evolution are verified by numerical simulations in a very wide spatial dynamic range. Actual astrophysical reconnection in solar flares and geomagnetospheric substorms can be treated as an evolutionary process in free space, because the resultant scale is much larger than the initial scale. In spite of this fact, most of the previous numerical works focused on the evolutionary characters strongly affected by artificial boundary conditions on the simulation boundary. Our new model clarifies a realistic evolution for such cases. The characteristic structure around the diffusion region is quite similar to the Petschek model which is characterized by a pair of slow-mode shocks and the fast-mode rarefaction-dominated inflow. However, in the outer region, a vortex-like return flow driven by the fast-mode compression caused by the piston effect of the plasmoid takes place. The entire reconnection system expands self-similarly.Comment: 17 Pages, 17 Figure

Estimation of the correlation structure of crustal velocity heterogeneity from seismic reflection data

Numerous sources of evidence point to the fact that heterogeneity within the Earth's deep crystalline crust is complex and hence may be best described through stochastic rather than deterministic approaches. As seismic reflection imaging arguably offers the best means of sampling deep crustal rocks in situ, much interest has been expressed in using such data to characterize the stochastic nature of crustal heterogeneity. Previous work on this problem has shown that the spatial statistics of seismic reflection data are indeed related to those of the underlying heterogeneous seismic velocity distribution. As of yet, however, the nature of this relationship has remained elusive due to the fact that most of the work was either strictly empirical or based on incorrect methodological approaches. Here, we introduce a conceptual model, based on the assumption of weak scattering, that allows us to quantitatively link the second-order statistics of a 2-D seismic velocity distribution with those of the corresponding processed and depth-migrated seismic reflection image. We then perform a sensitivity study in order to investigate what information regarding the stochastic model parameters describing crustal velocity heterogeneity might potentially be recovered from the statistics of a seismic reflection image using this model. Finally, we present a Monte Carlo inversion strategy to estimate these parameters and we show examples of its application at two different source frequencies and using two different sets of prior information. Our results indicate that the inverse problem is inherently non-unique and that many different combinations of the vertical and lateral correlation lengths describing the velocity heterogeneity can yield seismic images with the same 2-D autocorrelation structure. The ratio of all of these possible combinations of vertical and lateral correlation lengths, however, remains roughly constant which indicates that, without additional prior information, the aspect ratio is the only parameter describing the stochastic seismic velocity structure that can be reliably recovered

Fast Collisionless Reconnection Condition and Self-Organization of Solar Coronal Heating

I propose that solar coronal heating is a self-regulating process that keeps the coronal plasma roughly marginally collisionless. The self-regulating mechanism is based on the interplay of two effects. First, plasma density controls coronal energy release via the transition between the slow collisional Sweet-Parker regime and the fast collisionless reconnection regime. This transition takes place when the Sweet--Parker layer becomes thinner than the characteristic collisionless reconnection scale. I present a simple criterion for this transition in terms of the upstream plasma density (n_e), the reconnecting (B_0) and guide (B_z) magnetic field components, and the global length (L) of the reconnection layer: L < 6.10^9 cm [n_e/(10^{10}/cm^3)]^(-3) (B_0/30G)^4 (B_0/B_z)^2. Next, coronal energy release by reconnection raises the ambient plasma density via chromospheric evaporation and this, in turn, temporarily inhibits subsequent reconnection involving the newly-reconnected loops. Over time, however, radiative cooling gradually lowers the density again below the critical value and fast reconnection again becomes possible. As a result, the density is highly inhomogeneous and intermittent but, statistically, does not deviate strongly from the critical value which is comparable with the observed coronal density. Thus, in the long run, the coronal heating process can be represented by repeating cycles that consist of fast reconnection events (i.e., nanoflares), followed by rapid evaporation episodes, followed by relatively long periods (1-hour) during which magnetic stresses build up and simultaneously the plasma cools down and precipitates.Comment: 17 pages, no figures; accepted to the Astrophysical Journal; replaced to match the accepted versio

Simulations of Electron Acceleration at Collisionless Shocks: The Effects of Surface Fluctuations

Energetic electrons are a common feature of interplanetary shocks and planetary bow shocks, and they are invoked as a key component of models of nonthermal radio emission, such as solar radio bursts. A simulation study is carried out of electron acceleration for high Mach number, quasi-perpendicular shocks, typical of the shocks in the solar wind. Two dimensional self-consistent hybrid shock simulations provide the electric and magnetic fields in which test particle electrons are followed. A range of different shock types, shock normal angles, and injection energies are studied. When the Mach number is low, or the simulation configuration suppresses fluctuations along the magnetic field direction, the results agree with theory assuming magnetic moment conserving reflection (or Fast Fermi acceleration), with electron energy gains of a factor only 2 - 3. For high Mach number, with a realistic simulation configuration, the shock front has a dynamic rippled character. The corresponding electron energization is radically different: Energy spectra display: (1) considerably higher maximum energies than Fast Fermi acceleration; (2) a plateau, or shallow sloped region, at intermediate energies 2 - 5 times the injection energy; (3) power law fall off with increasing energy, for both upstream and downstream particles, with a slope decreasing as the shock normal angle approaches perpendicular; (4) sustained flux levels over a broader region of shock normal angle than for adiabatic reflection. All these features are in good qualitative agreement with observations, and show that dynamic structure in the shock surface at ion scales produces effective scattering and can be responsible for making high Mach number shocks effective sites for electron acceleration.Comment: 26 pages, 12 figure

Analytic solution for nonlinear shock acceleration in the Bohm limit

The selfconsistent steady state solution for a strong shock, significantly modified by accelerated particles is obtained on the level of a kinetic description, assuming Bohm-type diffusion. The original problem that is commonly formulated in terms of the diffusion-convection equation for the distribution function of energetic particles, coupled with the thermal plasma through the momentum flux continuity equation, is reduced to a nonlinear integral equation in one variable. Its solution provides selfconsistently both the particle spectrum and the structure of the hydrodynamic flow. A critical system parameter governing the acceleration process is found to be $\Lambda = M^{-3/4}\Lambda_1$, where $\Lambda_1 =\eta p_1/mc$, with a suitably normalized injection rate $\eta$, the Mach number M >> 1, and the cut-off momentum $p_1$. We particularly focus on an efficient solution, in which almost all the energy of the flow is converted into a few energetic particles. It was found that (i) for this efficient solution (or, equivalently, for multiple solutions) to exist, the parameter $\zeta =\eta\sqrt{p_0 p_1}/mc$ must exceed a critical value $\zeta_{cr} \sim 1$ ($p_0$ is the injection momentum), (ii) the total shock compression ratio r increases with M and saturates at a level that scales as $ r \propto \Lambda_1 (iii) the downstream power-law spectrum has the universal index q=3.5 over a broad momentum range. (iv) completely smooth shock transitions do not appear in the steady state kinetic description.Comment: 39 pages, 3 PostScript figures, uses aasms4.sty, to appear in Aug. 20, 1997 issue ApJ, vol. 48