Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “time series analysis”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 991 records · Page 55

Dynamics of Fungal and Bacterial Biomass Carbon in Natural Ecosystems: Site-level Applications of the CLM-Microbe Model

Explicitly representing microbial processes has been recognized as a key improvement to Earth system models for the realistic projections of soil carbon (C) and climate dynamics. The CLM-Microbe model builds upon the CLM4.5 and explicitly represents two major soil microbial groups, fungi and bacteria. Based on the compiled time-series data of fungal (FBC) and bacterial (BBC) biomass C from nine biomes, we parameterized and validated the CLM-Microbe model, and further conducted sensitivity analysis and uncertainty analysis for simulating C cycling. The model performance was evaluated with mean absolute error (MAE), root mean square error (RMSE), and coefficient of determination (R2) for relative change in FBC and BBC. The CLM-Microbe model is able to reasonably capture the seasonal dynamics of FBC and BBC across biomes, particularly for tropical/subtropical forest, temperate broadleaf forest, and grassland, with MAE < 0.49 for FBC and <0.36 for BBC and RMSE <0.52 FBC and <0.39 for BBC, while R2 values are relatively smaller in some biomes (e.g., shrub) due to small sample sizes. We found good consistencies between simulated and observed FBC (R 2 =0.70, P<0.001) and BBC (R 2 =0.26, P<0.05) on average across biomes, but the model is not able to fully capture the large variation in observed FBC and BBC. Sensitivity analysis shows the most critical parameters are turnover rate, carbon-to-nitrogen ratio of fungi and bacteria, and microbial assimilation efficiency. This study confirms that the explicit representation of soil microbial mechanisms enhances model performance in simulating C variables such as heterotrophic respiration and soil organic C density. The further application of the CLM-Microbe model would deepen our understanding of microbial contributions to the global C cycle.

54 ENVIRONMENTAL SCIENCES↗

Combinatorial Evaluation of Physical Feature Engineering, Classical Machine Learning, and Deep Learning Models for Synchrophasor Data at Scale

A major objective of the project was to train and evaluate the effectiveness of multiple event and anomaly detection, identification and classification deep temporal learning models for processing of real-time phasor measurement unit (PMU) data streams. A vast dataset, consisting of two years of phasor measurements from all three U.S. Interconnections, was curated and released by the Department of Energy (DOE) through Pacific Northwest National Laboratory (PNNL). The dataset also included an event log that provided event times and types (e.g. generator trips, line trips, planned service events, transformer operations, etc.). Our analysis of this dataset addressed six (6) of the eleven (11) research priorities identified in Funding Opportunity Announcement (FOA) DE-FOA-0001861 “Big Data Analysis of Synchrophasor Data” (FOA 1861). Rather than being limited to pre-determined specific algorithms, this project relied on the uniquely structured, highly performant underlying time series database capabilities of the PredictiveGrid platform to assess the vast dataset utilizing a wide variety of algorithms.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Detection of Extreme Phenomena in the Stable Boundary Layer over the Amazonian Forest

We apply different methods for detection of extreme phenomena (EP) in air-turbulent time series measured in the nocturnal boundary layer above the Amazon forest. The methods used were: (a) a Morlet complex wavelet transform, which is often used in analysis of non-linear application processes. Through the use of the wavelet, it is possible to observe a phase singularity that involves a strong interaction between an extensive range of scales; (b) recurrence plot tests, which were used to identify a sudden change between different stable atmospheric states. (c) statistical analysis of early-warning signals, which verify simultaneous increases in the autocorrelation function and in the variance in the state variable; and (d) analysis of wind speed versus turbulent kinetic energy to identify different turbulent regimes in the stable boundary layer. We found it is adequate to use a threshold to classify the cases of strong turbulence regime, as a result of the occurrence of EP in the tropical atmosphere. All methods used corroborate and indicate synergy between events that culminate in what we classify as EP of the stable boundary layer above the tropical forest.

