Asian summer monsoon anticyclone: trends and variability

. The Asian summer monsoon anticyclone (ASMA) has been a topic of intensive research in recent times re-garding its variability in dynamics, chemistry and radiation. This work explores the spatial variability and the trends of the ASMA using observational and reanalysis data sets. Our analysis indicates that the spatial extent and magnitude of the ASMA is greater during July and August than in June and September. The decadal variability of the anticyclone is very large at the edges of the anticyclone compared with the core region. Signiﬁcant decadal variability is observed in the northeastern and southwestern parts of the ASMA with reference to the 1951–1960 period. The strength of the ASMA shows a drastic increase in zonal wind anomalies in terms of temporal variation. Furthermore, our results show that the extent of the anticyclone is greater during the active phase of the monsoon, strong monsoon years, and La Niña events. Signiﬁcant warming with strong westerlies is observed exactly over the Tibetan Plateau from the surface to the tropopause during the abovementioned periods. Our results support the existence of transport process over the Tibetan Plateau and the Indian region during active, strong monsoon years and during strong La Niña years. Therefore, it is recommended that the different phases of the monsoon be taken into account when interpreting the variability of pollutants and trace gases in the anticyclone.


Introduction
The Asian summer monsoon anticyclone (ASMA) is a dominant circulation in the Northern Hemisphere (NH) summer that extends from Asia to the Middle East in the upper troposphere and lower stratosphere (UTLS). The ASMA is bordered by the subtropical westerly jet in the north and easterly jets in the south. The Asian summer monsoon (ASM) dynamics act as a pathway for the transport of trace gases and pollutants both vertically (through convection) and horizontally (via the low-level jet and tropical easterly jet). The ASMA circulation responds to heating corresponding to the deep convection of the south Asian monsoon (Hoskins and Rodwell, 1995;Highwood and Hoskins, 1998). This strong anticyclone circulation isolates the air and is tied to the outflow of deep convection, which has distinct maximum characteristics in terms of dynamic and chemical variability (Randel and Park, 2006;Park et al., 2007). Recently, the anticyclone circulation in UTLS has been paid more attention by researchers in order to understand the dynamics, chemistry and radiation of the region. This problem has been discussed by several authors (e.g., Park et al., 2007;Fadnavis et al., 2014;Glatthor et al., 2015;Vernier et al., 2015;Santee et al., 2017). Deep convection during the monsoon can transport tropospheric tracers from the surface to the UTLS (Vogel et al., 2015;Tissier and Legras, 2016). The tracers that are transported are confined in the anticyclone and, consequently, affect the trace gas concentrations in the UTLS, resulting in significant changes in radiative forcing (Solomon et al., 2010;Riese et al., 2012;Hossaini et al., 2015). The center of the anticyclone is located either over the Iranian Plateau or over the Tibetan Plateau, where the distribution of pollutants and tracers vary significantly (Yan et al., 2011).
The spatial extent, strength and location of the anticyclone vary on several temporal scales due to the internal dynamic variability of the Asian monsoon (Zhang et al., 2002;Randel and Park, 2006;Garny and Randel, 2013;Vogel et al., 2015;Pan et al., 2016). However, the variability of the anticyclone structure and the response to the Indian monsoon activity are not understood and, consequently, neither is the variability of the tracers (e.g., O 3 and CO) trapped in the anticyclone during the same period in the UTLS region. Since the anticyclone extends from the Middle East to East Asia, trapped pollutants are expected to make a large radiative forcing contribution to the background atmosphere. Thus, it is essential to understand the variability of anticyclone structure itself in detail as well as its response to the Indian summer monsoon (ISM). Therefore, in the first part of the study, we investigate the spatial, interannual and decadal variations of the anticyclone. Since the Indian monsoon responds at different timescales, we also investigate the anticyclone variability with respect to the active and break phases of the Indian monsoon, strong and weak monsoon years, and the stronger El Niño-Southern Oscillation (ENSO) years. For this, we utilize the National Center for Environmental Prediction/National Center for Atmospheric Research (NCEP/NCAR) reanalysis geopotential height from 1951 to 2016. The structure of the paper is as follows: the data sets used in this study are described in Sect. 2; Sect. 3 contains the seasonal and decadal variation of the anticyclone; Sect. 4 shows the influence of the ISM on the anticyclone, i.e., active and break phases, strong and weak monsoon years, and ENSO's effects on the anticyclone; and, finally, the results are presented in Sect. 5.

