The effect of orbital evolution on the Haumea (2003 EL61) collisional family

The Haumea family is currently the only identified collisional family in the Kuiper belt. We numerically simulate the long-term dynamical evolution of the family to estimate a lower limit of the family’s age and to assess how the population of the family and its dynamical clustering are preserved over Gyr timescales. We find that the family is not younger than 100 Myr, and its age is at least 1 Gyr with 95% confidence. We find that for initial velocity dispersions of 50–400 m s^?1, approximately 20–45% of the family members are lost to close encounters with Neptune after 3.5 Gyr of orbital evolution. We apply these loss rates to two proposed models for the formation of the Haumea family, a graze-and-merge type collision between two similarly sized, differentiated KBOs or the collisional disruption of a satellite orbiting Haumea. For the graze-and-merge collision model, we calculate that >85% of the expected mass in surviving family members within 150 m s^?1 of the collision has been identified, but that one to two times the mass of the known family members remains to be identified at larger velocities. For the satellite-break-up model, we estimate that the currently identified family members account for ～50% of the expected mass of the family. Taking observational incompleteness into account, the observed number of Haumea family members is consistent with either formation scenario at the 1σ level, however both models predict more objects at larger relative velocities (>150 m s^?1) than have been identified.     

The planetary fourier spectrometer (PFS) onboard the European Venus Express mission

The planetary fourier spectrometer (PFS) for the Venus Express mission is an infrared spectrometer optimized for atmospheric studies. This instrument has a short wavelength (SW) channel that covers the spectral range from 1700 to 11400 cm⁻¹ (0.9–5.5 μm) and a long wavelength (LW) channel that covers 250–1700 cm⁻¹ (5.5–45 μm). Both channels have a uniform spectral resolution of 1.3 cm⁻¹. The instrument field of view FOV is about 1.6 ° (FWHM) for the short wavelength channel and 2.8 ° for the LW channel which corresponds to a spatial resolution of 7 and 12 km when Venus is observed from an altitude of 250 km. PFS can provide unique data necessary to improve our knowledge not only of the atmospheric properties but also surface properties (temperature) and the surface-atmosphere interaction (volcanic activity).PFS works primarily around the pericentre of the orbit, only occasionally observing Venus from larger distances. Each measurements takes 4.5 s, with a repetition time of 11.5 s. By working roughly 1.5 h around pericentre, a total of 460 measurements per orbit will be acquired plus 60 for calibrations. PFS is able to take measurements at all local times, enabling the retrieval of atmospheric vertical temperature profiles on both the day and the night side.The PFS measures a host of atmospheric and surface phenomena on Venus. These include the:(1) thermal surface flux at several wavelengths near 1 μm, with concurrent constraints on surface temperature and emissivity (indicative of composition); (2) the abundances of several highly-diagnostic trace molecular species; (3) atmospheric temperatures from 55 to 100 km altitude; (4) cloud opacities and cloud-tracked winds in the lower-level cloud layers near 50-km altitudes; (5) cloud top pressures of the uppermost haze/cloud region near 70–80 km altitude; and (6) oxygen airglow near the 100 km level. All of these will be observed repeatedly during the 500-day nominal mission of Venus Express to yield an increased understanding of meteorological, dynamical, photochemical, and thermo-chemical processes in the Venus atmosphere. Additionally, PFS will search for and characterize current volcanic activity through spatial and temporal anomalies in both the surface thermal flux and the abundances of volcanic trace species in the lower atmosphere.Measurement of the 15 μm CO₂ band is very important. Its profile gives, by means of a complex temperature profile retrieval technique, the vertical pressure-temperature relation, basis of the global atmospheric study.PFS is made of four modules called O, E, P and S being, respectively, the interferometer and proximity electronics, the digital control unit, the power supply and the pointing device.     

Study of large solar energetic particle events with halo coronal mass ejections and their associated solar flares

