Challenges and opportunities of sustainability, certifications and traceability in the Italian beekeeping sector

Giulia Mastromonaco; Simone Blanc; Antonina Sparacino; Chiara Medoro; Stefano Predieri; Marta Cianciabella; Giulia Mastromonaco; Simone Blanc; Antonina Sparacino; Chiara Medoro; Stefano Predieri; Marta Cianciabella

doi:10.3934/agrfood.2025016

AIMS Agriculture and Food

2025, Volume 10, Issue 2: 314-336. doi: 10.3934/agrfood.2025016

Previous Article Next Article

Research article

Challenges and opportunities of sustainability, certifications and traceability in the Italian beekeeping sector

1.
Department of Agricultural, Forest, and Food Sciences, University of Turin, Largo Paolo Braccini 2, 10095, Grugliasco, Turin, Italy
2.
Institute for Bioeconomy (IBE), Italian National Research Council, Via P. Gobetti 101, 40129 Bologna, Italy

Received: 22 January 2025 Revised: 23 March 2025 Accepted: 29 April 2025 Published: 14 May 2025

The impact of beekeeping on agricultural activity is significant in environmental and economic terms, supporting biodiversity through pollination and contributing to rural livelihoods. However, in recent years, this activity has been increasingly threatened by climate change, competition from imported products, and widespread honey fraud, which undermine the sector's sustainability and profitability. This study explores the perceptions of Italian beekeepers and their associations with sustainability, certification, and traceability in the sector. The aim was to identify the challenges and the opportunities for the Italian beekeeping sector in sustainability, traceability, and certification. To this end, two types of analysis were applied to combine the results of the focus groups with three association/consortia and a survey of 360 beekeepers. Through this combination, the research highlights the challenges facing the sector, particularly those related to climate change and biodiversity loss. Associations emphasized the critical role of technical support and legislative measures in promoting sustainable practices, including improved regulation of hive products imports and combating honey fraud. Beekeepers showed strong awareness of their environmental and social role, particularly in pollination and biodiversity conservation, but less engagement with broader systemic issues such as certified production and climate change mitigation. Amateur beekeepers showed greater sensitivity to sustainability issues than professionals. The findings highlight the importance of collaboration between associations, policy makers and stakeholders to address immediate operational needs and promote long-term sustainability in the beekeeping sector.

Keywords:

Citation: Giulia Mastromonaco, Simone Blanc, Antonina Sparacino, Chiara Medoro, Stefano Predieri, Marta Cianciabella. Challenges and opportunities of sustainability, certifications and traceability in the Italian beekeeping sector[J]. AIMS Agriculture and Food, 2025, 10(2): 314-336. doi: 10.3934/agrfood.2025016

Related Papers:

[1]	Muhammad Farman, Muhammad Azeem, M. O. Ahmad . Analysis of COVID-19 epidemic model with sumudu transform. AIMS Public Health, 2022, 9(2): 316-330. doi: 10.3934/publichealth.2022022
[2]	Yosef Mohamed-Azzam Zakout, Fayez Saud Alreshidi, Ruba Mustafa Elsaid, Hussain Gadelkarim Ahmed . The magnitude of COVID-19 related stress, anxiety and depression associated with intense mass media coverage in Saudi Arabia. AIMS Public Health, 2020, 7(3): 664-678. doi: 10.3934/publichealth.2020052
[3]	Srikanth Umakanthan, Anuradha Chauhan, Madan Mohan Gupta, Pradeep Kumar Sahu, Maryann M Bukelo, Vijay Kumar Chattu . COVID-19 pandemic containment in the Caribbean Region: A review of case-management and public health strategies. AIMS Public Health, 2021, 8(4): 665-681. doi: 10.3934/publichealth.2021053
[4]	Jennifer L Lemacks, Laurie S Abbott, Cali Navarro, Stephanie McCoy, Tammy Greer, Sermin Aras, Michael B Madson, Jacqueline Reese-Smith, Chelsey Lawrick, June Gipson, Byron K Buck, Marcus Johnson . Passive recruitment reach of a lifestyle management program to address obesity in the deep south during the COVID-19 pandemic. AIMS Public Health, 2023, 10(1): 116-128. doi: 10.3934/publichealth.2023010
[5]	Toshikazu Kuniya, Hisashi Inaba . Possible effects of mixed prevention strategy for COVID-19 epidemic: massive testing, quarantine and social distancing. AIMS Public Health, 2020, 7(3): 490-503. doi: 10.3934/publichealth.2020040
[6]	Sushant K Singh . COVID-19: A master stroke of Nature. AIMS Public Health, 2020, 7(2): 393-402. doi: 10.3934/publichealth.2020033
[7]	Carlee Bellapigna, Zornitsa Kalibatseva . Psychosocial risk factors associated with social anxiety, depressive and disordered eating symptoms during COVID-19. AIMS Public Health, 2023, 10(1): 18-34. doi: 10.3934/publichealth.2023003
[8]	Boris G Andryukov, Natalya N Besednova . Older adults: panoramic view on the COVID-19 vaccination. AIMS Public Health, 2021, 8(3): 388-415. doi: 10.3934/publichealth.2021030
[9]	Maria Mercedes Ferreira Caceres, Juan Pablo Sosa, Jannel A Lawrence, Cristina Sestacovschi, Atiyah Tidd-Johnson, Muhammad Haseeb UI Rasool, Vinay Kumar Gadamidi, Saleha Ozair, Krunal Pandav, Claudia Cuevas-Lou, Matthew Parrish, Ivan Rodriguez, Javier Perez Fernandez . The impact of misinformation on the COVID-19 pandemic. AIMS Public Health, 2022, 9(2): 262-277. doi: 10.3934/publichealth.2022018
[10]	Areeb Khalid, Muhammad Waqar Younas, Hashim Khan, Muhammad Sarfraz Khan, Abdur Rehman Malik, Adam Umair Ashraf Butt, Basit Ali . Relationship between knowledge on COVID-19 and psychological distress among students living in quarantine: an email survey. AIMS Public Health, 2021, 8(1): 90-99. doi: 10.3934/publichealth.2021007

Abstract

1. Introduction

In December 2019, the world is facing the emergence of a new pandemic, which is called coronavirus disease 2019 (COVID-19). Then, COVID-19 spreads to world widely over the first two months in 2020. There were 492,510 confirmed cases of COVID-19 infection and 22,185 dead cases in world [1], [2]. Therefore, it poses a continuing threat to human health because of its high transmission efficiency and serious infection consequences as well, it transmits by direct contact. Many researchers have tried to study and understand the dynamical behavior of COVID-19 through the transmission dynamics and calculate the basic reproduction number of COVID-19. It has become a key quantity to determine the spread of epidemics and control it. For example, in [3], Li et al. conducted a study of the first 425 confirmed cases in Wuhan, China, showing that the reproduction number of COVID-19 was 2.2, and revealed that person to person transmission occurred between close contacts. Other research [4] shows that the reproduction number of COVID-19 becomes 2.90, which is being increasing. In [5], Riou et al. studied pattern of early human to human transmission of COVID-19 in Wuhan, China. In [6], Hellewell et al. investigated the feasibility of controlling 2019-nCoV outbreaks by isolation of cases and contacts. Chen et al.[7], suggested mathematical model for simulation the phase-based transmissibility of novel coronavirus. Bentout et al. [8] developed an susceptible exposed infectious recovered model to estimation and prediction for COVID-19 in Algeria. Belgaid et al.[9] suggested and analysis of a model for Coronavirus spread. Owolabi et al. [10] proposed and analyzed a nonlinear epidemiological model for SARS CoV-2 virus with quarantine class. Flaxman et al. [11] suggested and estimating the effects of non-pharmaceutical interventions on COVID-19 in Europe. Kennedy et al. [12] suggested a mathematical model involving the effects of intervention strategies on COVID-19 transmission dynamics. Feng et al. [13] studied a COVID-19 model with the effects of media and quarantine in UK. In this present study, we will show effects of the quarantine strategy and media reports on the spread of COVID-19.

We propose a mathematical model for COVID-19 transmission dynamics with the quarantine strategy and media effects. We start the model formulation by denoting the total size of the population by N which is classified further into five classes, the susceptible S(t), the exposed E(t), the infected I(t), the hospital quarantined Q(t) and the recovery R(t) at any time t, So, $N = S + E + I + Q + R$ . The exposed class means low-level virus carrier, which is considered to be non infectious. The quarantined class in which the individual who is in the process in hospital, we suppose that only those who treat it will be in contact with the infected population. Accordingly, the flow of corona virus pandemic along with the above assumptions can be representing in the following block diagram:

