North-eastern Poland Research Articles

Long-Term Changes in the Thermal and Ice Regime of the Biebrza River (Northeastern Poland) in the Era of Global Warming

In the context of ongoing environmental changes, particularly against the backdrop of global warming, significant attention is being given to areas of exceptional natural value that, in many aspects, retain a pristine character. One such area is the Biebrza River in northeastern Poland, which, together with the wetlands in its basin, forms one of the most valuable ecosystems of its kind in Europe. This study analyses the changes in the thermal and ice regime for two hydrological stations, Sztabin and Burzyn, in the period from 1959 to 2023. It was found that the average annual water temperature in this period for the Biebrza River increased by 0.28 °C/decade, and in the case of ice phenomena, statistically significant changes for both stations showed a decline, with an acceleration of the ice cover disappearance by an average of 3 days/decade. These recorded changes should be considered unfavourable, as they will affect the transformation of both the biotic and abiotic characteristics of the river itself, as well as the natural elements associated with it.

Does Every Strain of Pseudomonas aeruginosa Attack the Same? Results of a Study of the Prevalence of Virulence Factors of Strains Obtained from Different Animal Species in Northeastern Poland

Background: Pseudomonas aeruginosa is a pathogen that causes infections in animals and humans, with veterinary implications including ear infections in dogs, respiratory diseases in cats, and mastitis in ruminants. In humans, it causes severe hospital-acquired infections, particularly in immunosuppressed patients. This study aimed to identify and assess the prevalence of specific virulence factors in Pseudomonas aeruginosa isolates. Methods: We analyzed 98 Pseudomonas aeruginosa isolates from various animal samples (dogs, cats, ruminants, fowl) from northeastern Poland in 2019–2022 for virulence-related genes (toxA, exoU, exoT, exoS, lasB, plcN, plcH, pldA, aprA, gacA, algD, pelA, endA, and oprF) by PCR and assessed biofilm formation at 48 and 72 h. Genomic diversity was assessed by ERIC-PCR. Results: The obtained results showed that all strains harbored the pelA gene (100%), while the lowest prevalence was found for pldA (24%) and exoU (36%). Regardless of the animal species, strong biofilm forming ability was prevalent among the strains after both 48 h (75%) and 72 h (74%). We obtained as many as 87 different genotyping profiles, where the dominant one was profile ERIC-48, observed in four strains. Conclusions: No correlation was found between presence or absence of determined genes and the nature of infection. Similarly, no correlation was found between biofilm-forming genes and biofilm strength. The high genetic diversity indicates challenges for effective prevention, emphasizing the need for ongoing monitoring and research.

Spatial patterns of vascular plant species richness in Poland: Relations among species group richness and hot spot distribution

Knowledge of spatial patterns of species richness (SR) is highly relevant for theoretical research in ecology and the development of conservation plans. In Poland, despite a long tradition of botanical surveys, vascular plant SR has not been mapped, nor have the correlations in richness among different plant species groups been explored at the entire country scale. Here, we used a recently published data set to examine spatial patterns and relationships among the joined SR of vascular plant species, including native species, archaeophytes, neophytes, and species with high conservation value (i.e., red-list species). The basic spatial unit employed was a 10 × 10 km grid from the Atlas of Distribution of Vascular Plants in Poland (ATPOL). We found that SR of native species and archaeophytes, neophytes, and red-list species were positively correlated. The main patterns of SR and the percentage of particular groups in the joined SR were based on three components: (1) gradient of overall SR, (2) invasion level, and (3) peculiarity of flora in some regions resulting from the high number and proportion of rare species that often have high conservation value (red-list species). In general, northeastern Poland was species-poor, while the Carpathian Mountain range, the uplands in southern Poland, and some parts of Wisła River valley had the highest SR concentrations. e location of SR hotspots usually did not overlap with the existing national parks system. The correlations among native SR, high conservation value species, and neophyte SR suggest that biological invasions are among the most important threats to vascular plant diversity in Poland. Finally, despite likely biases in SR assessments, we demonstrated that the presented maps seem to reflect general ecological gradients influencing vascular plant distribution in Poland.

Role of Lake Morphometric and Environmental Drivers of Ice Cover Formation and Occurrence on Temperate Lakes: A Case Study from the Eastern Baltic Lakeland, Poland

