Laboratory calibration measurements of a piezoelectric lead zirconate titanate cosmic dust detector at low velocities

A cosmic dust monitor for use onboard a spacecraft is currently being developed using a piezoelectric lead zirconate titanate element (PZT). Its characteristics of the PZT sensor is studied by ground-based laboratory impact experiments using hypervelocity particles supplied by a Van de Graaff accelerator. The output signals obtained from the sensor just after the impact appeared to have a waveform that was explicitly related to the particle’s impact velocity. For velocities less than ∼6 km/s, the signal showed an oscillation pattern and the amplitude was proportional to the momentum of the impacting particle. For higher velocities, the signal gradually changed to a single waveform. The rise time of this single waveform was proportional to the particle’s velocity for velocities above ∼6 km/s. The present paper reports on results for the low velocity case and especially discusses the effect of an outer coating of the sensor with a paint, which is used to reduce heating by solar radiation.     

Velocity dependent response of a piezoelectric element to hypervelocity microparticles

Piezoelectric lead–zirconate–titanate (PZT) elements were studied by directly bombarding hypervelocity silver and iron particles in the velocity range of 2–63 km/s. This report focuses on the concept of “immediately after collision”. This concept was realized in the leading part of a pulsed signal that exhibited no effects due to reflection. The signal form was discussed from the viewpoint that information on impact was recorded by the waveform. The results indicated that the output amplitude was independent of the thickness of a PZT element, and that the waveform was explicitly related to the velocity at collision. Further, it appeared as an oscillating form at velocities less than 6 km/s. As the velocity increased, it gradually changed to a solitary pulse. The rise time of the single solitary waveform was linearly related to velocities above 6 km/s. The potential of a single PZT element is discussed as a real-time detector for hypervelocity microparticles.     

Response of a pentagonal PZT element as a component of a 4π-real-time detector

The detector characteristics of a pentagonal element were studied by colliding it with hypervelocity micro-particles. A charge-sensitive amplifier was developed for the element of its capacitance ∼10 nF. The output amplitudes were expressed as a linear function of the momentum at collision. Empirical formulas obtained from on-ground experiments could be used for the calibration of the detector. The pentagonal element was potential to measure the momentum during collision from the output amplitude. A set of electrodes on the surface was used to confirm the measurement of the coordinates at collision. A possible application of this pentagonal element on a real-time dust detector was discussed.     

Characteristics of piezoelectric lead zirconate titanate multilayered detector bombarded with hypervelocity iron particles

A cosmic dust detector is currently being developed using a piezoelectric lead zirconate titanate (PZT) element. The characteristics of the multilayered detector (MD), which was composed of one hundred PZT disks, were investigated by bombarding it with hypervelocity iron particles supplied by a Van de Graaff accelerator. It was confirmed that there was a linear relationship between the signal amplitude observed from MD and the momentum of the particles. As compared with the single-layered detector (SD) that was composed of one PZT disk, it was found that the sensitivity of MD was ∼3 times higher than that of SD within the limits of the experimental conditions.     

Systematic examination of the geomagnetic storm sudden commencement using multi resolution analysis

Geomagnetic storm sudden commencements (SSC) contain a wealth of information which is useful in many applications. It is important to point out that the SSCs used in this study are sampled at the rate of one sample per second in order get use of such high resolution data. In this paper, two studies are made on the geomagnetic SSC and an SSC onset automatic detection algorithm is introduced. The first study is about finding the relationship between the SSC rise time and its amplitude. Where it is found that there is a positive correlation between the amplitude and the amplitude gradient which is the amplitude divided by rise time. The second study is checking the spectrum of the SSC, starting from its onset until the end of the SSC rise time. This check had proved that the SSC contains both low and high frequency regions. This led us to introduce a new term, namely the SSC variation rate (VR). This VR is defined as the maximum rate of change of the field in the higher-frequency region of the SSC. These two studies were the guide to build an SSC automatic detector of one sample per second data using multi resolution analysis (MRA) of the discrete wavelet transform (DWT). The data set contains 134 SSCs with different VRs that were collected from the Circum-pan Pacific Magnetometer Network (CPMN). It is found that the standard deviation of the detection error is 41 s and that the average error is 9 s. From the calculated error distribution function, it is found that the detection error is within the range of −1.5 to 3 min. The detection process, as will be shown in the article, takes 70 s for one station and 3 min if the decision is related to the detection(s) of other stations. These results demonstrate the superiority of the proposed algorithm over other algorithms in which the detection error ranges between −8 and 5 min and the detection process takes 2–10 min. In addition to being faster and more accurate than the other algorithms, the proposed algorithm is the first algorithm that automatically detects the SSC onset times from high-resolution data unlike previous studies that focused on determining the SSC times automatically using one-minute resolution data.     

