Data Mining Applied to the Electrochemical Noise Technique in the Time/Frequency Domain for Stress Corrosion Cracking Recognition

: In this paper, time/frequency domain data processing was proposed to analyse the EN signal recorded during stress corrosion cracking on precipitation-hardening martensitic stainless steel in a chloride environment. Continuous Wavelet Transform, albeit with some limitations, showed a suitable support in the discriminatory capacity among transient signals related to the different stress corrosion cracking mechanisms. In particular, the aim is to propose the analysis of electro-chemical noise signals under stress corrosion cracking conditions in the time–frequency domain by using the Hilbert–Huang approach. The Hilbert–Huang Transform (performed by the Empirical Mode Decomposition approach) was ﬁnally proposed to carry out an identiﬁcation of the corrosion mechanisms in comparison to conventional data processing methods. By using this approach, a detailed simultaneous decomposition of the original electrochemical noise data in the time and frequency domain was carried out. The method gave useful information about transitions among different corrosion mechanisms, allowing us to (i) identify a speciﬁc characteristic response for each corrosion damaging phenomenon induced by stress corrosion cracking, (ii) time each corrosion of the damaging phenomenon, and (iii) provide a topological description of the advancing SCC damaging stages. This characteristic evidences that the Hilbert–Huang Transform is a very powerful technique to potentially recognize and distinguish the different corrosion mechanisms occurring during stress corrosion cracking.


Introduction
Stress corrosion cracking (SCC) represents one of the most crucial corrosion processes that can result in early catastrophic failures of structural parts or system, significantly reducing its life cycle.As a result, the adoption of a non-intrusive in situ assessment technique looks to be an intriguing way for obtaining knowledge and determining definite parameters in SCC procedures.Furthermore, various attempts are being made in the scientific community to develop NDT methods for monitoring stressed parts subjected to corrosive environments [1][2][3][4].
Among the several electrochemical techniques, electrochemical noise (EN) could potentially be considered most suitable for the assessment of damage evolution during SCC.Because this technology relies on spontaneous current and potential fluctuations created by corrosion processes, one of its primary advantages is its non-invasiveness [5,6].
Electrochemical noise (EN) can be applied to monitor and identify, in situ, the corrosion processes of metallic materials [7], and additionally monitor and assess the type and rate of corrosion in different harmful environments [8].It can also be used to monitor the durability of coated metals [6,9].Furthermore, this technique can detect localized corrosion, which affects small areas of the metal surface, as well as uniform corrosion, which affects the whole surface [10].
Indeed, various research efforts have shown that by using EN measurements, it is feasible to identify and discriminate between distinct forms of localized corrosion.[10], such as pitting [11], crevice [12], uniform corrosion [13], and stress-corrosion cracking [14].
In particular, as a reference, but not limited to, in [11], the stable and metastable pitting corrosion behaviour, pit initiation and growth in SLM 316L SS, via wavelet and stochastic analysis of ECN signals, were assessed.
Furthermore, in [12], localised crevice corrosion was investigated.Given that the mechanistic evaluation of localized corrosion mechanisms, with regard to electrochemical variables influencing the corrosion phenomena, is a critical point of discrimination, a combination of electrochemical impedance spectroscopy (EIS) and electrochemical noise (EN) approach using the Hilbert analysis was advantageously used.
Similarly, excellent potential for analysis of uniform corrosion using EN has been highlighted in several research activities [13], for which the uniform corrosion rate has been quantified by combining statistical characteristics and frequency domain analysis.
Analogously, the EN technique can be effectively applied to a complex phenomenon, such as stress corrosion cracking for a valuable detection of damage phenomena in SCC endorsing environment [14].
In this concern, studies performed on stainless steel under various conditions [15,16] showed that on SCC tests, it may be feasible to quantify electrochemical current and voltage noise variations caused by the cyclical activation and propagation of cracks.In a lot of research activities [17][18][19], hydrogen embrittlement and the impact of the related evolved hydrogen bubbles on the EN transient was accurately investigated.
SCC can form by either continuous and discontinuous mechanisms, characterized by quite distinct EN features.The discontinuous crack growth phenomenon is able to generate some distinctive current transients, whereas the continuous crack growth approach might generate low noise levels [20].Therefore, EN has a substantial advantage over other methods for understanding the mechanism of SSC.
Although, it is difficult to link the crack tip activity with the fracture propagation rate and the EN response data, as observed in [21], the technique verifies the estimation of stress-corrosion thresholds and, considering that EN signal increased with activity at the crack tip, can be used to analyse the influence factors have on cracking.
An EN signal could be analysed in the time domain (by means of visual analysis, statistical analysis, such as maxima and minima recognition, or identification of statistical moments and other statistical parameters) or in the frequency domain (by means of Fast Fourier Transform (FFT), Power Spectral Density (PSD), Hilbert Transform (HT), Wavelet Transform (WT) and many other techniques).
To know beforehand whether a processed signal is stationary and linear is generally difficult.Thus, there is a risk that not all appropriate signal analysis methods could be used, therefore obtaining non-reliable results.
To limit this problem, the use of transforms such as FFT and PSD, which in fact require stationary and linear conditions, are applied for the limited signal time window, whose size is chosen to ensure the hypothesis of constant signal amplitude and linearity [22,23].Therefore, the modified analysis becomes the so-called Short-term FT (STFT) [24][25][26], which shows critical issues on the choice of time window size, and, generally, provides a compromise between the assumption of stationariness and frequency resolution.Furthermore, typically acquired signals are a superposition of different signals from different sources (electrochemical dissolution, temperature drift, bubbles evolution...), each one with different intervals of stationary, so it is difficult to discriminate the correct values of each corrosion process.
Taking into account the aforementioned factors, traditional EN analysis approaches are insufficient to allow a clear discrimination among the different damage stages during the corrosion evolution in a component or system.Therefore, the use of a powerful timedependent signal processing technique for data analysis should be evaluated.
The constant wavelet transform (CWT) was introduced to overcome the principal FFT theoretical limitation [27]; it furnishes a native energy/time/frequency representation, and it could be applied to the non-stationary signals due to its regional convolution mechanism.This method was successfully applied to the load variation analysis [24] and to the corrosion field [27].However, it maintains a theoretical limitation on the applicability to non-linear signals, and, furthermore, there is an approximation related to the frequency/time resolution that should be taken into account.
A promising analysis method for non-stationary and non-linear signals of recorded EN data, which can be used to substitute or corroborate traditional time-dependent or frequency-dependent approaches, is the Hilbert-Huang Transform (HHT) [28].The ability to analyse non-stationary and non-linear signals is the main advantage of HHT compared to ordinary signal analysis methods [29].This relatively novel analysis technique was recently applied with promising results for electrochemical noise data to obtain trend removal and a characterization of EN data in some simple corrosion studies [30][31][32][33] With this in mind, the aim of this paper is to address the gap in knowledge about SCC mechanisms by integrating different methods of EN signal analysis.These methods include traditional time and frequency (PSD) analysis, Wavelet Transform analysis, and the HHT technique.This comprehensive analysis approach can distinguish each relevant transient event that corresponds to a specific SCC damage mechanism.In particular, the HHT approach was applied to electrochemical noise signals acquired during stress corrosion cracking tests of a martensitic stainless steel.Based on this technique, coupled with conventional electrochemical noise analysis techniques, it is possible to identify a more precise topological damage map which might be used to monitor the damage progression of a real structure over time, identifying or avoiding its eventual failure.