The presence of ice cover on temperate lakes is a crucial factor in determining the functioning of these ecosystems. The isolation of water from atmospheric influences significantly alters physical, chemical, and biological processes, and the intensity of this impact depends on the duration of the ice cover. This study analyzed the basic parameters of ice cover on several dozen lakes in Northeastern Poland. The aim of this study is to investigate the influence of morphometric parameters, alongside environmental factors, on the variation of ice cover characteristics in lakes located within the Eastern Baltic Lakeland. Characterization of ice conditions in the analyzed lakes was based on basic statistics such as minimum and maximum values, mean, standard deviation, coefficients of variation, skewness, and kurtosis. Given that the dataset contains variables describing ice phenomena in the studied lakes and data describing location, morphometric parameters, and land cover directly adjacent to the lake (treated as independent variables), a method of Spearman’s rank correlations and constrained ordination method were decided upon. Despite the relatively small study area, significant variability was observed, with average differences as follows: 26 days for the onset of ice cover, 17 days for the end date, 15 cm for ice thickness, and a 30-day difference in the average duration of ice cover. Key factors included parameters such as lake volume, average depth, and land use (urbanized and agricultural areas). Understanding parameters such as the onset and end of ice cover is essential for lake ecosystems, both from an ecological and economic perspective. This knowledge is crucial for interpreting the behavior of living organisms, water quality, and economic considerations.

Traffic Intensity as a Factor Influencing Microplastic and Tire Wear Particle Pollution in Snow Accumulated on Urban Roads

Traffic-related roads are an underestimated source of synthetic particles in the environment. This study investigated the impact of traffic volume on microplastic (MP) and tire wear particle (TWP) pollution in road snow. An examination was conducted in a medium-sized city situated in northeastern Poland, known for being one of the cleanest regions in the country. MPs and TWPs were found at all 54 sites, regardless of the intensity of traffic. The average concentration for all samples was 354.72 pcs/L. Statistically significant differences were found between the average values of the particle concentration on low, medium, and heavy traffic roads, amounting to 62.32 pcs/L, 335.97 pcs/L, and 792.76 pcs/L, respectively. Within all three studied groups of roads, MPs and TWPs with the smallest size, ranging from 50 to 200 μm, were prevalent. In all of the studied groups of roads, four analyzed shapes of particles were found, with irregular fragments being the most abundant form (89.23%). The most frequently recorded color among the collected samples was black (99.85%), and the least frequently recorded color was blue, constituting only 0.01%. This study suggests that snow cover on the roads may act like a temporary storage of pollutants during winter particularly in the temperate climate zone and, after thawing can significantly increase the concentration of MPs and TWPs in surface waters. Possible measures to decrease the release of MPs and TWPs into the environment in the city may include reducing the traffic volume and speed, implementing street sweeping, utilizing filtration chambers, and installing stormwater bioretention systems or settling ponds.

Prevalence of Aleutian Mink Disease Virus (AMDV) in Free-Ranging American Mink from Biebrza and Narew National Parks (Poland)-An Epidemiological Concern.

Aleutian Mink Disease Virus (AMDV) is the causative agent of Aleutian disease (AD). This progressive and chronic disorder significantly impacts the mink breeding industry, affecting farmed and free-ranging American and European mink. This study investigated AMDV variants isolated from free-ranging American mink in northeastern Poland. Between 2018 and 2019, 26 spleen samples were collected from mink in Narew National Park (NNP) and Biebrza National Park (BNP). DNA was extracted and subjected to PCR to amplify the NS1 gene, followed by sequencing and phylogenetic analysis. The NS1 gene was detected in 50% of samples from NNP minks and in 30% of samples from BNP minks, with an overall prevalence of 42.31%; these findings align with global data and indicate serious ecological and health concerns. Ten closely related AMDV variants and one distinct variant were identified. The grouped variants exhibited high genetic homogeneity, closely related to strains found in mink from the USA, Germany, Greece, Latvia, and Poland; meanwhile, the distinct variant showed similarities to strains found in mink from Finland, Denmark, China, Poland, and Latvia, suggesting multiple infection sources. These findings, consistent with data from Polish mink farms, indicate significant genetic similarity between farmed and wild mink strains, suggesting potential bidirectional transmission. This underscores the importance of a One Health approach, emphasizing the interconnectedness of human, animal, and environmental health. Continuous surveillance and genetic studies are crucial for understanding AMDV dynamics and mitigating their impacts. Measures to reduce transmission between farmed and wild mink populations are vital for maintaining mink health and ecosystem stability.

