\section{Introduction}

Radio Frequency Interference (RFI) has become a significant factor in our ability to observe natural phenomena
on Earth, in our solar system, and the nearby and distant universe.
From the radio scientist's point of view, RFI is considered to be any unwanted addition to the wanted signal that
has the potential to degrade or prevent the successful conduct of an observation.
The chance for encountering RFI has been increasing steadily because of the ability to observe increasingly
weaker signals with more sensitive detection systems that operate over larger bandwidths.
In addition, because spectral lines are fixed fundamental (rest) frequencies, and radio astronomers are looking at weaker
emission lines from more distant sources in the universe that are redshifted outside science service bands into bands allocated
to active radio services.
And while the active spectral usage is intensifying and the used active service bands are moving to higher frequency,  the
chance for finding strong interfering signals in the science data is becoming larger.

 Unlike thermal noise, which has stable temporal stochastic properties (white noise) and can be dealt with through radiometric
 detection (i.e. longer integration times and on-source minus off-source subtraction), an RFI signal is temporally, spatially or
 spectrally structured and can obscure a wanted signal or produce a false positive detection.

The Radio Regulations issued by the Radiocommunication Sector of the International Telecommunication Union (ITU-R) assign Primary and secondary users to allocated spectral bands, such that secondary users are not allowed to interfere with the Primary use of the band \citep{ITU-R}. This regulatory system works well in general, except that some users do not respect these rules.
In 'passive' bands that are allocated to the science services are Primary and all emissions are prohibited.
However, because this requirement is 'too difficult to adhere to' for neighbouring active services, RFI may need to be tolerated in those bands for 2\% of time for a single network or 5\% for an aggregate situation.
In addition, one may encounter out-of-band and spurious emissions resulting from emissions from adjacent bands and there may also be some illegal transmissions.
In bands shared between science services and other services, it depends on the status of the allocated services how RFI issues need to be coordinated.
Outside the bands allocated to the science services, science users have no rights of protection and may operate in these bands on a 'non-interference' basis relative to the Primary users of the band.

Traditionally much of the data corrupted by RFI has been discarded. However, technology and computing advances have made it possible to mitigate - to lessen in intensity and damage - the effect of the RFI signals in scientific data.
Yet RFI mitigation has often been an afterthought and retro-active implementation of mitigation options at observatories has been very slow. Often the severity of the anticipated long-term RFI conditions has been underestimated by the users.
Furthermore, the complexity of integrating RFI mitigation systems into existing observing systems was too high because these systems were not yet fully digital.
Finally, the users preferred (old-fashioned) hands-on flagging methods over the adaptation of black-box procedures.

Changes in technology have made it possible to implement RFI mitigation in existing systems as well as new instrument projects.
Rather than simply complaining about the presence of RFI in the data and about the damage it does, it it possible to do something about the situation.
In this paper we discuss methods of mitigation that have been tried and successfully tested.

Previous reviews covering many of the methods and issues presented in this paper are found in \citep{BriggsK2005, Kesteven2005, BaanB2010,AnEA2018} and ITU-R Report 2126 \citep{Rep2126}

\section{Managing radio frequency interference}

RFI always degrades the quality of the data and the working assumption for most observations has been that RFI-corrupted
data are partially or completely unusable.
The most common method for dealing with RFI is to excise spectral or temporal segments of a data set that are known
to be corrupted.
There are powerful motivations to move beyond throwing away data. Methods that remove or mitigate RFI, thus enabling
scientific usage of data that would otherwise be discarded, are becoming more and more essential and feasible.
In addition, the dramatic increase in the size of data sets make automated procedures ever more necessary.

In practical terms and from the point of the science user, any external signal that disturbs the wanted scientific signals can be
considered as RFI, whether is it strong and showing in the instrument bandpass or weak and only showing in integrated data.
The ITU-R has defined the 'power flux density' (pfd) levels that are considered detrimental interference that are well within
the noise of the observing instruments.
For Radioastronomy they are given in Recommendation ITU-R RA.769 and are based on 10\% of the system noise of the instrument.
The percentage of permissible data loss in Passive bands resulting from emissions above these thresholds is specified in
Recommendation ITU-R RA.1513 \citep{ITU-RA}.
In the exclusive Primary bands for the Passive Radio Astronomy Service (RAS) and and the Earth Exploration Service (EESS) are listed in RR No. 5.340, all emissions are prohibited.
In the other radio astronomy bands, listed in RR No. 5.149, administrations are urged to take all practicable
steps to protect the radio astronomy service from harmful interference \citep{ITU-R}.
For Remote Sensing (RS) and the EESS the interference levels are discussed in
Recommendations ITU-R RS.1029 for passive sensors, ITU-R RS.1166 for active sensors and ITU-R RS.1263 for
MetAids \citep{ITU-RS}.
An important Resolution 750(WRC-15) presents compatibility  studies  between EESS(passive) and
relevant active services and discusses hard limits for OOB emissions into EESS bands \citep{Res750}.

The aim of mitigation efforts is to optimise the increased availability of signal processing hardware and algorithms
and to enable science observations to be conducted in densely occupied bands and heavily used radio environments.
As a result, RFI mitigation requires more complexity in order to deal with the growing number of wireless
communications services and the increased need for observations outside protected/allocated bands.

A guiding principle for the implementation of RFI Mitigation should be to remove interfering signals as early as
possible in the data chain and at the highest possible time and spectral resolution.
Since the removal of RFI signals is generally limited by the instantaneous signal-to-noise of the signal in the system,
the strongest and most damaging RFI signals may be mitigated at the sampling speed in the system.
Weaker signals that are still buried in the sampling noise can only be removed after time-integration later in the signal path.
Mitigation of the RFI as early as possible in the data chain minimises the damage to the data and
keeps the data loss at a minimum level.

Mitigation at high time resolution is particularly beneficial for time-variable RFI and will reduce data loss due to time-smearing.
Higher spectral resolution reduces the occupation of spectral channels by the RFI and limits the damage of RFI
to the data due to spectral-smearing.

In the following sections, methods of mitigation are discussed for a more general structure of receiving systems.
Most of this discussion is geared towards radio astronomy instruments but many methodologies will also be applicable
for EESS systems for both terrestrial and space applications. However, specific applications for space-based Remote Sensing applications will not be discussed in this review.

\begin{figure}
\begin{center}
\includegraphics[width= 10cm]{LayersOfMitigation.eps}
\end{center}
\caption{Stages and methods of RFI Mitigation along the signal path in a radio observing instrument.}
\label{fig:fig1}
\end{figure}

\section{Implementation in observing systems - a multi-layer strategy}

Effective mitigation of RFI signals needs to be multi-layered using different technical methods at different locations
in the signal path of the system and can be subdivided in three types of activity:
1) Preventing radio frequency interference (RFI) signals from entering the science data,
including both the reduction of the observatory's vulnerability to RFI signals and the elimination of in-house signals.
 2) Removing the strongest and/or time-variable RFI signals from the data in real time at one or more locations in the data chain.
 3) Off-line removal or reduction of the impact of RFI.

Multi-layer mitigation can be introduced at different stages in the data path and can be grouped into different
categories (Figure \ref{fig:fig1}).  Naturally, the ability to apply these mitigation protocols will depend on the structure of observing
system, and some methods cannot be applied to all systems.
Certain methods can be applied at and before each individual observing station, others in the signal chain of the stations before
station data processing or before array central processing, others apply to the individual/central processing site, and lastly
some methods may be applied after the data has been processed.
The ability of all methods to excise or subtract RFI signals from the incoming data is limited by the Interference to Noise Ratio (INR)
and they can reduce the RFI to a level of INR $\approx$1. Therefore strong RFI signals can be dealt with early in the data path,
while weaker RFI signals can only be addresses after integration.

The multi-layer mitigation scheme presented in Figure \ref{fig:fig1} includes an "adaptive control' component.
After a 'training process' this unit should be able to recognize the specific type of
interference encountered and adjust the system to invoke the optimal mitigation algorithm.
The long-term goal for this unit is to install an
autonomous RFI mitigation system that can handle all kinds of interference environments.
Another function to be integrated into this control unit is the accurate bookkeeping on 'what has been done where'.

\section{Station - Pro-active and pre-detection}

{\it Cleaning the Site} - Pro-active and pre-detection actions are aimed at preventing RFI signals from entering the
detector systems or reducing their strength.
In the first place, locally generated RFI from inside and outside the observatory needs to be eliminated or reduced to
levels that are below RA.769 values at the entry of the sensor or telescope.
This form of Electromagnetic Compatibility (EMS) should address all radio frequency noise generated by electrical and
electronic equipment and requires shielding/screening or placement in Faraday cages. This effort would also include laptops
and mobile communication devices in the building and on the site.
In the case of distributed buildings at an observatory site, also the radio noise generated in such buildings should be reduced
to RA.769 levels at the location of the telescope or sensor.
The observatory should not be its own (worst) enemy.