The purpose of the present study is to investigate the association of solar energetic particle (SEP) events with halo coronal mass ejections (CME) and with their associated solar flares during the period 1997–2014 (solar cycle 23 and 24). We have found that halo CMEs are more effective in producing SEP events. The occurrence probability and peak fluxes of SEPs strongly depend on the halo CMEs speed (V) as follows. The highest associations, 56% for occurrence probability and 90% for average peak fluxes, are found for the halo CMEs with V> 1400 km s⁻¹ but the lowest associations, 20% for occurrence probability and 5% for average peak fluxes, are found for halo CMEs with speed range 600 ≤ V ≤ 1000 km s⁻¹. We have also examined the relationship between SEP events and halo CME associated solar flares and found that 73% of events are associated with western solar flares while only 27% are with eastern solar flares. For longitudinal study, 0–20° belt is found to be more dominant for the SEP events. The association of SEP events with latitudinal solar flares is also examined in the study. 51% of events are associated with those halo CMEs associated solar flares which occur in the southern hemisphere of the Sun while 49% are with those solar flares that occur in the northern hemisphere of the Sun. Also, 10–20° latitudinal belt is found to be likely associated with the SEP events. Further, 45% of SEP events are associated with M-class solar flares while 44% and 11% are with X and C-class respectively. Maximum number of SEP events are found for the fast halo CME associated X- class solar flares (68%) than M and C- class solar flares.     

Descent motions of the Huygens probe as measured by the Surface Science Package (SSP): Turbulent evidence for a cloud layer

The Huygens probe underwent vigorous short-period motions during its parachute descent through the atmosphere of Saturn's moon Titan in January 2005, at least some of which were excited by the Titan environment. Several sensors in the Huygens Surface Science Package (SSP) detect these motions, indicating the transition to the smaller stabilizer parachute, the changing probe spin rate, aerodynamic buffeting, and pendulum motions. Notably, in an altitude range of about 20–30 km where methane drops will freeze, the frequency content and statistical kurtosis of the tilt data indicate excitation by turbulent air motions like those observed in freezing clouds on Earth, supporting the suggestion of Tokano et al. [Tokano, T., McKay, C.P., Neubauer, F.M., Atreya, S.K., Ferri, F., Fulchignoni, M., Niemann, H.B. (2006a). Methane drizzle on Titan. Nature 442, 432–435] that the probe passed through such a cloud layer. Motions are weak below 20 km, suggesting a quiescent lower atmosphere with turbulent fluctuations of nominally <0.15 m/s (to within a factor of ∼2) but more violent motions in the upper troposphere may have been excited by turbulent winds with amplitudes of 1–2 m/s. Descent in part of the stratosphere (150–120 km) was smooth despite strong ambient wind (∼100 m/s), but known anomalies in the probe spin prevent investigation of turbulence in the known wind-shear layer from 60 to 100 km.     

Interpretation of combined infrared,submillimeter, and millimeter thermal flux data obtained during the Rosetta fly-by of Asteroid (21) Lutetia

The European Space Agency’s Rosetta spacecraft is the first Solar System mission to include instrumentation capable of measuring planetary thermal fluxes at both near-IR (VIRTIS) and submillimeter–millimeter (smm–mm, MIRO) wavelengths. Its primary mission is a 1 year reconnaissance of Comet 67P/Churyumov–Gerasimenko beginning in 2014. During a 2010 close fly-by of Asteroid 21 Lutetia, the VIRTIS and MIRO instruments provided complementary data that have been analyzed to produce a consistent model of Lutetia’s surface layer thermal and electrical properties, including a physical model of self-heating. VIRTIS dayside measurements provided highly resolved 1 K accuracy surface temperatures that required a low thermal inertia, I < 30 J/(K m² s^0.5). MIRO smm and mm measurements of polar night thermal fluxes produced constraints on Lutetia’s subsurface thermal properties to depths comparable to the seasonal thermal wave, yielding a model of I < 20 J/(K m² s^0.5) in the upper few centimeters, increasing with depth in a manner very similar to that of Earth’s Moon. Subsequent MIRO-based model predictions of the dayside surface temperatures reveal negative offsets of ～5–30 K from the higher VIRTIS-measurements. By adding surface roughness in the form of 50% fractional coverage of hemispherical mini-craters to the MIRO-based thermal model, sufficient self-heating is produced to largely remove the offsets relative to the VIRTIS measurements and also reproduce the thermal limb brightening features (relative to a smooth surface model) seen by VIRTIS. The Lutetia physical property constraints provided by the VIRTIS and MIRO data sets demonstrate the unique diagnostic capabilities of combined infrared and submillimeter/millimeter thermal flux measurements.     

Solar modulation of galactic cosmic rays during 19–23 solar cycles