Materials and Testing Conditions
Stress Corrosion Cracking (SCC) specimens were manufactured of 17-4 PH (UNS S17400) stainless steel.17-4 PH is a martensitic stainless steel with more Cr and Ni than the standard AISI 4xx series for increased corrosion resistance, as well as an extra alloying element(s) for precipitation hardening.(Table 1).Table 1 summarizes the 17-4 PH stainless steel nominal chemical composition obtained by X-ray fluorescence (XRF, Rigaku Supermini 200, Tokyo, Japan) analysis.[14] to obtain the mechanical properties reported in Table 2. SCC tests, in agreement with the ASTM G36 standard, were performed in a water-based solution containing 30% magnesium chloride by weight at a temperature of 100 • C.
For the application of tensile stress to specimens during stress corrosion testing, a dead-weight type machine was used.The stress exerted to the samples during the SCC test was 90% of the material's actual yield strength (745 MPa).
To allow the arrangement of the electrodes necessary for the collection of current and potential noise, a suitable cell was created by adapting the well-known design typically used for stress corrosion tests using dog-bone form specimens (i.e., NACE TM0177 technique A).The cell was made up of a Pyrex ring and two Teflon flanges with a central hole for holding the specimen.A three-electrode system was used to measure current and potential fluctuations simultaneously during EN measurements.The electrodes in this system were all made of the same material [34].The electrodes were placed in a parallel pattern about 4 cm apart.The EN current between the working electrode and one of the cylindrical electrodes was measured using a zero-resistance ammeter (ZRA).Simultaneously, the EN voltage between the working electrode and the other cylindrical electrode was measured using a high-impedance voltmeter (1 Tera Ohm), via a PAR VersaSTAT 4 potentiostat.The stressed specimen's effective working surface (working electrode) in contact with the test solution was about 24.1 cm 2 .The sampling rate was 1 Hz.A scheme of the SCC test setup is reported in Figure 1.SCC tests, in agreement with the ASTM G36 standard, were performed in a waterbased solution containing 30% magnesium chloride by weight at a temperature of 100 °C.
For the application of tensile stress to specimens during stress corrosion testing, a dead-weight type machine was used.The stress exerted to the samples during the SCC test was 90% of the material's actual yield strength (745 MPa).
To allow the arrangement of the electrodes necessary for the collection of current and potential noise, a suitable cell was created by adapting the well-known design typically used for stress corrosion tests using dog-bone form specimens (i.e., NACE TM0177 technique A).The cell was made up of a Pyrex ring and two Teflon flanges with a central hole for holding the specimen.A three-electrode system was used to measure current and potential fluctuations simultaneously during EN measurements.The electrodes in this system were all made of the same material [34].The electrodes were placed in a parallel pattern about 4 cm apart.The EN current between the working electrode and one of the cylindrical electrodes was measured using a zero-resistance ammeter (ZRA).Simultaneously, the EN voltage between the working electrode and the other cylindrical electrode was measured using a high-impedance voltmeter (1 Tera Ohm), via a PAR VersaSTAT 4 potentiostat.The stressed specimen's effective working surface (working electrode) in contact with the test solution was about 24.1 cm 2 .The sampling rate was 1 Hz.A scheme of the SCC test setup is reported in Figure 1.The EN acquisition covered the full period of the tests.The performed experimental tests were phenomenologically comparable.The fracture times of the tests varied in the range of 45,000-68,000 s.A scheme of the proposed damage analysis approach of EN data is reported in Figure 2. Following the removal of corrosion products, the specimens were examined using a scanning electron microscope (SEM), specifically, the ZEISS crossbeam 540 FESEM (Zeiss Microscopy, Jena, Germany).The EN acquisition covered the full period of the tests.The performed experimental tests were phenomenologically comparable.The fracture times of the tests varied in the range of 45,000-68,000 s.A scheme of the proposed damage analysis approach of EN data is reported in Figure 2. Following the removal of corrosion products, the specimens were examined using a scanning electron microscope (SEM), specifically, the ZEISS crossbeam 540 FESEM (Zeiss Microscopy, Jena, Germany).

Hilbert-Huang Transform Approach
HHT is a way to analyse the signal through its Intrinsic Mode Functions (IMFs).IMFs run in a simple oscillatory mode with varying amplitude and frequency along a t axis [35,36], and are used to obtain instantaneous frequencies at each moment in time w a symmetrical wave profile.The decomposition of the signal into its IMFs is obtained ing the Empirical Mode Decomposition (EMD) method.
Empirical Mode Decomposition (EMD) [29,36] is an empirical method that expl the definition of IMF to carry out the decomposition: a function with an equal number with a difference of 1) of extremes and zero crossing, whose envelopes are symmetric w respect to zero (i.e., at any point the average of the signal envelope determined by the lo maxima and the envelope determined by the local minima is zero).It represents a gen ally simple oscillatory mode.
So, the EMD algorithm is composed by two indented processes: the main process a the sifting process.The main process is performed via the following steps, which are plied to the signal s(t): 1. Executing the sifting process on the s(t), obtaining an IMF denominated as dn(t).2. Calculating the residual r(t) = sn−1(t) − dn(t) (then, this residual is considered a n s(t)).
The sifting process consists of the following steps: 1. Extracting relative extremes, obtaining two functions: v1(t) is the interpolation of maxima point, and v2(t) is the interpolation of the minima point 2. Determining the m(t) function as m The sifting process is iterated i-times, until the number of extreme values equals number of zero-crossings.The main process is iterated N-times, until r(t) is a monoto function or r(t) is smaller than a predefined threshold.
Therefore, the original signal is decomposed as: By means of this method, each data set, however complex, can be decomposed in finite number of components, each of which meets the definition of IMF.In Figure 3,

Hilbert-Huang Transform Approach
HHT is a way to analyse the signal through its Intrinsic Mode Functions (IMFs).The IMFs run in a simple oscillatory mode with varying amplitude and frequency along a time axis [35,36], and are used to obtain instantaneous frequencies at each moment in time with a symmetrical wave profile.The decomposition of the signal into its IMFs is obtained using the Empirical Mode Decomposition (EMD) method.
Empirical Mode Decomposition (EMD) [29,36] is an empirical method that exploits the definition of IMF to carry out the decomposition: a function with an equal number (or with a difference of 1) of extremes and zero crossing, whose envelopes are symmetric with respect to zero (i.e., at any point the average of the signal envelope determined by the local maxima and the envelope determined by the local minima is zero).It represents a generally simple oscillatory mode.
So, the EMD algorithm is composed by two indented processes: the main process and the sifting process.The main process is performed via the following steps, which are applied to the signal s(t): 1.
Executing the sifting process on the s(t), obtaining an IMF denominated as d n (t).