NCEP/NCAR reanalysis
The National Center for Environmental Prediction (NCEP) in collaboration with the National Center for Atmospheric Research (NCAR) produces reanalysis data from a consistent assimilation and modeling procedure that incorporates all of the available observed conditions obtained from conventional and satellite information from 1951 to the present (Kalnay et al., 1996). We used NCEP/NCAR reanalysis daily geopotential height (GPH) and wind data from 1951 to 2016. The NCEP/NCAR data assimilation uses a 3D-variational analysis scheme with 28 pressure levels and triangular truncation of 62 waves (at a horizontal resolution of 200 m). Both GPH and temperature at the chosen standard levels are described as class output variables (Kalnay et al., 1996), i.e., they are strongly influenced by observed data. Only the Indian summer monsoon months (June, July, August and September) containing gridded daily data were considered in this study. The NCEP/NCAR reanalysis data have a spatial resolution of 2.5 • . The seasonal values are estimated from daily data. To identify the spatial and temporal variations of the anticyclone center, we used the monthly mean values of the GPH and the zonal wind component. The quality of NCEP GPH reanalysis data is discussed in Bromwich et al. (2007).

IMD gridded precipitation data
The India Meteorological Department (IMD) high-resolution (0.25 • × 0.25 • ) gridded precipitation data are used to identify the active and break phases during June, July and August from 1951 to 2016 (Pai et al., 2016). These precipitation data have been validated extensively with observational and reanalysis data sets and display a very good correlation (Kishore et al., 2016). We identified the active and break phases based on daily rainfall over the monsoon core zone of India (which roughly spans from 18 to 28 • N and from 65 to 88 • E) during July and August, as reported by Rajeevan et al. (2010). The normalized anomaly estimated from the averaged daily rainfall in the monsoon core zone is subtracted from its long-term  mean by dividing it by its daily standard deviation. The active (break) phases were identified from the normalized anomaly when rainfall is greater (less) than −1.0 (+1.0) for 3 consecutive days or more.