54 ENVIRONMENTAL SCIENCES↗

Multigrid one shot methods for optimal control problems: Infinite dimensional control

The multigrid one shot method for optimal control problems, governed by elliptic systems, is introduced for the infinite dimensional control space. ln this case, the control variable is a function whose discrete representation involves_an increasing number of variables with grid refinement. The minimization algorithm uses Lagrange multipliers to calculate sensitivity gradients. A preconditioned gradient descent algorithm is accelerated by a set of coarse grids. It optimizes for different scales in the representation of the control variable on different discretization levels. An analysis which reduces the problem to the boundary is introduced. It is used to approximate the two level asymptotic convergence rate, to determine the amplitude of the minimization steps, and the choice of a high pass filter to be used when necessary. The effectiveness of the method is demonstrated on a series of test problems. The new method enables the solutions of optimal control problems at the same cost of solving the corresponding analysis problems just a few times.

Arian, Eyal↗

Phytoplankton production in the Sargasso Sea as determined using optical mooring data

Optical measurements from an untended mooring provide high-frequency observations of in-water optical properties and permit the estimation of important biological parameters continuously as a function of time. A 9-month time series, composed of three separate deployments, of optical data from the BIOWATT 1987 deep-sea mooring located in the oligotrophic waters of the Sargasso Sea at 34 deg N, 70 deg W are presented. These data have been tested using several bio-optical models for the purpose of providing a continuous estimate of phytoplankton productivity. The data are discussed in the context of contemporaneous shipboard observations and for future ocean color satellite observations. We present a continuous estimation of phytoplankton productivity for the 9-month time series. Results from the first 70-day deployment are emphasized to demonstrate the utility of optical observations as proxy measures of biological parameters, to present preliminary analysis, and to compare our bio-optical observations with concurrent physical observations. The bio-optical features show variation in response to physical forcings including diel variations of incident solar irradiance, episodic changes corresponding to wind forcing, variability caused by advective mesoscale eddy events in the vicinity of the mooring, and seasonal variability corresponding to changes in solar radiation, shoaling of the mixed layer depth, and succession of phytoplankton populations.

Waters, K. J.↗

Session on High Speed Civil Transport Design Capability Using MDO and High Performance Computing

Since the inception of CAS in 1992, NASA Langley has been conducting research into applying multidisciplinary optimization (MDO) and high performance computing toward reducing aircraft design cycle time. The focus of this research has been the development of a series of computational frameworks and associated applications that increased in capability, complexity, and performance over time. The culmination of this effort is an automated high-fidelity analysis capability for a high speed civil transport (HSCT) vehicle installed on a network of heterogeneous computers with a computational framework built using Common Object Request Broker Architecture (CORBA) and Java. The main focus of the research in the early years was the development of the Framework for Interdisciplinary Design Optimization (FIDO) and associated HSCT applications. While the FIDO effort was eventually halted, work continued on HSCT applications of ever increasing complexity. The current application, HSCT4.0, employs high fidelity CFD and FEM analysis codes. For each analysis cycle, the vehicle geometry and computational grids are updated using new values for design variables. Processes for aeroelastic trim, loads convergence, displacement transfer, stress and buckling, and performance have been developed. In all, a total of 70 processes are integrated in the analysis framework. Many of the key processes include automatic differentiation capabilities to provide sensitivity information that can be used in optimization. A software engineering process was developed to manage this large project. Defining the interactions among 70 processes turned out to be an enormous, but essential, task. A formal requirements document was prepared that defined data flow among processes and subprocesses. A design document was then developed that translated the requirements into actual software design. A validation program was defined and implemented to ensure that codes integrated into the framework produced the same results as their standalone counterparts. Finally, a Commercial Off the Shelf (COTS) configuration management system was used to organize the software development. A computational environment, CJOPT, based on the Common Object Request Broker Architecture, CORBA, and the Java programming language has been developed as a framework for multidisciplinary analysis and Optimization. The environment exploits the parallelisms inherent in the application and distributes the constituent disciplines on machines best suited to their needs. In CJOpt, a discipline code is "wrapped" as an object. An interface to the object identifies the functionality (services) provided by the discipline, defined in Interface Definition Language (IDL) and implemented using Java. The results of using the HSCT4.0 capability are described. A summary of lessons learned is also presented. The use of some of the processes, codes, and techniques by industry are highlighted. The application of the methodology developed in this research to other aircraft are described. Finally, we show how the experience gained is being applied to entirely new vehicles, such as the Reusable Space Transportation System. Additional information is contained in the original.

