Assessing the Value of Agroecosystem Services in Warmia and Mazury Province Using Choice Experiments

: Due to the often-intangible beneﬁts of agroecosystem services (AES), their cultural, social, and economic value is very difﬁcult to quantify and integrate into decision-making and policy-making processes. The aim of this study is to assess the non-market values of AES in an exceptionally environmentally rich area of the Warmia and Mazury region (Poland), identifying consumers’ preferences for them using the choice experiment (CE) method. Four AES attributes were selected for the research: (i) water quality; (ii) wildlife populations; and (iii) agricultural landscape. The study has revealed that the residents of Warmia and Mazury region were concerned about environmental issues that may be caused by agriculture. There was a demand for the provision of AES. The application of multinomial logit (MNL) model has revealed that marginal willingness to pay (MWTP) values were the highest for water quality (EUR 1.94), followed by wildlife population (EUR 1.02) and agricultural landscape (EUR 0.85). The ﬁndings have provided quantitative information related to the demand for improvements in AES through agri-environmental protection programs.


Introduction
Historically, agricultural systems have been primarily created for the production of food and fiber raw materials [1]. Nevertheless, they have a significant impact on the formation of cultural ecosystems and provision of different ecosystem services (ES) in rural areas: healthy and safe food, quality of surface water, air and soil, rich landscape, and biodiversity [2]. Agriculture is a managed or transformed or cultured ecosystem. Most food and other biological products are produced from agriculture as well as forestry, fishing, hunting, and harvesting activities [3]. Here, production is the use of man-made ecosystems and the use of natural ecosystem elements. In this context, according to Christiansen (1979) [3], agricultural systems are referred to as the agroecosystem, and the outputs delivered by them are referred to as the AES. The agroecosystem is the only ecosystem directly managed by people in order to meet own needs [4].
Therefore, farmers remain the key agroecosystem managers globally, and are responsible for the provision of multiple AES, in particular, food provisioning services. Besides the intensification of food production, it is very important to include farmers in the development of land management strategies [1]. On the other hand, farmers do not always have economic incentives to provide AES. When farmers receive income from the products sold, they receive no compensation for AES, so provision of AES may not be a part of their decision-making process.
Agriculture is an integral part of an agroecosystem; therefore, it is particularly important to analyze the functions of agroecosystems, reveal the concept of agroecosystem services, their social and economic benefits, and the need for special conditions for their provision.
Agriculture 2021, 11, 4 3 of 14 In this context, the aim of this study was to assess the non-market values of agroecosystem services in an exceptionally environmentally rich area of the Warmia and Mazury region (Poland), identifying consumers' preferences for them using the CE method. The paper is structured as follows: Section 2 describes the main characteristics of agriculture in Warmia and Mazury region; Section 3 explains the methodology applied in the research; the results of the empirical application are discussed in Section 4; and Section 5 ends the paper by presenting its main conclusions.