Figure 1. Flow diagram of the compartmental model of COVID-19.

DownLoad: Full-Size Img PowerPoint

And the corresponding dynamical model has formulated through the nonlinear differential equations as follows,

$\begin{array}{l} \frac{d S (t)}{d t} = A - (β_{1} - β_{2} \frac{I}{m + I}) S I - d S, \\ \frac{d E (t)}{d t} = (1 - c) (β_{1} - β_{2} \frac{I}{m + I}) S I - k (β_{1} - β_{2} \frac{I}{m + I}) E I - d E, \\ \frac{d I (t)}{d t} = c (β_{1} - β_{2} \frac{I}{m + I}) S I + k (β_{1} - β_{2} \frac{I}{m + I}) E I - (ϵ + γ_{1} + d + µ) I, \\ \frac{d Q (t)}{d t} = ϵ I - (d + γ_{2}) Q, \\ \frac{d R (t)}{d t} = γ_{1} I + γ_{2} Q - d R, \end{array}$ (1)

with initial conditions

$S (0) > 0, E (0) > 0, I (0) > 0, Q (0) > 0, R (0) > 0.$ (2)

In model (1), the birth rate A is taken into susceptible class and natural death rate of population is given by the parameter d. The susceptible will be infected through sufficient direct contacts with infected people in the absence of media alerts by $β_{1}$ , with fraction parameter c, where $c \in [0,1] .$ The term $β_{2} \frac{I}{m + I}$ reduce the transmission as media continuously alert the susceptible and exposed regarding infected cases and possible preventive measures. Usually, we assume that $β_{1} \geq β_{2}$ . As well, we consider the media awareness cannot stop the outbreak of COVID-19 but can aware the population to minimize the transmission risk through half saturation of media constant m. The death due to disease rate $µ$ affecting from infected class only. k represent to a fraction denoting the level of exogenous re-infection. The quarantined rate is given by $ϵ$ . And the mean recovery rates of class I, Q are $γ_{i}, i = 1,2$ , respectively.

It is easy see that the 4^th and 5^th equations are a linear differential equation with respect to variables I(t) and R(t), which are not appear in the other equations of model (1). Hence model (1) can be reduced to the following model:

In this paper, we will discuss the dynamics of model (3) with initial conditions

$S (0) > 0, E (0) > 0, I (0) > 0.$ (4)

This paper is organized as follows. In section 2, we will build the basic properties of model such as (positivity, boundedness of solutions and basic reproduction number). Existence of equilibrium points is presented in section 3. In section, the phenomenon of backward bifurcation is considered. The local and global stability of equilibrium points are studied in sections 4. In section 5, numerical simulation results are given. We conclude this paper with a brief conclusion.

2. Basic properties of model (3)

2.1. Positivity of solutions

On the positivity of solutions for model (3), we have the following result.

Theorem 2.1 Every solution of (3) with initial values (4) is positive as t > 0.

Proof. Let $t_{1} = s u p {t > 0 : S (t) > 0, E (t) > 0, I (t) > 0} > 0$ . It follows (3) that

$\frac{d S (t)}{d t} = A - (β_{1} - β_{2} \frac{I}{m + I}) S I - d S,$ (5)

which can be written as

$\begin{array}{l} \frac{d}{d t} {S (t) e x p [d t + \int_{0}^{t} (β_{1} - β_{2} \frac{I (τ)}{m + I (τ)}) S (τ) I (τ) d τ]} \\ = A e x p [d t + \int_{0}^{t} (β_{1} - β_{2} \frac{I (τ)}{m + I (τ)}) S (τ) I (τ) d τ] . \end{array}$ (6)

thus,

$\begin{array}{l} S (t_{1}) e x p [d t_{1} + \int_{0}^{t} (β_{1} - β_{2} \frac{I (τ)}{m + I (τ)}) S (τ) I (τ) d τ] - S (0) \\ = \int_{0}^{t_{1}} A e x p [d y + \int_{0}^{y} (β_{1} - β_{2} \frac{I (τ)}{m + I (τ)}) S (τ) I (τ) d τ] d y, \end{array}$ (7)

so that

$\begin{array}{l} S (t_{1}) = S (0) e x p [- d t_{1} - \int_{0}^{_{1}} (β_{1} - β_{2} \frac{I (τ)}{m + I (τ)}) S (τ) I (τ) d τ] \\ + e x p [- d t_{1} - \int_{0}^{_{1}} (β_{1} - β_{2} \frac{I (τ)}{m + I (τ)}) S (τ) I (τ) d τ] \\ \times \int_{0}^{t_{1}} A e x p [d y + \int_{0}^{y} (β_{1} - β_{2} \frac{I (τ)}{m + I (τ)}) S (τ) I (τ) d τ] d y > 0. \end{array}$ (8)

Similarly, it can be shown that E(t) > 0 and I(t) > 0 for all time t > 0. Hence all solutions of the model (3) remain positive for all non-negative initial conditions, as required.

2.2. Boundedness

Theorem 2.2 All solutions of model (1) which initiate in $ℝ_{+}^{5}$ are uniformly bounded.

Proof. Define the function $N (t) = S (t) + E (t) + I (t) + Q (t) + R (t)$ and then take the time derivative of N(t) along the solution of model (1) gives $\frac{d N}{d t} \leq A - L N .$ Then, $\frac{d N}{d t} + L N \leq A,$ where $L = min {d, d + µ}$ .

Now, it is easy to verify that the solution of the above linear differential inequalities can be written as

$N (t) \leq \frac{A}{L} + (N_{0} - \frac{A}{L}) e^{- L t},$ (9)

where $N_{0} = (S (0), E (0), I (0), Q (0), R (0))$ . Hence,

$\underset{t \to \infty}{l i m s u p} N (t) \leq \frac{A}{L} .$ (10)

and $N (t) \leq \frac{A}{L}$ for $\forall t > 0.$ Thus all solutions are uniformly bounded and the proof is complete.

2.3. Basic reproduction number

It is easy to see that model (3) always has a disease-free equilibrium $P_{0} (S_{0},0,0)$ , where $S_{0} = \frac{A}{d}$ . We can calculate the reproduction number $ℛ_{0}$ of model (3) by using the next-generation matrix method illustrated by van den Driessche and Watmough in [14].

$ℛ_{0} = \frac{c β_{1} A}{d (ϵ + γ_{1} + d + µ)} .$ (11)

Consequently, from Theorem 2 of [14], we have the following result.

Theorem 2.3 The disease-free equilibrium $P_{0}$ of the model (3) is locally asymptotically when $ℛ_{0} < 1$ and $P_{0}$ is unstable when $ℛ_{0} > 1$ .

The basic reproduction number for COVID-19 infection $ℛ_{0}$ measures the average number of new COVID-19 infections generated by a single infected individual in a completely susceptible population [14], [15]. Theorem 2.3 implies that COVID-19 can be eliminated from the community (when $ℛ_{0} < 1$ ) if the initial sizes of the sub-populations of the model (3) are in the attraction basin of the disease-free equilibrium $P_{0}$ . To ensure that COVID-19 elimination is independent of the initial sizes of the sub-populations, it is necessary to show that the disease-free equilibrium $P_{0}$ is globally asymptotically stable when $ℛ_{0} < 1$ .

3. Existence the COVID-19 equilibria point and backward bifurcation

In this section, we consider the number of equilibrium solutions the model (3). To do so, let $P^{*} (S^{*}, E^{*}, I^{*})$ be any arbitrary equilibrium of the model (3). Setting the right sides of the model (3) to zero gives

$\begin{array}{l} S^{*} = \frac{A}{X^{*} I^{*} + d}, \\ E^{*} = \frac{(1 - c) A X^{*} I^{*}}{(k X I^{*} + d) (X^{*} I^{*} + d)} . \end{array}$ (12)

here

$X^{*} = β_{1} - β_{2} \frac{I^{*}}{m + I^{*}} .$ (13)

Since we assume $β_{1} > β_{2}$ , S* and E* are positive. now, substituting (12) in 3^rd equation of the model (3) and simplifying it, we get

$D_{1} I^{* 4} + D_{2} I^{* 3} + D_{3} I^{* 2} + D_{4} I^{*} + D_{5} = 0,$ (14)

where

