A multi-objective approach to the design of low thrust space trajectories using optimal control

In this article, we introduce a novel three-step approach for solving optimal control problems in space mission design. We demonstrate its potential by the example task of sending a group of spacecraft to a specific Earth L ₂ halo orbit. In each of the three steps we make use of recently developed optimization methods and the result of one step serves as input data for the subsequent one. Firstly, we perform a global and multi-objective optimization on a restricted class of control functions. The solutions of this problem are (Pareto-)optimal with respect to ΔV and flight time. Based on the solution set, a compromise trajectory can be chosen suited to the mission goals. In the second step, this selected trajectory serves as initial guess for a direct local optimization. We construct a trajectory using a more flexible control law and, hence, the obtained solutions are improved with respect to control effort. Finally, we consider the improved result as a reference trajectory for a formation flight task and compute trajectories for several spacecraft such that these arrive at the halo orbit in a prescribed relative configuration. The strong points of our three-step approach are that the challenging design of good initial guesses is handled numerically by the global optimization tool and afterwards, the last two steps only have to be performed for one reference trajectory.     

Quantification of the surface brackish water layer and frontal zones in southern Chilean fjords between Boca del Guafo (43°30′S) and Estero Elefantes (46°30′S)

The southern Chilean region between the Boca del Guafo passage and Estero Elefantes contains one of the estuarine zones with the greatest freshwater influence on the planet. At the surface, plumes of freshwater from the fjord heads to their mouths, emptying into the Moraleda–Costa–Elefantes channel system and then the coastal ocean. The influence of this freshwater on the region’s estuarine dynamics, coastal ecology, and biogeochemical processes has only recently begun to be elucidated.Using hydrographic data from the CIMAR-Fiordos cruises (1998–2001), this study quantifies the equivalent height of freshwater, emphasizing the role it plays in the potential energy anomaly and front locations, as well as its relationship with river discharges. Using a criterion of equivalent height of freshwater >15% (density <1021 kg/m³ and salinity <28), the brackish layer was found to be 1–15 m thick (except in Estero Elefantes), with horizontal extensions on the order of 100 km. The limits of this layer tended to coincide with frontal zones having potential energy anomaly gradients >0.005 J/m⁴. The frontal zones were located in the extreme southeast of Jacaf Channel, at the head of Ventisquero Sound, in the central part of the Puyuguapi and Moraleda channels, and at the head and mouth of Aysén Fjord. The equivalent height of freshwater and potential energy anomaly showed a good correlation with the accumulated (5-day) river discharges (r²=0.87), which were greatest toward the fjord heads in spring. The brackish surface water had short residence times (3.5 days) in Aysen Fjord, unlike the deep layer, which other authors report to have a longer residence time (near 1 year).     

Scaling Issues in Forest Succession Modelling

This paper reviews scaling issues in forest succession modelling, focusing on forest gap models. Two modes of scaling are distinguished: (1) implicit scaling, i.e. taking scale-dependent features into account while developing model equations, and (2) explicit scaling, i.e. using procedures that typically involve numerical simulation to scale up the response of a local model in space and/or time. Special attention is paid to spatial upscaling methods, and downscaling is covered with respect to deriving scenarios of climatic change to drive gap models in impact assessments. When examining the equations used to represent ecological processes in forest gap models, it becomes evident that implicit scaling is relevant, but has not always been fully taken into consideration. A categorization from the literature is used to distinguish four methods for explicit upscaling of ecological models in space: (1) Lumping, (2) Direct extrapolation, (3) Extrapolation by expected value, and (4) Explicit integration. Examples from gap model studies are used to elaborate the potential and limitations of these methods, showing that upscaling to areas as large as 3000 km² is possible, given that there are no significant disturbances such as fires or insect outbreaks at the landscape scale. Regarding temporal upscaling, we find that it is important to consider migrational lags, i.e. limited availability of propagules, if one wants to assess the transient behaviour of forests in a changing climate, specifically with respect to carbon storage and the associated feedbacks to the atmospheric CO₂ content. Regarding downscaling, the ecological effects of different climate scenarios for the year 2100 were compared at a range of sites in central Europe. The derivation of the scenarios is based on (1) imposing GCM grid-cell average changes of temperature and precipitation on the local weather records; (2) a qualitative downscaling technique applied by the IPCC for central and southern Europe; and (3) statistical downscaling relating large-scale circulation patterns to local weather records. Widely different forest compositions may be obtained depending on the local climate scenario, suggesting that the downscaling issue is quite important for assessments of the ecological impacts of climatic change on forests.     