GNSS radio occultation (RO) data
We also used the Global Navigation Satellite System (GNSS) radio occultation (RO) data to investigate the temperature anomaly. The basic measurement principle of RO exploits the atmosphere-induced phase delay in the GNSS signals, which are recorded in the low Earth-orbiting satellite. This technique provides vertical profiles of refractivity, density, pressure, temperature and water vapor (Kursinski et al., 1997). The temperature profiles from this technique are available with low horizontal (∼ 200-300 km) and high vertical resolutions (10-35 km) with an accuracy of < 0.5 K. We used CHAllenging Minisatellite Payload (CHAMP) and Constellation Observing System for Meteorology, Ionosphere, and Climate (COSMIC) data covering the period from 2002 to 2016.
The CHAMP satellite was launched into a circular orbit by Germany on 15 July 2000 to measure the Earth's gravity and magnetic field and to provide global RO soundings (Wickert et al., 2001). About ∼ 230 RO profiles per day have been measured by the CHAMP payload since 2002. The CHAMP payload was solely designed to track the setting occultations, and the RO event is terminated when the signal is lost -this results in a decrease in the number of occultations with decreasing altitude (Beyerle et al., 2006). This receiver measures the phase delay of radio wave signals that are occulted by the Earth's atmosphere. From this phase delay, it is possible to retrieve the bending angle and refractivity vertical profiles.
COSMIC consists of a constellation of six satellites, which were launched into a circular, 72 • inclination orbit at a 512 km altitude in April 2006 and are capable of receiving signals from the Global Positioning System (GPS; Anthes et al., 2008). Compared with previous satellites, COS-MIC satellites employ an open-loop mode, which can track both the rising and setting of occultations (Schreiner et al., 2007). The open-loop tracking technique significantly reduces the GPS RO inversion biases by eliminating tracking errors (Sokolovskiy et al., 2006). The COSMIC temperature profiles display a very good agreement with radiosonde data, reanalyses and models (Rao et al., 2009;Kishore et al., 2011Kishore et al., , 2016. The CHAMP and COSMIC GPS RO data were interpolated to 200 m from their native resolution. We derived the cold-point tropopause altitude/temperature over the ASMA region, as discussed by Ratnam et al. (2014) and Ravindra Babu et al. (2015). Both the CHAMP and COSMIC data were obtained from the COSMIC Data Analysis and Archive Center (CDAAC; https://cdaac-www.cosmic.ucar.edu/cdaac/ products.html, last access: 30 June 2019).

Variability of the anticyclone
The climatological spatial variability of the GPH and wind vectors at 100 hPa during June, July, August and September from NCEP reanalysis data is shown in Fig. 1a-d. The anticyclone circulation is clearly depicted by wind vectors during the abovementioned months ( Fig. 1). During the months of September and June, the GPH values, which represent the spatial extent of the anticyclone, are low compared with July and August. Thus, the spatial extent and intensity of the anticyclone are greater during July than in the other months. During July and August, the anticyclone extends from East Asia to the Middle East. The spatial extent of the anticyclone circulation is clearly evident in the grid from 15 to 45 • N and from 30 to 120 • E at 100 hPa, and the climatological averaged GPH values vary from 16.5 to 17 km in the NCEP reanalysis between 1951 and 2016. Using the modified potential vorticity equation, Randel et al. (2006) showed the spatial variation of the anticyclone where GPH values are stationary in the range from 16.75 to 16.9 km. Similarly, Park et al. (2007) showed the anticyclone structure from the strongest wind at 100 hPa using the streamline function. Bian et al. (2012) reported the spatial variability of the anticyclone using 16.77 and 16.90 km in the GPH contour as the lower and the upper boundaries, respectively; thus, these empirically selected GPH values represent the anticyclone boundaries. Therefore, in this present study, we have chosen the values from 16.75 to 16.9 km to investigate the spatial features of the anticyclone, and the results are depicted in Fig. 1e-h. The spatial extent and existence of the anticyclone are highly prominent during July and August compared with June. During September, very low GPH values are seen compared with July and August. Therefore, we considered the average of the July and August GPH values from 1951 to 2016 for further analysis, as shown in Fig. S1 in the Supplement. The core region and the spatial extent of the anticyclone are clearly evident in Fig. S1. The core region of the anticyclone shows a bimodal distribution, i.e., one core located at the southwestern flank of the Himalayas and another over Iran. The core region over the southwestern flank of the Himalayas is due to large-scale updraft, which is caused by the moist energy over the Indo-Gangetic Plain, heating of the Tibetan Plateau and the orographic forcing of the Himalayas. Severe heating over the Arabian Peninsula supports the formation of the mid-tropospheric anticyclone in the west. This anticyclone can merge intermittently with the ASMA. It is also observed that the spatial extent of the anticyclone varies drastically at different temporal scales. Therefore, seasonal variation is much more pronounced.
The decadal variation of the anticyclone is studied with respect to the spatial variability. Figure 2 shows the decadal spatial variation of the anticyclone with reference to the 1951-1960 period. The significant difference in the decadal variation is noticed in Fig. 2. The edges (east, north and west) of the anticyclone undergo drastic changes during the 1961-1970 period. With respect to the 1971-1980 period, except for a small portion in the east, the whole anticyclone shows drastic changes. During the decade from 1971 to 1980, the recorded GPH values in the anticyclone are also ∼ 25 m lower compared with the values from 1951 to 1960. The inverse is seen from 1981 to 1990, as high GPH values (∼ 30 m) are observed compared with those from the reference period. The GPH difference is significant over the western, northeastern and southern regions of the anticyclone during the 1991-2000 period. Similar changes are observed from 2001 to 2010. Compared with all of the decadal differences, 2011-2016 shows a completely different picture: changes are only seen in the western and northeastern corners, whereas other parts of the anticyclone do not show any change. From this analysis, we observed significant changes in the anticyclone even from one decade to another, which can result in a change in chemical and dynamic variability over this region.
Furthermore, the spatial distribution of the trend is estimated from 1951 to 2016 using robust regression analysis at a 95 % confidence interval, as displayed in Fig. 3. The edges on all sides of the anticyclone undergo noticeable changes compared with the core region. The eastern and northwestern sides of the anticyclone show an increasing trend compared with the other regions. The trends at the northern end are more significant than at the southern end. A few areas on the northern side of the anticyclone show a reduction in strength. Therefore, in order to understand the asymmetry in the anticyclone variability, we divided the anticyclone region into four different sectors, as shown in  It is well known that the Indian monsoon rainfall varies at different timescales, i.e., daily, sub-seasonal, interannual, decadal and centennial scales (Rajeevan et al., 2010). Precipitation during the monsoon varies from intra-seasonal scales between active (good rainfall) and break (less rainfall) phases. Any small change in the precipitation pattern will affect the anticyclone due to the thermodynamics in-volved in rainfall. In this study, we also investigated the anticyclone variability during the active and break phases of the Indian monsoon. The active and break periods were identified in July and August using the high-resolution gridded (0.25 • × 0.25 • ) rainfall data from 1951 to 2016 as defined by Pai et al. (2016).
The number of active and break days is derived from the precipitation data shown in Fig. S2a and b. Daily GPH, temperature and zonal wind are taken from the NCEP reanalysis, whereas the tropopause altitude is derived from the GNSS RO data for active and break days. The anticyclone structure during active (red line) and break (blue line) days is shown in Fig. 6a. Two interesting aspects of the anticyclone variability can be noticed between active and break days: (1) the extent of the anticyclone is large during active days compared with break days, and (2) there are two-cell structures in the anticyclone core region during active days. The extent of the anticyclone is large in the eastern and northern sectors on active days. The zonal (meridional) cross section of temperature (colored shading), zonal wind (contour lines) and the difference between the active and break phases averaged in the longitude band from 80 to 90 • E (latitude band from 30 to 40 • N) along with the cold-point tropopause for active and break days are illustrated in Fig. 6b and c. During     active days, the temperature shows cooling at tropical latitudes, whereas it shows warming in the midlatitudes from the surface to the tropopause. Significant warming is observed during the active days in the mid-troposphere over the Tibetan Plateau and its northern side. Westerly (easterly) winds exist over the cooler (warmer) regions. The warm temperature anomalies stretch from 1.5 to 12 km between 25 and 60 • N. The tropopause altitude is low (high) during the active (break) phases of the Indian monsoon, as shown in Fig. 6b. The meridional cross section of temperature anomalies displays significant warming from ∼ 1.5 to 8 km over the Indian region. The tropopause altitude exemplifies random variability in the meridional cross section.
As discussed previously, the anticyclone circulation is significant during the months of July and August when most of the precipitation occurs over India (Basha and Ratnam, 2013;Basha et al., 2015;Kishore et al., 2015;Narendra Reddy et al., 2018). Thus, the influence of strong and weak monsoon years will have a drastic impact on anticyclone circulation. In order to understand these changes, we have divided the years into strong and weak monsoon years based on gridded precipitation data over the domain from 5 to 30 • N and from 70 to 95 • E between 1951 and 2016. This region is known to have heavy precipitation and orographic forcing, which helps transport water vapor via deep convection to the UTLS (Houze et al., 2007;Medina et al., 2010;Pan et al., 2016). The detrended precipitation represents the strong and weak monsoon years. Years with positive (negative) precipitation values correspond to strong (weak) monsoon years as shown in Fig. S2b. The composite of the mean distribution of the anticyclone circulation during strong and weak monsoon years is shown in Fig. 7a based on GPH values at 100 hPa from NCEP reanalysis data. The circulation expands on the eastern and western sides of the anticyclone during the strong monsoon years (red line). The core of the anticyclone is significant during strong monsoon years. A clear eye structure is observed in the core of the anticyclone on left (right) during the strong (weak) monsoon years. The composite mean difference of temperature and zonal wind between the strong and weak monsoon years along with the tropopause altitude averaged over the longitude range from 80 to 85 • E is shown in Fig. 7b. The warmest temperature anomalies are observed over the Tibetan Plateau. Positive (warm) temperature anomalies exactly above the Tibetan Plateau (11 km) and negative (cooling) anomalies on both sides are seen in the lower troposphere in Fig. 7b. Strong easterly (west- erly) winds are observed on the left (right) side of the Tibetan Plateau. The whole Tibetan Plateau acts as a barrier that drives cold air to upper altitudes during strong monsoon years. Strong anticyclone circulation with strong westerlies at 35 • N and easterlies on both sides of the Tibetan Plateau along with an elevated tropopause represent the impacts of the strong monsoon vertically above the anticyclone. The rising motion over East Asia is excited by the local heating of the Tibetan Plateau and links to the single stretch vertically. The longitude-altitude cross section of temperature and wind anomalies shown in Fig. 7c is averaged between a latitude band from 35 to 40 • N. Positive temperature anomalies are observed from the surface to 12 km in the longitude band from 60 to 80 • E and stretch towards the west. This process clearly demonstrates that a large-scale ascent develops over the Asian monsoon region. The tropopause altitude is high (low) during strong vertical motion, and heavy precipitation is observed over a region similar to that reported by Lau et al. (2018). The transport processes from the boundary layer to the tropopause occur on the eastern side of the anticyclone, i.e., the southern flank of Tibetan Plateau, northeast India and the head of the Bay of Bengal. This result is consistent with a previous study by Bergman et al. (2013).
ENSO typically shows the strongest signal in boreal winter, but it can affect the atmospheric circulation and constituent distributions until the next fall (McPhaden et al., 2006). It is well known that strong ENSO events have a significant influence on tropical upwelling and stratospheretroposphere exchange (Yan et al., 2018). This change can impact the distribution of the composition and structure of the UTLS region. In the UTLS, the tropopause responds to the annual and interannual variability associated with ENSO (Trenberth, 1990) and the quasi-biennial oscillation (Baldwin et al., 2001). Several studies have focused on the effects of the different impacts of El Niño in the tropopause and lower stratosphere (Hu and Pan, 2009;Zubiaurre and Calvo, 2012;Xie et al., 2012). In the present study, we investigated the changes in the anticyclone circulation and tropical upwelling associated with strong ENSO events during July and August. Therefore, we have also separated the GPH for the strongest El Niño (1958, 1966, 1973, 1983, 1988, 1992, 1998) and La Niña (1974, 1976, 1989, 1999, 2000, 2008 years to verify the change in the circulation pattern of the anticyclone. For this, we chose July and August GPH data at 100 hPa, as shown in Fig. 8. The red and blue colors indicate the composite of the La Niña and El Niño circulation, respectively. During La Niña years, the anticyclone circulation expands compared with El Niño years at 100 hPa, as shown in Fig. 8a. On the eastern and southern sides of the anticyclone, the expansion is greater during La Niña years. Warm temperatures and strong westerlies are observed in the latitude band from 43 to 55 • N during La Niña, as shown in Fig. 8b (Lau et al., 2018). The cooling impact is significant over the Tibetan Plateau during La Niña events compared with El Niño events. The significant cooling observed over the Tibetan Plateau distributes towards tropical latitudes between 600 and 100 hPa. The zonal wind shows a convergence of easterly winds over the Tibetan Plateau from the mid to the upper tropospheric region. On the right side of the Tibetan Plateau, strong westerly winds exist from the surface to tropopause altitudes along with strong warming. The meridional cross section of temperature and the zonal wind difference between La Niña and El Niño is shown in Fig. 8c. Significant cooling is observed during La Niña in the longitude band from 80 to 100 • E with strong easterlies from the surface to the tropopause. From this analysis, it is clear that the Indian summer monsoon variability has a significant impact on the ASMA; thus, it is necessary to consider the different phases of the monsoon while dealing with UTLS pol-lutants. In addition, we investigated the zonal mean vertical cross section in the longitude band from 50 to 60 • E, which represents the Iranian mode. Figure S3 depicts the difference between the active and break phases, strong and weak monsoon years, and La Niña and El Niño years along with the tropopause altitude. Significant warming is observed during La Niña years and strong monsoon years compared with the active phase of the Indian monsoon in the troposphere. Compared with the Tibetan mode, the Iranian mode warming is lower. The tropopause altitude is slightly higher during the active phase of the Indian monsoon, strong monsoon years and La Niña years. A moderate increase in the tropopause from the Equator to 40 • N is observed, and it decreases drastically afterward.

Summary and conclusions
Several authors have discussed the interannual and decadal variability of pollutants and tracers in the ASMA region from model, observational and reanalysis data sets (Kunze et al., 2016;Santee et al., 2017;Yuan et al., 2019). In this study, we have investigated the spatial variability, trends of the anticyclone and the influence of Indian monsoon activity, i.e., active and break days, strong and weak monsoon years, and strong La Niña and El Niño years, on the ASMA using longterm reanalysis, satellite and observational data sets that were not investigated earlier. In this study, we have considered the GPH values from 16.75 to 16.9 km, which represent the spatial structure of the anticyclone at 100 hPa. Our analysis shows that the spatial extent (magnitude) of the anticyclone structure is very large (strong) during July and August, whereas it is very weak in June at 100 hPa. The bimodal distribution (Tibetan and Iranian modes) of the anticyclone is clearly observed during the month of July, but it is absent during other months (June and August). The anticyclone variability undergoes significant decadal variation from one decade to another. The edges of the ASMA change drastically compared with the core of the anticyclone. However, there are significant spatial differences in the structure of the anticyclone at 100 hPa. The anticyclone undergoes a decreasing trend on the northern side, whereas an increasing trend is observed on the western side. A significant increasing trend is observed in the spatially averaged zonal wind in the four different sectors (Fig. 5) Furthermore, we investigated the Indian monsoon influence on the anticyclone region. Our results reveal that the spatial extent of the anticyclone expands during the active phase of the Indian monsoon, strong monsoon years and during strong La Niña years on the northern and eastern sides. A similar expansion of the anticyclone was noticed by Yuan et al. (2019) during strong monsoon years when considering MERRA 2 data. However, the ASMA boundaries are not always well defined in all the events. The zonal mean cross section of temperature shows significant warming over the Tibetan Plateau and from the surface to 12 km during the active phase of the Indian monsoon, strong monsoon years and strong La Niña years. Similarly, the rise of the tropopause during the abovementioned periods is also noticed. Since the Tibetan Plateau acts as a strong heat source in summer with the strongest heating layer located in the lower layers, thermal adaptation results in a shallow and weak cyclonic circulation near the surface and a deep and strong anticyclonic circulation above it. During summer, the Tibetan Plateau acts as a strong heat source that influences the whole UTLS region. The warm ascending air above pulls the air from below; the surrounding air in the lower troposphere converges towards the Tibetan Plateau area and climbs up the heated sloping surfaces (Bergman et al., 2013;Garny and Randel, 2016). Thus, significant warming is observed over the Tibetan Plateau and causes the strong transport of pollutants into the tropopause during the active phase of the Indian monsoon, the strong monsoon years and the strong La Niña years. Pan et al. (2016) reported the transport of carbon monoxide through the southern flank of the Tibetan Plateau from the model analysis. The abovementioned results indicate that the high mountain regions play a significant role as elevated heat sources during the formation and maintenance of the anticyclones over Asia. This also emphasizes the role of the thermal forcing of the Tibetan Plateau on the temporal and the spatial evolution of the South Asian high. Lau et al. (2018) reported the transport of the dust and pollutants from the Himalayas-Gangetic Plain and the Sichuan Basin.
Overall, we demonstrate the ASMA variability during different phases of the Indian monsoon. The uplifting of boundary layer pollutants to the tropopause primarily occurs on the eastern side of the anticyclone, centered near the southern flank of the Tibetan Plateau, northeastern India, Nepal and north of the Bay of Bengal. The variability in the tropopause altitude and temperature as well as trace gases (water vapor, ozone, carbon monoxide and aerosols) shows a distinct behavior in ASMA region. The ASMA itself is highly dynamic in nature and the confinement of tracers and aerosols results in changes in its chemistry and radiation . However, a more detailed and higher quality data set is needed to further understand the effects of the Tibetan Plateau on the transport of different tracers and pollutants to the UTLS region .
Data availability. The NCEP/NCAR reanalysis data are available from the NOAA website: https://www.esrl.noaa.gov/psd/ data/gridded/data.ncep.reanalysis.pressure.html (last access: 30 June 2019; Kalnay et al., 1996). The COSMIC and CHAMP data are available from the COSMIC CDAAC website. IMD gridded precipitation data are available from the National Climate Data Management Group, Pune, India. All the data used in the present study are freely available from the respective websites.
Author contributions. GB and MVR conceived and designed the scientific questions investigated in the study. PK estimated the active and break phases of the Indian monsoon. GB performed the analysis and wrote the first draft of the paper in close association with MVR. All authors edited the paper.
Competing interests. The authors declare that they have no conflict of interest.
Special issue statement. This article is part of the special issue "Interactions between aerosols and the South West Asian monsoon". It is not associated with a conference.