Deriving remote sensing reflectance from turbid Case II waters using green-shortwave infrared bands based model

The objectives of this study are to validate the applicability of a shortwave infrared atmospheric correction model (SWIR-based model) in deriving remote sensing reflectance in turbid Case II waters, and to improve that model using a proposed green-shortwave infrared model (GSWIR-based model). In a GSWIR-based model, the aerosol type is determined by a SWIR-based model and the reflectance due to aerosol scattering is calculated using spectral slope technology. In this study, field measurements collected from three independent cruises from two different Case II waters were used to compare models. The results indicate that both SWIR- and GSWIR-based models can be used to derive the remote sensing reflectance at visible wavelengths in turbid Case II waters, but GSWIR-based models are superior to SWIR-based models. Using the GSWIR-based model decreases uncertainty in remote sensing reflectance retrievals in turbid Case II waters by 2.6–12.1%. In addition, GSWIR-based model’s sensitivity to user-supplied parameters was determined using the numerical method, which indicated that the GSWIR-based model is more sensitive to the uncertainty of spectral slope technology than to that of aerosol type retrieval methodology. Due to much lower noise tolerance of GSWIR-based model in the blue and near-infrared regions, the GSWIR-based model performs poorly in determining remote sensing reflectance at these wavelengths, which is consistent with the GSWIR-based model’s accuracy evaluation results.     

Short-arc tracklet association for geostationary objects

Measurement association and initial orbit determination is a fundamental task when building up a database of space objects. This paper proposes an efficient and robust method to determine the orbit using the available information of two tracklets, i.e. their line-of-sights and their derivatives. The approach works with a boundary-value formulation to represent hypothesized orbital states and uses an optimization scheme to find the best fitting orbits. The method is assessed and compared to an initial-value formulation using a measurement set taken by the Zimmerwald Small Aperture Robotic Telescope of the Astronomical Institute at the University of Bern. False associations of closely spaced objects on similar orbits cannot be completely eliminated due to the short duration of the measurement arcs. However, the presented approach uses the available information optimally and the overall association performance and robustness is very promising. The boundary-value optimization takes only around 2% of computational time when compared to optimization approaches using an initial-value formulation. The full potential of the method in terms of run-time is additionally illustrated by comparing it to other published association methods.     

The ionosphere and the Latin America VLF Network Mexico (LAVNet-Mex) station

In order to detect and study the ionospheric response to solar flares (transient high energy solar radiation), we have constructed a radio receiver station at Mexico City, which is part of the “Latin American Very low frequency Network” (LAVNet-Mex). This station extends to the northern hemisphere the so called “South American VLF Network”.     

First results of operational ionospheric dynamics prediction for the Brazilian Space Weather program

It is shown the development and preliminary results of operational ionosphere dynamics prediction system for the Brazilian Space Weather program. The system is based on the Sheffield University Plasmasphere–Ionosphere Model (SUPIM), a physics-based model computer code describing the distribution of ionization within the Earth mid to equatorial latitude ionosphere and plasmasphere, during geomagnetically quiet periods. The model outputs are given in a 2-dimensional plane aligned with Earth magnetic field lines, with fixed magnetic longitude coordinate. The code was adapted to provide the output in geographical coordinates. It was made referring to the Earth’s magnetic field as an eccentric dipole, using the approximation based on International Geomagnetic Reference Field (IGRF-11). During the system operation, several simulation runs are performed at different longitudes. The original code would not be able to run all simulations serially in reasonable time. So, a parallel version for the code was developed for enhancing the performance. After preliminary tests, it was frequently observed code instability, when negative ion temperatures or concentrations prevented the code from continuing its processing. After a detailed analysis, it was verified that most of these problems occurred due to concentration estimation of simulation points located at high altitudes, typically over 4000 km of altitude. In order to force convergence, an artificial exponential decay for ion–neutral collisional frequency was used above mentioned altitudes. This approach shown no significant difference from original code output, but improved substantially the code stability. In order to make operational system even more stable, the initial altitude and initial ion concentration values used on exponential decay equation are changed when convergence is not achieved, within pre-defined values. When all code runs end, the longitude of every point is then compared with its original reference station longitude, and differences are compensated by changing the simulation point time slot, in a temporal adjustment optimization. Then, an approximate neighbor searching technique was developed to obtain the ion concentration values in a regularly spaced grid, using inverse distance weighting (IDW) interpolation. A 3D grid containing ion and electron concentrations is generated for every hour of simulated day. Its spatial resolution is 1° of latitude per 1° of longitude per 10 km of altitude. The vertical total electron content (VTEC) is calculated from the grid, and plotted in a geographic map. An important feature that was implemented in the system is the capacity of combining observational data and simulation outputs to obtain more appropriate initial conditions to the ionosphere prediction. Newtonian relaxation method was used for this data assimilation process, where ionosonde data from four different locations in South America was used to improve the system accuracy. The whole process runs every day and predicts the VTEC values for South America region with almost 24 h ahead.     