Remote sensing of stream depths with hydraulically assisted bathymetry (HAB) models

This article introduces a technique for using a combination of remote sensing imagery and open-channel flow principles to estimate depths for each pixel in an imaged river. This technique, which we term hydraulically assisted bathymetry (HAB), uses a combination of local stream gage information on discharge, image brightness data, and Manning-based estimates of stream resistance to calculate water depth. The HAB technique does not require ground-truth depth information at the time of flight. HAB can be accomplished with multispectral or hyperspectral data, and therefore can be applied over entire watersheds using standard high spatial resolution satellite or aerial images. HAB also has the potential to be applied retroactively to historic imagery, allowing researchers to map temporal changes in depth.We present two versions of the technique, HAB-1 and HAB-2. HAB-1 is based primarily on the geometry, discharge and velocity relationships of river channels. Manning's equation (assuming average depth approximates the hydraulic radius), the discharge equation, and the assumption that the frequency distribution of depths within a cross-section approximates that of a triangle are combined with discharge data from a local station, width measurements from imagery, and slope measurements from maps to estimate minimum, average and maximum depths at a multiple cross-sections. These depths are assigned to pixels of maximum, average, and minimum brightness within the cross-sections to develop a brightness–depth relation to estimate depths throughout the remainder of the river.HAB-2 is similar to HAB-1 in operation, but the assumption that the distribution of depths approximates that of a triangle is replaced by an optical Beer–Lambert law of light absorbance. In this case, the flow equations and the optical equations are used to iteratively scale the river pixel values until their depths produce a discharge that matches that of a nearby gage.R² values for measured depths versus depths estimated by HAB-1 and HAB-2 are 0.51 and 0.77, respectively, in the relatively simple Brazos River, Texas. R² values for HAB-1 and HAB-2 are 0.46 and 0.26, respectively, in the Lamar River, a complex mountain river system in Yellowstone National Park. Although the R² values are moderate, depth maps and cross-sections derived from the HAB techniques are consistent with typical stream geomorphology patterns and provide far greater spatial coverage and detail than could be achieved with ground-based survey techniques. Improved depth estimates can be achieved by stratifying the river into different habitat types that normalize for differences in turbulence and substrate.     

Sources,transport, and partitioning of organic matter at a highly dynamic continental margin

Continental shelves play a major role as transition zone during transport of multiply-sourced organic matter into the deep sea. In order to obtain a comprehensive understanding of the origin and fractionation processes of organic matter at the NW Iberian margin, 40 surface sediment samples were analyzed for a structurally diverse range of lipid biomarkers, lignin phenols, grain size distribution, organic carbon content (TOC), its stable carbon isotopic composition (δ¹³CTOC), and the organic carbon to nitrogen ratio (TOC/TN). The biomarker inventory reflected a heterogeneous mixture of organic matter from various marine and terrestrial sources. Soil- and vascular plant-derived continental organic matter, indicated by lignin phenols and plant-derived triterpenoids, was primarily associated with the silt fraction and transported by river run-off. The spatial distribution patterns of higher plant-derived waxes, long-chain n-alkanes, n-alcohols, and n-fatty acids suggested distinct different transport mechanisms and/or sources. The branched tetraether index, a molecular proxy expressing the relative abundance of branched dialkyl tetraethers vs. crenarchaeol and considered to signal soil-derived organic matter, was not as sensitive as the other molecular indicators in detecting continental organic matter. Hydrodynamic sorting processes on the shelf resulted in a separation of different types of terrestrial organic matter; grass and leaf fragments and soil organic matter were preferentially transported offshore and deposited in areas of lower hydrodynamic energy. Algal lipid biomarker distributions indicated a complex community of marine plankton contributing to organic matter. Spatial and seasonal patterns of phytoplankton growth primarily controlled the distribution of algal organic matter components. The interplay of all of these processes controls production, distribution, and deposition of organic matter and results in three distinct provinces at the Galicia–Minho shelf: (I) fresh marine organic matter dominated the inner shelf region; (II) high inputs of terrestrial organic matter and high TOC content characterized the mid-shelf deposited mudbelt; (III) lower concentrations of relatively degraded organic matter with increased proportions of refractory terrestrial components dominated the outer shelf and continental slope.     