Rehder, Joe↗

Aqua-planet simulations of the formation of the South Atlantic convergence zone

The impact of Amazon Basin convection and cold fronts on the formation and maintenance of the South Atlantic convergence zone (SACZ) is studied using aqua-planet simulations with a general circulation model. In the model, a circular patch of warm sea-surface temperature (SST) is used to mimic the effect of the Amazon Basin on South American monsoon convection. The aqua-planet simulations were designed to study the effect of the strength and latitude of Amazon Basin convection on the formation of the SACZ. The simulations indicate that the strength of the SACZ increases as the Amazon convection intensifies and is moved away from the equator. Of the two controls studied here, the latitude of the Amazon convection exerts the strongest effect on the strength of the SACZ. An analysis of the synoptic-scale variability in the simulations shows the importance of frontal systems in the formation of the aqua-planet SACZ. Composite time series of frontal systems that occurred in the simulations show that a robust SACZ occurs when fronts penetrate into the subtropics and become stationary there as they cross eastward of the longitude of the Amazon Basin. Moisture convergence associated with these frontal systems produces rainfall not along the model SACZ region and along a large portion of the northern model Amazon Basin. Simulations in which the warm SST patch was too weak or too close to the equator did not produce frontal systems that extended into the tropics and became stationary, and did not form a SACZ. In the model, the SACZ forms as Amazon Basin convection strengthens and migrates far enough southward to allow frontal systems to penetrate into the tropics and stall over South America. This result is in agreement with observations that the SACZ tends to form after the onset of the monsoon season in the Amazon Basin.

South American Monsoon↗

SAI VLBI Analysis Center Report 2012

This report presents an overview of the SAI VLBI Analysis Center activities during 2012 and the plans for 2013. The SAI AC analyzes all IVS sessions for computations of the Earth orientation parameters (EOP) and time series of the ICRF source positions and performs research and software development aimed at improving the VLBI technique.

Zharov, Vladimir↗

Effect of MODIS Terra Radiometric Calibration Improvements on Collection 6 Deep Blue Aerosol Products: Validation and Terra/Aqua Consistency

The Deep Blue (DB) algorithm's primary data product is midvisible aerosol optical depth (AOD). DB applied to Moderate Resolution Imaging Spectroradiometer (MODIS) measurements provides a data record since early 2000 for MODIS Terra and mid-2002 for MODIS Aqua. In the previous data version (Collection 5, C5), DB production from Terra was halted in 2007 due to sensor degradation; the new Collection 6 (C6) has both improved science algorithms and sensor radiometric calibration. This includes additional calibration corrections developed by the Ocean Biology Processing Group to address MODIS Terra's gain, polarization sensitivity, and detector response versus scan angle, meaning DB can now be applied to the whole Terra record. Through validation with Aerosol Robotic Network (AERONET) data, it is shown that the C6 DB Terra AOD quality is stable throughout the mission to date. Compared to the C5 calibration, in recent years the RMS error compared to AERONET is smaller by approximately 0.04 over bright (e.g., desert) and approximately 0.01-0.02 over darker (e.g., vegetated) land surfaces, and the fraction of points in agreement with AERONET within expected retrieval uncertainty higher by approximately 10% and approximately 5%, respectively. Comparisons to the Aqua C6 time series reveal a high level of correspondence between the two MODIS DB data records, with a small positive (Terra-Aqua) average AOD offset <0.01. The analysis demonstrates both the efficacy of the new radiometric calibration efforts and that the C6 MODIS Terra DB AOD data remain stable (to better than 0.01 AOD) throughout the mission to date, suitable for quantitative scientific analyses.