2.
Calculating the residual r(t) = s n−1 (t) − d n (t) (then, this residual is considered a new s(t)).
The sifting process consists of the following steps: 1.
Extracting relative extremes, obtaining two functions: v 1 (t) is the interpolation of the maxima point, and v 2 (t) is the interpolation of the minima point 2.
Determining the m(t) function as m(t The sifting process is iterated i-times, until the number of extreme values equals the number of zero-crossings.The main process is iterated N-times, until r(t) is a monotonic function or r(t) is smaller than a predefined threshold.
Therefore, the original signal is decomposed as: By means of this method, each data set, however complex, can be decomposed into a finite number of components, each of which meets the definition of IMF.In Figure 3, the developed EMD algorithm is schemed.All algorithms have been implemented in Python language (Python 3.11.4).The EMD method operates in the time domain, and it is highly efficient and very adaptive [37].It furnishes a holomorphic function that is suitable to be used with the wellknown Hilbert Transform (HT) to produce an accurate description of the signal.Hilbert Spectral Analysis (HSA) is then applicable to examine each IMFs as functions of time.The result is a frequency-time distribution of signal amplitude, named the Hilbert Spectrum, permitting the analysis of localized features [38].The Hibert Spectrum thus provides an accurate resolution in punctual time, instead of averaging over the entire dataset of a timewindow.

Time Domain Analysis
In Figure 4, the time series of the electrochemical noise potential and current, after a polynomial de-trend, are shown (Figure 4a,b, respectively).Polynomial de-trending is a type of filtering method that is widely used to remove trends in ENC and ENP signals [39].It entails approximating the trend using a polynomial of a given order m.The coefficients of the polynomial were identified using the least squares regression method.An m = 5 polynomial order was chosen.According to [40], this choice is a good computational compromise, since the change of noise resistance in m > 4 can be considered negligible, implying that a high m order of the polynomial is not required.Furthermore, m > 6 values are advised to avoid large computed mistakes and relevant variations at the start and end of the data.Moreover, with high m, Rn exhibits extremely irregular variations after drift removal.
EN potential, after polynomial de-trend, shows a low-frequency fluctuation at increasing time, with amplitude ranging from 10 mV to 20 mV.The potential oscillations are relatively large and wide with usually about 4000 s periodicity (point A in Figure 4a), although a constancy of this parameter is not clearly identifiable.Occasionally, mainly at longer times, some large fluctuations can be also observed (point B in Figure 4a).The EMD method operates in the time domain, and it is highly efficient and very adaptive [37].It furnishes a holomorphic function that is suitable to be used with the wellknown Hilbert Transform (HT) to produce an accurate description of the signal.Hilbert Spectral Analysis (HSA) is then applicable to examine each IMFs as functions of time.The result is a frequency-time distribution of signal amplitude, named the Hilbert Spectrum, permitting the analysis of localized features [38].The Hibert Spectrum thus provides an accurate resolution in punctual time, instead of averaging over the entire dataset of a time-window.

Time Domain Analysis
In Figure 4, the time series of the electrochemical noise potential and current, after a polynomial de-trend, are shown (Figure 4a,b, respectively).Polynomial de-trending is a type of filtering method that is widely used to remove trends in ENC and ENP signals [39].It entails approximating the trend using a polynomial of a given order m.The coefficients of the polynomial were identified using the least squares regression method.An m = 5 polynomial order was chosen.According to [40], this choice is a good computational compromise, since the change of noise resistance in m > 4 can be considered negligible, implying that a high m order of the polynomial is not required.Furthermore, m > 6 values are advised to avoid large computed mistakes and relevant variations at the start and end of the data.Moreover, with high m, Rn exhibits extremely irregular variations after drift removal.Current noise signal, after polynomial de-trend, is described by low-frequency broad fluctuations overlaid to high-frequency abrupt fluctuation with quick decay.At a low testing time, a series of small transient current peaks (repeated periodically in the whole transient) are observed (point C in Figure 4b); the majority of them are defined by a wide increase followed by a sharp reduction, which is said to be particularly prevalent during metastable pitting of stainless steels alloys [41][42][43] because of the slow growth of the pit, and prompt depassivation.
After about 1500 s, new transient phenomena, characterized by a sharp rise and a relatively short decay, were recorded (point D in Figure 4b).These events showed very high currents at, approximately, every 10-20 s, which suggests they were linked to the damage process of pits turning into cracks [44].
The third stage, at a longer testing time, is distinguished by a substantially large transient with a very sluggish decline.(point E in Figure 4b).
Considering the long duration of the acquired data and the significant fluctuations in frequency and amplitude ranges both on voltage and current transients, a visual timedomain analysis could not achieve a simple discrimination of several corrosion stages that take place during a SCC phenomenon [45].
However, in order to be able to acquire information related to the mode and stages of damaging mechanisms on SCC tests, electrochemical and statistical parameters from EN potential, after polynomial de-trend, shows a low-frequency fluctuation at increasing time, with amplitude ranging from 10 mV to 20 mV.The potential oscillations are relatively large and wide with usually about 4000 s periodicity (point A in Figure 4a), although a constancy of this parameter is not clearly identifiable.Occasionally, mainly at longer times, some large fluctuations can be also observed (point B in Figure 4a).
Current noise signal, after polynomial de-trend, is described by low-frequency broad fluctuations overlaid to high-frequency abrupt fluctuation with quick decay.At a low testing time, a series of small transient current peaks (repeated periodically in the whole transient) are observed (point C in Figure 4b); the majority of them are defined by a wide increase followed by a sharp reduction, which is said to be particularly prevalent during metastable pitting of stainless steels alloys [41][42][43] because of the slow growth of the pit, and prompt depassivation.
After about 1500 s, new transient phenomena, characterized by a sharp rise and a relatively short decay, were recorded (point D in Figure 4b).These events showed very high currents at, approximately, every 10-20 s, which suggests they were linked to the damage process of pits turning into cracks [44].
The third stage, at a longer testing time, is distinguished by a substantially large transient with a very sluggish decline.(point E in Figure 4b).
Considering the long duration of the acquired data and the significant fluctuations in frequency and amplitude ranges both on voltage and current transients, a visual timedomain analysis could not achieve a simple discrimination of several corrosion stages that take place during a SCC phenomenon [45].
However, in order to be able to acquire information related to the mode and stages of damaging mechanisms on SCC tests, electrochemical and statistical parameters from electrochemical noise data can be extrapolated.Several variables can be used to make a description of the corrosion phenomena, some of which derive from a time domain analysis, such as the standard deviation of noise potential and current, noise resistance, localization index, etc [46].
Noise Resistance, R N , is described as the ratio of the standard deviation of potential and current variation.
The Localization Index, LI, [47,48] is a variable that is related to the type of corrosion taking place: generalized (value close to 0) or localized (value close to 1).It has been demonstrated that this variable is very sensitive to electrode asymmetry [49].It is calculated as: where I noise corresponds to the electrochemical noise current signal, I to the mean value, and σ I to the standard deviation of the I noise signal.
Another approach aimed to analyse EN signal, based on the shot noise theory, is based on the concept that current noise signals may be thought of as charge packets (with certain limitations [50]).Using the equations that come next, the charge associated with every electrochemical event, q, and the event frequency, f n , may be calculated from the potential and current noise signals, according to [51]: where ψ E and ψ I are the low-frequency PSD values of the potential and current noise, respectively, B is the Stern-Geary coefficient (0.026 V [50]), and A the exposed electrode area (24.1 cm 2 ).The charge q gives an indication of the lost metal mass in the event, while f n provides information about the frequency of these events [51].The charge, q, is considered to be an affordable index of the localized corrosion activity [48,52], so the cumulative sum of the charge could be used as an index to evaluate the overall rate of electrochemical noise activities that take place on the specimen surface.
In order to use the previously mentioned set of variables, a sliding window method was used.Each window has a width of 10 min and an overlap of 60 s.Window dimension and overlapping were chosen in such a way to ensure a sufficient resolution to the variable q.
The cumulative charge-time evolution is reported in Figure 5.By analysing this trend, it is possible to detect a first region, up to 1500 s of stability, where no significant electrochemical activity was observed.This stage could be related with the pre-activation corrosion phase.After that, the cumulative charge increases, showing a progressive enhancement of the electrochemical activity on the specimen surface.Figure 6 confirms this by showing the evolution during the time of the current standard deviation.The spikes represent an increased variability in current signal, related with electrochemical event occurrence.The first significant peak, point A in Figure 6, occurs around 1000 s.This local current peak could be ascribed to the activation of metastable pits.Therefore, this phase can be related to electrochemical phenomena activation on the specimen surface.Therefore, the relevant step in the cumulative charge rate observed at about 1500 s may be associated with the activation and propagation of surface pits (point B in Figure 5).plot (point B' in Figure 5).In this stage, in the range 1500-10,000s, a relevant variability of the current standard deviation takes place, with intermediate peaks and valley regions (point C Figure 6).After approximately 10,000 s, a progressive rise of the cumulative charge trend takes place until sample failure occurs.Continuous fluctuations in the standard deviation trend with relevant peak values occur (points C in Figure 6), and are related to short to long range cracks propagation [53].In this stage, analysing the cumulative charge trend, continuous steps of increase and stabilization of the electrochemical activity can be identified (points C in Figure 5).Compatible with that in [14], this behaviour could be associated with the SCC cracks propagation.Due to the increase in crack length, the stress concentration increases at the crack tip, causing it to undergo more widespread blunting.A bigger plastic area means a longer period for the fracture to be reshaped by dissolution in order for a new crack to propagate [54].Hence, crack stabilization periods, where metal dissolution (and, therefore, electrochemical activity) is relevant, are coupled cyclically with crack evolution (related to low electrochemical activity) [55].
With the advancing of crack length, the steps, point B and B' in Figure 5, increase in height as a consequence of the progressive transition from short to long range crack propagation (magnitude from 0.6 nC to 2.7 nC). Afterward, for times above 40,000-45,000 s, the cumulative charge increases monotonically with short deviations from linearity (points C' in Figure 5).This period is the prelude of the final rupture of the sample that takes place after about 58,000 s.After approximately 10,000 s, a progressive rise of the cumulative charge trend takes place until sample failure occurs.Continuous fluctuations in the standard deviation trend with relevant peak values occur (points C in Figure 6), and are related to short to long range cracks propagation [53].In this stage, analysing the cumulative charge trend, continuous steps of increase and stabilization of the electrochemical activity can be identified (points C in Figure 5).Compatible with that in [14], this behaviour could be associated with the SCC cracks propagation.Due to the increase in crack length, the stress concentration increases at the crack tip, causing it to undergo more widespread blunting.A bigger plastic area means a longer period for the fracture to be reshaped by dissolution in order for a new crack to propagate [54].Hence, crack stabilization periods, where metal dissolution (and, therefore, electrochemical activity) is relevant, are coupled cyclically with crack evolution (related to low electrochemical activity) [55].
With the advancing of crack length, the steps, point B and B' in Figure 5, increase in height as a consequence of the progressive transition from short to long range crack propagation (magnitude from 0.6 nC to 2.7 nC). Afterward, for times above 40,000-45,000 s, the cumulative charge increases monotonically with short deviations from linearity (points C' in Figure 5).This period is the prelude of the final rupture of the sample that takes place after about 58,000 s.Then, at an increased time, a reduction of the electrochemical activity occurred, leading to the observation of a plateau with only some smooth steps in the cumulative charge plot (point B' in Figure 5).In this stage, in the range 1500-10,000s, a relevant variability of the current standard deviation takes place, with intermediate peaks and valley regions (point C Figure 6).
After approximately 10,000 s, a progressive rise of the cumulative charge trend takes place until sample failure occurs.Continuous fluctuations in the standard deviation trend with relevant peak values occur (points C in Figure 6), and are related to short to long range cracks propagation [53].In this stage, analysing the cumulative charge trend, continuous steps of increase and stabilization of the electrochemical activity can be identified (points C in Figure 5).Compatible with that in [14], this behaviour could be associated with the SCC cracks propagation.Due to the increase in crack length, the stress concentration increases at the crack tip, causing it to undergo more widespread blunting.A bigger plastic area means a longer period for the fracture to be reshaped by dissolution in order for a new crack to propagate [54].Hence, crack stabilization periods, where metal dissolution (and, therefore, electrochemical activity) is relevant, are coupled cyclically with crack evolution (related to low electrochemical activity) [55].
With the advancing of crack length, the steps, point B and B' in Figure 5, increase in height as a consequence of the progressive transition from short to long range crack propagation (magnitude from 0.6 nC to 2.7 nC). Afterward, for times above 40,000-45,000 s, the cumulative charge increases monotonically with short deviations from linearity (points C' in Figure 5).This period is the prelude of the final rupture of the sample that takes place after about 58,000 s.
Concerning this, different EN features (e.g., noise resistance, localization index, etc.), are usually applied as key parameters to detect localized corrosion [48,56].This is also applicable to a complex corrosive mechanism such as SCC [18].Details of the main statistical parameter of the ENC and ENP time series investigated in this paper are reported in [14]; however, usually it is required to take a multiple variable approach on EN features, e.g., by using a neural network method [14], to discriminate the fundamental parameters and their correlation during the different and synergistic phases of the process.Surely, the integration of all this information can represent a profitable approach to investigate the phenomena which evolves during the SCC.However, this requires support validation of frequency or time-frequency domain analysis.