\vspace{2mm}{\it \bf Quiet and Coordination Zones} - Working close with governmental spectrum and regulatory agencies
it is necessary to establish a Quiet Zone (QZ) around the observatory.
It is generally not possible to make such a quiet zone actually radio quiet, but it will prevent new transmitters to be located
inside the zone and help remove some old ones.
In addition, a Coordination Zone (CZ) should be established that facilitates coordination with the operators about
nearby transmitters in order to reduce or eliminate the signal strength in the direction towards the observatory.
Many radio astronomy observatories have established quiet/coordination zones.
Large-scale zones have been documented for several sites, such as the Mid West Radio
Quiet Zone in Western Australia \citep{ARQZWA2007}, the National Radio Quiet Zone around Green Bank
\citep{NRQZ1958} and the Puerto Rico Coordination Zone around the Arecibo Observatory \citep{PRCZ1998}.
Information about the structure and characteristics have been described in ITU-R Report RA.2259.
Quiet zones for Earth exploration systems are realistic for terrestrial applications.
However, RS space applications will mostly depend on coordination with terrestrial interferers.

Transmitters existing inside a Quiet Zone at the time of the establishment of the zone or after need
to be identified and need to be quieted or relocated with the help of the national spectrum agency.
Friendly coordination with transmitting stations existing in a Coordination zone may result in lowering their power or
changing their emission pattern or change channel in order to reduce the signal strength at the observing frequencies
of the observatory. In general, the operators of existing stations are sensitive to the needs of an observatory.
New stations in a coordination zone need to be coordinated before the installation takes place.

\section{Station Mitigation at the detection system}
Mitigation at the detection system should be applied to standalone stations as well as to all stations that are part of an array.
The use of robust receivers for observations that can withstand the presence and effects of strong RFI is essential.
The suppression of second-order intermodulation products and enough 'headroom' to avoid gain suppression are
guiding principles for the construction of such receivers.

Spectral filtering and band-edge filtering have been a standard method to avoid strong signals entering the band.
However, filtering comes at a price in gain performance. For in-band filtering, super-conducting notch filters have been
developed. These costly devices will slightly degrade the performance of the receiver but their tuning is not very flexible.

Edge filtering or band selection to avoid band sections with strong signals may be applied for continuum studies, but
observers enjoy less flexibility with regard to spectral line studies and redshifted surveys. Application of software
defined radio techniques will aid in finding clean spectrum and avoid RFI signals.

\section{Station Mitigation before correlation and processing}

\subsection{Excision at baseband}

\vspace{2mm}{\it \bf Time and frequency analysis} -
A powerful way of removing strong RFI signals that have entered the receiver is real-time excision at the Intermediate
Frequency (IF) (at baseband) at the sampling speed.
Signals that have significant SNR at the sampling frequency are also most damaging to the surrounding spectral data and
should be removed from the data before further processing.
Successful applications of excising RFI from the data in the IF stage have been based on thresholding and kurtosis schemes.
Real-time thresholding simply applies a threshold of $n$ times the root mean square (rms) of the noise for removing IF signals
above this level in the frequency and/or time domains.
A successful implementation of this and other methods has been documented for the Westerbork Synthesis Radio Telescope \citep[WSRT;][]{BaanFM2004,BaanEA2010}.
Recently, a real-time application of the Median Absolute Deviation (MAD) estimator \citep{Fridman2008} for RFI thresholding
has been implemented (pre-correlation) within the GMRT wideband backend in order to target  time-variable power-line RFI
particularly at frequencies less than 700 MHz \citep{BuchEA2018}.

Similarly a dynamical excision threshold may be based on a {\it Kurtosis analysis} in order to identify and excise statistical outliers.
Kurtosis is a measure of the 'tailedness' of the probability distribution using a scaled version of the fourth moment of the data such
that a higher kurtosis number results from infrequent extreme deviations.
This statistical method has been applied successfully for removing signals with kurtosis $> 3$ from solar data \citep{NitaEA2007,GaryLN2010}.
Kurtosis excision at baseband is also applied to pulsar observations at the Parkes telescope \citep{Hobbs2018} and works
well for RFI signals that are significantly stronger than the pulsed signals of the pulsar.
Caution should be applied for pulsar searching observations because the pulses for these sources are still unknown.

In particular, time-variable interference can be addressed optimally by these methods, although they will result in a certain percentage of data loss that depends on the bandwidth and duty cycle of the RFI signal.
However, not removing these signals would result in even more damage and data loss after integration at a later stage.
As an example, the out-of-band emissions from the uplink-downlink of a Mobile Satellite System with a 50\% duty cycle would
result in data loss in the adjacent RAS 1612 MHz band in the range range from 50-100\%.
Real-time thresholding at baseband with the RFIMS observing system at Westerbork Observatory made
it possible to only sample the uplink time segment of Iridium with good results, but still with a data loss of about 50\% (Figure \ref{fig:fig2}).
For pulsar observations with the new Apertif {\it phased-array feed} (PAF) system at WSRT, thresholding in time and frequency
is applied on one second integrated samples after the beamforming process together with coincidence detection for all beams
\citep{vLeeuwen2018}.

\begin{figure*}
\begin{center}
\includegraphics[scale=0.45]{pulsar1.eps}
\includegraphics[scale=0.45]{pulsar2.eps}
\end{center}
\caption{WSRT observations of pulsar PSR\,0332+5434 under severe RFI conditions at 1625 MHz. (Left) A period folded
spectrum without RFI mitigation. (Right) The period folded spectrum with thresholding excision clearly shows the pulsar profile.
Figures are courtesy of Petr Fridman and Ben Stappers. }
\label{fig:fig2}
\end{figure*}

Thresholding and kurtosis analysis may be used to clip the power of strong RFI signals and replace these with some
smaller values. This 'data loss' changes the noise properties in the affected data and may also affect the data contained wanted
signals in (sometimes) unknown locations in the spectrum. Accurate bookkeeping of the 'when and where' is therefore an
important part of the process.

Implementation of mitigation and excision at baseband requires the (simple) inclusion of processing power in the IF path of the
instrument before further processing or correlation. For IF processing in PAFs the compute power required to process all elements
will be significant, and processing after beamforming can suffice if the RFI does not affect the beamforming process.

\vspace{2mm}{\it \bf Time Blanking} - In the presence of strong periodic signals such as radar systems, a blanking scheme
may be adopted where the backend and data taking is interrupted starting at the leading edge (in time) of the radar pulse.
The blanking window will have to take into account both the pulse and the multi-path signals caused by the intervening terrain.
This method has been applied at the Arecibo Observatory for certain radars at the San Juan (SJU)
airport, which allows observations inside the radar band with a certain data loss that would not have been possible otherwise.
Thresholding at baseband (discussed below) would achieve the same results as blanking the radar at the backend and with a similar amount of data loss.

\begin{figure}[h]
\begin{center}
\includegraphics[width= 8cm]{Fig7-FB-2001.eps}
\end{center}
\caption{Methods of RFI cancellation using autoclean filtering (frame a) and RFI excision by filtering
using a reference channel (frame b). Image obtained from Fridman \& Baan (2001).}
\label{fig:fig3}
\end{figure}

\subsection{Removing the RFI signature}

A number of {\it waveform subtraction} methods and methods that employ the statistical properties of the data have
been used to cancel and remove the added RFI power from the spectral data in the time domain with minimal or
no damage to the data itself.  However, the effectiveness and the achieved suppression are limited by the quality of
the estimate of the interference received by the instrument.

\vspace{2mm}{\it \bf Autoclean} - Temporally spread and strongly correlated RFI can be suppressed using cancellation techniques
based on estimating the RFI waveform and subtracting it from the {\it signal + RFI} mixture.
The RFI waveform can be estimated using available filtering technique (spine-smoothing, wavelet analysis, Wiener filtering,
parametric estimation). In this {\it autoclean procedure} the estimate is then subtracted from the input data in the temporal or
frequency domain (see Fig. \ref{fig:fig3}a). In principle, these methods do little or no damage to the data and the
signal of interest.

An application of the autoclean procedure has been described employing a WSRT telescope and comparing the result with
that of an 'uncleaned' antenna \citep{FridmanB2001}.
An example of a parametric approach to this type of RFI cancellation has been presented for the strong
interfering signal of a GLONASS satellite operating within the OH 1612 MHz band \citep{EllingsonEA2001}.
The RFI parameters Doppler frequency,
phase code and complex amplitude were determined for each data segment, to be used to calculate the RFI waveform for
subtraction from the {\it RFI + noise signal} mixture.
A recent parametric application is described for the removal of a broadband digital video broadcast (DVB) signal from the data \citep{SteebDW2018}.