Collection 6 Deep Blue↗

Sensitivity of Inherent Optical Properties from Ocean Reflectance Inversion Models to Satellite Instrument Wavelength Suites

The Earth science community seeks to develop climate data records (CDRs) from satellite measurements of ocean color, a continuous data record that now exceeds 20 years. Space agencies will launch additional instruments in the coming decade that will continue this data record, including the NASA PACE spectrometer. Inherent optical properties (IOPs) quantitatively describe the absorbing and scattering constituents of seawater and can be estimated from satellite-observed spectroradiometric data using semi-analytical algorithms (SAAs). SAAs exploit the contrasting optical signatures of constituent matter at spectral bands observed by satellite sensors. SAA performance, therefore, depends on the spectral resolution of the satellite spectroradiometer. A CDR spanning SeaWiFS, MODIS, OLCI, and PACE, for example, would include IOPs derived using varied wavelength suites if all available wavelengths were considered. Here, we explored differences in derived IOPs that stem simply from the use of (eight) different wavelength suites of input radiometric measurements. Using synthesized data and SeaWiFS Level-3 mission-long composites, we demonstrated equivalent SAA performance for all wavelength suites, but that IOP retrievals vary by several percent across wavelength suites and as a function of water type. The differences equate to roughly ≤ 6, 12, and 7% for a(sub dg)(443), a(sub ph)(443), and b(sub bp)(443), respectively, for waters with C(sub a) ≤ 1 mg -cu m. These values shrink for sensors with similar wavelength suites (e.g., SeaWiFS, MODIS, and MERIS) and rise to substantially larger values for higher C(sub a) waters. Our results also indicate that including 400 nm (in the case of OLCI) influences the derived IOPs, using longer wavelengths (>600 nm) influences the derived IOPs when there is a red signal, and, including additional spectral information shows potential for improved IOP estimation, but not without revisiting SAA parameterizations and execution. While modest in scope, we believe this study contributes to the knowledge base for CDR development. The implication of ignoring such an analysis as CDRs continue to be developed is a prolonged inability to distinguish between algorithmic and environmental contributions to trends and anomalies in the IOP time-series.

ocean color satellites↗

Initial Validation of NDVI time seriesfrom AVHRR, VEGETATION, and MODIS

The paper will address Theme 7: Multi-sensor opportunities for VEGETATION. We present analysis of a long-term vegetation record derived from three moderate resolution sensors: AVHRR, VEGETATION, and MODIS. While empirically based manipulation can ensure agreement between the three data sets, there is a need to validate the series. This paper uses atmospherically corrected ETM+ data available over the EOS Land Validation Core Sites as an independent data set with which to compare the time series. We use ETM+ data from 15 globally distributed sites, 7 of which contain repeat coverage in time. These high-resolution data are compared to the values of each sensor by spatially aggregating the ETM+ to each specific sensors' spatial coverage. The aggregated ETM+ value provides a point estimate for a specific site on a specific date. The standard deviation of that point estimate is used to construct a confidence interval for that point estimate. The values from each moderate resolution sensor are then evaluated with respect to that confident interval. Result show that AVHRR, VEGETATION, and MODIS data can be combined to assess temporal uncertainties and address data continuity issues and that the atmospherically corrected ETM+ data provide an independent source with which to compare that record. The final product is a consistent time series climate record that links historical observations to current and future measurements.

Morisette, Jeffrey T.↗

Assessing Relay Communications for Mars Sample Return Surface Mission Concepts