We make a detailed analysis of cross-correlation and time-lag between monthly data of galactic cosmic rays (GCRs) intensity and different solar activity indices (e.g., sunspot number, sunspot area, green coronal Fe line and 10.7 cm solar radio flux) during 19–23 solar cycles. GCRs time-series data from Kiel neutron monitor station and solar data from the last 50 years period, covering five solar cycles (19–23), and alternating solar polarity states (i.e., five A < 0 and four A > 0) have been investigated. We find a clear asymmetry in the cross-correlation between GCRs and solar activity indicators for both odd and even-numbered solar cycles. The time-lags between GCRs and solar parameters are found different in different solar cycles as well as in the opposite polarity states (A < 0 and A > 0) within the same solar cycle. Possible explanations of the observed results are discussed in light of modulation models, including drift effects.     

Characteristics of the Expansion of Magnetic Funnels in Solar Quiet Regions

Via the potential field extrapolation of the observed photospheric magnetic field, the structure of the photospheric magnetic fields above solar quiet regions is renewed. As revealed by the result, below 20 Mm the open magnetic lines exhibit many obvious small funnel structures. These funnels expand with height and at the height of about 20 Mm they combine into large funnel structures. By a systematic study of the tendency of change of the cross section areas of funnels, it is discovered that the cross section areas of funnels in solar quiet regions expand approximately linearly. The velocity of expansion of magnetic funnels at rather low altitudes (< 20 Mm) is larger than that at high altitudes (> 20 Mm). This phenomenon has important significance for the two-dimensional numerical simulations of the origin of solar wind and the mass flow in magnetic loops. At the same time it is found that the number of closed magnetic lines decreases in the form of exponential function.     

The O2 nightglow in the martian atmosphere by SPICAM onboard of Mars-Express

We present observations of the O₂(a¹Δ_g) nightglow at 1.27 μm on Mars using the SPICAM IR spectrometer onboard of the Mars Express orbiter. In contrast to the O₂(a¹Δ_g) dayglow that results from the ozone photodissociation, the O₂(a¹Δ_g) nightglow is a product of the recombination of O atoms formed by CO₂ photolysis on the dayside at altitudes higher than 80 km and transported downward above the winter pole by the Hadley circulation. The first detections of the O₂(a¹Δ_g) nightglow in 2010 indicate that it is about two order of magnitude less intense than the dayglow (Bertaux, J.-L., Gondet, B., Bibring, J.-P., Montmessin, F., Lefèvre, F. [2010]. Bull. Am. Astron. Soc. 42, 1040; Clancy et al. [2010]. Bull. Am. Astron. Soc. 42, 1041). SPICAM IR sounds the martian atmosphere in the near-IR range (1–1.7 μm) with the spectral resolution of 3.5 cm^?1 in nadir, limb and solar occultation modes. In 2010 the vertical profiles of the O₂(a¹Δ_g) nightside emission have been obtained near the South Pole at latitudes of 82–83°S for two sequences of observations: L_s = 111–120° and L_s = 152–165°. The altitude of the emission maximum varied from 45 km on L_s = 111–120° to 38–49 km on L_s = 152–165°. Averaged vertically integrated intensity of the emission at these latitudes has shown an increase from 0.22 to 0.35 MR. Those values of total vertical emission rate are consistent with the OMEGA observations on Mars-Express in 2010. The estimated density of oxygen atoms at altitudes from 50 to 65 km varies from 1.5 × 10¹¹ to 2.5 × 10¹¹ cm^?3. Comparison with the LMD general circulation model with photochemistry (Lefèvre, F., Lebonnois, S., Montmessin, F., Forget, F. [2004]. J. Geophys. Res. 109, E07004; Lefèvre et al. [2008]. Nature 454, 971–975) shows that the model reproduces fairly well the O₂(a¹Δ_g) emission layer observed by SPICAM when the large field of view (>20 km on the limb) of the instrument is taken into account.     

Mercury’s radius change estimates revisited using MESSENGER data

Mapping tectonic features using MESSENGER data mainly acquired at high Sun incidence angle (>50°) reveals previously undetected structures. The analysis of the latter features determines an upward revision of measurements of density and spatial distribution of tectonism and thus of estimates of average contractional strain and planetary radius decrease. We calculated an average surface contraction of ～0.23–0.30% (～0.28% for fault dip angle θ = 30°) within an area corresponding to 21% of the planet. This strain, extrapolated to the entire surface, corresponds to a decrease in radius of about 2.4–3.6 km (～3.0 km for θ = 30°). These values are three–four times higher with respect to previous estimates and are compatible with results from thermomechanical models.     