Examining the Sensitivity of Satellite-Derived Vegetation Indices to Plant Drought Stress in Grasslands in Poland.

In this study, the emphasis is on assessing how satellite-derived vegetation indices respond to drought stress characterized by meteorological observations. This study aimed to understand the dynamics of grassland vegetation and assess the impact of drought in the Wielkopolskie (PL41) and Podlaskie (PL84) regions of Poland. Spatial and temporal characteristics of grassland dynamics regarding drought occurrences from 2020 to 2023 were examined. Pearson correlation coefficients with standard errors were used to analyze vegetation indices, including NDVI, NDII, NDWI, and NDDI, in response to drought, characterized by the meteorological parameter the Hydrothermal Coefficient of Selyaninov (HTC), along with ground-based soil moisture measurements (SM). Among the vegetation indices studied, NDDI showed the strongest correlations with HTC at r = -0.75, R2 = 0.56, RMSE = 1.58, and SM at r = -0.82, R2 = 0.67, and RMSE = 16.33. The results indicated drought severity in 2023 within grassland fields in Wielkopolskie. Spatial-temporal analysis of NDDI revealed that approximately 50% of fields were at risk of drought during the initial decades of the growing season in 2023. Drought conditions intensified, notably in western Poland, while grasslands in northeastern Poland showed resilience to drought. These findings provide valuable insights for individual farmers through web and mobile applications, assisting in the development of strategies to mitigate the adverse effects of drought on grasslands and thereby reduce associated losses.

Petrogenesis of microgranular enclaves in the A-type granitoid Krasnopol intrusion (Mazury Complex, northeastern Poland): Evidence of magma mixing

The origin of magmatic microgranular enclaves has been investigated in the Mesoproterozoic granitoid Krasnopol intrusion (1.5 Ga), part of the AMCG (anorthosite–mangerite–charnockite–granite) Mazury Complex in the East European Craton (NE Poland). The granitoids are ferroan and metaluminous, and display the typical characteristics of A-type granites, with high contents of Zr, Nb, Ga and rare earth elements (REEs). The enclaves are metaluminous and have a broad compositional range with two groups distinguished: silica-poor (45–50 wt% SiO2) and silica-rich (54 to 59 wt% SiO2), the latter overlapping in composition with the granitoid samples. The silica-poor enclaves are enriched in REEs compared to the silica-rich type, while the silica-rich enclaves exhibit trace-element patterns similar to those of the granitoids. Initial whole rock εNd values range between -3.8 and -4.0 for the granitoids and give a slightly wider range from -2.6 to -3.8 for the enclaves. The 87Sr/86Sr initial values vary from 0.7084 to 0.7138 for the granitoids and between 0.7052 and 0.7075 for the enclaves and indicate that the granitoids and enclaves are not isotopically identical. These may suggest that the two magmatic systems represented by the granitoid host rock and the enclaves, were probably derived from different sources, but with sufficient interaction, which led to a progressive change in the composition of the enclaves towards intermediate composition. We suggest that the mafic melts of the enclaves were generated at the base of the thickened crust through partial melting of the lower crustal source, with a significant contribution from mantle material. The increase in temperature resulted in anatexis of the lower crust and the formation of the granitoid parental magma.

The Quality of Jerusalem Artichoke Biomass Harvested Twice during the Growing Season in North-Eastern Poland

The Quality of Jerusalem Artichoke Biomass Harvested Twice during the Growing Season in North-Eastern Poland

Changes in Hospital Admissions of Patients with COPD in Poland: A 14-Year Nationwide Analysis (2006–2019)

