A hyperspectral experiment over tropical forests based on the EO-1 orbit change and PROSAIL simulation

ABSTRACT

We designed a unique hyperspectral experiment from the Earth Observing One (EO-1) orbit change to evaluate solar illumination effects over tropical forests in Brazil. Ten nadir-viewing Hyperion images collected over a fixed site and period of the year (July to August) were selected for analysis. We evaluated variations in reflectance and in 16 narrowband vegetation indices (VIs) with increasing solar zenith angle (SZA) from the pre-drift (2004–2008) to the EO-1 drift period (2011–2016). To detect changes in reflectance and shadows, we applied spectral mixture analysis (SMA) and principal component analysis (PCA) and calculated the similarity spectral angle (θ) between the vegetation spectra measured with variable SZA. The magnitude of the illumination effects was also evaluated from change-point analysis and nonparametric Mann-Whitney U tests applied over the time series. Finally, we complemented our experiment using the PROSAIL model to simulate the VIs variation with increasing SZA resultant from satellite drift. The results showed significant changes in Hyperion reflectance and VIs, especially when the EO-1 crossed the study area at earlier times and larger SZA in 2015 (9:05 a.m.; SZA = 59°) and 2016 (8:30 a.m.; SZA = 67°). Compared to the pre-drift period (10:30 a.m.; SZA = 45°), the SZA differences of 14° (2015) and 22° (2016) increased the shade fractions and decreased the vegetation brightness. PCA separated the pre-drift and drift reflectance datasets, showing shifts in scores due to changes in brightness. θ increased with SZA, indicating changes in the shape of the vegetation spectra with drift. For most VIs, the change-point analysis indicated 2015 (SZA = 59°) as the predominant year of detected changes. Compared to the EO-1 original orbit, the Plant Senescence Reflectance Index (PSRI), Anthocyanin Reflectance Index (ARI) and Structure Insensitive Pigment Index (SIPI) presented the largest positive changes during drift, while the Photochemical Reflectance Index (PRI), Visible Atmospherically Resistant Index (VARI) and Enhanced Vegetation Index (EVI) had the largest negative changes. The effect size of the illumination geometry on these VIs was large, as indicated by increasing values of the Cohen’s r metric toward 2016. The anisotropy of the Hyperion VIs was generally consistent with that from PROSAIL in the simulated pre-drift and drift periods. Focusing on structural indices, it affected the relationships between VIs and simulated leaf area index (LAI) at large SZA.     

Focal mechanisms of small earthquakes in the southeastern Brazilian shield: a test of stress models of the South American plate

Focal mechanisms of small earthquakes with magnitudes of about 3 in the SE Brazilian shield are calculated using S / P amplitude ratios. Low attenuation ( Q _p from 400 to 800) in the shield upper-crustal layers allowed sharp S arrivals to be recorded up to distances of 100 km. Besides P -wave polarities, SH -wave first motions were also used to constrain the nodal-plane orientations. Normal and reverse faulting mechanisms with strike-slip components were found. The inversion of four mechanisms to estimate the stress tensor indicated a strike-slip stress regime with roughly E–W-orientated σ ₁ and N–S σ ₃. Both the orientations and the shape factor ( φ =0.7) of the inverted stress are in excellent agreement with theoretical predictions for that part of Brazil from the driving-force model of Coblentz & Richardson (1996) . Good agreement with the nature of the stress, as well as its orientation, was also found for the model of Meijer (1995) . Both of these theoretical models include spreading stresses along the continent/ocean lithospheric transition. Because the earthquakes are more than 300 km from the continental shelf they should not be affected by the local flexural forces caused by sediment load in the marginal basins. The agreement between observed and theoretical stresses then confirms the importance of continental spreading forces in modelling intraplate stresses.     

Natural entropy production in an inflationary model for a polarized vacuum

Though entropy production is forbidden in standard FRW Cosmology, Berman and Som presented a simple inflationary model where entropy production by bulk viscosity, during standard inflation without ad hoc pressure terms can be accommodated with Robertson–Walker’s metric, so the requirement that the early Universe be anisotropic is not essential in order to have entropy growth during inflationary phase, as we show. Entropy also grows due to shear viscosity, for the anisotropic case. The intrinsically inflationary metric that we propose can be thought of as defining a polarized vacuum, and leads directly to the desired effects without the need of introducing extra pressure terms.     

The two pioneers anomalies and universal rotation

We show how to prove the two Pioneers Anomalies by means of the Godlowski et al. (, 2004) idea for a rotating General Relativistic Universe. The so-called clock effect is calculated.     

Combined influence of PDO and ENSO on northern Andean glaciers: a case study on the Cotopaxi ice-covered volcano,Ecuador

This paper describes the application of remote sensing in monitoring the fluctuations in one of the mountain glaciers in the Ecuadorean Andes during the past few decades using ASTER, EO-1 ALI, Landsat MSS, TM and ETM + images. Satellite images were used to calculate the snow line altitudes (SLAs) during the period 1979–2013. Cotopaxi ice covered volcano was studied as representative of Ecuadorian glaciers in the eastern cordillera. Precipitation and air temperature data from various gauging stations within the range of 30 km from the study site and monthly discharge and water level data from a gauging station were also utilized in this study. Anomalies in precipitation and temperature were found to be slightly different in the Cotopaxi region compared to nearby Antizana in the same cordillera and Chimborazo region in the western cordillera. An attempt to correlate the El Niño—southern oscillation phenomenon with the glacier fluctuations in Ecuadorian Andes was done successfully. Cold and warm regimes of Pacific Decadal Oscillation is also considered. The calculated glacier fluctuations obtained were similar to that performed on the nearby Antizana 15 in the eastern cordillera during 1995–2002. Precipitation and temperature anomalies were similar with Antizana 15. It is evident from the research that SLAs were highly fluctuated between the period of occurrence of El Niño and La Niña events. It is also seen that the glacier fluctuations show a negative mass balance trend in during the warm regime of Pacific Decadal Oscillation during the past three decades. Glaciated areas were advanced during the La Nina events in the cold regime of PDO during 1998–2002.     

