The effect of economic institutional based community, social network establishment, re-trust of social affiliation, and regrouping based on capital goods on economic welfare

Abstract The purpose of this paper is to determine the effect of the Economic Institutional Based Community, Social Network Establishment, Re-trust of Social Affiliates, and Regrouping Based on Capital Goods on Economic Welfare through the Moderation of Rubber Sector Coordination, and Collaboration on the Rubber Sector. This research utilized qualitative explanatory associative. The study was conducted on the rubber plantation community in Kuantan Singingi Regency, Riau Province, Indonesia. The research approach used was quantitative with an analysis tool Generalized Structure Component Analysis (GSCA). The results of the analysis, there is an influence of the Economic Institutional Based Community, Social Network Establishment, Re-trust of Social Affiliates, and Regrouping Based on Capital Goods and Rubber Sector Coordination and Collaboration Rubber Sector on Economic Welfare. The higher value of the Economic Institutional Based Community, Social Network Establishment, Re-trust of Social Affiliates, and Regrouping Based on Capital Goods and Rubber Sector Coordination and Collaboration on the Rubber Sector will have an impact on the increasing economic welfare. In addition, it was also found that the Rubber Sector Coordination and Collaboration moderated the influence of the Economic Institutional Based Community variable, Social Network Establishment, Re-trust of Social Affiliates, and quasi Regrouping Based on Capital Goods moderator and enrichment. The existence of a moderating effect in this study is a distinction between this research and previous studies. Previous research merely described the problem qualitatively. Similar research has never been carried out in the Kuantan Singingi area.

PUBLIC INTEREST STATEMENT
Indonesia is an agricultural country because most of its population has a livelihood in the agricultural sector, one of which is in the production of rubber. Riau is one of the provinces in Indonesia that has the largest rubber production. However, over time, rubber production in Riau Province has stagnated. This is partly due to the lack of empowerment of human resources, especially in small and medium social groups. Therefore, in this study, the development of small and medium-scale socio-economic groups in rubber plantations is the main target that in turn will provide opportunities for the development of economic activities.

Introduction
Rubber production in Riau Province is one of the largest in Sumatra. In the period 2008 to 2012, rubber production in Riau Province experienced stagnancy. The production rate was in the range of 325,000 tons to 396,000 tons in 500,949 hectares of land use (Statistics, 2013). Based on land use, the highest rubber commodity development potency is in Kuantan Singingi Regency possessing an area of 146,215 hectares. The second largest is Kampar Regency with an area of 101,597 hectares, and Indragiri Hulu Regency with an area of 61,372 hectares (Statistics, 2013). Furthermore, the rubber plantations' stagnant productivity is counterproductive to macro development scenarios compiled at both the central and regional levels. Master Plan for the Acceleration of Indonesian Economic Development (MP3EI) explicitly stated that there is a need for accelerated economic development in Sumatra through the rubber industry development. The output is to build\ a downstream industry for rubber commodities, which will eventually bring multiplier effects in other fields, especially the community economy.
The urgency to revitalize rubber commodities in Kuantan Singingi Regency increased when people who had been concerned about managing rubber plantations eventually moved to palm oil management as it is considered more profitable in the short/medium term. The decline in the community interest to manage rubber indicates a decline in the level of pull and push factors in rubber commodities development in the Kuantan Singingi Regency. Hence, the solution to this issue needs to be immediately sought.
Human resources on rubber commodities management in Kuantan Singingi Regency is concentrated in self-management management patterns. The self-reliant (Swadaya) pattern tends to position the farmer as subordinate. For example, in terms of the rubber selling price. The other four patterns have more bargain and bargaining power in the context of managing rubber commodities, as these patterns move their business in more managerial and systematic ways compared to selfreliant patterns. For example, with the awareness to form an institution. The latest data exhibited that out of 62,975 farmers, 35 rubber plantation business institutions have been formed with a total membership of around 1,012. There are a huge number of farmers managing rubber plantations. Nevertheless, these have not been developed and/or empowered according to their strengths.
The development of small and medium scale socio-economic groups on rubber plantations needs to be the main target in community-based development activities. Such community development will, in turn, provide opportunities for the development of local economic activities and other productive endeavors (Nasdian, 2014). Nasdian further argued that institutional development is crucial for the development of productive economic efforts as this development requires high transaction costs. Also, through institutional development, especially customary institutions, it was found to be able to overcome various social problems, such as child neglect and poverty (Hikmat, 2006) The development of rubber plantation institutions reflects the extent to which social networks have been built into these communities. When these two aspects (institutions and social networks) are synergized with various stakeholders functions of the rubber plantation business, will manifest into the development of social capital. Sujianto's (2008) study exhibited that social capital can provide community development and empowerment on the economic and social aspects through existing social networks. It cannot be denied that community development based on social capital can create a prosperous and empowered society (NurhNurhadi, 2008). The challenge is how to develop and utilize social capital. Identifying whether social capital possesses potency or not is important because it can not only support the ongoing development process but may weaken the development process itself (Adi, 2008).
Based on the description, a study was conducted which aims to determine the effect of the Economic Institutional Based Community (X1), Social Network Establishment (X2), Re-trust of Social Affiliates (X3), and Regrouping Based on Capital Goods (X4) on Economic Welfare (Y) through the Moderation of Rubber Sector Coordination (M1), and Collaboration on the Rubber Sector (M2). The existence of a moderating effect in this study is a distinction between this research and previous studies. Previous research merely described the problem qualitatively. Similar research has never been carried out in the Kuantan Singingi area, therefore this research is worthy of publication due to its originality. Nasdian (2014) stated that universities, NGOs, and other stakeholders could participate through institutional relations and social networking approaches in an endeavor to develop social groups. The network was formed for socio-economic group development by synergizing the functions of various stakeholders as a form of social capital development. Also, institutional development is crucial in the development of productive economic enterprises as it requires a high transaction cost. According to Nasdian, the development of small and medium-scale socio-economic groups needs to be the main target in community-based development activities. Through the development of such groups, it is hoped that it will be able to reduce unemployment, increase people's purchasing power, and in turn be able to have a dual impact, especially providing opportunities to develop local economic activities and productive efforts at the community level.