The gravitropic and phototropic responses of wheat grown in a space greenhouse prototype with hemispherical planting surface

The time course of gravicurvature of 3-day-old wheat (Triticum aestivum L., cv. Apogee) coleoptiles and 7-day-old wheat stems were studied in darkness and under red and red-blue light illumination after declination from the vertical at various angles. The experiments showed that the shortest gravitropic curvature corresponded to 30° initial angle of gravistimulation (IAG). The time course became longer as the IAG increased and with plant age. The effects of unilateral red (660 nm) and red-blue light (660 nm; 470 nm) at photosynthetic photon flux (PPF) of 30 μmol m⁻² s⁻¹ on the curvature of 3-day-old coleoptiles were evaluated. Red light did not produce phototropic bending of wheat coleoptiles in contrast with red-blue light. The analysis of experimental data showed that the curvature in response to a gravitropic stimulus or to combined gravity-light stimuli were not statistically different. Time course of gravitropic curvature were used to determine the acceptable crop rotation rate around the horizontal axis. Approximation of stem bending to a linear dynamic system described by a first-order aperiodic element with a lag allowed the determination of the dependence of the amplitude of apex oscillations on the rate of horizontal rotation under 1-g conditions. The calculated lowest minimal rotation rate (MRR) minimizing the gravitropic effects on wheat was about 1 revolution per hour (rph). Rotating the plant growth chamber (PGC) at a rate of more than MRR eliminated the effect of gravitropic curvature.     

溶胶-水热法合成PZT纳米粉体及性能研究

结合溶胶-凝胶和水热法的优点,采用改进的溶胶-水热复合法在较低温下制备出纯相PZT纳米粉体,并对粉体的烧结性能进行了研究,分别讨论了烧结温度、保温时间等工艺参数对烧结陶瓷密度、微观结构和压电性能的影响.270℃保温热处理2h可合成出粒径约为20~30nm,具有良好分散性的钙钛矿型PZT纳米粉体,且具有良好的烧结活性.在1150℃烧结保温2h,压电性能达到最优(机电耦合系数: 0.50,机械品质因数: 410,压电常数: 220pC/N,介电常数: 1060).结果表明,溶胶水热复合法具有合成温度低、组分易于控制、粉体烧结活性高的优点.      

谐振式液体密度传感器压电激励与检测的实现

本文利用正压电效应和逆压电效应,实现谐振式液体密度传感器的压电激励与压电检测,原理简单、功耗低、检测方便。设计了基于该激励和检测方式的自动增益控制(Automatic Gain Control,AGC)电路,当传感器的输出信号发生较大变化时,激励信号的增益可随输出信号的幅值进行自动调节,使输出信号的幅值保持稳定。该电路具有结构简单,适应性好,可靠性高,调节速度快等优点。同时,对所设计的谐振式液体密度传感器进行了标定实验,该密度传感器的精度约为±1.0kg/m³,重复性约为±0.05%,可以实现液体密度的高精度实时在线测量。     

Properties of traveling atmospheric disturbances (TADs) inferred from CHAMP accelerometer observations

Densities derived from accelerometer measurements on the CHAMP satellite near 400 km are used to statistically establish characteristics of large-scale (>1000 km) traveling atmospheric disturbances (TADs). Only TADs that at least propagate from the auroral zone to the equator are analyzed here, and a total of 21 identifiable events are found over the years 2001–2007. The average speed of all TADs, regardless of local time, is 646 ± 122 ms⁻¹. The average speeds on the dayside and nightside are 595 ± 127 ms⁻¹ and 685 ± 106 ms⁻¹, respectively, i.e., the speed appears to be 10% higher on average on the nightside. On six occasions TADs were only detected on the night side; however, TADs on the dayside often appear more distinctly in the data. Moreover, contrary to some theoretical expectations, dayside TADs do not dissipate more readily than night side TADs, although much less are detected between 8–20 solar local time. No clear dependence of TAD amplitude or phase speed with respect to Kp, or rate of increase of Kp, is found.     

Geomagnetic pulsation over conjugate locations during geomagnetic storms and substorm