Poland is a country with a high burden of COPD and its risk factors. This population-wide and geospatial study on COPD in Poland, aimed to analyze changes in hospitalizations of patients with a diagnosis of COPD between 2006 and 2019 as well as to identify changes in the demographic characteristics of patients hospitalized with COPD. This is a retrospective analysis of hospital discharge records of patients with a diagnosis of COPD (ICD-10 code J44 and age 40 and over), hospitalized in Poland between 2006 and 2019. Data were analyzed separately for all patients hospitalized with COPD (primary and secondary diagnosis), for patients admitted due to COPD (primary diagnosis) and patients hospitalized with COPD as a comorbidity (secondary diagnosis). Between 2006 and 2019, a total of 1,663,420 hospital admissions of patients with COPD were reported. Between 2006 and 2019, the annual number of patients hospitalized due to COPD decreased by 53.8%, including 57.3% among men and 46.6% among women. The number of patients hospitalized with COPD as comorbidity increased by 35.8%. The lowest hospital admissions rate was in north-western Poland and the highest in north-eastern Poland. Among patients admitted due to COPD, the percentage of women increased from 32.6% to 37.7%. The percentage of patients admitted due to COPD and living in rural areas decreased from 51.2% in 2006 to 40.8%in 2019. This study provided data on changes in sociodemographic characteristics of hospitalizations of patients with a diagnosis of COPD, including overall reduction in COPD admissions but an increase in COPD burden among women.

Transformation of the Cultural Landscape in the Central Part of the Historical Region of Warmia in Poland

The study analyzed the cultural landscape of Dobre Miasto municipality, which is situated in north-eastern Poland, in the heart of the Warmia region. Cultural traditions that have been passed down from generation to generation and are presently being revived constitute the backbone of Warmia’s cultural heritage. The Warmian landscape comprises several distinctive components, including traditional settlement patterns dating back to the 13th and 14th centuries, as well as cultural values (historical and architectural values, and the level of their preservation) that are being increasingly recognized and contribute to the region’s tourism potential. The growing awareness that the cultural landscape must be preserved for the future generations implies that landscape protection is no longer regarded as a niche problem that is of interest only to experts in the field, and that activities should also be initiated at the local level. The study aimed to specify measures that should be undertaken to protect the components of the cultural landscape and to propose strategies that will not only contribute to the appreciation of cultural values, but will also facilitate the promotion and preservation of cultural landscape components that are gradually disappearing from the rural fabric. This complex approach to the rational management of environmental resources can set an example for other lake regions with clearly defined historical roots and well-preserved cultural values. The resulting guidelines and conclusions can be useful for developing new planning documents and updating and validating the existing legal regulations on planning, development, land use, and protection of the cultural landscape.

Spelt in Diversified and Spelt-Based Crop Rotations: Grain Yield and Technological and Nutritional Quality

A properly designed crop rotation contributes to the equilibrium of the agro-ecosystem and the volume and quality of the yield. The cultivation of spelt in crop rotations enriches its biodiversity and provides grains with many different types of nutritional value. The aim of this current study was to investigate how the distribution of winter spelt in different positions and after different forecrops in four-field crop rotations would affect the technological quality of the grain, the nutrient content of the grain, and the grain yield. A 6-year field experiment, designed in a randomised block, was conducted from 2012 to 2018 in north-eastern Poland (53°35′47″ N, 19°51′20″ E). This study provides the results from a 6-year (2013–2018) field experiment. The spelt was cultivated in four crop rotations: CR1—winter rape + catch crop (blue tansy), spring barley, field pea and winter spelt; CR2,—winter rape, winter spelt + catch crop (blue tansy), field pea and winter spelt; CR3—winter rape + catch crop (blue tansy), field pea, winter spelt and winter spelt; and CR4—winter rape, winter spelt + catch crop (blue tansy), spring barley and winter spelt. This study evaluated grain yield and the following grain parameters: the total protein, wet gluten and starch contents, the Zeleny index, the falling number, the weight of 1000 grains, the N, P, K, Mg, Ca, Cu, Fe, Zn and Mn contents, and the grain yield. The results were assessed at the significance level p < 0.05. It was demonstrated that the cultivation of spelt in all four crop rotations after winter rape and after field pea was characterised by higher protein and wet gluten contents, Zeleny index value and falling number, a greater weight of 1000 grains, higher N, P, Fe and Zn contents, and greater grain yield than those harvested from the crop rotations CR3 and CR4 after spelt and after barley. It was demonstrated that the cultivation of spelt in crop rotations CR3 and CR4, in succession after spelt and after barley, caused deterioration in grain quality (lower protein and gluten contents, a lower Zeleny index value, a lower falling number, and a smaller weight of 1000 grains, and the N, P, Fe and Zn contents). In addition, a smaller grain yield was obtained from these crop rotation fields. Regardless of the type of crop rotation, the cultivation of spelt after winter rape and after pea produced a high yield and a good quality yield of this cereal. Due to the lower yield of grain and its lower quality, it is not recommended that winter spelt is grown after each other or after spring barley.