\vspace{2mm}{\it \bf Adaptive Noise Cancellation and Reference antennas} -
A separate dedicated reference channel may be used to obtain an independent estimate of the RFI signal.
This technique of {\it adaptive noise cancellation} (ANC) is used actively in digital processing for communication and military
technology.
Two data channels as shown in the block diagram of Figure \ref{fig:fig3}b are a 'main channel' radio telescope and an
'auxiliary or reference channel' radio telescope both pointing 'on-source' and containing the RFI signal.
Since the RFI signal in both channels is not exactly identical because of different propagation paths and radio receivers,
an adaptive filter is used to reduce the error signal and apply the result to the main channel.
This basic principle of this procedure both in the temporal domain ({\it adaptive filtering}) and also the frequency domain
is to make an fast Fourier transform (FFT) from the incoming data, perform an adaptation operation on the frequency bins,
and then return to the
frequency domain via an inverse FFT. This method, based on Wiener filtering, works for interfering signals with a
significant INR, i.e. when the RFI dominates the system noise, and the suppression of the interfering signal can be
equal to its instantaneous INR. Examples of the use of reference channels are described in \citep{FridmanB2001}.

A variation on adaptive filtering is to subtract a reference data-channel from a signal data-channel by comparing {\it on-source
plus RFI} and {\it off-source plus RFI} signals. Separate reference antennas pointed at the interfering source have been
successfully used at the GBT and at Parkes \citep{BarnbaumB1998,BriggsBK2000,NigraEA2010}.

Adaptive filters are effective for spectral line observations, where the RFI and the signal-of-interest occupy the same
frequency domain, and when spectral information is unimportant, such as in pulsar \citep{Kesteven2005}
and continuum studies.
However, these methods should preferable be applied before data processing but may also be applied after depending
on the SNR of the external signal.
In principle, the signatures for multiple interfering sources may be added before subtraction from the wanted signal.

\vspace{2mm}{\it \bf Higher-order statistics and Probability Analysis} -
The voltages of the system noise and a radio source generally have a Gaussian signal distribution with a zero mean.
Fourier transformation of such an ideal signal also gives real and imaginary components in every spectral bin, that are Gaussian random values with zero means.
On the other hand, the instantaneous power spectrum (the square of the magnitude of the complex spectrum) has an exponential
distribution that is described as a chi-squared distribution with two degrees of freedom.
The presence of an RFI signal modifies the ideal input signal and yields a change of its statistics, giving a power spectrum with
 a non-central chi squared distribution with two degrees of freedom. Real-time analysis and DSP processing of
 the distribution before any averaging will allow a separation of the two signal components such as an RFI signal superposed
 on a known spectral line \citep{FridmanB2001}.

Higher-order statistics methods and determining the moments of the probability distribution of the signal power can be
used both for spectral line and continuum observations, and it will be useful to introduce this procedure into new detection
systems \citep{Fridman2001}.
However, reliable estimates of these higher moments (or cumulants) will require more averaging than for the first moment (the mean), which is necessary for mitigation of the weak RFI. On the other hand, large averaging intervals will
smooth the variability of the RFI and will yield estimates with a considerable bias.
Therefore, there are limitations on the detection and excision of weak RFI signals, which is equally true for all other RFI
mitigation methods.

\vspace{2mm}{\it \bf Multi-Element and Phased-Array Systems} -
In addition to station processing using thresholding or kurtosis techniques, adaptive filtering is particularly useful for
multi-feed single-dish radio telescopes and array instruments, where one array
element is being used as a reference channel. For an array instrument each strong and distinct source of RFI requires a
separate reference antenna.

{\it Spatial filtering} is widely used in {\it smart antennas} in (multi-element) radar and communication systems.
The applicability of spatial filtering for sparse arrays may be limited because of their sparseness and because of offline processing of
long integrations during which RFI sources may also be moving.
A variety of specific algorithms including maximum
SNR, subspace projection, Wiener filtering, and multiple sidelobe cancelling have been studied for application to radio astronomical observing \citep{Boonstra2005, BoonstraT2005,Ellingson2003,EllingsonH2002,JeffsEA2005,FridmanB2001,LeshemVB2000, vdTolV2005}.

Multi-layer real-time {\it adaptive/spatial-nulling} techniques may be applied for telescopes with phased-array technology
but also for densely populated sensor stations, such as for recent low-frequency systems or new generation instruments \citep{vArdenneSH2000,BentumEA2008,BlackEA2015,Boonstra2005,Bregman2000}.
The use of a reference antenna with a direct look at the interferer with a higher SNR will greatly improve the nulling performance \citep{BriggsBK2000,JeffsEA2005,SardarabadiEA2016}

In order to fully exploit these techniques, multi-element systems should have computer control of the antenna phases
and their amplitudes \citep{FridmanB2001}. Adaptive filtering using a beam-forming algorithm requires a high INR and
is limited to a small number of RFI targets to be tracked during an observation. The RFI sources also need to
remain stable and predictable through an observation. Spatial filtering in beam-forming mode for a
limited number of RFI sources generally does not degrade the image generated by the main beam.

\vspace{2mm}{\it \bf Subspace Projections for Multi-Element Systems} -
{\it Subspace projection} for array null-formation identifies the interference in terms of correlations between
array elements, which can be used to determine beamforming coefficients that result in patterns that reject the
interference with little or no effect on the main lobe characteristics.
Subspace projection using outstanding RFI properties has significant advantages for radio astronomy \citep{RazaEA2002}
but does not help with poor detection and localisation performance for the interference \cite{EllingsonH2002, LeshemVB2000}.

Beam pattern distortions when subspace nulling a (rapidly) moving source of interference and for narrowband signals
can be reduced with deeper nulls when time-integrated data is stored and processed  \citep{JeffsW2008, LandonJW2011}.
However, the cost of digital/transport infrastructure and a possible auxiliary antenna is an obstacle for implementation
and there are limits in nulling below the noise floor \citep{JeffsW2013}. There is also the need for a good science case.
In general, null-forming is most applicable to mitigation of RFI from satellites, and can be expected
to be somewhat less effective against terrestrial RFI with intervening terrain.

\section{Station and Central - Mitigation at correlation and post-correlation}

As part of the correlation process, digitised data are generally integrated over time intervals ranging
from the sampling time up to seconds, which significantly raises the INR. In consequence,
persistent but weak RFI, that could not be treated in real-time, and weak (spectral) remnants of
earlier mitigation operations become accessible for processing.
Therefore, a second layer of (automatic) real-time excision can be applied on the accumulated data
records that would be complementary to off-line flagging of the data.

Flagging and excising interference in the frequency or time domain has been the standard procedure used by
radio scientists during post-correlation processing.
Flagging of baselines and antennas in array observations also identifies and eliminates system problems.
Post-correlation excising is performed on integrated/averaged and correlated data, and can result in considerable
data loss because even for time- and frequency-variable RFI whole time-slots, whole baselines, and/or whole
antennas need to be flagged.
This differs from antenna-based flagging/excising of IF baseband data, which results in less data loss overall.

\vspace{2mm}{\it \bf On-line or off-line processing} - Anti-coincidence protocols may be incorporated at the processing
stage in order to identify
the RFI components, as well as digital mitigation processing and the integration of a reference antenna during
(software) correlation. However, the implementation of these algorithms into pre-existing hardware backends requires the
addition of both special hardware and software.

Automated flagging and excision of calibrated/processed and integrated data records has been proposed and implemented for
single-dish systems or for each baseline of an interferometry system
\citep{Kalberla2010,KeatingEA2010,Middelberg2006, OffringaEA2010, SirothiaEA2009}.
Recent developments on Recurrent Neural Network algorithms and Deep Learning systems may provide different options
for recognizing/identifying RFI in the data \citep{AkeretEA2017,BurdEA2018}.
Systems incorporating the SumThreshold algorithm have been implemented for the low-frequency LOFAR observatory
and higher frequency MERLIN interferometer \citep{OffringaEA2012,PeckF2013}.
New generation software correlators permit the integration of threshold or kurtosis-based
flagging applications before and after FX (Fourier Transform before multiplication) correlation and
stacking protocols \citep{Deller2010}.
A reference antenna implemented at the post-correlation stage can remove the signal from a well-defined
RFI source using the available closure relations \citep{BriggsBK2000}.

\vspace{2mm}{\it \bf Fringe stopping and delay compensation} - Array instruments employ fringe-stopping and delay-compensation techniques to keep a zero fringe
rate at the central observing position during observations. As a result the stationary (terrestrial) and
satellite RFI components in data distinguish themselves by fringing faster than components from
astronomical sources. This distinctive (relative) motion allows the off-line identification and
elimination of stationary RFI sources from both the correlated data and the image plane without
causing data loss \citep{CornwellEA2004,WijnholdsEA2004}. The coding for this
operation originating at the GMRT is now incorporated into AIPS task UVRFI \citep{Athreya2009}.