Frequency Domain Analysis
With the purpose of better highlighting the differences in the EN signals during the SCC damage steps, frequency domain analysis was performed.
Due to the application hypotheses at the basis of the Fourier Transform, for a correct use of the FFT algorithm, a signal framing was performed.The size of the signal framing windows was chosen in order to have signal portions that are stationary and linear.This approach has the benefit to obtain combined time and frequency information, i.e., a nonnative representation for the Fourier Transform.
Due to the non-deterministic nature of the EN signals, based on FFT spectra, the PSD graph results in a more suitable analysis.This approach uses the auto-correlation function of the recorded signal, so it is capable to extract information about the signal power, but, at the same time, loses all the information based on the phase.In Figure 7, the time evolution of PSD of the current signal is reported.As shown in the graph, different phases, each of which could be related to the corrosion phases, are difficult to detect.Several fluctuations in the initial and medium period of the test can be identified.Therefore, it is possible to perceive different activities which take place, but it is difficult to discriminate and relate each different stage to a specific corrosion phenomenon.This is due to the multiple variation in PSD parameters (roll-off slope, white noise, cut-off frequency, stop frequency, magnitude), which simultaneously occurs during time.Therefore, this discrimination activity is very difficult to be argued against via the conventional amplitude graph, due to the multivariate nature of the corrosion phenomena.