The influence of forest site types on the distribution of moose Alces alces in north-eastern Poland

ABSTRACT Movements of moose are caused by variations in resource availability such as food, shelter and the energetic cost of foraging. These variables are important in forest and population management. Therefore, the main objective of the present study was to estimate the influence of forest site types on the distribution of the non hunted moose population in February and March 2020 in the Augustów Forest, a 1102.0 km2 forest complex in north-eastern Poland. In the study area, 20 sampling plots were distributed over a total area of 91.4 km2, where six teams, each consisting of 30 people, observed moose to estimate population numbers. The average population density was 1.22 ± 0.19 moose per km2 and ranged from 0.21 to 3.40 animals/km2. The total population number estimated this way amounted to 1344 animals. Non-linear regression demonstrated a positive relationship between the percentage of hygric-wet forest site types in the sampling plot area and the moose population distribution (r = 0.735, n = 20, p = 0.0096). In winter, these habitats are a critical refuge for moose. For this reason, tree harvesting and hunting in hygric-wet forest site types should be limited during this period.

Loose human bones as evidence of the multi-step burial rite: Case study of the Stone Age hunter-gatherer sites at Dudka and Szczepanki, Masuria (northeastern Poland)

Loose human bones as evidence of the multi-step burial rite: Case study of the Stone Age hunter-gatherer sites at Dudka and Szczepanki, Masuria (northeastern Poland)

Assessment of the impact of setting the heating curve on reducing gas consumption in a residential building while ensuring thermal comfort

Due to the large share of the residential building sector in CO2 emissions, it is necessary to look for quick-to-implement and cost-free solutions for heating systems aimed at increasing their efficiency and reducing fuel consumption. The aim of the study is a multi-aspect assessment of the proper setting of the heating curve in the weather controller of a low-temperature floor heating system with a condensing gas boiler, and an analysis of the impact of the heating curve slope on the natural gas consumption, the emission of gaseous and dust pollutants and on maintaining appropriate indoor air parameters: air temperature, operative temperature and PMV and PPD. A wide range of simulation tests were carried out using the TRNSYS software for an existing energy-efficient single-family residential building, with mechanical ventilation with heat recovery, located in north-eastern Poland in Bialystok. Changing the heating curve slope by 0.1 can reduce gas consumption and pollutant emissions by 1.6–10.3 %. Low heating curve slopes may result in failure to reach the temperature set by users (difference between the set and achieved indoor temperature up to 1.9 K), while maintaining standard thermal comfort parameters (−0.5<PMV<+0.5; PPD<10 %). Reducing the indoor air temperature by 2 K can contribute to reducing gas consumption by 11 % while maintaining standard thermal comfort parameters. Changing the heating curve slope is one of the effective methods of reducing fuel consumption in the building. Results are useful for users of buildings equipped with low-temperature heating systems with floor heating and installers of these systems.

Quality of Winter Wheat Flour from Different Sowing and Nitrogen Management Strategies: A Case Study in Northeastern Poland

Quality of Winter Wheat Flour from Different Sowing and Nitrogen Management Strategies: A Case Study in Northeastern Poland

Determination of the rabies virus variants circulating in Poland in 2021-2023 and their phylogeny with analysis of the strains in the Mazowieckie and Podkarpackie voivodeships.

Rabies is endemic in Europe and red foxes are the vector and reservoir of the rabies virus (RABV). Based on classification established in the early 1990s, four variants of the rabies virus have been distinguished in Europe. Rabies broke out in January 2021 in the Mazowieckie voivodeship in central north-eastern Poland. The virus spread rapidly, reaching the Świętokrzyskie voivodeship in the central southern part and the Lubelskie voivodeship in the eastern part in the next months. Nine rabies cases were reported in the Podkarpackie voivodeship in south-eastern Poland between 2021 and 2023, mainly in red foxes but also in dogs and wildcat. The aim of the study was the identification of RABV variants in wildlife and domestic animals in Poland between 2021 and 2023. The study involved 157 animal brains tested positive for rabies using a fluorescent antibody test. From 10% w/v brain homogenates, RNA was isolated and full-length RABV genomes were high-throughput sequenced with an RABV-enriched approach. Complete genomes of RABV isolates were phylogenetically analysed and the variants were estimated. Molecular and phylogenetic studies revealed 147 (93.6%) of the RABV strains out of 157 which had rapidly spread in the wildlife of the Mazowieckie, Świętokrzyskie and Lubelskie voivodeships to be Central European strains. Nine RABVs (5.7%) detected in foxes, a wildcat and a dog in the Podkarpackie voivodeship were identified as North-Eastern European. A vaccine-induced rabies case was detected in a red fox in the Lubelskie voivodeship in May 2023. Central European and North-Eastern European RABVs were circulating in Poland between 2021 and 2023.