Geomagnetic data collected during magnetic storm over magnetically conjugate pair (according to IGRF 2000) of high latitude stations viz., Maitri (70° 45′ S, 11° 42′ E) and Tromso (69° 40′ N, 18° 56′ E) reveal that amplitudes of Pc6 pulsation characteristically differ. The amplitude obtained from horizontal magnetic field for the Pc6 pulsation frequency range between 0.6 and 1.6 mHz significantly differs in time corresponding to peak amplitude. The relative differences in its time of occurrence found to gradually increase around initial phase of storm and remain exactly out of phase at peak amplitude of storm. Thence, it is found to be in agreement in phase gradually until storm unwinds. This indicates that simultaneous amplitude of Pc6 pulsation at conjugate pair of stations and its time of occurrence could be a key factor to infer storm arrivals somewhat prior to its peak effects. The emphasize remains on prediction of storm arrivals only by utilizing ground based magnetometer observations. However, it is necessary to understand differences on the basis of weak, moderate, strong, and super strong cases and more exactly how they behave along the line of magnetic Meridian. Nevertheless, the analysis implies that geo-effective magnetic ejecta/clouds/CIRs/sheaths/CMEs/ICMEs giving rise to geomagnetic storm can be predicted ahead of its peak effects by having magnetometer data over conjugate locations.     

Latitudinal variation of F2-region response to geomagnetic disturbance

The effect of geomagnetic storms on the F2 region was studied by calculating the deviation, ΔfoF2, of foF2 during 40 magnetic storms, ranging from moderate (Dst < −50 nT) to very intense (Dst < −200 nT) of the 21st solar cycle. In order to study the variation of storm-time foF2 with latitude, season and storm strength, ionosonde data were obtained from eight stations spanning a latitudinal range of +60–−60°. The stations chosen lay in a narrow longitudinal range of 140–151°, so that local time difference between the stations is practically negligible. The features exhibited by positive and negative phases were essentially different. The storm time ΔfoF2 clearly exhibited a latitudinal variation and this variation were found to be coupled with the seasonal variation. As for the variation with storm intensity, though ΔfoF2 was found to vary even between two storms of almost equal intensity, the amplitude of a positive or negative phase, |ΔfoF2max| showed a distinct upper limit for each intensity category of storms.     

Ground-based and satellite DC-ULF electric field anomalies around Wenchuan M8.0 earthquake

The ground-based and satellite DC-ULF electric field data were analyzed around Wenchuan M8.0 earthquake on May 12, 2008 in China. The results show that ground electric field anomalies occurred at 3 stations located to the north and south of the epicenter with the amplitude of 3–100 mV/km. The change shapes and their amplitude of ground electric field anomalies are different largely due to their individual underground layer conductivity, water level and so on. The analysis of long time series illustrates that the abnormal geoelectric field started since March 2008. Onboard the DEMETER satellite, the ULF waveforms of electric field were collected and processed by wavelet transform method. The disturbances in the ionosphere were about 3–5 mV/m at a frequency band lower than 0.5 Hz. When the ground and space electric field anomalies were compared, their occurrence time and spatial distribution points are consistent with each other, including the long time anomalies from March 2008 and the short term ones 1–2 days before the Wenchuan earthquake. Finally, the coupling mechanism was discussed.     

New flare model using recent measurements of the solar ultraviolet irradiance

The solar photon output from the Sun, which was once thought to be constant, varies considerably over time scales from seconds during solar flares to years due to the solar cycle. This is especially true in the wavelengths shorter than 190 nm. These variations cause significant deviations in the Earth and space environment on similar time scales, which then affects many things including satellite drag, radio communications, atmospheric densities and composition of particular atoms, molecules, and ions of Earth and other planets, as well as the accuracy in the Global Positioning System (GPS). The Flare Irradiance Spectral Model (FISM) is an empirical model that estimates the solar irradiance at wavelengths from 0.1 to 190 nm at 1 nm resolution with a time cadence of 60 s. This is a high enough temporal resolution to model variations due to solar flares, for which few accurate measurements at these wavelengths exist. This model also captures variations on the longer time scales of solar rotation (days) and solar cycle (years). Daily average proxies used are the 0–4 nm irradiance, the Mg II c/w, F10.7, as well as the 1 nm bins centered at 30.5 nm, 121.5 (Lyman Alpha), and 36.5 nm. The GOES 0.1–0.8 nm irradiance is used as the flare proxy. The FISM algorithms are given, and results and comparisons are shown that demonstrate the FISM estimations agree within the stated uncertainties to the various measurements of the solar Vacuum Ultraviolet (VUV) irradiance.     

Numerical simulation of equatorial plasma bubbles over Cachimbo: COPEX campaign