The Mars Sample Return (MSR) Campaign would be a 3-mission campaign concept supported by NASA and ESA to return samples from the Mars surface. MSR would, for the first time ever, present a need to communicate with multiple surface assets that are co-located on Mars in a coordinated effort to accomplish the unified objective of fetching, transporting, and returning samples from Mars. Currently, Mars surface assets relay data to and from Earth using a number of orbiters in what’s known as the Mars Relay Network (MRN). This network is characterized by a small number of surface assets distributed across the Martian globe and a larger number of orbiters to provide relay services. As of June 2020, there are two surface assets for which five orbiters are providing relay. During the MSR Campaign, there would be two rovers and a lander that all would require relay communication from a small number of Mars orbiters to meet the aggressive MSR timeline. The inversion of the current MRN paradigm, a system of many surface assets requiring relay and few orbiters to provide relay, necessitates the unique challenge of optimally allocating relay passes to maximize the operational capability of all assets. The allocation must consider a large number of trade variables including Mars asset operational requirements and Earth ground system constraints, including staffing schedules, operations planning across time zones, and more. To address these telecommunication challenges, the Mars Asset Relay Mission Link Allocation Design Environment (MARMLADE) tool was developed. It is a MATLAB-based tool to assign orbiter passes or Direct-From-Earth (DFE) links to each of the three surface assets and quantify the operational efficiency of each surface asset.MARMLADE uses a data set of simulated Mars relay orbiter geometry and telecommunication capabilities provided by JPL’s Telecom Orbit Analysis and Simulation Tool (TOAST) software to compute which asset should get each pass based on a series of heuristics and predictions of all assets’ states. Within MARMLADE, the user can provide inputs including the option for time-based pass splitting, fixed FWD data rate capabilities, DFE communication capabilities, and link parameters allowing for the assessment of complex operations and hardware trades using surface mission operational efficiency as a primary figure of merit. As the MSR mission concepts continue to mature, MARMLADE is being used to assess ability of all MSR elements to meet the surface mission timeline requirements and to provide relay link allocations to each of the MSR surface assets.This paper will describe the motivation and design of the MARMLADE tool and how it is being used to perform campaign and mission level trades, generate requirements, and support development of the MSR surface mission scenarios.

Lee, Charles↗

Automated detection of photovoltaic cleaning events: A performance comparison of techniques as applied to a broad set of labeled photovoltaic data sets

Extracting accurate soiling loss information from photovoltaic (PV) production data first requires segmenting the time series data per natural or manually occurring cleaning events. Maintenance logs are often incomplete, rain data are often unavailable, and the debate on rain thresholds for cleaning and dew or wind cleanings is still ongoing. The present work aims to overtake these issues by improving automated methods to detect these cleaning events and therefore improve extraction of soiling loss information. Time series power production data from 22 PV inverters were labeled for natural or manually occurring cleaning events. The data sets were carefully selected to include varying degrees of soiling, cleaning events, and noise. Several algorithms, including filtering logic and change point detection, were examined for efficacy at detecting the labeled cleanings. All the methods introduced except for changepoint detection showed significant improvement at detecting the labeled cleaning events per the mean F 1 score. Furthermore, the highest performing cleaning detection algorithm achieved an absolute increase in the mean F 1 score of 43% over the default version of the RdTools stochastic rate and recovery (SRR) algorithm. The highest performing algorithm included irradiance filtering and a cleaning detection threshold, adjusted based on the 40-day centered rolling median of the absolute day-to-day deviations in the daily performance index (PI). Furthermore, these improvements are promising as cleaning detection is an essential step in the automated analysis of PV soiling.

14 SOLAR ENERGY↗

MOBY, A Radiometric Buoy for Performance Monitoring and Vicarious Calibration of Satellite Ocean Color Sensors: Measurement and Data Analysis Protocols