Circulation of the Venus upper mesosphere/lower thermosphere: Doppler wind measurements from 2001–2009 inferior conjunction,sub-millimeter CO absorption line observations

Sub-millimeter ¹²CO (346 GHz) and ¹³CO (330 GHz) line absorptions, formed within the mesospheric to lower thermospheric altitude (70–120 km) region of the Venus atmosphere, have been mapped across the nightside disk of Venus during 2001–2009 inferior conjunctions, employing the James Clerk Maxwell Telescope (JCMT). Radiative transfer analysis of these thermal line absorptions supports temperature and CO mixing profile retrievals, as described in a companion paper (Clancy et al., 2012). Here, we consider the analysis of the sharp line absorption cores of these CO spectra in terms of accurate Doppler wind profile measurements at 95–115 km altitudes versus local time (～8 pm–4 am) and latitude (～60N–60S). These Doppler wind measurements support determinations of the nightside zonal and subsolar-to-antisolar (SSAS) circulation components over a variety of timescales. The average behavior fitted from 21 retrieved maps of ¹²CO Doppler winds (obtained over hourly, daily, weekly, and interannual intervals) indicates stronger average zonal (85 m/s retrograde) versus SSAS (65 m/s) circulation at the 1 μbar pressure (108–110 km altitude) level. However, the absolute and relative magnitudes of these circulation components exhibit extreme variability over daily to weekly timescales. Furthermore, the individual Doppler wind measurements within each nightside mapping observation generally show significant deviations (20–50 m/s, averaged over 5000 km horizontal scales) from the simple zonal/SSAS solution, with distinct local time and latitudinal characters that are also time variable. These large scale residual circulations contribute 30–70% of the observed nightside Doppler winds at any given time, and may be most responsible for global variations in nightside lower thermospheric trace composition and temperatures, as coincidentally retrieved CO abundance and temperature distributions do not correlate with solution retrograde zonal and SSAS winds (see companion paper, Clancy et al., 2012). Limited comparisons of these nightside submillimeter results with dayside infrared Doppler wind measurements suggest distinct dayside versus nightside circulations, in terms of zonal winds in particular. Combined ¹²CO and ¹³CO Doppler wind mapping observations obtained since 2004 indicate that the average zonal and SSAS wind components increase by 50–100% between altitudes of 100 and 115 km. If gravity waves originating from the cloud levels are responsible for the extension of zonal winds into the thermosphere (Alexander, M.J. [1992]. Geophys. Res. Lett. 19, 2207–2210), such waves deposit substantial momentum (i.e., break) in the lower nightside thermosphere.     

An investigation of the field-aligned currents associated with a large-scale ULF wave in the morning sector

Previous work by Scoffield, H.C., Yeoman, T.K., Wright, D.M., Milan, S.E., Wright, A.N., Strangeway, R.J. [2005. An investigation of the field aligned currents associated with a large scale ULF wave using data from CUTLASS and FAST. Ann. Geophys. 23, 487–498) investigated a large-scale ULF wave, occurring in the dusk sector (∼1900 MLT). The wave had a period of ∼800 s (corresponding to 1.2 mHz frequency), an azimuthal wave number of ∼7 and a full-width at half-maximum (FWHM) across the resonance of 350 km. IMAGE ground magnetometer and SuperDARN radar observations of the wave's spatial and temporal characteristics were used to parameterise a simple, two-dimensional field line resonance (FLR) model. The model-calculated field-aligned current (FAC) was compared with FACs derived from the FAST energetic particle spectra and magnetic field measurement. Here the authors use the same method to investigate the FAC structure of a second large-scale ULF wave, with a period of ∼450 s, occurring the dawn sector (∼0500 MLT) with an opposite sense background region 1–region 2 current system. This wave has a much larger longitudinal scale (m∼4.5) and a smaller latitude scale (FWHM=150 km). Unlike the dusk sector wave, which was dominated by upward FAC, FAST observations of the dawn sector wave show an interval of large-scale downward FAC of ∼1.5 μA m⁻². Downgoing magnetospheric electrons with energies of a few keV were observed, which are associated with upward FACs of ∼1 μA m⁻². For both wave studies, downward currents appear to be carried partially by upgoing electrons below the FAST energy detection threshold (5 eV), but also consist of a mixture of hotter downgoing magnetospheric electrons and upgoing ionospheric electrons of energies 30 eV–1 keV. Strong intervals of upward current show that small-scale structuring of scale ∼50 km has been imposed on the current carriers. In general, this study confirms the findings of Scoffield, H.C., Yeoman, T.K., Wright, D.M., Milan, S.E., Wright, A.N., Strangeway, R.J. [2005. An investigation of the FACs associated with a large-scale ULF wave using data from CUTLASS and FAST. Ann. Geophys. 23, 487–498).     