\vspace{2mm}{\it \bf Sub-space processing} -
In addition, more sophisticated statistical or sub-space processing can be implemented to remove the
RFI component with a minimum of data loss.
{\it Subspace filtering methods} may also be implemented in a digital correlation system to search for a particular signature
in the RFI power component of data in order to identify and remove it.
A particularly successful application is the search for cyclo-stationarity within the data, which works well
for digitally modulated RFI signals \citep{FeliachiEA2009, FeliachiEA2010,WeberEA2007}.
For array applications these methods depend on the estimation of the RFI spatial signature using the diagonalization of
the correlation matrix or the cyclic correlation matrix of the array \citep{HellbourgEA2012}.
In the case of a reference antenna, the evolution with time for a multiple RFI scenario requires a subspace tracking
approach using the covariance data where the reference antenna supports a faster convergence \citep{HellbourgEA2014}.

\section{Conclusions}
Both on-line and off-line data processing has been successful in mitigating the RFI environment of radio astronomy
observatories.
While there are an increasing variety of viable mitigation options, the choice of method depends
strongly on the RFI characteristics, the type of radio instrument, and the type of observation.
In particular, on-line real-time data processing may be preferred for a variable RFI environment, while special
measures such as reference antennas and spatial filtering may be preferred for known and fixed sources of RFI.
In addition to these factors, the absence of human involvement may also render automated on-line processing a
more attractive option.

 No universal method exists for mitigating RFI in astronomical data and no method can identify or remove RFI within
 the noise of the system.
 The effective suppression of RFI depends on the INR and its temporal and spectral characteristics.
 A quantitative evaluation of the method used is not always possible because mitigation algorithms are generally
 non-linear processes that also affect the noise characteristics and the gain calibration.
 The toxicity of the method used, i.e. the negative effect invoked on the data by the deployed method, and the amount
 of data loss resulting from the method are other guiding factors for evaluating the various methods.

 Multiple methods need to be applied to deal with a more general RFI environment. Because RFI characteristics
 change after each mitigation step and with increasing integration of the data, the cumulative effect of RFI mitigation
 at subsequent stages is not a linear sum of what each method can do, but rather the sum of what is practically
 possible at each step.

 The cost of computing hardware capabilities and digital applications at radio astronomy observatories are rapidly
 changing parameters. Upgrades of existing observing facilities and newly constructed instrumentation provide ample
 opportunity to implement and use automated RFI mitigation algorithms.
 These changing capabilities also result in increased observing bandwidths, higher time resolution, and higher spectral
 resolution.
 These increasingly large data volumes will force the introduction of automated data reduction pipelines and automated
 mitigation procedures.
 While the traditional user community of radio observatories does not or only reluctantly accepts automated mitigation
 implementations, they will be forced to do so in the future by the sheer volume of the data being handled.

 New telecommunication and broadcasting technologies are reaching the market place, and many of these involve
 unlicensed mobile devices.
 Their movement is difficult to control and they will rapidly affect observatory operations.
 Algorithmic research is needed in order to eliminate the signals of these devices from the data.
 In particular, spread spectrum (ultra-wide band) devices will pose problems for passive services because their digital
 modulation schemes do not respect the boundaries of spectrum allocations.
 Current estimates suggest that the number of transmitting devices used by each person is set to increase dramatically
 and many of these devices will rely on dynamic spectrum access.
The discovery space for radio astronomy is determined to a significant degree by the technical characteristics of the
observing system and by limiting factors such as the RFI environment.
While new generation instruments seek out the most pristine environments, existing facilities need to coexist with their
local environment.
In order to prevent the RFI environment becoming the limiting factor for each of the existing observatories, spectrum
management, both internal and external, as a method to control this environment must remain a very high priority.
Both observatory management and astronomers should be taking RFI issues seriously.


\section{Introduction}

Teaching Software Architecture (SA) in a meaningful and effective way is difficult. Many reasons have been posed to explain this reality, although there seems to be some consensus around the idea that SA exposes students to concepts that have significantly greater scale than handled before, and which are not easy to reproduce within the classroom-and-academic-term context~\cite{VanDeursen2017591}. These concepts require the development of correspondingly increased abilities of abstraction, being also often the first time issues like performance or security are addressed~\cite{Galster2016356}. Last but not least, there are non-technical skills at play, such as communication~\cite{Lago200535} and decision-making~\cite{Capilla2020231}, that need to be fostered, too~\cite{Ferrari200932}.

In this paper, we propose a SA teaching approach driven by the co-existence of different learning profiles, as defined by David Kolb~\cite{kolb2007kolb,McLeod2017}, in the classroom. Taking into account that traditional activities such as lectures and development projects are more aligned with some of those learning profiles, we contribute a set of role-playing activities designed to improve student's understanding of different software architectures characteristics, strengths and weaknesses, that target those learning profiles which are usually forgotten.

\section{Role-playing as a teaching tool}

The use of participatory teaching methods features a wide range of options which are applicable in the context of computer science~\cite{Jones1987155}: brainstorming, directed dialogues, small discussion groups, debates, panel discussions, or role-playing, amongst others. From this range of choices, for this work we have decided to use role-playing as a teaching tool.

Role-playing, sometimes referred to as \emph{role-play simulation} in educational settings, is an experiential learning method in which learners are involved in a proposed scenario by representing an interacting part in it. The scenario is outlined by the teacher or professional, and while it must allow improvisation, it represents a safe and supportive environment where students will develop their own meaningful first-person experience.

Role-playing is widely acknowledged as a powerful technique across multiple avenues of training and education. There is previous experience in software engineering featuring role-playing, most often to simulate the conditions of an industrial environment~\cite{Ohlsson1995291}, of software product lines~\cite{Zuppiroli201213}, or for requirements engineering~\cite{Zowghi2003233}. In these cases, students are asked to play the part of stakeholders, from clients to the different technical roles within the development team.

However, the way we have chosen to use role-playing to teach SA is not by assigning human roles to students, but by assigning them the roles of software components, as we will discuss in Section~\ref{sec:roleplaying}. While ``traditional'' (i.e.~human-focused) role-playing has been used before in the context of SA teaching~\cite{Montenegro20171} (as a trade-off analysis method), we have only found one previous similar approach to ours (i.e.~software-focused): specifically, in the context of a programming course where students play the part of objects in order to grasp the concepts of object-orientation~\cite{Andrianoff2002121}.

\section{Role-playing architecture models}
\label{sec:roleplaying}

SA at University of A Coruña is a 6-ECTS (\emph{European Credit Transfer and accumulation System}~\cite{ects}) course for Software Engineering students. This translates in practice into 150 hours of work, typically spread amongst the 15 weeks of a term. Of the 10 weekly hours, 3 correspond to class and lab hours, and 7 to student autonomous work. In this course, students are exposed to architectural concepts and practices: from architectural patterns to non-functional requirements' analysis, the impact of the latter in the former, architectural representation and modelling, etc. The goal is to offer students a learning environment in which they can acquire the skills that allow them to carry out architectural tasks in the context of software development: component identification and characterisation, assignment of responsibilities, motivated election of communication and integration alternatives, and architectural evolution and maintenance.

The teaching methodologies that were being used in the SA course at UDC were clearly aligned with some of Kolb's learning styles (cfg.~\ref{fig:kolblearningstylesinAS}). Lectures, and the materials provided (bibliographic references, articles, videos, assignments) work well for people which more prominent learning traits are watching and thinking, i.e.~those students that are reflective and take the most out of observation, and conceptualisation. Complementarily, practical work during lab sessions gave additional coverage to those most prominently guided by active experimentation. However, those with \emph{accommodating} and \emph{diverging} learning styles, in which experimentation or observation would most benefit from combination with concrete experience of the subject at hand, are entirely left up to self-manage their needs in their autonomous work.

\begin{figure}
  \centering
  \includegraphics[width=0.9\columnwidth]{learning-styles-kolb-as.pdf}
  \caption{SA course at UDC: link between methodologies/learning styles}
  \label{fig:kolblearningstylesinAS}
\end{figure}

Our proposal to address this situation was to role-play
different architectures of reference,
meaning that \textbf{scenarios were designed to provide students with a concrete experience of the structure, advantages and disadvantages of different architectural models}.

\subsection{Role-playing a layered architecture}

The instructions given for the first role-playing scenario are shown in \Cref{tab:layers}. Students were made to form systems with different number of layers, but the same functionality. The expected outcomes are first-hand experiences of:

\begin{itemize}
  \item Layers reuse from one system to the next, increased when responsibilities are specific rather than broad.
  \item Limited impact of changes, affecting to single layers (or adjacent, at most).
  \item Request processing in two different moments at each layer (on the way \emph{in}, on the way \emph{out}, on both), and consequences (i.e.,~security).
  \item Structural limitations to performance, related to the number of layers.
\end{itemize}

