IRI the International Standard for the Ionosphere

This paper gives a brief overview over the International Reference Ionosphere (IRI) project and model. IRI is recognized as the official standard for the ionosphere by the International Standardization Organization (ISO), the International Union of Radio Science (URSI), the Committee on Space Research (COSPAR), and the European Cooperation for Space Standardization (ECCS). Of great importance are the external drivers of the model that help IRI to represent ionospheric conditions as realistically as possible. The paper discusses the drivers currently used presents recent improvements and changes. Besides the standard solar, magnetic, and ionospheric indices the paper also reports on the adjustment of the model with data and equivalent indices and on the progress towards a Real-Time IRI using data assimilation. IRI has been widely validated with many different data sources and has fared very well in community wide assessment studies. We present some of these studies and document the wide usages of the model in the scientific literature. Finally, we present an outlook on things to come in IRI-2018 and thereafter.


Introduction
The International Reference Ionosphere (IRI) project started in the late sixties when the Committee on Space Research (COSPAR) and the International Union of Radio Science (URSI) recognizing the need for an international standard set up an Inter-Union Working Group with the goal of developing an international standard ionosphere.The charter given to the group was as follows: -IRI is an international project jointly sponsored by the Committee on Space Research (COSPAR) and the International Union of Radioscience (URSI) to develop and improve a reference model for the most important plasma parameters in the Earth ionosphere -COSPAR's prime interest is in a general description of the ionosphere as part of the terrestrial environment for the evaluation of environmental effects on spacecraft and experiments in space.
-URSI's prime interest is in a reference model for defining the background ionosphere for radiowave propagation studies and applications.
-The model should be primarily based on experimental evidence using all available ground and space data sources and should not depend on the evolving theoretical understanding of ionospheric processes.
-As new data become available and as older data sources are fully evaluated and exploited, the model should be revised in accordance with these new results.
-Where discrepancies exist between different data sources the IRI team should facilitate critical review and discussions to determine the reliability of the different data sets and to establish guidelines on which data should be used for IRI modelling.
Figure 1 gives a brief overview of the most important highlights and milestones of the IRI project from the beginning in 1968 to becoming the international ISO standard for the ionosphere in 2014, and the most recent publication of the latest version of the IRI model (Bilitza et al., 2017).New versions were released as new data became available and as old data were fully exploited for the modelling process.Better mathematical techniques and appropriate functional representation also contributed to the steady increase in accuracy of the model.From early on IRI was made openly available and was distributed first on magnetic tape, then punched IRI is the synthesis of almost all available and reliable ionospheric data.The IRI Working Group's (60+ experts) global distribution has helped to gain access to the global data base.The different data sources that have contributed to the built-up of IRI are depicted in Fig. 2 including data from ground-based observatories as well as spacecraft instruments.Combining the global picture recorded by satellites for different local times and levels of solar activity with the 24/7 365 analysis provided by ground stations is one of the core principle of the IRI development.An important first step here is to confirm that both techniques (from the ground and from space) get the same result when observing the same volume in space at the same time.IRI helped to resolve conflicting measurements by bringing together the ionospheric ground and space communities during biannual IRI workshops.Two examples of still existing conflicts are shown in Figs. 3 and 4. The data points in Fig. 3 show simultaneous ground and space measurements of the electron temperature at 840 km altitude.These measurements were acquired while one of the DMSP satellites flew over the Incoherent Scatter Radar (ISR) station at Millstone Hill.Large discrepancies are found at very low and very high temperatures.This may be due to data analysis problems both techniques encounter at the very low electron densities in the topside ionosphere.In Fig. 4 we compare different models for the F-peak height hmF2, that were developed with different data sources.All three model options are included in the latest version of IRI, IRI-2016.The IRI-M(3000)F2, IRI-AMTB, and IRI-Shubin models were developed with ionosonde data for the propagation factor M(3000)F2 (CCIR, 1967;Bilitza et al., 1979), with ionosonde data for hmF2 (Altadill et al., 2013), and with COSMIC radio occultation data (Shubin et al., 2013;Shubin, 2015), respectively.Figure 4 shows the variation of these three models from 1958 to 2015 during winter noon conditions at mid-latitudes.While the three models agree well during the solar cycle minimum, large and consistent discrepancies are found at high solar activities, pointing to potential discrepancies between these different techniques.
IRI describes monthly averages of electron density, electron temperature, ion temperature and ion composition (O + , H + , He + , N + , NO + , O + 2 , Cluster ions) in the altitude range 50-2000 km on the whole globe.Being a data-based model IRI's reliability depends on the spatial and temporal coverage provided by the underlying data base.High data density at northern mid-latitudes guarantees good accuracy whereas the performance is not as good in data sparse regions like the high (auroral and polar) and low (equatorial) latitudes.IRI uses a modular approach which allows to replace a specific sub-model without affecting the other model  parts, e.g., global models at profile anchor points are connected with profile functions.In the case of the electron density these anchor points are the E and F-peaks while in the case of the electron temperature the anchor points are at fixed heights (300, 400, 600, 1400, and 3000 km).44 switches allow the user to customize their application of IRI, e.g., turn certain features (storm model) on/off or choose between different model options for a specific parameter in a specific altitude range.Allowing different options for the same parameter simplifies the phasing in of new models with validation help from the users.  .IRI-M(3000)F2 is based on ionosonde measurements of propagation factor M(3000)F2 and the relationship between hmF2 and M(3000)F2 (Bilitza et al., 1979); AMTB-2013 is based on digisonde measurements of hmF2 (Altadill et al., 2013); Subin-2015 is the model developed by Shubin et al. (2013) and Shubin (2015) based primarily on COSMIC radio occultation data.
The improvements implemented with the latest version of the IRI model, IRI-2016, were introduced and discussed by Bilitza et al. (2017) and will not be repeated here.In the next chapter we will discuss the drivers that are being used in IRI including indices as well as actual measurements.The third chapter highlights some of the IRI successes and we close with an outlook at future improvements that are expected from ongoing projects.