The problem of day-to-day variability in onset of equatorial spread F (ESF) is addressed using data from the 2002 COPEX observational campaign in Brazil and numerical modeling. The observational results show that for values of virtual height of the F layer base less than 355 km at around 18:35 LT, and for the prereversal peak enhancement of the vertical plasma drift (Vp) less than 30 m/s, the spread-F (ESF) was absent on four nights over Cachimbo (9.5°S, 54.8°W, dip latitude = −2.1°). In this work we analyze the geophysical conditions for the generation of the irregularities by comparing the nights with and without the ESF. In the comparison a numerical code is used to simulate plasma irregularity development in an extended altitude range from the bottom of the equatorial F   layer. The code uses the flux corrected transport method with Boris–Book’s flux limiter for the spatial integration and a predictor–corrector method for the direct time integration of the continuity equation for O⁺${\mathrm{O}}^{+}$ and the SOR (Successive-Over-Relaxation) method for electric potential equation. The code is tested with different evening eastward electric fields (or vertical drifts Vp < 30 m/s and Vp > 30 m/s) in order to study the influence of the prereversal enhancement in the zonal electric field on plasma bubble formation and development. The code also takes into account the zonal wind, the vertical electric field and the collision frequency of ions with neutrals and the amplitude of initial perturbation. The simulation shows a good agreement with the observational results of the ESF. The results of the code suggest that the instability can grow at the F layer bottomside by the Rayleigh–Taylor mechanism only when the Vp > 30 m/s. In the analyzed cases we have considered the competition of other geophysical parameters in the generation of plasma structures.     

Monitoring the South Atlantic Anomaly using ATSR instrument series

The South Atlantic Anomaly (SAA) has been monitored for 19 years using the Along Track Scanning Radiometer (ATSR) series of instruments onboard the ERS-1, ERS-2 and ENVISAT ESA satellites. The time evolution of the night-time particle induced noise in the short wavelength infrared (SWIR, 1.6 μm) and visible (VIS, 0.55 μm) channels of the ATSR instrument series have been analysed. The monthly location and extension of the SAA are inferred by fitting a two-dimensional, elliptical Gaussian function to the coordinates of the night-time hot spots detected over the SAA region. The location of the centre of the SAA is found to drift westwards with an average drift rate of about 0.24 deg/year and northward with an average drift rate of about 0.12 deg/year. Irregularities are found where the drift speed is inverted and the SAA moves eastward and southward. Results indicate that, as expected, the retrieved values of SAA’s strength and extension are anti-correlated with the solar activity expressed by the solar flux at 10.7 cm (F10.7). Finally, the peak-to-peak amplitude of the seasonal variation of the SAA strength, estimated from monthly VIS data, is found to be 30% of the average value with the annual to semiannual amplitude ratio of 1.38.     

Global structure,seasonal and interannual variability of the eastward propagating tides seen in the SABER/TIMED temperatures (2002–2007)

The present paper is focused on the global spatial (altitude and latitude) structure, seasonal and interannual variability of the most stable in amplitude and phase eastward propagating diurnal and semidiurnal tides with zonal wavenumbers 2 and 3 derived from the SABER/TIMED temperatures for full 6 years (January 2002–December 2007). The tidal results are obtained by an analysis method where the tides (migrating and nonmigrating) and the planetary waves (zonally travelling, zonally symmetric and stationary) are simultaneously extracted from the satellite data. It has been found that the structures of the eastward propagating diurnal tides with zonal wavenumbers 3 and 2 change from antisymmetric with respect to the equator below ∼85 km height, to more symmetric above ∼95 km. The seasonal behavior of the DE3 is dominated by annual variation with maximum in August–September reaching average (2002–2007) amplitude of ∼15 K, while that of the DE2 by semiannual variation with solstice maxima and with average amplitude of ∼8 K. These tides revealed some interannual variability with a period of quasi-2 years. The seasonal behavior of the eastward propagating semidiurnal tide with zonal wavenumber 2 in the southern hemisphere (SH) is dominated by annual variation with maximum in the austral summer (November–January) while that in the northern hemisphere (NH) by semiannual variation with equinoctial maxima. The SE2 maximizes near 115 km height and at latitude of ∼30° reaching an average amplitude of ∼6 K. The seasonal behavior of the eastward propagating semidiurnal tide with zonal wavenumber 3 in both hemispheres indicates a main maximum during June solstice and a secondary one during December solstice. The tide maximizes near 110–115 km height and at a latitude of ∼30° reaching an average amplitude of ∼4.8 K in the SH and ∼4 K in the NH. The tidal structures of the two eastward propagating semidiurnal tides are predominantly antisymmetric about the equator.     

On the spatial extent of the proton radiation belt from solar cell output variation of the Akebono satellite