Literature review and conseptual framework
The development of social institutions is an alternative that could be developed with a collaborative institutional network approach from the community level to the locality level, exhibiting that the implementation of equality principles is more informal, participatory, has a strong commitment, and synergizes existing strengths. It is very helpful in solving problems and finding solutions in efforts to develop productive businesses at the community level.
The development of rubber plantation institutions reflects the extent to which social networks have been built into these communities. Should institutions and social networks aspects are synergized with the functions of various stakeholders of the rubber plantation business, will manifest into the development of social capital. Sujianto's (2008) study exhibited that social capital can provide community development and empowerment on the economic and social side through existing social networks.
The theory of social capital was originally developed by French sociologist Pierre Bourdieu, and by US sociologist James Coleman. Bourdieu states that there are three types of capital, namely money capital, social capital, and cultural capital. These three capitals would be more effectively used when there are social interactions or social relations. Social capital can be used for all purposes, but without physical resources and cultural knowledge, it will be difficult for individuals to build a social relationship. Social relations will be strong if the three elements above exist (Coleman, 1988). James Coleman defines social capital as a structure of relationships between individuals that enables them to create new values. According to Coleman, social capital is weak by processes that destroy kinships, such as divorce and separation, or migration. When families leave their existing kinship networks, friends, and other contacts, the value of their social capital will fall (Field, 2005, p. 140).
Social Capital is a resource seen as an investment to obtain new resources. The resources used for investment are called capital. Social capital is extensive and complex. Social capital is not interpreted based on material, but as a capital contained within someone. For example, in groups of family institutions, organizations, and all manners that lead to cooperation. Social capital places more emphasis on the potential of groups and relationship patterns between individuals in a group and between groups, with attention to beliefs, networks, norms, and values born from group members and become group norms.
Social affiliation thus refers to behaviors, such as the reduction in the physical distance from others, postural orientation towards others, looks at others, and any verbal or non-verbal behaviors (facial expressions, gestures) that initiate or maintain visual, physical, and/or verbal contact with another (see Delelis, 2002;Cottrell & Epley, 1977;Kulik & Mahler, 2000). The notion that social affiliation may be involved in emotion regulation appeared in Schachter's seminal work (Schachter, 1959) on the motives of social affiliation in anxiety-producing situations. What remains unclear is how, and under what conditions social affiliation contributes to emotion regulation.
Economic welfare is economic wellbeing expressed in terms of the sum of consumer and producer surplus-also known as community surplus. The concept of a welfare economy was coined by the famous economist Arthur Cecil Pigou in his work 'Economics of Welfare.' Welfare economics is a branch of economics that evaluates the distribution of goods and resources, the structure of markets and guides the making of public policies for the betterment of society. This theory is directly related to the distribution of income, economic efficiency, and the relation of these aspects with the well-being of society. The concept of welfare economics is a subjective study based on the assumption of how society and individuals in it can be benefitted or, in other words, how assumptions regarding policies can be constructed efficiently such that they can benefit the overall society. Famous economists who professed this theory are Vilfredo Pareto and Arthur Cecil Pigou.
Based on the definitions given by experts on social capital which broadly exhibited that social capital is an element crucial for the establishment of cooperation between individuals or groups or the establishment of collective cooperative behavior. Social capital is inseparable from the three main elements that exist in social capital which are used as elements of assessment on each type of social capital including (a) Trust (honesty, fairness, egalitarian attitude, tolerance, and generosity); (b) Social Networks (participation, reciprocity, solidarity, cooperation); (c) Norms (shared values, norms and sanctions, rules). The three elements of social capital above and its aspects are essential elements that exist or should exist in the life of a social group, whether the group is called community, ethnic group, and so on. In other words, social capital elements are necessary to run social structure engine well.

Material and method
Based on the research problems, this research utilized qualitative explanatory associative. It endeavored to explain causality relationships between exogenous and endogenous variables. The study was conducted on the rubber plantation community in Kuantan Singingi Regency, Riau Province. The sampling technique used in this study was non-probability sampling. Purposive sampling is a sampling method based on certain criteria. Researchers obtained a sample of 90 respondents. The research approach used was quantitative with an analysis tool Generalized Structure Component Analysis (GSCA) (Solimun, 2012).

Measurement model
The following

Assessing of linearity assumptions
Based on GSCA analysis, there is one assumption that must be fulfilled before an analysis is carried out, namely the linearity assumption, which requires a relationship between variables that are linear. The linearity assumption using the Curve Fit method is the relationship between variables declared linear when it meets one of the following two possibilities: (1) a linear model is significant (linear model sig <0.05), (2) a non-significant linear model and all possible non-significant models (sig linear model >0.05, and sig models other than linear >0.05). The assessment result exhibited the value of the linear model, which is <0.05, therefore, the model is linear and meet the assumptions set.

Goodness of fit
Based on the results of the structural feasibility assessment models measured using FIT and AFIT, FIT values of 0.871 and AFIT values of 0.795 were obtained. The FIT value indicated that the total diversity explained by the model is 87.1%. It indicates the formed model can explain all existing variables at 871.1%. The variables variety: Institutional Based Community, Social Network The overall feasibility of the model was measured using GFI and SRMR, GFI values of 0.941 and SRMR values of 0.031 were obtained. The GFI value is greater than 0.900 and the SRMS value is less than 0.08. It indicates that the model used is a good fit.

GSCA Analysis
Inner model (structural model) test assessed the research hypothesis. Hypothesis assessment was conducted using t-test (T-statistic) on each direct influence path in the partial method. GSCA and hypothesis assessment result is exhibited in the following Table 2 The table is presented in the following figure: Based on Table 2 and Figure 2, the inner model test results can be presented as follows: (1) Assessing the direct effect of Economic Institutional Based Community (X1) on Economic Welfare (Y), obtained a structural coefficient of 0.645, with a p-value of 0.000 < 0.05, it indicated a significant direct effect between Economic Institutional Based Community (X1) towards Economic Welfare (Y). The structural coefficient is positive, indicating that the relationship is positive. The higher the Institutional Based Community (X1), the higher the Economic Welfare will be (Y).
(2) Assessing the direct effect of Social Network Establishment (X2) on Economic Welfare (Y), obtained a structural coefficient of 0.558, with a p-value of 0.000. Because the p-value is <0.05, there is a significant direct effect between Social Network Establishment (X2) on Economic Welfare (Y). The structural coefficient is positive, indicating that the relationship is positive. the higher the value of the Social Network Establishment (X2), the higher the economic welfare will be (Y).
(3) Assessing the direct effect of Re-trust of Social Affiliation (X3) on Economic Welfare (Y), obtained a structural coefficient of 0.545, with a p-value of 0.000 < 0.05, then there is a significant direct effect between Re-trust of Social Affiliation (X3) against Economic Welfare (Y). The structural coefficient is positive, indicating that the relationship is positive. the higher the Re-trust of Social Affiliation (X3), the higher the Economic Welfare will be (Y).
(4) Assessing the direct effect of Regrouping Based on Capital Goods (X4) on Economic Welfare (Y), obtained a structural coefficient of 0.568, with a p-value of 0.000. Because the p-value <0.05, there is a significant direct effect between Regrouping Based on Capital Goods (X4) on Economic Welfare (Y). The structural coefficient is positive, indicating that the relationship is positive. the higher the value of the Regrouping Based on Capital Goods (X4), the higher the economic welfare will be (Y).
(5) Assessing the direct effect of Rubber Sector Coordination (M1) on Economic Welfare (Y), obtained a structural coefficient of 0.489, with a p-value of 0.000. Because the p-value is <0.05, there is a significant direct effect between Rubber Sector Coordination (M1) on Economic Welfare (Y). The structural coefficient is positive, indicating that the relationship is positive. the higher the Rubber Sector Coordination (M1), the higher the Economic Welfare will be (Y).
(6) Assessing the direct effect of the Rubber Sector Collaboration (M2) on Economic Welfare (Y), obtained a structural coefficient of 0.379, with a p-value of 0.000. Because the p-value is <0.05, there is a significant direct effect between the Rubber Sector Collaboration (M2) on Economic Welfare (Y). The structural coefficient is positive, indicating that the relationship is positive. the higher the Rubber Sector Collaboration (M2), the higher the Economic Welfare will be (Y).
GSCA analysis result obtained an interaction coefficient of 0.485 and P of 0.000. P-value <0.05 indicates that Rubber Sector Coordination is a moderator variable between the influence of Institutional Based Community (X1) on Economic Welfare (Y). The direct influence and the effect     GSCA analysis obtained an interaction coefficient of 0.473 and P of 0.000. P < 0.05 indicates that the Rubber Sector Collaboration is a moderator variable between the influence of Social Network Establishment (X2) on Economic Welfare (Y). Because the direct influence and the effect of interaction are equally significant for company performance (Y), the Rubber Sector Collaboration variable is a quasi-moderator. While the coefficient of interaction effect is positive, the Collaborative Rubber Sector (M) variable is reinforcing. Higher Rubber Sector Collaboration (M), increases the influence of the Social Network Establishment (X2) on Economic Welfare (Y). Figure 8.