IRI Drivers
There are a number of different ways in which measurements are being used to drive IRI towards describing conditions as realistically as possible: (1) solar, magnetic, and ionospheric indices that are deduced from solar, magnetic, and ionospheric measurements; (2) equivalent solar or ionospheric indices obtained by adjusting the model to ionospheric measurements; (3) direct updates locally if measurements of characteristic parameters are available; (4) assimilation of available data into the IRI background model.

Solar and ionospheric indices
The traditional method is to use indices which are themselves deduced from measurements.In IRI we use solar, magnetic, and ionospheric indices.Table 1 lists the solar and ionospheric indices applied in the different IRI sub-models.
Sunspot number R, which is related to the number of spots and spot groups observed on the sun, was the index of choice for many of the early ionospheric models and is still being used by some of the older IRI sub-models.Most often the models depend on the 12-month running mean of R (dewww.adv-radio-sci.net/16/1/2018/Adv.Radio Sci., 16, 1-11, 2018 ).Statistical studies had found much higher correlation between R 12 and ionosonde-measured foF2 (the F2 peak plasmas frequency) than between the daily value or the monthly average of R and foF2.This is understandable because daily and monthly variation of ionospheric parameters are not only affected by the solar input but also by the forcing from below (atmosphere) and above (magnetosphere).R is obtained and distributed by the Solar Influences Data analysis Center (SIDC, http://sidc.oma.be/, last access: 20 July 2018) in Brussels, Belgium.A recent review of the way R is determined from observations of spots and groups of spots on the Sun resulted in a recalibration of the R computation (Clette et al., 2014).The new method is now internationally accepted and the new sunspot number is distributed by SIDC.
The new index exceeds the old one by about a factor of 1.4.Models like IRI that were developed with the old R index need to account for this transition to a new index.In the case of IRI this is done by using the new index with a scaling factor of 0.7.Gulyaeva (2014) has shown that 0.7 is a very good estimate for the ratio of the 12 month running mean of these two indices.Future modelling efforts are strongly encouraged to use the solar radio flux F10.7 instead of the sunspot number for the representation of variations with solar activity.F10.7 is the solar radio flux at 10.7 cm wavelength that has been measured consistently since 1947, first at Ottawa, Ontario and then at the Penticton Radio Observatory in British Columbia both in Canada.Besides the daily F10.7 index IRI uses the 81-day (3 solar rotations) and the 365-day running mean of F10.7 as indicated in Table 1.For the topside electron temperature and ion composition a combination of daily and 81-day indices, PF10.7 = (F10.7_daily+ F10.7_81)/2, generated the best results.Figure 5 illustrates the differences between these indices.Only the daily index follows the 27-day solar rotation.With the 365-day smoothing the dominant pattern is the 11-year cycle.The IRI indices file APF107.DAT includes the F10.7 index adjusted to 1 AU.This is the index consistently provided by the distribution sites.In IRI-2016 the adjusted flux is back-ward corrected for the changing Earth-Sun distance to obtain the F10.7 index observed at the ground.It is not always clear which F10.7 indices ("adjusted" or "observed") were used in developing a specific IRI sub-model.In IRI-2016 the observed flux is used for NRLMSIS00 (computation of neutral temperature) and IRI-FLIP (computation of ion composition in the bottomside ionosphere).The difference of observed minus adjusted flux varies from −7 flux units in December to +7 units in June.
The Ionosonde Global (IG) index is obtained by adjusting the solar cycle variation of the CCIR (1967) foF2 model to ionosonde measurements of foF2.IG was introduced by Liu et al. (1983) and is being produced and distributed by the UK Solar System Data Centre (https://www.ukssdc.ac.uk/) with values going back to 1932.The CCIR model provides monthly averages of foF2 and accordingly IG is a monthly index and not a daily index.Indices are calculated at each ionosonde station and then IG is obtained as the global median.Originally 13 stations were used but the current index is obtained with only 4 stations (Chilton/U.K., Port Stanley/U.K., Kokubunji/Japan, and Canberra/Australia). IG has the advantage over solar indices that it is obtained from a ionospheric measurement and so includes dynamic effects not covered by a solar index.Disadvantages are (i) the averaging over several globally distributed station thereby smoothing out any regional differences; (ii) the change in number of stations over time affecting regional representation; (iii) only noontime data are used to produce the index; (iv) should be only used with the CCIR (1967) foF2 model not with other foF2 models.Because the global index smooths out regional differences it was found that foF2 is more closely correlated to the 12-month running mean of IG (IG 12 )than to the monthly IG and therefore IG 12 is used with CCIR-foF2 in IRI. Brown et al. (2018) found that with hemispheric IG indices the monthly indices perform better than the 12-month running means.They showed that use of these monthly hemispheric indices improves IRI model performance by almost a factor of 2. This is largely due to the fact that there are significant hemispheric differences between the station-specific IG values that are smoothed out with a global index.Themens and Jayachandran (2016) underscore this point in their study of IRI performance in the polar cap and subauroral region noting that the use of a monthly IG index, in place of the IRI's standard IG 12 index, leads to an improvement in TEC specification by up to 6 TECU.
In IRI these indices can be obtained from an internal file based on the date or the user can enter his/her own indices values for R 12 , IG 12 , F10.7_daily, and/or F10.7_81.Following user request starting with IRI-2016 the code now automatically adjusts IG 12 if only R 12 is provided by the user and vice versa using the approximate functional relationship between the two indices IG 12 = (−0.0031•R12 +1.5332)•R12 − 11.5634.For F10.7 all F107 indices (F10.7_daily,F10.7_81, and F10.7_365) are set to the same user-specified value if only one of them is specified.This was done to facilitate solar activity studies because a user may not always know which solar index affects which parameter.

Equivalent indices
An equivalent index is obtained by varying the model-driving index or indices until agreement is achieved between data and model.This is a method widely used for station-specific modelling.For regional or global modelling an averaging or interpolation process has to be applied to cover the region in between stations.Different data sources have been used for this purpose in conjunction with IRI.Most efforts have concentrated on obtaining equivalent R 12 or IG 12 indices because of the dominant role these two indices play in determining the F2 peak height and density, respectively.Bilitza et al. (1997) used worldwide ionosonde data from 1986 to 1989 to obtain equivalent IG indices with the goal of improving IRI ionospheric corrections for the altimeter measurements aboard the Geosat satellite.Komjathy et al. (1998) described a method of updating IRI with GPS total electron content (TEC) maps based on the grid of equivalent IG values corresponding to the TEC grid.Hernandez-Pajares et al. (2002) illustrated how the equivalent IG method can be applied directly to GPS slant TEC measurements.More recently this method has been expanded by using several data sources and/or by applying it iteratively to the IG and R indices.Migoya-Orué et al. (2015) adjusted the F2 peak parameter with ionosonde data and the topside profile with GPS data and Habarulema and Ssessanga (2016) (1998) did for the IG 12 variation in IRI.The results show that this is a very promising direction and future improvements are expected for IRI with the fine-tuning of these methods.

Local Updates with measured characteristic parameters
A direct update of IRI with measurements is possible if a user has data for the location in which s/he is interested in.
The IRI electron density profile is normalized to the F2, F1, and E peak parameters, so updating the IRI model with the input of these parameters is straight forward and includes the peak densities NmF2, NmF1, and NmE and the peak heights hmF2 and hmE.Instead of the F2, F1, and E peak densities (NmF2, NmF1, and NmE) a user can provide the corresponding plasma frequencies (foF2, foF1, and foE) and instead of the F2 peak height hmF2 one can enter the propagation factor M(3000)F2 that is inversely correlated to hmF2.This updating method has been widely used because routine ionosonde measurements of foF2 and M(3000)F2 can with the help of IRI predict the whole electron density profile for the ionosonde location.New with IRI-2016 a user now has the option to also enter the parameters B0 and B1, if measurements for these parameters are available.B0 determines the thickness of the IRI bottomside electron density profile and B1 the shape of the bottomside profile.This new capability is of particular interest for Digisonde users because the Automatic Real-Time Ionogram Scaler with True height (ARTIST) ionogram autoscaling software also determines B0 and B1 from the measured ionogram.In fact, the IRI Real-Time Assimilative Mapping (IRTAM) IRTAM makes use of this new option to assimilate measured B0 and B1 values into the background IRI (see next chapter).It should be pointed out that because of the special character of these parameters and to avoid unrealistic profiles B1 is constrained to the range 0.6-6 and B1 input is only accepted in conjunction with B0 input.

Data assimilation into IRI
Data assimilation into IRI and the development of a Real-Time IRI were discussed during several dedicated IRI Task Force meetings: 2009 Air Force Academy, Colorado Springs, Colorado, USA; 2012 Institute for Atmospheric Physics, Prague, Czech Republic; 2014 University of Massachusetts, Lowell, USA.A wide range of data assimilation techniques have been applied to IRI with many different data sources both in real-time and in retrospective mode.Schmidt et al. (2008) and Wenjing et al. (2015) represent the difference between GPS data and IRI with a multidimensional expansion in B-spline functions and have applied their technique successfully with GPS data from the South-American sector.Other examples are the GPS Ionospheric Inversion (GPSII) model of Fridman et al. (2006) using a Tikhonov methodology and the Electron Density Assimilative Model (EDAM) approach of Angling et al. (2009) based on a weighted, damped least mean squares estimation with IRI and GPS data.One of the most extensive efforts was undertaken by Yue et al. ( 2012) using a Kalman filter technique to assimilate GPS data, radio occultation data (CHAMP, GRACE, COSMIC, SAC-C, Metop-A, and TerraSAR-X), and Jason-1 and 2 altimeter TEC measurements into IRI.They did this for the time period 2002-2011 thus creating a valuable data source for space weather studies.Great progress towards a Real-Time IRI was achieved by Galkin et al. (2012) with the IRI Real-Time Assimilative Mapping (IRTAM, http: //giro.uml.edu/RTAM, last access: 20 July 2018).IRTAM is based on foF2 measurements from the worldwide network of Digisonde stations (the Global Ionospheric Radio Observatory -GIRO) and uses the CCIR (1967) model functions to represent the difference between data and model.A neural network interpolator is used for regions not covered by the GIRO network and IRTAM model coefficients for foF2 are generated every 15 min.The IRTAM data flow is illustrated in Fig. 6.The IRTAM project has made good progress and the foF2 modelling approach is now also applied to the F2 peak height hmF2 and the bottomside profile parameters B0 and B1, parameters that are routinely measured by the GIRO Digisondes.For these parameters the first step was the development of a CCIR-type model because CCIR (1967) only provides foF2.For hmF2 this work had been done already by Brunini et al. (2013).Another recent improvement is the inclusion of a linear trend term so the model parameter does not have to return to the same value after 24 h.
The latest version of IRI now also includes the necessary subroutines and changes to the computer code to facilitate the use of the IRTAM coefficients in the standard IRI code.

IRI successes
IRI is recognized and recommended as the standard for the ionosphere by almost all of the international bodies that in one way or other deal with the ionosphere.COSPAR and URSI, the two international unions that had initiated the IRI project, have both recognized IRI as the international standard for the ionosphere and recommended its usage to its member states.The European Cooperation for Space Standardization (ECSS) has elected IRI as its model of choice for the ionosphere.But the most important endorsement came with the election of IRI to become the ISO standard for the ionosphere in April 2014 (ISO 16457, 2014).ISO is the International Standardization Organization and ISO standards undergo a long vetting process whereby each country asks its ionospheric experts to review the standard proposal and to comment regarding the adoptions of the proposed standard.Final voting by country then decides about the acceptance or rejection of a standard.
IRI did very well in model assessment studies performed by independent/impartial assessors.Shim et al. (2011Shim et al. ( , 2012) ) undertook such an independent assessment of different models in the framework of the CEDAR Electrodynamics Thermosphere Ionosphere (ETI) Challenge for systematic assessment of ionosphere/thermosphere models; the Coupling, Energetics, and Dynamics of Atmospheric Regions (CEDAR) Program is sponsored by the National Science Foundation (NSF).Of the 8 models included in the study, IRI was the only empirical model.The others were theoretical models and two of them involved assimilation of event data.The data included NmF2 and hmF2 from ionosondes and COSMIC radio occultation and electron density along the CHAMP orbit.The assessment was done for 9 event periods: 2 strong, 4 moderate storm events, and 3 quiet periods.Metrics used included RMS, and the model-data ratio of maximum and minimum values and maximum-minimum range.IRI was the clear overall winner.It was, if not first than second or third for almost all events and metrics parameters.Shim et al. (2017)  Table 2 shows the statistics regarding the usage of IRI in papers published in the most prominent science journals cov-   1 in Bilitza et al. (2017).The numbers speak for themselves, per year IRI contributes to 5 %-8 % of JGR papers and even more impressive to 10 %-14 % of RS papers made use of IRI for their study.Figure 7, finally, displays the monthly accesses to the IRIweb online system for calculating IRI parameters.Accesses reached the 1 000 000 mark in 2014 and have continued to increase since then.

Summary and outlook
The International Reference Ionosphere (IRI) is the international standard for Earth's ionosphere and is recognized as such by COSPAR and URSI and most importantly by the organization in charge of international standards the International Standardization Organization (ISO).In April 2014 IRI was elected to become the ISO standard 16457 by unanimous vote by the ISO member states.This article focusses on the external drivers of the IRI model and discusses some of the changes introduced with IRI-2016.Besides the traditional solar, magnetic, and ionospheric indices, IRI can be also driven by characteristic profile parameters either by direct input if data are available locally or by model adjustment and data assimilation for regional or global modelling.Significant progress has been achieved in the area of data assimilation.Yue et al. ( 2012) have produced a valuable data set for the years 2002-2011 by assimilating a large amount of GNSS and altimeter data into IRI.Galkin et al. (2012) continue to improve their IRI Real-Time Assimilative Mapping (IRTAM) now including the parameters foF2, hmF2, B0, and B1.IRI has been widely used as documented by the large number of citations (4236) since 1990 in a wide range of journals (25) and the significant percentage of papers each year in the Journal of Geophysical Research (5 %-8 %) and Radio Science (10 %-14 %) that have used the IRI model in their scientific studies.Usage of the online IRI interface has now reached the level of 8 000 000 accesses per month.
Future improvements of the IRI model are expected in a number of areas.Several new models are already scheduled for inclusion in the next release of the IRI model.This includes a description of storm effects on hmF2 (Blanch and Altadill, 2012), a scintillation occurrence probability model (Chen et al., 2017).Other efforts are close to completion including an updated model for the D-region electron density based on a simple ion-chemical model that is corrected with rocket radio propagation data (Martin Friedrich, Graz University of Technology, Austria) and the Vary-Chap formalism for the representation of the topside electron density profile using Alouette and ISIS topside sounder data and the TOP-IST software (Bodo Reinisch, Ivan Galkin, and Osman Mohamed, University of Massachusetts Lowell, Lowell, USA).Of special importance is a better description of IRI densities and temperatures at very low solar activities in view of the discrepancies found during the last unusually low and extended solar minimum (2007)(2008)(2009) and because the current solar cycle may reach again very low solar activities.The topside electron density, in particular, has been shown to be overestimated by IRI in comparisons with C/NOFS, CHAMP, and GRACE satellite measurements (Lühr and Xiong, 2010;Bilitza et al., 2012).An effort to correct this problem is nearing completion using also Swarm in situ measurements and Alouette and ISIS topside sounder data (Di-eter Bilitza with Chao Xiong, Geo Research Centre, GFZ, Potsdam, Germany).Intense efforts are underway to enable the use of GNSS data in the Real-Time IRI (Ivan Galkin, University Massachusetts Lowell, Lowell, USA and Andrzej Krankowski, University of Warmia and Mazury, Olsztyn, Poland).As pointed out in this article the use of monthly hemispheric ionospheric indices may lead to a significant performance increase for foF2 model predictions.This will require setting up a routine computation of these indices.A model for the occurrence probability of sporadic E based on COSMIC radio occultation data is in the very early stages of development (Dieter Bilitza with Christina Arras, Geo Research Centre, GFZ, Potsdam, Germany).These are just a few examples of studies that are ongoing with the goal of making the IRI model an even better standard for the Earth's ionosphere.
Competing interests.The author declares that there is no conflict of interest.
Special issue statement.This article is part of the special issue "Kleinheubacher Berichte 2017".It is a result of the Kleinheubacher Tagung 2017, Miltenberg, Germany, 25-27 September 2017.

Figure 1 .
Figure1.Brief history of the IRI project with the most important milestones.Note: IRI-79 is also described inRawer et al. (1978b).

Figure 2 .
Figure 2. The data sources used for the development of IRI with depictions of the instruments that accumulated these data sources: ionosondes (a), incoherent scatter radars (b), rockets (c), topside sounder and in situ satellites (d), Global Navigation Satellite Systems (e), and radio occultation (f).Data sources highlighted in green have been used for validation but not yet for model development.

Figure 3 .
Figure 3.Comparison of simultaneous measurements of electron temperature over Millstone Hill at 840 km with the Millstone Hill incoherent scatter radar (Te_MH) and with the DMSP Langmuir Probe instrument (Te_DMSP) separated into daytime and nighttime values.Also show is the Te_MH = Te_DMSP line (solid line) and the least-square fits to the day (dotted lines) and night values (dashed lines).

Figure 6 .
Figure 6.Illustration of the data flow in the IRI Real-Time Assimilative Mapping (IRTAM) system from the Global Ionospheric Radio Observatory (GIRO) network of digisondes to the assimilation into the IRI electron density profile every 15 min (http://giro.uml.edu/RTAM, last access: 20 July 2018).
continued their model assessment effort with the study "CEDAR-GEM Challenge for Systematic Assessment of Ionosphere/Thermosphere Models in Predicting TEC during the 2006 December Storm Event".The event was the storm on 14-15 December 2006 (AGU storm event) and the data used in this case were GPS total electron content (TEC) measurements in eight selected longitude sectors.The result show IRI on par with the best theoretical models.IRI is highly cited in a wide range of scientific journals.Since 1990 the most important IRI papers have been cited 4236 times in 25 different journals including Advances in Space Research, Journal of Asian Earth Science, Advances in Radio Science, Journal of Atmospheric and Solar-Terrestrial Physics, Annales Geophysicae, Applied Optics, Journal of Geodesy, Astrophysics and Space Science, Journal of Geophysical Research, Chinese Journal of Aeronautica, Journal Space Weather & Space Climate, Computer Physics Communications, Planetary and Space Science, Computers & Geosciences, Plasma Science and Technology, Cosmic Research, Radio Science, Geochimica et Cosmochimica Acta, Solar Physics, Geodesy and Geodynamics, Space Science Review, Geophysical Research Letters, Space Weather, GPS Solutions, and Surveys in Geophysics.

Table 1 .
Indices used in the different IRI sub-models.

Table 2 .
Percentage of papers per year that acknowledge usage of the IRI model in the AGU journals Journal of Geophysical Research (JGR), Geophysical Research Letters (GRL), Space Weather (SW), and Radio Science (RS) for the years 2009 to 2017.ionospheric physics including the Journal of Geophysical Research (JGR), the Geophysical Research Letter, the Space Weather journal, and the journal Radio Science (RS).We first search for citations of IRI papers and then determine if IRI actually contributed to the paper or was just referenced.Please note that this extends and corrects a similar table shown as Table ering