Output current of silicon solar cells of Akebono satellite orbiting in the inner magnetosphere decreased from 13 A in 1989 to about 7 A in 2009, due to accumulated damage by energetic particles. A fair correlation between the monthly decrease rate of solar cell output current and trapped proton flux was found between 1989 and 1992, indicating that trapped energetic protons are responsible for the solar cell degradation. Simple variation from month to month was dominated by a fluctuation with one-year period and the correlation was not discernible after 1993. On the other hand, during orbiting the earth, another kind of decrease of output current emerged at low altitude above the sub-solar point, i.e., in the midst of sunlit condition. The decrease was larger in later years. The fact indicates more prominent temperature effect in the later years because of progress of the degradation. By removing data affected by the orbit condition above, variation from the same month in the previous year shows a clear oscillation due to orbit precession correlated with trapped energetic proton flux up to 1996. The amplitude of the annual variation oscillation tends to be larger than that expected from a degradation model based on energetic proton distribution of the NASA’s AP8 model. The larger oscillation amplitude suggests that the proton radiation belt was more sharply localized than given by the AP8 model throughout the early half of 1990s.     

Use of two-line element data for thermosphere neutral density model calibration

Traditional empirical thermospheric density models are widely used in orbit determination and prediction of low-Earth satellites. Unfortunately, these models often exhibit large density errors of up to around 30% RMS. Density errors translate into orbit errors, adversely affecting applications such as re-entry operations, manoeuvre planning, collision avoidance and precise orbit determination for geodetic missions. The extensive database of two-line element (TLE) orbit data contains a wealth of information on satellite drag, at a sufficiently high spatial and temporal resolution to allow a calibration of existing neutral density models with a latency of one to two days. In our calibration software, new TLE data for selected objects is converted to satellite drag data on a daily basis. The resulting drag data is then used in a daily adjustment of density model calibration parameters, which modify the output of an existing empirical density model with the aim of increasing its accuracy. Two different calibration schemes have been tested using TLE data for about 50 objects during the year 2000. The schemes involve either height-dependent scale factors to the density or corrections to CIRA-72 model temperatures, which affect the density output based on a physical model. Both schemes have been applied with different spherical harmonic expansions of the parameters in latitude and local solar time. Five TLE objects, varying in perigee altitude between 280 and 530 km, were deliberately not used during calibration, in order to provide independent validation. Even with a single daily parameter, the RMS density model error along their tracks can already be reduced from the 30% to the 15% level. Adding additional parameters results in RMS errors lower than 12%.     

Radiometric estimation of water vapor content over Brazil

A multi-channel microwave radiometre (make: Radiometrics Corporation) is installed at Instituto Nacional de Pesquisas Espaciais–INPE, Brazil (22°S). The radiometric output of two channels of the radiometer in the form of brightness temperature at 23.834 GHz and 30 GHz, initially, were used to find out the ambient water vapor content and the non-precipitable cloud liquid water content. The necessary algorithm was developed for the purpose. The best results were obtained using the hinge frequency 23.834 GHz and 30 GHz pair having an r.m.s. error of only 2.64. The same methodology was then adopted exploiting 23.034 GHz and 30 GHz pair. In that case the r.m.s. error was 3.42. These results were then compared with those obtained over Kolkata (22°N), India, by using 22.234 GHz and 31.4 GHz radiometric data. This work conclusively suggests the use of a frequency should not be at the water vapor resonance line. Instead, while measuring the vapor content for separation of vapor and cloud liquid, one of them should be a few GHz left or right from the resonance line i.e., at 23.834 GHz and the other one should be around 30 GHz.     

Improved orbit predictions using two-line elements

The density of orbital space debris constitutes an increasing environmental challenge. There are two ways to alleviate the problem: debris mitigation and debris removal. This paper addresses collision avoidance, a key aspect of debris mitigation. We describe a method that contributes to achieving a requisite increase in orbit prediction accuracy for objects in the publicly available two-line element (TLE) catalog. Batch least-squares differential correction is applied to the TLEs. Using a high-precision numerical propagator, we fit an orbit to state vectors derived from successive TLEs. We then propagate the fitted orbit further forward in time. These predictions are validated against precision ephemeris data derived from the international laser ranging service (ILRS) for several satellites, including objects in the congested sun-synchronous orbital region. The method leads to a predicted range error that increases at a typical rate of 100 m per day, approximately a 10-fold improvement over individual TLE’s propagated with their associated analytic propagator (SGP4). Corresponding improvements for debris trajectories could potentially provide conjunction analysis sufficiently accurate for an operationally viable collision avoidance system based on TLEs only.