The importance of vertical velocity variability for estimates of the indirect aerosol effects

The activation of aerosols to form cloud droplets is dependent upon vertical velocities whose local variability is not typically resolved at the GCM grid scale. Consequently, it is necessary to represent the subgrid-scale variability of vertical velocity in the calculation of cloud droplet number concentration. This study uses the UK Chemistry and Aerosols community model (UKCA) within the Hadley Centre Global Environmental Model (HadGEM3), coupled for the first time to an explicit aerosol activation parameterisation, and hence known as UKCA-Activate. We explore the range of uncertainty in estimates of the indirect aerosol effects attributable to the choice of parameterisation of the subgrid-scale variability of vertical velocity in HadGEM-UKCA. Results of simulations demonstrate that the use of a characteristic vertical velocity cannot replicate results derived with a distribution of vertical velocities, and is to be discouraged in GCMs. This study focuses on the effect of the variance (σw2) of a Gaussian pdf (probability density function) of vertical velocity. Fixed values of σw (spanning the range measured in situ by nine flight campaigns found in the literature) and a configuration in which σw depends on turbulent kinetic energy are tested. Results from the mid-range fixed σw and TKE-based configurations both compare well with observed vertical velocity distributions and cloud droplet number concentrations. The radiative flux perturbation due to the total effects of anthropogenic aerosol is estimated at −1.9 W m−2 with σw = 0.1 m s−1, −2.1 W m−2 with σw derived from TKE, −2.25 W m−2 with σw = 0.4 m s−1, and −2.3 W m−2 with σw = 0.7 m s−1. The breadth of this range is 0.4 W m−2, which is comparable to a substantial fraction of the total diversity of current aerosol forcing estimates. Reducing the uncertainty in the parameterisation of σw would therefore be an important step towards reducing the uncertainty in estimates of the indirect aerosol effects. Detailed examination of regional radiative flux perturbations reveals that aerosol microphysics can be responsible for some climate-relevant radiative effects, highlighting the importance of including microphysical aerosol processes in GCMs.

Comment on “B_y fluctuations in the magnetosheath and azimuthal flow velocity transients in the dayside ionosphere” by Newell and Sibeck

Newell and Sibeck [1993] (hereafter N&S) list some objections to our interpretation of dayside auroral transients and associated azimuthal flow bursts in terms of pulsed reconnection [e.g. Lockwood et al., 1989; 1993a]. They present what they term an “apparently overlooked” alternative explanation in terms of steady reconnection and fluctuations in the magnitude of the By component of the magnetosheath field. The objections of N&S can all be answered by reference to our previous publications and their alternative explanation was only “overlooked” in so far as it fails to explain the observations. Here we discuss just some of the reasons why the objections of N&S are invalid, and then give reasons why the events are not simply due to magnetosheath |By| changes.

Dayside ionospheric convection changes in response to long-period interplanetary Magnetic field oscillations: Determination of the ionospheric phase velocity

Ground magnetic field perturbations recorded by the CANOPUS magnetometer network in the 7 to 13 MLT sector are used to examine how reconfigurations of the dayside polar ionospheric flow take place in response to north-south changes of the IMF. During the 6-hour interval in question IMF Bz oscillates between ±7 nT with about a 1-hour period. Corresponding variations in the ground magnetic disturbance are observed which we infer are due to changes in ionospheric flow. Cross correlation of the data obtained from two ground stations at 73.5° magnetic latitude, but separated by ∼2 hours in MLT, shows that changes in the flow are initiated in the prenoon sector (∼10 MLT) and then spread outward toward dawn and dusk with a phase speed of ∼5 km s−1 over the longitude range ∼8 to 12 MLT, slowing to ∼2 km s−1 outside this range. Cross correlating the data from these ground stations with IMP 8 IMF Bz records produces a MLT variation in the ground response delay relative to the IMF which is compatible with these deduced phase speeds. We interpret these observations in terms of the ionospheric response to the onset, expansion and decay of magnetic reconnection at the dayside magnetopause.

Incoherent scatter radar observations of non-Maxwellian ion velocity distributions in the auroral F-region

