LONG-TERM RIGOROUS NUMERICAL INTEGRATION OF NAVIER-STOKES EQUATION BY NEWTON-GMRES ITERATION

The recent result of an orbit continuation algorithm has provided a rigorous method for long-term numer- ical integration of an orbit on the unstable manifold of a periodic solution. This algorithm is matrix-free and em- ploys a combination of the Newton-Raphson method and the Krylov subspace method. Moreover, the algorithm adopts a multiple shooting method to address the problem of orbital instability due to long-term numerical integra- tion. The algorithm is described through computing the extension of unstable manifold of a recomputed Nagata~s lower-branch steady solution of plane Couette flow, which is an example of an exact coherent state that has recently been studied in subcritical transition to turbulence.     

Canadian advanced life support capacities and future directions

Canada began research on space-relevant biological life support systems in the early 1990s. Since that time Canadian capabilities have grown tremendously, placing Canada among the emerging leaders in biological life support systems. The rapid growth of Canadian expertise has been the result of several factors including a large and technically sophisticated greenhouse sector which successfully operates under challenging climatic conditions, well planned technology transfer strategies between the academic and industrial sectors, and a strong emphasis on international research collaborations. Recent activities such as Canada’s contribution of the Higher Plant Compartment of the European Space Agency’s MELiSSA Pilot Plant and the remote operation of the Arthur Clarke Mars Greenhouse in the Canadian High Arctic continue to demonstrate Canadian capabilities with direct applicability to advanced life support systems. There is also a significant latent potential within Canadian institutions and organizations with respect to directly applicable advanced life support technologies. These directly applicable research interests include such areas as horticultural management strategies (for candidate crops), growth media, food processing, water management, atmosphere management, energy management, waste management, imaging, environment sensors, thermal control, lighting systems, robotics, command and data handling, communications systems, structures, in-situ resource utilization, space analogues and mission operations. With this background and in collaboration with the Canadian aerospace industry sector, a roadmap for future life support contributions is presented here. This roadmap targets an objective of at least 50% food closure by 2050 (providing greater closure in oxygen, water recycling and carbon dioxide uptake). The Canadian advanced life support community has chosen to focus on lunar surface infrastructure and not low Earth orbit or transit systems (i.e. microgravity applications). To advance the technical readiness for the proposed lunar missions, including a lunar plant growth lander, lunar “salad machine” (i.e. small scale plant production unit) and a full scale lunar plant production system, a suite of terrestrial developments and analogue systems are proposed. As has been successfully demonstrated by past Canadian advanced life support activities, terrestrial technology transfer and the development of highly qualified personnel will serve as key outputs for Canadian advanced life support system research programs. This approach is designed to serve the Canadian greenhouse industry by developing compliance measures for mitigating environmental impact, reducing labour and energy costs as well as improving Canadian food security, safety and benefit northern/remote communities.     

A gravity wave analysis near to the Andes Range from GPS radio occultation data and mesoscale numerical simulations: Two case studies

Global maps of potential wave energy per unit mass, recently performed with the Global Positioning System (GPS) Radio Occultation (RO) technique and different satellite missions (CHAMP and SAC-C since 2001, GRACE and COSMIC since 2006) revealed in Argentina, at the eastern side of the highest Andes Mountains, a considerable wave activity (WA) in comparison with other extra-tropical regions. The main gravity wave (GW) sources in this natural laboratory are deep convection (mainly during late Spring and Summer), topographic forcing and geostrophic adjustment.     

Linkage of the ionospheric peak electron density and height deduced from the topside sounding data

Topside sounding electron density profiles are analyzed to explore interrelations of the F2 layer critical frequency and the peak height for a representative set of conditions provided by ISIS1, ISIS2, IK19 and Cosmos-1809 satellites for the period of 1969–1987. The foF2 and hmF2 are delivered with exponential extrapolation of electron density profile to zero of its 1st derivative. It is shown that the linear regression exists between foF2 and hmF2 under different conditions. The linkage between the two parameters amended to the empirical model of the peak height [Gulyaeva, T.L., Bradley, P.A., Stanislawska, I., Juchnikowski, G. Towards a new reference model of hmF2 for IRI. Adv. Space Res. 42, 666–672, doi:10.1016/j.asr.2008.02.021, 2008] results in an empirical model of the both foF2 and hmF2 expressed by superposition of functions in terms of local-time, season, geodetic longitude, modified dip latitude and solar activity. For the solar activity we use a proxy Fsp index averaged from the mean solar radio flux F10.7s for the past 81 days (3 solar rotations) and F10.7 value for 1 day prior the day of observation. Impact of geomagnetic activity is not discernible with the topside sounding data due to mixed positive and negative storm-time effects. Appreciable differences have been revealed between IRI-CCIR predictions and outcome of the new model which might be attributed to the different techniques of the peak electron density and height derivation, different epochs and different global distribution of the source data as well as the different mathematical functions involved in the maps and the model presentation.     

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.     

Hybrid simulations of the O ion escape from Venus: Influence of the solar wind density and the IMF x component

