
Purpose of the paper: In a society facing epochal challenges, all educational levels are called to propose new answers and strategies. This paper suggests the collaboration between Universities and Primary Schools as a means to spread entrepreneurial mindsets in children and to conceive useful answers to social challenges. Methodology: This research is based on an in-depth analysis of an ongoing project that engages one University and three Primary Schools, involving 175 children aged from 6 to 10 years. To evaluate the project's efficacy, families and teachers were invited to fill out a questionnaire, and a statistical analysis has been done. Results: Findings show the efficacy of the project in terms of the improvement of social orientation, creativity, self-esteem, perseverance, autonomy, and proactivity of children. These skills have been registered within the educational context and the social life of children, underlining a positive effect on society. Research limitations: The main drawbacks refer to two distinct sides: the study and the evaluation of the project's efficacy. Regarding the study, focusing on a single project inhibits the possibility of generalizing results. Further studies can consider enlarging the sample, by analyzing similar projects oriented to primary schools, or propose comparative studies focused on different school levels, to deepen the knowledge of the best age to launch entrepreneurial projects. Concerning the evaluation of the project's efficacy, as it has been done in the short term, results could be affected by the initial enthusiastic appreciation of the proposed activities. However, to assess the long-term efficacy, further studies can replicate the analysis to observe the project's effects on the children's life, jointly considering school and family. Academic and practical implications: Theoretically, this study contributes to the literature on entrepreneurial education in primary schools by introducing the influential role of University. For practitioners, this research offers best practices to inspire other successful cases to spread an early entrepreneurial mindset. Originality of the paper: The study elucidates the relevance of the early entrepreneurial mindset as a means to face the current social challenges and suggests developing that entrepreneurial mindset by building a bridge between University and Primary schools. Specifically, through the in-depth analysis of a project, this study spotlights how Universities can support primary schools in overcoming difficulties due to resources constraint and lack of adequate entrepreneurial and managerial competencies by promoting extracurricular projects managed by researchers in Management and Entrepreneurship.
Citation: Michela Floris, Angela Dettori, Elisabetta Reginato. Budding entrepreneurs. The role of University in spreading early entrepreneurial mindset in school kids[J]. National Accounting Review, 2023, 5(1): 86-107. doi: 10.3934/NAR.2023006
[1] | Renjie Chu, Peiyuan Jin, Hanli Qiao, Quanxi Feng . Intrusion detection in the IoT data streams using concept drift localization. AIMS Mathematics, 2024, 9(1): 1535-1561. doi: 10.3934/math.2024076 |
[2] | Xiangjun Dai, Suli Wang, Baoping Yan, Zhi Mao, Weizhi Xiong . Survival analysis of single-species population diffusion models with chemotaxis in polluted environment. AIMS Mathematics, 2020, 5(6): 6749-6765. doi: 10.3934/math.2020434 |
[3] | Xiangyun Shi, Xiwen Gao, Xueyong Zhou, Yongfeng Li . Analysis of an SQEIAR epidemic model with media coverage and asymptomatic infection. AIMS Mathematics, 2021, 6(11): 12298-12320. doi: 10.3934/math.2021712 |
[4] | Paul Bosch, Héctor J. Carmenate, José M. Rodríguez, José M. Sigarreta . Generalized inequalities involving fractional operators of the Riemann-Liouville type. AIMS Mathematics, 2022, 7(1): 1470-1485. doi: 10.3934/math.2022087 |
[5] | Asma Hanif, Azhar Iqbal Kashif Butt . Atangana-Baleanu fractional dynamics of dengue fever with optimal control strategies. AIMS Mathematics, 2023, 8(7): 15499-15535. doi: 10.3934/math.2023791 |
[6] | Lijuan Chen, Tingting Liu, Fengde Chen . Stability and bifurcation in a two-patch model with additive Allee effect. AIMS Mathematics, 2022, 7(1): 536-551. doi: 10.3934/math.2022034 |
[7] | Miled El Hajji, Mohammed Faraj S. Aloufi, Mohammed H. Alharbi . Influence of seasonality on Zika virus transmission. AIMS Mathematics, 2024, 9(7): 19361-19384. doi: 10.3934/math.2024943 |
[8] | Yousef Alnafisah, Moustafa El-Shahed . Deterministic and stochastic model for the hepatitis C with different types of virus genome. AIMS Mathematics, 2022, 7(7): 11905-11918. doi: 10.3934/math.2022664 |
[9] | Fahad Al Basir, Konstantin B. Blyuss, Ezio Venturino . Stability and bifurcation analysis of a multi-delay model for mosaic disease transmission. AIMS Mathematics, 2023, 8(10): 24545-24567. doi: 10.3934/math.20231252 |
[10] | Jinsheng Guo, Shuang-Ming Wang . Threshold dynamics of a time-periodic two-strain SIRS epidemic model with distributed delay. AIMS Mathematics, 2022, 7(4): 6331-6355. doi: 10.3934/math.2022352 |
Purpose of the paper: In a society facing epochal challenges, all educational levels are called to propose new answers and strategies. This paper suggests the collaboration between Universities and Primary Schools as a means to spread entrepreneurial mindsets in children and to conceive useful answers to social challenges. Methodology: This research is based on an in-depth analysis of an ongoing project that engages one University and three Primary Schools, involving 175 children aged from 6 to 10 years. To evaluate the project's efficacy, families and teachers were invited to fill out a questionnaire, and a statistical analysis has been done. Results: Findings show the efficacy of the project in terms of the improvement of social orientation, creativity, self-esteem, perseverance, autonomy, and proactivity of children. These skills have been registered within the educational context and the social life of children, underlining a positive effect on society. Research limitations: The main drawbacks refer to two distinct sides: the study and the evaluation of the project's efficacy. Regarding the study, focusing on a single project inhibits the possibility of generalizing results. Further studies can consider enlarging the sample, by analyzing similar projects oriented to primary schools, or propose comparative studies focused on different school levels, to deepen the knowledge of the best age to launch entrepreneurial projects. Concerning the evaluation of the project's efficacy, as it has been done in the short term, results could be affected by the initial enthusiastic appreciation of the proposed activities. However, to assess the long-term efficacy, further studies can replicate the analysis to observe the project's effects on the children's life, jointly considering school and family. Academic and practical implications: Theoretically, this study contributes to the literature on entrepreneurial education in primary schools by introducing the influential role of University. For practitioners, this research offers best practices to inspire other successful cases to spread an early entrepreneurial mindset. Originality of the paper: The study elucidates the relevance of the early entrepreneurial mindset as a means to face the current social challenges and suggests developing that entrepreneurial mindset by building a bridge between University and Primary schools. Specifically, through the in-depth analysis of a project, this study spotlights how Universities can support primary schools in overcoming difficulties due to resources constraint and lack of adequate entrepreneurial and managerial competencies by promoting extracurricular projects managed by researchers in Management and Entrepreneurship.
In water resource planning, to quantify the flow requirements in space and time that are necessary to sustain desired ecosystem services [8,40] is referred to as instream flow need (IFN) assessment [3]. Ecologists, mathematicians and environment managers have paid increasing attention to the investigation of spatial population dynamics and invasions of stream or river species, which suggests solutions to IFNs and provides evidence on the influence of flows on ecosystems.
River ecologists are interested in the "drift paradox" problem, which asks how river organisms can persist in a river with unidirectional water flow [34,35]. Population persistence over large spatial and temporal scales is an important component of instream flow need assessment studies [3]. Krkosěk and Lewis [25] proposed three measures for population persistence that relate to lifetime reproductive output of a species in a spatially variable environment and they previously have been applied to river population models [20,31]. These measures are connected through the next-generation operator, which maps the population forward in time from one generation to the next. The first measure, Rloc(x), describes the fundamental niche of the species. It represents the local persistence at any location x in the absence of dispersal, depending only on reproduction and survival at the location x. The second measure, Rδ(x), describes the realized niche. It is the lifetime reproductive output of an individual, initially introduced at x, undergoing survival, reproduction, and dispersal, and hence, its value, larger than one or less than one, also determines the source and sink regions in the habitat. The third measure is the net reproductive rate R0. It represents the average number of offspring produced by a single individual over its lifetime where the spatial distribution of the individual is given by the dominant eigenfunction of the linearized next-generation operator; see theories for R0 in e.g., [9,47]. The net reproductive rate has a long history as measures for population persistence/extinction in ecological modeling. Intuitively, if R0>1 the population will grow, but if R0<1 the population will likely become extinct. While this has the potential to make R0 a powerful tool for studying population persistence, further mathematical work is required if R0 is to be linked rigorously to behavior of the associated nonlinear model. Existing studies of the three measures for river species have revealed the influence of hydrodynamics on long term behaviors of river populations, especially of the net reproductive rate R0, but it is only for idealized spatially homogeneous or one-dimensional (longitudinal) varying rivers (see e.g., [20,31]). Therefore, it is necessary and important to establish persistence theories in more realistic two or three-dimensional rivers.
Most population models in river ecology have been investigated under simple assumptions for hydrodynamics, such as a constant water depth or a constant flow velocity throughout the river (see e.g., [19,29,38,39,41]). A few have taken into account spatial or temporal variations in riverbed structure or the flow regime, but only for simplified and essentially one-dimensional structures (see e.g., [20,23,24,28]). The results from these models are interesting but only provide limited information about the influence of realistic hydrodynamics on spatial population dynamics. There also have been methodologies relying on habitat suitability models and physical habitat availability as a proxy for population status (see e.g., [3,14,21,26,40]). They link flow regimes with ecology qualitatively, but lack a mechanistic foundation, and hence they can well predict habitat quality rather than species long term behavior in the habitat. Thus, direct integration of ecological systems into water flow dynamics and quantitative analysis of the impact of flows on ecology are still rare but in urgent need. A hybrid modeling approach was proposed to simulate the spatial dynamics of macroinvertebrates in a section of the Merced River in central California in [2]. The hydrodynamic model in MIKE 21 Flow Model FM [11] was used to characterize the two dimensional flow field in the longitudinal and lateral directions through the Robinson Reach of this river. The hydraulic prediction from the two dimensional model was then coupled with a particle tracking algorithm [12] to compute the drift dispersal. While the work focused on the effect of hydraulic flow on population dispersal in a two dimensional river, the main work of integration of hydraulic dynamics and ecology was simplified to the study of the steady state of a one dimensional representation model. In [22], a hybrid physical-biological modeling approach was presented to directly link river hydrology with river population models. A governing equation for the gradually varied flow was first coupled to a single population equation in one-dimensional rivers and then in two-dimensional rivers. Moreover, population spread and invasion ratchet phenomenon were studied for these models.
By coupling hydrodynamics and ecological dynamics, one can explicitly investigate the influence of physical factors and hydraulic conditions on spatiotemporal dynamics of a population in a river, which strengthens the ecological component of environmental flow assessments that is currently still lacking. By using the hybrid physical-biological modeling approach coupled to River2D [22], one can observe the evolution of the density of a single stage drifting population in a two-dimensional river. While such calculation can certainly predict the long-term dynamics of persistence or extinction of a population in a specific river, the calculation itself is expensive and the transient running time may be long. Therefore it is advantageous to develop qualitative and quantitative methods that can directly determine population persistence in a two dimensional river.
In this paper, we are interested in species living both in the flowing water and on the benthos of a river (see e.g., [1,13,20,28,38]). We extend the results of persistence measures (Rloc, Rδ and R0) for a benthic-drift population [20] in a one-dimensional river to a benthic-drift model in a two-dimensional depth-averaged river environment, in order to see how different factors, especially hydraulic factors, influence population persistence, and to distinguish how the sedentary stage helps populations persist in a two dimensional river. Our hybrid biological-physical modeling approach, coupled to the River2D program, provides a method to analyze the impact of river morphology on population persistence in a realistic river. It can be used to calculate not only the benthic and drift population densities but also the net reproductive rate of a specific species and its source/sink regions in any depth-averaged river model, given the river morphologic information and population demographic information. To our knowledge, this is the first study that quantitatively describes population persistence in a two-dimensional river environment based on the explicit bio-physical coupling of hydrology and benthic and drift compartments. It is expected that this work could provide an important methodology for ecologists or environment managers to determine population dynamics of a specific species in a specific river.
This paper is organized as follows. In Section 2, we introduce a benthic-drift model for a population in a two-dimensional depth-averaged model of a river. In Section 3, we define the next generating operator for the linearized system and derive three measures of persistence, Rloc(x,y), Rδ(x,y), and R0. We show that R0 can be used as a threshold to determine global population persistence or extinction for the nonlinear model. In Section 4, we present a numerical method to calculate R0 and Rδ when coupling the population model and hydrodynamic model in River2D. In Section 5, we choose two typical rivers and investigate how the water flow, birth rate, transfer rates, diffusion rate, river bottom slope, and river bottom roughness influence population persistence via numerical simulations for R0 and Rδ under different conditions. Finally, a short section of discussion completes the paper.
Denote the two dimensional river region by a bounded and smooth domain Ω in R2 with boundary ∂Ω. Consider the following two dimensional benthic-drift population model, a derivation of which is provided in Appendix A:
{∂Nd∂t=μ(x,y)h(x,y)Nb⏟transfer from Nb−σ(x,y)Nd⏟transfer to Nb−md(x,y)Nd⏟death−1h(x,y)▽⋅(v(x,y)h(x,y)Nd)⏟advection+1h(x,y)▽⋅(D(x,y)h(x,y)▽Nd)⏟diffusion,(x,y)∈Ω,t>0,∂Nb∂t=f(x,y,Nb)Nb⏟reproduction−mb(x,y)Nb⏟death−μ(x,y)Nb⏟transfer to Nd+σ(x,y)h(x,y)Nd⏟transfer from Nd,(x,y)∈Ω,t>0,Nd(x,y,0)=N0d(x,y),Nb(x,y,0)=N0b(x,y),(x,y)∈Ω,a(x,y)Nd+b(x,y)∂Nd∂→n=0,(x,y)∈∂Ω,t>0, | (2.1) |
where Nd(x,y,t) is the population density in the drifting water (unit: 1/volume) at location (x,y) and time t, Nb(x,y,t) is the population density on the benthos (unit: 1/area) at location (x,y) and time t, N0d and N0b are initial population densities, f is the reproduction rate of the population (unit: 1/time), md(x,y) is the mortality rate of the drift population at location (x,y) (unit: 1/time), mb(x,y) is the mortality rate of the benthic population at location (x,y) (unit: 1/time), D(x,y) is the diffusion rate (unit: area/time), μ(x,y) is the per capita rate at which individuals on the benthos enter the drift (unit: 1/time), σ(x,y) is the per capita rate at which individuals return to the benthos from the drift (unit: 1/time), h(x,y) is the river depth at location (x,y) (unit: length), v(x,y)=(v1(x,y),v2(x,y)) is the depth-averaged flow velocity at location (x,y) with v1(x,y) the flow velocity (unit: length/time) in the x direction and v2(x,y) the flow velocity (unit: length/time) in the y direction, ▽=(∂/∂x,∂/∂y), →n is the unit outward normal vector on the boundary, and a(x,y) and b(x,y) are nonnegative bounded functions on ∂Ω satisfying a2+b2≠0. For each t≥0, the solutions Nb(⋅,⋅,t) and Nd(⋅,⋅,t) are in the function space X, which is the continuous function space C(ˉΩ,R) or the subspace of C(ˉΩ,R) consisting of continuously differentiable functions vanishing on the boundary if Dirichlet conditions are applied to some part of the boundary.
The boundary conditions of (2.1) can be Dirichlet, Neumann or Robin conditions. The following typical conditions could be assumed in different portions of the boundary of a river.
(BI) At the river inflow (upstream boundary condition): zero-flux or zero density.
(ⅰ) Zero-flux (individuals cannot enter or leave at the source). The condition is
a(x,y)Nd+b(x,y)∂Nd∂→n=0 with a(x,y)=−v(x,y)⋅→n>0 since v and →n have opposite directions and b(x,y)=D(x,y).
(ⅱ) Zero density (hostile condition) (e.g., the living condition at the upstream end is extremely bad so that individuals die whenever they reach there): Nd=0.
(BB) Along the bank: zero-flux a(x,y)Nd+b(x,y)∂Nd∂→n=0 with a(x,y)=−v(x,y)⋅→n and b(x,y)=D(x,y).
(ⅰ) If no discharge flows into the river through the bank, then v and →n are perpendicular along the bank, hence a=0. Thus, the boundary condition becomes ∂Nd∂→n=0.
(ⅱ) If water flows into the river from the bank of the river (e.g., small streams or tributaries or groundwater flow entering the larger main river), then we have a(x,y)Nd+b(x,y)∂Nd∂→n=0 with a≥0 since v and →n have opposite directions.
(BO) At the river outflow (downstream boundary condition): free-flux or zero density.
(ⅰ) Free-flux (Danckwert's condition): ∂Nd∂→n=0. Individuals can freely leave the river with water flow.
(ⅱ) Zero density (hostile condition) (e.g., all individuals die at the downstream boundary): Nd=0.
We make the following assumptions throughout this paper.
(H1) D,h,v1,v2,∈C2(ˉΩ,(0,∞)).
(H2) μ,σ,md,mb are positive continuous functions in Ω.
(H3) The function f(x,y,Nb) is Lipschitz continuous with respect to Nb with Lipschitz constant c, f(x,y,0)−mb(x,y)<∞, f(x,y,Nb) is monotonically decreasing in Nb and for each (x,y) there exists a unique value K(x,y)>0 such that f(x,y,K(x,y))−mb(x,y)=0, f(x,y,0)−mb(x,y)−μ(x,y)≤0 on Ω.
(H4) μ,σ,md,mb, D, f, a and b do not depend on the vertical height of the location in the river.
In this section, we extend the theories of persistence measures for a one-dimensional benthic-drift model in [20] to obtain three river metrics for population persistence based on the linearized system of (2.1) at the trivial solution (0,0). To this end, we first introduce the next generation operator.
Linearizing system (2.1) at the trivial steady state (N∗d,N∗b)≡(0,0) yields
{∂Nd∂t=μ(x,y)h(x,y)Nb−σ(x,y)Nd−md(x,y)Nd+LNd,(x,y)∈Ω,t>0,∂Nb∂t=f(x,y,0)Nb−mb(x,y)Nb−μ(x,y)Nb+σ(x,y)h(x,y)Nd,(x,y)∈Ω,t>0,Nd(x,y,0)=N0d(x,y),Nb(x,y,0)=N0b(x,y),(x,y)∈ˉΩ,a(x,y)Nd+b(x,y)∂Nd∂→n=0,(x,y)∈∂Ω, | (3.1) |
where L is the linear partial differential operator defined as
LNd:=−1h(x,y)▽⋅(v(x,y)h(x,y)Nd)+1h(x,y)▽⋅(D(x,y)h(x,y)▽Nd), |
for any Nd∈C2(Ω,R).
Remark 1. If f(x,y,0)−mb(x,y)−μ(x,y)>0 on Ω, then
∂Nb∂t=(f(x,y,0)−mb(x,y)−μ(x,y))Nb+σ(x,y)h(x,y)Nd≥(f(x,y,0)−mb(x,y)−μ(x,y))Nb, |
which implies that Nb grows exponentially as t→∞ and (0,0) is unstable for (3.1). Hence, population persists for (2.1) regardless of other conditions. For this reason, we impose the last condition in (H3) in the rest of the paper for interesting results.
Suppose that a population is introduced into the river environment ˉΩ with initial distribution (ψ0d,ψ0b)∈X×X. The individuals in this population then experience dispersal, transfer between mobile and stationary classes and reproduction until they die, so the distribution of these initially introduced individuals at time t, denoted by (ψd(x,y,t),ψb(x,y,t)), is governed by the following system
{∂ψd∂t=μ(x,y)h(x,y)ψb−σ(x,y)ψd−md(x,y,0)ψd+Lψd,(x,y)∈Ω,t>0,∂ψb∂t=−mb(x,y,0)ψb−μ(x,y)ψb+σ(x,y)h(x,y)ψd,(x,y)∈Ω,t>0ψd(x,y,0)=ψ0d,ψb(x,y,0)=ψ0b(x,y),(x,y)∈Ω,a(x,y)ψd+b(x,y)∂ψd∂→n=0,(x,y)∈∂Ω. | (3.2) |
Then f(x,y,0)ψb(x,y,t) is the rate of reproduction by the initially introduced individuals at location (x,y) at time t, and therefore the total reproduction by the initially introduced individuals during their lifetime at (x,y) is given by ∫∞0f(x,y,0)ψb(x,y,t)dt. Noting that offsprings are reproduced only on the benthos, we have the following definition.
Definition 2. The next generation operator Γ:X×X→{0}×X associated with (3.1) is defined by
Γ(ψ0dψ0b)(x,y)=∫∞0(000f(x,y,0))(ψd(x,y,t)ψb(x,y,t))dt=(0f(x,y,0)∫∞0ψb(x,y,t)dt), | (3.3) |
where (ψd(x,y,t),ψb(x,y,t)), the distribution of the initially introduced individuals in the river at time t, is the solution of (3.2).
This next generation operator Γ maps an initial population distribution to its "next generation" (offspring) distribution.
Now we define three river metrics for population persistence as those in [20,31].
River metric 1: Fundamental niche Rloc(x,y) – a local persistence metric that describes the full range of environmental conditions and resources (biological and physical) that the organism can possibly occupy and use, especially when limiting factors are absent.
Assume that an individual only experiences birth and death after being introduced into the river but does not disperse during its lifetime. Define Rloc(x,y) as the number of offspring produced over its lifetime by an individual introduced at location (x,y) in the benthic zone. For (x,y)∈Ω,
Rloc(x,y)=f(x,y,0)∫∞0nb(x,y,t)dt=f(x,y,0)mb(x,y), | (3.4) |
where nb(x,y,t)=e−mb(x,y)tdt is the solution of
{dnbdt=−mb(x,y)nb(x,y,t),x∈Ω,t>0,nb(x,y,0)=1,x∈Ω. | (3.5) |
Then Rloc(x,y)>1 implies that an individual introduced at (x,y) will produce more than one offspring at (x,y) in the next generation, hence the population size at (x,y) will increase over generations. Therefore, locations with Rloc(x,y)>1 correspond to the fundamental niche of the species.
River metric 2: Source-sink distribution Rδ(x,y) – a local persistence metric that determines source and sink regions in the river.
Assume that an individual undergoes birth, death and dispersal after being introduced into the river and we use Rδ(x,y) to describe the number of offspring produced over its lifetime by an individual introduced at location (x,y) in the benthic zone.
Note that the next generation distribution from a single individual introduced at location (x0,y0) is defined by
f(x,y,0)∫∞0ψb(x,y,t)dt, |
where (ψd(x,y,t),ψb(x,y,t)) is the solution of (3.2) with initial conditions ψd(x,y,0)=0, ψb(x,y,0)=δ(x−x0)δ(y−y0) and δ(⋅) is the Dirac delta distribution. Then Rδ(x0,y0) is defined as
Rδ(x0,y0)=∫Ωf(x,y,0)∫∞0ψb(x,y,t)dtdxdy. | (3.6) |
If Rδ(x0,y0)>1, then an individual introduced at location (x0,y0) will produce more than one offspring in the whole river, we call location (x0,y0) a source. If Rδ(x0,y0)<1, then an individual introduced at location (x0,y0) will produce less than one offspring in the whole river, we call location (x0,y0) a sink.
River metric 3: Net reproductive rate R0 – a global persistence metric that determines population persistence or extinction in the whole river.
For any nonnegative initial distribution (N0d(x,y),N0b(x,y)) of the model (3.1), the associated next generation distribution is
Γ(N0dN0b)(x,y)=(0f(x,y,0)∫∞0Nb(x,y,t)dt), | (3.7) |
where (Nd(x,y,t),Nb(x,y,t)) solves (3.2) with initial condition (N0d,N0b).
Define
R0:=r(Γ), | (3.8) |
where r(Γ) is the spectral radius of the linear operator Γ on X×X. Then R0 represents the average number of offspring an individual may produce during its lifetime. We call R0 the net reproductive rate. Similarly as we did in [20], we can show the following result. A sketch of the proof is given in Appendix B.
Theorem 3. Assume (H1)–(H4) for (2.1). The following statements are valid:
(i) If R0<1, then the extinction equilibrium (0,0) is asymptotically stable for (2.1) for all nonnegative initial value conditions.
(ii) If R0>1, then there exists ϵ0>0 such that any nonnegative solution (Nd(x,y,t),Nb(x,y,t)) of (2.1) with (N0d,N0b)≥,≢0 satisfies
lim supt→∞‖(Nd(⋅,⋅,t),Nb(⋅,⋅,t)‖∞=lim supt→∞max(x,y)∈ˉΩ{Nd(x,y,t),Nb(x,y,t)}≥ϵ0. | (3.9) |
Remark 4. Theorem 3 implies that R0 is a threshold for global dynamics of the population in the whole river. Population persists in the river in the sense that the inequality (3.9) is valid if R0>1; population will be extinct in the river if R0<1.
River2D is a hydrodynamic and fish habitat model developed specifically for use in natural streams and rivers. We implemented the population model into River2D to calculate river metrics R0 and Rδ for a species in a river with two dimensional depth-averaged velocities. Starting with creating a preliminary bed topography file from the raw river field data, we edited and refined the data, developed a computational discretization of the river, and then solved for the water depths and velocities throughout the discretization under given inflow and outflow conditions. To obtain steady river flow conditions the River2D model is run from an initial set of conditions with constant inflow and outflow conditions until a steady state of the river is obtained. We then implemented the two dimensional population model equations into the steady river flow in River2D to calculate R0 (see (3.8)) as well the next generation distribution and Rδ (see (3.6)). River2D is a finite element model, based on a conservative Petrov-Galerkin upwinding formulation (see e.g., [7]). The hydrodynamic component of the River2D model is based on the two-dimensional, depth averaged Saint Venant Equations expressed in conservative form. See [15,16,43] for more details about River2D.
In what follows, we present the methods to calculate Rδ and R0 for a benthic-drift model in a spatially two-dimensional (longitudinal-lateral) river in River2D.
For numerical calculation convenience, we introduce a new variable Nv for the offspring that the initially introduced individuals reproduce and consider the following system based on the linearization of model (2.1):
∂Nd(x,y,t)∂t=−σ(x,y)⋅Nd(x,y,t)+μ(x,y)h(x,y)⋅Nb(x,y,t)−md(x,y)⋅Nd(x,y,t)−1h(x,y)∂∂x[v1(x,y)h(x,y)Nd(x,y,t)]+1h(x,y)∂∂x[D(x,y)h(x,y)∂Nd(x,y,t)∂x]−1h(x,y)∂∂y[v2(x,y)h(x,y)Nd(x,y,t)]+1h(x,y)∂∂y[D(x,y)h(x,y)∂Nd(x,y,t)∂y]∂Nb(x,y,t)∂t=σ(x,y)⋅Nd(x,y,t)⋅h(x,y)−μ(x,y)⋅Nb(x,y,t)−mb(x,y)⋅Nb(x,y,t)∂Nv(x,y,t)∂t=f(x,y,0)Nb(x,y,t)Nd(x,y,0)=0,Nb(x,y,0)=N0b,Nv(x,y,0)=0,a(x,y)Nd+b(x,y)∂Nd∂→n=0,(x,y)∈∂Ω,t>0. | (4.1) |
Note that individuals can stay on the benthos or drift in water but the recruitment only occurs when they stay on the benthos. Here Nb and Nd are densities of benthic population and drift population of the first generation (i.e., initially introduced population), respectively; Nv is the population density of the second generation (i.e., offspring of the first generation population). Individuals are only introduced on the benthos. The boundary condition along the river bank is chosen as the no-flux condition and the conditions at the upstream and downstream ends can be chosen as any from boundary assumptions (BI) and (BO) in Section 2. In numerical simulations in the next section, we use the following hostile upstream and Dankwert's downstream conditions,
Nd|upstream=0,∂Nd∂→n|downstream=0, | (4.2) |
or zero flux upstream and Dankwert's downstream conditions,
−(v(x,y)⋅→n)Nd+D(x,y)∂Nd∂→n|upstream=0,∂Nd∂→n|downstream=0. | (4.3) |
For a well-defined river, we first run the river in River2D to the steady flow to obtain v(x,y) and h(x,y). Then we implement the above population model into the steady flow to calculate river metrics. The algorithms for computing Rδ and R0 are given as below.
Note that Nv(x,y,∞) is the total number of new offspring produced at location (x,y) and ∫ΩNv(x,y,∞)dxdy is the total number of offspring produced by initially introduced individuals in their lifetime in the whole river. If a single individual is introduced at location (x0,y0)∈Ω, then
Rδ(x0,y0)=∫ΩNv(x,y,∞)dxdy, | (4.4) |
where Nv(x,y,∞) is the third compartment of the solution to system (4.1) as t→∞, with initial distribution Nb(x0,y0,0)=N0b(x,y)=δ(x−x0,y−y0)=δ(x−x0)⋅δ(y−y0), where δ is the Dirac delta function. Numerically, in River2D, (4.4) is approximated by
Rδ(x0,y0)=∑(xi,yj)∈ˉΩNv(xi,yj,∞), | (4.5) |
where Nv(xi,yj,∞) is approximated by Nv(xi,yj,t1) for sufficiently large t1 and {(xi,yj)i∈I,j∈J} is a division of the river region ˉΩ in the simulation.
Consider system (4.1). For any initially introduced population distribution (0,N0b,0) with ∫ΩN0b(ξ,η)dξdη=1, define
Nl+1v(x,y)=N(l)v(x,y,∞),l=0,1,2,⋯,Nlb(x,y)=Nlv(x,y)∫ΩNlv(ξ,η)dξdη,l=1,2,3,⋯, | (4.6) |
where N(l)v(x,y,∞) is the third compartment of the solution to system (4.1) with the normalized initial condition (0,Nlb,0) as t→∞. Note that the total population size of Nlb is 1 since ∫ΩNlb(ξ,η)dξdη=1. Equations in (4.6) indicate that N(l)v(⋅,⋅,∞) is the final distribution of offspring of the initial population Nlb of total size 1.
In the numerical program in River2D, we approximate R0 by
R0≈∑(xi,yj)∈ˉΩNlv(xi,yj)=∑(xi,yj)∈ˉΩN(l−1)v(xi,yj,∞), | (4.7) |
and the next generation function by
ϕ∗(xi,yj)≈Nlv(xi,yj),(xi,yj)∈ˉΩ, | (4.8) |
for sufficiently large positive integer l, where {(xi,yj)i∈I,j∈J} is a division of the river region ˉΩ in simulation. Intuitively, the above approximations in (4.7) and (4.8) illustrate that R0 represents the average number of offspring a single individual produces in the whole river in its lifetime and that the distribution of all offspring of a single individual is ϕ∗n(x,y). Mathematically, the reason why the net reproductive rate R0 defined in (3.8) can be approximated by (4.7) is given in Appendix C.
Remark 5. The approximations (4.5), (4.7) and (4.8) for Rδ, R0 and the eigenfunction are derived for the general model (4.1); nevertheless, our current numerical program in River2D only assumes constant parameters for (4.1).
We consider a river with an island (see Figure 1) and find persistence metrics of a population in the river. The river is in a region with dimensions 458 m by 247 m (unit: m). The bed elevation and effective bed roughness height are defined as in Figure 1. We consider two flow cases: the instream flow of 200 m3/s and the outflow water depth of 0.5 m; the instream flow of 400 m3/s and the outflow water depth of 0.5 m. Under such flow boundary conditions, the steady flow conditions can be calculated in River2D and the plan view distributions of the water depth and flow velocity are shown in Figure 2.
We then introduce a population into the river and calculate Rδ and R0 as well as the corresponding next generation function under these two different flow conditions. The zero-flux condition is chosen along the bank and hostile conditions (4.2) are chosen at the upstream and downstream ends. We choose the following model parameters: D(x,y)≡1 m2/s, σ(x,y)≡0.75/s, μ(x,y)≡1.8/s, mb(x,y)≡md(x,y)≡0.1/s, and f(x,y,0)≡r=0.5 or 0.322/s.
Two different cases are considered. Case 1: Q=200 m3/s and r=0.5/s; case 2: Q=400 m3/s and r=0.322/s. We obtain that R0=1.54978 in case 1 and R0=0.986818 in case 2. The Rδ distribution and the next generation function are shown in Figure 3, respectively. It follows from Theorem 3 that the population persists in the river in case 1 (when Q=200 m3/s) and will be extinct in case 2 (when Q=400 m3/s). When the flow discharge is low (i.e., in case 1), the Rδ distribution shows that most parts of the river are source regions (Rδ(x,y)>1) that are suitable to introduce the population if it is expected to grow in the river. The next generation function distribution shows that the next generation offspring will not be distributed near the upstream end. They occupy other parts of the river, especially in the middle parts, and this ensures the persistence of the population in the whole river. However, when the flow discharge is large (i.e., in case 2), the population will go extinct even though there are some source regions in the river. Note that we choose different birth rate r(x,y) in the examples here to easily illustrate the phenomenon of persistence and extinction. In fact, even with the same birth rate, the population may persist under low flow but will be extinct under high flow.
We consider a two dimensional periodic meandering river and study the influence of ecological factors and flow characteristics on population persistence metrics in such a river.
The channel of a meandering river is represented by a sine generated curve θ=θmsin(2πs/Ls) (see Figure 4), where θ represents the angle between the channel and the longitudinal line with the maximum θm, s is the distance from the upstream end of the river, Ls is the length of one period of the channel, and the lateral bed slope is tanα where α=αmsin(2πs/Ls) with maximum |αm|. See also, e.g., [27]. We selected the following parameter values to build the meandering river profile for input to the River2D environment: the river length L=800m, the channel period Ls=200 m, the width 20 m, the longitudinal slope of the river bottom S0=0.0002, θm=π/4, αm=−π/1000, the bottom roughness n=0.1m. See Figure 5. In particular, most nodal points inside the numerical domain are uniformly distributed with a 5 m spacing, but nodes near the boundary are slightly shifted and the space discretization on the boundary is smaller than 5m, for better numerical results.
With a specific constant flow discharge at the upstream boundary, River2D can simulate the steady state flow for this river, where the water depths and flow velocities vary periodically. See Figure 6 for the steady flow when Q=0.1m3/s and Q=10m3/s.
1. The co-influence of the flow discharge and the population growth rate. Consider model (2.1) in the meandering river with boundary conditions (4.2). We calculate the net reproductive rate R0 and obtain the next generation function as well as the source/sink distribution Rδ under conditions of different birth rates, death rates, and flow discharges. Figure 7 shows the relationship between R0 and the flow discharge Q under different growth rate conditions. R0 decreases as the flow increases, but increases as the birth rate increases. When the flow only carries organisms to the downstream but does not help them grow, it has only a negative effect on population persistence in the river; a higher growth rate clearly helps population persistence. This coincides with the analysis and simulations for persistence conditions for single population models and benthic-drift population models; see [20,31,38]. Figures 8 and 9 show the next generation function distribution and the distribution of Rδ under different flow conditions. While the next generation function values are almost zero in the first three periods of the river, Figures 8 only shows the next generation function distribution in the last period of the river. Note that the next generation function corresponds to the steady distribution of offspring. This indicates that if the upstream boundary condition is hostile, then in a long run, the next generation concentrates near the downstream end, especially in the region with low flow velocity at the downstream end. Figure 9 shows that the source regions are located in the shallow areas with relatively low flow velocities and are more likely in the upstream. When the flow discharge becomes larger, the downstream contains more sink regions and hence becomes a worse environment for the population to live. Note that here we assume that the flow only has a negative effect on the population. In some rivers, the flow may help persistence by bringing nutrients to organisms; see [19] for drift feeders. In our simulations, the birth rate is assumed to be constant. If the birth rate increases in the downstream direction, then the increased flow plays a trade-off role (see [31] and [20] for models in one-dimensional rivers). Moreover, higher flow rate may be associated with changes in dissolved oxygen [50], hence it would be interesting to investigate how the flow rate affects the concentration of dissolved oxygen in a river, and thus affects the population survival.
2. The effect of the transmission rates. We consider the influence of individuals' transfer rates between the benthos and the water column on population persistence. Figure 10 shows the net reproductive rate and corresponding next generation function under boundary conditions (4.3) for different transfer rates. We choose the ratio between σ (the transfer rate from benthos to water column) and μ (the transfer rate from water column to benthos) to be 1/4. It shows that the smaller the transfer rate is, the larger the net reproductive rate is, which coincides with the result from the benthic-drift model in a one-dimensional river (see Figure 5 in [20]). Therefore, when the ratio between transfer rates is fixed, more frequent transfer between the benthos and the water column does not help population persistence in the whole river. Note that Remark 1 implies that if f(x,y,0)≡r>μ(x,y)+mb(x,y), then population persistence is guaranteed. In our situation, this means that if σ(x,y)<(f(x,y,0)−mb(x,y))/4, then there is unconditional persistence, regardless of the flow rate or other parameter conditions. This is consistent with the result for the benthic-drift model in a one-dimensional river in [38] and again indicates that the lower the transfer is, the better it would be for population persistence.
3. The effect of the diffusion rate. The net reproductive rate and corresponding next generation function with boundary conditions (4.3) for different diffusion rates are shown in Figure 11. The net reproductive rate slightly increases with the diffusion rate here. Hence, under the flow condition given in the example, diffusion in the water column helps population persistence in the river. This is because in a river with low flow more random movement of the population increases the chance of individuals' upstream dispersal and hence helps persistence since the upstream end is not lethal. This is consistent with the results from the one-dimensional models; see e.g., Figure 3.2 in [23].
4. The effect of the river bottom roughness height. The net reproductive rate and next generation function with boundary conditions (4.2) corresponding to different river bottom roughness heights are shown in Figure 12. The net reproductive rate increases with the bottom roughness height n of the river. As the bottom roughness indicates the living condition on the benthos, this implies that a river with a rough bottom (e.g., a river with rocks or stones on the bottom) is better for a population to live and persist than a river with smooth bottom (e.g. sand). It has been shown in [22] (see Figure 3(d) in [22]) that the river bottom roughness increases the upstream speed of a population in a one-dimensional river. Since the conditions for positive upstream speed in an infinitely long river and that for population persistence in a bounded river have been proved to be the same in many literatures (see e.g., [23,24,29]), our result that the river bottom roughness helps population persistence in a two-dimensional river is essentially consistent with that in [22].
5. The effect of the river bottom slope. The net reproductive rate and next generation function with boundary conditions (4.2) corresponding to different river bottom slopes are shown in Figure 13. The net reproductive rate decreases as the river bottom slope increases. Therefore, it is easier for the population to persist in a river with flatter bottom than in a river with steeper bottom. This is also consistent with the result in [22] that the river bottom slope decreases the upstream speed of a population in a one-dimensional river. In a previous work [20] for a one-dimensional benthic-drift model, we observed that when deep pools are close to the upstream, more offspring are produced and live upstream to avoid being washed out and hence deep pools close to upstream help population persistence. In a meandering river constructed as in our examples, deep pools and shallow regions may appear on the same cross section, so the effect of pools on persistence may not be as clear as in the one-dimensional case.
6. The effect of boundary conditions. We use boundary conditions (4.2) (with hostile upstream boundary condition) and (4.3) (with zero-flux upstream boundary conditions) in above simulations. It can been seen from Figures 8 to 15 that boundary conditions greatly influence the value of the net reproductive rate and the next generation function distribution as well as the source/sink regions. When a hostile boundary condition is applied at the upstream end, the upstream part of the river becomes a poor region for the distribution of the next generation function; the offspring occupy only the region near the downstream end and the density near the downstream end becomes high. See Figures 8, 12, and 13. By way of contrast, when the zero-flux boundary condition is applied at the upstream end, the entire river becomes a more suitable region for the next generation distribution, and hence, the next generation offspring are more evenly distributed in the river, although the population density is higher in deep water than in shallow water. See Figures 10, 11, and 14. Moreover, the zero-flux boundary condition at the upstream end yields a much larger net reproductive rate than a hostile boundary condition does. The source/sink regions in the river are also very different under the two different boundary conditions. Under the hostile boundary condition at the upstream, the upstream end is the sink region and the source regions are concentrated in the shallow areas in the river (see Figure 9). However, under the zero-flux condition at the upstream, Rδ attains its maximum at the upstream end and the upstream end becomes the most important source region for the population; see Figure 15. The different scales in the lower part of Figure 15 also illustrate why R0 is higher when Q=10 than when Q=0.1; Rδ is an order of magnitude larger when Q=10 compared to when it is Q=0.1. We also see a strange phenomenon, unlike in Figure 7, about the next generation function distribution for system (4.1) under the zero-flux upstream boundary condition (4.3), when varying the flow discharge; see Figure 14. When the flow discharge increases, the net reproductive rate first increases and then decreases. This seems to be hard to understand. When the water discharge increases, the water depth and flow velocity increase throughout the river. It is hard to see which parameter causes the drop of the net reproductive rate when the discharge increases from 5 to 10.
Understanding hydrodynamics and ecological dynamics is crucial in stream and river management. While there is a need to integrate hydraulic and biological features to discover how river morphology and water flows affect the ecological status of rivers [8,40], a reasonable and efficient method for such integration becomes critical. Based on the fruitful development of habitat models in river habitat assessment, most existing methods (e.g., the physical habitat simulation model (PHABSIM), ecological limits of hydrologic alteration (ELOHA) framework, Software for Assisted Habitat Modeling (SAHM), etc.) link habitat suitability index of river species (e.g., fish) to physical conditions, focusing primarily on habitat suitability and availability; see e.g., [4,32,33,36,40,43,44,46,48].
A few works have attempted to incorporate ecological factors explicitly into hydrodynamic modeling analysis; see e.g., [2,6,18,22,30]. The modelling framework in this paper extends the approach in [22]. The results therein have made it possible to directly analyze how river morphology influences short and long term behaviors of a population in a river. In the current work, by analyses and calculations for the coupled hydrodynamics and population models in partial differential equations, we are able to determine the source/sink regions and global dynamics of a specific population in a two-dimensional depth-averaged river model. Hydraulic, physical and demographic features are explicitly incorporated into the model and their effects on suitable habitat or population persistence/extinction can be numerically analyzed.
Our results showed that the increase of the growth rate and the diffusion rate yields larger net reproductive rate and helps population persistence while larger transfer rate decreases the net reproductive rate and drives population to be extinct. Moreover, higher river bottom roughness height helps population persist but higher bottom slope leads to population extinction more quickly, which might be due to the fact that changes to these parameters change the flow velocity, to a lower level in the former case and to a higher level in the latter case. When the upstream end is imposed with hostile boundary condition, the offspring concentrate in the downstream region and the source regions are in the shallow areas throughout the river. When the upstream end is imposed with the zero-flux boundary condition, the offspring concentrate in deep water regions and the source regions are mainly in the shallow water regions with the most significant source at the upstream end. We also saw a challenging phenomenon which is very hard to understand: When the upstream end is imposed with the zero-flux conditions, when the flow discharge increases, the net reproductive rate may increase first and then decrease.
In order to use the River2D program to calculate Rδ and R0 for a population in a river, the following data need to be collected by ecologists and river managers: Biological data including the birth rate, the death rate, the diffusion rate, transfer rates between the benthos and the drifting water, and population boundary conditions; river bed topography data including bed location (x-coordinate, y-coordinate), bed elevation, bed roughness height, and boundary type; and flow data such as the upstream inflow and the downstream water surface level. It is our hope that this work would be a useful tool for ecologists to predict or control long term behaviors of river and stream species and for water resources managers in identifying more accurately the targets for flow regulation.
In the current River2D program, the biological parameters (birth, death, diffusion, and transfer rates) are assumed to be constants. We will further incorporate spatial heterogeneities of these parameters in the program so that it can be widely used to investigate population persistence and source-sink dynamics in more realistic scenarios. Moreover, since the living conditions for aquatic species in rivers can vary seasonally, the theory developed in this work could be extended to more general situations by incorporating temporal variations in population demography features and flow regime. As a future work, we may consider dynamics of interacting species models coupled with hydrodynamic equations and hence provide better river ecology management strategies.
The authors thank Peter M. Steffler for valuable discussions and comments. The authors also wish to thank Lewis Lab for fruitful discussions. Y. J. gratefully acknowledges an AMS-Simons Travel grant and NSF Grant (DMS 1411703). Q. H. gratefully acknowledges support from NSF of China (11871060), the Fundamental Research Funds for the Central Universities (XDJK2018B031), and the faculty startup fund from Southwest University (20710948). M. A. L. also gratefully acknowledges a Canada Research Chair, NSERC Discovery and Acceleration grants, and a Killam Research Fellowship.
The authors declare no conflict of interest.
We first develop a single-compartment model, and then extend it to a benthic-drift model by dividing the river channel into two zones and the population into two corresponding groups.
Let N(x,y,t) be the population density (unit: quantity/volume) at location (x,y) and time t. Let h(x,y) be the water depth (unit: length) at location (x,y) and v(x,y)=(v1(x,y),v2(x,y)) be the depth averaged flow velocity (unit: length/time) at (x,y) with v1 and v2 being the flow velocities in the x and y directions, respectively.
Consider a water flow box containing the cross-section through (x,y) with the volume estimated by ΔV=ΔxΔyh(x,y) (see Figure 16). We want to understand how the population density in the box changes over time due to flows of population into and out of the box. In the x-direction, the population enters the box with a flux Jx and leaves the box with a flux Jx+Δx. In the y-direction, the population enters the box with a flux Jy and leaves the box with a flux Jy+Δy. The flux expresses the density of population that passes a unit area per unit of time. It has dimension density/(area⋅time). The population density that flows into or out of the box in the x direction per unit time can be estimated by the product of the flux in the x direction and the area of the surface over which the flux occurs, Jx⋅Ax, or Jx+Δx⋅Ax+Δx. Similar formulas apply for the population inflow and outflow in the y direction. We then have the balance equation:
∂N∂t=Jx⋅Ax−Jx+Δx⋅Ax+Δx+Jy⋅Ay−Jy+Δy⋅Ay+ΔyΔV=JxΔyh(x,y)−Jx+ΔxΔyh(x+Δx,y)+JyΔxh(x,y)−Jy+ΔyΔxh(x,y+Δy)ΔxΔyh(x,y)+o(Δx)+o(Δy). | (A.1) |
Letting Δx,Δy→0, we have
∂N∂t=−1h[∂(hJx)∂x+∂(hJy)∂y]. | (A.2) |
Microscopically, the flux due to advection can be written as
Jx|advection=v1(x,y)N,Jy|advection=v2(x,y)N. | (A.3) |
The flux due to dispersion can be written as
Jx|dispersion=−D(x,y)∂N∂x,Jy|dispersion=−D(x,y)∂N∂y, | (A.4) |
where D>0 is the diffusion coefficient. Moreover, Jx=Jx|advection+Jx|dispersion and Jy=Jy|advection+Jy|dispersion.
By combining the flux divergence formula (A.2) with the microscopic formula for advection (A.3) and dispersive flux (A.4) and also considering a reproduction or decay process of the population, we obtain the following general 2-dimensional model for a population that only lives in the drifting water.
∂N∂t=g(x,y,N)N−1h(x,y)[∂∂x(v1(x,y)h(x,y)N(x,y,t))+∂∂y(v2(x,y)h(x,y)N(x,y,t))]+1h(x,y)[∂∂x(D(x,y)h(x,y)∂N(x,y,t)∂x)+∂∂y(D(x,y)h(x,y)∂N(x,y,t)∂y)], | (A.5) |
where g(x,y,N) is the population growth rate.
Hydrologically, the presence of free-flowing water zones on the top and transient storage zones along the bottom in rivers is an important hydraulic characteristics in the ecology of streams. In the storage zones, water movement can be approximated as zero flow [5,10]. Ecologically, many aquatic organisms reside mainly in the storage zone but occasionally jump into the free-flowing zone and drifting downstream until they settle down on the benthos again [1]. Motivated by these facts, we extend the single-compartment model (A.5) to the following benthic-drift model by partitioning the river into two zones, drift zone and benthic zone, and dividing the population into two interacting compartments, individuals residing on the benthos and individuals dispersing in the drift zone:
{∂Nd∂t=μ(x,y)h(x,y)Nb⏟transfer from Nb−σ(x,y)Nd⏟transfer to Nb−md(x,y)Nd⏟death−1h(x,y)▽⋅(v(x,y)h(x,y)Nd)⏟advection+1h(x,y)▽⋅(D(x,y)h(x,y)▽Nd)⏟diffusion,∂Nb∂t=f(x,y,Nb)Nb⏟reproduction−mb(x,y)Nb⏟death−μ(x,y)Nb⏟transfer to Nd+σ(x,y)h(x,y)Nd⏟transfer from Nd, | (A.6) |
where Nd and Nb represent the population density in the drifting water and the population density on the benthos, respectively, f is the reproduction rate of the population md and mb are the mortality rates of individuals in the drift and individuals on the benthos, respectively, and ▽=(∂/∂x,∂/∂y).
Theorem 3 can be proved by following the same process as in Section 3.3 in [20]. Substituting Nd(x,y,t)=eλtϕ1(x,y) and Nb(x,y,t)=eλtϕ2(x,y) into (3.1), we obtain the associated eigenvalue problem
{Lϕ1−(σ(x,y)+md(x,y))ϕ1+μ(x,y)h(x,y)ϕ2=λϕ1,(x,y)∈Ω,(f(x,y,0)−mb(x,y)−μ(x,y))ϕ2+σ(x,y)h(x,y)ϕ1=λϕ2,(x,y)∈Ω,a(x,y)ϕ1+b(x,y)∂ϕ1∂→n=0,(x,y)∈∂Ω. | (B.1) |
By applying similar arguments as in Theorem 3 and Lemma 4 in [20], we can prove the following results.
Lemma 6. (a) The eigenvalue problem (B.1) has a simple principal eigenvalue λ∗ with a positive eigenfunction if f(x,y,0)−mb(x,y)−μ(x,y)<0.
(b) R0−1 and λ∗ have the same sign, where λ∗ is the principal eigenvalue of (B.1).
Then by using similar arguments as in the proof of Theorem 5 in [20], we can obtain the following results.
(ⅰ) If λ∗<0, then the extinction equilibrium (0,0) is asymptotically stable for (2.1) for all nonnegative initial value conditions.
(ⅱ) If λ∗>0, then there exists ϵ0>0 such that any positive solution of (2.1) satisfies
lim supt→∞‖(Nd(⋅,⋅,t),Nb(⋅,⋅,t)−(0,0)‖∞=lim supt→∞max(x,y)∈ˉΩ{Nd(x,y,t),Nb(x,y,t)}≥ϵ0. |
These results together with Lemma B.1 complete the proof of Theorem 3.
Let (Nd(x,y,t),Nb(x,y,t),Nv(x,y,t)) be the solution of system (4.1) with initially introduced population distribution (0,N0b,0). Denote the final distribution of offsprings by N1v, i.e.,
N1v(x,y)=Nv(x,y,∞)=limt→∞Nv(x,y,t). |
Let ˉΓ:X→X be defined as
ˉΓ(N0b)(x,y):=N1v(x,y)=Nv(x,y,∞). | (C.1) |
By using similar derivation as in Appendix A.5 in [20], we can obtain the formula of the operator ˉΓ:
ˉΓ(φ)(x,y)=f(x,y,0)φ(x,y)mb(x,y)+μ(x,y)+σ(x,y)Ad(x,y)f(x,y,0)(mb(x,y)+μ(x,y))Ab(x,y)∫L0∫L0k(x,y,ξ,η)μ(ξ,η)Ab(ξ,η)φ(ξ,η)(mb(ξ,η)+μ(ξ,η))Ad(ξ,η)dξdη | (C.2) |
for any φ∈X, where k(x,y,ξ,η) is the solution of the ordinary boundary value problem
{(L−σ(x,y)−md(x,y)+σ(x,y)μ(x,y)mb(x,y)+μ(x,y))k(x,y,ξ,η)=−δ(x−ξ)δ(y−η),(x,y)∈Ωa(x,y)Nd+b(x,y)∂Nd∂→n=0,(x,y)∈∂Ω,t>0, | (C.3) |
for a fixed point (ξ,η)∈Ω. Note that the solution to (C.3) is a Green's function (see (2–11) in Chapter 7 in [17], (2.9) in Chapter 3 in [42], or Appendix B in [31]). Furthermore, we obtain the following result.
Proposition 7. The spectral radius of Γ is equal to the spectral radius of ˉΓ, i.e.,
R0=r(Γ)=r(ˉΓ). | (C.4) |
Proof. Define an operator ˆΓ:X→X as
ˆΓ(φ)(x,y)=f(x,y,0)φ(x,y)mb(x,y)+μ(x,y)+σ(x,y)Ad(x,y)(mb(x,y)+μ(x,y))Ab(x,y)∫L0∫L0k(x,y,ξ,η)f(ξ,η,0)μ(ξ,η)Ab(ξ,η)φ(ξ,η)(mb(ξ,η)+μ(ξ,η))Ad(ξ,η)dξdη, | (C.5) |
for all φ∈X, where k(x,y,ξ,η) is the solution of the ordinary boundary value problem (C.3). By similar proof as for Theorem 6 in [20], we obtain that the spectral radius of Γ is equal to the spectral radius of ˆΓ, i.e.,
r(Γ)=r(ˆΓ). |
Define two operators T1,T2:X→X:
T1(φ)(x,y)=f(x,y,0)φ(x,y),T2(φ)(x,y)=φ(x,y)mb(x,y)+μ(x,y)+σ(x,y)Ad(x,y)(mb(x,y)+μ(x,y))Ab(x,y)∫L0∫L0k(x,y,ξ,η)μ(ξ,η)Ab(ξ,η)φ(ξ,η)(mb(ξ,η)+μ(ξ,η))Ad(ξ,η)dξdη, | (C.6) |
for all φ∈X, where k is defined in (C.3). Then both T1 and T2 are bounded and linear, and we further have ˆΓ=T2∘T1, ˉΓ=T1∘T2. This yields the result that the spectral radii of these two operators are the same, i.e.,
r(ˆΓ)=limn→∞||(T2∘T1)n||1/n=limn→∞||(T1∘T2)n||1/n=r(ˉΓ). |
Hence, we obtain
R0=r(Γ)=r(ˆΓ)=r(ˉΓ). |
To numerically approximate ˉΓ, we can discretize (C.2) in the river region ˉΩ. Divide ˉΩ into a grid (uniform except in the area near the boundary). Let {(xi,yj)i∈I,j∈J} represent all the nodes of the grid. Let n be the total number of the nodes and rearrange all nodes such that they are represented by a sequence of points {Ph,1≤h≤n}. By applying a numerical quadrature (e.g., Newton-Cotes, Gauss, Simpson, etc.) to (C.2), we can approximate (C.2) by
ˉΓ(φ)(Ph)≈n∑s=1ηh,sφ(Ps),1≤h≤n, |
where ηh,s>0 depends on parameter functions in (4.1) and the numerical quadrature (see e.g., [37,45]). Let
ˉΓn=(ηh,s). |
Then we can approximate the operator ˉΓ by the positive matrix operator ˉΓn and hence approximate r(ˉΓ) by r(ˉΓn), i.e.,
r(ˉΓ)≈r(ˉΓn), | (C.7) |
for large positive integer n.
The Perron-Frobenius Theorem implies that ˉΓn admits a positive principal eigenvalue λ∗n which is equal to r(ˉΓn) and is associated with a positive eigenvector ϕ∗n. Since ˉΓn is a positive matrix, we can use the power method to numerically calculate λ∗n. For any nonnegative function N0b with ∑nh=1N0b(Ph)=1, let
Nlb(Ph)=ˉΓnNl−1b(Ph)∑nh=1ˉΓnNl−1b(Ph),1≤h≤n,l=1,2,3,⋯. | (C.8) |
We then have
(Nlb)T(ˉΓnNlb)(Nlb)TNlb→λ∗n and Nlb→ϕ∗n, as l→∞. | (C.9) |
Note that we use ˉΓn to approximate ˉΓ. The definition of ˉΓ in (C.1) implies that ˉΓnNl−1b approximates the offspring distribution of system (4.1) with the initial population Nl−1b. Denote
Nlv(Ph)=ˉΓnNl−1b(Ph),1≤h≤n,l=1,2,⋯. |
(C.9) implies that
n∑h=1Nlv(Ph)=n∑h=1ˉΓnNl−1b(Ph)→λ∗n,asl→∞ | (C.10) |
and
Nlv=ˉΓnNl−1b→λ∗nϕ∗n, as l→∞. | (C.11) |
Following this idea, we do not directly find ˉΓn to estimate λ∗n (and hence R0) in our River2D simulations. Instead, we solve (4.1) numerically for any nonnegative function N0b satisfying ∑(xi,yj)∈ˉΩN0b(xi,yj)=1, and define the following sequence by following the idea as in (C.1) and (C.8):
Nl+1v(xi,yj)=N(l)v(xi,yj,∞),l=0,1,2,⋯,Nlb(xi,yj)=Nlv(xi,yj)∑(xi,yj)∈ˉΩNlv(xi,yj),l=1,2,3,⋯, | (C.12) |
where N(l)v(xi,yj,∞) is the offspring compartment of the numerical solution to system (4.1) with the normalized initial condition (0,Nlb,0) as t→∞. Then by (C.4), (C.7), (C.10) and (C.11), we approximate R0 by
R0=r(Γ)=r(ˉΓ)≈r(ˉΓn)=λ∗n≈∑(xi,yj)∈ˉΩNlv(xi,yj), | (C.13) |
and the next generation function by
ϕ∗(xi,yj)≈ϕ∗n(xi,yj)≈Nlv(xi,yj), | (C.14) |
for sufficiently large positive integer l.
[1] |
Acosta Castellanos PM, Queiruga-Dios A (2022) From environmental education to education for sustainable development in higher education: a systematic review. Int J Sust Higher Ed 23: 622–644. https://doi.org/10.1108/IJSHE-04-2021-0167 doi: 10.1108/IJSHE-04-2021-0167
![]() |
[2] |
Adey P (2005) Issues arising from the long-term evaluation of cognitive acceleration programs. Res Sci Educ 35: 3–22. https://doi.org/10.1007/s11165-004-3430-5 doi: 10.1007/s11165-004-3430-5
![]() |
[3] | Adey P (2006) Thinking in Science-Thinking in General, In: Asia-Pacific Forum on Science Learning and Teaching 7: 1–6. |
[4] | Allen IE, Seaman CA (2007) Likert scales and data analyses. Qual Prog 40: 64–65. |
[5] |
Alvesson M, Sandberg J (2011) Generating research questions through problematization. Acad Manage Rev 36: 247–271. https://doi.org/10.5465/amr.2009.0188 doi: 10.5465/amr.2009.0188
![]() |
[6] |
Amirkhanyan AA, Holt SB, Mccrea AM, et al. (2019) Managing Racial Diversity: Matching Internal Strategies with Environmental Needs. Public Admin Rev 79: 69–81. https://doi.org/10.1111/puar.12977 doi: 10.1111/puar.12977
![]() |
[7] |
Armstrong JS, Overton TS (1977) Estimating nonresponse bias in mail surveys. J Mark Res 14: 396–402. https://doi.org/10.2307/3150783 doi: 10.2307/3150783
![]() |
[8] | Audretsch DB, Grilo I, Thurik AR (2007a) Explaining entrepreneurship and the role of policy: a framework. The handbook of research on entrepreneurship policy, 1–17. |
[9] | Audretsch DB, Grilo I, Thurik AR (2007b) Handbook of research on entrepreneurship policy, Edward Elgar Publishing. |
[10] | Avenier MJ (2000) Ingénierie des pratiques collectives: la cordée et le quatuor, Editions L'Harmattan. |
[11] | Bacigalupo M, Kampylis P, Punie Y, et al. (2016) EntreComp: The entrepreneurship competence framework, Publication Office of the European Union, Luxembourg. |
[12] |
Bell R, Liu P (2019) Educator Challenges in the Development and Delivery of Constructivist Active and Experiential Entrepreneurship Classrooms in Chinese Vocational Higher Education. J Small Bus Enterp Dev 26: 209–227. https://doi.org/10.1108/JSBED-01-2018-0025 doi: 10.1108/JSBED-01-2018-0025
![]() |
[13] |
Bentler PM, Chou CP (1987) Practical issues in structural modeling. Sociol Methods Res 16: 78–117. https://doi.org/10.1177/0049124187016001 doi: 10.1177/0049124187016001
![]() |
[14] |
Berkovich I, Eyal O (2019) Transformational Leadership, Transactional Leadership, and Moral Reasoning. Leadersh Policy Sch 20: 131–148. https://doi.org/10.1080/15700763.2019.1585551 doi: 10.1080/15700763.2019.1585551
![]() |
[15] | Beyes T, Parker M, Steyaert C (2016) Introduction: why does management education need reinventing? In: Steyaert, C., Beyes, T., Parker, M., The Routledge Companion to Reinventing Management Education, Routledge, 17–36. https://doi.org/10.4324/9781315852430 |
[16] | Black TR (1999) Doing quantitative research in the social sciences: An integrated approach to research design, measurement, and statistics, Sage. |
[17] | Boffo S, Cocorullo A (2019) University Fourth Mission, Spin-Offs and Academic Entrepreneurship: Connecting Public Policies with New Missions and Management Issues of Universities. High Educ Forum 16: 125–142. |
[18] | Bourgeois A (2011) Entrepreneurship Education at School in Europe: National Strategies, Curricula and Learning Outcomes. |
[19] | Bourgeois A, Balcon MP (2016) Entrepreneurship Education at School in Europe. |
[20] | Bruner JS (1966) Toward a theory of instruction, Harvard University Press. |
[21] | Burns P (2012) Corporate entrepreneurship: innovation and strategy in large organizations. |
[22] |
Cai Y, Lattu A (2022) Triple helix or quadruple helix: which model of innovation to choose for empirical studies? Minerva 60: 257–280. https://doi.org/10.1007/s11024-021-09453-6 doi: 10.1007/s11024-021-09453-6
![]() |
[23] |
Campbell DF, Carayannis EG (2016) The academic firm: a new design and redesign proposition for entrepreneurship in innovation-driven knowledge economy. J Innov Entrep 5: 12. https://doi.org/10.1186/s13731-016-0040-1 doi: 10.1186/s13731-016-0040-1
![]() |
[24] |
Canrinus ET, Klette K, Hammerness K (2019) Diversity in coherence: Strengths and opportunities of three programs. J Teach Educ 70: 192–205. https://doi.org/10.1177/0022487117737305 doi: 10.1177/0022487117737305
![]() |
[25] |
Carayannis EG, Grigoroudis E, Campbell DF, et al. (2018) The ecosystem as helix: an exploratory theory-building study of regional co-opetitive entrepreneurial ecosystems as Quadruple/Quintuple Helix Innovation Models. R D Manage 48: 148–162. https://doi.org/10.1111/radm.12300 doi: 10.1111/radm.12300
![]() |
[26] |
Cascavilla I, Hahn D, Minola T (2022) How You Teach Matters! An Exploratory Study on the Relationship between Teaching Models and Learning Outcomes in Entrepreneurship Education. Adm Sci 12: 12. https://doi.org/10.3390/admsci12010012 doi: 10.3390/admsci12010012
![]() |
[27] |
Chilenga N, Dhliwayo S, Chebo AK (2022) The entrepreneurial mindset and self-employment intention of high school learners: The moderating role of family business ownership. Front Educ 7: 1–13. https://doi.org/10.3389/feduc.2022.946389 doi: 10.3389/feduc.2022.946389
![]() |
[28] |
Churchill Jr GA (1979) A paradigm for developing better measures of marketing constructs. J Mark Res 16: 64–73. https://doi.org/10.2307/3150876 doi: 10.2307/3150876
![]() |
[29] | Collins A (1996) Design issues for learning environments, In: Vosniadou, S., De Corte, E., Glaser, R., et al., International perspectives on the design of technology-supported learning environments, Lawrence Erlbaum Associates, Inc, 347–361. |
[30] | Commission E (2006) The Oslo agenda for entrepreneurship education in Europe, In: Entrepreneurship education in Europe: Fostering entrepreneurial mindsets through education and learning conference. |
[31] | Commission E (2013) Thematic Report PLA on Stakeholder Engagement in Entrepreneurship Education, Retrieved from Copenhagen. |
[32] |
Cooper D (2009) University-Civil Society (U-CS) research relationships: The importance of a 'fourth helix'alongside the 'triple helix'of University-Industry-Government (UIG) relations. S Afr Review Sociol 40: 153–180. https://doi.org/10.1080/21528586.2009.10425106 doi: 10.1080/21528586.2009.10425106
![]() |
[33] |
Cooper D (2015) Social Justice and S outh A frican University Student Enrolment Data by 'Race', 1998–2012: From 'Skewed Revolution'to 'Stalled Revolution. High Educ Quart 69: 237–262. https://doi.org/10.1111/hequ.12074 doi: 10.1111/hequ.12074
![]() |
[34] |
Cooper D (2017) Concepts of "applied and public sociology": Arguments for a bigger theoretical picture around the idea of a "university third mission. J Appl Soc Sci 11: 141–158. https://doi.org/10.1177/1936724417722580 doi: 10.1177/1936724417722580
![]() |
[35] |
Cuesta-Claros A, Malekpour S, Raven R, et al. (2022) Understanding the roles of universities for sustainable development transformations: A framing analysis of university models. Sustain Dev 30: 525–538. https://doi.org/10.1002/sd.2247 doi: 10.1002/sd.2247
![]() |
[36] | Cunha F, Heckman JJ (2007) The evolution of inequality, heterogeneity and uncertainty in labor earnings in the US economy. National Bureau of Economic Research. https://doi.org/10.3386/w13526 |
[37] | Das SC (2014) Learning management from Aesop's Fables and panchatantra or Jatak tales. Bus Manag Rev 5: 434. |
[38] |
Daspit JJ, Fox CJ, Findley SK (2023) Entrepreneurial mindset: An integrated definition, a review of current insights, and directions for future research. J Small Bus Manag 61: 12–44. https://doi.org/10.1080/00472778.2021.1907583 doi: 10.1080/00472778.2021.1907583
![]() |
[39] | Drucker PF (1963) Twelve fables of research management. Harvard Business Review. |
[40] | Dubey P (2022) Entrepreneurial characteristics, attitude and self-employment intention: A study of public and private technical institutes' undergraduates. J Small Bus Manag 24: 21–30. |
[41] |
Eisenhardt KM, Graebner ME (2007) Theory building from cases: Opportunities and challenges. Acad Manage J 50: 25–32. https://doi.org/10.5465/amj.2007.24160888 doi: 10.5465/amj.2007.24160888
![]() |
[42] |
Eldridge K, Larry L, Baird J, et al. (2018) A collaborative governance approach to improving tertiary education in Papua New Guinea. Asia Pac J Educ 38: 78–90. https://doi.org/10.1080/02188791.2018.1423949 doi: 10.1080/02188791.2018.1423949
![]() |
[43] |
Elmore R (1996) Getting to scale with good educational practice. Harvard Educ Rev 66: 1–27. https://doi.org/10.17763/haer.66.1.g73266758j348t33 doi: 10.17763/haer.66.1.g73266758j348t33
![]() |
[44] | Falcone D, Silvestri A, Cerbaso C, et al. (2015) Proposal of a Methodology For Non-Formal Competences Certification. Online J New Horiz Educ 122. |
[45] | Falk-Lundqvist Å, Hallberg PG, Leffler E, et al. (2011) Entrepreneurial pedagogy in schools: driving forces for pupils' learning, Liber, Stockholm. |
[46] | Farrington CA, Roderick M, Allensworth E, et al. (2012) Teaching Adolescents to Become Learners: The Role of Non-cognitive Factors in Shaping School Performance—A Critical Literature Review. |
[47] | Filion LJ (1994) Ten steps to entrepreneurial teaching. J Small Bus Entrep 11: 68–78. |
[48] | Fisher SL, Graham ME, Compeau M (2008) Starting from scratch: Understanding the learning outcomes of undergraduate entrepreneurship education, In: Harrison, R., Leitch, C.M., Entrepreneurial Learning, 335–362, Routledge. https://doi.org/10.4324/9780203931929 |
[49] | Floris M, Dettori A (2020) Fostering Early Entrepreneurial Competencies: An Action Research Approach, In: Sociological Perspectives on Educating Children in Contemporary Society, IGI Global, 259–281. https://doi.org/10.4018/978-1-7998-1847-2.ch011 |
[50] | Floris M, Floris S (2021) Soft skills and entrepreneurial competencies to promote equity in schools. In: Bending the arc towards justice: equity-focused practices for educational leaders, Information Age Publishing, 257–282. |
[51] |
Floris M, Pillitu D (2019) Improving entrepreneurship education in primary schools: a pioneer project. Int J Educ Manag 33: 1148–1169. https://doi.org/10.1108/IJEM-09-2018-0283 doi: 10.1108/IJEM-09-2018-0283
![]() |
[52] | Gasse Y (1985) A strategy for the promotion and identification of potential entrepreneurs at the secondary school level. Front Entrep, 538–559. |
[53] |
Gibb A (2002) In Pursuit of a New 'Enterprise' and 'Entrepreneurship' Paradigm for Learning: Creative Destruction, New Values, New Ways of Doing Things and New Combinations of Knowledge. Int J Manag Rev 4: 213–231. https://doi.org/10.1111/1468-2370.00086 doi: 10.1111/1468-2370.00086
![]() |
[54] | Gibb A (2008) Entrepreneurship and enterprise education in schools and colleges: Insights from UK practice. Int J Entrep Educ 6: 48. |
[55] |
Gilder G (1981) Moral sources of capitalism. Society 18: 24–27. https://doi.org/10.1007/BF03180594 doi: 10.1007/BF03180594
![]() |
[56] | Hair JF, Anderson RE, Tatham RL, et al. (1998) Multivariate data analysis 5: 87–135. |
[57] |
Handscombe RD, Rodriguez-Falcon E, Patterson EA (2008) Embedding Enterprise in Science and Engineering Departments. Educ Train 50: 615–625. https://doi.org/10.1108/00400910810909063 doi: 10.1108/00400910810909063
![]() |
[58] |
Hartikainen H, Ventä-Olkkonen L, Kinnula M, et al. (2023) We were proud of our idea: How teens and teachers gained value in an entrepreneurship and making project. Int J Child-Comput Interact 35: 100552. https://doi.org/10.1016/j.ijcci.2022.100552 doi: 10.1016/j.ijcci.2022.100552
![]() |
[59] |
Hassi A (2016) Effectiveness of early entrepreneurship education at the primary school level: Evidence from a field research in Morocco. Citizenship Soc Econ Educ 15: 83–103. https://doi.org/10.1177/2047173416650448 doi: 10.1177/2047173416650448
![]() |
[60] |
Heinonen J (2007) An entrepreneurial-directed approach to teaching corporate entrepreneurship at university level. Educ Train 49: 310–324. https://doi.org/10.1108/00400910710754453 doi: 10.1108/00400910710754453
![]() |
[61] |
Hjorth D (2017) Critique nouvelle–an essay on affirmative-performative entrepreneurship research. Revue de l'Entrepreneuriat 16: 47–54. https://doi.org/10.3917/entre.161.0047 doi: 10.3917/entre.161.0047
![]() |
[62] |
Höglund L, Linton G (2018) Smart specialization in regional innovation systems: a quadruple helix perspective. R D Manag 48: 60–72. https://doi.org/10.1111/radm.12306 doi: 10.1111/radm.12306
![]() |
[63] |
Ireland RD, Hitt MA, Sirmon DG (2003) A model of strategic entrepreneurship: The construct and its dimensions. J Manage 29: 963–989. https://doi.org/10.1016/S0149-2063(03)00086-2 doi: 10.1016/S0149-2063(03)00086-2
![]() |
[64] | Johannisson B (2010) The agony of the Swedish school when confronted by entrepreneurship. In: Creativity and Innovation: Preconditions for Entrepreneurial Education, 91–105. |
[65] | Johannisson B (2018) Limits to and prospects of entrepreneurship education in the academic context, In: Fayolle, A., A Research Agenda for Entrepreneurship Education, Edward Elgar Publishing, 139–163. https://doi.org/10.4337/9781786432919.00014 |
[66] |
Jones C (2019) A signature pedagogy for entrepreneurship education. J Small Bus Enterp Dev 26: 243–254. https://doi.org/10.1108/JSBED-03-2018-0080 doi: 10.1108/JSBED-03-2018-0080
![]() |
[67] | Kanselaar G (2002) Constructivism and socio-constructivism. Constructivism and socio-constructivism, 1–7. |
[68] |
Kautonen T, Van Gelderen M, Fink M (2015) Robustness of the theory of Planned Behavior in Predicting Entrepreneurial Intentions and Actions. Entrep Theory Pract 39: 655–674. https://doi.org/10.1111/etap.12056 doi: 10.1111/etap.12056
![]() |
[69] |
Kibler E, Kautonen T, Fink M (2014) Regional social legitimacy of entrepreneurship: Implications for entrepreneurial intention and start-up behaviour. Reg Stud 48: 995–1015. https://doi.org/10.1080/00343404.2013.851373 doi: 10.1080/00343404.2013.851373
![]() |
[70] | Komulainen K, Naskali P, Korhonen M, et al. (2011) Internal Entrepreneurship-a Trojan horse of the neoliberal governance of education? Finnish pre-and in-service teachers' implementation of and resistance towards entrepreneurship education. J Crit Educ Policy Stud 9: 341–374. |
[71] |
Kuratko DF (2005) The emergence of entrepreneurship education: Development, trends, and challenges. Entrep Theory Pract 29: 577–597. https://doi.org/10.1111/j.1540-6520.2005.00099.x doi: 10.1111/j.1540-6520.2005.00099.x
![]() |
[72] |
Lackéus M (2014) An emotion based approach to assessing entrepreneurial education. Int J Managt Educ 12: 374–396. https://doi.org/10.1016/j.ijme.2014.06.005 doi: 10.1016/j.ijme.2014.06.005
![]() |
[73] | Lackéus M (2015) Entrepreneurship in education: What, why, when, how. Available from: https://www.oecd.org/cfe/leed/BGP_Entrepreneurship-in-Education.pdf. |
[74] |
Lackéus M, Lundqvist M, Middleton KW (2016) Bridging the traditional-progressive education rift through entrepreneurship. Int J Entrep Behav Res 22: 777–803. https://doi.org/10.1108/IJEBR-03-2016-0072 doi: 10.1108/IJEBR-03-2016-0072
![]() |
[75] | Lackéus M, Lundqvist M, Williams Middleton K (2013) How can entrepreneurship bridge between traditional and progressive education? Available from: file: ///C: /Users/zhuan/Downloads/Howcanentrepreneurshipbridgebetweentraditionalandprogressiveeducation-FINAL.pdf. |
[76] |
Langley A (1999) Strategies for theorizing from process data. Acad Manage Rev 24: 691–710. https://doi.org/10.5465/amr.1999.2553248 doi: 10.5465/amr.1999.2553248
![]() |
[77] |
Lanivich SE, Bennett A, Kessler SR, et al. (2021) RICH with well-being: An entrepreneurial mindset for thriving in early-stage entrepreneurship. J Bus Res 124: 571–580. https://doi.org/10.1016/j.jbusres.2020.10.036 doi: 10.1016/j.jbusres.2020.10.036
![]() |
[78] |
Leonard-Barton D (1990) A dual methodology for case studies: Synergistic use of a longitudinal single site with replicated multiple sites. Organ Sci 1: 248–266. https://doi.org/10.1287/orsc.1.3.248 doi: 10.1287/orsc.1.3.248
![]() |
[79] | Lindner J (2020) Entrepreneurship Education by Youth Start—Entrepreneurial Challenge-Based Learning. In: Auer, M., Tsiatsos, T., The Challenges of the Digital Transformation in Education. ICL 2018. Advances in Intelligent Systems and Computing, Springer, Cham, 916: 866–875. https://doi.org/10.1007/978-3-030-11932-4_80 |
[80] | Löw M (2020) Social and Spatial Uncertainty and Inequality: The Refiguration of Spaces as Today's Challenge for Cities. In: Smagacz-Poziemska, M., Gómez, M., Pereira, P., et al., Inequality and Uncertainty, Singapore: Palgrave Macmillan, 23–41. https://doi.org/10.1007/978-981-32-9162-1_2 |
[81] |
McCann JE, SELSKY J (1984) Hyperturbulence and the emergence of type 5 environments. Acad Manage Rev 9: 460–470. https://doi.org/10.5465/amr.1984.4279678 doi: 10.5465/amr.1984.4279678
![]() |
[82] | McGrath RG, MACMILLAN IC (2000) The entrepreneurial mindset: Strategies for continuously creating opportunity in an age of uncertainty, Harvard Business Press. |
[83] | McGuinness C (1999) From thinking skills to thinking classrooms: A review and evaluation of approaches for developing pupils' thinking. Department for Education and Employment, London. https://doi.org/10.13140/RG.2.1.4000.1129 |
[84] | Moberg K (2014a) Assessing the impact of entrepreneurship education: From ABC to PhD, Frederiksberg: Copenhagen Business School (CBS). |
[85] |
Moberg K (2014b) Two approaches to entrepreneurship education: The different effects of education for and through entrepreneurship at the lower secondary level. Int J Manag Educ 12: 512–528. https://doi.org/10.1016/j.ijme.2014.05.002 doi: 10.1016/j.ijme.2014.05.002
![]() |
[86] |
Morawska-Jancelewicz J (2022) The role of universities in social innovation within quadruple/quintuple helix model: Practical implications from polish experience. J Knowl Econ 13: 2230–2271. https://doi.org/10.1007/s13132-021-00804-y doi: 10.1007/s13132-021-00804-y
![]() |
[87] | Motta VF, Galina SVR (2023) Experiential learning in entrepreneurship education: A systematic literature review. Teach Teach Educ 121: 103919. https://doi.org/10.1016/j.tate.2022.103919 |
[88] |
Mucci N (2019) Social Work and Social Welfare from a multidisciplinary perspective–A challenge for researchers. Soc Work Soc Welf 1: 1–3. https://doi.org/10.25082/SWSW.2018.01.002 doi: 10.25082/SWSW.2018.01.002
![]() |
[89] | Mulvey J (2016) Educating communities, In: Mulvey, J., Cooper, B.(Eds.), Understanding the power and politics of public education: Implementing policies to achieve equal opportunity for all, Rowman & Littlefield Publishers, 81–91. |
[90] |
Mwasalwiba ES (2010) Entrepreneurship education: a review of its objectives, teaching methods, and impact indicators. Educ Train 52: 20–47. https://doi.org/10.1108/00400911011017663 doi: 10.1108/00400911011017663
![]() |
[91] | Ngwoke AN, Aneke AO, Oraelosi CC (2020) Developing entrepreneurial mindset in the Nigerian child through quality early childhood education. Int J Stud Educ 16: 162–174. |
[92] |
Nicolaou N, Shane S (2009) Can genetic factors influence the likelihood of engaging in entrepreneurial activity? J Bus Ventur 24: 1–22. https://doi.org/10.1016/j.jbusvent.2007.11.003 doi: 10.1016/j.jbusvent.2007.11.003
![]() |
[93] |
Obschonka M, Silbereisen RK, Schmitt-Rodermund E, et al. (2011) Nascent entrepreneurship and the developing individual: Early entrepreneurial competence in adolescence and venture creation success during the career. J Vocat Behav 79: 121–133. https://doi.org/10.1016/j.jvb.2010.12.005 doi: 10.1016/j.jvb.2010.12.005
![]() |
[94] |
Olutuase SO, Brijlal P, Yan B (2023) Model for stimulating entrepreneurial skills through entrepreneurship education in an African context. J Small Bus Entrep 35: 263–283. https://doi.org/10.1080/08276331.2020.1786645 doi: 10.1080/08276331.2020.1786645
![]() |
[95] |
Oosterbeek H, Van Praag M, Ijsselstein A (2010) The impact of entrepreneurship education on entrepreneurship skills and motivation. Eur Econ Rev 54: 442–454. https://doi.org/10.1016/j.euroecorev.2009.08.002 doi: 10.1016/j.euroecorev.2009.08.002
![]() |
[96] | Patrício LD, Ferreira JJ (2023) Unlocking the connection between education, entrepreneurial mindset, and social values in entrepreneurial activity development. Rev Manag Sci, 1–23. https://doi.org/10.1007/s11846-023-00629-w |
[97] | Palmér H, Johansson M, Karlsson L (2018) Teaching for entrepreneurial and mathematical competencies: teachers stepping out of their comfort zone, In: Palmér, H., Skott, J., Students' and Teachers' Values, Attitudes, Feelings and Beliefs in Mathematics Classrooms, Springer, Cham, 13–23. https://doi.org/10.1007/978-3-319-70244-5_2 |
[98] |
Palumbo R, Manna R (2019) Making Educational Organizations Able to Change: A Literature Review. Int J Educ Manag 33: 734–752. https://doi.org/10.1108/IJEM-02-2018-0051 doi: 10.1108/IJEM-02-2018-0051
![]() |
[99] |
Pepin M (2012) Enterprise education: a Deweyan perspective. Educ Train 54: 801–812. https://doi.org/10.1108/00400911211274891 doi: 10.1108/00400911211274891
![]() |
[100] | Piaget J (1936) Origin of intelligence in the child, Penguin books, New York. |
[101] | Piaget J (1975) L'équilibration des structures cognitives: problème central du développement, Presses universitaires de France. |
[102] | Piaget J (1976) Piaget's theory, In: Inhelder, B., Chipman, H.H., Zwingmann, C. (eds) Piaget and His School. Springer Study Edition, Berlin Heidelberg: Springer, 11–23. https://doi.org/10.1007/978-3-642-46323-5_2 |
[103] |
Pittaway L, Cope J (2007) Simulating entrepreneurial learning: Integrating experiential and collaborative approaches to learning. Manage Learn 38: 211–233. https://doi.org/10.1177/1350507607075776 doi: 10.1177/1350507607075776
![]() |
[104] |
Rae D (2010) Universities and enterprise education: responding to the challenges of the new era. J Small Bus Enterp Dev 17: 591–606. https://doi.org/10.1108/14626001011088741 doi: 10.1108/14626001011088741
![]() |
[105] |
Rauch A, Hulsink W (2015) Putting entrepreneurship Education where the intention to Act lies: An investigation into the impact of entrepreneurship education on entrepreneurial behavior. Acad Manag Learn Educ 14: 187–204. https://doi.org/10.5465/amle.2012.0293 doi: 10.5465/amle.2012.0293
![]() |
[106] | Rebernik M (2009) Quadruple helix of entrepreneurship and management education. Revista de Management Comparat Internațional 10: 910–921. |
[107] |
Rodriguez S, Lieber H (2020) Relationship between entrepreneurship education, entrepreneurial mindset, and career readiness in secondary students. J Experient Educ 43: 277–298. https://doi.org/10.1177/1053825920919462 doi: 10.1177/1053825920919462
![]() |
[108] |
Rosendahl Huber L, Sloof R, Van Praag M (2014) The effect of early entrepreneurship education: Evidence from a randomized field experiment. Eur Econ Rev 72: 76–97. https://doi.org/10.1016/j.euroecorev.2014.09.002 doi: 10.1016/j.euroecorev.2014.09.002
![]() |
[109] | Rustiana Y, Mohd OB, Binti MN (2022) The distribution role of entrepreneurship mindset and task technology fit: an extended model of theory of planned behavior. J Distrib Sci 20: 85–96. http://dx.doi.org/10.15722/jds.20.05.202205.85 |
[110] |
Sahlberg P, Oldroyd D (2010) Pedagogy for economic competitiveness and sustainable development. Eur J Educ 45: 280–299. https://doi.org/10.1111/j.1465-3435.2010.01429.x doi: 10.1111/j.1465-3435.2010.01429.x
![]() |
[111] |
Sampene AK, Li C, Khan A, et al. (2022) Yes! I want to be an entrepreneur: a study on university students'entrepreneurship intentions through the theory of planned behavior. Curr Psychol 1: 1–19. https://doi.org/10.1007/s12144-022-03161-4 doi: 10.1007/s12144-022-03161-4
![]() |
[112] |
San Tan S, Ng CF (2006) A problem-based learning approach to entrepreneurship education. Educ Train 48: 416–428. https://doi.org/10.1108/00400910610692606 doi: 10.1108/00400910610692606
![]() |
[113] | Sánchez JC (2013) The Impact of an Entrepreneurship Education Program on Entrepreneurial Competencies and Intention. J Small Bus Manag 51: 447–465. |
[114] |
Sarasvathy SD, Venkataraman S (2011) Entrepreneurship as method: Open questions for an entrepreneurial future. Entrep Theory Pract 35: 113–135. https://doi.org/10.1111/j.1540-6520.2010.00425.x doi: 10.1111/j.1540-6520.2010.00425.x
![]() |
[115] |
Shayer M, Adhami M (2007) Fostering cognitive development through the context of mathematics: Results of the CAME project. Educ Stud Math 64: 265–291. https://doi.org/10.1007/s10649-006-9037-1 doi: 10.1007/s10649-006-9037-1
![]() |
[116] |
Short JC, Ketchen Jr DJ (2005) Teaching timeless truths through classic literature: Aesop's Fables and strategic management. J Manag Educ 29: 816–832. https://doi.org/10.1177/1052562904271450 doi: 10.1177/1052562904271450
![]() |
[117] |
Smith R, Neergaard H (2015) Telling business stories as fellowship-tales. Int J Gend and Entrep 7: 232–252. https://doi.org/10.1108/IJGE-08-2014-0026 doi: 10.1108/IJGE-08-2014-0026
![]() |
[118] |
Somprach K, Popoonsak P, Sombatteera S (2014) Soft skills development to enhance teachers' competencies in primary schools. Procedia-Soc behav Sci 112: 842–846. https://doi.org/10.1016/j.sbspro.2014.01.1240 doi: 10.1016/j.sbspro.2014.01.1240
![]() |
[119] | Surlemont B (2007) Promoting enterprising: a strategic move to get schools' cooperation in the promotion of entrepreneurship, In: Fayolle, A (ed.), Handbook of Research in Entrepreneurship Education, Edward Elgar Publishing, 2: 255–265. |
[120] |
Volkmann C, Fichter K, Klofsten M, et al. (2021) Sustainable entrepreneurial ecosystems: An emerging field of research. Small Bus Econ 56: 1047–1055. https://doi.org/10.1007/s11187-019-00253-7 doi: 10.1007/s11187-019-00253-7
![]() |
[121] |
Von Glasersfeld E (2001) Radical constructivism and teaching. Prospects 31: 161–173. https://doi.org/10.1007/BF03220058 doi: 10.1007/BF03220058
![]() |
[122] |
Vorley T, Nelles J (2008) (Re) conceptualising the academy. High Educ Manag Policy 20: 1–17. https://doi.org/10.1787/17269822 doi: 10.1787/17269822
![]() |
[123] |
Vorley T, Nelles J (2009) Building entrepreneurial architectures: A conceptual interpretation of the third mission. Policy Futures Educ 7: 284–296. https://doi.org/10.2304/pfie.2009.7.3.284 doi: 10.2304/pfie.2009.7.3.284
![]() |
[124] | Vygotsky L (1962) Thought and language, MIT Press. https://doi.org/10.1037/11193-000 |
[125] | Vygotsky L (1978a) Interaction between learning and development. Readings on the development of children 23: 34–41. |
[126] | Vygotsky L (1978b) Internalization of higher psychological functions. Mind in society: The development of higher psychological processes 37: 52–57. |
[127] | Woods P, Jeffrey B, Troman G, et al. (2019) Restructuring schools, reconstructing teachers: Responding to change in the primary school, Routledge. |
[128] |
Yan J, Huang T, Xiao Y (2022) Assessing the impact of entrepreneurial education activity on entrepreneurial intention and behavior: role of behavioral entrepreneurial mindset. Environ Sci Pollut Res 30: 26292–26307. https://doi.org/10.1007/s11356-022-23878-w doi: 10.1007/s11356-022-23878-w
![]() |
[129] | Yin RK (1994) Case study research (3 ed.), Sage, Thousand Oaks, CA. |
[130] |
Zupan B, Cankar F, Setnikar Cankar S (2018) The development of an entrepreneurial mindset in primary education. Eur J Educ 53: 427–439. https://doi.org/10.1111/ejed.12293 doi: 10.1111/ejed.12293
![]() |