Discussion
Analysis result indicated that Economic Institutional Based Community had an effect on Economic Welfare. A high Institutional Based Community Economic will have an impact on the high Economic Welfare. The direction expected from this economic informal institution is about how the newly formed informal institutions can play an active role in the rubber auction process which has been dominated by the Tauke and Joint Business Group (KUB). The output of the model offered for the short term is the realization of a community-based economic institution (based on capital owned) by making the rubber plantation community a pioneer (agent of change). At the next level, the new institution that has been formed is directed to be able to participate in the rubber auction system. The outcome to be achieved is the economic stability of the community as this is sought by the community (pragmatism mindset). Furthermore, at the mid-term level, this study concludes the importance of optimizing the function and role of formal institutions in the management locus of rubber commodities in Kuantan Singingi Regency as the rubber management in Kuantan Singingi is inadequate. There is a lack of coordination of authorized formal stakeholders and a grand design on a long-term basis. Therefore, when it comes to creating a rubber downstream business in Riau, it is absolutely necessary to have new rubber land that is more conducive and productive both in the environmental and socio-economic perspectives of the community.
Furthermore, this study also found that the Social Network Establishment had an effect on Economic Welfare. A high social network establishment will have an impact on the high economic welfare. The potential development of social networks is between rubber plantation communities and (especially) rubber processing companies and the Office of Industry and Trade of Kuantan Singingi Regency. In the context of social networking, it is also necessary to carry out strengthening of social networks, especially with the Office of Plantation in order to increase effectiveness and efficiency by providing rubber seeds. The development and strengthening of this social network will later be related to the immediate term stages in the form of strengthening the functions and roles of the formal institutions are also recommended in this study Re-trust of Social Affiliation has an effect on Economic Welfare. The high Re-trust of Social Affiliation will have an impact on the high economic welfare. The findings in the field indicate that in terms of a trust, the rubber plantation community in Kuantan Singingi was degraded after a rubber cooperative manager embezzled a certain amount of members' capital. This incident caused rubber farmers no longer believe in the presence of socio-economic institutions, which were based on cooperative systems. On the other hand, when it comes to developing social capital-based institutions, trust factors play a vital role. When the trust factor has been obtained, despite plummeting rubber price, it remained a mainstay of the community because the big name of Kuantan Singingi rubber has become popular. This is a form of positive capital not easy to build on national and international markets. Nevertheless, it remained that how the "human" or the rubber processing community itself maintain Kuantan Singingi rubber reputation. The Regrouping Based on Capital Goods has an effect on Economic Welfare. The high Regrouping Based on Capital Goods will have an impact on the high economic welfare. The urgency of regrouping rubber plantation communities is due to the diversity of capital goods ownership from the rubber plantation community itself. The results of the study exhibited that the general pattern found in the ownership structure of community capital exhibits that most rubber processing communities possess other capital. Capital goods ownership variations in the community structure of the Kuantan Singingi rubber plantation are as follows: (a) rubber and palm oil; (b) rubber and PETI (unauthorized mining); (c) rubber and cattle; (d) rubber, cattle, and palm oil. By conducting regrouping, it is expected that the people who initially have managed rubber will prevail using rubber commodity itself without having to hinder or prohibit them from managing/processing other commodities (livestock, PETI, and palm).

Conclusions and recommendations
Based on the results of the analysis, there is the influence of Society-Based Economic Institutions (X1), Social Networking (X2), Social Affiliation Trust (X3), and Capital Goods-Based Grouping (X4) and Rubber. Field Coordination (M1) and Rubber Sector Collaboration (M2) for Economic Welfare (Y). The higher the value of Community-Based Economic Institutions (X1), Social Networking (X2), Social Affiliation Trust (X3), and Capital Goods-Based Grouping (X4) and Rubber Sector Coordination (M1) and Collaboration in the Rubber Sector (M2) will have an impact on increasing economic welfare (Y). Also, it was found that the Rubber Sector Coordination and Collaboration moderated the influence of the variables of Society-Based Economic Institutions (X1), Social Networking (X2), Social Affiliation Trust (X3), and Quasi Regrouping Based on Moderators and Enrichment of Capital Goods (X4).
The results of this study provide support for the concepts offered by the Theory of Trust. If the application of this theory has been proven to apply to many public sector organizations in several countries, the results of this study provide additional understanding of its reliability to be applied to professional organizations in the mid-level rubber cultivation sector. Theoretically, the research results have implications for the development of social theory concepts, beliefs, and economic theories, particularly in public sector organizations. The results showed the domination of the influence of Community-Based Economic Institutions over other variables on the economic welfare of the community at the Rubber Plantation in Kuantan Singingi Regency which can be used as input and consideration for regulators (government) and the community in improving institutions. aspects of yield and efficiency of rubber production.
Based on the results of the study, it is suggested to increase the community's trust in rubber plantations to existing and/or established social institutions. If the trust aspect has been able to be returned to a better level, the next step is to build a joint commitment to the rubber sector in Kuantan Singingi. It must be maintained and enhanced both qualitatively and quantitatively. Furthermore, strengthening the awareness and ownership of the rubber plantation community on the existence of Kuantan Singingi rubber. This is considered important because the current trend for people to be indifferent and leave the rubber sector is simply due to a lack of innovation. On the other hand, it is necessary to take steps to integrate or collaborate rubber commodities with other commodities that have been "glimpsed" by the rubber plantation community in Kuantan Singingi. This study also recommends the need for a grand design related to land reclamation to increase the area of rubber land in Kuantan Singingi. This of course requires support from the central government level.
Likewise, to improve economic welfare, both the community and the government should increase the value aspects of the Community-Based Economic Institutions, Social Networking, Social Affiliation Trust, and Capital Goods-Based Grouping and Rubber Sector Coordination and Collaboration in the Rubber Sector with the hope of greater economic welfare good for the future.

Funding
The author received no direct funding for this research.

Citation information
Cite this article as: The effect of economic institutional based community, social network establishment, re-trust of social affiliation, and regrouping based on capital goods on economic welfare, Sujianto Sujianto, Cogent Business & Management (2021), 8: 1959008.


extfieldsofstudy:
Computer Science
Medicine

text:
A Comparative Study of BDS Triple-Frequency Ambiguity Fixing Approaches for RTK Positioning

Concerning the triple-frequency ambiguity resolution, in principle there are three different realizations. The first one is to fix all the ambiguities of the original frequencies together. However, it is also believed that fixing the combined integer ambiguities with longer wavelength, such as extra-wide-lane (EWL), wide-lane (WL), should be advantageous. Also, it is demonstrated that fixing sequentially EWL, WL and one type of original ambiguities provides better results, as the previously fixed ambiguities increase parameters’ precision for later fixings. In this paper, we undertake a comparative study of the three fixing approaches by means of experimental validation. In order to realize the three fixing approaches from the same information in terms of adjustment, we developed a processing strategy to provide fully consistent normal equations. We first generate the normal equation with the original undifferentiated carrier phase ambiguities, then map it into that with the combined and double-differenced ambiguities required by the individual approach for fixing. Four baselines of 258 m, 22 km, 47 km and 53 km are selected and processed in both static and kinematic mode using the three ambiguity-fixing approaches. Indicators including time of first fixed solution (TFFS), the correct fixing rate, positioning accuracy and RATIO are used to evaluate and investigate results. We also made a preliminary theoretical explanation of the results by looking into the decorrelation procedure of the ambiguity searching algorithm and the intermediate results. As conclusions, integrated searching of original ambiguities or combined ambiguities has almost the same fixing performance, whereas the sequential fixing of EWL, WL and B1 ambiguities overperforms the integrated searching. By the way, the third-frequency data can shorten the TFFS significantly but can hardly improve the positioning.