Case study: Warmia and Mazury
Warmia and Mazury province is located in the northeast part of Poland ( Figure 2) [16]. This area is one of the most attractive tourist, sightseeing, and recreational regions in Poland. Its population is almost 1.5 million. The capital of the region is Olsztyn. The Warmia and Mazury province is the fourth largest voivodeship in Poland (area over 24,000 km 2 ). The EU's Common Agricultural Policy has specific economic initiatives in order to ensure the provision of ecosystem services. It is focused on restoring, preserving, and enhancing ecosystems related to agriculture. The main measures are "Agri-Environment and Climate" measures, organic farming, payments to farmers in areas with handicaps (except for mountainous areas), greening and other environmental measures.
These agri-environmental programs are implemented by the means of monetary support from EU citizens who are also the consumers of such services. However, these measures are not always based on consumer expectations regarding agricultural activities and the AES they create. The analysis of consumer opinion and preferences towards AES creation and use should be included into the decision-making process of distribution of the agri-environmental payments. Literature [7,8] has suggested that the stated preference methods, in particular, choice experiments (CEs) are the best techniques of evaluation of the AES or agricultural public goods when analyzing both the overall product and each component or attribute.
Several researchers have analyzed consumer preferences for AES or multiple/several agricultural public goods in the EU countries [6,[9][10][11], focusing on particular elements or service, such as recreation [12] or agricultural landscape [13]. A few research studies were conducted in Poland, for example, Włodarczyk-Marciniak et al. (2020) implemented socio-cultural valuation of Polish agricultural landscape components by farmers [14], and Czyżewski (2020) valuated public goods from agricultural landscape [15]. However, no research was conducted about consumer preferences for AES in Poland or Polish regions. In this context, the aim of this study was to assess the non-market values of agroecosystem services in an exceptionally environmentally rich area of the Warmia and Mazury region (Poland), identifying consumers' preferences for them using the CE method. The paper is structured as follows: Section 2 describes the main characteristics of agriculture in Warmia and Mazury region; Section 3 explains the methodology applied in the research; the results of the empirical application are discussed in Section 4; and Section 5 ends the paper by presenting its main conclusions.

Case Study: Warmia and Mazury
Warmia and Mazury province is located in the northeast part of Poland ( Figure 2) [16]. This area is one of the most attractive tourist, sightseeing, and recreational regions in Poland.
Its population is almost 1.5 million. The capital of the region is Olsztyn. The Warmia and Mazury province is the fourth largest voivodeship in Poland (area over 24,000 km 2 ). Warmia and Mazury province is characterized by specific natural, economic, and social features. The most important of them include: -A relatively low population density per unit area (more than twice lower than the average in Poland); -Low pollution of the atmosphere of harmful emissions of gases and dust; -Lack of soil contamination with heavy metals and lack of potential premises for the occurrence of such threats on a larger scale; -Branched network of rivers and lakes, which causes a system of circulation of surface and deep ground water, requiring special protection, as opposed to elsewhere than in the country; -Unique and extremely beautiful landscape diversity, enabling the development of numerous forms of tourism, recreation, and treatment; -Diversified terrain configuration and soil mosaicism, justifying the use of appropriate directions of land use and maintaining the appropriate level of intensification of agricultural production, which results from the need to prevent erosion processes and leaching nutrients from the soil; -Low level of industrialization of the area, including lack of environmentally troublesome industry; -Relatively low share of people living from non-agricultural sources; -Favorable area structure of farms; -Significant area of forests [17].
Forest areas and lakes are characteristic for the nature of the region. Half of the ten largest lakes in Poland are located in Warmia and Mazury province, including the two largest ones: Śniardwy and Mamry. Forests in this climatic and geographical zone are the most natural formation. Woodiness in the Warmia and Mazury province was 31.4% (sixth in the country with an average for Poland-29.6%) in 2017 [17].
The leading branch of the economy of the voivodeship is agriculture. The intensity of agricultural production in Warmia and Mazury province is relatively small. Data from the Central Statistical Office show that the use of NPK fertilizers per 1 ha of arable land is at a relatively low level.
The crop structure of Warmia and Mazury province is dominated by cereals (65.5%) with the largest share being wheat, i.e., approx. one-quarter, then rape-7.2%, green maize-6.8%, legumes-2.8%, and root corn-2.4%. Agriculture and forestry account for almost 30% of the total water consumption (including most of the water in fishponds). The irrigated area of arable land and forest covers approximately three thousand hectares. Irrigation is implemented by means of the subsoil water. There is no sprinkling or flooding. The demand for water is estimated to increase; however, the risk of water shortages is not expected to increase significantly in Warmia and Mazury province. Warmia and Mazury province is characterized by specific natural, economic, and social features. The most important of them include: -A relatively low population density per unit area (more than twice lower than the average in Poland); -Low pollution of the atmosphere of harmful emissions of gases and dust; -Lack of soil contamination with heavy metals and lack of potential premises for the occurrence of such threats on a larger scale; -Branched network of rivers and lakes, which causes a system of circulation of surface and deep ground water, requiring special protection, as opposed to elsewhere than in the country; -Unique and extremely beautiful landscape diversity, enabling the development of numerous forms of tourism, recreation, and treatment; -Diversified terrain configuration and soil mosaicism, justifying the use of appropriate directions of land use and maintaining the appropriate level of intensification of agricultural production, which results from the need to prevent erosion processes and leaching nutrients from the soil; -Low level of industrialization of the area, including lack of environmentally troublesome industry; -Relatively low share of people living from non-agricultural sources; -Favorable area structure of farms; -Significant area of forests [17].
Forest areas and lakes are characteristic for the nature of the region. Half of the ten largest lakes in Poland are located in Warmia and Mazury province, including the two largest ones:Śniardwy and Mamry. Forests in this climatic and geographical zone are the most natural formation. Woodiness in the Warmia and Mazury province was 31.4% (sixth in the country with an average for Poland-29.6%) in 2017 [17].
The leading branch of the economy of the voivodeship is agriculture. The intensity of agricultural production in Warmia and Mazury province is relatively small. Data from the Central Statistical Office show that the use of NPK fertilizers per 1 ha of arable land is at a relatively low level.
The crop structure of Warmia and Mazury province is dominated by cereals (65.5%) with the largest share being wheat, i.e., approx. one-quarter, then rape-7.2%, green maize-6.8%, legumes-2.8%, and root corn-2.4%. Agriculture and forestry account for almost 30% of the total water consumption (including most of the water in fishponds). The irrigated area of arable land and forest covers approximately three thousand hectares. Irrigation is implemented by means of the subsoil water. There is no sprinkling or flooding. The demand for water is estimated to increase; however, the risk of water shortages is not expected to increase significantly in Warmia and Mazury province.
The area of Warmia and Mazury is one of the richest areas in Poland in terms of fauna. In addition to common animals, there are a lot of rare ones, including those that are at the risk of extinction.
Warmia and Mazury province has specific environmental conditions for economic development. About one-fifth of Poland's surface under water (including seawater) falls to Warmia and Mazury, which is the most in the country. Additionally, the share of this area in the total area of the province is the highest in Warmia and Mazury (about 6%). In a large part of the province municipalities, the share of surface waters in the area of the commune exceeds 10%, and forests occupy almost 10% of all such areas in the country, which is almost one-third of the province's area and is the sixth-most in Poland. Regions with special natural values which are legally protected occupy over 10% of these types of areas in Poland [18].
In addition, it should be noted that Warmia and Mazury region is an environmentally rich area, where agriculture is multifunctional, delivering multiple ecosystem services which are highly important for society. Therefore, the understanding and analysis of the preferences of residents of Warmia and Mazury concerning agroecosystem services will be useful for the improvement of agri-environmental policy in the attempt to change performance of farmers in term of sustainability. Payments to farmers have to reflect societal demand for agroecosystem services in Warmia and Mazury region.

Experimental Design
Lancaster's Theory of Value [19] and the Random Utility Theory [20] (Thurstone 1927) are the main theoretical basis of choice experiments. According to Lancaster's Theory, the utility for goods and services can be decomposed into individual utilities by their characteristics or attributes [9]. The function of the attributes of alternatives relevant to a given choice problem can model the choices of consumers [21]. Random Utility Theory explains the diversity of the opinions choosing the offered combinations. It is useful for analysis of the provision of attributes (and their levels) within goods and services and the relationships between these attributes [22]. The services subject to the estimation are decomposed into different attributes, which are presented at different levels in CE, thus creating different scenarios, which could be chosen by respondents during the survey [23]. In this study, CEs enabled valuation of different ES within an agroecosystem of Warmia and Mazury. Implementation of the CE for the respondents' willingness to pay, revealed as implicit price for estimation of the agroecosystem services, consists of four main stages, i.e., selection of the attributes of AES and their levels, construction of the experimental design, questionnaire design, and data collection. These stages are presented below.
The first and one of the essential stages of CEs is the selection of attributes suitable for the research. Following Novikova (2016) [24] and the analysis of environmental issues of agriculture in Warmia and Mazury [25], the attributes were selected for the present research. Lithuania and Warmia and Mazury were identified as having a similar geographical location, environmental conditions, and issues, suggesting the assumption that the attributes could be feasible for Warmia and Mazury case study. The three selected AES attributes were substantiated by extensive review of the previous research.
Each of the attributes was set to a different level. The lowest level corresponded to the status quo (current), the second level-10% improvement in AES, and the highest level-20% improvement achieved in AES (i.e., the best possible performance scenario). The last monetary attribute was personal contribution into environmental management, which was proposed to be considered on an annual basis, for a five-year period ( Table 1). The selected attributes and their levels were checked by the focus group and during the pilot survey. The survey was divided into the set of alternative options, consisting of two alternative options (option A and option B), representing improved situations and status quo situation. The unlabeled experimental design, representing combination of attributes, was used. Following Colombo et al. [30], definitions of the levels were presented using percentage.
D-efficient fractional factorial design excluding unrealistic cases was adapted to each of the choice questions to make the statistically efficient choice design for the main survey. Thirty-six choice sets were created with linear D-optimal using SAS© software, which were then distributed into four blocks. For 1-3 bocks four cards, and for the fourth block five cards were proposed. Table 2 illustrates an example of the choice cards shown to the respondents. be considered on an annual basis, for a five-year period ( Table 1). The selected attributes and their levels were checked by the focus group and during the pilot survey. The survey was divided into the set of alternative options, consisting of two alternative options (option A and option B), representing improved situations and status quo situation. The unlabeled experimental design, representing combination of attributes, was used. Following Colombo et al. [30], definitions of the levels were presented using percentage.
D-efficient fractional factorial design excluding unrealistic cases was adapted to each of the choice questions to make the statistically efficient choice design for the main survey. Thirty-six choice sets were created with linear D-optimal using SAS© software, which were then distributed into four blocks. For 1-3 bocks four cards, and for the fourth block five cards were proposed. Table 2 illustrates an example of the choice cards shown to the respondents. After the questionnaire had been developed, its feasibility was tested by the pilot survey in June-July, 2019. The results of the pilot survey showed that the selected attributes, as well as their levels and questions were understood and accepted by the respondents. The final questionnaire consisted of the following three parts: (1) questions about respondents' knowledge and attitudes on farming impacts on environment; (2) choice cards with different agri-environmental options for identification of the respondents' preferences; (3) questions related to respondent socioeconomic characteristics.

Modeling Framework
Agriculture 2021, 11, x FOR PEER REVIEW be considered on an annual basis, for a five-year period ( Table 1). The selected their levels were checked by the focus group and during the pilot survey. The survey was divided into the set of alternative options, consisting o tive options (option A and option B), representing improved situations an situation. The unlabeled experimental design, representing combination of a used. Following Colombo et al. [30], definitions of the levels were present centage.
D-efficient fractional factorial design excluding unrealistic cases was ad of the choice questions to make the statistically efficient choice design for the Thirty-six choice sets were created with linear D-optimal using SAS© sof were then distributed into four blocks. For 1-3 bocks four cards, and for the five cards were proposed. Table 2 illustrates an example of the choice cards respondents. After the questionnaire had been developed, its feasibility was tested survey in June-July, 2019. The results of the pilot survey showed that the s utes, as well as their levels and questions were understood and accepted by ents. The final questionnaire consisted of the following three parts: (1) qu respondents' knowledge and attitudes on farming impacts on environme cards with different agri-environmental options for identification of the preferences; (3) questions related to respondent socioeconomic characteristi

Modeling Framework
Agriculture 2021, 11, x FOR PEER REVIEW be considered on an annual basis, for a five-year period (T their levels were checked by the focus group and during th The survey was divided into the set of alternative o tive options (option A and option B), representing imp situation. The unlabeled experimental design, representi used. Following Colombo et al. [30], definitions of the centage.
D-efficient fractional factorial design excluding unr of the choice questions to make the statistically efficient c Thirty-six choice sets were created with linear D-optim were then distributed into four blocks. For 1-3 bocks fo five cards were proposed. Table 2 illustrates an example respondents. After the questionnaire had been developed, its f survey in June-July, 2019. The results of the pilot surve utes, as well as their levels and questions were understo ents. The final questionnaire consisted of the following respondents' knowledge and attitudes on farming imp cards with different agri-environmental options for id preferences; (3) questions related to respondent socioeco

Modeling Framework
After the questionnaire had been developed, its feasibility was tested by the pilot survey in June-July, 2019. The results of the pilot survey showed that the selected attributes, as well as their levels and questions were understood and accepted by the respondents. The final questionnaire consisted of the following three parts: (1) questions about respondents' knowledge and attitudes on farming impacts on environment; (2) choice cards with different agri-environmental options for identification of the respondents' preferences; (3) questions related to respondent socioeconomic characteristics.

Modeling Framework
According to the Random Utility Theory, the person chooses the alternative, which gives the highest utility. People choose among alternatives on the basis of the utility function with two components: a systematic (i.e., observable) component plus a random term (non-observable by the researcher) within this theoretical framework [31]: where U in is the utility provided by alternative i to subject n, V in is the systematic component of the utility, Z i is the vector of attributes of alternative i, S n is the vector of socioeconomic characteristics of the respondent n, and ε is the random term. The econometric analysis is based on a multinomial logit (MNL) model. The MNL is one of the available probabilistic choice models, mostly used in CEs [9,13,31,32]. According to the MNL model, the probability that an individual n will choose alternative i (P in ) among other alternatives j (j = 1 . . . J) of a set C n is expressed by the equation [31]: where V in is the systematic component of the utility provided by alternative i, and µ is a scale parameter which is inversely proportional to the standard deviation of the error terms and usually is assumed to be equal to one [32]. The maximum likelihood method could be used to evaluate the robustness of the econometric model [33]. The willingness to pay (WTP) is calculated by estimating the marginal rates of substitution (MRS) between attributes. The implicit price is calculated as follows: where ∆β i is the variation in utility due to the change in the level of the attribute i, while β bid is the estimation of the coefficient corresponding to the attribute, expressed in monetary terms [34].

Sampling Characteristics
The final survey was performed in September-November 2019; 550 questionnaires were distributed. We received 353 valid questionnaires, excluding protest responses and imperfectly filled-in questionnaires. The valid questionnaires were obtained, delivering 1505 choice observations. Warmia and Mazury citizens over 18 years old comprised the target population. The random sampling technique was applied during the survey. The survey was implemented using a face-to-face method, randomly selecting the respondents at different events, courses, and public places. The descriptive statistics of the respondents' main sociodemographic characteristics are presented in Table 3. The largest part of our sample consisted of young people aged 18-39. This did not correspond to the actual age distribution in Warmia and Mazury. Women were slightly overrepresented; the distribution of males and females in Warmia and Mazury is almost equal. Slightly more than half of the respondents were living in urban areas, and the remaining share in rural areas, which was almost in line with the Warmia and Mazury region data available. Therefore, our sample included more people of young age, with approximately EUR 460 monthly earnings. The main cause of this situation was a very high unwillingness of older people to participate in the survey. Unfortunately, the research was not provided with any additional/external funding to pay people for participation in the survey in order to ensure the representativeness of the sample. The respondents with secondary and fundamental education comprised the largest share of the sample, and the respondents with comparative and special secondary education accounted for the smallest share in the sample, which was almost in line with the distribution of the Warmia and Mazury population. The present study was primarily a demonstration study, and the first attempt to evaluate consumer preferences for AES in Warmia and Mazury. Therefore, the present paper does not consider reaching a study sample which could meet the general population characteristic to be fundamental for its purpose. It is a demonstration study. Table 3 shows the attitudes and awareness of the agricultural impacts on the environment among the residents of Warmia and Mazury. The majority (about 65%) of the respondents found the impacts of farming on natural environment positive and very positive. About 30% of them stated that impacts of agriculture in Warmia and Mazury were bad. Respondents were asked about their feelings concerning the environmental aspects from agriculture in Warmia and Mazury with the focus on water quality in rivers and lakes, drinking water quality, decline in flora populations, decline in fauna populations, and formation and maintenance of landscape. More than 40% of the respondents stated that they were concerned (always or very often) about water quality in the rivers and lakes. About 80% of the respondents were very concerned (always or very often) about drinking water quality. Respondents were less concerned about loss in biodiversity and agricultural landscape issues. Respondents were asked about the awareness of positive and negatives issues from agriculture. Most of them claimed that they were aware of the issues related to farming and agriculture in Warmia and Mazury (Table 4).

Modeling Results
Two MNL models were estimated using the data derived from the survey with Nlogit 6. The first model (Model I without additional variables) was a basic specification showing the importance of the choice attributes in explaining consumer preferences of the different options of environmental protection programs. The second model (Model II or with additional variables) refers to socioeconomics and knowledge about the impact of agricultural activity on environment variables in addition to the attributes of AES in the choice set. In Model I, the utility is determined by the levels of four attributes (water quality, wildlife populations, agricultural landscape, and personal contribution) in the choice sets. Thereby, Model I provides an estimate of the effect of a change in any of the attributes on the probability that the respondent chooses one of these options. Literature suggests that additional characteristics about the respondents' reveals their preferences and explains their choices better [7]. Therefore, additional socioeconomic and knowledge characteristics were included into Model II through the interaction with the water, landscape and wildlife attributes. Five additional characteristics were found to be statistically significant. Hence, they were included into Model II: − through the interaction with water quality: respondents' awareness about the environmental issues from agriculture and area of residence; − through the interaction with wildlife populations and income; − through the interaction with agricultural landscape: area of residence and gender (Table 5). Other socioeconomic characteristics and knowledge about the impact of agricultural activity on environment were not statistically significant. Therefore, they had no impact on respondents' choices and preferences towards the water quality, wildlife populations and agricultural landscape. The estimation results of the MNL model with and without additional variables are shown in the Table 6. The signs of the coefficients for all attributes in both models are as expected and statistically significant at 1% level. Water, wildlife, and landscape coefficients are positive, meaning improvement of these attributes leads to the gained utility of the respondents. Negative price coefficients suggest that the respondents were willing to accept the policy with lower cost to them.
The MNL results show that the socioeconomic characteristics and knowledge on the impact of agricultural activity variables influence the respondents' choices and preferences: (i) negative and significant coefficient of water quality and residence interaction suggests that the residents from rural areas would be less inclined to choose improvements in the water quality because of the agricultural activity compared to the residents from urban areas. The relationship between the agricultural landscape and place of residence has shown the opposite situation, when people from rural areas are more likely to choose the agri-ecosystem compared to those from urban areas; (ii) negative and significant coefficient of the water quality and gender relationship indicates higher willingness to pay for drinking water quality by women compared to men (men would choose improvements in the water quality more rarely); (iii) positive and significant sign of the relationship between wildlife and resident income suggests that the respondents gaining more income would choose the program for improvement of wildlife populations more frequently. Based on the MNL model II, the implicit prices/WTP for the attributes analyzed were calculated, as described in Equation (3). According to the results, the respondents' MWTP/person for water quality was EUR 1.94 for the reduction in underground water pollution due to agricultural activity by 1%; EUR 1.02 for improvement for protecting the diversity of wildlife; and EUR 0.85 for an improved agricultural landscape. The trade-off between attributes indicated that Warmia and Mazury region residents gave preference to water quality, while wildlife populations and agricultural landscape attained lower preferences.
Results for the WTP estimates for the different levels of provision considering the attributes (10% and 20% improvements) are presented in Table 7. Consumers' willingness to pay substantiates their higher concerns about the drinking water quality. It shows that residents of Warmia and Mazury were mostly concerned about the primary needs. The analysis of their opinions has shown that biodiversity and landscape services were less significant for them, which was substantiated with the lowest willingness to pay for these services.

Discussion
Findings of the current research about consumer preferences for AES are in line with the previous studies. Research has shown that the water quality is the most important attribute for consumers [11,30,34,35], which is in line with our findings. The respondents ranked biodiversity second by importance, as found previously by Faccioni et al. [11]. As highlighted by Soini and Aakkula [36], the main reason is that people could see and feel the direct effect of water quality on their welfare. They are not able to estimate the consequences of loss in biodiversity, because the overall concept and importance of biodiversity for people is not usually understood [6]. Being unconscious of the concept biodiversity is due to its intangibility [2,6,36]. Our research has revealed that the agricultural landscape the least important for the respondents, thus substantiating the conclusions of the previous research [11,26,30].
Bernués et al. [6], Alcon et al. [37] and Rodríguez-Ortega et al. [38] obtained biodiversity and landscape values higher than the results obtained by current research. This happened due to the lower of income level of Warmia and Mazury region residents. Similar results of willingness to pay were obtained by Novikova [24], who assessed willingness to pay for AES in Lithuania, where the level of subsistence was more similar to Poland.
Our results underline the importance of AES depending on respondents' socioeconomic characteristics and individuals' backgrounds and interests, the same as in the research by Faccioni et al. [11]. For example, women have a higher willingness to pay for AES than men do. These findings are in line with other research [38]. Similarly to other studies [34], the present research has revealed that higher income respondents have higher willingness to pay for agroecosystem services. Arriaza et al. [9] and Baskaran et al. [34] have found out that respondents from urban and rural areas have different opinions about agroecosystem services; urban residents show a higher willingness to pay than those from rural areas. The present research has found that urban residents were willing to pay less for drinking water quality than rural residents. Respondents who had more knowledge about the impact of agriculture on nature were willing to pay less for AES (drinking water quality in current study), as shown by other studies [13,24,28,30].

Conclusions
The present study was the first instance of CE application in the Warmia and Mazury region. The focus was to analyze the local residents' preferences towards the provision of AES. According to the characteristics of the Warmia and Mazury province, it is an agricultural region with high natural values. The obtained research results may have practical application in shaping agricultural policy at the regional level, in creating smart specializations by including them in strategic documents concerning the region.
The study has clearly shown that residents of Warmia and Mazury region were concerned about environmental problems caused by agriculture. The results of the study have revealed that, in environmentally rich areas such as Warmia and Mazury region, agriculture is multifunctional, delivering multiple AES highly important for the society. Particularly strong support was claimed towards the provision of AES, such as water quality and wildlife populations.
Three AES (water quality, wildlife populations, and agricultural landscape) were selected for evaluation for the Warmia and Mazury region. The choice was made according the analysis of the situation in the region and extensive literature research about the impacts of agricultural activity. In exploration of the choice experiment method, the value was calculated and expressed as average personal contribution (EUR per year for the next five years), expressed in implicit prices. The study has revealed that the socioeconomic characteristics and knowledge about agricultural activity impact influence consumers' choices and the WTP for AES. Women, higher income residents and residents from rural areas are more willing to pay for AES, in particular, for drinking water quality. The present study results are in line with similar research performed by other scientists. However, the values of WTP were lower compared to other studies. These results may have been influenced by lower income level of Warmia and Mazury region residents.
It should be noted that the present study has limitations, because it covered only three attributes (water, biodiversity and landscape). Current agri-environmental policy measures cover a lot more ES. Therefore, an even more detailed study could be carried out in the future, the results of which would be applied for the development of agrienvironmental policy measures, or as suggested by Niedermayr et al. (2018) [39], could aim to evaluate the willingness of farmers to accept for the implementation of environmentally friendly farming practice