Tuning a Kalman filter carrier tracking algorithm in the presence of ionospheric scintillation

GPS Solutions - Strong ionospheric electron content gradients may lead to fast and unpredictable fluctuations in the phase and amplitude of the signals from Global Navigation Satellite Systems...     

WIND observations of coherent electrostatic waves in the solar wind

The time domain sampler (TDS) experiment on WIND measures electric and magnetic wave forms with a sampling rate which reaches 120 000 points per second. We analyse here observations made in the solar wind near the Lagrange point L1. In the range of frequencies above the proton plasma frequency f_pi and smaller than or of the order of the electron plasma frequency f_pe, TDS observed three kinds of electrostatic (e.s.) waves: coherent wave packets of Langmuir waves with frequencies f ≃ f_pe, coherent wave packets with frequencies in the ion acoustic range f_pi ≥ f ≥ f_pe, and more or less isolated non-sinusoidal spikes lasting less than 1 ms. We confirm that the observed frequency of the low frequency (LF) ion acoustic wave packets is dominated by the Doppler effect: the wavelengths are short, 10 to 50 electron Debye lengths λ_D. The electric field in the isolated electrostatic structures (IES) and in the LF wave packets is more or less aligned with the solar wind magnetic field. Across the IES, which have a spatial width of the order of ≃25_D, there is a small but finite electric potential drop, implying an average electric field generally directed away from the Sun. The IES wave forms, which have not been previously reported in the solar wind, are similar, although with a smaller amplitude, to the weak double layers observed in the auroral regions, and to the electrostatic solitary waves observed in other regions in the magnetosphere. We have also studied the solar wind conditions which favour the occurrence of the three kinds of waves: all these e.s. waves are observed more or less continuously in the whole solar wind (except in the densest regions where a parasite prevents the TDS observations). The type (wave packet or IES) of the observed LF waves is mainly determined by the proton temperature and by the direction of the magnetic field, which themselves depend on the latitude of WIND with respect to the heliospheric current sheet.     

Groundwater recharge in a confined paleovalley setting,Northeast British Columbia,Canada

Hydrogeology Journal - Ancient river channels or subglacial drainage networks infilled with younger sediments can include significant deposits of highly permeable sands and gravels. Despite being...     

Towards generating highly detailed 3D CityGML models from OpenStreetMap

About one decade has passed since US vice president Al Gore articulated his vision of Digital Earth (DE). Within this decade, a global multi-resolution and three-dimensional (3D) representation of the Earth, which sums up the DE vision, increasingly gained interest in both public and science. Due to the desired high resolution of the available data, highly detailed 3D city models comprise a huge part of DE and they are becoming an essential and useful tool for a range of different applications. In the past as well as at present, 3D models normally come from a range of different sources generated by professionals, such as laser scans or photogrammetry combined with 2D cadaster data. Some models are generated with semi-automated or fully automated approaches, but in most cases manual fine tuning or even manual construction from architectural plans is required. Further beyond outdoor city models, DE additionally envisages the provision of indoor information. That is, the interior structure of public or publically accessible buildings, such as airports or shopping malls, is represented and made available in 3D; however, at the moment, such models are mostly created by hand and essentially based on professional data sources. In contrast to such professional data, which is mainly captured by surveyors or companies, the last few years revealed the phenomenon of crowdsourced geodata, which receives an increasing attractiveness as an alternative data source for many Geographic Information Systems (GIS). Former research already demonstrated the power and richness of such geodata – especially OpenStreetMap (OSM) – and it has also been proved that this non-standardized, crowdsourced geodata can be combined with international standards of the Open Geospatial Consortium (OGC). For example, CityGML Level-of-Detail 1 (LoD1) and LoD2 models have already been created automatically from OSM. The research presented in this article will further continue on the automated generation of CityGML models from OpenStreetMap. Essentially, a method for the creation of highly detailed CityGML LoD4 models with interior structures will be explained. By applying the invented approach on existing OSM data, limitations and restrictions of the IndoorOSM mapping proposal, the available data and the developed approach are revealed and discussed.