Observations by the EISCAT experiments “POLAR” and Common Programme CP-3 reveal non-Maxwellian ion velocity distributions in the auroral F-region ionosphere. Analysis of data from three periods is presented. During the first period, convection velocities are large (≈2 km s-1) and constant over part of a CP-3 latitude scan; the second period is one of POLAR data containing a short-lived (<1 min.) burst of rapid (>1.5 km s-1) flow. We concentrate on these two periods as they allow the study of a great many features of the ion-neutral interactions which drive the plasma non-thermal and provide the best available experimental test for models of the 3-dimensional ion velocity distribution function. The third period is included to illustrate the fact that non-thermal plasma frequently exists in the auroral ionosphere: the data, also from the POLAR experiment, cover a three-hour period of typical auroral zone flow and analysis reveals that the ion distribution varies from Maxwellian to the threshold of a toroidal form.

The effect of rapid changes in ionospheric flow on velocity vectors deduced from radar beam-swinging experiments

The effects on the horizontal ionospheric velocity vectors deduced from radar beam-swinging experiments, which occur when changes in the flow take place on short time scales compared with the experiment cycle time, are analysed in detail. The further complications which arise in the interpretation of beam-swinging data, due to longitudinal gradients in the flow and to field-aligned flows, are also considered. It is concluded that these effects are unlikely to seriously compromise statistical determinations of the response time of the flow, e.g. to changes in the north-south component of the IMF, such as have been recently reported by Etemadiet al. (1988, Planet. Space Sci.36, 471), using EISCAT ‘Polar’ data.

The influence of anisotropic F region ion velocity distributions on ionospheric ion outflows into the magnetosphere

The contribution to the field-aligned ionospheric ion momentum equation, due to coupling between pressure anisotropy and the inhomogeneous geomagnetic field, is investigated. We term this contribution the “hydrodynamic mirror force” and investigate its dependence on the ion drift and the resulting deformations of the ion velocity distribution function from an isotropic form. It is shown that this extra upforce increases rapidly with ion drift relative to the neutral gas but is not highly dependent on the ion-neutral collision model employed. An example of a burst of flow observed by EISCAT, thought to be the ionospheric signature of a flux transfer event at the magnetopause, is studied in detail and it is shown that the nonthermal plasma which results is subject to a hydrodynamic mirror force which is roughly 10% of the gravitational downforce. In addition, predictions by the coupled University College London-Sheffield University model of the ionosphere and thermosphere show that the hydrodynamic mirror force in the auroral oval is up to 3% of the gravitational force for Kp of about 3, rising to 10% following a sudden increase in cross-cap potential. The spatial distribution of the upforce shows peaks in the cusp region and in the post-midnight auroral oval, similar to that of observed low-energy heavy ion flows from the ionosphere into the magnetosphere. We suggest the hydrodynamic mirror force may modulate these outflows by controlling the supply of heavy ions to regions of ion acceleration and that future simulations of the effects of Joule heating on ion outflows should make allowance for it.

Non-Maxwellian ion velocity distributions observed using EISCAT

Recent observations from the EISCAT incoherent scatter radar have revealed bursts of poleward ion flow in the dayside auroral ionosphere which are consistent with the ionospheric signature of flux transfer events at the magnetopause. These bursts frequently contain ion drifts which exceed the neutral thermal speed and, because the neutral thermospheric wind is incapable of responding sufficiently rapidly, toroidal, non-Maxwellian ion velocity distributions are expected. The EISCAT observations are made with high time resolution (15 seconds) and at a large angle to the geomagnetic field (73.5°), allowing the non-Maxwellian nature of the distribution to be observed remotely for the first time. The observed features are also strongly suggestive of a toroidal distribution: characteristic spectral shape, increased scattered power (both consistent with reduced Landau damping and enhanced electric field fluctuations) and excessively high line-of-sight ion temperatures deduced if a Maxwellian distribution is assumed. These remote sensing observations allow the evolution of the distributions to be observed. They are found to be non-Maxwellian whenever the ion drift exceeds the neutral thermal speed, indicating that such distributions can exist over the time scale of the flow burst events (several minutes).

Affective trajectories: locating diegetic velocity in the cinema experience

Films that feature high-speed diegetic motion, and present those high speeds through fast mobile framing and fast cutting, are frequently charged with generating a sensory overload which empties out meaning or any sense of spatial orientation. Inherent in this discourse is a privileging of optical-spatial intelligibility that suppresses consideration of the ways cinema can represent diegetic velocity, and the spectator’s sensory experience of the same. This paper will instead highlight the centrality of the evocation of a trajectory for movement for the spectator’s experience of diegetic speed, an evocation that does not depend on optical-spatial legibility for its affective force.

DENSITY PROFILE, VELOCITY ANISOTROPY, AND LINE-OF-SIGHT EXTERNAL CONVERGENCE OF SLACS GRAVITATIONAL LENSES