Introduction
Along with the development of real-time kinematic (RTK) positioning, the distance between users and the reference station is always a challenge in order to enlarge its service coverage. The crucial issue is the Integer Ambiguity Resolution (IAR) which is the key for fast and precise global navigation satellite system (GNSS) positioning. The phase observations with fixed ambiguity is equal to ultra-precise range observations and improve the other parameters remarkably, so that positioning accuracy of centimeter-or even millimeter-level can be achieved [1]. This challenge still remains a research hotspot in network RTK or PPP-RTK, since on the one hand, it could be applied to the ambiguity resolution of reference networks, and on the other hand, there are still a lot of regions with only very sparse reference stations.
In principle IAR performance depends heavily on the observation noise and modelling errors, wavelength of the ambiguity parameters and the observing geometry. The smaller the observation noise or the larger the wavelength, the better is the fixing. A stronger observing geometry will provide a better float solution, and consequently a more reliable fixing. However, several factors which cannot be efficiently removed by between-station differencing will remain in the observations used in the adjustment. The most significant ones are the tropospheric and ionospheric delays which are increasing with the inter-station distance [2]. They can be corrected with a priori correction models [3,4], respectively, however these corrections are not accurate enough for precise positioning, especially, the feasibility of ambiguity resolution in precise positioning is still related to the state of the ionosphere and can be affected by ionospheric disturbances of different scales [5]. Thus, additional parameters must be introduced and estimated simultaneously with the other parameters. This degrades the strength of the float solution and makes the instantaneous ambiguity resolution very difficult.
Thereafter, on the one hand, a number of studies worked on providing possible accurate ionospheric and tropospheric delay information and imposing proper constraints on them to retain the strength of the solution and reduced observation residuals and demonstrated their efficiency via the resulting ambiguity fixing performance [6][7][8]. On the other hand, from the processing of short baselines where dual-frequency data can be applied as independent observations, due to the very strong constraint among ambiguities of different frequencies, ambiguities can be fixed very quickly and reliably [9]. Fortunately, the third-frequency observations available in the modernized GPS and the newly emerging system BDS and Galileo bring obviously additional information on ionospheric delays and will certainly enhance the constraint among ambiguities at different frequencies [10][11][12]. The ambiguity resolution of multi-frequency are investigated thoroughly together with the consideration of multi-constellations, as both sides rely tightly on each other through more information in ionospheric delays.
Beside strengthening the solution and reducing the errors to increase the fixing opportunity as discussed above, to quickly and precisely find out the correct integer solution is of the same importance, i.e., the fixing approaches. In the past, several approaches were developed, most of them are based on searching the optimal integer ambiguities statistically closest to the corresponding float solution [9][10][11][12][13][14][15][16][17]. The mostly used approach is the Least-Squares AMBiguity Decorrelation Adjustment (LAMBDA) method [1], for example, which is also applied in the recent contributions [18][19][20].
A large number of scholars have studied the long-range RTK positioning using triplefrequency observations with the concentration of triple-frequency ambiguity resolution and observation modelling as well, One group of scientists tried to resolve the ambiguities directly based on a set of optimized combinations of the phase and range observations, which is referred as triple-frequency carrier-phase ambiguity resolution (TCAR). It is proved that EWL and WL can be fixed reliably, however the last one with a rather shorter wavelength can only be fixed using geometrical observations and ambiguity searching approaches. It is also confirmed TCAR or the later improved method does not have obvious advantage to some of the searching schemes [21]. For example, [12,22] has studied the contribution of the third-frequency of GPS and Galileo on long-range RTK ambiguity resolution, the results of a cascading fixing scheme in the sequence of EWL, WL and L1/E5 ambiguities are very promising although only simulated observations were available and employed.
Overall, there are also several different schemes for multi-frequency ambiguity resolution: (1). to searching all ambiguities of the original frequencies, and (2) to map the ambiguities to extra-wide lane (EWL), wide-lane (WL) and one original frequency and to search the integer solution altogether, (3) similar to (2) but search the ambiguities of the same wavelength in batch with that of the longest wavelength first. Usually the second one is believed better than the first one as the EWL and WL ambiguities should be easily fixed due to their longer wavelength. However, there is also the opinion that LAMBDA should give the same result as the first two schemes are based on the same float solution information. A comparative study of the three schemes is the major goal of this contribution.
The Chinese BeiDou Navigation Satellite System (BDS) should be the first one of the GNSS which provides triple-frequency data across the constellation. In fact, it can provide signals at more than triple frequencies since 27 December 2018 [23]. Nowadays, Galileo and all GPS Block IIF and III satellites are also transmitting signals at third or even more frequencies. In this contribution, the comparative study on the above-mentioned three fixing schemes is carried out based on the BDS triple-frequency data. In order to do a fair comparison, all the fixing will start from exactly the same normal equation information. We use double-differenced but uncombined observations with undifferenced ambiguities at the original frequencies, which should be equivalent to that using undifferenced or single-differenced observations if their covariance matrix is rigorously considered [24]. Before the ambiguity resolutions these ambiguities are mapped to the combination among frequencies, i.e., EWL, WL etc. and between stations and satellites, i.e., double-differenced ambiguities for fixing. The time of first fixed solution (TFFS) and the correct fixing rate, and of course also the derived positioning deviations are employed as fixing performance indicators for evaluation. The result of the comparison will be further investigated by examining the theory and algorithm of the ambiguity resolution schemes. The paper is arranged as follows: the related functional and stochastically models will be introduced with the handling of atmosphere delays after this introduction. Then, ambiguity mapping method is presented for both among frequencies and between stations and satellites. Three ambiguity fixing schemes are defined for comparison based on the ambiguity searching approach. Afterwards, data processing strategy is presented, and results are compared and analyzed before conclusions are drawn.

Estimating Model
For a better understanding of the ambiguity resolution schemes, we briefly introduce the observation equations, and the handling of tropospheric and ionospheric delays, as well as ambiguity mapping for triple-frequency data processing.

Observation Equation
The observation equation of GNSS multi-frequency undifferenced and uncombined carrier phase can be expressed as: where p, r and i are the indicates of the satellite, receiver and carrier phase frequency, respectively; X r is the coordinate parameters vector after first order Taylor expression with design matrix H p r ; φ p i,r is the carrier phase observation, ρ p r is the geometric distance from the satellite to the receiver r; dt r and dt p are the clock error of the receiver and satellite, respectively; c represents the speed of light; T r is the zenith tropospheric delay; M p r is tropospheric mapping function; µ i is the ionospheric coefficient for frequency i, µ i = f 2 1 / f 2 i ; I p 1,r is the ionospheric delay of the first frequency, and λ i is the wavelength of carrier phase; B p i,r is the carrier phase ambiguity, ε p i,r is the measurement noise, P p i,r is the weight matrix of the observations and which value is calculated according to the altitude angle of the satellite.
In the baseline or RTK processing, the single-differenced (SD) observation between stations can be applied, as it can eliminate the satellite clock biases and reduce the effect of satellite orbit biases and atmospheric delays. In order to also eliminate the receiver clock bias, further difference between satellites is usually employed in practice, i.e., doubledifferenced (DD) observations are used. Suppose that two stations simultaneously observe the reference satellite q, and the station coordinate and zenith tropospheric delay at the reference stations are not estimated, only those at the moving station are estimated, the DD observation equations read as: where r, s represent rover station and reference station, respectively, v pq i,rs represents observation residual. As is well known, using DD or SD or undifferenced observations are equivalent if their covariance matrices are considered rigorously [24].
Usually, DD ambiguities should be involved in the DD observation equations, as they have the nature integer characteristics and can be fixed directly. However, the DD ambiguities may change along with satellite rises and sets, so the normal equation must be reformed accordingly. This makes the realization of the algorithm much more inconvenient. Therefore, SD ambiguity parameters are used in this study and then are mapped into DD ambiguities before integer ambiguity resolution.