Water vapor near the cloud tops of Venus from Venus Express/VIRTIS dayside data

Observations of the dayside of Venus performed by the high spectral resolution channel (–H) of the Visible and Infrared Thermal Imaging Spectrometer (VIRTIS) on board the ESA Venus Express mission have been used to measure the altitude of the cloud tops and the water vapor abundance around this level with a spatial resolution ranging from 100 to 10 km. CO₂ and H₂O bands between 2.48 and 2.60 μm are analyzed to determine the cloud top altitude and water vapor abundance near this level. At low latitudes (±40°) mean water vapor abundance is equal to 3 ± 1 ppm and the corresponding cloud top altitude at 2.5 μm is equal to 69.5 ± 2 km. Poleward from middle latitudes the cloud top altitude gradually decreases down to 64 km, while the average H₂O abundance reaches its maximum of 5 ppm at 80° of latitude with a large scatter from 1 to 15 ppm. The calculated mass percentage of the sulfuric acid solution in cloud droplets of mode 2 (～1 μm) particles is in the range 75–83%, being in even more narrow interval of 80–83% in low latitudes. No systematic correlation of the dark UV markings with the cloud top altitude or water vapor has been observed.     

Photometric,spectral phase and temperature effects on 4 Vesta and HED meteorites: Implications for the Dawn mission

Phase angle and temperature are two important parameters that affect the photometric and spectral behavior of planetary surfaces in telescopic and spacecraft data. We have derived photometric and spectral phase functions for the Asteroid 4 Vesta, the first target of the Dawn mission, using ground-based telescopes operating at visible and near-infrared wavelengths (0.4–2.5 μm). Photometric lightcurve observations of Vesta were conducted on 15 nights at a phase angle range of 3.8–25.7° using duplicates of the seven narrowband Dawn Framing Camera filters (0.4–1.0 μm). Rotationally resolved visible (0.4–0.7 μm) and near-IR spectral observations (0.7–2.5 μm) were obtained on four nights over a similar phase angle range. Our Vesta photometric observations suggest the phase slope is between 0.019 and 0.029 mag/deg. The G parameter ranges from 0.22 to 0.37 consistent with previous results (e.g., Lagerkvist, C.-I., Magnusson, P., Williams, I.P., Buontempo, M.E., Argyle, R.W., Morrison, L.V. [1992]. Astron. Astrophys. Suppl. Ser. 94, 43–71; Piironen, J., Magnusson, P., Lagerkvist, C.-I., Williams, I.P., Buontempo, M.E., Morrison, L.V. [1997]. Astron. Astrophys. Suppl. Ser. 121, 489–497; Hasegawa, S. et al. [2009]. Lunar Planet. Sci. 40. ID 1503) within the uncertainty. We found that in the phase angle range of 0° < α ? 25° for every 10° increase in phase angle Vesta’s visible slope (0.5–0.7 μm) increases 20%, Band I and Band II depths increase 2.35% and 1.5% respectively, and the BAR value increase 0.30. Phase angle spectral measurements of the eucrite Moama in the lab show a decrease in Band I and Band II depths and BAR from the lowest phase angle 13° to 30°, followed by possible small increases up to 90°, and then a dramatic drop between 90° and 120° phase angle. Temperature-induced spectral effects shift the Band I and II centers of the pyroxene bands to longer wavelengths with increasing temperature. We have derived new correction equations using a temperature series (80–400 K) of HED meteorite spectra that will enable interpretation of telescopic and spacecraft spectral data using laboratory calibrations at room temperature (300 K).     