The Marine Optical Buoy (MOBY) is the centerpiece of the primary ocean measurement site for calibration of satellite ocean color sensors based on independent in situ measurements. Since late 1996, the time series of normalized water-leaving radiances L(sub WN)(lambda) determined from the array of radiometric sensors attached to MOBY are the primary basis for the on-orbit calibrations of the USA Sea-viewing Wide Field-of-view Sensor (SeaWiFS), the Japanese Ocean Color and Temperature Sensor (OCTS), the French Polarization Detection Environmental Radiometer (POLDER), the German Modular Optoelectronic Scanner on the Indian Research Satellite (IRS1-MOS), and the USA Moderate Resolution Imaging Spectrometer (MODIS). The MOBY vicarious calibration L(sub WN)(lambda) reference is an essential element in the international effort to develop a global, multi-year time series of consistently calibrated ocean color products using data from a wide variety of independent satellite sensors. A longstanding goal of the SeaWiFS and MODIS (Ocean) Science Teams is to determine satellite-derived L(sub WN)(labda) with a relative combined standard uncertainty of 5 %. Other satellite ocean color projects and the Sensor Intercomparison for Marine Biology and Interdisciplinary Oceanic Studies (SIMBIOS) project have also adopted this goal, at least implicitly. Because water-leaving radiance contributes at most 10 % of the total radiance measured by a satellite sensor above the atmosphere, a 5 % uncertainty in L(sub WN)(lambda) implies a 0.5 % uncertainty in the above-atmosphere radiance measurements. This level of uncertainty can only be approached using vicarious-calibration approaches as described below. In practice, this means that the satellite radiance responsivity is adjusted to achieve the best agreement, in a least-squares sense, for the L(sub WN)(lambda) results determined using the satellite and the independent optical sensors (e.g. MOBY). The end result of this approach is to implicitly absorb unquantified, but systematic, errors in the atmospheric correction, incident solar flux, and satellite sensor calibration into a single correction factor to produce consistency with the in situ data.

Clark, Dennis K.↗

Data-Driven Learning for the Mori--Zwanzig Formalism: A Generalization of the Koopman Learning Framework

A theoretical framework which unifies the conventional Mori--Zwanzig formalism and the approximate Koopman learning of deterministic dynamical systems from noiseless observation is presented. In this framework, the Mori--Zwanzig formalism, developed in statistical mechanics to tackle the hard problem of construction of reduced-order dynamics for high-dimensional dynamical systems, can be considered as a natural generalization of the Koopman description of the dynamical system. We next show that, similar to the approximate Koopman learning methods, data-driven methods can be developed for the Mori--Zwanzig formalism with Mori's linear projection operator. We have developed two algorithms to extract the key operators, the Markov and the memory kernel, using time series of a reduced set of observables in a dynamical system. We have adopted the Lorenz `96 system as a test problem and solved for the above operators. These operators exhibit complex behaviors, which are unlikely to be captured by traditional modeling approaches in Mori--Zwanzig analysis. The nontrivial generalized fluctuation-dissipation relationship, which relates the memory kernel with the two-time correlation statistics of the orthogonal dynamics, was numerically verified as a validation of the solved operators. Here we present numerical evidence that the generalized Langevin equation, a key construct in the Mori--Zwanzig formalism, is more advantageous in predicting the evolution of the reduced set of observables than the conventional approximate Koopman operators.

97 MATHEMATICS AND COMPUTING↗

Atmospheric condition identification in multivariate data through a metric for total variation