Data from 58 strong-lensing events surveyed by the Sloan Lens ACS Survey are used to estimate the projected galaxy mass inside their Einstein radii by two independent methods: stellar dynamics and strong gravitational lensing. We perform a joint analysis of these two estimates inside models with up to three degrees of freedom with respect to the lens density profile, stellar velocity anisotropy, and line-of-sight (LOS) external convergence, which incorporates the effect of the large-scale structure on strong lensing. A Bayesian analysis is employed to estimate the model parameters, evaluate their significance, and compare models. We find that the data favor Jaffe`s light profile over Hernquist`s, but that any particular choice between these two does not change the qualitative conclusions with respect to the features of the system that we investigate. The density profile is compatible with an isothermal, being sightly steeper and having an uncertainty in the logarithmic slope of the order of 5% in models that take into account a prior ignorance on anisotropy and external convergence. We identify a considerable degeneracy between the density profile slope and the anisotropy parameter, which largely increases the uncertainties in the estimates of these parameters, but we encounter no evidence in favor of an anisotropic velocity distribution on average for the whole sample. An LOS external convergence following a prior probability distribution given by cosmology has a small effect on the estimation of the lens density profile, but can increase the dispersion of its value by nearly 40%.

Moho map of South America from receiver functions and surface waves

We estimate crustal structure and thickness of South America north of roughly 40 degrees S. To this end, we analyzed receiver functions from 20 relatively new temporary broadband seismic stations deployed across eastern Brazil. In the analysis we include teleseismic and some regional events, particularly for stations that recorded few suitable earthquakes. We first estimate crustal thickness and average Poisson`s ratio using two different stacking methods. We then combine the new crustal constraints with results from previous receiver function studies. To interpolate the crustal thickness between the station locations, we jointly invert these Moho point constraints, Rayleigh wave group velocities, and regional S and Rayleigh waveforms for a continuous map of Moho depth. The new tomographic Moho map suggests that Moho depth and Moho relief vary slightly with age within the Precambrian crust. Whether or not a positive correlation between crustal thickness and geologic age is derived from the pre-interpolation point constraints depends strongly on the selected subset of receiver functions. This implies that using only pre-interpolation point constraints (receiver functions) inadequately samples the spatial variation in geologic age. The new Moho map also reveals an anomalously deep Moho beneath the oldest core of the Amazonian Craton.

Seismic activity triggered by water wells in the Parana Basin, Brazil

Triggered seismicity is commonly associated with deep water reservoirs or injection wells where water is injected at high pressure into the reservoir rock. However, earth tremors related solely to the opening of groundwater wells are extremely rare. Here we present a clear case of seismicity induced by pore-pressure changes following the drilling of water wells that exploit a confined aquifer in the intracratonic Parana Basin of southeastern Brazil. Since 2004, shallow seismic activity, with magnitudes up to 2.9 and intensities V MM, has been observed near deep wells (120-200 m) that were drilled in early 2003 near the town of Bebedouro. The wells were drilled for irrigation purposes, cross a sandstone layer about 60-80 m thick and extract water from a confined aquifer in fractured zones between basalt flow layers. Seismic activity, mainly event swarms, has occurred yearly since 2004, mostly during the rainy season when the wells are not pumped. During the dry season when the wells are pumped almost continuously, the activity is very low. A seismographic network, installed in March 2005, has located more than 2000 microearthquakes. The events are less than 1 km deep (mostly within the 0.5 km thick basalt layer) and cover an area roughly 1.5 km x 5 km across. The seismicity generally starts in a small area and expands to larger distances with an equivalent hydraulic diffusivity ranging from 0.06 to 0.6 m(2)/s. Geophysical and geothermal logging of several wells in the area showed that water from the shallow sandstone aquifer enters the well at the top and usually forms waterfalls. The waterfalls flow down the sides of the wells and feed the confined, fractured aquifer in the basalt layer at the bottom. Two seismic areas are observed: the main area surrounds several wells that are pumped continuously during the dry season, and a second area near another well (about 10 km from the first area) that is not used for irrigation and not pumped regularly. The main area displays cyclic annual activity, but the second area does not. We explain the earthquake swarms as being triggered by pore pressure diffusion in the fractured basalt layer due to additional pressure from the newly connected surface aquifer. This reaches critically prestressed areas up to a few kilometers away from the wells. During periods of continuous pumping, the reduction of pore pressure in the confined aquifer stops the seismic activity. Our study suggests that this kind of activity may be more common than previously thought and implies that many other cases of small tremors associated with the drilling of water wells may have gone unnoticed.

The intracratonic Caraibas-Itacarambi earthquake of December 09, 2007 (4.9 m(b)), Minas Gerais State, Brazil

On December 9, 2007, a 4.9 m(b) earthquake occurred in the middle of the Sao Francisco Craton, in a region with no known previous activity larger than 4 m(b). This event reached intensity VII MM (Modified Mercalli) causing the first fatal victim in Brazil. The activity had started in May 25, 2007 with a 3.5 magnitude event and continued for several months, motivating the deployment of a local 6-station network. A three week seismic quiescence was observed before the mainshock. Initial absolute hypocenters were calculated with best fitting velocity models and then relative locations were determined with hypoDD. The aftershock distribution indicates a 3 km long rupture for the mainshock. The fault plane solution, based on P-wave polarities and hypocentral trend, indicates a reverse faulting mechanism on a N30 degrees E striking plane dipping about 40 degrees to the SE. The rupture depth extends from about 0.3 to 1.2 km only. Despite the shallow depth of the mainshock, no surface feature could be correlated with the fault plane. Aeromagnetic data in the epicentral area show short-wavelength lineaments trending NNE-SSW to NE-SW which we interpret as faults and fractures in the craton basement beneath the surface limestone layer. We propose that the Caraibas-Itacarambi seismicity is probably associated with reactivation of these basement fractures and faults under the present E-W compressional stress field in this region of the South American Plate. (c) 2009 Elsevier B.V. All rights reserved.

Deep crustal structure of the Parana Basin from receiver functions and Rayleigh-wave dispersion: Evidence for a fragmented cratonic root

The deep crustal structure of the Parana Basin of southern Brazil is investigated by analyzing P- and PP-wave receiver functions at 17 Brazilian Lithosphere Seismic Project stations within the basin. The study area can be described as a typical Paleozoic intracratonic basin that hosts one of the largest Large Igneous Province of the world and makes a unique setting for investigating models of basin subsidence and their interaction with mantle plumes. Our study consists of (1) an analysis of the Moho interaction phases in the receiver functions to obtain the thickness and bulk Vp/Vs ratio of the basin`s underlying crust and (2) a joint inversion with Rayleigh-wave dispersion velocities from an independent tomographic study to delineate the detailed S-wave velocity variation with depth. The results of our analysis reveal that Moho depths and bulk Vp/Vs ratios (including sediments) vary between 41 and 48 km and between 1.70 and 1.76, respectively, with the largest values roughly coinciding with the basin`s axis, and that S-wave velocities in the lower crust are generally below 3.8 km/s. Select sites within the basin, however, show lower crustal S-wave velocities slightly above 3.9 km/s suggestive of underplated mafic material. We show that these observations are consistent with a fragmented cratonic root under the Parana basin that defined a zone of weakness for the initial Paleozoic subsidence of the basin and which allowed localized mafic underplating of the crust along the suture zones by Cenozoic magmatism.

Genetic algorithm inversion of the average 1D crustal structure using local and regional earthquakes

Knowing the best 1D model of the crustal and upper mantle structure is useful not only for routine hypocenter determination, but also for linearized joint inversions of hypocenters and 3D crustal structure, where a good choice of the initial model can be very important. Here, we tested the combination of a simple GA inversion with the widely used HYPO71 program to find the best three-layer model (upper crust, lower crust, and upper mantle) by minimizing the overall P- and S-arrival residuals, using local and regional earthquakes in two areas of the Brazilian shield. Results from the Tocantins Province (Central Brazil) and the southern border of the Sao Francisco craton (SE Brazil) indicated an average crustal thickness of 38 and 43 km, respectively, consistent with previous estimates from receiver functions and seismic refraction lines. The GA + HYPO71 inversion produced correct Vp/Vs ratios (1.73 and 1.71, respectively), as expected from Wadati diagrams. Tests with synthetic data showed that the method is robust for the crustal thickness, Pn velocity, and Vp/Vs ratio when using events with distance up to about 400 km, despite the small number of events available (7 and 22, respectively). The velocities of the upper and lower crusts, however, are less well constrained. Interestingly, in the Tocantins Province, the GA + HYPO71 inversion showed a secondary solution (local minimum) for the average crustal thickness, besides the global minimum solution, which was caused by the existence of two distinct domains in the Central Brazil with very different crustal thicknesses. (C) 2010 Elsevier Ltd. All rights reserved.