Constraints on mantle plumes on Venus: Implications for volatile history

The analysis of Venus’ gravity field and topography suggests the presence of a small number of deep mantle plumes (～9). This study predicts the number of plumes formed at the core–mantle boundary, their characteristics, and the production of partial melt from adiabatic decompression. Numerical simulations are performed using a 3D spherical code that includes large viscosity variations and internal heating. This study investigates the effect of several parameters including the core–mantle boundary temperature, the amount of internal heating, and the mantle viscosity. The smallest number of plumes is achieved when no internal heating is present. However, scaling Earth’s radiogenic heating to Venus suggests a value of ～16 TW. Cases with internal heating produce more realistic lid thickness and partial melting, but produce either too many plumes or no plumes if a high mantle temperature precludes the formation of a hot thermal boundary layer. Mantle viscosity must be reduced to at least 10²⁰ Pa s in order to include significant internal heating and still produce hot plumes. In all cases that predict melting, melting occurs throughout the upper mantle. Only cases with high core temperature (>1700 K) produce dry melting. Over time the upper mantle may have lost significant volatiles. Depending on the water content of the lower mantle, deep plumes may contribute to present-day atmospheric water via volcanic outgassing. Assuming 50 ppm water in mantle, 10 plumes with a buoyancy flux of 500 kg/s continuously erupting for 4 myr will outgas an amount of water on the order of that in the lower atmosphere. A higher level of internal heating than achieved to date, as well as relatively low mantle viscosity, may be required to achieve simulations with ～10 plumes and a thinner lid. Alternatively, if the mantle is heating up due to the stagnant lid, the effect is equivalent to having lower rates of internal heating. A temperature increase of 110 K/byr is equivalent to ?13 TW. This value along with the internal heating of 3 TW used in this study may represent the approximate heat budget of Venus’ mantle.     

Valleys,paleolakes and possible shorelines at the Libya Montes/Isidis boundary: Implications for the hydrologic evolution of Mars

We describe the results of our morphologic, stratigraphic and mineralogic investigations of fluvial landforms, paleolakes and possible shoreline morphologies at the Libya Montes/Isidis Planitia boundary. The landforms are indicative of aqueous activity and standing bodies of water, including lakes, seas and oceans, that are attributed to a complex hydrologic cycle that may have once existed on Mars in the Noachian (>3.7 Ga) and perhaps also in the Hesperian (>3.1 Ga). Our observations of the Libya Montes/Isidis Planitia boundary between 85°/86.5°E and 1.8°/5°N suggest, that (1) the termination of valley networks between roughly ?2500 and ?2800 m coincide with lake-size ponding in basins within the Libya Montes, (2) an alluvial fan and a possible delta, layered morphologies and associated Al-phyllosilicates identified within bright, polygonally fractured material at the front of the delta deposits are interpreted to be the results of fluvial activity and discharge into a paleolake, (3) the Arabia “shoreline” appears as a series of possible coastal cliffs at about ?3600 and ?3700 m indicating two distinct still stands and wave-cut action of a paleosea that temporarily filled the Isidis basin the Early Hesperian, and (4) the Deuteronilus “shoreline” appears at ?3800 m and is interpreted to be a result of the proposed sublimation residue of a frozen sea that might have filled the Isidis basin, similar to the Vastitas Borealis Formation (VBF) identified in the northern lowlands. We interpret the morphologic–geologic setting and associated mineral assemblages of the Libya Montes/Isidis Planitia boundary as results of fluvial activity, lake-size standing bodies of water and an environmental change over time toward decreasing water availability and a cold and dry climate.     

The heat flow during the formation of ribbon terrains on Venus