Atmospheric Delay Estimation
The atmospheric delay is one of the major error sources that limits the accuracy and reliability of RTK positioning and the impact becomes more serious as the inter-station distance increases, since both troposphere and ionosphere delays can be neither cancelled by between stations differencing nor precisely corrected. For tropospheric delays large height differences between stations can also cause similar inaccurate modelling due to the different path of the transmitting signals [25].
In practice, the dry component of the tropospheric delay can be corrected by a model, for example the Saastamoinen model using pressure and temperature from GPT2 [26,27], whereas the wet component is parameterized as either piece-wise liner function or randomwalk process to be estimated. Regarding to the mapping functions, the mostly used are GMF and VMF1 and the later one is used in this study [28,29].
Ionosphere delays can be eliminated by forming the ionosphere-free combination from observations of at least two frequencies. This is very often applied in long-range baseline data processing and the ambiguity of ionosphere-free observations must be separated into wide-lane and narrow-lane and fixed sequentially [30,31]. A more advantageous strategy is to use uncombined observations with slant ionospheric delays as unknown where the ambiguity of original frequencies can be estimated by imposing constraints on spatial and/or temporal ionospheric delay variations [32]. The later one is extremely suitable for triple-frequency data processing to avoid the selection of ionosphere-free combinations to involved in the data processing. Therefore, the later one is employed in this study with temporal constraint, i.e., inter-epoch ionospheric delay parameter constraint.
The random walk process for both ionospheric and tropospheric delays can be mathematically expressed as: where E(·), Cov(·) are the functions for expectation and, variance, respectively, and t is epoch time. The symbols d, w, σ 2 , and q denote residual value, change of residual value, variance and power density of the corresponding atmospheric parameter, respectively. For RTK positioning, the remaining delay in the SD observations between stations depends on the inter-station distance, so the power density should be fine-tuned according to the baseline length. For example, both ionospheric and zenith tropospheric delays are almost fully eliminated for short baselines, the constraint of the initial state and the power density should be very close to zero, so that all estimates of state parameter are zero. In practice, the power density of the random walk processes can be sophistically adjusted in order to optimize the positioning performance, for example, fine-tuning the power density during the initial period to shorten the ambiguity convergence time and to improve positioning accuracy and reliability [33]. It can also be adapted according to the ionospheric variation calculated from ionosphere combinations of phase observations or the estimated delay parameters in the previous epochs [34].

Parameter Estimation
The DD observation equation system of Equation (3) is abbreviated for the presentation of parameter estimation and ambiguity resolution as: where x represent position parameters; y includes the residual zenith tropospheric delay after applying corrections from VMF1 model and the slant ionospheric delay parameters each for one satellite, i.e., (T r , I 1 1 , I 2 1 , · · · I nsat 1 ) T with nsat as the number of satellites observed at the epoch, and will be estimated as random walk process as described in Section 2.2; and b represents SD ambiguity parameters; and A, D, C are the corresponding coefficient matrices for parameters x, y and b, respectively; l is pre-fit observation residuals. It should be mentioned that for kinematic positioning the position parameters are time-varying instead of constant parameters all over the time for the static positioning. According to the least square's criterion, the normal equation of Equation (5) is which is abbreviated as: This normal equation system is estimable even if a loose constraint on the reference ambiguity parameters is applied.

Ambiguity Mapping
As is mentioned, the ambiguities in Equations (5)-(7) are SD-ambiguities between the two stations at the original signal frequencies in order to avoid the complicated reformation due to satellites rises or sets. Since only DD-ambiguities have integer characteristics because of the existence of uncalibrated phase delays, the SD ambiguities must be mapped into DD ambiguities for integer ambiguity resolution. This re-parameterization can be realized by the normal equation mapping based on a defined parameter transformation proposed by Blewitt and Dong and Bock [30,31].
From a given set of SD ambiguities, a number of sets of independent DD-ambiguities can be defined. There are various criteria to select the ones to be fixed the most easily, for example according to their elevations and length of continuous observations etc. Obviously from SD to DD the mapping is rank-defect, independent SD ambiguities must be further added to complete the transformation to full-rank. Assume that a full rank transformation then the corresponding normal equation system can be derived from Equation (6) as: Besides the mapping from SD to DD ambiguities, ambiguities combination of various frequencies can also be implemented with this mapping approach. For the triple-frequency data processing, the mapping can be generalized defined as: where (i 1 , j 1 , k 1 ), (i 2 , j 2 , k 2 ) and (i 3 , j 3 , k 3 ) are the combination coefficients and they must be integer values, so that the mapped ambiguities also have integer characteristics for fixing. For example, with (0, 1, −1), (1, −1, 0), and (1, 0, 0). we obtain the EWL, WL and B1 as mapped ambiguities. The normal equation system with the mapped ambiguities can be derived in the same way as for the mapping from SD to DD ambiguities. Usually, it is believed that ambiguities with longer wavelength can be fixed more easily. It might be true for most of the searching approaches but not for the integrated searching with LAMBDA method. This will be demonstrated and investigated in this study.
In principle, the former mapping, from SD to DD ambiguities, is definitely needed for integer ambiguity resolution, whereas the later one is carried out on demand.

Ambiguity Resolution
There are two widely used ambiguity fixing methods: the rounding method and the searching method. The rounding method, as is named, directly rounds the estimated float ambiguity to its nearest integer. This is the simplest method to obtain the integer ambiguity, however it requires the float estimate with a certain accuracy so that the rounding is reliable. Therefore, it is mostly used in static post processing with long-term observations. The searching method is to find the integer ambiguities mostly compatible to the float solution according to the float estimates and their covariance. The most famous one is the LAMBDA method proposed by Teunissen [1].

Ambiguity Searching Approaches
Statistically speaking, the float estimates and their covariance matrix of the ambiguities define the confidence ellipsoid to which the integer solution belongs. However, how to find out the integer candidate vectors from this confidence ellipsoid has always been a key issue, especially in the case that the quantity of the observations is rather small, and the observation geometry is very poor.
Before LAMBDA the searching methods define the search areas according to the STD of each individual ambiguity [1]. Due to high correlation of ambiguities, the searching area defined this way can be extremely large than the confidence ellipsoid, so that thousands of candidates are included and must be checked and results into very long computation time. Although it was improved by introducing further constraint on ambiguity differences [9], but it is still not good enough for ambiguity fixing using observations of a single epoch or few epochs [35]. The key aspect of LAMBDA method is to find out the best integer linear combination of original ambiguities, so that the correlation among them is at the lowest. This means the searching region can be rather precisely defined by the STD of each individual ambiguity, and only few candidates are identified and checked for fixing.
The LAMBDA method is widely used for integer ambiguity searching, especially for fast ambiguity fixing with short-term observations. It can make full use of the original observation information, as is proved, it can maximize the success rate of ambiguity fixing [14]. We use LAMBDA to search the best ambiguity in the study.
Instead of directly searching integer ambiguity, Gaussian integer transformation is carried out on the ambiguity matrix to reduce the correlation of the ambiguities in the LAMBDA method. Afterwards, the integer ambiguities are converted to the original to obtain the parameters of the fixed solution. Since the task of the LAMBDA approach is to find out the best integer combination among the original ambiguities, therefore it should reach the same integer combination in principle no matter the searching starts with ambiguities of the original frequencies or that containing EWL and WL in spite of that the later ones are usually believed to be fixed more easily. More specifically, transforming original ambiguities into them before applying LAMBDA is not necessary which will be investigated and validated later on.