Analysis of the Hybrid Power-Heating System in a Single-Family Building, along with Ecological Aspects of the Operation

This study evaluates a hybrid heating system in a single-family building in northeastern Poland, which has a temperate continental climate. The analysis covers two heating seasons in 2021/2022 and 2022/2023. The hybrid heating system includes an air heat pump HPA–08 CS Plus with a heating power of 8.2 kW (AHP), a condensing gas boiler VC146/5–5 with a power of 14 kW (GB–Condens.), and a solid fuel boiler with a power of 11 kW for central heating. Additionally, hot water is heated by a Basic 270 (DHW’s AHP) air–water heat pump with a power of 2 kW, utilizing a tank with a capacity of 270 dm3 equipped with two heating coils. The building’s average electricity consumption is around 5400 kWh/year. A 4.96 kWp photovoltaic installation is installed on the building’s roof at a 40° angle towards the south to supplement the hybrid system. The study aims to assess whether the PV installation can adequately cover the energy needs of the hybrid heat source for heating and hot water. Furthermore, the study calculates the emission of pollutants (CO2, SOx, NOx, CO, and PM10) into the atmosphere. The total annual electricity production from PV installations was 5444.9 kWh in 2021/2022 and 5684.8 kWh in 2022/2023. The excess electricity was stored in the PGE power grid as per the Prosumer settlement rules. The installed PV installation is sufficient to power the following devices annually: AHP, DHW’s AHP, and GB–Condens. However, the daily electricity production from the PV installation is not enough to cover the energy needs of the heat pump for heating during the cold months in Poland (I–III, XI–XII). It can meet the power needs of a PC all year round and can also be stored during the summer months, for example, in energy warehouses or by directly storing it in the PGE power grid. The use of the PV installation resulted in an average reduction in pollutant emissions into the atmosphere: CO2—94.1%, SOx—91.8%, NOx—95.6%, CO—9.7%, and PM10—32.1%.

Jerusalem Artichoke: Energy Balance in Annual and Perennial Cropping Systems—A Case Study in North-Eastern Poland

This article presents the results of a three-year experiment (2018–2020) conducted at the Agricultural Experiment Station in Bałcyny (north-eastern Poland) with the aim of determining Jerusalem artichoke (JA) yields and the energy balance of biomass production in (i) a perennial cropping system (only aerial biomass was harvested each year) and (ii) an annual cropping system (both aerial biomass and tubers were harvested each year). When JA was grown as a perennial crop, the demand for energy reached 25.2 GJ ha−1 in the year of plantation establishment and 12.3–13.4 GJ ha−1 in the second and third year of production. The energy inputs associated with the annual cropping system were determined in the range of 31.4–37.1 GJ ha−1. Biomass yields were twice as high in the annual than in the perennial cropping system (20.98 vs. 10.30 Mg DM ha−1). Tuber yield accounted for 46% of the total yield. The energy output of JA biomass was 1.8 times higher in the annual than in the perennial cropping system (275.4 vs. 157.3 GJ ha−1). The average energy gain in JA cultivation ranged from 140 (perennial crop) to 241 GJ ha−1 (annual crop). The energy efficiency ratio of JA biomass production reached 7.7–13.3 in the perennial cropping system, and it was 20% lower in the annual cropping system. These results imply that when JA was grown as an annual crop, an increase in energy inputs associated with plantation establishment (tillage and planting) and the harvest and transport of tubers was not fully compensated by the energy output of tubers.

Oilseed Radish: Nitrogen and Sulfur Management Strategies for Seed Yield and Quality—A Case Study in Poland

Oilseed Radish: Nitrogen and Sulfur Management Strategies for Seed Yield and Quality—A Case Study in Poland