Time-Frequency Domain Analysis
The Wavelet transform has a strong theoretical base as in the case of FFT/PSD theory, but it adds native support to the non-stationary signals.In addition, it furnishes a native energy-time-frequency representation.

Time-Frequency Domain Analysis
The Wavelet transform has a strong theoretical base as in the case of FFT/PSD theory, but it adds native support to the non-stationary signals.In addition, it furnishes a native energy-time-frequency representation.
The mother wavelet function, applied here, is a Ricker wavelet (also known as the Mexican hat wavelet) [57], with a scale depth of 30.The scale parameter is inversely proportional to the frequency parameter.
Figure 8 shows the Continuous Wavelet Transform applied to the EN current signal.A mixed time scale was performed.In particular, in order to better evidence all corrosion evolution stages during the SCC test, a log-time scale at a short time, and a linear scale at a longer time were utilised.According to the stages detected by the univariate analysis, the graph was divided in three stages.At first, up to about 1500 s, a growing electrochemical activity can be noticed.This area of activity shows a peak at the medium frequency range, and, as time increases, a widening of the peak towards the low frequencies is shown.A second phase, in the range 1500-5000 s, is characterized by the presence of competitive mechanisms, respectively, at medium-high and medium-low frequencies, followed by, at around 5000-10,000 s, an electrochemical quiescence zone.This region can be related to the triggering and mechanical propagation of the short-range cracks [14].
Afterwards, starting from 10,000 s and coded as the SCC short to long range crack propagation phase, according to Figure 5, more relevant electrochemical activities take place.This region is characterized by several subclusters of current transient signals, which are depicted by high magnitude and a scale frequency of about 10 3 Hz.These regions identify that, at this stage, the propagation phenomenon by SCC is due to electrochemical contribution.Furthermore, the continuous alternation between high-intensity peaks and areas of lower relevance, compatible with the presence of transition mechanisms between stabilization due to local plastic deformation and long range crack propagation, are slightly visible.This last phase turns out to be preparatory to the final failure.Consequently, although this technique is able to discriminate the pit triggered SCC phase with relative sub-steps, it does not allow effective discrimination of the long-range crack propagation stages by providing partial information on the damage evolution in this phase.
A further improvement could be provided by the fractional energy of smooth and Similarly to PSD analysis, the graph analysis shows the presence of different signal stages (detected by different magnitude at different scale).This analysis is advantageous, and enables us to obtain a better time resolution due to the absence of the approximation induced by the critical choice of the temporal window size which afflicts the STFT algorithm.
The Wavelet Transform cannot overcome the limitations related to the non-linearity of the signal.In addition, it has some inaccuracies related to the frequency resolution.Moreover, the choice of the mother wavelet turns out to be very critical.Despite these limitations, it is possible to make an interpretation of the electrochemical phenomena in progress during the SCC test.
According to the stages detected by the univariate analysis, the graph was divided in three stages.At first, up to about 1500 s, a growing electrochemical activity can be noticed.This area of activity shows a peak at the medium frequency range, and, as time increases, a widening of the peak towards the low frequencies is shown.A second phase, in the range 1500-5000 s, is characterized by the presence of competitive mechanisms, respectively, at medium-high and medium-low frequencies, followed by, at around 5000-10,000 s, an electrochemical quiescence zone.This region can be related to the triggering and mechanical propagation of the short-range cracks [14].
Afterwards, starting from 10,000 s and coded as the SCC short to long range crack propagation phase, according to Figure 5, more relevant electrochemical activities take place.This region is characterized by several subclusters of current transient signals, which are depicted by high magnitude and a scale frequency of about 10 3 Hz.These regions identify that, at this stage, the propagation phenomenon by SCC is due to electrochemical contribution.Furthermore, the continuous alternation between high-intensity peaks and areas of lower relevance, compatible with the presence of transition mechanisms between stabilization due to local plastic deformation and long range crack propagation, are slightly visible.This last phase turns out to be preparatory to the final failure.Consequently, although this technique is able to discriminate the pit triggered SCC phase with relative sub-steps, it does not allow effective discrimination of the long-range crack propagation stages by providing partial information on the damage evolution in this phase.
A further improvement could be provided by the fractional energy of smooth and low-frequency crystals in wavelet analysis to better identify the corrosion mode (general, metastable or stable pitting) and its start time in the sample during the corrosion phenomenon [58].By comparing the energy ratio of each crystal in the data sets, it could be possible, based on the power of WA, to distinguish corrosion types and assess their severity [59].
The Hilbert-Huang Transform (HHT), performed via the empirical mode decomposition (EMD) approach, was proposed to investigate the damage corrosion steps from the acquired EN non-stationary and non-periodic data.The EMD algorithm decomposes the original signal in its Intrinsic Mode Functions (IMFs): functions that respect the application condition of the Hilbert Transform.
The crack advancement rate, according to the SCC phenomenon, should be constant during the subcritical crack advancement phase.Certainly, given the interdependence of stress corrosion cracking growth rates, as determined by an integration of crack tip mechanics and crack tip oxidation kinetics [60], a suitable forecast for the influence of K on stress corrosion cracking growth rate is possible on these phenomena, and therefore could provide an improvement in knowledge regarding this concern.At first, the current approach is to evaluate, based on the ENC and ENP transient, how the electrochemical events are modified and how using a time-frequency technique (such as HHT), it is possible to identify macro-features of the process that are able to better discriminate, compared to conventional techniques, the phenomena for their correlation with the damage evolution steps in SCC.
In Figure 9, the results of the EMD decomposition for the current noise signal are shown.Only the first eight IMFs are presented.The current electrochemical noise was used for the HHT analysis, due to its powerful representation of corrosion phenomena compared to the potential signal.
The EN current signal has been decomposed into eleven IMFs.The IMFs that are characterized by a non-flat phase and an amplitude with a constant mean value are the most significant ones.IMFs that are characterized by a higher amplitude (and consequently greater power amount) includes large energy information of the source signal.A large number of IMFs (from the 6th to the 11th ones) are characterized by a very low amplitude and frequency (very close to the DC component) spectrum.Thus, these IMFs can probably be related to slow fluctuation phenomena (e.g., temperature, solution conductivity, etc).Consequently, only the first five IMFs can be considered relevant to the study of the corrosion phenomenon.
In Figure 10, the first IMF of the electrochemical noise current data is represented.The first IMF represents the most significant contributor to the whole signal and to the Hilbert Spectrum analysis, either in amplitude and in frequency content contribution.
sible to identify macro-features of the process that are able to better discriminate, compared to conventional techniques, the phenomena for their correlation with the damage evolution steps in SCC.
In Figure 9, the results of the EMD decomposition for the current noise signal are shown.Only the first eight IMFs are presented.The current electrochemical noise was used for the HHT analysis, due to its powerful representation of corrosion phenomena compared to the potential signal.The EN current signal has been decomposed into eleven IMFs.The IMFs that are characterized by a non-flat phase and an amplitude with a constant mean value are the most significant ones.IMFs that are characterized by a higher amplitude (and consequently greater power amount) includes large energy information of the source signal.A large number of IMFs (from the 6th to the 11th ones) are characterized by a very low amplitude and frequency (very close to the DC component) spectrum.Thus, these IMFs can probably be related to slow fluctuation phenomena (e.g., temperature, solution In Figure 10, the first IMF of the electrochemical noise current data is represented.The first IMF represents the most significant contributor to the whole signal and to the Hilbert Spectrum analysis, either in amplitude and in frequency content contribution.
It is possible to identify regions in which the activity is remarkably clear, interspersed with other regions where the quiescence of the electrochemical activity is noted.In particular, identifying the stages in which the amplitude magnitude of the first current IMF has major peaks, four time steps can be identified: 1.At first, up to 1500 s, a low electrochemical activity region can be identified.As previously discussed, this stage could be related with the pre-activation corrosion phase.At very short times the electrochemical activity is not representative.However, with It is possible to identify regions in which the activity is remarkably clear, interspersed with other regions where the quiescence of the electrochemical activity is noted.
In particular, identifying the stages in which the amplitude magnitude of the first current IMF has major peaks, four time steps can be identified: At first, up to 1500 s, a low electrochemical activity region can be identified.As previously discussed, this stage could be related with the pre-activation corrosion phase.At very short times the electrochemical activity is not representative.However, with increasing time, the intensity of the IMF progressively increases.The first intensity and wide peak, point A in Figure 10, occurs at about 1000 s, in correspondence with the triggering of metastable pits [61].Consequently, this phase can be associated with the superficial electrochemical activation of the specimen [62].

2.
The pre-activation stage continues to evolve during time, with the presence of several sub-peaks (point B in Figure 10), up to about 2000 s, in correspondence with the activation and propagation of superficial pits.Subsequently the IMF signal undergoes a sharp reduction in intensity, during which there are only occasional signal spikes (point B' in Figure 10).During this time, in the time interval 2000-10,000 s, predominantly, the pit-to-crack transition takes place.The triggering of the SCC phenomenon obviously stimulates the defect evolution partly by mechanical contribution making therefore the electrochemical one less significant.The occasional secondary IMF peaks, identified by point B', may be associated with the local formation and propagation of additional surface pits that evolve simultaneously during the stress corrosion triggering phase.

3.
After about 10,000 s a third stage was defined as short to long range crack propagation can be identified.During this stage several fluctuations in IMS signal can be highlighted.In particular After a long electrochemical quiescence step (point C in Figure 10) a significant peak can be identified at about 16,000-18,000 s (point C' in Figure 10).This trend related to stabilization and abrupt increase steps in IMF noise current signal is periodically identifiable.This behaviour is caused by the crack evolution by SCC [36].At this stage, the size of the crack can reach a length sufficient to assume that the stress concentration at the crack tip becomes significant, thus implying that there is a large area of plastic deformation, leading to a more extensive blunt at the crack tip.Due to the crack reshape, the mechanical crack evolution mechanism is inhibited, despite electrochemical dissolution [48].Consequently, a phase of electrochemical stabilization associated with the mechanical evolution of SCC damage is complementary to a phase of metal dissolution at the crack tip, and, therefore, complementary to an electrochemical activity phase.Depending on various electrochemical dissolution factors, the current transient generated by the micro-crack opening and re-passivation can last for a few seconds [63].These sub-steps alternate cyclically during this stage (points C* in Figure 10), where the sub-critical propagation of the medium and long-range crack takes place.4.
After about 45,000 s, after a transient region (coded as pre-quiescence), a large electrochemical plateau in the IMF signal can be identified (above 48,000 s).During this step, no relevant electrochemical vents can be identified.This stage is the prelude for the critical failure of the sample that took place after about 58,000s.
In order to be able to better discriminate the different time steps with significant electrochemical activity, it may be desirable to relate each time step with both amplitude and frequency of each IMF.
In Figure 11, the Hilbert transform of the entire electrochemical noise current signal is shown.This graph represents the trend of the frequency content vs. time.
is shown.This graph represents the trend of the frequency content vs. time.
The Hilbert Spectra allow to discriminate events occurring in different time instants with high precision, as well as facilitating the identification and characterization of these phenomena.Figure 11 can be used as a topological map by highlighting different clustered areas of the map, and relating them to specific degradation phenomena that take place during the SCC test.Each cluster in the map is defined as a homogeneous area characterized by similar frequency contents.Based on these considerations, four main cluster areas can be identified: The Hilbert Spectra allow to discriminate events occurring in different time instants with high precision, as well as facilitating the identification and characterization of these phenomena.Figure 11 can be used as a topological map by highlighting different clustered areas of the map, and relating them to specific degradation phenomena that take place during the SCC test.Each cluster in the map is defined as a homogeneous area characterized by similar frequency contents.
Based on these considerations, four main cluster areas can be identified: 1. Stabilization: at first, in the rage 0-100 s (step I in Figure 11), the noise signal is not characterised by any significant fluctuations.EN transient related to this phenomenon is characterised by the absence of high-frequency events: only some events at about 10 −3 Hz can be observed.In this phase, the interaction between the electrolyte and specimen surface takes place, with an unstable alternation of general corrosion and passivation phases.The aggressive ions locally destabilize the passive oxide layer of the sample, causing a local thinning of the oxide.This process has a short duration (approximately 100 s) due to the high aggressiveness of the environment conditions of the test.

2.
Electrochemical activation: in the range 100-1500 s (step II in Figure 11), in correspondence with depassivated surface regions, the electrochemical activity is enhanced and the pit initiation stage extends to micrometre-sized, large dissolved holes.In this phase, an increase in the cumulative charge trend was also identified.This substage is identified by a high activity at medium and low frequencies.In particular, at increasing time (after about 500 s) and with a progressive increment in the number of pits on the surface, the signal activity evolves towards a higher frequency (from 2 × 10 −2 to 5 × 10 −1 Hz, defining two substep: IIa and Iib).This stage can be related with the electrochemical activation phase (Iia) and the subsequent triggering of metastable pits (Iib) on the sample surface [14].

3.
SCC Activation and Propagation: This region, ranging from 1500 to 10,000 s, can be related mainly with the SSC triggering by short range crack activation.At this stage, the electrochemical activity becomes significant.Some sub-clusters can be identified.
At first, region IIIa, as shown in Figure 11, suggests that high-frequency contribution can be identified.At this stage, the aggressive ions locally destabilize the passive oxide layer of the sample, causing a local thinning of the oxide.In correspondence with depassivated surface regions, the electrochemical activity is enhanced and the pit initiation stage extends to micrometre-sized, large dissolved holes (see pit in Figure 12d).This process leads to the formation of preferential areas of localized attacks by pitting, identifiable by a high magnitude of the IMF.Progressively, the signal activity evolves towards a lower frequency (from 5 × 10 −3 to 5 × 10 −2 Hz) and magnitude (from 10 0.4 to 10 −0.4 dB), as shown via region IIIb in Figure 11.In this low frequency phase, the mechanical contribution is more relevant than the electrochemical one, as evidenced by the low EN signal magnitude, indicating that SCC short-range activation of the cracks was triggered [14].4.
SCC propagation.At about 10,000 s, a new step begins.During this step, the SCC damage evolution phase of the sample occurs (in Figure 12c, it is possible to note a secondary crack originating from a pit).The progressive crack growth (region III) leads to an increase in the region of plastic deformation at the crack tip [64].This implies that a more severe anodic dissolution is needed by the crack to re-sharpen before inducing a further propagation.At this phase, in fact, the dissolution within the crack is still taking place.When increasing the crack length, an increase in stress concentration occurs, which induces a larger plastic zone ahead of the crack tips.As a result, the crack point becomes more blunted.A larger plastic zone suggests that it will take more time for the fracture to be re-sharpened via dissolving for future crack propagation.This phase was indeed indicated as a mechanical quiescence [53] and the driving force in crack propagation.This stage shows metal dissolution at the crack tip.This region could correspond to the high-frequency activity in the first IMF.The electrochemical dissolution (sub-steps IVa in Figure 11), identifiable via the step-wise increase in cumulative charge (Figure 5), stimulates the mechanical instability of the crack (sub-steps IVb in Figure 11), primed for future subsonic propagation.Consequently, this stage identifies the prelude to the catastrophic failure of the specimen.In the topological map, this region is identified by medium-low amplitude valleys alternating with high-frequency peaks (5 × 10 −2 -3 × 10 −1 Hz).Afterwards, for times above 45,000 s, the Hilbert spectrum shows a low magnitude region where an electrochemical quiescence occurs.Finally, for catastrophic failure, the sample takes place after about 58,000s.
From a morphological point of view, the specimen's fracture surface was irregular, showing mixed ductile/brittle propagation modes (Figure 12a), and confirmed the previous argument in relation to analysing the data in Figure 11.Cracks began with a cleavage-like manner from the pit's bottom (point A, B and C in Figure 12a).A detail of a referencepitted cracked area is shown in Figure 12b, where the short-range crack propagation by cleavage mechanisms is identifiable.In Figure 12d, a surface illustration of a pit to crack progression is depicted (stage III in Figure 11).The corrosion is caused by random germinating open-hemispherical holes with diameters of around 20 µm on the surface of the samples, as well as the gradual initiation and propagation of river-shaped fractures with the pit as their origin.Moreover, for long-distance cracks, ductile crack propagation can occur in conjunction with the crack apex due to localized plastic deformation at the crack tip (dotted ellipse in Figure 12d).Figure 12c depicts further information on crack propagation mechanism in relation to the dog-bone gauge cross section.The crack started from superficial pit and then propagated in a branched trans-granular mode.
On the basis of the above reported considerations, it can be possible to divide the HHT spectrum into specific damage mechanism areas, according to the schematic representation shown in Figure 13, in which we can distinguish SCC short-and long-range crack evolution and failure stages.The possibility to have topological maps is a further advantage of the HHT method, providing a graphical support in order to better discriminate the significant evolving damage phases.
nating open-hemispherical holes with diameters of around 20 µm on the surface of the samples, as well as the gradual initiation and propagation of river-shaped fractures with the pit as their origin.Moreover, for long-distance cracks, ductile crack propagation can occur in conjunction with the crack apex due to localized plastic deformation at the crack tip (dotted ellipse in Figure 12d).Figure 12c depicts further information on crack propagation mechanism in relation to the dog-bone gauge cross section.The crack started from superficial pit and then propagated in a branched trans-granular mode.On the basis of the above reported considerations, it can be possible to divide the HHT spectrum into specific damage mechanism areas, according to the schematic representation shown in Figure 13, in which we can distinguish SCC short-and long-range crack evolution and failure stages.The possibility to have topological maps is a further advantage of the HHT method, providing a graphical support in order to better discriminate the significant evolving damage phases.Concerning the SCC short-range evolution stage, a crack activation cluster can be identified.This cluster is characterized by sharp peaks due to local corrosion attacks on the sample surface (pitting).Instead, the low-frequency region is related to short-range mechanical crack growth.Analogously, several sub-steps can be identified via analysing the SCC short-to-long-range evolution stage.EN events characterized by high magnitude and high frequency are clustered, and able to identify the dissolution of the crack tip's metal.All events at low frequency identify the transition of the crack growth from shortto-long-range propagation.This transition is also related to a modification from brittle (cleavage) to ductile in crack evolution (as valuable in Figure 12).This region is usually related with a mechanistic SCC quiescence phase [65].Afterwards, at a longer time, the electrochemical quiescence phase can be considered as a precondition for the final catastrophic failure of the sample.In conclusion, it is the author's opinion that HHT algorithms and the related Hilber-Huang Spectrum are powerful instruments for analysing non-stationary and non-linear signals of an EN data set, which overcomes the theoretical limits of Fast Fourier and Wavelet transforms.In these particular instances, the use of EN measurements not only vali- Concerning the SCC short-range evolution stage, a crack activation cluster can be identified.This cluster is characterized by sharp peaks due to local corrosion attacks on the sample surface (pitting).Instead, the low-frequency region is related to short-range mechanical crack growth.Analogously, several sub-steps can be identified via analysing the SCC short-to-long-range evolution stage.EN events characterized by high magnitude and high frequency are clustered, and able to identify the dissolution of the crack tip's metal.All events at low frequency identify the transition of the crack growth from short-to-long-range propagation.This transition is also related to a modification from brittle (cleavage) to ductile in crack evolution (as valuable in Figure 12).This region is usually related with a mechanistic SCC quiescence phase [65].Afterwards, at a longer time, the electrochemical quiescence phase can be considered as a precondition for the final catastrophic failure of the sample.
In conclusion, it is the author's opinion that HHT algorithms and the related Hilber-Huang Spectrum are powerful instruments for analysing non-stationary and non-linear signals of an EN data set, which overcomes the theoretical limits of Fast Fourier and Wavelet transforms.In these particular instances, the use of EN measurements not only validated earlier ideas established from time and frequency analysis, but also allowed the identification of multiple SCC damage processes happening over time on stainless steel samples.This separation was accomplished by the establishment of correlations between key EN properties and the relevant damage mechanisms.
Considering that, in SCC, the effects of the electrochemical reaction and mechanical stress occur simultaneously, a combined approach of the two techniques can be considered effective.In several papers [14,66,67], it was possible to detect how the combination of an electrochemical technique (EN) with a mechanistic one (e.g., acoustic emission, AE), by means of neural network analysis, can be profitable in order to determine the key features of the SCC process and the evolutionary states of damage.With this in mind, a future scope will be the possibility to synergistically combine the EN technique, analysed by means of HHT, with the acoustic emission technique by integrating them through neural network analysis.

Conclusions
The EN technique was used to investigate the SCC damage evolution on 17-4PH stainless steel samples.The conventional methods (time and frequency domain analysis) showed promising results on predicting the damage phase.However, they require a complex interpretation of the results, due to the system's multivariate nature.Furthermore, the theoretical assumptions involve some approximations when applied to the EN signals.
For this purpose, the HHT method, based on an empirical mode decomposition (EMD) approach, allows calculating and analysing the EN signal accurately.This ability to handle and process non-stationary and non-linear data results in a detailed representation of the signal contents in the time-frequency domain.
Using the technique, clear information was obtained about the transitions between different corrosion phases caused by the SCC test.It was possible (i) to identify a specific response to each corrosion damage phenomenon, (ii) determine the timeline with accuracy, and (iii) also provide a topological representation of the SCC damage stages as they evolved.
In particular, the following stages were detected: • Stabilization: at first, in the rage 0-100 s.The EN transient is characterised by the absence of high-frequency events.Only some events at about 10 −3 Hz were observed.

•
Electrochemical activation: in the range 100-1500 s.This stage was identified by a high activity at medium and low frequencies (from 2 × 10 −2 to 5 × 10 −1 Hz,).This stage can be related with electrochemical activation and the subsequent triggering of metastable pits on the sample surface.

•
SCC Activation and Propagation: ranging from 1500 to 10,000 s.This is mainly related to the SSC triggering by short-range crack activation.At this stage, the electrochemical activity became significantly identifiable by a high magnitude of IMF.Then, pro-gressively, the signal activity evolved towards a lower frequency (from 5 × 10 −3 to 5 × 10 −2 Hz) and magnitude (from 10 0.4 to 10 −0.4 dB), because the mechanical contribution was more relevant than the electrochemical.• SCC propagation.Above 10,000 s.This is related to progressive crack growth.This region was identified by medium-low amplitude valleys alternating with high-frequency peaks (5 × 10 −2 -3 × 10 −1 Hz).For times above 45,000 s a low magnitude region, related to an electrochemical quiescence, was identified before catastrophic failure of the samples.
These results indicate the HHT method as a potentially powerful technique for EN signal processing.Additionally, HHT is able to recognize and distinguish different damage mechanism stages that occur during degradation.

Figure 2 .
Figure 2. Scheme of EN damage analysis procedure.

Figure 2 .
Figure 2. Scheme of EN damage analysis procedure.

Figure 3 .
Figure 3. Schematic representation of EMD algorithm.# indicates the analysed numbers.

Figure 3 .
Figure 3. Schematic representation of EMD algorithm.# indicates the analysed numbers.

Figure 5 .
Figure 5. Normalized Charge Cumulative trend.A coupled log-linear time scale is performed, in order to evidence the short time segment.

Figure 5 .
Figure 5. Normalized Charge Cumulative trend.A coupled log-linear time scale is performed, in order to evidence the short time segment.

Figure 5 .
Figure 5. Normalized Charge Cumulative trend.A coupled log-linear time scale is performed, in order to evidence the short time segment.

Figure 6 .
Figure 6.Current Standard Deviation evolution during time.A coupled log-linear time scale is performed in order to evidence the short time segment.

Figure 7 .
Figure 7. Time evolution of the Power Spectral Density related to the current EN signal.The coloured bar indicates PSD amplitude [A 2 Hz −1 ].

Figure 7 .
Figure 7. Time evolution of the Power Spectral Density related to the current EN signal.The coloured bar indicates PSD amplitude [A 2 Hz −1 ].

12 Figure 8 .
Figure 8. Continuous Wavelet Transform applied to the EC current signal.A non-linear time scale is performed in order to evidence the low time characteristics.A coupled log-linear time scale is performed in order to evidence the short time segment.

Figure 8 .
Figure 8. Continuous Wavelet Transform applied to the EC current signal.A non-linear time scale is performed in order to evidence the low time characteristics.A coupled log-linear time scale is performed in order to evidence the short time segment.

Figure 10 .
Figure 10.Hilbert-Huang Transform applied to the electrochemical noise current data: only the first IMF is shown.A coupled log-linear time scale is performed in order to evidence the short time segment.

Figure 10 .
Figure 10.Hilbert-Huang Transform applied to the electrochemical noise current data: only the first IMF is shown.A coupled log-linear time scale is performed in order to evidence the short time segment.

Figure 11 .
Figure 11.Hilbert Spectrum representation: up: 3D view, down: Frequency/Time contour plot.A coupled log-linear time scale is performed in order to evidence the short time components.

Figure 11 .
Figure 11.Hilbert Spectrum representation: up: 3D view, down: Frequency/Time contour plot.A coupled log-linear time scale is performed in order to evidence the short time components.

Figure 12 .
Figure 12.SEM-SE morphology of the fracture surface, details (a-c).Cracks began in a cleavage pattern from the pit's bottom.The specimen's final rupture was identified by ductile propagation.Furthermore, gauge cross section displaying a pit-originating fracture (Vilella's reagent etch) (d).

Figure 12 .
Figure 12.SEM-SE morphology of the fracture surface, details (a-c).Cracks began in a cleavage pattern from the pit's bottom.The specimen's final rupture was identified by ductile propagation.Furthermore, gauge cross section displaying a pit-originating fracture (Vilella's reagent etch) (d).

Figure 13 .
Figure 13.Damage steps in Hilbert Spectrum Map.

Figure 13 .
Figure 13.Damage steps in Hilbert Spectrum Map.

Table 1 .
PH stainless steel nominal chemical composition.