Assessment of troposphere mapping functions using three-dimensional ray-tracing

The troposphere delay is an important source of error for precise GNSS positioning due to its high correlation with the station height parameter. It has been demonstrated that errors in mapping functions can cause sub-annual biases as well as affect the repeatability of GNSS solutions, which is a particular concern for geophysical studies. Three-dimensional ray-tracing through numerical weather models (NWM) is an excellent approach for capturing the directional and daily variation of the tropospheric delay. Due to computational complexity, its use for positioning purposes is limited, but it is an excellent tool for evaluating current state-of-the-art mapping functions used for geodetic positioning. Many mapping functions have been recommended in the past such as the Niell Mapping Function (NMF), Vienna Mapping Function 1 (VMF1), and the Global Mapping Function (GMF), which have been adopted by most IGS analysis centers. A new Global Pressure Temperature model (GPT2) has also been developed, which has been shown to improve upon the original atmospheric model used for the GMF. Although the mapping functions mentioned above use the same functional formulation, they vary in terms of their atmospheric source and calibration approach. A homogeneous data set of three-dimensional ray-traced delays is used to evaluate all components of the mapping functions, including their underlying functional formulation, calibration, and compression method. Additionally, an alternative representation of the VMF1 is generated using the same atmospheric source as the truth data set to evaluate the differences in ray-tracing methods and their effect on the end mapping function. The results of this investigation continue to support the use of the VMF1 as the mapping function of choice when geodetic parameters are of interest. Further support for the GPT2 and GMF as reliable back-ups when the VMF1 is not available was found due to their high consistency with the NWM-derived mapping function. Additionally, a small latitude-dependent bias in station height was found in the current mapping functions. This bias was identified to be due to the assumption of a constant radius of the earth and was largest at the poles and at the equator. Finally, an alternative version of the VMF1 is introduced, namely the UNB-VMF1 which provides users with an independent NWM-derived mapping function to support geodetic positioning.     

Center-to-Limb Variation of Solar Microwave Bursts

Statistical studies of hard X-ray flares position on the solar disk have shown that the more energetic hard X-rays have a tendency to be more concentrated near the limb rather than at disk center, whereas lower-energy hard X-ray emission seems isotropic. Since the high-frequency radio emission is believed to be produced by the same energetic electron population responsible for the high-energy hard X-rays, we searched the microwave/millimeter emitting bursts for center-to-limb variation in their emission. A total of 499 bursts observed by the radio telescopes in Bern at the frequencies of 3.1, 5.2, 8.4, 11.8, 19.6, 35.0, and 50.0 GHz were analyzed. Simultaneous Hα flares were used for determination of the radio burst position on the solar disk. For each of the 7 frequencies, the peak flux and duration were studied as a function of heliocentric position. For 312 bursts, spectral parameters such as spectral index, peak frequency, and flux at spectral maximum were analyzed. For a subset of 43 bursts with emission at all frequencies, the emission and spectral parameters were analyzed. Center-to-limb variations of the spectral parameters for all bursts were sought. In order to interpret the observational results, we have performed a numerical simulation of gyrosynchrotron spectra. We find that high-frequency events, which are also the more energetic ones, have larger center-to-limb variations in their parameters than do the overall flares. Moreover, this behavior agrees with theoretical predictions.     

From desert to Mediterranean rangelands: will increasing drought and inter-annual rainfall variability affect herbaceous annual primary productivity?

Climate change is predicted to alter the rainfall regime in the Eastern Mediterranean Basin: total annual rainfall will decrease, while seasonal and inter-annual variation in rainfall will increase. Such changes in the rainfall regime could potentially lead to large-scale changes in aboveground net primary productivity (ANPP) in the region. We conducted a data-driven evaluation of herbaceous ANPP along an entire regional rainfall gradient, from desert (90 mm MAR [Mean Annual Rainfall]) to Mesic-Mediterranean (780 mm MAR) ecosystems, using the largest database ever collated for herbaceous ANPP in Israel, with the aim of predicting consequences of climate change for rangeland productivity. This research revealed that herbaceous ANPP increases with increasing rainfall along the gradient, but strong dependence on rainfall was only apparent within dry sites. Rain Use Efficiency peaks at mid-gradient in Mediterranean sites without woody vegetation (560 and 610 mm MAR). Inter-annual coefficients of variation in rainfall and herbaceous ANPP decrease along the rainfall gradient up to ca. 500 mm MAR. Climate change is more likely to affect herbaceous ANPP of rangelands in the arid end of the rainfall gradient, requiring adaptation of rangeland management, while ANPP of rangelands in more mesic ecosystems is less responsive to variation in rainfall. We conclude that herbaceous ANPP in most Mediterranean rangelands is less vulnerable to climate change than generally predicted.     

Study of bi-orthogonal modes in magnetic butterflies

We present a bi-orthogonal decomposition of the temporal and latitudinal distribution of solar magnetic fields from synoptic magnetograms. Results are compared with a similar decomposition of the distribution of sunspots since 1874. We show that the butterfly diagrams can be interpreted as the result of approximately constant amplitudes and phases of two oscillations with periods close to 22 years. A clear periodicity of 7 years can also be identified in the most energetic modes of both spatio-temporal series. These results can be used to obtain relevant information concerning the physics of the solar dynamo.