Ambiguity Resolution Strategy
For the ambiguity resolution of a special positioning solution, the searching approach can be applied in two different ways: the first one is to search all ambiguities simultaneously, whereas the second is to do the searching in sequential batch, also referred as to stepwise approach or cascading approach. In the second approach, the set of ambiguities which can be fixed most likely is selected and fixed, and then the solution is updated with the fixed ambiguities. Obviously, the remaining ambiguities in the updated solution will have a better quality, so that they can be more easily and reliably fixed now. This procedure can be repeated until no more ambiguities can be fixed. Usually, taking the advantage of the long wavelength of EWL and WL ambiguities, they should be fixed firstly to improve the solution and consequently to assist in resolving the short-wavelength ambiguities. In fact, the partial fixing can be considered as a special case of this approach, in which difficultly to be fixed ambiguities, for example, ambiguities of a rising satellites with very few observations or at a very low elevation, or those with the largest variance, will be excluded in the first round of searching, and will be fixed after their estimates are improved by the already fixed ambiguities [36]. From this point of view, this approach is superior to the first one, especially in the case that the first one does not work properly, i.e., cannot fix all ambiguities at once.
We adapt the stepwise fixing approach to the ambiguity resolution of triple-frequency GNSS observations to resolve the integer ambiguities step-by-step according to their wavelength, and to investigate its fixing efficiency.
After applying the ambiguity mapping from the original ambiguity (b 1 , b 2 , b 3 ) to the EWL, WL and L1 ambiguities (b EW , b W , b 1 ), we obtain the normal equation for fixing. In order to facilitate the distinction between the DD-ambiguity parameters and other non-DDambiguity parameters, x, y and other non-ambiguity parameters are collectively referred to as X. The normal equation system for ambiguity fixing reads as The ambiguity searching should be carried out at least three times, i.e., for EWL, WL and B1, respectively, if no partial fixing is required. After each successful fixing, the solution is updated by inserting the fixed integer values into the normal equations and the updated ambiguities and their Q matrix are employed in the next searching.
Taking the first step as example, applying LAMBDA searching approach to Q b EW b EW and b EW of the float solution, we obtain the integer EWL ambiguity vector b EW . Inserting the integer into the normal equations, the normal equations with WL and B1 ambiguities can be expressed as: Solve the normal equation system, the updated solution with fixed EWL ambiguities is derived. Then we can repeat the above process for WL ambiguities by applying LAMBDA searching to the estimates and their co-variance matrix.
There is also algorithm for updating the solution with fixed ambiguities without resolving the normal equations again, which can be find in a number of references, i.e., in Teunissen [21]. Considering the computer capacity nowadays and for better understanding, we resolve the normal equations in this study.

Ambiguity Resolution Scheme
Based on the above discussion, there could be three ambiguity fixing schemes for triple-frequency data processing: integrated fixing of DD ambiguities of original signal frequencies, integrated fixing of DD ambiguities of combined frequencies with longer wavelength, and stepwise/cascading fixing of DD ambiguities of combined frequencies, referred as to Scheme 1 (S1), Scheme 2 (S2) and Scheme 3 (S3), respectively. Here the combined frequencies are EWL, WL and one original frequency B1.
As uncombined observations are mostly suggested for multi-frequency data processing for making use of ionosphere constraints and to avoid the selection of ionosphere-free combinations for triple-frequency data, S1 is the natural way for ambiguity fixing. However, mapping to EWL, WL and B1 for fixing is believed more efficient because of the longer wavelength of EWL and WL ambiguities. As already mentioned, the LAMBDA method tries to reach the best linear integer combination of ambiguities to reduce the correlation to obtain exactly the confidence ellipsoid containing the best integer candidate. The combination includes also those among the ambiguities of different frequencies. This means that LAMBDA method may automatically achieve the same or similar combination no matter whether original or mapped ambiguities are employed. Therefore, S2 is designed for a numerical clarification of such a misunderstanding. S3 is actually to improve the fixability of ambiguities sequentially by fixing the most easily to be fixed ones in advance. It should be the best strategy if enough satellites are observed.
For all the three fixing schemes, DD observations of the signal frequencies are used and normal equation with SD ambiguities are formed as in Equation (6). Then the SD ambiguities are mapped into DD ambiguities using Equation (8). For Schemes 2 and 3 the DD ambiguities of the signal frequencies are further mapped into EWL, WL and B1 with Equation (10). In this way, the float solution of the three schemes are exactly the same, so that we can evaluate the efficiency of the ambiguity fixing.

Experimental Validation
In order to validate the implemented algorithm for medium-range RTK positioning, especially the comparison of the three ambiguity resolution schemes, data of baselines with various inter-station distances are collected and processed using the three fixing scheme with dual-and triple-frequency data of GPS and BDS and in static and kinematic mode, respectively. the atmosphere after double difference can effectively eliminate for short baseline, so the short baseline cases get conclusion is almost not affected by the atmospheric. However, atmosphere is the main influence factors for the long baseline. Whether conclusions change as baseline length raising, this is also one of the research objects in this paper. In this paper, the conclusions of the above three schemes were first obtained by statistical analysis of the short baseline data, and then further verified in the long baseline, and the results are analyzed and investigated.

Experimental Data
In total four baselines are selected with from different networks with baseline length of 258 m, 22 km, 47 km and 53 km representing short and middle-range baselines, respectively. Baseline A is the data collected using the receiver, both stations are equipped with SinoGNSS K708. Baseline B between station CASM and LFSL stations are equipped with UNICORE UB4B0 receiver. Baseline A and B are measured data. Baseline C1 between LNSB and LNTL are from the Liaoning CORS network, all with Trimble Net R9. Baseline C2 of Guangdong CORS network data, the receivers are SinoGNSS K708.
All the receivers are capable of tracking GPS and BDS triple-frequency data and the sampling rate is 1 Hz. However, as not all GPS satellites are transmitting L5 signals only dual-frequency data is utilized for GPS in this study. The information of the baselines is listed in Table 1. Although baselines C1 and C2 are of similar length, they are in the north and south of China, respectively, and therefore may have quite different ionosphere impacts. Among the four baselines, the average number of GPS observables is about 7. Baseline B and C1 have almost the same number of BDS and GPS satellites, while Baseline A and C2 have at least two more BDS satellites on average than GPS.

Data Processing
The data of each baseline is divided into 10 min sessions and processed independently in static and kinematic mode, respectively. For each 10-min session, data of dual-frequency of GPS and both dual-and triple-frequency of BDS are processed with the three ambiguity fixing schemes. The solutions are named according to the systems (BDS, GPS), frequencies (TF, DF) and fixing schemes (S1, S2, S3) accordingly. In total, for each session we have nine solutions for either static or kinematic processing, namely BDS-TF-S1, BDS-TF-S2, BDS-TF-S2; BDS-DF-S1, BDS-DF-S2, BDS-DF-S3; and GPS-DF-S1, GPS-DF-S2 and GPS-DF-S3. For all the solutions, the epoch interval is 1 s and cutoff elevation is 15 degree. In the data preprocessing, outliers are detected according to the posterior phase residuals and residuals larger than 0.5 cycles is considered as cycle slips which will not be repaired but modelled with a new ambiguity parameter. For making the fixing decision, RATIO-test is employed. For the float ambiguity vector b, and ambiguity variance matrix Q bb , two integer ambiguity candidates with the smallest quadratic form (n − b) T Q −1 bb (n − b) are selected. The proportion of the value of the second smallest quadratic form with respect to the smallest one is defined as RATIO. If RATIO is larger than a threshold, the integer candidate with the smallest quadric form is the fixed integer vector. This threshold is 2.0 in this study. In order to verify the ambiguity resolution performance, RATIO, TFFS and the rate of correct fixings are used as indices for its convergence and fixing reliability, respectively. The correct fixing rate is defined as the percentage of the correctly fixed epochs as: P = N suc N all (13) where N suc represents the number of epochs where ambiguities are correctly fixed, and N all the total number of epochs. The correctness of a fixing is confirmed by full consistence of the integer ambiguities with that derived from all data of the baseline in static mode using the high-precision data processing software package PANDA [37]. Atmospheric errors have a significant impact on the success of ambiguity resolution, especially in the case of long baselines. The influence of atmospheric errors on the short baseline, for example shorter than 10 km, is negligible due to the strong spatial correlation between the two stations. For middle and long baselines, the troposphere and ionospheric delays are estimated with constraints, the degree of improvement depends on the accuracy of the imposed tropospheric and ionospheric constraints [8]. The tropospheric zenith delay is constrained with an initial STD of 2 mm, and a power density of 10 mm/sqrt(hour) for Baseline B, C1 and C2. For the ionospheric delay parameters, the a priori value for the zenith direction is selected according to the linear function: σ I [cm] = 0.04 · length[km] [34,38], then the power density of the temporal ionospheric variation is calculated from ionosphere combinations of phase observations of the previous epochs and a simple mapping function 1/sin(E) is applied to obtain the values for the slant delay parameters of each satellite.