Identification of atmospheric conditions within a multivariable atmospheric data set is a necessary step in the validation of emerging and existing high-fidelity models used to simulate wind plant flows and operation.Atmospheric conditions relevant for wind energy research include stationary conditions, given the need for well-converged statistics for model validation, as well as conditions observed less frequently, such as extreme atmospheric events, which are used in wind turbine and wind plant design.Aggregation of observations without regard to covariance between time series discounts the dynamical nature of the atmosphere and is not sufficiently representative of atmospheric conditions.Identification and characterization of continuous time periods with atmospheric conditions that have a high value for analysis or simulation set the stage for more advanced model validation and the development of real-time control and operational strategies.The current work explores a single metric for variation in a multivariate data sample that quantifies variability within each channel as well as covariance between channels.The total variation is used to identify conditions of interest that conform to desired objective functions, such as stationary conditions, ramps or waves of wind speed, and changes in wind direction.Total variation is somewhat sensitive to the presence of outliers in the input data, and the method is best complemented by quality-control procedures to ensure reliable results.The direct detection and classification of events or conditions of interest within atmospheric data sets is vital to developing our understanding of wind plant response and to the formulation of forecasting and control models.

17 WIND ENERGY↗

Analysis of Atmospheric Mesoscale Models for Entry, Descent and Landing

Each Mars Exploration Rover (MER) is sensitive to the martian winds encountered near the surface during the Entry, Descent and Landing (EDL) process. These winds are strongly influenced by local (mesoscale) conditions. In the absence of suitable wind observations, wind fields predicted by martian mesoscale atmospheric models have been analyzed to guide landing site selection. Two different models were used, the MRAMS model and the Mars MM5 model. In order to encompass both models and render their results useful to the EDL engineering team, a series of statistical techniques were applied to the model results. These analyses cover the high priority landing sites during the expected landing times (1200 to 1500 local time). The number of sites studied is limited by the computational and analysis cost of the mesoscale models.

D M Kass↗

Time and Measurement Days

Questions in data analysis involving the concepts of time and measurement are often pushed into the background or reserved for a philosophical discussion. Some examples are: a) Is causality a consequence of the laws of physics, or can the arrow of time be reversed? b) Can we determine the arrow of time of an event? c) Do we need the continuum hypothesis for the underlying function in any measurement process? d) Can we say anything about the analyticity of the underlying process of an event? e) Would it be valid to model a non-analytical process as function of time? f) What are the implications of all these questions for classical Fourier techniques? However, in the age of big data gathered either from space missions supplying ultra-precise long time series, or e.g. LIGO data from the ground, the moment to bring these questions to the foreground seems arrived. The limitations of our understanding of some fundamental processes is emphasized by the lack of solution for problems open for more than 2 decades, such as the non-detection of solar g-modes, or the modal identification of main sequence stellar pulsators like delta Scuti stars. Flicker noise or 1/f noise, for example, attributed in solar-like stars to granulation, is analyzed mostly only to apply noise reduction techniques, neither considering the classical problem of 1/f noise that was introduced a 100 years ago, nor taking into account ergodic or non-ergodic solutions that make inapplicable spectral analysis techniques in practice. This topic was discussed by Nicholas W. Watkins during the ITISE meeting held in Granada in 2016. There he presented preliminary results of his research on Mandelbrot's related work. We reproduce here his quotation of Mandelbrot (1999) "There is a sharp contrast between a highly anomalous ("non-white") noise that proceeds in ordinary clock time and a noise whose principal anomaly is that it is restricted to fractal time", suggesting a connection with the above proposed topics that could be phrased as the following additional questions:a) Is self-organized criticality (SOC) frequent in astrophysical phenomena? b) Could all fractals in nature be considered stochastic? c) Could we establish mathematical/physical relationships between chaotic and fractal behaviors in time series? d) Could the differences between fractals and chaos in terms of analyticity be used to understand the residuals of the fitting of stellar light curves? In this meeting we would like to approximate these problems from a holistic and multidisciplinary perspective, taking into account not only technical issues but also the deeper implications. In particular the concept of connectivity (introduced in Pascual-Granado et al. A&A, 2015) could be used to implement, within the framework of ARMA processes, an "arrow of time" (see attached document), and so studying the possible implications in the concept of time as envisaged by Watkins.

data analysis↗