Ribbons are regularly spaced, between 2 and 6 km, troughs that exist on venusian tesserae, which are mainly located in, and characterize to, venusian crustal plateaus. Independent of the geological or temporal relations with other features, regularly and similarly spaced ribbons on several tesserae strongly suggest a thermal control on the thickness of the deformed layer. This can be used to constraint the heat flow at the time of ribbon formation, which holds important implications for the viability of the hypotheses that address the origin and evolution of crustal plateaus. For a brittle–ductile transition ∼1–3 km deep (as proposed from ribbon spacing), realistic strain rates, and a present-day surface temperature of 740 K, the implied heat flow is very high, 130–780 mW m⁻². If Venus has experienced higher surface temperatures due to climate forcing by massive volcanism, then the heat flow could be greatly reduced. For surface temperatures of 850–900 K the heat flow is 190–560, 60–230 and 20–130 mW m⁻² for brittle–ductile transition depths of 1, 2 and 3 km, respectively. Heat flow values around 80–100 mW m⁻² are reasonable for venusian hotspots, based on terrestrial analogs, but hardly consistent with coldspot settings. High surface temperatures are also required to maintain the crustal solidus deeper than a few kilometers during the formation of ribbon terrains. For the obtained heat flows, a solidus deeper than ∼30 km (the likely mean value for the crustal thickness) is difficult to achieve. This suggests that a substantial proportion of the crust beneath crustal plateaus was emplaced subsequently to the time when ribbon terrains were formed. Alternatively, at that time a magma reservoir inside the crust could have existed.     

Absolute magnitudes of asteroids and a revision of asteroid albedo estimates from WISE thermal observations

We obtained estimates of the Johnson V absolute magnitudes (H) and slope parameters (G) for 583 main-belt and near-Earth asteroids observed at Ond?ejov and Table Mountain Observatory from 1978 to 2011. Uncertainties of the absolute magnitudes in our sample are <0.21 mag, with a median value of 0.10 mag. We compared the H data with absolute magnitude values given in the MPCORB, Pisa AstDyS and JPL Horizons orbit catalogs. We found that while the catalog absolute magnitudes for large asteroids are relatively good on average, showing only little biases smaller than 0.1 mag, there is a systematic offset of the catalog values for smaller asteroids that becomes prominent in a range of H greater than ～10 and is particularly big above H ～ 12. The mean (H_catalog ? H) value is negative, i.e., the catalog H values are systematically too bright. This systematic negative offset of the catalog values reaches a maximum around H = 14 where the mean (H_catalog ? H) is ?0.4 to ?0.5. We found also smaller correlations of the offset of the catalog H values with taxonomic types and with lightcurve amplitude, up to ～0.1 mag or less. We discuss a few possible observational causes for the observed correlations, but the reason for the large bias of the catalog absolute magnitudes peaking around H = 14 is unknown; we suspect that the problem lies in the magnitude estimates reported by asteroid surveys. With our photometric H and G data, we revised the preliminary WISE albedo estimates made by Masiero et al. (Masired, J.R. et al. [2011]. Astrophys. J. 741, 68–89) and Mainzer et al. (Mainzer, A. et al. [2011b]. Astrophys. J. 743, 156–172) for asteroids in our sample. We found that the mean geometric albedo of Tholen/Bus/DeMeo C/G/B/F/P/D types with sizes of 25–300 km is p_V = 0.057 with the standard deviation (dispersion) of the sample of 0.013 and the mean albedo of S/A/L types with sizes 0.6–200 km is 0.197 with the standard deviation of the sample of 0.051. The standard errors of the mean albedos are 0.002 and 0.006, respectively; systematic observational or modeling errors can predominate over the quoted formal errors. There is apparent only a small, marginally significant difference of 0.031 ± 0.011 between the mean albedos of sub-samples of large and small (divided at diameter 25 km) S/A/L asteroids, with the smaller ones having a higher albedo. The difference will have to be confirmed and explained; we speculate that it may be either a real size dependence of surface properties of S type asteroids or a small size-dependent bias in the data (e.g., a bias towards higher albedos in the optically-selected sample of asteroids). A trend of the mean of the preliminary WISE albedo estimates increasing with asteroid size decreasing from D ～ 30 down to ～5 km (for S types) showed in Mainzer et al. (Mainzer, A. et al. [2011a]. Astrophys. J. 741, 90–114) appears to be mainly due to the systematic bias in the MPCORB absolute magnitudes that progressively increases with H in the corresponding range H = 10–14.     

Measurement of the radial density gradient of cosmic ray in the heliosphere by the GRAPES-3 experiment