Statistics of Ambiguity Resolution
We first have a close look at the statistics of ambiguity fixing, namely the TFFS and the correct fixing rate for both static and kinematic processing. As the statistics for both processing modes have very similar behavior, we will present them together and more detailed analysis will be carried out for the kinematic mode since we are more interested in medium-range RTK positioning. Table 2 shows the statistics of TFFS and the correct fixing rate of the nine solutions in both the static and kinematic mode for Baseline A, in which the TFFS are in unit of epochs and it is consistent with the length of time in seconds as the sampling rate is 1 s.

BDS-TF
BDS-DF GPS-DF S1 S2 S3 S1 S2 S3 S1 S2 S3 From Table 2, for the short baseline, the ambiguity convergence times of the fixing scheme S1 and S2 are the same for both GPS and BDS, and DF and TF observations. The stepwise fixing scheme S3 can further improve the fixing by reducing the ambiguity convergence time for BDS-DF and GPS-DF data, from 4.0 s to 3.0 s and 5.0 s to 4.0 s, respectively and increasing the correct fixing rates, although the magnitude is rather small. From the statistics of the kinematic positioning shown on the second half of the table, only few more epochs are needed for the first fixing and the correct fixing rate is slight decreased, and the overall behavior of the three processing scheme is almost the same compared with that of the static mode.

Static
The same statistics for the middle range baseline B are listed in Table 3. From Table 3 we also observed the similar phenomena as for Baseline A that the ambiguity convergence of S1, S2 are the same and S3 brings additional improvement compared to S1 and S2 BDS. The triple-frequency data shorten the TFFS very significantly for this middle range baseline which is hardly noticeable for the short baseline. It should be pointed out that the contribution of S3 is obvious, which is rather small in the short baseline. Compared with the static mode, the convergence gets clearly slower and the correct fixing rate decreases of several points in percentage as well for the kinematic mode, except the triple-frequency BDS data, especially the processing scheme S3 where fixing can be achieved again within a few epochs on average as for the short baseline.

BDS-TF
BDS-DF GPS-DF S1 S2 S3 S1 S2 S3 S1 S2 S3  Tables 4 and 5 show the same statistics for the two medium-range baselines: Baseline C1 and C2 where remaining ionospheric and tropospheric delays in differenced observations get larger and must be handled properly as above described. From Tables 4 and 5 we can again obtain the similar conclusions as that for the medium baseline B that S1 and S2 have the same fixing performance for both GPS and BDS and for all data types DF and TF, and the impact of the third-frequency data is even stronger, whereas S3 can reduce the TFFS for DF data markable and keep the TFFS for TF of few seconds as for the short baseline. Table 4. Statistics of ambiguity fixing performance of Baseline C1, Sampling rate 1 s.

TFFS(s)
BDS-TF BDS-DF GPS-DF S1 S2 S3 S1 S2 S3 S1 S2 S3 From the overall statistics of ambiguity resolution, S1 and S2 have the same performance in terms of TFFS and the correct fixing rate for all the solutions: DF and TF of BDS data and DF of GPS, while S3 overperforms S1 and S2 obviously except the short baseline. When we have a close look at the statistics of the S3 solutions, the performance seems correlated with the baseline length, the advantage of S3 is better over longer baselines. One possible reason could be that over short baselines the float solution is already very strong and precise, so that sequential fixing could not further improve the solution.

Static
To confirm this possible reason, we furthermore calculated the RMS of posterior phase residuals of the float solutions of the baseline A and C1 and depicted on Figure 1. The RMS of the residuals of the short Baseline A shown in the top panel is only about 0.003 cycles which is consistent with the observation noise as most of the biases are well cancelled by the between-station differencing. The RMS of the two systems are comparable and BDS has a slight larger RMS. For the medium-range Baseline C1 the RMS increases to 0.02 cycles because of inaccurately modeled biases, such as atmospheric delays, multi-path and even satellite orbits. As is well known, larger phase residuals will make the B1 ambiguity fixing difficult, therefore fixing EWL and WL in advance will improve the solutions for B1 fixing. That is why the advantage of the stepwise fixing S3 over S1 and S2 can be better reflected over medium-range baseline. The atmospheric delay errors could be one of the major factors affecting the ambiguity resolution, especially in the case of long baselines. Taking Baseline C1 processed in kinematic mode and in 10-min sessions as an example, the estimated tropospheric zenith delays (top-left panel) and the slant ionospheric slant delays of all satellites are illustrated in Figure 2. The Ionospheric delays derived from the B1 and B2 phase observations, referred to as "observed" are also shown in the corresponding panels. The estimates before the ambiguity-fixing are not plotted, as there is usually a big jump in the estimated time series due to the big change of ambiguities from their float values to the integer ones. Therefore, there is usually a short gap between two adjacent sessions which is hardly visible due the high sampling rate.
As can be seen from the Figure 2, the remaining tropospheric zenith delay is rather stable and fluctuates within (−0.01 m, 0.01 m) with a slight increasing trend. The estimated ionospheric slant delays are of the similar amplitude, but the temporal changing is much fast. More important is that both the estimated and observed agree with each other very well after aligned by a constant. From this example, the ionospheric and tropospheric delays should be estimated in order to reduce their impact on the other parameters, especially the ambiguity parameters and their fixing.

Comparison of Fixing Index RATIO
To further compare the three fixing schemes, it is worth to check the fixing index RATIO at each epoch. Taking the BDS-TF solutions of Baseline A and Baseline C1 as example, the RATIO time series of S1 and S2 and S3 of Baseline A are depicted in Figure 3, where S3 has three RATIOs with respect to the fixing of EWL, WL and L1 ambiguities, in order to express the state of RATIO time series more clearly, we set an upper limit of 1500 for RATIO value. From the top panel of Figure 3, the RATIO time series of S1 and S2 are almost the same, their differences shown in the middle panel are within ±0.01. This clearly proves that the preset integer transformation, i.e., from raw to EWL and/or WL ambiguities, has almost no influence on the ambiguity fixing using LAMBD method. Theoretically, LAMBDA algorithm can find out the best Z transformation to narrow the searching space. However, when we compare the transformation matrix, the mapped ambiguities for S1 and S2 are not exactly the same. The reason is most likely the numerical calculation and round from real-value to integer value of the z-transformation. However, it should be further investigated in future. In the processing with S3, EWL, WL and B1 are fixed sequentially and the solution is updated after each fixing. The RATIO of EWL and WL fixing of S3 shown in the bottom panel is on average about 2-3 times larger than that of S1 and S2. This reflects that EWL as well as WL can be fixed much more easily than fixing all ambiguities together. After both EWL and WL are fixed, the updated solution with B1 ambiguities is improved, the related ambiguities can be fixed with a larger RATIO overall than that of S1 and S2 which can be seen from the top panel of Figure 3.
Correspondingly, the RATIO time series of S1 and S2 and S3 are depicted in Figure 4 for the BDS-TF solutions of Baseline C1 in static mode. From Figure 4, the RATIO time series show very similar characteristics except the RATIO values get a little bit smaller due to estimating more parameters in static mode. The same studying is carried out to the other three baselines and we reveal the same fact. For the above four baselines, S1 and S2 are equivalent in the ambiguity convergence measured by the fixing index RATIO in this experiment, and S3 performs usually better than S1 and S2, regardless of GNSS systems, data types and the positioning modes, especially for medium-range baselines.
In summary, regardless of the length of baselines and the use of DF or TF observations, the ambiguity convergence time indicators of S1 and S2 are the same, and the preset Z transform has almost no effect on the LAMBDA method, while the scheme S3 is superior to S1 and S2 for multi-frequency data.

