This work studies a nonlinear inverse problem of reconstructing the diffusion coefficient in a parabolic-elliptic system using the final measurement data, which has important application in a large field of applied science. Being different from other works, which are governed by single partial differential equations, the underlying mathematical model in this paper is a coupled parabolic-elliptic system, which makes theoretical analysis rather difficult. On the basis of the optimal control framework, the identification problem is transformed into an optimization problem. Then the existence of the minimizer is proved, and the necessary condition that must be satisfied by the minimizer is also given. Since the optimal control problem is nonconvex, one may not expect a unique solution universally. However, the local uniqueness and stability of the minimizer are deduced in this paper. DENG ET AL 3415 coefficient problem of identification of the diffusion coefficient from the final overspecified data has been considered by several authors (see, for instance, Demir and Hasanov 22,23 ).Compared with other papers, which also deal with parameter identification problems, the results to be given in this paper have the following features. Firstly, the mathematical model proposed in this paper is a coupled parabolic-elliptic system rather than a single parabolic or elliptic equation. So this kind of problem indeed belongs to the boundary value problems of mixed type. Secondly, our task is to identify the diffusion coefficient in the coupled system, and this inverse problem generally belongs to the mathematical class of severely ill-posed problems. Furthermore, though the underlying mathematical model is linear, the corresponding inverse problem is completely nonlinear. Thirdly, since the exchange of matter and energy occurs on the interface, and after all, the connective condition is concerned with the diffusivity, the underlying mathematical model can be viewed as a nonlinear boundary-value problem combined with a nonlinear parameter-identification problem. From this point of view, such problem is rather difficult than those of single equation. Finally, the extra condition used in the paper is only imposed on the solution in the small domain rather than on the whole domain Ω. Because of the lack of measurement data, we shall carefully treat every integral to avoid using more information outside . In summary, the complexity of the inverse problem 1.2 to be introduced in following Section 3 is higher than those that are governed by single equations.In this work, we use an optimal control framework (see, eg, recent works 7,21,24,25 ) to discuss Problem 1.1 mainly from the theoretical analysis angle. As mentioned above, since the mathematical model in this paper is a coupled parabolic-elliptic system rather than the single parabolic equation in previous studies, 7,21,24,25 the major differences between them should be noted. It is well known that the properties of the solution for elliptic or parabolic system are quite different f...