A radial anisotropy in the flux of cosmic rays in heliosphere was theoretically predicted by Parker and others within the framework of the diffusion–convection mechanism. The solar wind is responsible for sweeping out the galactic cosmic rays, creating a radial density gradient within the heliosphere. This gradient coupled with the interplanetary magnetic field induces a flow of charged particles perpendicular to the ecliptic plane which was measured and correctly explained by Swinson, and is hereafter referred as ‘Swinson flow’. The large area GRAPES-3 tracking muon telescope offers a powerful probe to measure the Swinson flow and the underlying radial density gradient of the galactic cosmic rays at a relatively high rigidity of ∼100 GV. The GRAPES-3 data collected over a period of six years (2000–2005) were analyzed and the amplitude of the Swinson flow was estimated to be (0.0644 ± 0.0008)% of cosmic ray flux which was an ∼80σ effect. The phase of the maximum flow was at a sidereal time of (17.70 ± 0.05) h which was 18 min earlier than the expected value of 18 h. This small 18 min phase difference had a significance of ∼6σ indicating the inherent precision of the GRAPES-3 measurement. The radial density gradient of the galactic cosmic rays at a median rigidity of 77 GV was found to be 0.65% AU⁻¹.     

Rotation periods of binary asteroids with large separations – Confronting the Escaping Ejecta Binaries model with observations

Durda et al. (Durda, D.D., Bottke, W.F., Enke, B.L., Merline, W.J., Asphaug, E., Richardson, D.C., Leinhardt, Z.M. [2004]. Icarus 170, 243–257), using numerical models, suggested that binary asteroids with large separation, called Escaping Ejecta Binaries (EEBs), can be created by fragments ejected from a disruptive impact event. It is thought that six binary asteroids recently discovered might be EEBs because of the high separation between their components (～100 > a/R_p > ～20).However, the rotation periods of four out of the six objects measured by our group and others and presented here show that these suspected EEBs have fast rotation rates of 2.5–4 h. Because of the small size of the components of these binary asteroids, linked with this fast spinning, we conclude that the rotational-fission mechanism, which is a result of the thermal YORP effect, is the most likely formation scenario. Moreover, scaling the YORP effect for these objects shows that its timescale is shorter than the estimated ages of the three relevant Hirayama families hosting these binary asteroids. Therefore, only the largest (D ～ 19 km) suspected asteroid, (317) Roxane, could be, in fact, the only known EEB.In addition, our results confirm the triple nature of (3749) Balam by measuring mutual events on its lightcurve that match the orbital period of a nearby satellite in addition to its distant companion. Measurements of (1509) Esclangona at different apparitions show a unique shape of the lightcurve that might be explained by color variations.     

Cluster observations of Pc 1–2 waves and associated ion distributions during the October and November 2003 magnetic storms

Unusual wave activity in the Pc 1–2 frequency band (0.1–5 Hz) was observed by the Cluster spacecraft in association with the two large geomagnetic storms of late 2003. During the onset of the Hallowe’en storm on October 29, 2003, intense broadband activity between ∼0.1 and 0.6 Hz appeared at all 4 spacecraft on both sides of the magnetic equator at perigee (near 1400 UT and 08:45 MLT). Power was especially strong and more structured in frequency in the compressional component: a minimum in wave power was observed at 0.38 Hz, corresponding to the oxygen ion cyclotron frequency. Poynting vector calculations indicated that wave power was primarily directed radially inward rather than along the magnetic field. Narrowband purely compressional waves near 0.15 Hz appeared at higher dayside latitudes in the southern hemisphere. CIS ion spectrometer data during this pass revealed that O⁺ was the dominant energetic ion. During the recovery phase of the November storm, on November 22, 2003, predominantly transverse 1.8 Hz waves with peak-to-peak amplitude of 10 nT were observed by all four spacecraft near perigee at L=4.4. During this more typical Pc 1 event, wave power was directed along B, toward the northern ionosphere. An unusually polarized 2.3 Hz emission (with power in the radial and compressional, but not azimuthal directions) was observed at L=5.4–5.9, 10–15° south of the magnetic equator. We infer that this wave event may have been generated on lower L shells and propagated obliquely to Cluster's location. Consistent with other recent observations, anisotropic plasma sheet/ring current proton distributions appeared to be a necessary condition for occurrence of waves during both passes, but was not always a sufficient condition. The transverse waves of November 22 occurred in regions which also contained greatly increased fluxes of cool ions (E<1 keV). On both days, Cluster observed features not previously reported, and we note that the purely compressional nature of the October 29 events was not anticipated in previous theoretical studies. The fact that these unusually polarized waves occurred in association with very intense geomagnetic storms suggests that they are likely to be extremely rare.