Spectral transformation of ASTER and Landsat TM bands for lithological mapping of Soghan ophiolite complex,south Iran

Spectral transformation methods, including correlation coefficient (CC) and Optimum Index Factor (OIF), band ratio (BR) and principal component analysis (PCA) were applied to ASTER and Landsat TM bands for lithological mapping of Soghan ophiolitic complex in south of Iran. The results indicated that the methods used evidently showed superior outputs for detecting lithological units in ophiolitic complexes. CC and OIF methods were used to establish enhanced Red–Green–Blue (RGB) color combination bands for discriminating lithological units. A specialized band ratio (4/1, 4/5, 4/7 in RGB) was developed using ASTER bands to differentiate lithological units in ophiolitic complexes. The band ratio effectively detected serpentinite dunite as host rock of chromite ore deposits from surrounding lithological units in the study area. Principal component images derived from first three bands of ASTER and Landsat TM produced well results for lithological mapping applications. ASTER bands contain improved spectral characteristics and higher spatial resolution for detecting serpentinite dunite in ophiolitic complexes. The developed approach used in this study offers great potential for lithological mapping using ASTER and Landsat TM bands, which contributes in economic geology for prospecting chromite ore deposits associated with ophiolitic complexes.     

Langley method applied in study of aerosol optical depth in the Brazilian semiarid region using 500, 670 and 870 nm bands for sun photometer calibration

Due to the high costs of commercial monitoring instruments, a portable sun photometer was developed at INPE/CRN laboratories, operating in four bands, with two bands in the visible spectrum and two in near infrared. The instrument calibration process is performed by applying the classical Langley method. Application of the Langley’s methodology requires a site with high optical stability during the measurements, which is usually found in high altitudes. However, far from being an ideal site, Harrison et al. (1994) report success with applying the Langley method to some data for a site in Boulder, Colorado. Recently, Liu et al. (2011) show that low elevation sites, far away from urban and industrial centers can provide a stable optical depth, similar to high altitudes. In this study we investigated the feasibility of applying the methodology in the semiarid region of northeastern Brazil, far away from pollution areas with low altitudes, for sun photometer calibration. We investigated optical depth stability using two periods of measurements in the year during dry season in austral summer. The first one was in December when the native vegetation naturally dries, losing all its leaves and the second one was in September in the middle of the dry season when the vegetation is still with leaves. The data were distributed during four days in December 2012 and four days in September 2013 totaling eleven half days of collections between mornings and afternoons and by means of fitted line to the data V₀ values were found. Despite the high correlation between the collected data and the fitted line, the study showed a variation between the values of V₀ greater than allowed for sun photometer calibration. The lowest V₀ variation reached in this experiment with values lower than 3% for the bands 500, 670 and 870 nm are displayed in tables. The results indicate that the site needs to be better characterized with studies in more favorable periods, soon after the rainy season.     

Development of in-situ Space Debris Detector

Due to high relative velocities, collisions of spacecraft in orbit with Space Debris (SD) or Micrometeoroids (MM) can lead to payload degradation, anomalies as well as failures in spacecraft operation, or even loss of mission. Flux models and impact risk assessment tools, such as MASTER (Meteoroid and Space Debris Terrestrial Environment Reference) or ORDEM (Orbital Debris Engineering Model), and ESABASE2 or BUMPER II are used to analyse mission risk associated with these hazards. Validation of flux models is based on measured data. Currently, as most of the SD and MM objects are too small (millimeter down to micron sized) for ground-based observations (e.g. radar, optical), the only available data for model validation is based upon retrieved hardware investigations e.g. Long Duration Exposure Facility (LDEF), Hubble Space Telescope (HST), European Retrievable Carrier (EURECA). Since existing data sets are insufficient, further in-situ experimental investigation of the SD and MM populations are required. This paper provides an overview and assessment of existing and planned SD and MM impact detectors. The detection area of the described detectors is too small to adequately provide the missing data sets. Therefore an innovative detection concept is proposed that utilises existing spacecraft components for detection purposes. In general, solar panels of a spacecraft provide a large area that can be utilised for in-situ impact detection. By using this method on several spacecraft in different orbits the detection area can be increased significantly and allow the detection of SD and MM objects with diameters as low as 100 μm. The design of the detector is based on damage equations from HST and EURECA solar panels. An extensive investigation of those panels was performed by ESA and is summarized within this paper. Furthermore, an estimate of the expected sensitivity of the patented detector concept as well as examples for its implementation into large and small spacecraft are presented.     

Leading-edge flow separation control over an airfoil using a symmetrical dielectric barrier discharge plasma actuator