Positioning Accuracy
The position differences of the fixed solutions with respect to the static solution using all the observations can be used to qualify the positioning accuracy. S3 scheme is superior to S1 and S2 scheme in terms of ambiguity resolution, in order to more clearly and intuitively reflect the positioning result accuracy of multi-frequency data, and in view of the equivalence between S1 and S2 scheme, we use the positioning result of S1 scheme which is superior to the suboptimal performance to calculate the RMS value. The RMS of the position differences of the S1 is calculated for the four baselines using BDS-TF, BDS-DF and GPS-DF observations and listed in Table 6 for static and kinematic processing mode in unit of m.
It can be seen from the table that in static mode the four baselines are better than 2 cm in the north and east directions, and 3 cm in the up direction; and in the kinematic mode similar accuracy can be accomplished in N, E direction and better than 5 cm in up directions.
The positioning results using BDS triple-frequency and dual-frequency data are also compared on Baseline A and B. They are very similar and comparable with GPS dualfrequency. That means the advantage of triple-frequency data in positioning accuracy is hardly visible.

Further Comparison of Ambiguity Resolution Strategies
The core of the LAMBDA method is the reduction process of the correlation of the ambiguity variance covariance matrix through the Z transformation. The final result of Z transformation is the diagonalization of the ambiguity factor matrix. Taking the static processing of BDS triple-frequency positioning of Baseline A as example, the influence of S1, S2 and S3 schemes on Z transformation is discussed for a better understanding of their performance. Figure 5 is the eigenvalue obtained after reducing the correlation. It can be inferred that S1 and S2 are similar in ellipse flatness from the definition of the condition number after Z transformation in Figure 5. This indicates that the external preset integer transformation applied in S2 has little effect on the correlation reduction of the entire ambiguity search space. However, they are not exactly the same, that means the realization of the best Z-transformation is not unique although theoretically it should be. One possibility could be that the numerical computation from a real-valued to an integer-valued transformation may depends on the starting ambiguity types. This should be further studied in the future. Of course, the final estimated integer ambiguities of original frequencies must be the same for all the three strategies, so that the positioning accuracy of the fixing solutions can be guaranteed. Figure 5. Example of the Eigenvalues of the covariance matrix after Z transformation of S1 and S2 for Baseline A.
In Figure 6, from left to right are the two-dimensional diagram of the ambiguity cofactor matrix after Z-transformation of the 33rd epoch obtained from the three schemes: S1, S2 and S3. For S3, there are three blocks with respect to the EWL, WL and B1 ambiguities, as they are fixed in three sequential steps, and the cofactor matrices are based on the updated solution by introducing the integer ambiguities fixed in the previous blocks. The colors represent the values of the corresponding matrix elements and non-diagonal elements illustrate the correlations of the ambiguities. Therefore, after Z-transformation, the smaller the non-diagonal elements, the higher the degree of diagonalization of the matrix, and the better the effect of correlation reduction. Figure 6. Element values of the covariance matrix after Z-transformation for S1 and S2 and S3 from left to right, respectively. S3 has three blocks for EWL, WL and B1 sequentially.
From Figure 6, although the cofactor matrix elements of the EWL ambiguities are larger than that of S1 and S2, the element values of the WL ambiguities and the B1 ambiguities are obviously smaller than those of S1 and S2. This is clearly the contribution of previously fixed integer ambiguities which enhances the solutions and reduces the correlations among the parameters. The theoretical principle can be found in [21] in which it is pointed out that the conditional standard deviations of the optimal combination of ambiguities are very small, and it is feasible to use this advantage to calculate the correct ambiguity within a few epochs. That is why in this study we have designed the S3 scheme to verify this conclusion from the experimental point of view.
The condition number is the ratio of the maximum eigenvalue to the minimum eigenvalue of the normal equation (≥1). The closer it is to 1, the closer the major and minor axes of the ellipsoid are, and the more they tend to be a sphere. Figure 7 is the number of ambiguity resolution conditions for each of S1, S2 and S3. The condition numbers of the S1 and S2 schemes have a relatively small difference, which means that there is little difference between the two in the searching ellipsoid. The condition number obtained by the S3 scheme to solve the ambiguity of the extra-wide lane is the smallest overall. The search ellipse of the wide lane ambiguity search is still smaller than S1 and S2. For the search interval of the original ambiguity, after the first two steps of correction, the accuracy of floating solution of the original ambiguity has been well corrected. Compared with S1 and S2, although the original ambiguity condition number in S3 is larger, it is smaller than S1 and S2 in the initial stage. It can be seen from Figure 3 that in terms of the Ratio value, the ambiguity obtained by the S3 scheme is more reliable, and the experimental results show that this conclusion is valid. Since ambiguities of long wavelength have the advantage to be fixed and integrated searching of very high dimension of ambiguities could result into computational difficulty, so S3 should be a better choice to effectively shorten the fixed time of multi-frequency ambiguity resolution.

Conclusions
In the comparative study of the three triple-frequency ambiguity resolution schemes, medium-range RTK is realized by applying adapted temporal constraint to ionospheric slant delay parameters and tropospheric zenith delay parameters. Starting from the same normal equation as the uncombined and undifferenced carrier phase ambiguities, the normal equation with the desired double-differenced ambiguity parameters are generated through ambiguity mapping in order to have a fair comparison. Three fixing schemes, namely S1 fixing all ambiguities of original frequencies simultaneously, S2 fixing all ambiguities of the EWL, WL and one original frequency together, and S3 fixing the same ambiguities as S2 but in a cascading way, are applied to four baselines of 258 m, 22 km, 47 km and 53 km for using BDS triple-frequency and dual-frequency data and GPS dualfrequency data and in both static and kinematic mode. The data is divided into 10 min pieces or sessions and processed independently in order to obtain the TFFS and the correct fixing rate as evaluation indicators.
From the statistics of the ambiguity fixing, S1 and S2 have almost the same fixing performance for all the baselines and sessions and their fixing RATIO time series agree with each other. That means for the integrated searching mapping to ambiguities of longer wavelength does not improve the correct fixing rate, as LAMBDA method searches for the best integer combination for reliable fixing. We also found that the integer transformation matrices from LAMBDA method are not exactly the same for S1 and S2 which should be investigated further.
The cascading scheme S3 overperforms S1 and S2 for all baselines and sessions and its advantage is more significant over longer baselines due to the existence of inaccurately modelled biases and possibly also large noises caused by station environment. The RATIO of S3 of EWL and WL are several ten times larger than that of S1 and S2 and even the RATIO of S3 B1 fixing is also obviously larger than that of S1 and S2 because the float solution of S3 B1 is already improved by the fixed EWL and WL ambiguities.
The efficiency of BDS triple-frequency carrier phase ambiguity resolution is significantly better than that using dual-frequency observation data. For the medium-range baselines in the experiment, the TFFS is always shorter than 10 s on average. This could further extend inter-station distance of medium-range RTK, so that a single reference station could cover a larger area.
In terms of positioning results, the positioning accuracy, measured by the RMS of the position differences of the solutions using dual-frequency data and triple-frequency data is almost the same and the improvement is hardly visible.