Natural rivers connect to each other to form networks. The geometric structure of a river network can significantly influence spatial dynamics of populations in the system. We consider a process-oriented model to describe population dynamics in river networks of trees, establish the fundamental theories of the corresponding parabolic problems and elliptic problems, derive the persistence threshold by using the principal eigenvalue of the eigenvalue problem, and define the net reproductive rate to describe population persistence or extinction. By virtue of numerical simulations, we investigate the effects of hydrological, physical, and biological factors, especially the structure of the river network, on population persistence. physical and hydrological features in a river network, can greatly influence the spatial distribution of the flow profile (including the flow velocity, water depth, etc.) in branches of the network. Moreover, the population dispersal vectors may be constrained by the network configuration and the flow profile, and species life history traits may depend on varying habitat conditions in the network. As a result, population distribution and persistence in river systems can be significantly affected by the network topology or structure, see e.g., [8,11,12,[14][15][16]46]. Then there arise interesting questions such as whether a population can persist in the desert streams of the southwestern United State while the streams are experiencing substantial natural drying trends [11], or whether dendritic geometry enhances dynamic stability of ecological systems [15] etc. Furthermore, other related dynamics in the network, such as the dynamics of water-born infectious diseases like cholera may also be greatly affected by the river network geometry (see e.g., [6]).Branches in a river network have been modeled as point nodes in a network of habitats in individual based models [11,16] and matrix population models for stage-structured populations [14,39]. However this oversimplifies the spatial heterogeneity of river networks. In a real river ecosystem, organisms mainly live in the branches of the network and the connections between branches (e.g., the network nodes) are mainly for population transitions from one branch to another. To take into account this realistic situation, in recent works [48-50], integro-differential equations and reactiondiffusion-advection equations were used to model population dynamics in river networks where the network branches, instead of the network nodes, are the main habitat for organisms. Here the river networks are modeled under the framework of metric graphs (or metric networks). A metric graph is a graph G = (V, E) with a set V of vertices and a set E of edges, such that each edge e ∈ E is associated with either a closed bounded interval. Mathematic notion of metric graphs was first introduced in the context of wave propagation on thin graph-like domains [5,25], and they are also called quantum graphs.The theories of parabolic and elliptic equations as well as the corr...