\begin{table*}
  \caption[Roles to play in a layered architecture]{Roles to play in a layered architecture\textsuperscript{a}}
  \label{tab:layers}
  \footnotesize
  \centering
  \begin{tabular}{cp{0.6\textwidth}}
  \textsf{Role: layers} & \textsf{Instructions} \\\hline \hline
  \textit{You are the input layer} & Receive requests from the client (teacher) and pass them along after processing. Do not take another request until you have received the reply to the last. Give the reply back to the client. \\
  \textit{You are an intermediate layer} & Take requests only from your previous, pass them along after processing to your next. Do not take another request until you have received the reply to the last. Give the reply back only to your previous. \\
  \textit{You are the last layer} & Take requests only from your previous, and process it. Do not take requests until you are done. Give the reply back only to your previous. \\
  \end{tabular}
  \mbox{} \\[0.25cm]
  {\scriptsize\textsuperscript{a} Due to space limitations, these role character description cards are not as detailed as the ones students' received, but instead some abstracted, summarized versions.}
\end{table*}

\subsection{Role-playing a pipe \& filter architecture}

The instructions given for role-playing a pipeline are shown in \Cref{tab:pipelines} (\cpageref{tab:pipelines}). Groups of students were asked to define their own filters, and design their own pipeline. Expected outcomes were first-hand experience of:

\begin{itemize}
  \item Filters reuse from one pipeline to the next, increased when responsibilities are specific rather than broad.
  \item Limited impact of changes, affecting to single filters (or adjacent, at most).
  \item Advantage to performance given by structural parallelism, increased by the possibility of duplicating filter stages when one is a bottleneck.
  \item Differentiated input and output points in the system.
  \item One-way path for requests and consequences (i.e.~error handling).
\end{itemize}

\begin{table*}
  \caption{Roles to play in a pipeline architecture}
  \label{tab:pipelines}
  \footnotesize
  \centering
  \begin{tabular}{cp{0.6\textwidth}}
  \textsf{Role: filters} & \textsf{Instructions} \\\hline\hline
  \textit{You are the input filter} & Receive requests from the client (teacher) and pass them along after processing. Take a new request as soon as you have passed along the previous. \\
  \textit{You are an intermediate filter} & Take requests only from your previous, pass them along after processing to your next. Take new requests only after passing along the previous. \\
  \textit{You are the last filter} & Take requests only from your previous, process and reply to the expected output. \\
  \end{tabular}
\end{table*}

\subsection{Role-playing a client-server architecture}

\Cref{tab:client-server} (\cpageref{tab:client-server}) shows the instructions given for the client-server scenario, the first of the distributed architectural models. Students showcased the creation and evolution of several systems. The fist-hand experiences outcomes were:

\begin{itemize}
  \item Independence of services, in terms of development, availability, etc.
  \item Robustness: to the ability to have a working system as soon as the first service, and the directory, are operational, service independence adds the ability to keep in operation even if one or several services are down.
  \item Directory component as single point of failure, that can be mitigated by having different directories (opens the door for offering tailored service catalogues to distinct client profiles).
\end{itemize}

\begin{table*}
    \caption{Roles to play in a client-server architecture}
    \label{tab:client-server}
    \footnotesize
    \centering
    \begin{tabular}{cp{0.75\textwidth}}
    \textsf{Role: clients} & \textsf{Instructions} \\\hline\hline
    \textit{You are a client} & Your job is to submit requests for the services offered by the directory. You may submit one or many requests at a time, to the same or different services, whenever you want. \\
    \textsf{Role: directory} & \\\hline
    \textit{You are the directory} & Your job is to receive the requests from the clients and pass them along to one of the instances of the corresponding service. As soon as you have passed a request along, you can take the next one right away. \\
    \textsf{Role: services} & \\\hline
    \textit{You are a service} & Your job is to receive client requests forwarded by the directory and produce a reply for them. Do not take requests from clients directly. Do not take a new request until you are done with processing the previous. \\
    \end{tabular}
\end{table*}

\subsection{Role-playing a leader-follower architecture}

The instructions given for role-playing a leader-follower architecture are shown in \Cref{tab:leader-follower}. Students played the dynamic evolution during operation of several systems. Expected outcomes were first-hand experience of:

\begin{itemize}
  \item Independence of workers in terms of availability.
  \item Robustness: ability to stay operational even if some workers go down.
  \item Scalability: ability to elastically react to demand by starting more workers.
  \item Resource consumption optimization: ability to elastically react to demand by stopping idle workers.
  \item Leader as single point of failure.
\end{itemize}

\begin{table*}
    \caption{Roles to play in a leader-follower}
    \label{tab:leader-follower}
    \footnotesize
    \centering
    \begin{tabular}{cp{0.75\textwidth}}
    \textsf{Role: clients} & \textsf{Instructions} \\\hline\hline
    \textit{You are a client} & Your job is to submit requests to the system, via the leader. You may submit one or many requests at a time, whenever you want. \\
    \textsf{Role: leader} & \\\hline
    \textit{You are the leader} & Your job is to receive the requests from the clients and pass them along to one of the workers. You have to decide how to choose such worker (randomly, round-robin, etc). You may create a new worker if you find that all of them are busy. As soon as you have passed a request along, you can take the next one right away. You may stop a number of workers if you find that most of them are often idle. You may send the request to several workers. You may split the request into several workers. \\
    \textsf{Role: worker} & \\\hline
    \textit{You are a worker} & Your job is to process leader' requests. Do not take requests from anyone else. Do not take a new request until you are done with processing the previous. \\
    \end{tabular}
\end{table*}

\subsection{Role-playing a peer-to-peer architecture}

The instructions given for the last role-playing scenario are shown in \Cref{tab:p2p}. Part of the students were made to form one network of peers, while the rest were to perform as clients. Expected experiences-as-outcomes were:

\begin{itemize}
  \item Need to embed both management logic and business logic into each peer (as opposed to the leader-follower architecture, in which these responsibilities are assigned to distinct components).
  \item Need to manage life of requests, to ensure clients get some reply even when the network is under heavy load.
  \item Difficulties of implementing security measures to differentiate malicious from faulty or unavailable peers.
  \item Absence of coordination (impact on scalability, self-load balance).
  \item No single point of failure: maximum availability.
\end{itemize}

\begin{table*}
    \caption{Roles to play in a peer-to-peer architecture}
    \label{tab:p2p}
    \footnotesize
    \centering
    \begin{tabular}{cp{0.75\textwidth}}
    \textsf{Role: clients} & \textsf{Instructions} \\\hline\hline
    \textit{You are a client} & Your job is to submit requests to the system, using any of the peers. You may submit one or many requests at a time, whenever you want. \\
    \textsf{Role: peer} & \\\hline
    \textit{You are a peer} & Your job is to receive the requests (which may come from other peers or from clients) and either handle them and produce a reply, or pass them along to one or several of your neighbours. If you produce a reply, give it back to whoever sent it to you. As soon as you have processed/passed a request along, you can take the next one right away. Periodically, you need to revisit your neighbourhood, look for new peers, dismiss others. \\
    \end{tabular}
\end{table*}

\section{Discussion}

Literature shows that multi-role projects are an effective teaching strategy~\cite{Warin2016137}. It also reveals that role-playing affects three areas that support student learning and engagement, namely: identified personalized learning, deepened content understandings, and enhanced collaboration skills~\cite{Toth2015386}.

However, empirical evaluation also reports that while many students feel that the infusion of role-playing aspects into the courses supported their learning and engagement, some other students do not~\cite{Toth2015386}. This is coherent with the fact that certain teaching activities resonate better with certain learning profiles, and may also explain our quantitative evaluation.

In any case, we consider this does not pose a threat to the validity of our efforts, since our role-playing activities aim precisely to provide an advantage to those learning profiles that benefit the less from the rest of teaching methodologies already present. There are a couple of factors that may influence how impactful this is when evaluating an initiative such as ours in the context of a whole class or course:

\begin{itemize}
    \item On the one hand, given that experimental learning experiences that address all learning styles are not commonplace, it makes sense to assume that students, especially by the time they reach the university level, have already adapted to the lack of methods and activities that best resonate with them, in the case of the less-often addressed learning profiles.
    \item On the other hand, given that we do not know how the distribution of different learning profiles amongst students in general, and neither amongst our CS students at UDC, it may also be the case their numbers are not significant when considered as part of the whole class.
\end{itemize}

\section{Conclusions}

In this paper, we have faced the challenge of extending the teaching methodologies that were being applied in an undergraduate SA course. The goal was to take into account the different learning profiles that students have. In doing so, we aim to provide a more fair learning environment, where every individual has opportunities to connect with the learning experiences in the most meaningful and effective way to them.

Given that the teaching methods that were in place were found to be miss-aligned when it came to students with \emph{accommodating} or \emph{diverging} learning styles (following Kolb's nomenclature), the way we have addressed this challenge has been by incorporating role-playing in an innovative way. We have designed a set of scenarios where students play the part of components in different architectural models. Such role-playing game provides concrete experiences of the structure, advantages, and disadvantages of different architectural models.


\section{Introduction}\label{sec1}
\begin{figure}
\begin{center}
\scalebox{0.38}{\epsfig{file=fig1.eps}}
\caption{The jumping chain.}
\label{fig1}
\end{center}
\end{figure}
There is a very interesting physical phenomenon
as shown in FIG. \ref{fig1}: a coiled but not entangled chain is put
in a cup. When the end of the chain is pulled out of the cup,
it will fall freely along the wall of the cup.
At a certain point, usually in a sudden moment, the chain can jump up
and remain in a steady arc shape which can last for a certain
duration in the air. This experiment is often demonstrated  in college
physics classes as an example  of an
intriguing classical mechanical problem.
The reader can also observe the details of the jumping in the
supplementary video or on the internet website~\cite{vedio}.
In this paper, we attempt to analyze and study the pattern and
related topics such as the height of the jump and the height
of the fall, {\it etc.}  Experiments will be conducted and
the data will be analyzed to verify our theoretical analysis~\cite{wxxy}.

In fact, the phenomenon is an application of classical
mechanics in string, rope or fluid physics~\cite{jordan}. The jumping rope is a familiar example of a slender structure that interacts
with the fluid through which it moves~\cite{aerodynamics}.
The chain is usually  suspended and swung by a mechanical apparatus,
and the pattern is well studied in the literature~\cite{hanging}.
This is because important class of truly 1D problems involves the motion of lines or filaments embedded in 3D space, such as vortex lines in
a defect lines in liquid crystals, type-II superconductor,
or polymers in imposed flows.  Nonlinear partial differential
equations describes the dynamics, and the twisting, breaking, reconnecting, or the knotting,  are the
essential topological parts of the chain dynamics~\cite{topology}.
However, the jumping chain studied in this paper,
is not suspended but jump from a higher platform, falling
freely in the space under the attraction of gravity.
The damping from the viscosity of
the air seems to be a sub-dominant effect in the whole process.
The gravity and the tension inside the chain play an important role.
The subtle point is that, although the
matter is moving from the higher position to the lower
ground, the system is not a fluid or a falling stream. A detailed
study will reveal the different dynamics.

In all, the paper is organized as follows:
The phenomenology and dynamics of the jumping chain
will be analyzed in Sec.~\ref{sec2}; The
experimental verification of  our analysis will be presented
in Sec.~\ref{sec3};
The conclusion is given in Sec. \ref{sec4}.

\section{Phenomenology and Theoretical Analysis}\label{sec2}
\subsection{Phenomena}
First, let's  describe  how the jumping phenomenon takes place.
As shown in FIG.~\ref{fig1}, the chains are coiled layer
by layer inside a cup, and it is kept from becoming  entangled
 (It is difficult to keep all steel balls untangled
in the preparation, and it is one of the dominant disturbing
factors  of the experiment.).
Then the cup is placed at a certain height.
a small length of chain is taken out of the cup
and allowed to  fall freely.
In the actual operation,  we also
tried to use a soft and non-contractible rope to do the same experiment.
The rope can barely jump up but the phenomenon is rather obscure,
so the thin metal chain composed of steel balls is adopted
in our experiment. The steel balls are loosely  connected, and the
chain is  non-retractable. It
means that the  chain do not have flexural capacity.
This is the key difference between the chain and the soft rope.
The reason will be discussed in the following context.

The experiments show that the chain jumps
to appropriate higher position from the cup if the cup is set at
some specific heights.
Due to the complex way of chain winding in the cup,
the starting point of the jump keeps  changing, so the
shape of the jumping from cup to the highest point is very complicated.
As for the pattern between the highest point and the falling off point
on the ground, the chain can maintain an almost stable shape
near the highest point on the top.
The shape keeps  shaking near the ground, but this
does not affect the stable shape on the top.

\begin{figure}
\begin{center}
\scalebox{1.2}{\epsfig{file=fig2.eps}}
\caption{The sketch map of our analysis on the jumping chain.}
\label{fig2}
\end{center}
\end{figure}

Summarizing our observations, we think this physical phenomenon
can be divided into two sub-processes: the accelerating
process and the stable suspension process.
Thus our analysis on the phenomenon is set out as
shown in FIG.~\ref{fig2}, in which point A is the beginning
of the suspension and the end is point B, which
is the collision point on the ground.
The point C is the starting point of the jumping
in the cup where the elements of the chain are static.
The part from C to A is the accelerating process.
Note that, as discussed above, point B keeps shaking,
and point C also varies quickly.  At the same time,
the height of Point C is decreasing. The definitions of
the variables such as $h_j,~h_i$ {\it etc.} are given
in the following context. One can see that
the stable suspension is obviously of the most concern,
and according to following analysis, we find that this is
actually the key point to comprehend the phenomenon.
Thus we analyze  the stable suspension in the following subsection.

\subsection{Analysis on the shape of the stable Suspension}
At first glance, the jumping chain appears like a parabolic curve.
However, the truth is not so simple as expect.
A free-falling mass point follows a parabolic trajectory.
However, for the non-retractable chain,
each element on the suspended chain falls
at the same velocity, only the directions vary.
Based on our observations, the coiled chain in the cup
almost moves to the ground at a constant rate. thus it is
natural to assume that the suspended chain can maintain
a constant velocity $v_0$ for appropriate interval as
an ideal case. Then, in $\Delta t$ time, there is a part of the
chain with mass $\lambda v_0dt$  ($\lambda$ is denoted as
the line density) in the cup which is accelerated from 0 to $v_0$,
At the point A, the chain feels the dragging tension from
the chain at rest in the cup as follows:
\begin{equation}
T_0=\lim_{\Delta t \to 0}
 \frac{\lambda v_0^2 \Delta t}{\Delta t}=\lambda v_0^2\, .
\end{equation}
It should be noted that this is merely the ideal case, and the following study will show that the actual tension
at the end of A will be less than $\lambda v_0^2$.
At point B, the chain element with velocity $v_0$ collides with the
ground. If the collision is elastic, the chain will rebound from the
ground with velocity $v_0$. Similar analysis gives that
impulsive force to the ground is $2\lambda v_0^2$.
Of course, the collision is not completely elastic due to the
dragging. We shall provide the specific method for dealing with the
impulsive force at point B later.

We find that, very interestingly, the shape of the suspended chain
is actually an inverted catenary. The exploration of the catenary is an
important progress in the history of classical mechanics and mathematics.
Our discovery in fact compensates a block in the catenary theory.
Here we briefly show the derivation of the catenary and
explain the reason for the same shape of the suspended chain.
Suppose the shape of the chain is $y=y(x)$, the solution of the
catenary is to find a $y(x)$ such that the potential
or center of gravity of the chain is at its lowest point.
Using the variational principle, the action of the whole system is
the center of gravity in the vertical direction
\begin{equation}
  S=-\int \lambda g(y-c)\sqrt{1+{y^\prime}^{2}}dx=\int\mathcal{L}dx,
\end{equation}
where $g$ is the acceleration of gravity and $c$ is an arbitrary
constant, $y^\prime$ is the  derivative of $y$. $\mathcal{L}$
is the Lagrangian of the system
\begin{eqnarray}
&&\mathcal{L} = -\lambda g(y-c)\sqrt{1+{y^\prime}^{2}}\label{motion2}.
\end{eqnarray}
The motion for $\mathcal{L}$ gives the non-linear differential equation of
the catenary
\begin{equation}
 (y-c)y^{\prime\prime}-{y^\prime}^2-1=0.\label{motion3}
\end{equation}
The solution is
\begin{equation}
  y=c+\frac{a}{2}(e^{\frac{x}{a}+b}+e^{-\frac{x}{a}-b}).
\end{equation}
where $a,~b$ are the parameters determined by the inputs of the catenary.
The approach described above is the standard derivation of the catenary.
Note that here suspension chain are hanged at any height,
thus an additional constant $c$ are introduced.
In fact, we notice that, although the equation of motion
Eq.~\eqref{motion3} is nonlinear, it has a discrete symmetry and $y(x)$
can be transformed
\begin{equation}
  y(x)\rightarrow Y(x)=-y(x)+2c.
\end{equation}
$Y(x)$ is also a solution of the catenary equation.
Despite its lack of practical application,  the inverted
catenary solution is still a valid solution of the catenary equation.
This solution is an unstable extremum
point of variation, which has been largely overlooked in the literature.
Nevertheless, our research has revealed that this inverted solution is
indeed the exact shape of the suspended chain.

For the jumping chain, the potential energy is the same as
that of the catenary. Therefore only the kinetic term
needs to be included in the action
\begin{equation}
  S=\int \mathcal{L}dx=\int \left(\frac{1}{2}\lambda\sqrt{1+{y^\prime}^2}v_0^2
-\lambda gy\sqrt{1+{y^\prime}^2}\right)dx\, .
\end{equation}
 By comparing Eq.~\eqref{motion2} and assuming that $v_0^2$ is a constant,
 we can observe that the Lagrangian of the jumping chain
 is identical to that of the  catenary,
 with the only difference being that  $c$ is taken as a constant
\begin{equation}
  c=\frac{v_0^2}{2g}\, .
\end{equation}
Due to the given conditions, the solution will be
a convexity, specifically, an inverted  catenary
\begin{equation}
  y=c-\frac{a}{2}(e^{\frac{x}{a}}+e^{-\frac{x}{a}})\,.
\end{equation}

\begin{figure}
\begin{center}
\scalebox{0.42}{\epsfig{file=fig3.eps}}
\caption{The force analysis on the elements of the jumping chain.}
\label{fig3}
\end{center}
\end{figure}

The difference between the deviation of an ordinary catenary and
a jumping chain is that the chain is constantly
in motion. The elements of the chain is continually jumping out of the
cup and falling on to the ground, thus changing the mass of
system. This may lead some readers to be skeptical about the validity
of variational principle. However, the same differential
equation of motion can be derived by the force analysis.
As shown in FIG.~\ref{fig3},
the tension in the jumping chain is denoted as $T=T(x)$
and the convention of the forward direction is defined
as from point A to point B. Every element of the chain
is moving in a curvilinear motion with the same speed, so
the tangential acceleration is zero, and the infinitesimal
tension $dT$ cancels the tangential component of gravity
\begin{eqnarray}
    dT=\lambda g\sqrt{1+{y^\prime}^2}dx  \frac{y^\prime}{\sqrt{1+{y^\prime}^2}}
    =\lambda g\frac{dy}{dx} dx=\lambda g dy\,.
\end{eqnarray}
Thus the tension in the chain can be expressed as
\begin{equation}
	T=T_A+\lambda gy\,.\label{zhangli}
\end{equation}
This is a very interesting result that
the tension in the stable jumping chain changes along with the height of
the elements, which is  the same as the tension in the
ordinary catenary.
To analyze the normal component force of the infinitesimal
chain imposed by tension from the curvature of the chain,
we can look to the  catenary or hanged rope for comparison.
In the case of the ordinary catenary,
 the normal component of the tension is canceled by the gravity,
 but in the case of the jumping chain,
the normal component force provides the centripetal force
together with the gravity. The geometrical configuration
is shown in FIG.~\ref{fig3} where we
investigate the centripetal force acting on the two chain elements
$x-dx\to x$ and $x\to x+dx$. From this, we
 can easily derive  the  equation for the
centripetal force
\begin{eqnarray}
\frac{2\lambda \sqrt{1+y'^{2}}dx v^2}{\rho}=
	 2T d\theta+2\lambda g \sqrt{1+y'^{2}} dx\frac{1}{\sqrt{1+y'^{2}}}\, ,
  \nonumber
\end{eqnarray}
where $\rho$ is the radius of curvature at point $x$.
By substituting  $\rho$, $d\theta$, and tension $T$, we can get
\begin{eqnarray}
	\frac{2\lambda|y''|v_0^2 dx}{1+y'^2}=\frac{2(T_A+\lambda g y)|y''|dx}{1+y'^2}+2\lambda g dx\,. \label{ten1}
\end{eqnarray}
Due to the convexity of the jumping chain
\begin{eqnarray}
	|y^{\prime\prime}|=-y^{\prime\prime}\, .
\end{eqnarray}
Eq.~\eqref{ten1} are reduced to
\begin{eqnarray}
-(T_A+\lambda gy -\lambda v_0^2)y''+\lambda g y'^2+\lambda g=0\, .
\end{eqnarray}
Given that $T_A$ is the variable to be solved, it is convenient
to define
\begin{eqnarray}
	T_A-\lambda v_0^2=-\lambda gc \label{tac}\, .
\end{eqnarray}
Then the differential equation of $y$ is
\begin{eqnarray}
	(y-c)y^{\prime\prime}-{y^{\prime}}^2-1=0\, .\label{def-cat}
\end{eqnarray}
Proved the same catenary equation again.

Although the analysis of force is a bit more complicate
than that of the variational principle of mechanics,
the physical picture is easier to understand because
it is an analysis of forces.
We can also see that the forces at the two ends of the catenary
and the jumping chain  are different.
In case of the ordinary catenary, the tensions at the two ends
are dragging the chain. Whereas, in the
case of the jumping chain, the tension
at point A is dragging the chain, while the tension at
point B is pushing the chain. Of course, this is the ideal case;
the jumping process will be further analyzed in the following.

\subsection{Analysis on the jumping}
First, let's begin with the explanation for why $T_A\neq\lambda v_0^2$.
$T_A=\lambda v_0^2$ corresponds to the ideal case  that
each chain element instantaneously accelerates from 0 to $v_0$.
In contrast, the actual situation is that the chain is loosely coiled and
the chain element needs an interval for the  acceleration,
making $T_A$ smaller than $\lambda v_0^2 $.
Analytically,  if $c=0$, the shape of the suspended chain line is
\begin{eqnarray}
	y=-\frac{a}{2}(e^{\frac{x}{a}}+e^{-\frac{x}{a}})\,,
\end{eqnarray}
then
\begin{eqnarray}
	y\le -a\,.
\end{eqnarray}
It is evident that  $y>0$ solution does not exist, as indicated by
Eq.~\eqref{def-cat}, which a non-linear differential equation.
Therefore, if the coordinate system shown in FIG.~\ref{fig2} is used
a $y>0$ solution requires a non-zero $c$.
From Eq.~\eqref{tac}, it is clear that
$T_A$ must be less than $\lambda v_0^2$.
As shown in Fig.~\ref{fig2}, the acceleration process
for the upper elements in the stable suspended shape is added.
And the height of the acceleration jumping up process
from point C to point A is measured as $h_i$.
The jumping height of the stable suspended chain line is $h_r$.
The total jumping height is thus the sum of $h_i$ and $h_r$
\begin{eqnarray}
	h_j = h_i +h_r\, .
\end{eqnarray}
Along with the conventions mentioned above,
the falling height (the height of the cup) of the chain is denoted as $H$
and the width of the jump is $2x_0$ as shown in Fig.~\ref{fig2}.
We will  further investigate the relationship between these
quantities for jumping chain and attempt to
conduct  experiment to verify our analysis.

The tension at point A is
\begin{eqnarray}
	T_A=\lambda v_0^2-\lambda gc\,.
\end{eqnarray}
According to the Eq.~(\ref{zhangli}),
the tension $T$ from point A to point B is
\begin{eqnarray}
	T=\lambda v_0^2+\lambda gy -\lambda gc\,.
\end{eqnarray}
As discussed previously,  the elastic collision implies that
the force imposed from the ground to the chain at point
B will be $2\lambda v_0^2$. This may appear suspicious to some
readers, as a deceleration process is needed at point B.
When the chain finally fall on the ground, the tension should
be less than $2\lambda v_0^2$. However,
we found that the chain elements can slightly bounce on the ground,
thus the actual impulsion to the ground should be greater
than $\lambda v_0^2$. To account for the force with more precision,
a constant coefficient $\alpha$ can be introduced
\begin{eqnarray}
  -\alpha\lambda v_0^2=\lambda v_0^2-\lambda gH-\lambda gc\, .
\end{eqnarray}
The value of $\alpha$ would be greater than 1 and less than 2.
It should be noted that
tension at point B should be negative according the
direction defined above. Since falling height $H$ is much larger
than the jumping height $h_j$,  the tension fluctuation
at point B is not sensitive to stable suspension above. Then
\begin{eqnarray}
  v_0^2=\frac{1}{\alpha+1}g(H+c)\label{v0}\,.
\end{eqnarray}
This is the relation between the falling height and the velocity
of the suspended chain.
The jumping height $y(x=0)$ of the stable suspended chain is
\begin{eqnarray}
  h_r=c-a\,.
\end{eqnarray}
Then all the exploration  focuses on the pending
parameter $c$ and $\alpha$, which will be analyzed in the
next section when discussing the experimental verification.
Before that, let us analyze the acceleration process.

As the velocity of the chain element changes over time during
the acceleration process, it seems difficult to do the analysis. However,
the time parameter can be eliminated after some substitutions.
The velocity of accelerating elements is denoted as $v$. According to
Fig.~\ref{fig2} and Eq.~\eqref{zhangli}, the tangential acceleration
should be
\begin{eqnarray}
  \lambda vdt\frac{dv}{dt}=\lambda v dv=dT-\lambda gdy\,. \label{jia}
\end{eqnarray}
Since the chain elements are loosely connected, the normal
acceleration dynamical equation can be disregarded. It can be seen that
the time dependence $dt$ can be eliminated. This is the
mathematical formulation of the transformation from gravity potential to
kinetic energy. Carry out the  integration on both sides, we have
\begin{eqnarray}
  \int ^{v_0}_{0}\lambda v dv&=&\int_0^{T_A} dT-\int_{-h_i}^0\lambda gdy\nonumber\\
&=&T_A-\lambda gh_i\,.
\end{eqnarray}
At point A
\begin{eqnarray}
  \frac{1}{2}\lambda v_0^2 =\lambda v_0^2 -\lambda g c -\lambda  g h_i\,,\nonumber
\end{eqnarray}
namely
\begin{eqnarray}
  h_i = \frac{v_0^2}{2g} -c\,.
\end{eqnarray}
Substitute $v_0$ from Eq.~\eqref{v0}, we can get the
relation between $h_i$ and $H$
\begin{eqnarray}
  h_i =\frac{H-(2\alpha +1)c}{2(\alpha+1)}\,.
\end{eqnarray}
From above equation,
the requirement for the jumping can also be got. The jump
can happen only when $h_i$ greater than zero, thus the condition
for the jumping is
\begin{eqnarray}
  H>(2\alpha+1) c \label{qitiao}\,.
\end{eqnarray}
The total jumping height is
\begin{eqnarray}
 h_j =h_i+h_r =\frac{H+c}{2(\alpha+1)}-a\label{hjh}\,.
\end{eqnarray}
It is evident that the jumping height basically is linearly correlated
with the falling height $H$, with
additional parameter $c$, $a$ and $\alpha$ to be determined.

\section{Experimentation and Verification}\label{sec3}

\subsection{ Analysis and execution of the experiment}\label{sec3.1}

\begin{figure}
\begin{center}
\scalebox{0.42}{\epsfig{file=fig4.eps}}
\caption{The apparatus of the experiment: A, the height
calibration on the wall;
B, the cup; C, the chain; D, the chain elements, E, the screenshot.}
\label{fig4}
\end{center}
\end{figure}

Through the theoretical analysis in Sec.~\ref{sec2}, we can verify
the relation through experimentation. Nevertheless, it is important
to bear in mind that all the theoretical analysis are
based on ideal conditions, and  there are certain
limitations in the actual experiment, such as
\begin{enumerate}
    \item Ideally, the chain in the cup should be loosely
    connected and freely coiled. Unfortunately, the elements in chain
    inevitably become  tangled, which prevents the chain from
    maintaining homogeneous during the the acceleration
    and stable suspension
    processes. We discovered that this is the primary source of
    disturbance in the experiment. Even a small knot can cause
    a significant variation of the stable suspension.
    \item Ideally, the jumping point C and falling point B should remain
     fixed in order to obtain a  stable suspension. However, these
     two points tend to  shake during the falling process,
     making it difficult to accurately measure the shape.
     Nevertheless, Since the shape of catenary is an exponential function
     in which $y$ increases exponentially with $x$, the shaking chain
     almost forms a straight line at the bottom. This makes it
     challenging to measure the correlation between the shape and the exponential curve. Furthermore, due to the upper two points, it is difficult to precisely determine the
     positions of points A, B and C. Generally, only the total
     jumping height $h_j$ and falling height $H$ can be measured.
\end{enumerate}
Though there are these difficulties, it is sure that
\begin{eqnarray}
  H\gg c,~a\,.
\end{eqnarray}
and the jumping height $h_j$, $a$ and $c$ are at the same order. Based
on this, we can assume that the actual shape of the
jumping chain resembles and  inverted catenary.
The experiment should then test the linear
relation between $h_j$ and $H$ (Eq.~\eqref{hjh}) and
the linear relation between $v_0^2$ and $H$ (Eq.~\eqref{v0}).
Both relations could be used to evaluate the parameter $\alpha$.
which will validate our analysis in Sec.~\ref{sec2}.

\begin{figure*}
\begin{center}
\scalebox{0.52}{\epsfig{file=fig5-1.eps}}
\scalebox{0.52}{\epsfig{file=fig5-2.eps}}
\caption{Left: The recorded results between the jumping height
and falling height; Right:  the velocity $v_0^2$ of the stable
shape and the falling height}
\label{fig5}
\end{center}
\end{figure*}

The apparatus used in our experiment is illustrated
in Fig.~\ref{fig4}.
We drawn a height calibration on the wall,  which can be seen
in penal A of Fig.~\ref{fig4}. The chain and the cup are also depicted
in the figure, and the details of chain are provided
in the Tab.~\ref{tab1}.
  \begin{table}[htbp]
     \centering
     \begin{tabular}{cc}
     \hline
     \hline
   Diameter of the steal ball   & $5\pm0.5$mm\\
   Length of the steal ball & $3\pm0.5$mm\\
   Length of the chain     & 8.4$\pm$0.1m\\
   Density $\lambda$   & 0.036$\pm$0.03kg/m\\
   Height of the cup     & $95\pm0.5$mm\\
   Diameter of the cup     & $66\pm0.5$mm\\
     \hline
     \hline
     \end{tabular}
     \caption{The setup of the exepriment.}
     \label{tab1}
 \end{table}
The experiment was conducted as follows:
The chain was pulled out of the cup
from a height 200cm, then gradually lowered by  5cm
each time. When the height reached 100cm, it became
increasingly
difficult to lower the chain further. Therefore,
we chose to lower the chain by 10cm for heights lower than
100cm, until the height reached 40cm. At each specific
height, a cellphone was used to record the jumping chain,
while the camera were placed at the same height at the  top of
jumping chain (see penal E of Fig.~\ref{fig4}).
Due to the presence of various disturbances
in the experiment, multiple vedios were taken in order to select
the best one for the data collection.

 In Principle, catenaries with different $c$ and  $a$ can exist
 when the requirements discussed in Sec.~\ref{sec2} are satisfied.
 However, since the chain needs to jump from the cup, thus the shape
 of the jumping chain depends on specific initial conditions,
 namely, the diameter of the cup, the depth of the
 cup at the beginning of the jump, {\it etc}. Thus though
 the shape of the jumping at a specific height
 is not in a single parameters set,  We found
 that $c$, $a$ are at the same order of the diameter of the cup.
 thus the linear relation could be found in the following.

The data collection was carried out as follows. each frame
of the chosen vedio was examined, and the best moment with the highest
jumping height and stable shape was identified. The falling height
and jumping height were then recorded base on the frame at the best
moment. To measure $v_0$, the
time interval between a specific height of the chain left in the cup
and the final moment was recorded. With these measured raw data, we can
analyze on the relation between the physical quantities talked
in the Sec.~\ref{sec2}

\subsection{Analysis on the data}\label{sec3.2}
As discussed in the previous subseciton, linear
relation between $h_j$ and $H$ (Eq.~\eqref{hjh}) and
the linear relation between $v_0^2$ and $H$ (Eq.~\eqref{v0})
can be used to verify  our analysis.
The corresponding results are shown in Fig.~\ref{fig5},
with the left panel displaying the recorded jumping heights
and falling height, and the right panel shows the calculated
velocity $v_0$ of the stable suspension shape and falling height. Note
that due to presence of various disturbances in the experiment,
only data with falling higher than 100cm (black points) were chosen for
the analysis. As can be seen in Fig.~\ref{fig5}, the linear relations
are realized and the slopes give
the values of the $\alpha_1$ (left panel) and $\alpha_2$ (right panel)
\begin{eqnarray}
    \alpha_1 = 2.6^{+5}_{-1.3},~~~ \alpha_2 = 2.5^{+7.0}_{-0.7}\,.
\end{eqnarray}
We can see that the errors are quite large and the central values
exceed 2. However, given the difficult setup of the experiment,
these two measurements on $\alpha$ are acceptable.
Note that although $\alpha_1$ and $\alpha_2$ come from the same
series of derivations, Eq.~\eqref{hjh} and Eq.~\eqref{v0} are assumed
to independent with each other. Thus,  we can conclude that
$\alpha_1$ and $\alpha_2$ are in agreement with each other.
Another interest point is that the extension line in the left panel
intersects the $H$ coordinate axis at a point greater than 0.
This point should correspond to  the required height of the jumping
condition in Eq.~\eqref{qitiao},
Although  this is only a rough estimation.

Finally, we can explain  why the soft rope
cannot effectively jump. the elements of chain are connected loosely,
and the rope cannot achieve this condition. This is also
the reason why the catenary cannot be realized by the soft rope.
Specifically, the soft rope lacks the non-flexural capacity
and the absence of disturbance from normal
force, making the jumping impossible.

\section{Conclusion}\label{sec4}
The jumping of a free-falling chain is an intriguing  phenomenon,
and this paper provides a  detailed analysis of  the physics
behind it. The analysis
are divided into two parts: the stable suspension and the jumping.
Both variational principle and force analysis demonstrate that  the
stable suspension is an inverted catenary. The
parameters that describe the phenomenon
are studied and tested through experiments, The requirement for
the jumping are discussed.
The measurements of the parameter $\alpha$ meet  our expectations.
The physical picture  of the falling chain could
be useful in some falling systems.
and the inverted catenary may be a complement to classical mechanics.