In order to promote an in-depth understanding of the mechanism of leading-edge flow separation control over an airfoil using a symmetrical Dielectric Barrier Discharge(DBD) plasma actuator excited by a steady-mode excitation, an experimental investigation of an SC(2)-0714 supercritical airfoil with a symmetrical DBD plasma actuator was performed in a closed chamber and a low-speed wind tunnel. The plasma actuator was mounted at the leading edge of the airfoil.Time-resolved Particle Image Velocimetry(PIV) results of the near-wall region in quiescent air suggested that the symmetrical DBD plasma actuator could induce some coherent structures in the separated shear layer, and these structures were linked to a dominant frequency of f0= 39 Hz when the peak-to-peak voltage of the plasma actuator was 9.8 kV. In addition, an analysis of flow structures without and with plasma actuation around the upper side of the airfoil at an angle of attack of18° for a wind speed of 3 m/s(Reynolds number Re = 20000) indicated that the dynamic process of leading-edge flow separation control over an airfoil could be divided into three stages. Initially, this plasma actuator could reinforce the shedding vortices in the separated shear layer. Then, these vortical structures could deflect the separated flow towards the wall by promoting the mixing between the outside flow with a high kinetic energy and the flow near the surface. After that, the plasma actuator induced a series of rolling vortices in the vicinity of the suction side of the airfoil, and these vortical structures could transfer momentum from the leading edge of the airfoil to the separated region, resulting in a reattachment of the separated flow around the airfoil.     

Assessing and minimizing collisions in satellite mega-constellations

We aim to provide satellite operators and researchers with an efficient means for evaluating and mitigating collision risk during the design process of mega-constellations. We first introduce a novel algorithm for conjunction prediction that relies on large-scale numerical simulations and uses a sequence of filters to greatly reduce its computational expense. We then use this brute-force algorithm to establish baselines of endogenous (intra-constellation), or self-induced, conjunction events for the FCC-reported designs of the OneWeb LEO and SpaceX Starlink mega-constellations. We demonstrate how these deterministic results can be used to validate more computationally efficient, stochastic techniques for close-encounter prediction by adopting a new probabilistic approach from Solar-System dynamics as a simple test case. Finally, we show how our methodology can be applied during the design phase of large constellations by investigating Minimum Space Occupancy (MiSO) orbits, a generalization of classical frozen orbits that holistically account for the perturbed-Keplerian dynamics of the Earth-satellite-Moon-Sun system. The results indicate that the adoption of MiSO orbital configurations of the proposed mega-constellations can significantly reduce the risk of endogenous collisions with nearly indistinguishable adjustments to the nominal orbital elements of the constellation satellites.     

Examining LiDAR – WorldView-3 data synergy to generate a detailed stand map in a mixed forest in the north-west of Turkey

Although stand delineation approach based on aerial photographs and field survey produces high accuracy maps, it is labour-intensive and time consuming. Furthermore, conventional forest stand maps may have some uncertainties that can hardly be verified due to the experiments and skills of photo-interpreters. Therefore, researchers have been seeking more objective and cost-effective methods for forest mapping. LiDAR (Light Detection and Ranging) data have a high potential to automatically delineate forest stands. Unlike optical sensors, LiDAR height data provides information about both the horizontal and vertical structural characteristics of forest stands. However, it deprives of spectral data that may be successfully used in separating tree species. In this study, we investigate the potential of LiDAR – WorldView-3 data synergy for the automatic generation of a detailed forest stand map which can be used for a tactical forest management plan. Firstly, image segmentation was applied to LiDAR data alone and LiDAR/WorldView-3 data set in order to obtain the most suitable image objects representing forest stands. Visual inspection of the segmentation results showed that image objects based on the LiDAR/WorldView-3 data set were more compatible with the reference forest stand boundaries. After the segmentation process, the LiDAR and LiDAR/WorldView-3 data sets were independently classified using object-based classification method. We tested two levels of classification. The first was a detailed classification with 14 classes considering reference stand types. The second was the rough classification with 9 classes where some stand types were combined. The mean, standard deviation and texture features of LiDAR metrics and spectral information were used in the classification. The accuracy assessment results of LiDAR data showed that the Overall Accuracy (OA) was calculated as 0.31 and 0.43, and the Kappa Index (KIA) was calculated as 0.26 and 0.32 for the detailed and rough classifications, respectively. For the LiDAR/WorldView-3 data set, the OA values were calculated as 0.50 and 0.61, and the KIA were calculated as 0.46 and 0.55 for the detailed and rough classifications, respectively. These results showed that the forest stand map derived from the LiDAR/WorldView-3 data synergy is more compatible with the reference forest stand map. In conclusion, it can be said that the forest stand maps produced in this study may provide strategic forest planning needs, but it is not sufficient for tactical forest management plans.