$\begin{array}{l} D_{1} = - k (ϵ + γ_{1} + d + µ) (β_{1} - β_{2})^{2}, \\ D_{2} = k A {(β_{1} - β_{2})}^{2} - (ϵ + γ_{1} + d + µ) (β_{1} - β_{2}) [d (k + 1) + 2 k m β_{1}], \\ D_{3} = A (β_{1} + β_{2}) (2 k m β_{1} + c d) - (ϵ + γ_{1} + d + µ) [m d (k + 1) (2 β_{1} + β_{2}) + d^{2} + k β_{1}^{2} m^{2}], \\ D_{4} = k A β_{1}^{2} m^{2} + c d A [β_{1} (m + 1) - β_{2} - m d (ϵ + γ_{1} + d + µ) (2 d + m β_{1}), \\ D_{5} = d^{2} m^{2} (ϵ + γ_{1} + d + µ) (ℛ_{0} - 1) . \end{array}$ (15)

From (15), we can find that $D_{1} < 0$ . And $D_{5} > 0$ when $ℛ_{0} > 1$ , $D_{5} < 0$ when $ℛ_{0} < 1$ . Thus, the number of possible positive real roots the polynomial (12) can have depends on the signs of D₂, D₃ and D₄. Let $f (x) = D_{1} x^{4} + D_{2} x^{3} + D_{3} x^{2} + D_{4} x + D_{5}$ . The various possibilities for the roots of f(x) can be analyzed using the Descartes Rule of Signs. The various possibilities for the roots of f(x) are tabulated in Table 1.

Table 1. Number of possible positive real roots of equation (14).

Cases	D₁	D₂	D₃	D₄	D₅	R₀	Number of sign changes	Number of possible positive real roots
1	−	+	+	+	+	ℛ₀ > 1	1	1
	−	+	+	+	−	ℛ₀ < 1	2	0,2
2	−	+	+	−	+	ℛ₀ > 1	3	1,3
	−	+	+	−	−	ℛ₀ < 1	2	0,2
3	−	+	−	+	+	ℛ₀ > 1	3	1,3
	−	+	−	+	−	ℛ₀ < 1	4	0,2,4
4	−	+	−	−	+	ℛ₀ > 1	3	1,3
	−	+	−	−	−	ℛ₀ < 1	2	0,2
5	−	−	+	+	+	ℛ₀ > 1	1	1
	−	−	+	+	−	ℛ₀ < 1	2	0,2
6	−	−	+	−	+	ℛ₀ > 1	3	1,3
	−	−	+	−	−	ℛ₀ < 1	2	0,2
7	−	−	−	+	+	ℛ₀ > 1	1	1
	−	−	−	+	−	ℛ₀ < 1	2	0,2
8	−	−	−	−	+	ℛ₀ > 1	1	1
	−	−	−	−	−	ℛ₀ < 1	0	0

| Show Table

DownLoad: CSV

Theorem 3.1 The model (3)

(i) has a unique endemic equilibrium if $ℛ_{0} > 1$ and whenever Cases 1, 5, 7 and 8 are satisfied;

(ii) could have more than one endemic equilibrium if $ℛ_{0} > 1$ and Cases 2, 3, 4 and 6 are satisfied;

(iii) could have 2 or more endemic equilibria if $ℛ_{0} < 1$ and Cases 1–7 are satisfied.

From the 4^th and 5^th equations of model (1) we can determent the values of Q* and R* through

$\begin{matrix} Q^{*} = \frac{ϵ I^{*}}{d + γ_{2}}, \\ R^{*} = \frac{γ_{1} I^{*} + γ_{2} Q^{*}}{d} . \end{matrix}$ (16)

The existence of multiple endemic equilibria when $ℛ_{0} < 1$ suggests the possibility of backward bifurcation [16], where the stable disease-free equilibrium co-exists with a stable endemic equilibrium when $ℛ_{0} < 1$ . This is can be obtained using Centre Manifold Theory.

Theorem 3.2 The model (3) exhibits backward bifurcation whenever $m > \frac{{(1 - c)}^{2} A β_{2}}{d^{2}}$ and no backward bifurcation otherwise.

Proof. To prove existence of backward bifurcation in the model (3) the Center Manifold approach as outlined by Castillo-Chavez and Song in [17] is used.

Firstly, for clarity and understanding of the Center Manifold Theory the model (3) variables are transformed as follows $x_{1} = S, x_{2} = E, x_{3} = I$ . Define $X = (x_{1}, x_{2}, x_{3})^{Τ}$ ( $Τ$ denotes transpose), such that the model (3) can be rewritten as $\frac{d X}{d t} = F (X)$ where $F = (f_{1}, f_{2}, f_{3}) .$ Hence,

$\begin{array}{l} \frac{d x_{1} (t)}{d t} = f_{1} = A - (β_{1} - β_{2} \frac{x_{3}}{m + x_{3}}) x_{1} x_{3} - d x_{1}, \\ \frac{d x_{2} (t)}{d t} = f_{2} = (1 - c) (β_{1} - β_{2} \frac{x_{3}}{m + x_{3}}) x_{1} x_{3} - k (β_{1} - β_{2} \frac{x_{3}}{m + x_{3}}) x_{2} x_{3} - d x_{2}, \\ \frac{d x_{3} (t)}{d t} = f_{3} = c (β_{1} - β_{2} \frac{x_{3}}{m + x_{3}}) x_{1} x_{3} + k (β_{1} - β_{2} \frac{x_{3}}{m + x_{3}}) x_{2} x_{3} - (ϵ + γ_{1} + d + µ) x_{3} . \end{array}$ (17)

Now let $β_{1} = β_{1}^{*}$ be the bifurcation parameter. Observe that at $ℛ_{0} = 1$ ,

$β_{1} = β_{1}^{*} = \frac{d (ϵ + γ_{1} + d + µ)}{c A} .$ (18)

With $β_{1} = β_{1}^{*}$ the transformed model equation (17) has a simple eigenvalue with zero real part and all other eigenvalues are negative (that is has a hyperbolic equilibrium point). Thus, Center Manifold Theory can be applied to investigate the local dynamics of the transformed system (17) near $β_{1} = β_{1}^{*}$ . Now the Jacobian matrix of the transformed system evaluated at COVID-19 free equilibrium $P_{0}$ is obtained as

$J (P_{0}) = (\begin{matrix} - d & 0 & β_{1} S_{0} \\ 0 & - d & (1 - c) β_{1} S_{0} \\ 0 & 0 & c β_{1} S_{0} - (ϵ + γ_{1} + d + µ) \end{matrix}) .$ (19)

It is easy to obtain the right eigenvectors of this Jacobian matrix as $V = (v_{1}, v_{2}, v_{3})^{Τ}$ , where ${(v_{1}, v_{2}, v_{3})}^{Τ} = (\frac{β_{1} S_{0}}{d}, \frac{(1 - c) β_{1} S_{0}}{d},1)$ . Similarly, it is possible to obtain the left eigenvectors which are denoted by $W = (w_{1}, w_{2}, w_{3}) = (0,0,1)$ . Now proceeding to obtain the bifurcation coefficients a and b as defined in Theorem 4.1 in [17].

First the non-vanishing partial derivatives of the transformed model (17) evaluated at COVID-19 free equilibrium are obtained as

$\frac{\partial^{2} f_{1}}{\partial x_{1} \partial x_{3}} = \frac{\partial^{2} f_{1}}{\partial x_{3} \partial x_{1}} = \frac{m}{β_{2}} - β_{1}, \frac{\partial^{2} f_{1}}{\partial x_{3}^{2}} = \frac{2 β S_{0}}{m},$ (20)

$\begin{matrix} \frac{\partial^{2} f_{2}}{\partial x_{1} \partial x_{3}} = \frac{\partial^{2} f_{2}}{\partial x_{3} \partial x_{1}} = (1 - c) β_{1}, \frac{\partial^{2} f_{2}}{\partial x_{2} \partial x_{3}} = \frac{\partial^{2} f_{2}}{\partial x_{3} \partial x_{2}} = - k β_{1}, \\ \frac{\partial^{2} f_{2}}{\partial x_{3}^{2}} = - 2 (1 - c) \frac{β_{2} S_{0}}{m}, \frac{\partial^{2} f_{3}}{\partial x_{1} \partial x_{3}} = \frac{\partial^{2} f_{3}}{\partial x_{3} \partial x_{1}} = c β_{1}, \frac{\partial^{2} f_{3}}{\partial x_{3}^{2}} = - \frac{2 c β_{2}}{m}, \end{matrix}$ (21)

so that

$\begin{array}{l} a = \sum_{k, i, j = 1}^{3} w_{k} v_{i} v_{j} \frac{\partial^{2} f_{k}}{\partial x_{i} \partial x_{j}} \\ = 2 w_{3} v_{1} v_{3} \frac{\partial^{2} f_{3}}{\partial x_{1} \partial x_{3}} + w_{3} v_{3}^{2} \frac{\partial^{2} f_{3}}{\partial x_{3}^{2}} \\ = \frac{2 c β_{1}^{2} S_{0}}{d^{2}} (1 - \frac{{(1 - c)}^{2} β_{2} S_{0}}{d m}) . \end{array}$ (22)

The sign of the bifurcation parameter b is associated with the following non-vanishing partial derivatives of F(X), also evaluated at the disease free equilibrium $P_{0}$ :

$\frac{\partial^{2} f_{1}}{\partial x_{3} \partial β_{1}} = - S_{0}, \frac{\partial^{2} f_{2}}{\partial x_{3} \partial β_{1}} = (1 - c) S_{0}, \frac{\partial^{2} f_{3}}{\partial x_{3} \partial β_{1}} = c S_{0} .$ (23)

The bifurcation coefficient b is obtained as

$\begin{array}{l} b = \sum_{k, i = 1}^{3} v_{k} w_{i} \frac{\partial^{2} f_{k}}{\partial x_{i} \partial β_{1}} \\ = v_{1} w_{3} \frac{\partial^{2} f_{1}}{\partial x_{3} \partial β_{1}} + v_{2} w_{3} \frac{\partial^{2} f_{2}}{\partial x_{3} \partial β_{1}} + v_{3} w_{3} \frac{\partial^{2} f_{3}}{\partial x_{3} \partial β_{1}} \\ = c S_{0} (1 + \frac{(2 - c β_{1}) A}{d^{2}}) > 0. \end{array}$ (24)

Obviously, b is always positive. From Theorem 3.2 the system (17) will exhibit backward bifurcation phenomena if the bifurcation coefficient a is positive. The positivity of a in (22) gives the condition for backward bifurcation, which leads to

$m > \frac{{(1 - c)}^{2} A β_{2}}{d^{2}} .$ (25)

4. Stability analysis

In this section, the stability analysis of the all equilibrium points of model (3) studied as shown in the following theorems by used some criterion.

Theorem 4.1 The COVID-19 equilibrium point P* of the model (3) is locally asymptotically if the following conditions are hold

$β_{2} \frac{I^{*} (2 m + I^{*})}{{(m + I^{*})}^{2}} < β_{1},$ (26)

$[(β_{1} - β_{2} \frac{I^{*} (2 m + I^{*})}{{(m + I^{*})}^{2}}) ((1 - c) S^{*} - k E^{*})] + \frac{d (d - 1 + k)}{k} < X I^{*} < \frac{c d}{k (1 - c) - c} .$ (27)

Proof. The Jacobian matrix of model (3) at $P^{*}$ can be written as

$J (P^{*}) = (\begin{matrix} a_{11} & 0 & a_{13} \\ a_{21} & a_{22} & a_{23} \\ a_{31} & a_{32} & 0 \end{matrix}),$ (28)

here

$\begin{matrix} a_{11} = - (X I^{*} + d), a_{13} = \frac{β_{2} S^{*} I^{*} (2 m + I^{*})}{{(m + I^{*})}^{2}} - β_{1} S^{*}, \\ a_{21} = (1 - c) X I^{*}, a_{22} = - (k X I^{*} + d), \\ a_{23} = (1 - c) (β_{1} S^{*} - \frac{β_{2} S^{*} I^{*} (2 m + I^{*})}{{(m + I^{*})}^{2}}) - k (β_{1} E^{*} - \frac{β_{2} E^{*} I^{*} (2 m + I^{*})}{{(m + I^{*})}^{2}}), \\ a_{31} = c X I^{*}, a_{32} = k X I^{*} . \end{matrix}$ (29)

clearly, the characteristics equation of J(P*) is given by

$λ^{3} + B_{1} λ^{2} + B_{2} λ + B_{3} = 0,$ (30)

where

$\begin{matrix} B_{1} = - [a_{11} + a_{22}], \\ B_{2} = a_{11} a_{22} - a_{13} a_{31} - a_{23} a_{32}, \\ B_{3} = - [a_{11} (- a_{23} a_{32}) + a_{13} (a_{21} a_{32} - a_{22} a_{31})] . \end{matrix}$ (31)

furthermore, we have that

$\begin{matrix} Δ = B_{1} B_{2} - B_{3} \\ = - a_{11} a_{22} [a_{11} + a_{22}] + a_{11} a_{13} a_{31} + a_{22} a_{23} a_{32} + a_{13} a_{21} a_{32} . \end{matrix}$ (32)

Now, according to Routh-huewitz criterion P* will be locally asymptotically stable provided that $B_{1} > 0, B_{3} > 0$ and $Δ > 0.$ It is clear that if above conditions (26)–(27) hold.

The purpose of this section is to investigate the global stability by using Lyapunov function for COVID-19 free equilibrium point and COVID-19 equilibrium point respectively. We obtain the result in the following theorems

Theorem 4.2 The disease-free equilibrium $P_{0}$ is globally asymptotically stable provided that the following condition holds:

$\frac{ℛ_{0}}{c} < 1.$ (33)

Proof. Consider the following function

$V_{0} (S, E, I) = (S - S_{0} - S_{0} ln \frac{S}{S_{0}}) + E + I .$ (34)

clearly, $V_{0} : ℝ_{+}^{3} \to ℝ$ is a continuously differentiable function such that $V_{0} (S_{0}, 0, 0) = 0$ and $V_{0} (S, E, I) > 0, \forall (S, E, I) \neq (S_{0}, 0, 0)$ . Further, we have

$\begin{matrix} \frac{d V_{0}}{d t} = (\frac{S - S_{0}}{S}) [A - (β_{1} - β_{2} \frac{I}{m + I}) S I - d S] \\ + [(1 - c) (β_{1} - β_{2} \frac{I}{m + I}) S I - k (β_{1} - β_{2} \frac{I}{m + I}) E I - d E] \\ + [c (β_{1} - β_{2} \frac{I}{m + I}) S I + k (β_{1} - β_{2} \frac{I}{m + I}) E I - (ϵ + γ_{1} + d + µ) I] . \end{matrix}$ (35)

now, by doing some algebraic manipulation and using the condition (33), we get

$\frac{d V_{0}}{d t} \leq - \frac{d}{S} {(S - S_{0})}^{2} - β_{2} S_{0} \frac{I^{2}}{m + I} - d E - [(ϵ + γ_{1} + d + µ) - β_{1} S_{0}] I .$ (36)

Obviously, $d V_{0} / d t = 0$ at $P_{0} = (S,0,0)$ , moreover $d V_{0} / d t < 0$ otherwise. Hence $d V_{0} / d t$ is negative definite and then the solution starting from any initial point satisfy the condition (33), will approaches asymptotically to COVID-19 free equilibrium point. Hence the proof is complete.

Theorem 4.3 P* in case i of Th. (3.1) is globally asymptotically stable if $ℛ_{0} > 1$ .

Proof. At the COVID-19 equilibrium point $P^{*} = (S^{*}, E^{*}, I^{*}), S^{*}, E^{*}$ and I* satisfies the following equations

$\begin{array}{l} A - (β_{1} - β_{2} \frac{I}{m + I}) S I - d S = 0, \\ (1 - c) (β_{1} - β_{2} \frac{I}{m + I}) S I - k (β_{1} - β_{2} \frac{I}{m + I}) E I - d E = 0, \\ c (β_{1} - β_{2} \frac{I}{m + I}) S I + k (β_{1} - β_{2} \frac{I}{m + I}) E I - (ϵ + γ_{1} + d + µ) I = 0 \end{array}$ (37)

By above equations (4.4) and assumptions

$\frac{S}{S^{*}} = x, \frac{E}{E^{*}} = y, \frac{I}{I^{*}} = u$

we obtian

$\begin{array}{l} \dot{x} = x [\frac{A}{S^{*}} (\frac{1}{x} - 1) - β_{1} I^{*} (u - 1) + \frac{β_{2} I^{2}^{*}}{m + I^{*}} (\frac{u^{2} (m + I^{*})}{m + I} - 1)] \\ \dot{y} = y {(1 - c) [\frac{β_{1} S^{*} I^{*}}{E^{*}} (\frac{x u}{y} - 1) - \frac{β_{2} S^{*} I^{2}^{*}}{(m + I^{*}) E^{*}} (\frac{(m + I^{*}) x u^{2}}{(m + I) y} - 1)] - k β_{1} I^{*} (u - 1) + \frac{k β_{2} I^{2}^{*}}{m + I^{*}} (\frac{(m + I^{*}) u^{2}}{m + I} - 1)} \\ \dot{u} = u [c β_{1} S^{*} (x - 1) - \frac{c β_{2} S^{*} I^{*}}{m + I^{*}} (\frac{(m + I^{*}) x u}{m + I} - 1) + k β_{1} E^{*} (y - 1) - \frac{k β_{2} E^{*} I^{*}}{m + I^{*}} (\frac{(m + I^{*}) y u}{m + I} - 1)] \end{array}$ (38)

now, define the Lyapunov function

$V_{1} = S^{*} (x - 1 - l n x) + E^{*} (y - 1 - l n y) + I^{*} (u - 1 - l n u)$ (39)

clearly, by derivative of $V_{1}$ we get

$\frac{d V_{1}}{d t} = S^{*} \frac{x - 1}{x} \dot{x} + E^{*} \frac{y - 1}{y} \dot{y} + I^{*} \frac{u - 1}{u} \dot{u}$

$\begin{array}{l} \frac{d V_{1}}{d t} = (x - 1) [A (\frac{1}{x} - 1) - β_{1} S^{*} I^{*} (u - 1) + \frac{β_{2} S^{*} I^{2}^{*}}{m + I^{*}} (\frac{u^{2} (m + I^{*})}{m + I} - 1)] \\ + (y - 1) {(1 - c) [β_{1} S^{*} I^{*} (\frac{x u}{y} - 1) - \frac{β_{2} S^{*} I^{2}^{*}}{m + I^{*}} (\frac{(m + I^{*}) x u^{2}}{(m + I) y} - 1)] + \frac{k β_{2} E^{*} I^{2}^{*}}{m + I^{*}} (\frac{(m + I^{*}) u^{2}}{m + I} - 1)} \\ + (u - 1) [c β_{1} S^{*} I^{*} (x - 1) - \frac{c β_{2} S^{*} I^{2}^{*}}{m + I^{*}} (\frac{(m + I^{*}) x u}{m + I} - 1) - \frac{k β_{2} E^{*} I^{2}^{*}}{m + I^{*}} (\frac{(m + I^{*}) y u}{m + I} - 1)] \end{array}$

furthermore, by simplifying the resulting terms, we get that

$\begin{array}{l} = A [2 - x - \frac{1}{x}] + β_{1} S^{*} I^{*} [x + u - c (x + u) - (1 - c) (y + \frac{x u}{y})] \\ - \frac{β_{2} S^{*} I^{2}^{*}}{m + I^{*}} [x - y (1 - c) - c u + (u^{2} - (1 - c) \frac{x u^{2}}{y} - c x u) (\frac{m + I^{*}}{m + I})] \\ - \frac{k β_{2} E^{*} I^{2}^{*}}{m + I^{*}} [y - u + (u^{2} - u y) (\frac{m + I^{*}}{m + I})] \end{array}$

Since the arithmetical mean is greater than, or equal to the geometrical mean, then $2 - x - \frac{1}{x} \leq 0$ for x>0 and $2 - x - \frac{1}{x} = 0$ if and only if x=1; $x + u - c (x + u) - (1 - c) (y + \frac{x u}{y}) \leq 0$ for x,y,u>0 and $x + u - c (x + u) - (1 - c) (y + \frac{x u}{y}) = 0$ if and only if x=y=u=1; $y - u + (u^{2} - u y) (\frac{m + I^{*}}{m + I}) \leq 0$ for y,u>0 and $y - u + (u^{2} - u y) (\frac{m + I^{*}}{m + I}) = 0$ if and only if y=u=1. Therefore, ${\dot{V}}_{1} \leq 0$ for x,y,u>0 and ${\dot{V}}_{1} = 0$ if and only if $x = y = u = 1$ , the maximum invariant set of model (3) on the set ${(x, y, u) : {\dot{V}}_{1} = 0}$ is the singleton (1,1,1). Thus, the COVID-19 equilibrium point P* is globally asymptotically stable if $ℛ_{0} > 1$ ,by LaSalle Invariance Principle [18]. Hence, the proof is complete.

5. Numerical simulation

For the parameters values of model (1.1), we can chosen the parameters values from real data available sense Feb. 24 2020 to Apr. 5 2020. The total population of the Iraq for the year 2020 is approximately 40 × 10⁶ [19]. The life expectancy in Iraq is approximatily 71.08 [19]. Clearly, we can obtain that the natural death rate $d = 3.8545 \times 10^{- 5}$ per day. The birth rate is estimated from A/d=N, and assumed that this is to be the bound population in the disease absence. So, A=1541.8 per day and the other parameters of our model shows that in Table 2.

Table 2. Definitions and values of model parameters.

Parameter	Definition	Value	Source
A	Birth rate	1541.8	[19]
β₁	Transmission contact rate between S and I	0.5	Estimated
c	Fraction constant	[0–1]	Estimated
β₂	Awareness rate	0.1	Estimated
m	Half saturation of media constant	70	Estimated
d	Natural death rate	3.854510⁻⁵	[19],[20]
k	Fraction denoting the level of exogenous re-infection	0.05	Estimated
ϵ	Quarantined rate	1/7	[13]
γ₁	Recovery rate from infected wihout quarantin strategy	0.033	Estimated
γ₂	Recovery rate from quarantin class	1/18	[13]
µ	Death due to disease rate	0.38	[19]

| Show Table

DownLoad: CSV

We plot the solution trajectories of model (1) with initial point (15,20,500,1000,150) which converges to COVID-19 equilibrium point P*=(1,27,2773,5428,19371), shown that in Figure 2.

Figure 2. Solution trajectories converge to COVID-19 equilibrium point P*=(1,27,2773,5428,19371), by parameter value in Table 2.

DownLoad: Full-Size Img PowerPoint

Table 3. Different government control measures and corresponding β₁ values.

No.	Date	Government measures	β₁
1	Feb. 24 2020	(1) detection of the first case of COVID-19 in Iraq	0.3
1	Feb. 24 2020	(2) quarantined as preliminary control	0.3
2	Feb. 25 2020	(1) medical examination for all individuals who are in contact with the affected case	0.1
		(2) cancellation of some mass gatherings
		(3) increase the awareness programs about prevention measures
3	Feb. 25-Mar. 24 2020	(1) cancellation of all religious and social events throughout Iraq	0.09
		(2) preventing movement between all provinces
		(3) the suspension of attendance at universities and schools
		(4) providing a number of hospitals to be places for prevention confirmed cases
4	Mar. 24-Apr. 5 2020	(1) close all borders with neighboring countries	0.08
		(2) to declare a state of emergency and impose a curfew
		(3) medical support from the government
		(4) methodological improvement on the diagnosis and treatment strategy
		(5) spontaneous household quarantine by citizens
		(6) more newly-hospitals put into use
		(7) massive online teaching in postponed semester
		(8) addition of new diagnosis method clinically diagnosis in Baghdad and some provinces

| Show Table

DownLoad: CSV

In the face of the COVID-19 outbreak, many stringent measures were taken by Iraqi government will show in Table 3, to simulate the impact of different government control measures on the number of all $S (t), E (t), I (t), Q (t)$ and R(t). We assumed that some values to contact rates with awareness Table 3, with the other parameters in Table 2 staying still on the all stages.

The following Figure 3 shows the values of $S (t), E (t), I (t), Q (t)$ and R(t) under government measures that above to control of COVID-19 outbreak.

Figure 3. Time series to value simulation curve of different values of contact rates β₁=0.3,0.1,0.09,0.08 respectively with keeping other parameters values are taken in Table 2.

DownLoad: Full-Size Img PowerPoint

Clearly, from above figure for effect of contact rate Table 3, We obtain that in case decrease the contact rate (social isolation) the reproduction number less than one and the dynamical behavior of model (1.1) still approaches to COVID-19 equilibrium point. Hence, the backward bifurcation is occur. Now, to investigate the effect of the quarantined strategy it is given by ϵ on the dynamical behavior of model (1.1) and to control to COVID-19 outbreak in Iraq. We study the impact of this parameter on values of $S (t), E (t), I (t), Q (t)$ and R(t) in follows Table 4 and shows the results in Figure 4.

Table 4. Different government control measures and corresponding ϵ values.

No.	Date	Government measures	ϵ
1	Feb. 24 2020	(1) quarantined as preliminary control in Iraq	0.2
2	Feb. 25 2020	(1) medical examination for all individuals who are in contact with the affected case	0.4
		(2) cancellation of some mass gatherings
		(3) increase the awareness programs about prevention measures
3	Feb. 25-Mar. 24 2020	(1) direct the media to explain the symptoms of the epidemic	2.5
		(2) Preventing movement between all provinces
		(3) Providing a number of hospitals to be places for prevention confirmed cases
4	Mar. 24-Apr. 5 2020	(1) to declare a state of emergency and impose a curfew to reduce the contact between people	4.5
		(2) medical support from the government
		(3) methodological improvement on the diagnosis and treatment strategy
		(4) spontaneous household quarantine by citizens
		(5) addition of new diagnosis method clinically diagnosis in Baghdad and some provinces

| Show Table

DownLoad: CSV

The following Figure 4 shows the values of $S (t), E (t), I (t), Q (t)$ and R(t) under government measures that above to control of COVID-19 outbreak.

Figure 4. Time series to value simulation curve of different values of quarantined rates ϵ=0.2,0.4,2.5,4.5 respectively with keeping other parameters values are taken in Table 2.

DownLoad: Full-Size Img PowerPoint

Clearly, from above investigate to impact of the quarantined strategy Table 4, when the quarantine strategy increasing we get the number of infected is decrease and other classes are increase. Here, we ask whether the quarantine strategy is the best solution? The answer is possible, but for specific numbers. Whereas, if the quarantine is more than the capacity of the health institutions. We get the dynamical behavior of model (1.1) lose the stability as shown in Figure 5.

Figure 5. Time series to value simulation curve of different values of quarantined rates 20.5 ≤ ϵ ≤ 30.5. With keeping other parameters values are taken in Table 2.

DownLoad: Full-Size Img PowerPoint

6. Discussion and results

In this research, a mathematical model of COVID-19 transmission has been proposed by compartment the total population into five epidemiological status, namely, susceptible S(t), exposed E(t), infected I(t), quarantine Q(t) and recovered R(t). The model incorporates the impact of social awareness programs conducted by public health officials with quarantine strategy in hospital. It has been noticed that these awareness programs and quarantine strategy result in human behavioral changes in order to avoid risk of disease transmission. The model mainly accounts for the reduction in disease class due to awareness. While we can say the disease goes away due to applied the quarantine it well. The proposed model has two biological equilibrium points are COVID-19 free and COVID-19. The COVID-19 free has been local stability when $ℛ_{0} < 1$ . Otherwise when $ℛ_{0} > 1$ , the COVID-19 free point becomes unstable and the dynamical behavior of the model converges to COVID-19 equilibruim point. The backward bifurcation occur if $ℛ_{0} = 1$ at the parameter bifurcation $β_{1} = β_{1}^{*} = d (ϵ + γ_{1} + d + µ) / c A$ . As well as the different government control measures have been also discussed. Furthermore, to shown and understand the effect of quarantine rate of disease we have choosed many different value of it say parameter then we have obtained some different results see Table 4 and Figure 4.

References

[1]	ISMEA (2024) Api e miele—News e analisi—Tendenze. Available from: https://www.ismeamercati.it/flex/cm/pages/ServeBLOB.php/L/IT/IDPagina/13170.
[2]	Matias DMS, Leventon J, Rau A-L, et al. (2017) A review of ecosystem service benefits from wild bees across social contexts. Ambio 46: 456–467. https://doi.org/10.1007/s13280-016-0844-z doi: 10.1007/s13280-016-0844-z
[3]	Papa G, Maier R, Durazzo A, et al. (2022) The honey bee Apis mellifera: An insect at the interface between human and ecosystem health. Biology 11: 233. https://doi.org/10.3390/biology11020233 doi: 10.3390/biology11020233
[4]	Vinci G, Rapa M, Roscioli F (2018) Sustainable Development in Rural Areas of Mexico through Beekeeping. https://doi.org/10.23958/ijsei/vol04-i08/01
[5]	Chambó E (2016) Beekeeping and Bee Conservation: Advances in Research, BoD—Books on Demand. https://doi.org/10.5772/61424
[6]	Durazzo A, Lucarini M, Plutino M, et al. (2021) Bee products: A representation of biodiversity, sustainability, and health. Life 11: 970. https://doi.org/10.3390/life11090970 doi: 10.3390/life11090970
[7]	Gutiérrez-Briceño I, García-Llorente M, Ortega-Marcos JJ, et al. (2023) Exploring the effect of landscape composition and agroecological practices on wild bees in horticultural farms. Basic Appl Ecol 71: 33–44. https://doi.org/10.1016/j.baae.2023.05.003 doi: 10.1016/j.baae.2023.05.003
[8]	Winfree R (2010) The conservation and restoration of wild bees. Ann New York Acad Sci 1195: 169–197. https://doi.org/10.1111/j.1749-6632.2010.05449.x doi: 10.1111/j.1749-6632.2010.05449.x
[9]	Kim Y, Vonortas NS (2014) Cooperation in the formative years: Evidence from small enterprises in Europe. Eur Manage J 32: 795–805. https://doi.org/10.1016/j.emj.2014.02.003 doi: 10.1016/j.emj.2014.02.003
[10]	Cvitanovic C, Shellock RJ, Mackay M, et al. (2021) Strategies for building and managing 'trust' to enable knowledge exchange at the interface of environmental science and policy. Environ Sci Policy 123: 179–189. https://doi.org/10.1016/j.envsci.2021.05.020 doi: 10.1016/j.envsci.2021.05.020
[11]	Sparacino A, Merlino VM, Blanc S, et al. (2022) A choice experiment model for honey attributes: Italian consumer preferences and socio-demographic profiles. Nutrients 14: 4797. https://doi.org/10.3390/nu14224797 doi: 10.3390/nu14224797
[12]	Zanchini R, Blanc S, Pippinato L, et al. (2022) Consumers' attitude towards honey consumption for its health benefits: First insights from an econometric approach. Br Food J 124: 4372–4386. https://doi.org/10.1108/BFJ-09-2021-0992 doi: 10.1108/BFJ-09-2021-0992
[13]	Šedík P, Hudecová M, Predanócyová K (2023) Exploring consumers' preferences and attitudes to honey: Generation approach in Slovakia. Foods 12: 1941. https://doi.org/10.3390/foods12101941 doi: 10.3390/foods12101941
[14]	Sparacino A, Merlino VM, Borra D, et al. (2023) Web content analysis of beekeeping website companies: Communication and marketing strategies in the Italian context. J Mark Commun 0: 1–22.
[15]	Šedík P, Horská E, Adam Š, et al. (2020) Mineral content as an aspect of nutrition marketing: Case study of honey market in Slovakia. J Food Nutr Res 59: 185–192.
[16]	Pippinato L, Blanc S, Mancuso T, et al. (2020) A sustainable niche market: How does honey behave? Sustainability 12: 10678. https://doi.org/10.3390/su122410678 doi: 10.3390/su122410678
[17]	Covaci B (2020) Natural Honey Market: EU28 and BRICS Competitiveness Evidences Regarding Mountain Natural Honey.
[18]	Simone Blanc, Giulia Mastromonaco, Valentina Maria Merlino, et al. (2024) Indagine socioeconomica sull'apicolturain Piemonte: prospettive e sfide per il futuro del settore, DISAFA Editore.
[19]	Aguado-Gragera C, Sama-Berrocal C, Mesías FJ, et al. (2024) Challenges and opportunities of short food supply chains in Spain: A stakeholder participatory study. J Agric Food Res 18: 101276. https://doi.org/10.1016/j.jafr.2024.101276 doi: 10.1016/j.jafr.2024.101276
[20]	Tarnavölgyi G (2003) Analysis of consumersx attitudes towards food additives using focus group survey. Agric Conspectus Sci 68: 193–196.
[21]	Nyumba TO, Wilson K, Derrick CJ, et al. (2018) The use of focus group discussion methodology: Insights from two decades of application in conservation. Methods Ecol Evolution 9: 20–32. https://doi.org/10.1111/2041-210X.12860 doi: 10.1111/2041-210X.12860
[22]	Bilska B, Tomaszewska M, Kołożyn-Krajewska D (2022) The management of meals in food service establishments in the context of food waste—Results of focus group interviews with employees and owners. Sustainability 14: 9258. https://doi.org/10.3390/su14159258 doi: 10.3390/su14159258
[23]	Baudino C, Giuggioli NR, Briano R, et al. (2017) Integrated methodologies (SWOT, TOWS, LCA) for improving production chains and environmental sustainability of kiwifruit and baby kiwi in Italy. Sustainability 9: 1621. https://doi.org/10.3390/su9091621 doi: 10.3390/su9091621
[24]	King C, Gunton J, Freebairn D, et al. (2000) The sustainability indicator industry: Where to from here? A focus group study to explore the potential of farmer participation in the development of indicators. Aust J Exp Agric 40: 631–642. https://doi.org/10.1071/EA99148 doi: 10.1071/EA99148
[25]	Bruzzese S, Blanc S, Brun F (2020) Strategies for the valorisation of chestnut resources in Italian mountainous areas from a sustainable development perspective. Resources 9: 60. https://doi.org/10.3390/resources9050060 doi: 10.3390/resources9050060
[26]	Blanc S, Brun F, Mosso A (2019) Performance analysis of logging enterprises operating in the Western Italian Alps. Qual-Access Success 20: 74–81.
[27]	Weihrich H (1982) The TOWS matrix—A tool for situational analysis. Long Range Plann 15: 54–66. https://doi.org/10.1016/0024-6301(82)90120-0 doi: 10.1016/0024-6301(82)90120-0
[28]	Clermont A, Eickermann, Michael, Kraus, François, et al. (2014) A survey on some factors potentially affecting losses of managed honey bee colonies in Luxembourg over the winters 2010/2011 and 2011/2012. J Apic Res 53: 43–56. https://doi.org/10.3896/IBRA.1.53.1.04 doi: 10.3896/IBRA.1.53.1.04
[29]	R: The R Project for Statistical Computing. Available from: https://www.r-project.org/.
[30]	Nikolaou IE, Evangelinos KI (2010) A SWOT analysis of environmental management practices in Greek mining and mineral industry. Resour Policy 35: 226–234. https://doi.org/10.1016/j.resourpol.2010.02.002 doi: 10.1016/j.resourpol.2010.02.002
[31]	Van Espen M, Williams JH, Alves F, et al. (2023) Beekeeping in Europe facing climate change: A mixed methods study on perceived impacts and the need to adapt according to stakeholders and beekeepers. Sci Total Environ 888: 164255. https://doi.org/10.1016/j.scitotenv.2023.164255 doi: 10.1016/j.scitotenv.2023.164255
[32]	Vercelli M, Novelli S, Ferrazzi P, et al. (2021) A qualitative analysis of beekeepers' perceptions and farm management adaptations to the impact of climate change on honey bees. Insects 12: 228. https://doi.org/10.3390/insects12030228 doi: 10.3390/insects12030228
[33]	Remotti F, Mattalia G, Porporato M, et al. (2024) 'The rules of nature are changing; every year is unpredictable': perceptions of climate change by beekeepers of Liguria, NW Italy. Reg Environ Change 24: 79. https://doi.org/10.1007/s10113-024-02242-3 doi: 10.1007/s10113-024-02242-3
[34]	Cappellari A, Malagnini V, Fontana P, et al. (2024) Impact of landscape composition on honey bee pollen contamination by pesticides: A multi-residue analysis. Chemosphere 349: 140829. https://doi.org/10.1016/j.chemosphere.2023.140829 doi: 10.1016/j.chemosphere.2023.140829
[35]	Chiesa LM, Labella GF, Giorgi A, et al. (2016) The occurrence of pesticides and persistent organic pollutants in Italian organic honeys from different productive areas in relation to potential environmental pollution. Chemosphere 154: 482–490. https://doi.org/10.1016/j.chemosphere.2016.04.004 doi: 10.1016/j.chemosphere.2016.04.004
[36]	Tosi S, Costa C, Vesco U, et al. (2018) A 3-year survey of Italian honey bee-collected pollen reveals widespread contamination by agricultural pesticides. Sci Total Environ 615: 208–218. https://doi.org/10.1016/j.scitotenv.2017.09.226 doi: 10.1016/j.scitotenv.2017.09.226
[37]	El Agrebi N, Steinhauer N, Renault V, et al. (2022) Beekeepers perception of risks affecting colony loss: A pilot survey. Transboundary Emerging Dis 69: 579–590. https://doi.org/10.1111/tbed.14023 doi: 10.1111/tbed.14023
[38]	Fisher A Ⅱ, Tadei R, Berenbaum M, et al. (2023) Breaking the cycle: Reforming pesticide regulation to protect pollinators. BioScience 73: 808–813. https://doi.org/10.1093/biosci/biad088 doi: 10.1093/biosci/biad088
[39]	Prodanović R, Brkić I, Soleša K, et al. (2024) Beekeeping as a tool for sustainable rural development. J Agron, Technol Eng Manag 7: 1054–1066. https://doi.org/10.55817/IXVM2800 doi: 10.55817/IXVM2800
[40]	Al Naggar Y, Codling G, Giesy JP, et al. (2018) Beekeeping and the need for pollination from an agricultural perspective in Egypt. Bee World 95: 107–112. https://doi.org/10.1080/0005772X.2018.1484202 doi: 10.1080/0005772X.2018.1484202
[41]	Devkota K, Dhakal SC, Thapa RB (2016) Economics of beekeeping as pollination management practices adopted by farmers in Chitwan district of Nepal. Agric Food Secur 5: 6. https://doi.org/10.1186/s40066-016-0053-9 doi: 10.1186/s40066-016-0053-9
[42]	Richardson K (2023) Beekeeping role in enhancing food security and environmental public health. HEM 4: 69–79. https://doi.org/10.61093/hem.2023.4-06 doi: 10.61093/hem.2023.4-06
[43]	Entrena-Durán F, Muñoz-Sánchez V-M, Pérez-Flores A-M (2021) Sustainability and development: From agrarian development to the paradigm of sustainability. Sustainability 13: 6175. https://doi.org/10.3390/su13116175 doi: 10.3390/su13116175
[44]	Huet J-C, Bougueroua L, Kriouile Y, et al. (2022) Digital transformation of beekeeping through the use of a decision making architecture. Appl Sci 12: 11179. https://doi.org/10.3390/app122111179 doi: 10.3390/app122111179
[45]	Moore JC, Spink J, Lipp M (2012) Development and application of a database of food ingredient fraud and economically motivated adulteration from 1980 to 2010. J Food Sci 77: R118–R126. https://doi.org/10.1111/j.1750-3841.2012.02657.x doi: 10.1111/j.1750-3841.2012.02657.x
[46]	García NL (2018) The current situation on the international honey market. Bee World 95: 89–94. https://doi.org/10.1080/0005772X.2018.1483814 doi: 10.1080/0005772X.2018.1483814
[47]	Gustafson CR, Champetier A, Tuyizere O, et al. (2024) The impact of honey fraud information on the valuation of honey origin: Evidence from an incentivized economic experiment. Food Control 155: 110070. https://doi.org/10.1016/j.foodcont.2023.110070 doi: 10.1016/j.foodcont.2023.110070
[48]	Mascarello G, Pinto A, Crovato S, et al. (2024) Consumers' perceptions and behaviors regarding honey purchases and expectations on traceability and sustainability in Italy. Sustainability 16: 8846. https://doi.org/10.3390/su16208846 doi: 10.3390/su16208846
[49]	Schifani G, Romeo P, Dara Guccione G, et al. (2016) Conventions of quality in consumer preference toward local honey in southern Italy. Qual-Access Success 17: 92–97.

This article has been cited by:

1.	Daniel E. Platt, Laxmi Parida, Pierre Zalloua, Lies, Gosh Darn Lies, and not enough good statistics: why epidemic model parameter estimation fails, 2021, 11, 2045-2322, 10.1038/s41598-020-79745-6
2.	Chandan Maji, Mayer Humi, Impact of Media-Induced Fear on the Control of COVID-19 Outbreak: A Mathematical Study, 2021, 2021, 1687-9651, 1, 10.1155/2021/2129490
3.	M. M. Hikal, M. M. A. Elsheikh, W. K. Zahra, Stability analysis of COVID-19 model with fractional-order derivative and a delay in implementing the quarantine strategy, 2021, 1598-5865, 10.1007/s12190-021-01515-y
4.	Ana P. Lemos-Paião, Cristiana J. Silva, Delfim F.M. Torres, A New Compartmental Epidemiological Model for COVID-19 with a Case Study of Portugal, 2020, 44, 1476945X, 100885, 10.1016/j.ecocom.2020.100885
5.	Khalid Hattaf, Mly Ismail El Karimi, Ahmed A. Mohsen, Zakaria Hajhouji, Majda El Younoussi, Noura Yousfi, Mathematical Modeling and Analysis of the Dynamics of RNA Viruses in Presence of Immunity and Treatment: A Case Study of SARS-CoV-2, 2023, 11, 2076-393X, 201, 10.3390/vaccines11020201
6.	Ahmed A. Mohsen, Hassan F. AL-Husseiny, Raid Kamel Naji, The dynamics of Coronavirus pandemic disease model in the existence of a curfew strategy, 2022, 25, 0972-0502, 1777, 10.1080/09720502.2021.2001139
7.	Azhar Iqbal Kashif Butt, Muhammad Imran, Saira Batool, Muneerah AL Nuwairan, Theoretical Analysis of a COVID-19 CF-Fractional Model to Optimally Control the Spread of Pandemic, 2023, 15, 2073-8994, 380, 10.3390/sym15020380
8.	A.I.K. Butt, W. Ahmad, M. Rafiq, D. Baleanu, Numerical analysis of Atangana-Baleanu fractional model to understand the propagation of a novel corona virus pandemic, 2022, 61, 11100168, 7007, 10.1016/j.aej.2021.12.042
9.	Musadaq A. Hadi, Zainab M. Amean, New strategy to control covid-19 pandemic using lead/lag compensator, 2021, 68, 17468094, 102669, 10.1016/j.bspc.2021.102669
10.	Abdelhamid Ajbar, Rubayyi T. Alqahtani, Mourad Boumaza, Dynamics of a COVID-19 Model with a Nonlinear Incidence Rate, Quarantine, Media Effects, and Number of Hospital Beds, 2021, 13, 2073-8994, 947, 10.3390/sym13060947
11.	Srikanth Umakanthan, Anuradha Chauhan, Madan Mohan Gupta, Pradeep Kumar Sahu, Maryann M Bukelo, Vijay Kumar Chattu, COVID-19 pandemic containment in the Caribbean Region: A review of case-management and public health strategies, 2021, 8, 2327-8994, 665, 10.3934/publichealth.2021053
12.	Abdelhamid Ajbar, Rubayyi T. Alqahtani, Mourad Boumaza, Dynamics of an SIR-Based COVID-19 Model With Linear Incidence Rate, Nonlinear Removal Rate, and Public Awareness, 2021, 9, 2296-424X, 10.3389/fphy.2021.634251
13.	T. S. Faniran, E. A. Bakare, A. O. Falade, Ali R. Ashrafi, The COVID-19 Model with Partially Recovered Carriers, 2021, 2021, 1687-0042, 1, 10.1155/2021/6406274
14.	Rahat Zarin, Amir Khan, Ali Akgül, Esra Karatas Akgül, Usa Wannasingha Humphries, Fractional modeling of COVID-19 pandemic model with real data from Pakistan under the ABC operator, 2022, 7, 2473-6988, 15939, 10.3934/math.2022872
15.	Asma Hanif, Azhar Iqbal Kashif Butt, Waheed Ahmad, Numerical approach to solve Caputo‐Fabrizio‐fractional model of corona pandemic with optimal control design and analysis, 2023, 0170-4214, 10.1002/mma.9085
16.	Taye Samuel Faniran, Aatif Ali, Nawal E. Al-Hazmi, Joshua Kiddy K. Asamoah, Taher A. Nofal, Matthew O. Adewole, Dan Selişteanu, New Variant of SARS-CoV-2 Dynamics with Imperfect Vaccine, 2022, 2022, 1099-0526, 1, 10.1155/2022/1062180
17.	Khalid Hattaf, Ahmed A. Mohsen, Jamal Harraq, Naceur Achtaich, Modeling the dynamics of COVID-19 with carrier effect and environmental contamination, 2021, 12, 1793-9623, 2150048, 10.1142/S1793962321500483
18.	Taye Samuel Faniran, Leontine Nkague Nkamba, Thomas Timothee Manga, Qualitative and Quantitative Analyses of COVID-19 Dynamics, 2021, 10, 2075-1680, 210, 10.3390/axioms10030210
19.	Azhar Iqbal Kashif Butt, Doremure Badalge Duresha Chamaleen, Muhammad Imran, Saira Batool, Muneerah AL Nuwairan, A new design and analysis of optimal control problems arising from COVID‐19 outbreak, 2023, 46, 0170-4214, 16957, 10.1002/mma.9482
20.	Andrea Miconi, Simona Pezzano, Elisabetta Risi, Framing pandemic news. Empirical research on Covid-19 representation in the Italian TV news, 2023, 35, 0862397X, 65, 10.58193/ilu.1752
21.	Sonu Lamba, Prashant K. Srivastava, Cost-effective optimal control analysis of a COVID-19 transmission model incorporating community awareness and waning immunity, 2023, 11, 2544-7297, 10.1515/cmb-2023-0154
22.	Vasileios E. Papageorgiou, Georgios Vasiliadis, George Tsaklidis, Analyzing the Asymptotic Behavior of an Extended SEIR Model with Vaccination for COVID-19, 2023, 12, 2227-7390, 55, 10.3390/math12010055
23.	Safoura Rezaei Aderyani, Reza Saadati, Stability and controllability results by n–ary aggregation functions in matrix valued fuzzy n–normed spaces, 2023, 643, 00200255, 119265, 10.1016/j.ins.2023.119265
24.	Laura Francis, Mogens Steffensen, Individual life insurance during epidemics, 2024, 18, 1748-4995, 152, 10.1017/S1748499523000209
25.	A.I.K. Butt, W. Ahmad, M. Rafiq, N. Ahmad, M. Imran, Optimally analyzed fractional Coronavirus model with Atangana–Baleanu derivative, 2023, 53, 22113797, 106929, 10.1016/j.rinp.2023.106929
26.	Brett McKinney, Muhammad Imran, Azhar Iqbal Kashif Butt, SEIR Mathematical Model for Influenza-Corona Co-Infection with Treatment and Hospitalization Compartments and Optimal Control Strategies, 2025, 142, 1526-1506, 1899, 10.32604/cmes.2024.059552

Reader Comments

Your name:*

Email:*
© 2025 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)