As an initial effort to study the evolution of the Venus atmosphere, the influence of the solar wind density and the interplanetary magnetic field (IMF) x component (the x-axis points from Venus towards the Sun) on the O⁺ ion escape rate from Venus is investigated using a three-dimensional quasi-neutral hybrid (HYB-Venus) model. The HYB-Venus model is first applied to a case of the high-density (100 cm⁻³) solar wind interaction with Venus selected from the Pioneer Venus Orbiter observations to demonstrate its capability for the study. Two sets of simulations with a wide range of solar wind densities and different IMF x components are then performed. It is found that the O⁺ ion escape rate increases with increasing solar wind density. The O⁺ ion escape rate saturates when the solar wind density becomes high (above 100 cm⁻³). The results also suggest that the IMF x component enhances the O⁺ ion escape rate, given a fixed IMF component perpendicular to the x-axis. Finally, the results imply a higher ion loss rate for early-Venus, when solar conditions were dramatically different.     

Orbit determination for the GOCE satellite

Precise Orbit Determination (POD) for the Gravity field and steady-state Ocean Circulation Explorer (GOCE), the first core explorer mission by the European Space Agency (ESA), forms an integrated part of the so-called High-Level Processing Facility (HPF). Two POD chains have been set up referred to as quick-look Rapid and Precise Science Orbit determination or RSO and PSO, respectively. These chains make use of different software systems and have latencies of 1 day and 2 weeks, respectively, after tracking data availability. The RSO and PSO solutions have to meet a 3-dimensional (3D) position precision requirement of 50 cm and a few cm, respectively. The tracking data will be collected by the new Lagrange GPS receiver and the predicted characteristics of this receiver have been taken into account during the implementation phase of the two chains.     

Surface temperature estimation in Singhbhum Shear Zone of India using Landsat-7 ETM+ thermal infrared data

Land surface temperature (LST) is an important factor in global change studies, heat balance and as control for climate change. A comparative study of LST over parts of the Singhbhum Shear Zone in India was undertaken using various emissivity and temperature retrieval algorithms applied on visible and near infrared (VNIR), and thermal infrared (TIR) bands of high resolution Landsat-7 ETM+ imagery. LST results obtained from satellite data of October 26, 2001 and November 2, 2001 through various algorithms were validated with ground measurements collected during satellite overpass. In addition, LST products of MODIS and ASTER were compared with Landsat-7 ETM+ and ground truth data to explore the possibility of using multi-sensor approach in LST monitoring. An image-based dark object subtraction (DOS3) algorithm, which is yet to be tested for LST retrieval, was applied on VNIR bands to obtain atmospheric corrected surface reflectance images. Normalized difference vegetation index (NDVI) was estimated from VNIR reflectance image. Various surface emissivity retrieval algorithms based on NDVI and vegetation proportion were applied to ascertain emissivities of the various land cover categories in the study area in the spectral range of 10.4–12.5 μm. A minimum emissivity value of about 0.95 was observed over the reflective rock body with a maximum of about 0.99 over dense forest. A strong correlation was established between Landsat ETM+ reflectance band 3 and emissivity. Single channel based algorithms were adopted for surface radiance and brightness temperature. Finally, emissivity correction was applied on ‘brightness temperature’ to obtain LST. Estimated LST values obtained from various algorithms were compared with field ground measurements for different land cover categories. LST values obtained after using Valor’s emissivity and single channel equations were best correlated with ground truth temperature. Minimum LST is observed over dense forest as about 26 °C and maximum LST is observed over rock body of about 38 °C. The estimated LST showed that rock bodies, bare soils and built-up areas exhibit higher surface temperatures, while water bodies, agricultural croplands and dense vegetations have lower surface temperatures during the daytime. The accuracy of the estimated LST was within ±2 °C. LST comparison of ASTER and MODIS with Landsat has a maximum difference of 2 °C. Strong correlation was found between LST and spectral radiance of band 6 of Landsat-7 ETM+. Result corroborates the fact that surface temperatures over land use/land cover types are greatly influenced by the amount of vegetation present.     

22 Years of AJISAI spin period determination from standard SLR and kHz SLR data

Satellite Laser Ranging (SLR) is a powerful and efficient technique to measure spin parameters of satellites equipped with corner cube reflectors. We obtained spin period determination of the satellite AJISAI from SLR data only: 17246 pass-by-pass estimates from standard 1–15 Hz SLR data (14/Aug/1986–30/Dec/2008) and 1444 pass-by-pass estimates (9/Oct/2003–30/Dec/2008) from data of the first 2 kHz SLR system from Graz, Austria. A continuous history of the slowing down of AJISAI spin is derived from frequency analysis, and corrected for the apparent effects. The apparent corrections, elaborated here, allowed very accurate determination of AJISAI initial spin period: 1.4855 ± 0.0007 [s]. The paper identifies also non-gravitational effects as a source of the periodical changes in the rate of slowing down of the satellite.     

Radial propagation velocity of energetic particle injections according to measurements onboard the <Emphasis Type="Italic">Cluster</Emphasis> satellites

Injections of energetic electrons with a dispersion over energies were observed during the February 23, 2004 (at about 03:20 UT) substorm onboard the Cluster satellites in the vicinity of perigee near the midnight meridian. The delays in the particle observation caused by the energy dependence of the magnetic drift velocities made it possible to determine the position and time of the beginning of the drift, tracing the trajectories of the leading center of particles back in time in the magnetospheric model. The comparisons of the measurements of four satellites allowed us to determine the radial propagation of the injection front with a velocity of 100–150 km/s at a distance of 7–9 R _E. The comparison with a few previous measurements shows a substantial slowing down of injections as they approached the Earth, and this confirms the prospects of this method for more detailed study of propagation of plasma injection into the inner magnetosphere.