通讯作者: 陈斌, bchen63@163.com

1.
沈阳化工大学材料科学与工程学院沈阳 110142

AIMS Agriculture and Food

1.9 3.9

Metrics

Article views(246) PDF downloads(20) Cited by(0)

Preview PDF

Download XML

Export Citation

Article outline

Show full outline

AIMS Agriculture and Food

Challenges and opportunities of sustainability, certifications and traceability in the Italian beekeeping sector

Related Papers:

Abstract

1. Introduction

2. Basic properties of model (3)

2.1. Positivity of solutions

2.2. Boundedness

2.3. Basic reproduction number

3. Existence the COVID-19 equilibria point and backward bifurcation

4. Stability analysis

5. Numerical simulation

6. Discussion and results

References

This article has been cited by:

Reader Comments

通讯作者: 陈斌, bchen63@163.com

Metrics

Figures and Tables

Other Articles By Authors

Catalog

Abstract

1. Introduction

2. Basic properties of model (3)

2.1. Positivity of solutions

2.2. Boundedness

2.3. Basic reproduction number

3. Existence the COVID-19 equilibria point and backward bifurcation

4. Stability analysis

5. Numerical simulation

6. Discussion and results

References

AIMS Agriculture and Food

Challenges and opportunities of sustainability, certifications and traceability in the Italian beekeeping sector

Related Papers:

Abstract

1. Introduction

2. Basic properties of model (3)

2.1. Positivity of solutions

2.2. Boundedness

2.3. Basic reproduction number

3. Existence the COVID-19 equilibria point and backward bifurcation

4. Stability analysis

5. Numerical simulation

6. Discussion and results

References

This article has been cited by:

Reader Comments

通讯作者: 陈斌, bchen63@163.com

Metrics

Figures and Tables

Other Articles By Authors

Related pages

Tools

Export File

Citation

Format

Content

Catalog

Abstract

1. Introduction

2. Basic properties of model (3)

2.1. Positivity of solutions

2.2. Boundedness

2.3. Basic reproduction number

3. Existence the COVID-19 equilibria point and backward bifurcation

4. Stability analysis

5. Numerical simulation

6. Discussion and results

References